柱子稳定分析算

! 要点：预应力，特征值屈曲，添加初始缺陷，几何非线性分析

finish 
/CLEAR 
/UNITS,SI 
/PREP7 
!* 
FORCE=100 
OFFSET=0.1 !初始缺陷为 0.1 
! 建立模型 
ET,1,BEAM4 
ET,2,LINK10 
R,1,0.1*0.12,0.12*0.1**3/12,0.1*0.12**3/12,0.12,0.1, , 
R,2,0.002*0.002,2e-3, !预应力 
MPTEMP,1,0 
MPDATA,EX,1,,200e9 
MPDATA,PRXY,1,,0.27 
MPDATA,DENS,1,,7800 
k,1, 
k,2,0,0,5 
k,3,0,0,-5 
k,11,0.2,0,0 
k,12,-0.2,0,0 
k,13,0,0.2 
k,15,0,-0.2 
l,1,2 
l,1,3 
l,1,11 
l,1,12 
l,1,13 
l,1,15 270 
l,2,11 
l,2,12 
l,2,13 
l,2,15 
l,3,11 
l,3,12 
l,3,13 
l,3,15 
lsel,,,,1,6 
latt,1,1,1 
ALLSEL,ALL 
lsel,,,,7,14 
latt,1,2,2 
ALLSEL,ALL 
lsel,,,,1,6 
LESIZE,all,0.3, , , , , , ,1 
lsel,,,,7,14 
LESIZE,all, , ,1 , , , , ,1 
ALLSEL,ALL 
LMESH,ALL 
FINISH 
/ESHAPE,1.0 
!* 
! 求解特征值屈曲荷载 
*DO,I,1,100 
FINISH 
/SOLU 
DK,3, , , ,0,UX,UY,UZ, , ,ROTZ , 
DK,2, , , ,0,UX,UY, , , , , 
FK,2,FZ,-FORCE 
ANTYPE,0 
!设定时间步 
TIME,1 
A