地震工程学大作业

合集下载

大工19秋《工程抗震》大作业题目及要求

大工19秋《工程抗震》大作业及要求注意：从以下五个题目中任选两个进行解答（注意：从题目一、二中选择一道计算题，并从题目三、四、五中选择一道问答题，分别进行解答，不可同时选择两道计算题或者问答题）；解答前，需将所选题目复制（使老师明确你所选的题目）。

题目一：底部剪力法。

钢筋混凝土5层框架经质量集中后计算简图如下图所示，各层高均为3m ，集中于各楼层的重力荷载代表值分别为：1500kN G =，2550kN G =，3580kN G =，4600kN G =，5450kN G =。

结构阻尼比0.05ξ=，自振周期为10.55s T =，Ⅰ1类场地类别，设计地震分组为第一组，抗震设防烈度为8度（设计基本地震加速度为0.30g ）。

按底部剪力法计算结构在多遇地震时的水平地震作用及地震剪力。

3580kN=2550kN=1500kN=（a ）计算简图4600kN=5450kN=题目二：振型分解反应谱法。

某3层钢筋混凝土框架结构，阻尼比0.05ξ=，部分尺寸如下图所示。

各楼层重力荷载代表值为11200kN G =，21000kN G =，3650kN G =，场地土Ⅱ类，设防烈度8度（设计基本地震加速度为0.20g），地震分组在第二组。

现算得前3个振型的自振周期为10.68sT=，20.24sT=，30.16sT=，振型如下图所示。

试用振型分解反应谱法求该框架结构的层间地震剪力标准值。

3212.1481.7351.000 1.000 1.000-1.138 1.4670.139-1.316（a）框架（b）计算简图（c）第一振型（d）第二振型（e）第三振型题目三：怎样判断土的液化？如何确定土的液化严重程度，并简述抗液化措施。

题目四：什么叫鞭端效应？设计时如何考虑这种效应？题目五：多高层钢筋混凝土结构抗震设计的防震缝设计有哪些要求？作业具体要求：1. 封面格式（见文件最后部分）封面名称：大连理工大学《工程抗震》大作业，字体为宋体加黑，字号为小一；姓名、学号、学习中心等字体为宋体，字号为小三号。

地震工程学作业

地震动反应谱计算与分析廖志娟 0920020164一、概论地震动反应谱是将特定地震动输入固有周期和阻尼比各不相同的线性单自由度体系中，计算单自由度体系响应（一般是相对位移、相对速度、绝对加速度）的最大值，将响应最大值表达成固有周期和阻尼比的函数的形式。

地震动反应谱建立了地震动特性与结构动力反应之间的桥梁。

在本质上，地震动反应谱反映了地震动频谱特性，同时，它又描述了一般结构地震反映的某些基本特征。

