$lsel,u,loc,x, 0.161,2.25 $lsel,u,loc,y,0.125-0.015 $lsel,u,loc,y,-0.125+0.015 $lsel,u,loc,y,-0.125-0.06 $lsel,u,loc,y,0.125+0.06 $lsel,u,loc,y,-0.095,0.095 $lsel,u,loc,y,0.251,1.5 $lsel,u,loc,y,-1.5,-0.251 $lsel,u,loc,y, 0.125 $lsel,u,loc,y,-0.125 $Lsel,r,loc,x,0.16 $Lsel,r,loc,z,0,0.16
$cm,ssgj4,line$
$cmsel,s,ssgj1,line$$cmsel,a,ssgj2,line$
cmsel,a,ssgj3,line$cmsel,a,ssgj4,line$$cm,guj3,line$
$cmsel,s,guj1,line$$cmsel,a,guj3,line$$cm,guj,line$latt,5,4,3,0$lesize,all,elemsiz
$cmsel,s,guj,line$cmsel,a,zj,line$cm,gj,line
Cmsel,s,gzg,volu$VSEL,r,LOC,Y,-0.125+0.015,0.125-0.015 $VSEL,r,LOC,z,0,0.003 $cm,fb,VOLU
Cmsel,s,gzg,volu$VSEL,u,LOC,Y,-0.125,0.125$cm,GZhu,VOLU Cmsel,s,gzg,volu$$cmSEL,U,fb,VOLU$VSEL,r,LOC,Y,-0.125+0.015,0.125-0.015$cm,mianban,VOLU
Cmsel,s,gzg,volu$$cmSEL,U,fb,VOLU$$cmSEL,U,gzhu,VOLU$$cmSEL,U,mianban,VOLU$$cm,yy,VOLU
!!!!!!!!!!!!!!!!!!!!!!!! Mesh钢筋!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! cmsel,s,gj,line$lmesh,all $lsel,all $/eshape,1$eplot !钢筋划分完毕
!!!!!!!!!!!!!!!!!!!! mesh 钢筋混凝土柱!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Cmsel,s,zhuzi,volu $esize,50e-3$vatt,2,1,2,0 $MSHAPE,0,3D$MSHKEY,1$VMESH,all allsel,all
!!!!!!!!!!!!!!!!!! Mesh 翼缘!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Cmsel,s,yy,volu$ $LSEL,s,length,,1.8 LESIZE,all, , ,12
vatt,1,1,1,0$esize,50e-3$$MSHAPE,0,3D$MSHKEY,1 $VMESH,all$allsel,all
!!!!!!!!!!!!!!!!!! Mesh fb!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Cmsel,s,fb,volu$
vatt,1,1,1,0$esize,50e-3$$MSHAPE,0,3D$MSHKEY,1 $VMESH,all$allsel,all
!!!!!!!!!!!!!!!!!! Mesh mianban gzhu!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Cmsel,s,mianban,volu$Cmsel,a,gzhu,volu$
vatt,3,5,1,0$esize,50e-3$$MSHAPE,0,3D$MSHKEY,1 $VMESH,all$allsel,all
$nummrg,node$numcmp,all !定义远距离加载点 $ET,4,TARGE170
ET,5,CONTA174
R,6
keyopt,4,2,1$keyopt,5,2,2$keyopt,5,4,2$keyopt,5,12,5 asel,s,loc,x,2$nsla,s,1$real,6$type,5$esurf nmax=ndinqr(0,14)+1$n,nmax,2,, type,4$tshap,pilo$e,nmax R,7
keyopt,4,2,1$keyopt,5,2,2$keyopt,5,4,2$keyopt,5,12,5 asel,s,loc,x,-2$nsla,s,1$real,7$type,5$esurf nmax=ndinqr(0,14)+1$n,nmax,-2,, type,4$tshap,pilo$e,nmax R,8
keyopt,4,2,1$keyopt,5,2,2$keyopt,5,4,2$keyopt,5,12,5 asel,s,loc,y,1.5$nsla,s,1$real,8$type,5$esurf nmax=ndinqr(0,14)+1$n,nmax,0,1.5, type,4$tshap,pilo$e,nmax Allsel,all
!施加约束和荷载 Allsel,all Asel,s,loc,z,0 Da,all,symm $D,8665,uy,,,,,uz $D,8666,uy,,,,,uz $
Asel,s,loc,y,-1.5 Da,all,all Allsel,all
/solu$ANTYPE,0$OUTRES,ALL,ALL$$pred,off
$neqit,30$lnsrch,on$nlgenom,off$Cnvtol,f,,0.05$kbc,0$eqslv,sparse,1.0E-5,2 $p=0.01
$time,1$nsubst,50,,10$Asel,s,loc,y,1.5
f,8667,fy,-582000$$allsel,all$OUTRES,ALL,ALL $ALLSEL,ALL $LSWRITE,1
TIME,2$D,8667,Ux,19*p$$kbc,0$ NSUBST,100,1000000,1$AUTOTS,on$OUTRES,ALL,ALL $ALLSEL,ALL $LSWRITE,2 $Lssolve,1,2