地震工程学大作业.pdf

《地震工程学大作业.pdf》由会员分享，可在线阅读，更多相关《地震工程学大作业.pdf（26页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、目录目录一、二、题目.错误错误!未定义书签。未定义书签。解答.错误错误!未定义书签。未定义书签。框架设计.错误错误!未定义书签。未定义书签。钢框架内力计算.错误错误!未定义书签。未定义书签。钢梁验算.错误错误!未定义书签。未定义书签。钢柱验算.错误错误!未定义书签。未定义书签。钢框架杆件模型弹性分析.错误错误!未定义书签。未定义书签。钢框架弹性阶段分析.错误错误!未定义书签。未定义书签。简谐波下的结果分析.错误错误!未定义书签。未定义书签。钢框架弹塑性阶段分析.错误错误!未定义书签。未定义书签。弹塑性分析介绍.错误错误!未定义书签。未定义书签。恢复力模型介绍.错误错误!未定义书签。未定义书签。

2、克拉夫（Clough）退化双线模型.错误错误!未定义书签。未定义书签。纤维模型理论.错误错误!未定义书签。未定义书签。克拉夫模型地震波反应计算.错误错误!未定义书签。未定义书签。纤维模型弹塑性地震反应计算.错误错误!未定义书签。未定义书签。参考文献.错误错误!未定义书签。未定义书签。一、一、题目题目1.1.设计并确定结构模型：设计并确定结构模型：如图所示平面框架，跨度 L=6m，层高 H=，榀距 4m。请首先设计确定框架杆件截面，需考虑楼面恒荷载和楼面活荷载（活荷载可为 m2）。2.2.弹性阶段：观察共振现象弹性阶段：观察共振现象首先计算自振频率，然后采用简谐波进行弹性阶段计算。（输入波的幅值

3、保持不变，改变频率与自振频率的比值分别为：，）3.3.计算弹塑性地震反应计算弹塑性地震反应以 El-centro N-S 波为基本波，改变幅值，从弹性阶段计算至弹塑性阶段，直至破坏。分别采用以下两种模型，并比较两种模型的差异：1）杆件模型：恢复力模型，如混凝土结构采用武田三线性模型，钢结构可采用Clough 模型；2）截面模型：采用纤维模型。二、二、解答解答框架设计框架设计楼板采用混凝土板，板厚为 110mm。荷载计算：楼面荷载：30 厚水泥砂浆结合层：0.0320 0.60kN/m2110 厚现浇混凝土楼板：0.1125 2.8kN/m215 厚混合砂浆天棚抹灰：0.01517 0.26kN

4、/m2合计：3.66kN/m2本框架中采用现浇混凝土楼板，可视为刚性铺面，能够阻止主梁上翼缘的侧向失稳，而不需考虑其整体稳定，只需满足其强度、刚度和局部稳定要求。本框架的梁选用宽翼缘 H 型钢梁。框架梁的跨度为 6m，按高跨比1/101/20，选取HN2005001016。框架梁自重为0.77kN/m。相应的梁柱截面参数见表 1。表 梁柱构件截面特性构件梁柱截面尺寸HN5002001016HN400200813A(cm2)Ix(cm4)4780023700Wx(cm3)19101190自重(kN/m)作用在梁上的恒荷载为：恒荷载：g 3.6640.896 17.536kN/m活荷载：q 2.5

5、4 10kN/m作用在柱上的恒荷载：g1 0.660kN/mgg1图作用在框架上的恒载钢框架内力计算钢框架内力计算利用 MIDAS 计算平面框架的内力，得到其内力图如图所示。钢梁验算钢梁验算抗弯强度验算：抗弯强度验算：a)弯矩图（单位：kN*m）b)剪力图（单位：kN）-69c)轴力图（单位：kN）图 钢框架内力图满足要求抗剪强度验算：抗剪强度验算：满足要求刚度验算：刚度验算：均布荷载标准值g 3.6640.896 17.536kN/m跨中挠度与跨度的比值满足要求整体稳定验算：整体稳定验算：设主梁跨中位置有侧向支撑，主梁受压翼缘的自由长度L1等于 3m，L1与梁受压翼缘宽度 b1之比为因此，梁