一般来说，反应谱有如下形式：(,)S S T ξ=考察一受地面加速度)(t a 作用的线性单自由度体系，)(),(),(t x t xt x 分别为单自由度体系的相对位移、相对速度和相对加速度，)()()(t a t x t y += 为绝对加速度，ξ为体系阻尼比，ω为体系基本频率。

2/s cm设初始条件为0)0()0(====t xt x 的单自由度体系，由结构动力学基本原理，可得到关于)(t a 的相对位移、相对速度和绝对加速度反应谱(d ωω≈)：max)(max )](sin[)(1)(),(⎰-==--tt d d t e a t x T S ττωτωξτξω max)(max )](cos[)()(),(⎰-==--t t v d t ea t xT S ττωτξτξωmax 2max )()(),(t x t yT S a ωξ== 可以证明，对于中频及高频体系，拟速度差不多等于最大相对速度，或当体系具有中等阻尼时，在从很低频率到很高频率的整个频率范围内，拟加速度与最大加速度的差别都不大。

总之，在周期不是很长，即频率不是很小的情况下，可以用拟相对速度谱、拟绝对加速度谱和相对位移谱近似计算以上反应谱而不会产生很大的误差。

max)()](sin[)(1),(⎰-=--tt d d t e a T S ττωτωξτξω),(),(ξωξT S T S d v =),(),(2ξωξT S T S d a =根据上式所表示的关系，我们可以用对数坐标把位移、速度和加速度反应谱画在一张图上，通常称它为三联反应谱。

大工19秋《工程抗震》大作业题目及要求标准答案

大工19秋《工程抗震》大作业题目及要求标准答案XXX《工程抗震》大作业研究中心：姓名：学号：题目二：底部剪力法。

一个六层钢筋混凝土框架结构建造于基本烈度为7度的区域，场地为Ⅱ类，设计地震分组为第二组。

结构首层层高为4m，其余各层层高为3m。

经质量集中后，各楼层G1=900kN，G2=G3=G4=G5=700kN，G6=650kN。

重力荷载代表值分别为：该结构的基本周期为0.65s，阻尼比为ξ=0.05.求各层水平地震作用标准值。

解：根据题目可知，抗震设防烈度为7度（由于题目没有给出基本地震加速度的值，一般理解为0.15g），且为多遇地震。

查《建筑设计抗震规范》表5.1.4-1表可得：a max=0.08.由题目可知，设计地震分组为第二组，Ⅱ类场地类别，查《建筑设计抗震规范》表5.1.4-2表可得：Tg=0.40s。

钢筋混凝土结构阻尼比题目已给出，为ξ=0.05，则水平地震影响系数曲线衰减系数为：γ=0.9+（0.05-ξ）/(0.3+6ξ)=0.9阻尼比调整系数为：η2=1+(0.05-ξ)/(0.08+1.6ξ)=1则：α1=（Tg/T1）rη2αmax=（0.40/0.65）0.9×1.0×0.08≈0.05由于T1=0.65s＞1.4Tg=1.4×0.40=0.56s，0.35<Tg=0.40s≤0.55.查《建筑设计抗震规范》表5.2.1表可知，δn=0.08T1+0.01=0.08×0.65+0.01≈0.06总水平地震作用标准值为：FEk=α1Geq=0.05×（900+700+700+700+700+650）×85%=184.88KN各楼层水平地震作用标准值为：Fi=GiHiFEk(1-δn)/∑GjHj (i=1,2,3，n)GjHj=900×4+700×7+700×10+700×13+700×16+650×19=KN ·mF1k=[900×4×184.88×(1-0.06)]/=12.99KNXXX[700×7×184.88×(1-0.06)]/=17.69KNXXX[700×10×184.88×(1-0.06)]/=25.27KNXXX 32.84KN。

工程结构抗震大作业7页

工程结构抗震大作业专业：土木工程姓名：张卫东学号：124507891.结构场地简述设计类型：宿舍类型：钢筋混凝土框架结构（3跨3层）结构设防烈度为8度，设计基本地震加速度为0.20g。

2.结构平面图3结构横向立面图4.结构框架简述框架梁：b×h=300mm×60;框架柱B×H=600mm×600mm;层高1~3层：h=4.0m;混凝土采用C30，E=30×106kN/m25.计算框架梁内力我们选出一品中间框架进行计算。

计算步骤如下：（1）计算梁的线刚度横向框架梁线刚度计算表楼层层高（m ）H （m m ）Bbh部位跨度l (m )I O (m 4)I b (m 4)JE c (k N /m 2)i b =E c I b /l (k N /m )一层中间梁 4 600 600 300 600 边 4 0.0054 0.010830000000C 30;30.0×10681000 一层边梁4600600300600边40.00540.008130000000 60750 一层中间梁 4600600300600中20.00540.010830000000162000一层边梁 4 600 600 300 600 中 2 0.0054 0.0081 30000000 121500(2)计算柱的抗侧刚度（3）计算柱的D值ic1 (KN*m)(第一排) i2(KN*m)(第一排)KaCk D首层边柱（A3A4）16000 60750 3.796875 0.741239892 12000 8894.878706首层中间柱（B3B4）16000 121500 11.390625 0.88798133 12000 10655.77596ic1(KN *m)(第二排) i2(KN*m)(第二排)KaCk D首层边柱（A3A4）16000 81000 5.0625 0.787610619 12000 9451.327434首层中间柱（B3B4）16000 162000 15.1875 0.912727273 12000 10952.72727ic2(KN *m)(第一排) i2(KN*m)(第一排)KaCk D一般层边柱（A2A3) 16000 60750 3.796875 0.654986523 12000 7859.838275E(KN/m2) B1(mm) H1(mm) I(m4) H(m) i(KN*m) k(KN/m) 一层柱30000000 400 400 0.002133333 4 16000 12000 二层柱30000000 400 400 0.002133333 4 16000 12000 三层柱30000000 400 400 0.002133333 4 16000 12000一般层中间柱（B2B3)16000 121500 11.390625 0.850641774 12000 10207.70128 ic1(KN *m)(第二排) i2(KN*m)(第二排) K a Ck D一般层边柱（A2A3) 16000810005.0625 0.716814159 12000 8601.769912一般层中间柱（B2B3)16000 162000 15.1875 0.883636364 12000 10603.63636（4）计算框架梁的横向基本自震周期T 1G i (kN)Vi(kN)∑D(KN/m)△U(m)U(m)一层 400012100201926.9469 0.060310941 0.060310941二层 4000 8100 187502.5959 0.043199402 0.103510343三层4100 4000 104452.6233 0.039252245 0.142762588 T 1=1.7ΨT u T ^0.5=1.7×0.7×0.142763^0.5=0.449629497(s)(5)计算各楼层剪力V（6）多遇地震作用下的层间位移验算表多遇地震作用下的层间位移验算表G eq (kN )αmax Tg(s) α 1 T(s) F Ek (kN) 1.4Tgδ710285 0.16 0.25 0.0943401750.449629497 970.2886952 0.35 0.10597036ΔF 7(kN ) F1(k N)F2(kN) F3(kN) ΔF 3(kN) V1(kN) V2(kN) V3(kN) 102.8218421 142.792898285.5857953 439.0881602102.8218421 970.2886952 827.4957976 541.9100023楼层层间剪力Vi(KN/m) 层间刚度∑Di(K N/m)层间弹性位移Δu i =Vi/∑Di(m)顶点位移∑ui(m) 层高h(m) 层间弹性位移角Δui/h3 541.910002 104452.62330.0051880940.01443762741/7702 827.495798 187502.59590.00441325 0.00924953341/9061 970.288695 201926.94690.0048362830.00483628341/827（7）水平地震作用下中间框架柱地震剪力及柱端弯矩计算水平地震作用下中间框架柱地震剪力及柱端弯矩计算表楼层 ∑D i (K N /m ) 部位 D i k(K N /m )V i (k N ) V i k (k N )yh (m )M i k上(K N ·m ) M i k下(K N ·m )3 104452.623 边柱7859.838275 541.910040.778 0.5 2 81.55515568 81.55515568 中间柱 10603.63636 55.0130.52 110.0253191 110.0253191 2 187502.596边柱7859.838275827.495734.687 0.54 2.1663.82487093 74.92484848中间柱 10603.63636 46.7960.5293.59299255 93.592992551201926.947边柱9451.327434970.288695245.7090.552.282.27672976100.5604475中间柱10655.7759651.5340.552.292.76182689 113.3755662（8）反弯点高度计算y 0y 1y 2y 3yh ´(m)A1A2 0.5 0 0 0 0.5 2 A2A3 0.54 0 0 0 0.54 2.16 A3A4 0.55 0 0 0 0.55 2.2 B1B2 0.5 0 0 0 0.5 2 B2B3 0.5 0 0 0 0.5 2 B3B40.550 0 0 0.552.2（9）水平地震作用下中间框架梁端弯矩计算楼层节点i 节点l (kN ·m) i 节点r (kN ·m) M i+1,k (上)(kN ·m)(柱) M i+1,k (下)(kN·m)（柱）M 节点，r(kN ·m)3 边 —— 81000 81.55515568 —— 81.6 中 81000 162000 110.0253191 —— 73.4 2 边 —— 81000 63.82487093 81.555 145 中 81000 162000 93.59299255 110.03 136 1 边 —— 81000 82.27672976 74.925 157 中 8100016200092.76182689 93.593124（10）水平地震作用下中间框架梁梁端剪力及柱轴力计算楼层边跨梁走道梁柱轴力 Mb l Mb r l(m ) Vb(kN) Mb l (kN ·m) l(m ) Vb(k N) 边柱（kN ）中柱（kN ）(kN ·m)3 81.5551557 36.6754 29.55756551 73.35 2 73.4 -29.5575655 -43.79264722 145.380027 67.8734 53.31319929 135.75 2 136 -82.8707648 -126.2249891 157.201978 62.1184 54.82996285124.24 2 124 -137.700728-195.631572（11）框架梁内力图（12）。

大工19春《工程抗震》大作业题目及要求大作业答案可直接上传

⼤⼯19春《⼯程抗震》⼤作业题⽬及要求⼤作业答案可直接上传⼤连理⼯⼤学《⼯程抗震》⼤作业学习中⼼： XXX姓名： XXX学号： XXXXX题⽬⼀：底部剪⼒法。

钢筋混凝⼟框架结构共7层，经质量集中后得到集中于各楼层的重⼒荷载代表值分别为：1900kN G =，23456700kN G G G G G =====，7600kN G =。

