原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
作业4
结11 2011012315 朱尧于
三层钢筋混凝土平面框架，层高4.0m，8度设防(0.2g)，II类场地，设计地震分组为第一组，质量m1=1600kg、m2=1600kg、m3=1200kg，层刚度k1=1800kN/m、K2=1600kN/m、k3=1500kN/m。用振型分解反应谱法计算结构在多遇地震下各振型的地震作用、各振型的层间位移和总层间位移。
解：
II类场地，设计地震分组为第一组：特征周期为：Tg=0.35s；
8度设防（0.2g）αmax=0.16
默认阻尼比为0.05：衰减系数r=0.9；η1=0.02；η2=1
建立结构自由振动的动力方程如下：
m
m
m
其中：
k11=3400KN/m，k
k22=3100KN/m，k
由特征方程求解频率（周期）和振型：
设：X
则：X
代入原方程可得：
(
取系数行列式等于0得：
K=[k11-m1*om,k12,k13
k12,k22-m2*om,k23
k13,k23,k33-m3*om];
px=det(K)=0
px =
-3072000000 x^3 + 16320000000 x^2 - 21576000000 x + 4320000000=0
ω2=242.76 1739.34 3330.40
即：
ω=15.58，41.71，57.71
T=0.40s，0.15s，0.11s
当ω2=242.76,一阶主振型为：X1:X2:X3=1:
当ω2=1739.34,二阶主振型为：X1:X2:X3=1:
当ω2=3330.4,三阶主振型为：X1
计算各振型的振型参与系数：
由Excel表格计算可知：r1=0.536;r2=0.34;r3=0.119
计算各振型的地震影响系数：
第一振型：0.142
第二振型：0.160
第三振型：0.160
计算各振型各质点的水平地震作用和各振型水平地震作用效应：
Xji——j振型i质点的振幅系数
振型组合得到总效应：
总绝对位移：
总层间位移为：
q
附录：
计算表格：见附带Excel（电子版已经发到邮箱）
MATLAB计算源代码：
主过程
clear all;
clc;
k11=3400000;k12=-1600000;k13=0;
k22=3100000;k23=-1500000;k33=1500000;
m1=1600;m2=1600;m3=1200;
syms om;
K=[k11-m1*om,k12,k13
k12,k22-m2*om,k23
k13,k23,k33-m3*om];
b=det(K) %得到矩阵行列式
求解行列式为0 的多项式方程得到特征根ω2
p=[- 3072000000,16320000000000, - 21576000000000000,4320000000000000000]
px=poly2str(p,x)
format rat
r=roots(p)
分别把三个ω2代入原方程组的前两个方程中求解主振型
z=solve((3400000-1600*1739.34)*1-1600000*x2+0*x3=0,-1600000*1+(3100000-1600*1739.34)*x2-1500000*x3=0,x2,x3)
a2=z.x2
a3=z.x3
文档评论（0）