chinact6 发表于 2006-7-24 17:21

求助:求解不收敛

*set,h,100                  !板高
*set,b,560                  !板宽
*set,l,3270                   !板长
*set,lz,100                   !支座长
*set,zg,20                  !支座高
*set,a,25                     !保护层厚度
*set,dis,120                  !钢筋间距
*set,zj,8                                                                                        !钢筋直径
*set,pi,acos(-1)            
*set,sr,pi*(zj/2)**2          !钢筋面积   
*set,cb,300                   !碳纤维布宽度
*set,ch,0.103               !碳纤维布厚度
*set,cl,l-2*lz                !碳纤维布长度

*set,f1,3.14e-3       
*set,f2,5.58e-3

/prep7
et,1,solid65                  !设置混凝土单元
et,2,link8                  !设置钢筋单元
et,3,shell41                  !设置碳纤维布单元
keyopt,3,1,2                  
et,4,solid45                  !设置支座单元

mp,ex,1,3.1e4               
mp,prxy,1,0.2

tb,kinh,1,1,7               
tbpt,,0.0001,3.1
tbpt,,0.0004,11.477
tbpt,,0.0008,20.40
tbpt,,0.0012,26.79
tbpt,,0.0016,30.60
tbpt,,0.002,31.88
tbpt,,0.0033,31.88

tb,conc,1,1,9                                                                                             
tbdata,,0.3,0.5,1.75,-1

mp,ex,2,2e5                                 
mp,prxy,2,0.3
tb,bkin,2,1,2,1
tbdata,,265,0

mp,ex,3,2.8e5
mp,prxy,3,0

mp,ex,4,2.1e5
mp,prxy,4,0.3

r,1,2,0,0,0,,,                        
rmore,,,,,,,

r,2,sr,,
r,3,sr/2,,

r,4,ch,,,,,,
rmore,,

block,0,-b/2,0,h,(l/2-lz),-(l/2-lz)   


lgen,2,8,,,b/2-2*dis                  
lgen,2,13,,,dis
adrag,13,,,,,,9
adrag,14,,,,,,9
vsba,1,7
vsba,3,8
lsel,s,loc,z,(l/2-lz)
lsel,r,loc,y,0
lgen,2,all,,,,a
lsel,s,loc,z,(l/2-lz)
lsel,r,loc,y,a
adrag,all,,,,,,9
allsel,all
vsba,2,15
vsba,4,14
vsba,1,13

lsel,s,loc,y,0                     !划分混凝土单元
lsel,a,loc,y,a
lsel,a,loc,y,h
lsel,u,loc,z,(l/2-lz)
lsel,u,loc,z,-(l/2-lz)
lesize,all,,,50

lsel,s,loc,y,0               
lsel,a,loc,y,a
lsel,a,loc,y,h
lsel,u,loc,x,0
lsel,u,loc,x,-dis
lsel,u,loc,x,-2*dis
lsel,u,loc,x,-b/2
lsel,r,loc,x,-2*dis,0
lesize,all,30
allsel
lsel,s,loc,y,0               
lsel,a,loc,y,a
lsel,a,loc,y,h
lsel,u,loc,x,0
lsel,u,loc,x,-dis
lsel,u,loc,x,-2*dis
lsel,u,loc,x,-b/2
lsel,r,loc,x,-b/2,-2*dis
lesize,all,40

lsel,s,loc,z,(l/2-lz)
lsel,a,loc,z,-(l/2-lz)
lsel,u,loc,y,0
lsel,u,loc,y,a
lsel,u,loc,y,h
lesize,all,25
allsel,all

type,1
mat,1
real,1
esys,0
mshape,0,3d
mshkey,1
vmesh,all
allsel,all

asel,s,loc,z,(l/2-lz)         !生成左支座上混凝土单元
type,1
extopt,esize,3,0
extopt,aclear,1
extopt,attr,0,0,0
mat,1
real,1
esys,0
vext,all,,,0,0,lz

asel,s,loc,z,-(l/2-lz)         !生成右支座上混凝土单元
type,1
extopt,esize,3,0
extopt,aclear,1
extopt,attr,0,0,0
mat,1
real,1
esys,0
vext,all,,,0,0,-lz
allsel,all

lsel,s,loc,x,-2*dis                     !划分板计算长度内钢筋单元
lsel,a,loc,x,-dis
lsel,r,loc,y,a
lsel,r,loc,z,-(l/2-lz),(l/2-lz)
type,2
mat,2
real,2
esys,0
lesize,all,,,,50
lmesh,all
allsel,all

lsel,s,loc,x,0
lsel,r,loc,y,a
lsel,r,loc,z,-(l/2-lz),(l/2-lz)
type,2
mat,2
real,3
esys,0
lesize,all,,,,50
lmesh,all
allsel,all

lsel,s,loc,x,-2*dis                     !划分板支座上钢筋单元
lsel,a,loc,x,-dis
lsel,r,loc,y,a
lsel,r,loc,z,(l/2-lz),l/2
type,2
mat,2
real,2
esys,0
lesize,all,,,,3
lmesh,all
allsel,all

lsel,s,loc,x,0
lsel,r,loc,y,a
lsel,r,loc,z,-l/2,-(l/2-lz)
type,2
mat,2
real,3
esys,0
lesize,all,,,,3
lmesh,all
allsel,all

asel,s,loc,y,0                        !生成支座
asel,u,loc,z,-(l/2-lz),(l/2-lz)
type,4
extopt,esize,2,0
extopt,aclear,1
extopt,attr,0,0,0
mat,4
esys,0
vext,all,,,,-zg
allsel,all

allsel,all                                    !碳纤维布
wpstyle,,,,,,,,1
wprota,,-90
rectng,0,-cb/2,(l/2-lz),-(l/2-lz)
wpstyle,,,,,,,,0
allsel,all
*get,mian,area,,num,max                     !划分碳纤维布
asel,s,area,,mian
lsla,s
lsel,u,loc,z,(l/2-lz)
lsel,u,loc,z,-(l/2-lz)
lesize,all,,,50

allsel,all
asel,s,area,,mian
lsla,s
lsel,u,loc,x,0
lsel,u,loc,x,-cb/2
lesize,all,30

allsel,all                  
asel,s,area,,mian
type,3
mat,3
real,4
esys,0
mshape,1,2d                     
amesh,all

nummrg,node            
numcmp,node                     

/solu
nsel,s,loc,x,0
dsym,symm,x                     

allsel,all                     
nsel,s,loc,z,l/2-lz
nsel,r,loc,y,-zg
d,all,,,,,,uy,uz

allsel,all
nsel,s,loc,z,-(l/2-lz)
nsel,r,loc,y,-zg
d,all,,,,,,uy

allsel,all
asel,s,loc,y,h
asel,r,loc,z,-(l/2-lz),(l/2-lz)
sfa,all,1,pres,f1

antype,0
time,1
nlgeom,1
autots,1
nropt,full
eqslv,spar
outres,all,all
nsubst,20
cnvtol,f,,0.05,2         
solve

把子布等改了好多次,就是不收敛,请各位大侠帮帮我

xiven2008 发表于 2006-7-24 21:06

光改子步数是没用的,看看收敛准则里的允许误差、求解器等,反正能改的都试试

混凝土方面的我不懂的,爱莫能助!
页: [1]
查看完整版本: 求助:求解不收敛