该结构每层层⾼均为3.5⽶，结构阻尼⽐0.05ξ=，⾃振周期为10.55s T =，Ⅰ1类场地类别，设计地震分组为第⼀组，抗震设防烈度为8度（设计基本地震加速度为0.30g ）。

按底部剪⼒法计算结构在多遇地震时的⽔平地震作⽤及地震剪⼒。

解：由题⽬可知，抗震设防烈度为8度（设计基本地震加速度为0.30g ），且⼯程有多遇地震的特点，查表可得：αmax=0.24。

由题⽬可知，设计地震分组为第⼀组，Ⅰ1类场地类别，查下表可得：Tg=0.25s 。

钢筋混凝⼟结构阻尼⽐题⽬已给出，为ξ=0.05，则：⽔平地震影响系数曲线衰减系数：γ=0.9+（0.05-ξ）/(0.3+6ξ)=0.9阻尼⽐调整系数：η2=1+(0.05-ξ)/(0.08+1.6ξ)=1则：α1=（Tg/T1）rη2αmax=（0.25/0.55）0.9×1.0×0.24≈0.12因为：T1=0.55s＞1.4Tg=1.4×0.25=0.35s，Tg=0.25s≤0.35。

查下表可知，δn=0.08T1+0.07=0.08×0.55+0.07≈0.11总⽔平地震作⽤标准值：FEk=α1Geq=0.12×（900+700+700+700+700+700+600）×85%=510 KN △Fn =δn F EK = 0.11×510=56.1 KN各楼层⽔平地震作⽤标准值：Fi=GiHiFEk(1-δn)/∑GjHj (i=1,2,3，n)∑GjHj=900×3.5+700×7+700×10.5+700×14+700×17.5+700×21+600×24.5=66850KN?mF1k=[900×3.5×510×(1-0.11)]/66850=21.39KNF2k=[700×7×510×(1-0.11)]/66850=33.27KNF3k=[700×10.5×510×(1-0.11)]/66850=49.91KNF4k=[700×14×510×(1-0.11)]/66850=66.54KNF5k=[700×17.5×510×(1-0.11)]/66850=83.18KNF6k=[700×21×510×(1-0.11)]/66850=99.81KNF7k=[600×24.5×510×(1-0.11)]/66850=99.81KN题⽬三：简述建筑结构抗震计算应遵循的五个原则。

(完整版)大连理工大学(工程抗震)大作业

大连理工大学《工程抗震》大作业题目1：底部剪力法。

结构阻尼比0.05ξ=，自振周期为10.55s T =，Ⅰ1类场地类别，设计地震分组为第一组，抗震设防烈度为8度（设计基本地震加速度为0.30g ）。

按底部剪力法计算结构在多遇地震时的水平地震作用及地震剪力。

3580kN=2550kN =1500kN=（a ）计算简图4600kN=5450kN=解：查《建筑设计抗震规范》表5.1.4-1.2知，8度多遇地震，αmax=0.16 设计地震分组为第一组， Ι类场地，取Tg=0.25s Tg=0.25＜T1=0.55s ＜5Tg=1.25sα1=（Tg/T1）r η2αmax =（0.25/0.55）0.9×1.0×0.16=0.079≈0.08 查《建筑设计抗震规范》表5.2.1知，T 1=0.55s ＞1.4Tg=1.4×0.25=0.35s 取δn=0.08 T1+0.07=0.08×0.55+0.07=0.11 总水平地震作用标准值：F EK =α1Geq=0.08×（500+550+580+600+450）×85%=182.2KN各楼层水平地震作用标准值： Fi=G i H i F EK (1-δn)/∑G j H j (i=1,2,3…n)∑G j H j =500×3 +550×6+580×9+600×12+450 ×15=23970KN ·m F 1=[500×3×182.2×(1-0.11)]/23970=10.15KN F 2=[550×6×182.2×(1-0.11)]/23970=22.32KN F 3=[580×9×182.2×(1-0.11)]/23970=35.31KN F 4=[600×12×182.2×(1-0.11)]/23970=48.71KN F 5=[450×15×182.2×(1-0.11)]/23970=45.66KN 计算各楼层的层间地震剪力V 1= F 1+ F 2+ F 3+ F 4+ F 5=10.15+22.32+35.31+48.71+45.66=162.15KN V 2= F 2+ F 3+ F 4+ F 5=22.32+35.31+48.71+45.66=152KN V 3= F 3+ F 4+ F 5=35.31+48.71+45.66=129.68KN V 4= F 4+ F 5=48.71+45.66=94.37KN V 5=F 5=45.66KN题目3：怎样判断土的液化？如何确定土的液化严重程度，并简述抗液化措施。

大工17春《工程地震性》在线作业

大工17春《工程地震性》在线作业一、作业要求本学期《工程地震性》在线作业总共包含以下几个要求：1. 学生需要通过学校指定的在线平台完成作业；2. 作业内容涵盖了本学期研究的各个章节，包括地震基础知识、地震力计算、结构响应等；3. 每个章节的作业都包含了选择题、判断题和简答题，学生需要根据掌握的知识完成答题；4. 作业完成后，学生需要及时提交答案，截止日期为本学期末。

二、作业注意事项为了顺利完成本学期的《工程地震性》在线作业，请注意以下几个事项：1. 学生需要保证网络畅通，并及时登录学校指定的在线平台；2. 学生需要仔细阅读每个章节的作业要求，对于有疑问的地方可以向老师咨询；3. 学生需要遵守学校的学术诚信规定，不得抄袭他人作业；4. 学生需要按时完成作业，并按要求提交答案；5. 学生需要认真复每个章节的知识点，以便更好地完成作业。

三、作业评分标准本学期《工程地震性》在线作业的评分参考以下标准：1. 选择题和判断题：每题1分，错误不扣分；2. 简答题：根据回答的完整度和准确性进行评分，满分为5分；3. 作业完成情况：作业按时提交且答案正确，满分为10分。

四、作业指导与解答本学期《工程地震性》在线作业的指导与解答将由老师在每周课程结束后发布在学校指定的在线平台上。

学生可以参考指导和解答，提高作业完成的准确性。

五、作业完成后的操作学生在完成作业后，请进行以下操作：1. 确认答案正确无误；2. 及时提交答案，确保在截止日期之前完成。

如有其他疑问，请及时向老师咨询。

谢谢大家的配合！祝大家顺利完成本学期的《工程地震性》在线作业！。

大工13秋《工程抗震》大作业题目及要求

