马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
帮助看一下我的非线性屈曲分析为什么不收敛啊!!!好郁闷啊!!:'(
finish
/clear,start
/filname,pipeline-5,1
/title,buckling of pipeline
/com,structural
/unit,si
/config,nres,1100000
*set,max_deffect,0.004
*dim,x,,10
*dim,y,,10
*dim,z,,10
/PREP7
l=30
e=6
ET,1,BEAM189
MP,EX,1,2e11
MP,PRXY,1,0.3
!define section data
SECTYPE,1,BEAM, CTUBE,,0
SECOFFSET, CENT
SECDATA,0.295,0.299,,0,0,0,0,0,0,0
!define keypoints
/view,1,1,1,1
k,1,
k,2,0,0,l
l,1,2
LATT, 1, , 1, , , , 1
esize,e
lmesh,all
save
finish
/solu
d,1,ux
d,1,uy
d,1,rotx
d,2,ux
d,2,uy
d,2,uz
allsel
fk,1,fz,1
allsel
pstres,on
outrr,,1
solve
save
finish
/solu
antype,1
bucopt,subsp,1,0,0
subopt,1,
mxpand,1,0,0,0,0.001
solve
finish
/post1
set,list
set,first
*get,pressure1,mode,1,freq
nsort,u,sum
*get,max_disp_sum,sort,,max
*set,scale_factor,max_deffect/max_disp_sum
/prep7
f,1,fz,1.2*pressure1
upgeom,scale_factor,1,1,'pipeline-5','rst',''
finish
/solu
antype,0
nlgeom,on
nsubst,10000
arclen,on,2
arctrm,l
solve
finish
/post26
nsol,2,7,u,y,uy_2
prod,3,1,,,,,,1.2*pressure1,1,1
/axlab,x,Dis-x(m)
/axlab,y,Pres(Pa)
xvar,2
plvar,3
prvar,3
[ 本帖最后由 xinyuxf 于 2007-3-17 19:43 编辑 ] |