非线性屈曲的分析，后面那部分不太会，高手帮忙指教一下

分享蓝天

非线性屈曲的分析，后面那部分不太会，高手帮忙指教一下

finish
/clear
/filname,test02
!设置参数
a=1.5
b=1
cd=10
rad=cd/57.2956
r=0.1
c=b*cos(rad)-(a-b*sin(rad)-r+r*sin(rad))*tan(rad)-r*cos(rad)
!设置截面参数
b1=0.05
h1=0.02
a1=b1*h1
iy1=b1*b1*b1*h1/12
iz1=h1*h1*h1*b1/12
!前处理
/prep7
et,1,beam3
r,1,a1,iz1,h1
mp,ex,1,2e11
mp,nuxy,1,0.3
!建模
k,1000,0,0
k,1001,b*sin(rad),b*cos(rad)
circle,1000,,,1001,cd
k,1003,a-r,c
k,1004,a,c
circle,1003,,,1004,90-cd
l,4,1
k,5,a,0
l,5,3
!网格划分
lesize,1,,,b*rad*10
lesize,2,,,r1*(3.1416-rad)*10*2
lesize,3,,,a*10
lesize,4,,,c*10
lmesh,all
/eshape
/replot
save
!静力解
/solu
antype,static
pstres,on
nsel,s,,,2
dsym,symm,x  !在节点2施加x对称载荷
nsel,s,,,24
dsym,symm,y  !在节点24施加y对称载荷
sfbeam,all,,pres,-1e5
solve
finish
!特征值屈曲解
/solu
antype,buckle
bucopt,lanb,1
mxpand,1
outpr,all,all
solve
finish
!非线性分析设置
/config,nres,1000  !设置分析终止条件（该处不理解）
/prep7
tb,biso,1,1,2  !修正材料模型
tbtemp,0
tbdata,,2e8,0
upgeom,0.01,1,1,'test02','rst'  !引入几何缺陷
finish
!非线性分析
/solu
antype,static
nlgeom,1  !打开大变形
outres,all,all  !设置输出
arclen,1,0,0   !打开弧长法
arctrm,u,0.25,2,uy  !设置弧长法选项(该处不会设置)
nsubst,200,,,1   !设置子步数(该处不是很明确)
solve
finish
!后处理
/post26
nsol,2,2,u,y,uy0  !指定位移变量（该处不会设置，不知道如何选择位移变量)
esol,3,2,2,f,x,fx0  !指定力分量（同上，不会选择）
/axlab,y,fx0
xvar,2
plvar,3
finish