求教箱梁极限强度分析命令流
我做了箱型梁的一半，中间做对称条件，梁的两个端面上各作刚性梁耦合到一个节点上，然后加弯矩，想做出它的最大承载力，现在压根不收敛，不知道是个啥问题
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
页:
[1]