隨手做了個命令流，對平面杆進行特徵值屈曲和非線性屈曲
水平有限，算出來的非線性屈曲居然比特徵值屈曲大，不合理啊！
命令流如下：
FINISH
/CLEAR
/PREP7
ET,1,beam3
r,1,1,1/12,1
mp,ex,1,2.1e5
mp,prxy,1,0.3
n,1,
n,2,0,10
e,1,2
/solu
antype,static
alls
pstres,on
d,1,all
F,2,fy,-1000
ALLS
solve
fini
antype,buckle
bucopt,lanb,4
OUTPR,all,ALL
finish
/SOLU !!對解進行擴充套件
EXPASS,ON
MXPAND,,,,YES
OUTRES
SOLVE
!以下進行非線性屈曲分析
*get,freq1,mode,1,freq
force=-1000*freq1
f,2,fy,force
/post1
set,first
nsort,u,sum
*get,dmax,sort,,max
/prep7
!對模型加缺陷，file改為你的定義的工作檔名
upgeom,0.1*dmax,1,1,file,rst
antype,0
force1=2*force
f,2,fy,force1
nlgeom,on
nsubst,20,0,0
arclen,1,0,0
outres,basic,1
!繪出節點2載荷位移圖
/post26
nsol,2,2,u,y,N2_Y
prod,3,2,,,N2_Y1,,,-1
prod,4,1,,,LOAD,,,-2*force,
/AXLAB,X,LOAD(N)
/AXLAB,Y,U2_Y(MM)
XVAR,4
PLVAR,3
隨手做了個命令流，對平面杆進行特徵值屈曲和非線性屈曲
水平有限，算出來的非線性屈曲居然比特徵值屈曲大，不合理啊！
命令流如下：
FINISH
/CLEAR
/PREP7
ET,1,beam3
r,1,1,1/12,1
mp,ex,1,2.1e5
mp,prxy,1,0.3
n,1,
n,2,0,10
e,1,2
/solu
antype,static
alls
pstres,on
d,1,all
F,2,fy,-1000
ALLS
solve
fini
/solu
antype,buckle
bucopt,lanb,4
OUTPR,all,ALL
solve
finish
/SOLU !!對解進行擴充套件
EXPASS,ON
MXPAND,,,,YES
OUTPR,all,ALL
OUTRES
SOLVE
FINISH
!以下進行非線性屈曲分析
/solu
alls
*get,freq1,mode,1,freq
force=-1000*freq1
f,2,fy,force
solve
finish
/post1
set,first
nsort,u,sum
*get,dmax,sort,,max
finish
/prep7
!對模型加缺陷，file改為你的定義的工作檔名
upgeom,0.1*dmax,1,1,file,rst
finish
/solu
antype,0
force1=2*force
f,2,fy,force1
nlgeom,on
nsubst,20,0,0
arclen,1,0,0
outres,basic,1
alls
solve
finish
!繪出節點2載荷位移圖
/post26
nsol,2,2,u,y,N2_Y
prod,3,2,,,N2_Y1,,,-1
prod,4,1,,,LOAD,,,-2*force,
/AXLAB,X,LOAD(N)
/AXLAB,Y,U2_Y(MM)
XVAR,4
PLVAR,3