!EX8.6 受彎矩作用懸臂梁幾何非線性分析
finish
/clear
/prep7
cdl=12
gdh=1
kdb=1
txml=30e6
mja=gdh*kdb
gxji=kdb*gdh*gdh*gdh/12.0
hzm=acos(-1)*txml*kdb*gdh*gdh*gdh/12.0/cdl
et,1,beam3
mp,ex,1,txml
r,1,mja,gxji,gdh
k,1
k,2,cdl
l,1,2
lesize,all,,,10
lmesh,all
finish       
/solu
dk,1,all
antype,0
nlgeom,1
nsubst,100
outres,all,all
*do,i,1,40
fk,2,mz,i/10.0*hzm
lswr,i
*enddo
lssolve,1,40 
/post1
pldisp,0
/post26
nsol,2,2,u,y
XVAR,2
PLVAR,1