是一个网友做过的例子，和大家分享！
!题目：钢筋混凝土简支梁模拟计算
!方法：分离式；solid65和link8
!材料：混凝土采用concr和钢筋为弹性材料，但不考虑压碎
!---------------------------------------------------
!为方便，假定钢筋置于梁底两侧.
!===================================================
/config,nres,2000
/prep7
!定义单元及其材料特性等
rd0=20.0 !钢筋直径
et,1,solid65
et,2,link8
mp,ex,1,33e3
mp,prxy,1,0.20
r,1
hntra=28
hntrl=2.6
tb,concr,1
tbdata,,0.7,1.0,hntrl,-1
mp,ex,2,2.1e5
mp,prxy,2,0.30
r,2,acos(-1)*0.25*rd0*rd0
!定义梁体即单元划分
blc4, , ,100,200,3000
/view,1,1,1,1
/ang,1
gplot
!定义网分时边长控制
lsel,s,loc,z,1,2999
lsel,r,loc,y,0
latt,2,2,2
lesize,all,,,20 !钢筋网格数目
lmesh,all
lsel,s,loc,z,0
lesize,all,,,4 !截面上的网格数目4x4
vsel,all
vatt,1,1,1
mshape,0,3d
mshkey,1
vmesh,all
allsel,all
finish
/solu
!施加约束
lsel,s,loc,z,0
lsel,r,loc,y,0
dl,all,,uy
dl,all,,uz
lsel,all
lsel,s,loc,z,3000
lsel,r,loc,y,0
dl,all,,uy
lsel,all
ksel,s,loc,x,0
ksel,r,loc,y,0
dk,all,ux
allsel,all
!施加荷载
qmz=0.3
asel,s,loc,y,200
sfa,all,1,pres,qmz
allsel,all
nsubst,40
outres,all,all
time,qmz*10
neqit,40
solve
finish
/post1
pldisp,1
etable,zxyl,ls,1
plls,zxyl,zxyl,1
finish
/post26
nsol,2,33,u,y
prod,3,1,,,,,,1/100
prod,4,2,,,,,,-1
xvar,4
plvar,3
