钢筋混凝土简支梁.docx

!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=1et,2,link8et,3,link8mp,ex,1,13585mp,prxy,1,0.2fc=14.3ft=1.43tb,concr,1tbdata,,0.5,0.95,ft,-1 !定义混凝土材料及相关参数，关闭压碎tb,miso,1,,10 !MIso 模型tbpt,,0.0002,fc*0.19tbpt,,0.0004,fc*0.36tbpt,,0.0006,fc*0.51tbpt,,0.0008,fc*0.64tbpt,,0.001,fc*0.75tbpt,,0.0012,fc*0.84tbpt,,0.0014,fc*0.91tbpt,,0.0016,fc*0.96tbpt,,0.0018,fc*0.99tbpt,,0.002,fcmp,ex,2,2.0e5mp,prxy,2,0.3tb,biso,2tbdata,,360,0pi=acos(-1)r,1,0.25*pi*16*16r,2,0.25*pi*8*8mp,ex,3,2.1e5mp,prxy,3,0.3tb,biso,3tbdata,,270,0pi=acos(-1) r,3,0.25*pi*8*8 r,4,0.25*pi*8*8/2 !2.创建几何模型blc4,,,150,300,1400wpoff,,,1400!箍筋位置*do,i,1,8wpoff,,,-100 !平移工作平面vsbw,all*enddo!加载面位置wpcsys,-1!wpoff,,,-600wpoff,,,550vsbw,allwpoff,,,100vsbw,all!支撑面wpcsys,-1wpoff,,,1150vsbw,allwpoff,,,100vsbw,all!钢筋位置wpcsys,-1wprota,,-90wpoff,,,30vsbw,allwpoff,,,240vsbw,allwpcsys,-1wpoff,30wprota,,,90vsbw,allwpoff,,,90vsbw,allwpcsys,-1!划分钢筋网格elemsiz=50!受拉钢筋lsel,s,loc,x,30lsel,r,loc,y,30cm,lj1,linelatt,2,1,2lesize,all,elemsizlsel,s,loc,x,120lsel,r,loc,y,30cm,lj2,linelatt,2,1,2lesize,all,elemsiz!受压钢筋lsel,s,loc,x,30lsel,r,loc,y,270lsel,r,loc,z,600,1400cm,yj1,linelatt,2,2,2lesize,all,elemsizlsel,s,loc,x,120lsel,r,loc,y,270lsel,r,loc,z,600,1400cm,yj2,linelatt,2,2,2lesize,all,elemsiz顷筋lsel,s,tan1,z !选取所有与Z垂直的线lsel,r,loc,y,30,270lsel,r,loc,x,30,120lsel,u,loc,z,0 !删除位置lsel,u,loc,z,550lsel,u,loc,z,650lsel,u,loc,z,1400lsel,u,loc,z,1150lsel,u,loc,z,1250cm,gj,linelatt,2,3,2lesize,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,elemsizlsel,all cmsel,s,lj1cmsel,a,lj2cmsel,a,yj1cmsel,a,yj2cmsel,a,gj!cmsel,a,gjbcm,gj,linelmesh,alllsel,all!划分混凝土网格vatt,1,,1 !设置体单元属性mshkey,1esize,elemsizvmesh,allallsel,allnumcmp,allnummrg,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,,uzasel,s,loc,y,0asel,r,loc,z,1150,1250da,all,uxasel,s,loc,y,0asel,r,loc,z,1150,1250da,all,uyasel,s,loc,z,0da,all,symm! p0=100000q0=p0/150/100asel,s,loc,z,550,650asel,r,loc,y,300sfa,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/soluantype,0nsubst,80 !荷载子步80!NSUBST,NSBSTP,NSBMX,NSBMN,Carryoutres,all,allautos,on !打开自动时间步!如果使用了自动时间步（即AUTOTS,ON ）则该子步数仅用于第一子步，也即第一子步的荷 载增量用NSBSTP求得，其余子步的荷载增量由程序自动确定。

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