本帖最后由 西风独自凉 于 2011-8-19 13:03 编辑
最近一工程要计算地震作用,故先作一简单算例考题,贴之于大家共享
!本文计算例题取自:
!建筑结构抗震设计原理,同济大学出版社,朱伯龙,2000年12月第5次印刷,p68,例题3-7
发信人: gxn (老农), 信区: FEA
标 题: Re: 一个振型分解法算例
发信站: 同舟共济站 (2002年11月04日21:54:56 星期一), 站内信件
标 题: Re: 一个振型分解法算例
发信站: 同舟共济站 (2002年11月04日21:54:56 星期一), 站内信件
!---------------------------------------------------------------------------
FINISH
/CLEAR
/prep7
et,1,mass21 !质量单元
r,1,116620 !m1=116620kg
r,2,110850 !m2=110850kg
r,3,59450 !m3=59450kg
et,2,combin14 !定义弹簧单元模拟层间刚度
keyopt,2,3,2 !两维弹簧
r,4,30910000 !k1=30910kN/m
r,5,36900000 !k2=36900kN/m
r,6,37440000 !k3=37440kN/m
!---------------------------------------------------------------------------
n,1,0,0,0
n,2,6.5,0,0
n,3,12.5,0,0
n,4,17.5,0,0
type,1 !建立单元m1
real,1
e,2
type,1 !建立单元m2
real,2
e,3
type,1 !建立单元m3
real,3
e,4
type,2 !建立单元k1
real,4
e,1,2
type,2 !建立单元k2
real,5
e,2,3
type,2 !建立单元k3
real,6
e,3,4
d,1,all !约束节点1
!---------------------------------------------------------------------------
/solu
antype,modal !求振型和频率
modopt,subs,3 !用子空间迭代法求前3阶振型
solve
finish
!---------------------------------------------------------------------------
/solu
antype,spectr !求振型和频率
spopt,sprs,3,yes !单一反应谱
sed,1 !X方向地震作用
svtype,2 !加速度谱
!定义反应谱(抗震规范GBJ11-89):Tg=0.4s,amax=0.16(8度)
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,,disp !应用SRSS方法计算单元应力,该命令将写*.mcom文件
solve
finish
!---------------------------------------------------------------------------
/post1
set,,,,,,,1 !画出第1模态的弹簧单元轴力图
etable,N1,smisc,1
plls,N1,N1,1,0
set,,,,,,,2 !画出第2模态的弹簧单元轴力图
etable,N2,smisc,1
plls,N2,N2,1,0
set,,,,,,,3 !画出第3模态的弹簧单元轴力图
etable,N3,smisc,1
plls,N3,N3,1,0
/input,,mcom !输入*.mcom文件后可查看SSRS内力组合值
etable,N,smisc,1 !画出弹簧单元轴力图组合值
plls,N,N,1,0
!结果比较:
!各模态的振型参与系数(ansys中可查看*.mcom文件)
! mode1 mode2 mode3
!ansys: 6.28804 0.324219 0.0366454
!SSRS内力组合值:
! V1 V2 V3
!Ansys: 253.21 181.38 73.405
!参考文献: 253.563 181.578 73.445
!第1模态内力值(第1模态振型系数为6.28804):
! V1 V2 V3
!Ansys: 40.089 28.429 10.787
!Ansys乘以振型系数后: 252.081 178.763 67.829
!参考文献: 252.438 179.043 67.909
!第2模态内力值(第2模态振型系数为0.324219):
! V1 V2 V3
!Ansys: 73.689 -92.450 -84.405
!Ansys乘以振型系数后: 23.891 -29.974 -27.365
!参考文献: 23.831 -29.930 -27.324
!第3模态内力值(第3模态振型系数为0.0366454):
! V1 V2 V3
!Ansys: 33.994 -120.840 170.250
!Ansys乘以振型系数后: 1.246 -4.428 6.239
!参考文献: 1.192 -4.294 6.057
|