是一個網(wǎng)友做過的例子，和大家分享！
!題目：鋼筋混凝土簡支梁模擬計算
!方法：分離式；solid65和link8
!材料：混凝土采用concr和鋼筋為彈性材料，但不考慮壓碎
!---------------------------------------------------
!為方便，假定鋼筋置于梁底兩側(cè).
!===================================================
/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
!定義網(wǎng)分時邊長控制
lsel,s,loc,z,1,2999
lsel,r,loc,y,0
latt,2,2,2
lesize,all,,,20 !鋼筋網(wǎng)格數(shù)目
lmesh,all
lsel,s,loc,z,0
lesize,all,,,4 !截面上的網(wǎng)格數(shù)目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
