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

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

自己顶顶