6、的整体稳定性满足要求。钢柱验算钢柱验算强度验算：强度验算：惯性矩回转半径长细比查表得，则满足要求局部稳定验算：局部稳定验算：翼缘腹板满足要求刚度验算：刚度验算：满足要求钢框架杆件模型弹性分析钢框架杆件模型弹性分析钢框架弹性阶段分析钢框架弹性阶段分析首先计算自振频率，然后采用简谐波进行弹性阶段计算。（输入波的幅值保持不变，改变频率与自振频率的比值分别为：，）。由于是水平抗震分析，故在用 Midas-Gen 进行自振周期计算时仅需考虑 X 方向（结构分析类型为 X-Z 平面），将恒载+活载以及自重转化为质量，且将自重转化为集中质量，体系只有两个自由度，采用子空间迭代法，迭代振型为两阶，分析后表明第

7、一阶的振型参与质量就达到了 100%，故仅需考虑第一阶振型即该框架变成了单自由度体系，图 4 为该框架的振型图。弹塑性时程分析中静力荷载取+作用在结构上，在 midas-Gen 分析中将此荷载转化为质量，通过特征值分析，结构自振频率与周期见表，其中 1、2 阶模态如图所示。表 结构自振频率与周期阶数自振频率(Hz)自振周期(s)1 阶2 阶a)第一阶模态b)第二阶模态图结构模态图结构阻尼矩阵采用瑞雷阻尼假定，1、2 阶阻尼比均取 2%。图 倍基频度响应时程谐振力下左柱顶点加速图 倍基频谐振力下左柱顶点加速度响应时程图基频谐振力下左柱顶点加速度响应时程图 倍基频谐振力下左柱顶点加速度响应时程图2

8、 倍基频谐振力下左柱顶点加速度响应时程简谐波下的结果分析简谐波下的结果分析共振下结构仍处于弹性状态，从而保证了简谐波作用下的分析为弹性分析。由于均采用 Newmark时程分析方法，杆模型与纤维模型计算结果一致。从而可以看出弹性状态下，虽然纤维模型相比杆模型截面划分更精细，但对计算结果精度提高很小。观察共振现象，由于结构各阶频率相距较远，取激励频率在结构一阶自振频率附近变化观察共振反应。当激励频率小于结构一阶自振频率时，在阻尼作用下结构很快达到稳态反应；当激励频率等于结构一阶自振频率时，结构出现共振反应，结构位移逐渐增大，由于阻尼的存在，最后达到稳态反应。当激励频率大于结构一阶自振频率时，结构位

9、移反应幅值开始阶段较大，并逐渐减小达到稳态反应。在保持相同激励幅值时，结构共振下的位移反应相比非共振下大的多。钢框架弹塑性阶段分析钢框架弹塑性阶段分析弹塑性分析介绍弹塑性分析介绍结构弹塑性动力分析的基本动力方程为：&MUCUF MUg（）其中，F为非线性恢复力向量。上式的增量方程为：&MUjCUjKUj MUgj（）式（）可用动力分析的逐步积分法求解。全量形式的动力方程为FKpU。在采用数值分析技术的前提下，结构线性地震反应分析与非线性地震反应分析的主要差别在于刚度矩阵是否可变。对于弹塑性结构，在每一步增量反应计算之前，要先修正矩阵K或Kp中各元素的量值，即所谓的刚度修正技术。刚度修正过程实质

10、上是一个重新形成总刚度矩阵的过程。恢复力模型介绍恢复力模型介绍恢复力模型是根据大量的从试验中获得的恢复力与变形关系曲线经适当抽象和简化而得到的实用数学模型，是构件的抗震性能在结构弹塑性地震反应分析中具体体现。若仅用静力非线性分析，模型一般是指力与变形关系骨架曲线的数学模型；而如果是用于结构动力非线性时程分析，恢复力模型不仅包含骨架曲线，同时也包含各阶段滞回环的数学模型。常见的恢复力模型有兰伯格奥斯古德模型、克拉夫（Clough）模型和武田模型。兰伯格奥斯古德模型曾被广泛用于土体和包括钢结构在内的各种结构物的非线性反应分析，偶尔也用于钢筋混凝土弯曲构件。克拉夫模型主要是针对钢筋混凝土受弯构件的恢

11、复力特性提出的，而武田模型是从较多的钢筋混凝土构架试验所得的恢复力特性曲线抽象得出的，适用于以弯曲破坏为主的情况。钢筋混凝土结构构件的恢复力模型一般分为曲线型和折线型两种，其中曲线型比较接近结构的实际受力特性，结果比较精确，但是刚度计算比较复杂，因此，应用很少；折线型恢复力模型由若干直线段所构成，刚度变化不连续，存在拐点问题，但刚度计算比较简单，故在实际工程中得到广泛应用。本次分析中，将采用克拉夫（Clough）退化双线模型作为钢结构的恢复力模型。克拉夫（克拉夫（CloughClough）退化双线模型）退化双线模型初次加载时沿着双折线骨架曲线移动，屈服后卸载路径沿着退化后的斜率移动；当反向加载

