圆管受集中力作用下的非线性屈曲(第二类稳定)

本文链接：https://blog.youkuaiyun.com/yueliang2100/article/details/82190045

本文介绍了一个使用有限元分析软件进行的结构变形与载荷关系的研究案例。通过定义材料属性、几何参数及载荷条件，模拟了不同温度下结构的位移响应，并分析了总载荷与位移之间的变化趋势。

r1=2
l=2
pi=acos(-1)
theta=(0.5*pi)*180/pi
/config,nres,2000
/prep7

mp,ex,1,2.160687961e11
mp,nuxy,1,0.3
tb,MISO,1,1,23
tbtemp,0
tbPT,,0.000814,175.88e6
tbPT,,0.001531,324.98e6
tbPT,,0.002819,591.26e6
tbPT,,0.003846,801.457e6
tbPT,,0.004239,873.22e6
tbPT,,0.004594,924.399e6
tbPT,,0.004955,956.88e6
tbPT,,0.00531,976.99e6
tbPT,,0.005699,990.2138e6
tbPT,,0.006261,1000e6
tbPT,,0.006957,1008.34e6
tbPT,,0.007804,1015.5e6
tbPT,,0.008779,1021.185e6
tbPT,,0.009827,1026.47e6
tbPT,,0.010929,1031.289e6
tbPT,,0.012066,1035.444e6
tbPT,,0.0132,1039.693e6
tbPT,,0.014354,1044.886e6
tbPT,,0.015494,1049.796e6
tbPT,,0.016631,1053.951e6
tbPT,,0.017785,1058.106e6
tbPT,,0.019039,1062.827e6
tbPT,,0.020276,1067.643e6

et,1,93
keyopt,1,6,1
r,1,0.02
csys,1
n,1,r1,90
n,2,r1,90,l
k,1,r1,90
k,2,r1,90-theta
k,3,r1,90,l
k,4,r1,90-theta,l

a,1,3,4,2
esize,0,8
mshkey,1
amesh,all
nummrg,node,,,,low
/solu
nsel,s,loc,z,0
dsym,symm,z
nsel,s,loc,y,90
dsym,symm,x

nsel,s,loc,y,90-theta
d,all,,,,,,ux,uy,uz
alls
/view,1,1,1,1
antype,0
nlgeom,1
f,1,fy,-450e3
nsubst,100,,,1
outres,all,all
arclen,1,4,0
solve

/post26
nsol,2,1,u,y,n1_disp
nsol,3,2,u,y,n2_disp

prod,4,1,,,load,,,4*450e3,1,1
prod,5,2,,,,,,-1,1,1
prod,6,3,,,,,,-1,1,1

*get,uy1,vari,2,extrem,vmin
*get,uy2,vari,3,extrem,vmin
/axlab,x,delfection(mm)
/axlab,y,total load(N)
/xrange,0,5
/yrange,-2.1*450e3,4*450e3,1
/grid,1
xvar,5
plvar,4
xvar,6
plvar,4