实体结构施加弯矩非线性分析求助

cqwlyjdr

取了实体结构的1/4进行分析，按照论坛上介绍的刚域的方法施加弯矩，可是总是算不出结果，求问是哪里不对呢？
命令流如下：
！建模试验
FINISH$/CLEAR$/FILNAME,joint$/PREP7
ALLSEL
ET,1,SOLID92                          !定义材料
mp,ex,1,206000
mp,nuxy,1,0.3
tb,bkin,1
tbdata,1,235,0.0
et,2,mass21
r,2,2
!定义模型参数
R2=200
R1=190
R3=40
R4=50
L0=3*R4
L1=R4*3
!建立几何模型
n,1,R2+L1,R3,L0/2
CYL4,,,R1,40,R2,90,L0
CYL4,,,R1,,R2,40,L0                  !切分钢管
CYL4,,,R2,,R2,90,-(R2-R1)            !建立封板
WPOFF,,,L0
CYL4,,,R2,,R2,90,(R2-R1)
WPCSYS,-1
WPOFF,,,L0/2
WPROTA,,,90
CYL4,,,R3,0,R4,180,L1+R2
VPTN,ALL
VDELE,7,,,1
ALLSEL
WPCSYS,-1
WPOFF,R2+L1/3$WPROTA,,,90$VSBW,ALL   !切分钢板
WPCSYS,-1
VPTN,ALL

!几何模型完成，开始划分网格
VSEL,S,VOLU,,1,2
VSEL,A,VOLU,,10,11
CM,VOL1,VOLU
MSHAPE,1,3d
MSHKEY,0
ESIZE,(R2-R1)
VMESH,ALL
ALLSEL
CMSEL,U,VOL1
ESIZE,(R2-R1)/2
VMESH,ALL
!施加约束
nsel,S,LOC,X,0
D,all,ux
nsel,S,LOC,Y,0
D,ALL,UY
allsel
NSEL,S,LOC,X,R2+L1
NSEL,R,LOC,Y,0
DDELE,ALL,ALL
TYPE,2
REAL,2
e,1
NSEL,S,LOC,X,R2+L1-0.01,R2+L1+0.01
cerig,1,all,all
allsel
sbctran
!施加荷载
NSEL,S,LOC,X,R2+L1
N=NDINQR(0,13)
LOAD=-300000
LOADMX=LOAD*0.75
F,ALL,FX,LOAD/N
nsel,s,node,,1
F,1,my,-LOADMX
allsel
FINISH
/solu
allsel,all
antype,0
nlgeom,ON
outres,all,all
arclen,ON
nsubst,200
solve
finish

cqwlyjdr

自己顶顶