/prep7
et,1,beam188           !梁的單元類型
et,2,shell63           !筒體,樓板及外墻單元類型
mp,ex,1,3.25e10        !C40混凝土
mp,prxy,1,0.2          !泊松比
mp,dens,1,2700         !密度
mp,ex,2,3.0e10         !C30混凝土
mp,prxy,2,0.2          !泊松比
mp,dens,2,2700         !密度

sectype,1,beam,rect    !框架柱截面
secdata,1.1,1.1

sectype,2,beam,rect    !外環梁截面
secdata,0.4,0.6

sectype,3,beam,rect    !內框架梁截面
secdata,0.5,0.8

sectype,4,beam,rect    !次梁
secdata,0.3,0.5


r,2,0.2                !樓板厚度
r,3,0.3                !筒體厚度

k,5000,22,8,72  !設置beam188單元節點的方向
k,1 
k,12,44  
kfill,1,12  
kgen,5,1,12,,,4   
kgen,19,1,60,,,,3 
/view,1,1,1,1
/ANGLE,1,270,XM,0   
/REPLOT

!框架柱共288根
*do,i,1,1021,60
l,i,i+60
*enddo
lgen,3,1,18,1,8,,,2
lgen,2,1,54,1,,16,,48
lgen,2,1,18,1,,8,,24
lgen,2,37,54,1,,8,,24
lgen,2,1,144,1,28,,,7

latt,1,,1,,5000,,1
lesize,all,1.5
lmesh,all
lsel,u,,,all
!線單元的最大編號為288

!外環梁KL2共144根
l,61,65
l,65,89
l,89,92
l,92,68
l,68,72
l,72,120
l,120,109
l,109,61
lgen,18,289,296,1,,,3,60

latt,1,,1,,5000,,2
lesize,all,2
lmesh,all
lsel,u,,,all
!線單元的最大編號為432


!內框架梁KL1共108根
l,63,111
l,85,89
l,89,113
l,70,118
l,92,96
l,92,116
lgen,18,433,438,1,,,3,60

latt,1,,1,,5000,,3
lesize,all,2
lmesh,all
lsel,u,,,all
!線單元的最大編號為540

!次梁共162根
l,62,110
l,64,112
l,69,117
l,71,119
l,97,108
l,73,77
l,80,84
l,90,102
l,91,103
lgen,18,541,549,1,,,3,60

latt,1,,1,,5000,,4
lesize,all,2
lmesh,all
allsel
!線單元的最大編號為702

!樓層板
*do,i,61,64,1
a,i,i+1,i+13,i+12
*enddo
agen,4,1,4,1,,4,,12
*do,i,68,71,1
a,i,i+1,i+13,i+12
*enddo
agen,4,17,21,1,,4,,12
*do,i,89,91,1
a,i,i+1,i+13,i+12
*enddo
agen,18,1,35,1,,,3,60
a,1133,1136,1124,1121

!墻體
!632-635
*do,i,61,64,1
a,i,i+1,i+61,i+60
*enddo

!636-638
*do,i,89,91,1
a,i,i+1,i+61,i+60
*enddo

!639-642
*do,i,68,71,1
a,i,i+1,i+61,i+60
*enddo

!643-646
*do,i,109,112,1
a,i,i+1,i+61,i+60
*enddo

!647-650
*do,i,116,119,1
a,i,i+1,i+61,i+60
*enddo

!651-654
*do,i,61,97,12
a,i,i+12,i+72,i+60
*enddo


!655-656
*do,i,65,77,12
a,i,i+12,i+72,i+60
*enddo

!657-658
*do,i,68,80,12
a,i,i+12,i+72,i+60
*enddo

!659-662
*do,i,72,108,12
a,i,i+12,i+72,i+60
*enddo
agen,17,632,662,1,,,3,60

aatt,2,2,2
aesize,all,2
amesh,all
asel,u,mat,,2
aplot

!筒體

!*do,i,53,1124,60
!a,i,i+3,i+63,i+60
!*enddo
!*do,i,41,1112,60
!a,i,i+3,i+63,i+60
!*enddo
!*do,i,41,1112,60
!a,i,i+12,i+72,i+60
!*enddo
!*do,i,44,1115,60
!a,i,i+12,i+72,i+60
!*enddo
a,53,56,1136,1133
a,41,44,1124,1121
aatt,1,3,2
aesize,all,2
amesh,all

asel,u,mat,,1
aplot

a,41,53,1133,1121
a,44,56,1136,1124
aatt,1,3,2
aesize,all,1.5
amesh,all

nummrg,all
numcmp,all
allsel
finish

/solu
nsel,s,loc,z,0
d,all,all
allsel
antype,modal
modopt,subsp,5
mxpand,5,,,yes
solve
finish