85边界非线性分析.docVIP

8-5 边界非线性时程分析(Boundary Nonlinear Time History Analysis) 8-5-1 分析模型 边界非线性单元由两个节点构成，在土木结构中用于模拟隔震和消能装置。当结构中设置了隔震或消能装置后，结构中其它构件一般设计成弹性构件。 分析模型一般由弹性构件(弹性系)和非线性弹簧构件(非弹性系)构成。midas程序支持模型中有非线性构件，但是需要注意的是当动力分析中同时考虑材料的非线性(非线性构件)和边界的非线性(非线性弹簧)时，时程分析方法一定要选择直接积分法。时程分析方法中的振型分解法不能同时考虑材料非线性(非线性构件)和边界非线性(非线性弹簧)。 midas程序中的每个非线性连接单元由六个方向的弹簧(轴向、两个剪切方向、扭转、两个受弯方向)构成，用户可以选择方向输入该方向的线性和非线性特性，没有输入特性的方向认为该方向上没有刚度，只输入线性特性没有输入非线性特性时认为该方向只具有弹性刚度。每个非线性连接单元至少要定义一个方向的非线性特性。 弹簧的特性由线性特性和非线性特性两部分组成。只具有线性特性的弹簧为线弹性弹簧，具有非线性特性的弹簧一般称为非线性弹簧(在程序中称为一般连接)。非线性弹簧中需要输入线性特性和非线性特性，其线性特性用于静力分析和线性动力分析(包括反应谱分析)，非线性特性用于非线性动力分析。 程序中提供“单元”和“内力”两种恢复力类型的非线性弹簧。“单元”类型顾名思义就是将弹簧视为“单元”(非线性构件)通过定义“铰”的内力-变形(刚度)关系即材料非线性特性来反映弹簧的非线性特性。“单元”类型的非线性弹簧有弹簧、阻尼、弹簧+阻尼三种类型。“内力”类型的恢复力模型是通过定义内力与位移(内力)、速度、阻尼的关系来反映弹簧的非线性特性。程序中提供的“内力”类型的非线性弹簧有粘弹性阻尼器、间隙、钩、滞回系统、铅芯橡胶支座、摩擦摆隔震装置等六种非线性弹簧。 本章中将介绍“内力”类型的非线性特性的定义方法，“单元”类型的非线性特性请参考“材料非线性”章节的介绍。 8-5-2 边界非线性时程分析概要 如前所述，边界非线性模型由非线性系统和线性系统构成，边界非线性分析通过将非线性系统发生的内力转换成线性系统的外部动力荷载在进行结构分析的。 图2.8.8 边界非线性的分析方法 程序中提供的边界非线性方法有非线性振型叠加法时程分析和非线性直接积分法。隔震或消能装置的单元特性用常微分方程表现并使用龙格-库塔法(Runge-Kutta Method)数值分析方法计算，计算的结果作为有效荷载用于求解动力方程。 非线性振型叠加法时程分析方法仅用于边界非线性+结构线性的结构模型的动力分析，且必须事先做特征值分析获得结构的自振周期和振型。非线性振型叠加法时程分析方法与直接积分法相比具有计算时间短的优点。一般来说提高特征值分析中的振型数量可提高与直接积分法的结果误差。 当分析模型中需要考虑构件的材料非线性特性时，需要使用非线性直接积分法。 8-5-3 使用振型叠加法做边界非线性时程分析 振型叠加法时程分析概要 包含内力类型的非线性弹簧的结构动力方程如下。 在此， : 质量矩阵 : 阻尼矩阵 : 不包含内力类型非线性弹簧的其它构件的弹性刚度矩阵 : 内力类型非线性弹簧的有效刚度矩阵 , : 转换矩阵 , , : 节点的位移、速度、加速度 : 动力荷载 : 内力类型非线性弹簧的非线性成分由有效刚度计算的内力 : 内力类型非线性弹簧的非线性成分的实际内力 方程右侧的与左侧由计算的节点力中内力类型非线性弹簧的非线性成分相抵， 因此只有对动力分析结果有影响。使用有效刚度矩阵的原因是因为非线性弹簧所处的位置当仅考虑结构构件的刚度矩阵时结构可能成为不稳定结构。有效刚度矩阵是由用户直接输入非线性弹簧刚度的线性特性值。 质量矩阵可由一致质量法或集中质量法计算，结构的自振周期和振型可由结构的质量矩阵和刚度矩阵使用特征值向量法(子空间迭代法或兰佐斯法)或多重里兹向量法计算。结构的阻尼矩阵使用振型阻尼方法计算。振型阻尼可使用由用户直接输入各振型阻尼的方法，也可使用质量和刚度因子计算各振型阻尼的瑞利阻尼法。程序还提供由各构件阻尼比通过应变能方法计算各振型阻尼的“应变能因子”法。 考虑振型的正交性可将上面的动力方程转换为各振型的动力方程。 在此， : 第i振型的振型形状向量 : 第i振型的阻尼比 : 第i振型的固有圆频率 ,, : 第i振型的广义位移、广义速度、广义加速度 右侧的和由内力类型的非线性弹簧在单元坐标系下的实际变形和应变决定。因为单元的实际变形包含了所有振型的结果，所以上面基于振型的运动方程式中各振型间依然相互关联。为了利用振型叠加法的优点，可通过对各时间步骤的和进行一些假设获得