钢筋混凝土简支梁

!------------------------------------------------------
!EX8.26 钢筋混凝土简支梁数值分析
!分离式模型,关闭压碎,keyopt(1)=0,keyopt(7)=1
!力加载,位移收敛准则，误差1.5%,1/2模型分析
!--------------------------------------------
finish
/clear
/config,nres,2000
/prep7
!1.定义单元与材料性质--------------------
et,1,solid65,,,,,,,1 !K1=0,k7=1
et,2,link8
et,3,link8
mp,ex,1,13585
mp,prxy,1,0.2
fc=14.3
ft=1.43
tb,concr,1
tbdata,,0.5,0.95,ft,-1 !定义混凝土材料及相关参数，关闭压碎tb,miso,1,,10 !MIso模型
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.36
tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64
tbpt,,0.001,fc*0.75
tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91
tbpt,,0.0016,fc*0.96
tbpt,,0.0018,fc*0.99
tbpt,,0.002,fc
mp,ex,2,2.0e5
mp,prxy,2,0.3
tb,biso,2
tbdata,,360,0
pi=acos(-1)
r,1,0.25*pi*16*16
r,2,0.25*pi*8*8
mp,ex,3,2.1e5
mp,prxy,3,0.3
tb,biso,3
tbdata,,270,0
pi=acos(-1)
r,3,0.25*pi*8*8
r,4,0.25*pi*8*8/2
!2.创建几何模型
blc4,,,150,300,1400 wpoff,,,1400
!箍筋位置
*do,i,1,8
wpoff,,,-100 !平移工作平面vsbw,all
*enddo
!加载面位置
wpcsys,-1
!wpoff,,,-600
wpoff,,,550
vsbw,all
wpoff,,,100
vsbw,all
!支撑面
wpcsys,-1
wpoff,,,1150
vsbw,all
wpoff,,,100
vsbw,all
!钢筋位置
wpcsys,-1
wprota,,-90
wpoff,,,30
vsbw,all
wpoff,,,240
vsbw,all
wpcsys,-1
wpoff,30
wprota,,,90
vsbw,all
wpoff,,,90
vsbw,all
wpcsys,-1
!划分钢筋网格
elemsiz=50
!受拉钢筋
lsel,s,loc,x,30
lsel,r,loc,y,30
cm,lj1,line
latt,2,1,2
lesize,all,elemsiz
lsel,s,loc,x,120
lsel,r,loc,y,30
cm,lj2,line
latt,2,1,2
lesize,all,elemsiz
!受压钢筋
lsel,s,loc,x,30
lsel,r,loc,y,270
lsel,r,loc,z,600,1400
cm,yj1,line
latt,2,2,2
lesize,all,elemsiz
lsel,s,loc,x,120
lsel,r,loc,y,270
lsel,r,loc,z,600,1400
cm,yj2,line
latt,2,2,2
lesize,all,elemsiz
!箍筋
lsel,s,tan1,z !选取所有与Z垂直的线lsel,r,loc,y,30,270
lsel,r,loc,x,30,120
lsel,u,loc,z,0 !删除位置
lsel,u,loc,z,550
lsel,u,loc,z,650
lsel,u,loc,z,1400
lsel,u,loc,z,1150
lsel,u,loc,z,1250
cm,gj,line
latt,2,3,2
lesize,all,elemsiz
!lsel,s,loc,z,1350
!lsel,r,loc,y,30,270
!lsel,r,loc,x,30,120
!cm,gjb,line
!latt,2,4,2
!lesize,all,elemsiz
lsel,all
cmsel,s,lj1
cmsel,a,lj2
cmsel,a,yj1
cmsel,a,yj2
cmsel,a,gj
!cmsel,a,gjb
cm,gj,line
lmesh,all
lsel,all
!划分混凝土网格
vatt,1,,1 !设置体单元属性mshkey,1
esize,elemsiz
vmesh,all
allsel,all
numcmp,all
nummrg,all
!5.施加荷载和约束
!lsel,s,loc,y,0
!lsel,r,loc,z,1200
!dl,all,,uy
!lsel,s,loc,y,0
!lsel,r,loc,z,1200
!dl,all,,ux
!lsel,s,loc,y,0
!lsel,r,loc,z,1200
!dl,all,,uz
asel,s,loc,y,0
asel,r,loc,z,1150,1250 da,all,ux
asel,s,loc,y,0
asel,r,loc,z,1150,1250 da,all,uy
asel,s,loc,z,0
da,all,symm
!-----------------
p0=100000
q0=p0/150/100
asel,s,loc,z,550,650 asel,r,loc,y,300
sfa,all,1,pres,q0
!dl,all,,uy,-15
allsel,all
!6.求解控制设置
!/solu
!antype,0
!nsubst,100
!outres,all,all
!autos,on
!neqit,50
!cnvtol,u,,0.03
!solve
/solu
antype,0
nsubst,80 !荷载子步80
!NSUBST,NSBSTP,NSBMX,NSBMN,Carry
outres,all,all
autos,on !打开自动时间步
!如果使用了自动时间步（即AUTOTS,ON）则该子步数仅用于第一子步，也即第一子步的荷载增量用NSBSTP求得，其余子步的荷载增量由程序自动确定。

neqit,50
!这里采用默认的收敛准则
!cnvtol,f,,0.015 !位移收敛准则
solve
!8.进入时程后处理
!荷载-挠度曲线
/post26
nsol,2,997,u,y
prod,3,1,,,P,,,p0/1000
prod,4,2,,,f,,,-1
xvar,4
plvar,3
/XRANGE, 0,10
plvar,3
/post26
nsol,2,296,u,y
!prod,3,1,,,P,,,p0/1000
!prod,4,2,,,f,,,-1
xvar,1
plvar,2
/XRANGE, 0,1 plvar,2
!钢筋应力
/post26
esol,6,28,35,s,z xvar,1
plvar,6
/XRANGE, 0,1 plvar,6。