12、时，指向反向最大变形点；反向没有发生屈服时，屈服点为最大变形点。克拉夫模型中认为全截面处于开裂状态，截面的刚度由受拉钢筋的受弯屈服状态决定。对正向和负向可定义不同的屈服后的刚度折减系数，适用于梁、柱、支撑构件。P PK2K2P1P1(+)(+)(+)(+)D D()max(),P PmaxK K0 0KrKr(-)(-)D1D1(-)(-)D1D1(+)(+)KrKr(+)(+)D DD D()max(),P PmaxK2K2(-)(-)P1P1(-)(-)P PP1P1(+)(+)D1D1(-)(-)D1D1(+)(+)D DP1P1(-)(-)图克拉夫（Clough）退化双线模型克拉夫模型

13、的骨架曲线由下列参数决定：:Pl()、Pl()正向和负向的第一屈服强度；Dl()、Dl()正向和负向的第一屈服强度；Dl()、Dl()正向和负向的第一屈服强度；K0初始刚度；K2()、K2()正向和负向的第二条折线刚度，K2()1()K0、K2()1()K0；1()、1()正向和负向的第一屈服后刚度折减系数；Kr()、Kr()正向和负向卸载时的刚度：Kr()D1()K0()Dmax K0Kr()，D1()K0()Dmax K0其中Dmax()、Dmax()为正向和负向的最大变形，没有屈服的区段使用屈服变形；为计算卸载刚度的幂阶。克拉夫模型的路径移动规则：（1）Dmax Dl时沿斜率为K0的直线

14、移动。（2）变形D第一次超过Dl()时或者超过当前的最大变形时，沿着斜率为K2()、K2()的第二折线移动。（3）在Dl()D、D Dl()状态下卸载时，沿着卸载刚度Kr()、Kr()的斜率移动。（4）卸载过程中荷载的符号发生变化时，将沿着指向反向最大变形点的直线移动（如反向未屈服则指向反向屈服点）。Clough 模型一般只适用于具有梭形滞回曲线的单纯受弯构件。纤维模型理论纤维模型理论纤维模型是将构件单元分割为许多沿构件纵向的纤维的模型。纤维模型不仅可以准确模拟受弯构件的力学特性，而且可以考虑截面内纤维的局部损伤状态。使用纤维模型时可利用纤维材料的应力应变关系和截面应变的分布形状假定较为准确地

15、确定截面的弯矩曲率关系，特别是可以考虑轴力引起的中和轴变化。另外纤维模型同样可以考虑轴力和弯矩、两个弯矩之间的相互影响，但是因为不能反映剪切破坏，所以一般用于剪切变形不大的线单元。本次分析采用的分析软件 MIDAS/GEN 中的纤维模型使用了一下几个假定：1.截面的变形维持平截面与构件轴线垂直；2.不考虑钢筋与混凝土之间的滑移；（本次分析采用钢结构）3.杆件单元截面形心的连线为直线。（1）纤维模型的计算过程首先假设在各单元的积分点上存在用纤维模型定义的截面，积分点数量最多不超过 20 个，积分方法使用高斯-罗贝托方法，只有当积分点数量为 2 时才使用古典高斯方法。然后将前次时间步中计算的两端的

16、构件内力通过转换并排除刚体运动成分后获得构件两端的的广义单元内力(轴力和弯矩)。最后通过内力的内插函数计算单元内各位置的内力。图构件任意位置的构件内力和变形单元荷载向量：Q Q单元变形向量：q qQ1,Q2,Q3,Q4,Q5T Tq1,q2,q3,q4,q5T TT T截面内力向量：D D(x)Mz(x)My(x)N(x)截面变形向量：d d(x)z(x)y(x)(x)D Di i(x)b b(x)Q Qi iT T其中，b b(x)为内力内插函数(Force Interpolation Function)，公式如下。x x-1000 LL x xb b(x)00-1 0L L00001其中，

