求高手指点我的命令流有何问题

yanfengpang007

要用混凝土做一个沉井，底部有混凝土底板和深梁，单元采用solid65，考虑弥散钢筋，下面是命令流，但计算一直不收敛。希望有前辈可以详细的指点下，如果命令流里有什么问题，也希望可以指正。
finish
/clear
/filname,chenjing-solid
/title,analysis of chenjing-solid


!参数定义
t1=1.5    !沉井壁厚
t2=1      !混凝土底板厚
t3=0.6    !深梁宽

b1=9.5       !沉井外半径
b2=1.5       !开孔半径


!--------前处理---------------------------

/PREP7

!定义单元类型

ET,1,solid65

KEYOPT,1,1,0
KEYOPT,1,5,0
KEYOPT,1,6,0
KEYOPT,1,7,1

!定义实常数

!实参数1：沉井不同方向配筋
R,1,2,0.004,90,,2,0.04,
RMORE,90,90,,,,,             ! 环向配筋率为0.4%，纵向配筋率为4%

R,2,2,0.04,,,2,0.04,
RMORE,90,,,,,,               ! x向配筋率为4%，y向配筋率为4%               

!显示控制

/PNUM,KP,0  
/PNUM,LINE,0
/PNUM,AREA,0
/PNUM,VOLU,1


!定义材料属性
!混凝土基本材料属性
fc=14.3e6
ft=1.43e6
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,3e10
MPDATA,PRXY,1,,0.2
!破坏准则
TB,CONC,1,1,9,
TBTEMP,0
TBDATA,,0.3,0.5,ft,fc,,
TBDATA,,,,1,,,
MPTEMP,,,,,,,,
MPTEMP,1,0
!钢材基本属性
MPDATA,EX,2,,2.06e11
MPDATA,PRXY,2,,0.3
!屈服准则
TB,BISO,2,1,2,
TBTEMP,0
TBDATA,,235e6,,,,,

!创建关键点,创建沉井
k,1,0,0,0
k,2,0,0,30.8
l,1,2
PCIRC,8,9.5,0,360
VDRAG,1,,,,,,1
allsel
vplot


!开孔1
k,19,0,0,8
k,20,15,0,8
l,19,20
wpoffs,,,8
wprota,,,90
PCIRC,1.5,,0,360
VDRAG,11,,,,,,26
vptn,all
vdele,3,5,,1
numcmp,all

!开孔2
k,29,-13.703,-6.101,8
l,19,29
wprota,,156
PCIRC,1.5,,0,360
VDRAG,15,,,,,,41
vptn,all
vdele,3,5,,1
allsel
numcmp,all

!创建混凝土底板
wprota,,-156,-90
wpoffs,,,-3
cyl4,0,0,8,,,,1
allsel

!深梁
wpstyl,defa
a,8,39,40,9
voffst,23,0.3
vsymm,y,3
a,7,38,41,10
voffst,35,0.3
vsymm,x,5
vsel,s,,,2,6
vadd,all
allsel
numcmp,all

vovlap,all
numcmp,all
allsel
vdele,1,8,,1
numcmp,all
vglue,all
allsel

!网格划分
wpoffs,,,12
vsbw,2
numcmp,all

CSWPLA,11,1,1,1
vsel,s,,,2,3
VATT,1,1,1,11  
allsel
vsel,s,,,1
vatt,1,2,1,0
allsel
csys,0

MSHAPE,1,3D
mshkey,0
vsel,s,,,1
esize,0.6
vmesh,all
allsel

vsel,s,,,2
esize,1.5
vmesh,all
allsel

vsel,s,,,3
esize,0.6
vmesh,all
vimp,all,1,1
allsel




!-----------施加约束--------------------
asel,s,,,23,26,1
asel,a,,,59,66,1
da,all,ux
da,all,uy
da,all,uz
allsel
csys,0
SBCTRAN
                    

!--------------求解---------------------
/solu
csys,0
acel,,,9.8
antype,0
!NLGEOM,ON
nropt,full
!arclen,on
outres,all,all
AUTOTS,ON
nsubst,200,800,200

/PSF,PRES,NORM,2,0,1      !主动土压力
csys,1
nsel,s,loc,x,9.5
nsel,r,loc,y,-90,90
nsel,r,loc,z,0,4.8
sfgrad,pres,0,z,0,-16.47e3
sf,all,pres,329.96e3
allsel

csys,1
nsel,s,loc,x,9.5
nsel,r,loc,y,-90,90
nsel,r,loc,z,4.8,9.1
sfgrad,pres,0,z,4.8,-18.16e3
sf,all,pres,250.92e3
allsel

csys,1
nsel,s,loc,x,9.5
nsel,r,loc,y,-90,90
nsel,r,loc,z,9.1,16
sfgrad,pres,0,z,9.1,-11.17e3
sf,all,pres,172.82e3
allsel

csys,1
nsel,s,loc,x,9.5
nsel,r,loc,y,-90,90
nsel,r,loc,z,16,25.48
sfgrad,pres,0,z,16,-10.1e3
sf,all,pres,95.75e3
allsel

!被动土压力
csys,1
nsel,s,loc,x,9.5
nsel,r,loc,y,90,270
nsel,r,loc,z,0,4.8
sfgrad,pres,0,z,0,-21.44e3
sf,all,pres,477.81e3
allsel

csys,1
nsel,s,loc,x,9.5
nsel,r,loc,y,90,270
nsel,r,loc,z,4.8,9.1
sfgrad,pres,0,z,4.8,-21.19e3
sf,all,pres,374.91e3
allsel

csys,1
nsel,s,loc,x,9.5
nsel,r,loc,y,90,270
nsel,r,loc,z,9.1,16
sfgrad,pres,0,z,9.1,-20.20e3
sf,all,pres,283.78e3
allsel

csys,1
nsel,s,loc,x,9.5
nsel,r,loc,y,90,270
nsel,r,loc,z,16,25.48
sfgrad,pres,0,z,16,-15.23e3
sf,all,pres,144.39e3
allsel

!底板浮力
csys,0
asel,s,area,,19,22,1
sfa,all,1,pres,156.55e3
allsel
SBCTRAN

solve
finish

senja

建议楼主看看别的帖子。
有关发帖规则的。
有关不收敛的可能原因分析的，然后自己尝试调