【轉帖】鋼筋混凝土模擬計算
!題目：鋼筋混凝土簡支梁模擬計算 
!方法：分離式；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