fini
/cle
/config,nres,20000
/prep7
/title,buckling of a bar with hinged solves
et,1,beam3
r,1,0.25,52083e-7,.5
mp,ex,1,30e6
mp,prxy,1,0.3
n,1
n,11,,100
fill
e,1,2
egen,10,1,1 !复制单元，10为复制次数，包括本身；1节点增量；1复制元素
finish
!静力求解
/solu
antype,static
pstres,on !打开预应力效应
d,1,all
f,11,fy,-1
solve
finish
!特征值屈曲分析（线性）
/solu
antype,buckle !屈曲分析
bucopt,lanb,1 !使用莱布尼兹法扩展模态，提取特征值数为1阶
mxpand,1 !扩展一阶模态
solve
finish
!施加初始扰动
/post1
set,first
pldisp,1
upcoord,0.1,on !失稳模态变形乘以0.1，更新节点坐标系
save,pert,db !保存有初始扰动的模型
finish
resu,pert,db !调入有初始扰动的模型
!非线性屈曲分析
/solu
antype,static
nlgeom,on !打开大变形
kbc,0 !斜坡载荷
nsubst,100 !定义子步数
arclen,on,2 !采用弧长法，半径最大为2
outres,,1 !控制输出所有子步结果
f,11,fy,-40 !施加载荷40，线性屈曲分析得到预测载荷38.553
solve
fini
/post1
pldisp,1 !显示屈曲后的变形
/post26
nsol,2,11,u,y,U1 !11节点的y位移作为2号变量
rforce,3,1,f,y,FY !1节点的y方向的支座反力作为3号变量
prod,4,2,,,U2,,,-1,1,1 !2号变量乘-1作为4号变量
/AXLAB,X,U !x轴用U标注，默认为TIME
/AXLAB,Y,F(N) !x轴用F(N)标注，默认为VALUE
/xrange,4 !指定x轴的范围
xvar,4 !以4号变量为x轴
plvar,3 !以3号变量为y轴，显示载荷-位移曲线
fini