求教箱梁极限强度分析命令流

xie364532537

我做了箱型梁的一半，中间做对称条件，梁的两个端面上各作刚性梁耦合到一个节点上，然后加弯矩，
想做出它的最大承载力，现在压根不收敛，不知道是个啥问题

finish
/clear
M=2E8!弯矩,梁的极限弯矩在1.6e9以下

/PREP7  
ET,1,SHELL181   
  
!ET,2,BEAM188
ET,2,184
KEYOPT,2,1,1
!*  
R,1,4.88, , , , , ,
RMORE, , , , , , ,  
!*  
R,2,3.38,0,0,0,0,0,
RMORE,0,0,0,0,0,0,  
!*  

MPTEMP,,,,,,,,         !第一个材料
MPTEMP,1,0  
MPDATA,EX,1,,208500
MPDATA,PRXY,1,,0.3  

TB,MISO,1,1,3,  

TBTEMP,0
TBPT,,0,0   
TBPT,,298/208500,298
TBPT,,1,298   

MPTEMP,,,,,,,,         !第二个材料
MPTEMP,1,0  
MPDATA,EX,2,,216200
MPDATA,PRXY,2,,0.3  

TB,MISO,2,1,3,  

TBTEMP,0
TBPT,,0,0   
TBPT,,211.6/216200,211.6
TBPT,,1,211.6

MPTEMP,,,,,,,,         !第三个材料
MPTEMP,1,0  
MPDATA,EX,3,,191500
MPDATA,PRXY,3,,0.3  

TB,MISO,3,1,3,  

TBTEMP,0
TBPT,,0,0   
TBPT,,276.5/191500,276.5
TBPT,,1,276.5


K,1,0,0,0,
K,2,609.5,0,0
K,3,609.5,0,914
K,4,0,0,914
K,5,0,787.4,0


L,1,2
L,2,3
L,3,4
L,1,5

LPLOT
/PNUM,LINE,1
/REPLOT


NUMMRG,ALL, , , ,LOW

ADRAG,4,,,,,,1,2,3
NUMMRG,ALL, , , ,LOW
/REPLOT

LESIZE,1, , ,10, , , , ,1  
LESIZE,2, , ,14, , , , ,1  
LESIZE,3, , ,10, , , , ,1  

LESIZE,7, , ,10, , , , ,1  
LESIZE,10, , ,14, , , , ,1  
LESIZE,13, , ,10, , , , ,1  

LESIZE,4, , ,9, , , , ,1  
LESIZE,5, , ,9, , , , ,1  
LESIZE,8, , ,9, , , , ,1  
LESIZE,11, , ,9, , , , ,1  


TYPE,   1   
MAT,    1
REAL,   1   
AMESH,1

AMESH,3


TYPE,   1   
MAT,    2
REAL,   2   
AMESH,2

nsel,s,loc,x,0
nsel,r,loc,y,1,780
D,ALL, , , , , ,UY,ROTX,ROTZ, , ,

allsel

nsel,s,loc,Y,0
n,1000,0,-100,462.732
type,2

*get,nnum,node,0,count
*get,ND,node,0,num,min
*do,i,2,nnum
!!!!nnum个mpc单元
E, 1000,ND
ND=NDNEXT(ND)
*enddo
allsel


nsel,s,loc,Y,787.4
n,1001,0,887.4,462.732
type,2

*get,nnum,node,0,count
*get,ND,node,0,num,min
*do,i,2,nnum
!!!!生成nnum个mpc单元
E, 1001,ND
ND=NDNEXT(ND)
*enddo
allsel

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

/SOLU
ANTYPE,0
NSEL,S,NODE,,1000
D,ALL, , , , , ,UX,UZ,ROTY,ROTZ, ,
F,NODE(0,-100,462.732),MX,-M

NSEL,S,NODE,,1001
D,ALL, , , , , ,UX,UZ,ROTY,ROTZ, ,
F,NODE(0,887.4,462.732),MX,M
nlge,on

arcl,on

nsub,100

outr,all,all

allsel
SOLVE
FINISH  


/POST1
SET,LAST
!PLDISP,1
PLNSOL, S,EQV, 0,1.0
FINISH


/POST26
！载荷时间曲线
ESOL,2,177,136,S,x!梁中间截面上的一个节点
prod,2,1,,,load,,,M/462.732*4.88/1219

prvar,2
/axlab,x,TIME(S)
/axlab,y,total load(Pa)
/grid,1
/xrange,0,1
/yrange,0,5000
plvar,2