大工13秋《工程抗震》大作业题目及要求注意：从以下五个题目中任选两个进行解答（注意：从题目一、二中选择一道计算题，并从题目三、四、五中选择一道问答题，分别进行解答，不可同时选择两道计算题或者问答题）；解答前，需将所选题目复制（使老师明确你所选的题目）。

题目一：求自振圆频率和振型。

三层剪切型结构如下图所示，求该结构的自振圆频率和振型。

3m 3m 4m198MN m=2196MN m=3245MN m=180tm =题目二：底部剪力法。

某四层钢筋混凝土框架结构，建造于基本烈度为8度区域，设计基本地震加速度为0.30g ，场地类别为Ⅱ类，设计地震分组为第二组，结构层高和重力代表值见下图。

取一榀典型框架进行分析，结构的基本周期为0.5s 。

要求：用底部剪力法计算在多遇地震时的水平地震作用及地震剪力。

5m4m 4m4m4500kNG =3800kNG =2800kN G =11000kNG =题目三：简述抗震设防中“三水准两阶段设计”的内容。

题目四：简述地基基础抗震验算的原则。

哪些建筑可不进行天然地基及基础的抗震承载力验算？题目五：多层砌体结构房屋的抗震构造措施包括哪些方面，请逐一详细阐述。

作业具体要求：1. 封面格式（见文件最后部分）封面名称：大连理工大学《工程抗震》大作业，字体为宋体加黑，字号为小一；姓名、学号、学习中心等字体为宋体，字号为小三号。

2. 文件名大作业上交时文件名写法为：[姓名奥鹏卡号（学号）学习中心]（如：戴卫东101410013979浙江台州奥鹏学习中心[1]VIP）；离线作业需用word文档写作，不要提交除word文档以外的文件，word文档无需压缩。

以附件形式上交离线作业（附件的大小限制在10M以内），选择已完成的作业（注意命名），点提交即可。

如下图所示。

截止时间：2014年3月5日前。

3. 正文格式作业正文内容统一采用宋体，字号为小四。

注意：作业应该独立完成，不准抄袭其他网站或者请人代做，如有雷同作业，成绩以零分计。

大工15秋《工程抗震》大作业答案

大连理工大学《工程抗震》大作业学习中心：姓名：学号：题目一：底部剪力法钢筋混凝土7层框架结构各层高均为3m ，经质量集中后，各楼层的重力荷载代表值分别为：11000kN G =，234567750kN G G G G G G ======，。

结构阻尼比0.05ξ=，自振周期为10.80s T =，Ⅱ类场地类别，设计地震分组为第二组，抗震设防烈度为7度。

按底部剪力法计算结构在多遇地震时的水平地震作用及地震剪力。

解：设防烈度对应的多遇地震和罕遇地震的αmax 值场地特征周期Tg （S ）（1）抗震设防烈度为7度，多遇地震下，查表得αmax=0.12，Ⅱ类场地类别，设计地震分组为第二组，查表得Tg=0.40s ，所以Tg=0.40<T 1=0.80<5Tg=2.00 （2）计算地震影响系数αiKNG G G T T E g i 55257006800100085.085.007.012.0180.040.0)eq eq9.0max 21=+⨯+⨯==≈⨯⨯==）（计算结构等效总重力）（（αηαγ（3）计算结构底部剪力F EKKN G 8.386552507.0F eq 1EK =⨯==α （4）计算各质点的水平地震作用F i已知Tg=0.40s ，T 1=0.80s >1.4Tg=0.56s 。

该结构为钢筋混凝土房屋结构，需要考虑结构顶部附加集中作用。

m*84600247002118151296800310004336.327162.288.386074.0074.001.080.008.001.008.0n1i ii n n n 1n KN HG KNF KN F F T EK EK =⨯++++++⨯+⨯==-=⨯==∆=+⨯=+=∑=）（）（则δδδKNF HG H G F KN F HG H G F KN F HG H G F KN F HG H G F KNF HG H G F KNF HG H G F KNF HG H G F KNF HG H G F EK EK EK EKEK EK EK EK 431.138.38684600310001946.218.3868460068001919.328.3868460098001892.438.38684600128001865.548.38684600158001838.658.38684600188001811.768.38684600218001784.878.38684600247001n n1i ii 111n n1i ii 222n n1i ii 333n n1i ii 444n n1i ii 555n n1i ii 666n n1i ii 777n n1i ii 888=⨯⨯=-==⨯⨯=-==⨯⨯=-==⨯⨯=-==⨯⨯=-==⨯⨯=-==⨯⨯=-==⨯⨯=-=∑∑∑∑∑∑∑∑========）（）（）（）（）（）（）（）（δδδδδδδδ（5）计算各楼层地震剪力KN V F V KN V F V KN V F V KN V F V KN V F V KN V F V KN V F V KN F F V 652.423221.410431.13221.410275.388946.21275.388356.355919.32356.355464.311892.43464.311599.256865.54599.256761.190838.65761.19095.113811.7695.113166.26784.87211322433544655766877n 88=+=+==+=+==+=+==+=+==+=+==+=+==+=+==+=∆+=题目四：多层小砌块房屋的芯柱应符合哪些构造要求？答：1、小砌块房屋芯柱截面不宜小于120mm×120mm.2、芯柱混凝土强度等级，不应低于Cb20.3、芯柱的竖向插筋应贯通墙身且与圈梁连接；插筋不应小于1φ12，6、7度时超过五层，8度时超过四层和9度时，插筋不应小于1φ14.4、芯柱应伸入室外地面下500mm或与埋深小于500mm的基础圈梁相连。

抗震大作业

1、2、建筑概况五层钢筋混凝土框架结构，设防烈度7度，Ⅰ类场地，设计地震分为第一组。

混凝土等级：梁为C25，柱为C30。

主筋采用HRB335钢，箍筋采用HPB235级钢，抗震等级为二级框架。

梁截面尺寸为250mm*600mm,柱截面尺寸为500mm*500mm 。

3、重力荷载代表值计算重力荷载代表值时，按《抗震规范》规定，楼面均布荷载的组合系数取0.5，屋面活荷载不予考虑。

