马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
小弟对一个薄壁墩顺利进行特征值稳定分析之后求解非线性屈曲极限荷载,一直循环不止,求教!命令流如下:
hm=1
w1=0.60
w2=0.40
w3=0.60
m=50
/filnam,sheel-pole
/view,1,1,1,1
/prep7
!国际单位制
!设计常量
h=50
wh=12 !横向宽
wz=2.5 !侧面顶宽
wzb=wz+2*h/m !侧面底宽
k,1,0,h-hm/2,0
k,2,-wz/2+w3/2,h-hm/2,0
k,3,-wz/2-(h-hm/2)/m+w3/2,0,0
kgen,2,1,3,1,,,wh/2-w1/2
k,7,0,0,0
k,8,0,0,wh/2-w1/2
a,1,2,5,4
a,2,5,6,3
a,4,5,6,8
/PREP7
et,1,shell63,,1 !keopt(2)=1
mp,dense,1,2650
mp,ex,1,3.45e10
mp,prxy,1,0.1667
r,1,hm/2
r,2,w3
r,3,w2
r,4,w1
et,2,shell63 !刚性板
mp,dense,2,7850
mp,ex,2,3.45e12
mp,prxy,2,0.3
r,5,3
smrt,8 !应该指定分网精细度
allsel
asel,s,,,1
aatt,1,1,1
mshape,0,2d
mshkey,2
amesh,all
asel,s,,,3
aatt,1,4,1
mshkey,2
amesh,all
asel,s,,,2
aatt,1,2,1
mshkey,2
amesh,all
allsel
nsym,x,210,all
allsel
esym,,210,all
nummrg,node
allsel
!nsel,u,loc,z,0,0.01
nsym,z,420,all
allsel
esym,,420,all
nummrg,node
numcmp,node
nummrg,elem
numcmp,elem
a,1,2,3,7
asel,s,,,4
aatt,1,3,1
mshkey,2
amesh,all
!施加墩底约束
esel,s,type,,1
nsle,s,all
csys,0
nsel,r,loc,y,0
d,all,all
allsel
!各种工况下的荷载
/solu
antype,static
pstres,on
allsel
nsel,s,loc,y,h-hm/2,,
!acel,,9.8
f,all,fy,-100
allsel
solve
finish
/solu
antype,buckle
bucopt,subsp,5
!expass,on
!mxpand,2
solve
finish
/post1
set,first
*get,eigenload,mode,1,freq
pldisp,2
set,next
pldisp,2
finish
/config,nres,100
/prep7
tb,biso,1,1,2
tbtemp,0
tbdata,,2e8,0
upgeom,0.1,1,1,sheel-pole,rst
finish
/solu
antype,0
nlgeom,1
outres,all,all
arclen,1,0,0
arctrm,u,0.25e-5,node(0,h-hm/2,0),ux
nsubst,20,,,1
solve
finish
/post26
nsol,2,node(0,h-hm/2,0),u,x,deflection
rforce,3,node(0,0,0),f,y,load
/axlab,x,deflection
/axlab,y,load
xvar,2
plvar,3 |