地震工程学大作业

目录目录 一 题目 1 二 解答 1 2 1 框架设计 1 2 1 1 钢框架内力计算 2 2 1 2 钢梁验算 3 2 1 3 钢柱验算 4 2 2 钢框架杆件模型弹性分析 5 2 2 1 钢框架弹性阶段分析 5 2 2 2 简谐波下的结果分析 7 2 3 钢框架弹塑性阶段分析 8 2 3 1 弹塑性分析介绍 8 2 3 2 恢复力模型介绍 8 2 3 3 克拉夫 Clough 退化双线模型 9 2 3 4 纤维模型理论 10 2 3 5 克拉夫模型地震波反应计算 14 2 3 5 纤维模型弹塑性地震反应计算 17 参考文献 22 1 22 一 一 题目题目 1 设计并确定结构模型 设计并确定结构模型 如图所示平面框架 跨度 L 6m 层高 H 3 6m 榀距 4m 请首先设计确定框架杆件截面 需考虑楼面恒荷载 和楼面活荷载 活荷载可为 2 5kN m2 2 弹性阶段 观察共振现象弹性阶段 观察共振现象 首先计算自振频率 然后采用简谐波进行弹性阶段计算 输入波的幅值保持不变 改变频率与自振频率的比 值分别为 0 25 0 5 1 0 1 5 2 0 3 计算弹塑性地震反应计算弹塑性地震反应 以 El centro N S 波为基本波 改变幅值 从弹性阶段计算至弹塑性阶段 直至破坏 分别采用以下两种模型 并比较两种模型的差异 1 杆件模型 恢复力模型 如混凝土结构采用武田三线性模型 钢结构可采用 Clough 模型 2 截面模型 采用纤维模型 二 二 解答解答 2 1 框架设计框架设计 楼板采用混凝土板 板厚为 110mm 荷载计算 楼面荷载 30 厚水泥砂浆结合层 2 0 03 200 60kN m 110 厚现浇混凝土楼板 2 0 11 252 8kN m 2 22 15 厚混合砂浆天棚抹灰 2 0 015 170 26kN m 合计 2 3 66kN m 本框架中采用现浇混凝土楼板 可视为刚性铺面 能够阻止主梁上翼缘的侧向失稳 而不需考 虑其整体稳定 只需满足其强度 刚度和局部稳定要求 本框架的梁选用宽翼缘 H 型钢梁 框架梁 的跨度为 6m 按高跨比 1 10 1 20 选取 HN200 500 10 16 框架梁自重为 相应的梁柱0 77kN m 截面参数见表 1 表 2 1 梁柱构件截面特性 构件截面尺寸 A cm2 Ix cm4 Wx cm3 自重 kN m 梁 HN500 200 10 16114 24780019100 896 柱 HN400 200 8 1384 12237001190 0 660 作用在梁上的恒荷载为 恒荷载 0 893 66417 5 6kN m63g 活荷载 2 5410kN mq 作用在柱上的恒荷载 1 0 660kN mg g g 3q 3 图 2 1 作用在框架上的恒载 2 1 1 钢框架内力计算钢框架内力计算 利用 MIDAS 计算平面框架的内力 得到其内力图如图 2 2 所示 1 g 3 22 40 840 8 62 7 40 8 40 8 17 6 17 6 a 弯矩图 单位 kN m 16 2 16 2 16 2 16 2 69 0 b 剪力图 单位 kN 16 2745 69 03 69 0 c 轴力图 单位 kN 图 2 2 钢框架内力图 2 1 2 钢梁验算钢梁验算 抗弯强度验算 抗弯强度验算 62 7 106 1 05 1910 103 52 83 2 215 2 满足要求 抗剪强度验算 抗剪强度验算 69 0 69 69 090 4 22 69 0 10 3 1 05 106 4 78 108 10 25 63 2 125 2 满足要求 刚度验算 刚度验算 均布荷载标准值 0 893 66417 5 6kN m63g 跨中挠度与跨度的比值 5 3 384 5 17 536 60003 384 2 06 105 4 78 108 5 97 10 4 1 400 满足要求 整体稳定验算 整体稳定验算 设主梁跨中位置有侧向支撑 主梁受压翼缘的自由长度 L1等于 3m L1与梁受压翼缘宽度 b1之比 为 1 1 3000 200 15 16 因此 梁的整体稳定性满足要求 2 1 3 钢柱验算钢柱验算 强度验算 强度验算 惯性矩 2 37 108 4 1 74 107 4 回转半径 168 45 4 长细比 0 3600 168 21 43 5 22 0 3600 45 4 79 30 查表得 则 0 582 69 03 103 0 582 8 412 103 24 41 2 215 2 满足要求 局部稳定验算 局部稳定验算 翼缘 96 13 7 38 10 0 1 235 10 0 1 79 30 17 93 腹板 0 374 8 46 75 25 0 5 235 25 0 5 79 30 64 65 满足要求 刚度验算 刚度验算 79 30 150 满足要求 2 2钢框架杆件模型弹性分析钢框架杆件模型弹性分析 2 2 1 钢框架弹性阶段分析钢框架弹性阶段分析 首先计算自振频率 然后采用简谐波进行弹性阶段计算 输入波的幅值保持不变 改变频率与自振频率的比 值分别为 0 25 0 5 1 0 1 5 2 0 由于是水平抗震分析 故在用 Midas Gen 进行自振周期计算时仅需考虑 X 方向 结构分析 类型为 X Z 平面 将 1 0 恒载 0 5 活载以及自重转化为质量 且将自重转化为集中质量 体系只有 两个自由度 采用子空间迭代法 迭代振型为两阶 分析后表明第一阶的振型参与质量就达到了 100 故仅需考虑第一阶振型即该框架变成了单自由度体系 图 4 为该框架的振型图 弹塑性时程分析中静力荷载取 1 0DL 0 5LL 作用在结构上 在 midas Gen 分析中将此荷载转化 为质量 通过特征值分析 结构自振频率与周期见表 2 2 其中 1 2 阶模态如图 2 3 所示 表 2 2 结构自振频率与周期 6 22 阶数1 阶2 阶 自振频率 Hz 5 21593451 925757 自振周期 s 0 1917200 019258 a 第一阶模态 b 第二阶模态 图 2 3 结构模态图 结构阻尼矩阵采用瑞雷阻尼假定 1 2 阶阻尼比均取 2 图 2 4 0 25 倍基频谐振力下左柱顶点加速度响应时程 图 2 5 0 5 倍基频谐振力下左柱顶点加速度响应时程 7 22 图 2 6 基频谐振力下左柱顶点加速度响应时程 图 2 7 1 5 倍基频谐振力下左柱顶点加速度响应时程 图 2 8 2 倍基频谐振力下左柱顶点加速度响应时程 2 2 2 简谐波下的结果分析简谐波下的结果分析 共振下结构仍处于弹性状态 从而保证了简谐波作用下的分析为弹性分析 由于均采用 Newmark 时程分析方法 杆模型与纤维模型计算结果一致 从而可以看出弹性状态下 虽然纤维模 型相比杆模型截面划分更精细 但对计算结果精度提高很小 观察共振现象 由于结构各阶频率相距较远 取激励频率在结构一阶自振频率附近变化观察共 8 22 振反应 当激励频率小于结构一阶自振频率时 在阻尼作用下结构很快达到稳态反应 当激励频率 等于结构一阶自振频率时 结构出现共振反应 结构位移逐渐增大 由于阻尼的存在 最后达到稳 态反应 当激励频率大于结构一阶自振频率时 结构位移反应幅值开始阶段较大 并逐渐减小达到 稳态反应 在保持相同激励幅值时 结构共振下的位移反应相比非共振下大的多 2 3 钢框架弹塑性阶段分析钢框架弹塑性阶段分析 2 3 1 弹塑性分析介绍弹塑性分析介绍 结构弹塑性动力分析的基本动力方程为 2 3 1 g MUC UFMU 其中 为非线性恢复力向量 上式的增量方程为 F 2 3 2 jjjgj MUCUKUMU 式 2 3 1 可用动力分析的逐步积分法求解 全量形式的动力方程为 p FKU 在采用数值分析技术的前提下 结构线性地震反应分析与非线性地震反应分析的主要差别在于 刚度矩阵是否可变 对于弹塑性结构 在每一步增量反应计算之前 要先修正矩阵或中各 K p K 元素的量值 即所谓的刚度修正技术 刚度修正过程实质上是一个重新形成总刚度矩阵的过程 2 3 2 恢复力模型介绍恢复力模型介绍 恢复力模型是根据大量的从试验中获得的恢复力与变形关系曲线经适当抽象和简化而得到的实 用数学模型 是构件的抗震性能在结构弹塑性地震反应分析中具体体现 若仅用静力非线性分析 模型一般是指力与变形关系骨架曲线的数学模型 而如果是用于结构动力非线性时程分析 恢复力 模型不仅包含骨架曲线 同时也包含各阶段滞回环的数学模型 常见的恢复力模型有兰伯格 奥斯古德模型 克拉夫 Clough 模型和武田模型 兰伯格 奥 斯古德模型曾被广泛用于土体和包括钢结构在内的各种结构物的非线性反应分析 偶尔也用于钢筋 混凝土弯曲构件 克拉夫模型主要是针对钢筋混凝土受弯构件的恢复力特性提出的 而武田模型是 从较多的钢筋混凝土构架试验所得的恢复力特性曲线抽象得出的 适用于以弯曲破坏为主的情况 钢筋混凝土结构构件的恢复力模型一般分为曲线型和折线型两种 其中曲线型比较接近结构的 实际受力特性 结果比较精确 但是刚度计算比较复杂 因此 应用很少 折线型恢复力模型由若 干直线段所构成 刚度变化不连续 存在拐点问题 但刚度计算比较简单 故在实际工程中得到广 9 22 泛应用 本次分析中 将采用克拉夫 Clough 退化双线模型作为钢结构的恢复力模型 2 3 3 克拉夫 克拉夫 Clough 退化双线模型 退化双线模型 初次加载时沿着双折线骨架曲线移动 屈服后卸载路径沿着退化后的斜率移动 当反向加载时 指向反向最大变形点 反向没有发生屈服时 屈服点为最大变形点 克拉夫模型中认为全截面处于 开裂状态 截面的刚度由受拉钢筋的受弯屈服状态决定 对正向和负向可定义不同的屈服后的刚度 折减系数 适用于梁 柱 支撑构件 Kr Kr P1 P1 D1 D1 P D K0 K2 K2 maxmax DP maxmax DP P1 P1 D1 D1 P D 图 2 9 克拉夫 Clough 退化双线模型 克拉夫模型的骨架曲线由下列参数决定 正向和负向的第一屈服强度 Pl Pl 10 22 正向和负向的第一屈服强度 Dl Dl 正向和负向的第一屈服强度 Dl Dl 初始刚度 0 K 正向和负向的第二条折线刚度 2K 2K 0 21KK 0 21KK 正向和负向的第一屈服后刚度折减系数 1 1 正向和负向卸载时的刚度 r K r K 0 max 0 1 K D D KKr 0 max 0 1 K D D KKr 其中 为正向和负向的最大变形 没有屈服的区段使用屈服变形 为 max D max D 计算卸载刚度的幂阶 克拉夫模型的路径移动规则 1 时沿斜率为的直线移动 max DDl 0 K 2 变形第一次超过时或者超过当前的最大变形时 沿着斜率为 的第二折线D Dl 2K 2K 移动 3 在 状态下卸载时 沿着卸载刚度 的斜率移动 DlD DDl r K r K 4 卸载过程中荷载的符号发生变化时 将沿着指向反向最大变形点的直线移动 如反向未屈 服则指向反向屈服点 Clough 模型一般只适用于具有梭形滞回曲线的单纯受弯构件 2 3 4 纤维模型理论纤维模型理论 纤维模型是将构件单元分割为许多沿构件纵向的纤维的模型 纤维模型不仅可以准确模拟受弯 构件的力学特性 而且可以考虑截面内纤维的局部损伤状态 使用纤维模型时可利用纤维材料的应 力 应变关系和截面应变的分布形状假定较为准确地确定截面的弯矩 曲率关系 特别是可以考虑 轴力引起的中和轴变化 另外纤维模型同样可以考虑轴力和弯矩 两个弯矩之间的相互影响 但是 因为不能反映剪切破坏 所以一般用于剪切变形不大的线单元 本次分析采用的分析软件 MIDAS GEN 中的纤维模型使用了一下几个假定 1 截面的变形维持平截面与构件轴线垂直 11 22 2 不考虑钢筋与混凝土之间的滑移 本次分析采用钢结构 3 杆件单元截面形心的连线为直线 1 纤维模型的计算过程 首先假设在各单元的积分点上存在用纤维模型定义的截面 积分点数量最多不超过 20 个 积分 方法使用高斯 罗贝托方法 只有当积分点数量为 2 时才使用古典高斯方法 然后将前次时间步中计 算的两端的构件内力通过转换并排除刚体运动成分后获得构件两端的的广义单元内力 轴力和弯矩 最后通过内力的内插函数计算单元内各位置的内力 图 2 10 构件任意位置的构件内力和变形 单元荷载向量 12345 Q Q Q Q Q T Q 单元变形向量 12345 q q q q q T q 截面内力向量 zy x M x M x N x T D 截面变形向量 zy x x x x T d x x ii DbQ 其中 为内力内插函数 Force Interpolation Function 公式如下 x b xx 1000 LL xx x 00 10 LL 00001 b 其中 为单元长度 L 12 22 通过截面的柔度可计算截面的变形 通过各截面位置的轴向 弯曲变形可计算各纤维的轴向应 变 x ECS x axis ECS y axis ECS z axis zi yi i th fiber ECS y axis ECS z axis 图 2 11 纤维模型的截面分割示意图 y iiiz x x z y1 x x 其中 截面的位置 y x 在 x 处截面的绕单元坐标系 y 轴的曲率 z x 在 x 处截面的绕单元坐标系 z 轴的曲率 x x 在 x 处截面的轴向应变 yi 截面上第 i 个纤维的 y 轴位置 zi 截面上第 i 个纤维的 z 轴位置 i 第 i 个纤维的应变 各纤维的轴向应变 i 对应的纤维的应力和纤维的切线刚度可通过纤维材料的本构关系 constitutive relation 计算 并由此判断纤维的状态 将一个截面内所有纤维的应力进行积分可获得 截面的轴力和弯矩 对各纤维的切线刚度进行积分可获得截面的柔度矩阵 对单元内所有积分点上 的截面的柔度进行积分可得单元的柔度 n x n x n x j2jj iiiiiiiiii i 1i 1i 1 n x n x n x jj2j iiiiiiiiii i 1i 1i 1 n x n x n x jjj iiiiiiii i 1i 1i 1 EAy EAyz EAy x EAyzEAzEAz EAyEAzEA j k 1 x x jj fk j 步骤时位置 x 处截面的切线刚度 x j k j 步骤时位置 x 处截面的柔度 x j f 13 22 位置 x 处截面内总纤维数量 n x j 步骤时位置 x 处截面内第 i 个纤维的切线刚度 x j i E 第 i 个纤维的截面面积 i A 第 i 个纤维在截面内的位置 ii y z 计算截面不平衡力所需的截面内力按下式计算 n x n x n x jjj ii iii iii i 1i 1i 1 x A y A z A T j R D 在计算过程中使用牛顿 拉普森迭代方法计算至满足收敛条件 纤维模型中单元的非线性特性表 现在纤维的非线性应力 应变关系 材料本构关系 上 2 采用的钢材纤维本构模型 MIDAS 软件中的钢材纤维本构模型 修正梅内戈托与平托模型 Modified Menegotto Pinto Steel Model 本构模型形状为逐渐逼近按随动硬化 Kinematic Hardening 准则定义的双折线的曲线 即本构模 型在卸载路径和应变硬化区段之间的两个渐近线之间的转换区段为曲线 两个渐近线的交点和卸载 方向的最大变形点之间的距离越远 转换区段的形状越光滑 可利用这样的特性模拟包辛格效应 1 1 1 R R b b 其中 1 0 002 rr rr a RR a 钢纤维的应变 钢纤维的应力 r r 卸载点位置 在弹性状态时假定为 0 0 0 0 定义当前加载或卸载路径的两个渐近线的交点 b 刚度折减率 R0 a1 a2 常数 决定曲线状态的参数 程序使用了试验获得的优化值 在加载 卸载的方向上最大应变与 0的差 绝对值 最大应变的初始值假定与 Fy E 相同 参见 图 2 12 14 22 b E E 1 r r 1 0 0 1 2 r r 2 0 0 2 Fy 图 2 12 钢材纤维本构模型 双折线钢材纤维模型 常用的双折线本构模型 屈服前加载和卸载时使用弹性刚度 屈服后加载时使用屈服后刚度 屈服后卸载 再加载时使用弹性刚度 三折线钢材纤维模型 常用的三折线本构模型 可定义第一屈服和第二屈服刚度折减率 受拉区和受压区的刚度折减 率可以不同 屈服前的加载和卸载使用弹性刚度 屈服后加载时使用折减的刚度 屈服后的卸载和 再加载使用弹性刚度 折线形状也可以通过输入应力 应变定义 非对称双折线钢材纤维模型 该本构模型可以模拟钢筋的非线性特性 在受拉区可考虑屈服和拉断 在受压区可考虑屈服和 失稳后的破坏 2 3 5 克拉夫模型弹塑性地震波反应计算克拉夫模型弹塑性地震波反应计算 不同幅值的 El centro N S 波下的柱底弯矩 曲率关系图和柱顶位移时程图如下所示 图 2 13 幅值为 1 0g 的 El centro N S 波下的柱底弯矩 曲率关系图 15 22 图 2 14 幅值为 1 0g 的 El centro N S 波下的柱顶位移时程图 图 2 15 幅值为 1 6g 的 El centro N S 波下的柱底弯矩 曲率关系图 图 2 16 幅值为 1 6g 的 El centro N S 波下的柱顶位移时程图 16 22 图 2 17 幅值为 3g 的 El centro N S 波下的柱底弯矩 曲率关系图 图 2 18 幅值为 3g 的 El centro N S 波下的柱顶位移时程图 图 2 19 幅值为 5g 的 El centro N S 波下的柱底弯矩 曲率关系图 17 22 图 2 20 幅值为 3g 的 El centro N S 波下的柱顶位移时程图 2 3 5 纤维模型弹塑性地震反应计算纤维模型弹塑性地震反应计算 本次分析采用双折线钢材纤维模型 以 El centro N S 波为基本波 改变幅值 从弹性阶段计算 至弹塑性阶段 直至破坏 以节点 2 的 X 方向位移为时间的函数 画出的曲线如图 2 21 2 24 所示 图 2 20 纤维模型中梁柱纤维划分 18 22 图 2 21 增幅系数为 0 2 时节点 2 的 X 方向位移时程曲线 纤维模型 图 2 22 增幅系数为 1 0 时节点 2 的 X 方向位移时程曲线 纤维模型 图 2 1 增幅系数为 2 0 时节点 2 的 X 方向位移时程曲线 纤维模型 图 2 2 增幅系数为 10 时节点 2 的 X 方向位移时程曲线 纤维模型 右柱 单元 2 的弯矩 曲率曲线如图 2 25 2 28 所示 19 22 图 2 3 增幅系数为 0 2 时的右柱的弯矩 转角曲线 图 4 增幅系数为 1 0 时右柱的弯矩 转角曲线 20 22 图 2 5 增幅系数为 2 0 时的右柱的弯矩 转角曲线 当增幅系数为 0 2 1 0 2 0 时 结构未出现塑性铰 增幅系数为 10 时 结构进入弹塑性阶段并破坏 图 2 6 增幅系数为 10 时的右柱的弯矩 转角曲线 21 22 图 2 7 增幅系数为 10 时柱顶 节点 2 的弯矩 转角曲线 纤维模型 2 3 6 弹塑性地震分析评价弹塑性地震分析评价 1 从上述图中可以看出 无论是纤维截面模型还是杆件模型 随着地震激励的增加 结构的反 应也随之增大 当地震动幅值加大时 结构位移值也增大 且反应曲线变的稀疏 反应峰值减少 峰值变的平坦 这主要是幅值大 时间间隔短 构件来不及恢复变形 又遭受一个强大脉冲 所以 反应波谷变高