结果不收敛怎么办

光语者

算出的结果不收敛怎么办？


命令流如下：




/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