屋面恒荷载标准值计算：20mm 厚水泥砂浆找平20.02200.4/kN m ⨯=100～140mm 厚（2%找坡）膨胀珍珠岩20.140.1070.84/2kN m +⨯= 100mm 厚现浇钢筋混凝土楼板 20.1025 2.5/kN m ⨯=15mm厚纸筋石灰抹底20.015160.24/kN m ⨯=屋面恒荷载标准值：23.98/kN m楼面恒荷载标准值计算：25mm 厚水泥砂浆面层20.025200.50/kN m ⨯=100mm 厚现浇钢筋混凝土楼板20.1025 2.5/kN m ⨯=15mm厚纸筋石灰抹底20.015160.24/kN m ⨯=楼面恒荷载标准值：23.24/kN m框架梁线荷载标准值（包括粉刷层）：0.250.6252(0.60.1)0.0217 4.09/kN m ⨯⨯+⨯-⨯⨯=框架柱的线荷载标准值（包括粉刷层）：0.50.525(0.5440.24)0.0217 6.61/kN m ⨯⨯+⨯-⨯⨯⨯=填充墙（包括粉刷）的恒荷载标准值：20.24190.02217 5.24/kN m ⨯+⨯⨯=楼面活荷载标准值：21.5/kN m 则楼面重力荷载代表值：5(3.987810 3.98285)(4.09820 4.09710 4.0925)(6.61424)/2 5.24(80.5)(40.6)20 5.24(70.5)(40.6)12 5.24(20.5)(40.6)6/25957G kN=⨯⨯⨯+⨯⨯⨯+⨯⨯+⨯⨯+⨯⨯+⨯⨯+⨯-⨯-⨯+⨯-⨯-⨯+⨯-⨯-⨯=4(3.247810 3.24285)(4.09820 4.09710 4.0925)(6.61424)/22 5.24(80.5)(40.6)20 5.24(70.5)(40.6)12 5.24(20.5)(40.6)6/22(0.5 1.578100.5 1.5285)8392G kN=⨯⨯⨯+⨯⨯⨯+⨯⨯+⨯⨯+⨯⨯+⨯⨯⨯+⨯-⨯-⨯+⨯-⨯-⨯+⨯-⨯-⨯⨯+⨯⨯⨯⨯+⨯⨯⨯⨯=38392G kN =28392G kN =1(3.247810 3.24285)(4.09820 4.09710 4.0925)(6.61424 6.61524)/2 5.24(80.5)(40.6)20 5.24(70.5)(40.6)125.24(20.5)(40.6)6 5.24(80.5)(50.6)20 5.24(70.5)(50.6)G =⨯⨯⨯+⨯⨯⨯+⨯⨯+⨯⨯+⨯⨯+⨯⨯+⨯⨯+⨯-⨯-⨯+⨯-⨯-⨯+⨯-⨯-⨯+⨯-⨯-⨯+⨯-⨯-⨯125.24(20.5)(50.6)6/2(0.5 1.578100.5 1.5285)9713kN+⨯-⨯-⨯+⨯⨯⨯⨯+⨯⨯⨯⨯=Gi 为楼层重力荷载代表值总重力荷载代表值G=40846Gi kN =∑4、框架刚度（1）梁的刚度计算结果列于下表部位断面b*h跨度矩形截面惯性边跨梁中跨梁Ib=1.5I ib=EIb/Ib=2I0 ib=EIb/L4)强度等级为C25的混凝土弹性模量：622810/c E kN m =⨯ 强度等级为C30的混凝土弹性模量：623010/c E kN m =⨯（2）、柱的刚度计算上部各层柱：3634410.50.5/41230101/120.50.5/4 3.9110c c i E kN m =⨯⨯⨯=⨯⨯⨯⨯=⨯⋅ 底层柱： 3634410.50.5/51230101/120.50.5/5 3.1210c c i E kN m =⨯⨯⨯=⨯⨯⨯⨯=⨯⋅5、自振周期计算假想顶点侧移T u ，计算结果列于下表结构基本自震周期考虑非结构墙影响的折减系数00.6α=1.7 1.70.63l T s α==⨯=6、多遇水平地震作用计算设防烈度7度，Ⅰ类场地时，设计地震分组为第一组，max0.08,0.25g T α==，0.90.9max 0.25()0.080.02770.65g l l T T αα⎛⎫==⨯= ⎪⎝⎭由于 1.4l g T T >，需附加顶部集中力，0.080.070.080.650.070.122n l T δ=+=⨯+=结构总水平地震作用效应标准值为0.02770.8540846962EK l eq F G kN α==⨯⨯=附加顶部集中力为0.122962117.4n n EK F F kN δ∆==⨯= 1(1)i ii EK n njjj G H F F G Hδ==-∑844.6EK n F F kN -∆=7、框架各层地震力及弹性位移计算结果列于下表多遇地震下楼层剪力和楼层弹性位移8、在水平地震作用下框架内力分析取出一榀中间框架，计算结果如下图9、框架重力荷载效应计算恒荷载作用于顶层的框架梁上的线荷载为：51515252 4.09/3.98727.86/3.9827.96/AB BC AB BC g g kN m g kN m g kN m===⨯==⨯= 作用于中间层框架梁上的线荷载为：1122 4.09 5.24(40.6)21.91/4.09/3.24722.68/3.242 6.48/AB BC AB BC g kN m g kN mg kN m g kN m=+⨯-===⨯==⨯=屋面框架节点集中荷载标准值为：顶层边节点集中荷载：55174.098(887) 3.982295.41A D G G kN==⨯+⨯+-⨯⨯= 顶层中节点集中荷载：55124.098(882) 3.982217(887) 3.98123.2722B C G G kN ==⨯+⨯+-⨯⨯+⨯+-⨯⨯=楼面框架节点集中荷载标准值(包括钢窗自重及门重)为：中间层边节点集中荷载：[]174.098(887) 3.24222.432 5.24(40.6)(80.5)23 1.86.614192.08A D G G kN==⨯+⨯+-⨯⨯+⨯+⨯-⨯--⨯⨯+⨯= 中间层中节点集中荷载：[]174.098(887) 3.242212(882) 3.24225.24(40.6)(80.5) 2.1 1.00.2 2.1 1.06.614255.91A D G G kN==⨯+⨯+-⨯⨯+⨯+-⨯⨯+⨯-⨯--⨯+⨯⨯+⨯= 活荷载顶层节点：5555555 1.5710.5/1.523/17(887) 1.523.632217(887) 1.52212(882) 1.534.1322AB CD BC A D B C p p kN m p kN mP P kN P P kN ==⨯==⨯===⨯+-⨯⨯===⨯+-⨯⨯+⨯+-⨯⨯=中间层节点： 1.5710.5/2.024/17(887) 1.523.632217(887) 1.52212(882)237.6322AB CD BC A D B C p p kN m p kN mP P kN P P kN ==⨯==⨯===⨯+-⨯⨯===⨯+-⨯⨯+⨯+-⨯⨯=荷载示意图用结构力学求解器求解的恒荷载作用下内力图如下：荷载图弯矩图剪力图轴力图轴力和剪力单位为kN；弯矩单位为kN m 。

大连理工14春《工程抗震》大作业答案