17、L为单元长度。通过截面的柔度可计算截面的变形，通过各截面位置的轴向、弯曲变形可计算各纤维的轴向应变。ECS z-axisECS z-axisECS y-axisECS x-axisyii-th fiberziECS y-axisx图纤维模型的截面分割示意图y(x)1z(x)(x)xizi-yi其中，：截面的位置y(x)：在 x 处截面的绕单元坐标系 y 轴的曲率z(x)：在 x 处截面的绕单元坐标系 z 轴的曲率x(x)：在 x 处截面的轴向应变yi：截面上第 i 个纤维的 y 轴位置zi：截面上第 i 个纤维的 z 轴位置i：第 i 个纤维的应变各纤维的轴向应变i对应的纤维的应力和纤维的切线

18、刚度可通过纤维材料的本构关系(constitutiverelation)计算，并由此判断纤维的状态。将一个截面内所有纤维的应力进行积分可获得截面的轴力和弯矩，对各纤维的切线刚度进行积分可获得截面的柔度矩阵，对单元内所有积分点上的截面的柔度进行积分可得单元的柔度。n(x)EjiAiy2ii=1n(x)j jk k(x)-EjiAiyizii=1n(x)-EjiAiyii=1n(x)-EjiAiyiziAiz2iAizii=1n(x)Ei=1n(x)i=1jiEji-Aiyii=1n(x)jEiAizii=1n(x)jEiAii=1n(x)Ejif fj j(x)k kj j(x)1k kj j(

19、x):j 步骤时位置 x 处截面的切线刚度f fj j(x):j 步骤时位置 x 处截面的柔度n(x):位置 x 处截面内总纤维数量E Ej ji i(x):j 步骤时位置 x 处截面内第 i 个纤维的切线刚度A Ai i:第 i 个纤维的截面面积y yi i,z zi i:第 i 个纤维在截面内的位置计算截面不平衡力所需的截面内力按下式计算。n(x)j jD DR R(x)-jiAiyii=1 A zjii ii=1n(x)n(x)jiAii=1T T在计算过程中使用牛顿-拉普森迭代方法计算至满足收敛条件。纤维模型中单元的非线性特性表现在纤维的非线性应力-应变关系(材料本构关系)上。（2）采

20、用的钢材纤维本构模型MIDAS 软件中的钢材纤维本构模型：修正梅内戈托与平托模型（Modified Menegotto&Pinto Steel Model）本构模型形状为逐渐逼近按随动硬化(Kinematic Hardening)准则定义的双折线的曲线，即本构模型在卸载路径和应变硬化区段之间的两个渐近线之间的转换区段为曲线。两个渐近线的交点和卸载方向的最大变形点之间的距离越远，转换区段的形状越光滑。可利用这样的特性模拟包辛格效应。b(1b)1R1/R 其中，rra ,R R010r0ra2:钢纤维的应变:钢纤维的应力(r,r):卸载点位置，在弹性状态时假定为(0,0)(0,0):定义当前加载或

