计算拱肋稳定性的
finish
/cle
! 不加耳朵
/prep7
et,1,solid45
mp,ex,1,2.1e11
mp,prxy,1,0.28
mp, dens, 1, 7800 ! 材料
tb,biso,1
tbdata,,235e6
BLC4,0,0,0.022,0.12,0.3 ! 建模
vsweep,all
nsel,,loc,y,0
d,all,all
allsel
save,a1-model,db
finish
! 线性稳定
/solu
antype,static
OUTPR,all,all
OUTRES,all,all
! NLGEOM,on ! 不考虑几何非线性
! sstif,on
pstres,on
acel,,9.8
asel,,,,4
SFA,all,1,PRES,1e8 ! 基准荷载q=1e8 N/(m*m) P=0.022*0.3*q=66 t
allsel
solve !静力求解
finish
!稳定求解
/solu
antype,buckle
outpr,nsol,all
bucopt,lanb,5
mxpand,5
solve
finish
/post1
SET,LIST
!! 非线性稳定
/eof
/solu
antype,static
OUTPR,all,all
OUTRES,all,all
NLGEOM,on ! 考虑几何非线性
sstif,on
AUTOTS,0
NSUBST,1
acel,,9.8
allsel
lswrite,1 ! 自重
NSUBST,500,0,0
asel,,,,4
SFA,all,1,PRES,5e9 ! 1e8*50
nsel,,loc,x,0
nsel,r,loc,y,0.12
f,all,fx,10000 ! 加初始缺陷 干挠力
allsel
lswrite,2
lssolve,1,2,1
finish |