框架模态分析final.txt

是用apdl语言编的

/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