结果不收敛怎么办
算出的结果不收敛怎么办？命令流如下：
/filname,TYX
/prep7
*afun,deg !定义角度单位为度
!确定基坑参数
ar_h=0.98
ar_d=12
h1=-1.5
h2=-4.5
h3=-16.5
sh=h1+h2+h3+ar_d
o=18.4548
l=25.831
h=5.7863
R_i=l/2/sin(o)
a=l/2-h*sin(o)/cos(o)/2+h/2/cos(o)
b=R_i-(l/2-h*sin(o)/cos(o)/2)*cos(o)/sin(o)
x1=l/2-h*sin(o)/(2*cos(o))
hkw=-0.75
!定义材料参数，三层土，一种水泥墙
et,1,solid45
et,101,45
et,102,45
et,2,solid45
et,201,solid45
et,202,solid45
et,203,solid45
et,204,solid45
et,205,solid45
et,206,solid45
et,3,solid45
et,301,solid45
et,302,solid45
!第一层
mp,ex,1,10e6
mp,nuxy,1,0.4
mp,dens,1,19572
tb,dp,1
tbdata,1,20.1e3,16.8,30
!第二层
mp,ex,2,10e6
mp,nuxy,2,0.4
mp,dens,2,17839
tb,dp,2
tbdata,1,13.2e3,11.3,30
!第三层
mp,ex,3,10e6
mp,nuxy,3,0.4
mp,dens,3,17431.2
tb,dp,3
tbdata,1,11.3e3,9.6,30
!拱
et,4,solid45
mp,ex,4,6e8
mp,prxy,4,0.25
mp,dens,4,21600
!面单元
et,5,plane42
!建立关键点
csys,0
k,1,0,0
k,3,l/2,h/2
k,5,0,b
k,6,0,b+ar_h
k,7,a,0
k,8,a+ar_h,0
k,9,0,3*b
k,10,3*a,3*b
k,11,3*a,0
k,13,a+5*ar_h,0
k,14,a+5*ar_h,b+5*ar_h
k,15,0,b+5*ar_h
k,16,3*a,b+5*ar_h
k,17,a+5*ar_h,3*b
LOCAL,11,1,0,-(R_i-b),0,,,,,,
csys,11
k,12,R_i+ar_h,90-o
LOCAL,12,1,x1,0,0,,,,,,
csys,12
l,3,7
l,12,8
csys,11
l,5,3
l,6,12
csys,0
l,1,5
l,1,7
l,7,8
l,8,13
l,13,11
l,5,6
l,6,15
l,15,9
l,13,14
l,15,14
l,11,16
l,16,10
l,9,17
l,16,10
l,14,16
l,14,17
l,17,10
l,3,12
al,1,7,2,21
al,3,21,4,10
al,6,1,3,5
al,13,14,11,4,2,8
al,15,18,13,9
al,16,20,19,18
al,19,17,12,14
lsel,s,,,9
lsel,a,,,18
lsel,a,,,20
lesize,all,,,15,4
lsel,s,,,14,17,3
lesize,all,,,12
lsel,s,,,13,15,2
lesize,all,,,10
lsel,s,,,11
lesize,all,,,5
lsel,s,,,7,10,3
lsel,a,,,21
lesize,all,,,3
lsel,s,,,5
lesize,all,,,6
lsel,s,,,6
lesize,all,,,8
lsel,s,,,8
lesize,all,,,5
lsel,s,,,1,2
lesize,all,,,4
lsel,s,,,3,4
lesize,all,,,12
mshape,0,2D
mshkey,2
type,5
amesh,1,7
eplot
extopt,esize,h1/hkw,0
type,101
mat,1
vext,3,,,,,h1/2
asel,s,loc,z,h1/2
extopt,esize,h1/hkw,0
type,102
mat,1
vext,all,,,,,h1/2
asel,s,loc,z,h1
extopt,esize,h1/hkw,0
type,201
mat,2
vext,all,,,,,h2/6
asel,s,loc,z,h1+hkw
extopt,esize,h1/hkw,0
type,202
mat,2
vext,all,,,,,h2/6
asel,s,loc,z,h1+2*hkw
extopt,esize,h1/hkw,0
type,203
mat,2
vext,all,,,,,h2/6
asel,s,loc,z,h1+3*hkw
extopt,esize,h1/hkw,0
type,204
mat,2
vext,all,,,,,h2/6
asel,s,loc,z,h1+4*hkw
extopt,esize,h1/hkw,0
type,205
mat,2
vext,all,,,,,h2/6
asel,s,loc,z,h1+5*hkw
extopt,esize,h1/hkw,0
type,206
mat,2
vext,all,,,,,h2/6
asel,s,loc,z,h1+h2
extopt,esize,h1/hkw,0
type,301
mat,3
vext,all,,,,,hkw
asel,s,loc,z,h1+h2+hkw
extopt,esize,h1/hkw,0
type,302
mat,3
vext,all,,,,,hkw
asel,s,loc,z,h1+h2+2*hkw
extopt,esize,12,0
type,3
mat,3
vext,all,,,,,-4.5
asel,s,,,1,2
extopt,esize,-2*ar_d/hkw,0
type,4
mat,4
vext,1,2,,,,-ar_d
asel,s,,,4,7
extopt,esize,2*h1/hkw,0
type,1
mat,1
vext,4,7,,,,h1
asel,s,,,72,79,7
asel,a,,,83,87,4
extopt,esize,2*h2/hkw,0
type,2
mat,2
vext,all,,,,,h2
asel,s,,,90,97,7
asel,a,,,101,105,4
extopt,esize,16,0
type,3
mat,3
vext,all,,,,,2*hkw-4.5
asel,s,loc,z,h1+h2+2*hkw-4.5
extopt,esize,(h3-(h1+h2+2*hkw-4.5))/hkw,0
type,3
mat,3
vext,all,,,,,h3-(h1+h2+2*hkw-4.5)
nummrg,node
/solu
DOFSEL,S,DISP !约束叠加,以免后面的约束把前面的取消掉了
DCUM,ADD,
DOFSEL,ALL
!**************设置约束***************
!选择对称的面。
nsel,s,loc,y,0
dsym,symm,y
nsel,s,loc,x,0
dsym,symm,x !定义上面选择的面为对称约束
allsel
!**************设置固定面边界***************
!选择固定的面
nsel,s,loc,x,3*a
nsel,a,loc,y,3*b
nsel,a,loc,z,h1+h2+h3
d,all,all
allsel
ANTYPE,0
NLGEOM,1
/AUTO,1
/REP,FAST
NSUBST,300,1000,50
AUTOTS,1
EQSLV,PCG,1E-4
LNSRCH,1
NEQIT,50
PRED,ON,,ON
TIME,1
acel,,,1
NLGEOM,1
NROPT,FULL, ,ON
LUMPM,0
EQSLV,PCG, ,3,
PRECISION,0
MSAVE,0
PIVCHECK,1
SSTIF
PSTRES
TOFFST,0,
allsel
kbc,0
allsel
allsel
allsel
EPLOT
allsel
kbc,0
solve
页:
[1]