大工14春《工程抗震》大作业及要求注意：从以下五个题目中任选两个进行解答（注意：从题目一、二中选择一道计算题，并从题目三、四、五中选择一道问答题，分别进行解答，不可同时选择两道计算题或者问答题）；解答前，需将所选题目复制（使老师明确你所选的题目）。

题目一：底部剪力法。

结构阻尼比0.05ξ=，自振周期为10.55s T =，Ⅰ1类场地类别，设计地震分组为第一组，抗震设防烈度为8度（设计基本地震加速度为0.30g ）。

按底部剪力法计算结构在多遇地震时的水平地震作用及地震剪力。

3580kN=2550kN =1500kN=（a ）计算简图4600kN =5450kN=题目二：振型分解反应谱法。

某3层钢筋混凝土框架结构，阻尼比0.05ξ=，部分尺寸如下图所示。

各楼层重力荷载代表值为11200kNG=，21000kNG=，3650kNG=，场地土Ⅱ类，设防烈度8度（设计基本地震加速度为0.20g），地震分组在第二组。

现算得前3个振型的自振周期为10.68sT=，20.24sT=，30.16sT=，振型如下图所示。

试用振型分解反应谱法求该框架结构的层间地震剪力标准值。

题目四：什么叫鞭端效应？设计时如何考虑这种效应？答：地震作用下突出建筑物屋面的附属小建筑物，如电梯间、女儿墙、附墙烟囱等由于重量和刚度突然变小，高振型影响较大，震害较为严重，称为鞭端效应。

结构按底部剪力法计算时，只考虑了第一振型的影响，突出屋出的小建筑物在地震中相当于受到从屋面传来的放大了的地面加速度，采用基底剪力法计算这类小建筑的地震作用效应时应乘以放大系数3。

抗震大作业

目录作业题目........................................................................................................................................... 1 1、求解结构的固有频率和固有阵型 (2)1.1. 单元离散及节点编号 ................................................................................................... 2 1.2. 形成整体刚度矩阵 ....................................................................................................... 3 1.3. 形成整体阻尼矩阵 ....................................................................................................... 4 1.4. 求固有频率及固有振型 ............................................................................................... 4 1.5. 求节点等效荷载 ........................................................................................................... 5 2、求解地震加速度作用下节点5的响应位移时程曲线 . (6)2.1 中心差分法 (6)2.2 Newmark 法（14β=）............................................................................................ 7 2.3 Wilson θ-法（ 1.42θ=） .......................................................................................... 9 2.4 三种计算方法结果对比如下: ......................................................................................... 10 2.5 使用Mdias 软件计算 ...................................................................................................... 11 附录： .. (14)作业题目：平面梁（每节点3个自由度）{}1122Y y y θθT=[]2232212612664621261266264L L L L L L EI K L L L L L L L --=----⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦ []2222156225413224133541315622420133224L L LL LL mL M L L L LL L --=----⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦0.0480.045C M K =+ 2A=125cm 4I 4150cm =2E 2100kN/cm = m 1.0e-4KN/cm =请采用自编程序和通用程序进行对比求：1、结构的固有频率和振型2、图示地震加速度作用下节点响应位移时程曲线。

地震工程学大作业

目录一、题目 .............................................................................................................. 错误!未定义书签。

二、解答 .............................................................................................................. 错误!未定义书签。

框架设计 ......................................................................................................... 错误!未定义书签。

钢框架内力计算........................................................................................ 错误!未定义书签。

钢梁验算 ................................................................................................. 错误!未定义书签。

钢柱验算.................................................................................................... 错误!未定义书签。

钢框架杆件模型弹性分析................................................................................ 错误!未定义书签。

钢框架弹性阶段分析 ............................................................................. 错误!未定义书签。

