!----- Scrip for the model of Concealed bracing RC shear wall using ANSYS-APDL
! Unit:N-mm-MPa
! written by blackeagle,2006.1.20
! blog: http://wangjian.host.9466.cn
! Email: wang.jian@gerb.com.cn
/prep7
!-----------element--------------
et,1,beam188
et,2,solid65
et,3,solid45
KEYOPT,1,1,0
KEYOPT,1,3,2
!-----------section--------------
SECTYPE,1,BEAM,CSOLID,F4,0 !stirrups of concealed column
SECDATA,2
SECTYPE,2,BEAM,CSOLID,F6,0 !distributing reinforces,stirrups of beam
SECDATA,3
SECTYPE,3,BEAM,CSOLID,F10,0 !reinforces of column
SECDATA,5
SECTYPE,4,BEAM,CSOLID,F12,0 !reinforces of beam
SECDATA,6
!-----------material--------------
!--reinforce
EX,1,2E5
NUXY,1,0.3
TB,BISO,1,1,2,
TBDATA,,350,2e3
!--concrete
ex,2,25806.46
nuxy,2,0.2
TB,CONC,2,1
TBDATA,,0.125,0.9,2.7,25
ex,3,3e5
nuxy,3,0.2
!----------------model--------------------
!---stirrups of concealed column
n,1,-500,-45
ngen,5,1,1,1,,30
ngen,5,5,1,5,1,0,90/4
ngen,56,25,1,25,1,0,0,20
type,1 $ mat,1 $ secnum,1
*do,i,1,54,3
*do,j,1,4,1
e,j+i*25,j+1+i*25
e,j+20+i*25,j+21+i*25
*enddo
*do,j,1,4,1
e,j*5-4+i*25,(j+1)*5-4+i*25
e,j*5+i*25,(j+1)*5+i*25
*enddo
*enddo
!---reinforces of column
type,1 $ mat,1 $ secnum,3
*do,i,1,55,1
e,i*25+1-25,(i+1)*25+1-25
e,i*25+5-25,(i+1)*25+5-25
e,i*25+21-25,(i+1)*25+21-25
e,i*25+25-25,(i+1)*25+25-25
*enddo
!---distributing reinforces into column
type,1 $ mat,1 $ secnum,2
*do,i,1,39,4
e,i*25+11,i*25+15
*enddo
nsel,none
!---stirrups of beam
n,2001,-500,-80,800
ngen,11,1,2001,2001,1,0,0,30
ngen,2,11,2001,2011,1,0,35,
ngen,5,11,2012,2022,1,0,90/4
ngen,2,11,2056,2066,1,0,35
ngen,26,77,2001,2077,1,20,
*do,i,1,25,3
*do,j,1,10,1
e,j+i*77+2000-77,j+1+i*77+2000-77
e,j+i*77+66+2000-77,j+1+i*77+66+2000-77
*enddo
*do,j,1,6,1
e,j*11-10+i*77+2000-77,(j+1)*11-10+i*77+2000-77
e,j*11+i*77+2000-77,(j+1)*11+i*77+2000-77
*enddo
*enddo
!---reinforces of beam
secnum,4
*do,i,1,25,1
e,i*77-77+1+2000,(i+1)*77+1+2000-77
e,i*77-77+11+2000,(i+1)*77+11+2000-77
e,i*77-77+67+2000,(i+1)*77+67+2000-77
e,i*77+2000,(i+1)*77+2000
*enddo
!---distributing reinforces
nsel,none
n,5001,-380,-45
ngen,20,1,5001,,,20
ngen,5,20,5001,5020,1,0,90/4
ngen,56,100,5001,5100,1,0,0,20
secnum,2
*do,i,1,40,4
*do,j,1,19,1
e,i*100+5041+j-1,i*100+5041+j
*enddo
*enddo
*do,i,1,20,4
*do,j,1,55,1
e,(j-1)*100+5041+i,j*100+5041+i
*enddo
*enddo
!---concrete of column
block,-500,-380,-45,45,0,800
lsel,s,loc,z,0
lsel,a,loc,z,800
lesize,all,,,4
lsel,s,loc,z,1,800-1
lesize,all,,,40
type,2 $ mat,2
vmesh,1
!---concrete of beam
lsel,none
block,-500,-380,-80,-45,800,1100
lsel,s,length,,120
lesize,all,,,4
lsel,s,length,,35
lesize,all,,,1
lsel,s,length,,300
lesize,all,,,10
vmesh,2
asel,s,loc,x,-380
asel,r,loc,z,800,1100
num=arnext(0)
extopt,attr,0,0,0
extopt,esize,19,0,
voffst,num,380
vgen,2,2,3,1,0,125,0
!---corner of beam and column
asel,s,loc,z,800
asel,r,loc,y,-45,45
num=arnext(0)
extopt,esize,10
voffst,num,300
!---beam
asel,s,loc,z,800,1100
asel,r,loc,y,-45,45
asel,r,loc,x,-380
num=arnext(0)
extopt,esize,19
voffst,num,380
!---plate
asel,s,loc,z,0,800
asel,r,loc,x,-380
num=arnext(0)
voffst,num,380
allsel,all
nsym,x,20000,all
esym,,20000,all
nummrg,all, !节点合并
ndele,all
!--bottom
nsel,s,loc,z,800,890
esln,s,1
esel,u,type,,1
ngen,2,50000,all,,,,,-890
egen,2,50000,all,,,1,1
allsel,all
nummrg,all
ndele,all
!--concealed reinforces
nsel,s,loc,z,0,800
nsel,r,loc,x,-380,380
nsel,r,loc,y,0
n,100001,-380,-45,
ngen,39,1,100001,,,20,,20
n,200001,-380,-45,120
ngen,35,1,200001,,,20,,20
n,300001,380,-45,
ngen,39,1,300001,,,-20,,20
n,400001,380,-45,120
ngen,35,1,400001,,,-20,,20
nsel,u,loc,y,0
ngen,2,400000,all,,,,90
ksel,none
k,100001,-380,-45,
kgen,39,100001,,,20,,20,1
k,200001,-380,-45,120
kgen,35,200001,,,20,,20,1
k,300001,380,-45,
kgen,39,300001,,,-20,,20,1
k,400001,380,-45,120
kgen,35,400001,,,-20,,20,1
ksel,u,loc,y,0
kgen,2,all,,,,90,,400000
esel,none
lsel,none
type,1 $ mat,1 $ secnum,2
*do,i,1,38,1
e,100000+i,100001+i
e,300000+i,300001+i
e,500000+i,500001+i
e,700000+i,700001+i
*enddo
*do,i,1,34,1
e,200000+i,200001+i
e,400000+i,400001+i
e,600000+i,600001+i
e,800000+i,800001+i
*enddo
secnum,1
*do,i,6,36,2
l,100000+i,200000+i-3
l,300000+i,400000+i-3
l,500000+i,600000+i-3
l,700000+i,800000+i-3
*enddo
*do,i,6,36,2
l,100000+i,500000+i
l,200000+i-3,600000+i-3
l,300000+i,700000+i
l,400000+i-3,800000+i-3
*enddo
lesize,all,,,4
lmesh,all
allsel,all
nummrg,all
ndele,all
numcmp,all
!--------------------loading-------------------------
/solu
!---loading parameter----step 1
Fh=350e3
Fv=270e3
!---vertical load
nsel,s,loc,z,-90
d,all,all
allsel,all
nsel,s,loc,z,1100
sf,all,pres,Fv/(1000*160)
allsel,all
!solution control
NLGEOM,on
NSUBST,2,0,0
OUTRES,ALL,ALL
AUTOTS,1
NEQIT,100
TIME,1
solve
!---horizonal load---step 2
nsel,s,loc,z,800,1100
nsel,r,loc,x,-500
sf,all,pres,Fh/(160*300)
allsel,all
CNVTOL,,-1,1 !delete force criterion
CNVTOL,u,,0.05,0
NSUBST,150,0,0
OUTRES,ALL,ALL
AUTOTS,1
NEQIT,100
TIME,350
solve |