21、卸载路径的两个渐近线的交点b:刚度折减率R0,a1,a2:常数(决定曲线状态的参数，程序使用了试验获得的优化值):在加载/卸载的方向上最大应变与0的差(绝对值)，最大应变的初始值假定与(Fy/E)相同(参见图2(0,0)2FybEE(r,r)1(r,r)2(0,0)11图钢材纤维本构模型双折线钢材纤维模型常用的双折线本构模型，屈服前加载和卸载时使用弹性刚度，屈服后加载时使用屈服后刚度，屈服后卸载、再加载时使用弹性刚度。三折线钢材纤维模型常用的三折线本构模型，可定义第一屈服和第二屈服刚度折减率，受拉区和受压区的刚度折减率可以不同。屈服前的加载和卸载使用弹性刚度，屈服后加载时使用折减的刚度，屈服后

22、的卸载和再加载使用弹性刚度。折线形状也可以通过输入应力-应变定义。非对称双折线钢材纤维模型该本构模型可以模拟钢筋的非线性特性，在受拉区可考虑屈服和拉断，在受压区可考虑屈服和失稳后的破坏。克拉夫模型弹塑性地震波反应计算克拉夫模型弹塑性地震波反应计算不同幅值的 El-centro N-S 波下的柱底弯矩-曲率关系图和柱顶位移时程图如下所示：图 幅值为的 El-centro N-S波下的柱底弯矩-曲率关系图图 幅值为的 El-centro N-S波下的柱顶位移时程图图 幅值为的 El-centro N-S波下的柱底弯矩-曲率关系图图幅值为的 El-centro N-S波下的柱顶位移时程图图 幅值为

23、3g 的 El-centro N-S 波下的柱底弯矩-曲率关系图图幅值为 3g 的 El-centro N-S 波下的柱顶位移时程图图幅值为 5g 的 El-centro N-S 波下的柱底弯矩-曲率关系图图幅值为 3g 的 El-centro N-S 波下的柱顶位移时程图纤维模型弹塑性地震反应计算纤维模型弹塑性地震反应计算本次分析采用双折线钢材纤维模型。以El-centro N-S 波为基本波，改变幅值，从弹性阶段计算至弹塑性阶段，直至破坏。以节点 2 的 X 方向位移为时间的函数，画出的曲线如图所示。图纤维模型中梁柱纤维划分图增幅系数为时节点 2 的 X 方向位移时程曲线（纤维模型）图增幅

24、系数为时节点 2 的 X 方向位移时程曲线（纤维模型）图 2.1增幅系数为时节点 2 的 X 方向位移时程曲线（纤维模型）图 2.2增幅系数为 10 时节点 2 的 X 方向位移时程曲线（纤维模型）右柱（单元 2）的弯矩曲率曲线如图所示。图 2.3增幅系数为时的右柱的弯矩转角曲线图 4增幅系数为时右柱的弯矩转角曲线图 2.5增幅系数为时的右柱的弯矩转角曲线当增幅系数为、时，结构未出现塑性铰。增幅系数为10 时，结构进入弹塑性阶段并破坏。图 2.6增幅系数为 10 时的右柱的弯矩转角曲线图 2.7增幅系数为 10 时柱顶（节点 2）的弯矩转角曲线（纤维模型）弹塑性地震分析评价弹塑性地震分析评价（

25、1）从上述图中可以看出，无论是纤维截面模型还是杆件模型，随着地震激励的增加，结构的反应也随之增大。当地震动幅值加大时，结构位移值也增大，且反应曲线变的稀疏，反应峰值减少，峰值变的平坦。这主要是幅值大，时间间隔短，构件来不及恢复变形，又遭受一个强大脉冲，所以反应波谷变高，曲线变得稀疏。前后峰值影响的范围加大，削减了许多波峰，所以波峰变少。(2)以杆件模型为例，当 1 倍 El Centro 波 N-S 分量作用于结构时，从柱顶剪力-位移关系图中看出，柱顶剪力与柱顶位移呈线性关系，结构仍处于弹性阶段；将 El Centro 波 N-S 分量增大至 5 倍时，柱顶剪力-位移滞回曲线呈现捏拢现象，表明

26、结构已进入弹塑性阶段；当El Centro 波 N-S 分量增大至 10倍时，滞回曲线变的稀疏很多，且捏拢现象更加明显，即进入塑性反应的阶段明显增多，而柱顶位移已达 40cm，可认为结构已发生破坏。(3)将纤维截面模型与杆件模型进行对比，从图中可以看出，杆件模型分析 5 倍 El Centro 波 N-S分量反应时，柱顶剪力-位移滞回曲线呈现捏拢现象，而10 倍 El Centro 波 N-S 分量作用时，捏拢现象更加明显；而纤维截面模型分析 5 倍 El Centro 波 N-S 分量反应时，柱顶剪力-位移滞回曲线较为丰满，El Centro 波 N-S 分量增大到 10 时，较杆件模型而言

27、，滞回曲线虽呈现捏拢现象，但不明显。这主要是因为纤维截面模型与杆件模型相比，纤维截面模型分析时，构件的恢复力特性为截面上纤维本构关系的积分结果，因此纤维模型法具有较高的计算精度。当用于承受双向弯曲和变轴力耦合作用下构件的非线性分析时，纤维截面模型有着明显的优势。参考文献参考文献1 李杰，李国强.地震工程学导论M.北京：地震出版社，1992.2 胡聿贤.地震工程学（第二版）M.北京：地震出版社，2006.3 Roy，Jr 着，常岭，李振邦译.结构动力学M.人民交通出版社，1996.4 同济大学计算数学教研室.现代数值计算M.北京：人民邮电出版社，2009.5 中华人民共和国国家标准.建筑抗震设计规范（GB500011-2010）S.北京：中国建筑工业出版社，2010.6 杨志勇,黄吉锋,邵弘.弹性与弹塑性动力时程分析方法中若干问题探讨J.建筑结构学报,2009.7 裴星洙.建筑结构抗震分析与设计M.北京：北京大学出版社，2013.