! 书中例题
FINISH$/CLEAR$/CONFIG,NRES,5000$/PREP7 et,2,link8 et,1,solid65 keyopt,1,1,1 keyopt,1,3,2 keyopt,1,7,1 mp,ex,1,19095 mp,prxy,1,0.2 fck=20.1 ftk=2.01 tb,concr,1
tbdata,,0.5,0.95,ftk,-1 tb,miso,1,,11
tbpt,,0.0002,fck*0.19 tbpt,,0.0004,fck*0.36 tbpt,,0.0006,fck*0.51 tbpt,,0.0008,fck*0.64 tbpt,,0.0010,fck*0.75 tbpt,,0.0012,fck*0.84 tbpt,,0.0014,fck*0.91 tbpt,,0.0016,fck*0.96 tbpt,,0.0018,fck*0.99 tbpt,,0.0020,fck*1.00 tbpt,,0.0033,fck*0.85 !定义钢筋材料参数 mp,ex,2,2e5 mp,prxy,2,0.3 tb,biso,2 tbdata,,335,0 mp,ex,3,2e5
mp,prxy,3,0.3 tb,biso,3 tbdata,,235,0 pi=acos(-1) !π
r,1,0.25*pi*22*22 !22直径圆 r,2,0.25*pi*22*22/2 r,3,0.25*pi*10*10 r,4,0.25*pi*10*10/2 r,5
!建立几何模型并切分 blc4,,,150/2,300,2000/2
*do,i,1,2$wpoff,,,100$vsbw,all$*enddo !沿z向移动两个100并切分 *do,i,1,4$wpoff,,,50$vsbw,all$*enddo
*do,i,1,5$wpoff,,,100$vsbw,all$*enddo
wpcsys,-1$wprota,,-90$wpoff,,,30 !坐标轴回并旋转、移动、切分 vsbw,all
wpoff,,,240 vsbw,all wpcsys,-1 wpoff,30 wprota,,,90 vsbw,all wpcsys,-1 elemsiz=50
!附钢筋单元属性并划分网格
lsel,s,loc,x,30$lsel,r,loc,y,30$cm,zj,line$latt,2,1,2$lesize,all,elemsiz !选中钢筋位置线并定义名称，附钢筋材质
lsel,s,loc,x,75$lsel,r,loc,y,30$cm,zjb,line$latt,2,2,2$lesize,all,elemsiz lsel,s,loc,x,30$lsel,r,loc,y,270$cm,jlj,line$latt,3,3,2$lesize,all,elemsiz lsel,s,tan1,z$lsel,r,loc,y,30,270$lsel,r,loc,x,30,70$lsel,u,loc,z,250 lsel,u,loc,z,350$cm,gj,line$latt,3,3,2$lesize,all,elemsiz lsel,s,loc,z,0$lsel,r,loc,y,30,270$lsel,r,loc,x,30,70 cm,gjb,line$latt,3,4,2$lesize,all,elemsiz lsel,all
cmsel,s,zj$cmsel,a,zjb$cmsel,a,jlj$cmsel,a,gj cmsel,a,gjb$cm,gj,line lmesh,all lsel,all
!划分混凝土单元网格 vatt,1,5,1 mshkey,1 esize,elemsiz vmesh,all allsel,all !施加约束
lsel,s,loc,y,0$lsel,r,loc,z,900 dl,all,,uy asel,s,loc,z,0 da,all,symm asel,s,loc,x,75 da,all,symm p0=180e3
q0=p0/150/100 asel,s,loc,z,250,350 asel,r,loc,y,300 sfa,all,1,pres,q0 allsel
!求解控制设置 /solu antype,0 outres,all,all autots,on neqit,40
nsubst,400,,200 cnvtol,u,,0.05 cnvtol,f,,0.05 solve
!后处理 /post26
node1=node(75,0,0) nsol,2,node1,u,y prod,3,1,,,,,,p0/1000 pprod,4,2,,,,,,-1 xvar,4 plvar,3
某钢筋混凝土分离式模型计算命令流
FINISH$/CLEAR$/CONFIG,NRES,5000$/PREP7 et,2,link8 et,1,solid65 keyopt,1,1,1 keyopt,1,3,2 keyopt,1,7,1 mp,ex,1,11305 mp,prxy,1,0.2 fck=11.9 ftk=1.27 tb,concr,1
tbdata,,0.5,0.95,ftk,-1 tb,miso,1,,11
tbpt,,0.0002,fck*0.19 tbpt,,0.0004,fck*0.36 tbpt,,0.0006,fck*0.51 tbpt,,0.0008,fck*0.64 tbpt,,0.0010,fck*0.75 tbpt,,0.0012,fck*0.84
tbpt,,0.0014,fck*0.91 tbpt,,0.0016,fck*0.96 tbpt,,0.0018,fck*0.99 tbpt,,0.0020,fck*1.00 tbpt,,0.0033,fck*0.85 !定义钢筋材料参数 mp,ex,2,2e5 mp,prxy,2,0.3 tb,biso,2
tbdata,,335,0 pi=acos(-1) !π
r,1,0.25*pi*8*8 !8直径圆 r,2,0.25*pi*8*8/2 r,3
!建立几何模型并切分 blc4,,,7800,160,7200
*do,i,1,17$wpoff,,,400$vsbw,all$*enddo !沿z向移动200并切分 wprota,,-90
wpoff,,,40$vsbw,all wpoff,,,80$vsbw,all wpcsys,-1$wprota,,-90 wprota,,,-90
*do,i,1,19$wpoff,,,-7800/20$vsbw,all$*enddo elemsiz=100 wpcsys,-1
!附钢筋单元属性并划分网格
lsel,s,loc,y,40$cm,bj,line$latt,2,1,2$lesize,all,elemsiz !选中钢筋位置线并定义名称，附钢筋材质 lsel,s,loc,y,120$cm,dj,line$latt,2,1,2$lesize,all,elemsiz cmsel,s,bj
cmsel,a,dj$cm,gj,line lmesh,all
!划分混凝土单元网格 lsel,all vatt,1,3,1 mshkey,1 esize,elemsiz vmesh,all allsel,all !施加约束 lsel,s,loc,x,0 dl,all,,all
lsel,s,loc,x,7800 dl,all,,all lsel,s,loc,z,0
dl,all,,all lsel,s,loc,z,7200 dl,all,,all allsel,all !加载 q0=2e-3
asel,s,loc,y,160 sfa,all,1,pres,q0 allsel
!求解控制设置 /solu antype,0 outres,all,all autots,on neqit,40
nsubst,400,,200 cnvtol,u,,0.05 cnvtol,f,,0.05 solve !后处理
!将背景设为白色 jpgprf,500,100,1 /rep
/post26
node1=node(75,0,0) nsol,2,node1,u,y prod,3,1,,,,,,p0/1000 pprod,4,2,,,,,,-1 xvar,4 plvar,3