地震工程学大作业代码(求ELCENTRO波的位移

M=10^5*diag([3.4 3.4 3.2 3.2 3.0 2.8 2.7 2.6]); %输入质量矩阵K=10^8*(diag([4.4 4.2 4 3.8 3.6 3.4 3.2 1.6])+...diag([-2.2 -2.0 -2.0 -1.8 -1.8 -1.6 -1.6],1)+...diag([-2.2 -2.0 -2.0 -1.8 -1.8 -1.6 -1.6],-1)); %刚度矩阵[fai,w]=eig(inv(K)*M); %特征值问题求解omiga=1./diag(w.^0.5); %得到周期for i=1:8 %归一化振型[val,poi]=max(abs(fai(:,i)));fai(:,i)=fai(:,i)/fai(poi,i);endgama=fai'*M*ones(8,1)./diag(fai'*M*fai); %计算振型参与系数kesai=(0.1347+0.006306*omiga.^2)/2./omiga; %计算各阶阻尼file_id='D:\地震动响应图\';mkdir(file_id); %创建文件夹保存生成的所有图A= textread('NS.txt','%n'); %读取地震波数据N=length(A);mint=0.01; %求时间间隔T=(N-1)*mint;B=xlsread('T');allmotion=zeros(N,8); %用来存储所有时刻的位移allacceleration=zeros(N,8); %用来存储所有时刻的加速度值allacceleration(1,:)=A(1); %对0时刻加速度矩阵初始化allvelocity=zeros(N,8); %用于存储所有时刻速度vg(1)=0;xg(1)=0;for m=1:N-1 %计算场地速度、位移用于对加速度进行基线修正vg(m+1)=vg(m)+A(m)*mint+(A(m+1)-A(m))*mint/2;xg(m+1)=xg(m)+vg(m)*mint+A(m)*mint^2/2+mint^2*(A(m+1)-A(m))/6;endc1=28/13/T^2*(2*vg(N)-15/T^5*xg*(3*T*B(:,1).^2-2*B(:,1).^3)*mint); %计算系数c1c0=(vg(N)-c1*T^2/2)/T; %计算系数c0A=max(A)/max(abs(A-(c0+c1*B(:,1))))*(A-(c0+c1*B(:,1))); %修正后的地震波加速度xg=xg'-(0.5*c0*B(:,1).^2+1/6*c1*B(:,1).^3); %修正后的场地位移for i=1:8 %第三层循环用于对不同周期值计算反应S=1+kesai(i)*omiga(i)*mint+((omiga(i)*mint)^2)/6;for j=2:N %对每个周期利用NewMark方法进行数值积分得到对应值（具体请查看word文档）Q=2*omiga(i)*kesai(i)*allacceleration(j-1,i)*mint+...(omiga(i)^2)*(allvelocity(j-1,i)*mint+0.5*allacceleration(j-1,i)*mint^2);deltaa=-(gama(i)*(A(j)-A(j-1))+Q)/S;allacceleration(j,i)=allacceleration(j-1,i)+deltaa;allvelocity(j,i)= allvelocity(j-1,i)+allacceleration(j-1,i)*mint+deltaa*mint/2;allmotion(j,i)= allmotion(j-1,i)+allvelocity(j-1,i)*mint+0.5*allacceleration(j-1,i)*mint^2+...deltaa*(mint^2)/6;end %积分计算反应并存放在前面的空矩阵里endmotion1=allmotion*fai'; %计算位移acceleration=allacceleration*fai'; %计算加速度velocity=allvelocity*fai'; %计算速度for n=1:9figure('color','white') %利用figure函数准确的控制画图窗口if n~=9axis([0,50,-0.03,0.03]) %控制坐标轴范围，便于比较endgrid onhold onbox off%annotation('arrow',[0.132 0.132],[0.8 1]);%annotation('arrow',[0.8 1],[0.108 0.108]); %产生坐标轴箭头（可能不够美观）if n==1 %计算相关数据和得到图名name='底层层间位移时程';motion(:,n)=motion1(:,n);elseif n==9name='顶层位移时程';motion(:,9)=motion1(:,8)+xg;elsemotion(:,n)=motion1(:,n)-motion1(:,n-1);name=strcat('第',num2str(n-1),'-',num2str(n),'层间位移时程');endplot(B(:,1),motion(:,n),'linewidth',2) %画相对位移时程xlabel('时间（s）','FontName','宋体','FontSize',16); %x轴标注ylabel('位移(m)','FontName','宋体','FontSize',16); %y轴标注title(name,'FontName','宋体','FontSize',20) %标注图名set(gcf,'position',get(0,'screensize')); %图形全屏，便于查看shgF=getframe(gcf);imwrite(F.cdata,[file_id,strcat(name,'.png')]); %存储得到的时程图endGK=repmat(-[-2 -2 -1.8 -1.8 -1.8 -1.8 -1.6 -1.6]*10^8,N,1);%输入刚度矩阵motion(:,8)=motion1(:,8)-motion1(:,7); %计算顶层的层间相对位移Fq=motion(:,1:8).*GK; %计算剪力figure('color','white') %生成图形窗口，背景白色plot(B(:,1),Fq(:,1)) %画底层剪力时程图xlabel('时间（s）','FontName','宋体','FontSize',16); %x轴标注ylabel('剪力（N）','FontName','宋体','FontSize',16); %y轴标注title('底层剪力时程','FontName','宋体','FontSize',20); %标注图名set(gcf,'position',get(0,'screensize')); %图形全屏，便于查看F=getframe(gcf); %存储生成的图形imwrite(F.cdata,[file_id,strcat('底层剪力时程图','.png')]);maxFq=zeros(1,16);for f=1:8 %计算所需剪力和位移包络图数据，并得到便于画图的矩阵maxmotion(2*f-1)=max(abs(motion(:,f)));maxmotion(2*f)=maxmotion(:,2*f-1);maxFq(2*f-1)=max(abs(Fq(:,f)));maxFq(2*f)=max(abs(Fq(:,f)));endfigure('color','white')plot(maxmotion,[0 0.5 0.5 1 1 1.5 1.5 2 2 2.5 2.5 3 3 3.5 3.5 4],'k-');%生成层间位移包络图title('层间位移绝对值最大值包络图','FontName','宋体','FontSize',20); %标注图名text(maxmotion(1:2:16),0.25:0.5:4,num2str(maxmotion(1:2:16)')) %标注位移数据set(gcf,'position',get(0,'screensize')); %图形全屏，便于查看F=getframe(gcf);imwrite(F.cdata,[file_id,strcat('层间位移绝对值最大值包络图','.png')]); %存储图形figure('color','white')plot(maxFq/1000,[0 0.5 0.5 1 1 1.5 1.5 2 2 2.5 2.5 3 3 3.5 3.5 4],'k-');title('层间剪力绝对值最大值包络图','FontName','宋体','FontSize',20); %标注图名text(maxFq(1:2:16)/1000,0.25:0.5:4,num2str(0.001*maxFq(1,1:2:16)'),'FontSize',10);set(gcf,'position',get(0,'screensize')); %图形全屏，便于查看F=getframe(gcf);imwrite(F.cdata,[file_id,strcat('层间剪力绝对值最大值包络图','.png')]);。

地震工程作业

作业1绘制1940 El Centro,N—S分量地震动的绝对加速度、相对速度和相对位移反应谱。

地震动：在PEER Ground Motion Database自行下载经典的1940 El Centro，N—S分量。

要求：在此模板内完成,A4纸打印.自编程序与软件（Bispec或Seismosigna等）计算反应谱进行对比。

提交自编写程序.Matlab程序：clearfid = fopen（'E:\Earthquake\El centro。

txt’)；[Accelerate,count］= fscanf(fid,’%g’)；％count 读入的记录的量Accelerate=9.8＊Accelerate'；％单位统一为m和stime=0：0。

02:(count-1)＊0。

02; ％单位sDisplace=zeros（1，count); %相对位移Velocity=zeros（1,count）；％相对速度AbsAcce=zeros(1，count); ％绝对加速度DampA=[0.00,0。

02，0。

05]; ％三个阻尼比TA=0。

0：0.02：4；％TA=0。

000001:0。

02:4；结构周期Dt=0。

02；%地震记录的步长％记录计算得到的反应，MDis为某阻尼时最大相对位移%MVel为某阻尼时最大相对速度，MAcc某阻尼时最大绝对加速度MDis=zeros(3，length（TA））;MVel=zeros（3，length（TA））；MAcc=zeros（3，length(TA）);j=1;for Damp=［0.00,0。

05,0。

1］t=1；for T=0.0:0。

02:4Frcy=2*pi/T ;DamFrcy=Frcy＊sqrt（1-Damp＊Damp）；e_t=exp(-Damp＊Frcy＊Dt);s=sin（DamFrcy＊Dt）;c=cos(DamFrcy*Dt）;A=zeros(2，2);A(1,1）=e_t＊（s＊Damp/sqrt(1-Damp*Damp)+c);A(1，2)=e_t＊s/DamFrcy；A（2,1）=—Frcy*e_t*s/sqrt(1—Damp*Damp）；A（2,2)=e_t*（—s*Damp/sqrt(1—Damp＊Damp）+c）;d_f=(2*Damp^2-1）/(Frcy^2*Dt）；d_3t=Damp/（Frcy^3＊Dt）；B=zeros（2，2）；B(1,1）=e_t*(（d_f+Damp/Frcy)＊s/DamFrcy+（2*d_3t+1/Frcy^2）＊c)-2*d_3t；B(1，2）=—e_t＊（d_f＊s/DamFrcy+2＊d_3t＊c)-1/Frcy^2+2＊d_3t；B（2，1）=e_t*（(d_f+Damp/Frcy)*（c-Damp/sqrt(1-Damp^2）*s)—（2＊d_3t+1/Frcy^2)*(DamFrcy＊s+Damp＊Frcy*c)）+1/（Frcy^2*Dt);B（2，2）=e_t＊（1/（Frcy^2*Dt）*c+s＊Damp/（Frcy＊DamFrcy＊Dt））—1/(Frcy^2*Dt）;for i=1：(count—1)Displace（i+1)=A(1,1)＊Displace(i)+A(1,2）＊Velocity(i)+B（1,1)＊Accelerate（i）+B（1,2)＊Accelerate（i+1）；Velocity（i+1)=A(2，1）*Displace（i）+A（2,2）＊Velocity（i)+B（2,1）*Accelerate(i）+B(2，2)*Accelerate(i+1)；AbsAcce(i+1)=-2＊Damp＊Frcy*Velocity（i+1)—Frcy^2*Displace(i+1）;endMDis（j,t）=max(abs（Displace）)；MV el(j,t）=max(abs(Velocity))；if T==0。

地震工程作业

作业1绘制1940 El Centro，N-S分量地震动的绝对加速度、相对速度和相对位移反应谱。

地震动：在PEER Ground Motion Database自行下载经典的1940 El Centro，N-S分量。

要求：在此模板内完成，A4纸打印。

自编程序与软件（Bispec或Seismosigna等）计算反应谱进行对比。

提交自编写程序。

Matlab程序：clearfid = fopen('E:\Earthquake\El centro.txt');[Accelerate,count] = fscanf(fid,'%g'); %count 读入的记录的量Accelerate=9.8*Accelerate'; %单位统一为m和stime=0:0.02:(count-1)*0.02; %单位sDisplace=zeros(1,count); %相对位移Velocity=zeros(1,count); %相对速度AbsAcce=zeros(1,count); %绝对加速度DampA=[0.00,0.02,0.05]; %三个阻尼比TA=0.0:0.02:4; %TA=0.000001:0.02:4;结构周期Dt=0.02; %地震记录的步长%记录计算得到的反应，MDis为某阻尼时最大相对位移%MVel为某阻尼时最大相对速度，MAcc某阻尼时最大绝对加速度MDis=zeros(3,length(TA));MVel=zeros(3,length(TA));MAcc=zeros(3,length(TA));j=1;for Damp=[0.00,0.05,0.1]t=1;for T=0.0:0.02:4Frcy=2*pi/T ;DamFrcy=Frcy*sqrt(1-Damp*Damp);e_t=exp(-Damp*Frcy*Dt);s=sin(DamFrcy*Dt);c=cos(DamFrcy*Dt);A=zeros(2,2);A(1,1)=e_t*(s*Damp/sqrt(1-Damp*Damp)+c);A(1,2)=e_t*s/DamFrcy;A(2,1)=-Frcy*e_t*s/sqrt(1-Damp*Damp);A(2,2)=e_t*(-s*Damp/sqrt(1-Damp*Damp)+c);d_f=(2*Damp^2-1)/(Frcy^2*Dt);d_3t=Damp/(Frcy^3*Dt);B=zeros(2,2);B(1,1)=e_t*((d_f+Damp/Frcy)*s/DamFrcy+(2*d_3t+1/Frcy^2)*c)-2*d_3t;B(1,2)=-e_t*(d_f*s/DamFrcy+2*d_3t*c)-1/Frcy^2+2*d_3t;B(2,1)=e_t*((d_f+Damp/Frcy)*(c-Damp/sqrt(1-Damp^2)*s)-(2*d_3t+1/Frcy^2)* (DamFrcy*s+Damp*Frcy*c))+1/(Frcy^2*Dt);B(2,2)=e_t*(1/(Frcy^2*Dt)*c+s*Damp/(Frcy*DamFrcy*Dt))-1/(Frcy^2*Dt); for i=1:(count-1)Displace(i+1)=A(1,1)*Displace(i)+A(1,2)*Velocity(i)+B(1,1)*Accelerate(i)+B(1,2)* Accelerate(i+1);Velocity(i+1)=A(2,1)*Displace(i)+A(2,2)*Velocity(i)+B(2,1)*Accelerate(i)+B(2,2)* Accelerate(i+1);AbsAcce(i+1)=-2*Damp*Frcy*Velocity(i+1)-Frcy^2*Displace(i+1);endMDis(j,t)=max(abs(Displace));MVel(j,t)=max(abs(Velocity));if T==0.0MAcc(j,t)=max(abs(Accelerate));elseMAcc(j,t)=max(abs(AbsAcce));endDisplace=zeros(1,count);Velocity=zeros(1,count);AbsAcce=zeros(1,count);t=t+1;endj=j+1;endclose allfigure %绘制位移反应谱plot(TA,MDis(1,:),'-b',TA,MDis(2,:),'-r',TA,MDis(3,:),':k')title('Displacement')xlabel('Tn(s)')ylabel('Displacement(m)')legend('ζ=0','ζ=0.02','ζ=0.05')gridfigure %绘制速度反应谱plot(TA,MVel(1,:),'-b',TA,MVel(2,:),'-r',TA,MVel(3,:),':k') title('Velocity')xlabel('Tn(s)')ylabel('velocity(m/s)')legend('ζ=0','ζ=0.02','ζ=0.05')gridfigure %绘制绝对加速度反应谱plot(TA,MAcc(1,:),'-b',TA,MAcc(2,:),'-r',TA,MAcc(3,:),':k') title('Absolute Acceleration')xlabel('Tn(s)')ylabel('absolute acceleration(m/s^2)')legend('ζ=0','ζ=0.02','ζ=0.05')grid由此得到下图反应谱与由Seismosigna 得到的反应谱作出比较相对位移反应谱（上）与Seismosigna 计算相对位移反应谱（下）的比较00.10.20.30.40.50.60.7DisplacementT n(s)D is p la c e m e n t (m )ζ=0ζ=0.02ζ=0.05相对速度反应谱（上）与Seismosigna 计算相对速度反应谱（下）的比较0.511.522.5VelocityT n(s)velocity(m/s)ζ=0ζ=0.02ζ=0.05绝对加速度反应谱（左图）与Seismosigna计算绝对加反应谱（右图）的比较5101520253035Absolute AccelerationT n(s)absoluteacceleration(m/s2)ζ=0ζ=0.02ζ=0.05作业2、使用中心差分方法和Newmark 方法分别计算下面单自由度结构的地震反应。