|
!地处Ⅱ类场地,抗震设防烈度为八度,近震
!反应谱最大地震影响系数αｍａｘ=0.16,
!特征周期Ｔｇ=0.30.结构阻尼比取为0.05.
/PREP7
BLOCK,0,0.15,0,2.85,0,-0.15,
BLOCK,0.15,3,0,2.85,0,-0.15,
BLOCK,3,6,0,2.85,0,-0.15
BLOCK,6,9,0,2.85,0,-0.15
BLOCK,0,0.15,0,2.85,-0.15,-4.5,
BLOCK,0,0.15,2.85,3,0,-0.15,
BLOCK,0.15,3,2.85,3,0,-0.15
BLOCK,3,6,2.85,3,0,-0.15
BLOCK,6,9,2.85,3,0,-0.15
BLOCK,0,0.15,2.85,3,-0.15,-4.5,
BLOCK,0.15,3,2.85,3,-0.15,-4.5
Block,3,6,2.85,3,-0.15,-4.5
Block,6,9,2.85,3,-0.15,-4.5
BLOCK,0,0.15,3,6,0,-0.15,
BLOCK,0.15,3,3,6,0,-0.15,
BLOCK,3,6,3,6,0,-0.15,
BLOCK,6,9,3,6,0,-0.15,
BLOCK,0,0.15,3,6,-0.15,-4.5,
vsel,s,,,all
vglue,all
/pnum,area,on
aplot
vsel,s,,,all
cm,f_sub,volum !前面部分
cmsel,s,f_sub
vsymm,z,all,,,,0,0
vsel,s,,,all
vsel,u,,,f_sub
vgen,,all,,,,,-9,,,1
vsel,s,,,all
vglue,all
!定义单元
et,1,solid45
mp,ex,1,2.5e8
mp,nuxy,1,0.2
mp,dens,1,2500
lsel,s,,,all
lsel,s,length,,,0.15
lesize,all,,3
lsel,s,,,all
lsel,u,length,,0.15
lesize,all,0.15
mshape,0,3D
mshkey,1
vsel,s,,,all
vmesh,all
asel,s,loc,y,0
cm,base_f,area !基础面
asel,s,,,all
asel,s,loc,x,9
cm,l_wall,area !x对称面
allsel
FINISH
/solu
antype,modal !求振型和频率
modopt,subsp,4 !用子空间迭代法求前3阶振型
asel,s,,,all
cmsel,s,base_f
nsla,s,1
d,all,all
asel,s,,,all
asel,s,,,l_wall
nsla,s
dsym,symm,x
allsel
mxpand,5
solve
FINISH
/post1
set,list,2
set,first
/view,1,1,1,1,
/ang,1
/rep,fast
pldisp,0
anmode,10,0.5e-1,20
/anfile,save,''1'',''avi'',''''
set,next
/view,1,1,1,1,
/ang,1
/rep,fast
pldisp,0
anmode,10,0.5e-1,20
/anfile,save,''2'',''avi'',''''
set,next
/view,1,1,1,1,
/ang,1
/rep,fast
pldisp,0
anmode,10,0.5e-1,20
/anfile,save,''3'',''avi'',''''
/solu
antype,spectr !求振型和频率
spopt,sprs,3,yes !单点谱分析
svtype,2 !加速度谱
sed,1,0,0 !X方向地震作用
!定义反应谱(抗震规范GBJ11-89)：Tg=0.4s,amax=0.16(8度) !F=a*G=a*m*g a值乘以9.82
freq,0.333333333,0.353773168,0.376883475,0.403224181,0.433523506,0.468746338,0.510198876,0.559694183,0.619824466,
sv,,0.256257352,0.270357153,0.286201463,0.304142684,0.324636126,0.348280823,0.375881551,0.408546201,0.447844563,
freq,0.694429977,0.789452909,0.914603474,1.086909265,1.339207992,1.744044089,2.5,10,1000
sv,,0.496079267,0.556774157,0.635616298,0.742436687,0.895877629,1.136284746,1.5712,1.5712,0.70704
solve
finish
/solu
antype,modal !求振型和频率
expass,on
mxpand,3,,,yes !扩展前3阶振型，计算单元解和反力
solve
finish
/solu
antype,spectr
srss,1,disp !应用SRSS方法计算单元应力,该命令将写*.mcom文件
solve
finish
/post1
set,list
/inp,,mcom
pldisp
pletab,avg
prnsol,dof
presol,elem
prrsol,f
|