有限元数值方法-讲稿9：动力学问题有限元法.pptVIP

State Key Laboratory of Structural Analysis for Industrial Equipment Dalian University of Technology, China 有限元与数值方法第九讲 第五章 弹性动力学问题的有限元法 结构动力有限元分析（三维实体）步骤 连续区域的离散化：网格，单元表示 构造插值函数： 形成系统的求解方程 结构动力有限元分析步骤 4. 求解动力方程 这是一个常系数的微分方程组 基本解法：Runge-Kutta：太费时 常用解法：2种： （1） 直接积分法：直接对运动方程进行积分 （2） 振型叠加法：首先求解无阻尼自由振动问题（特征值问题），将位移按振型展开，表示成振型的线性组合，然后代入控制方程，获得关于待定系数的方程。 优点：速度快！ 5. 计算应变；应力等 2. Raileigh商方法 3. 子空间迭代法 5.4 动力响应问题 单自由度系统的振动方程的求解 2. 基本的直接积分格式：线性加速度法 2. 基本的直接积分格式：线性加速度法 Wilson-?法 逐步积分法的精度与稳定性 则原振动方程就可转化如下： 其中： 原问题近似转化为求解q 阶方程组（降阶） 减缩自由度的方法 本节参考： The finite element method. Volume 1: the basis. O.C. Zienkiewicz. Chapter 19. 有限元分析的概念与应用。库克著。12.7节 特征值分析和子空间映射是计算力学最重要的概念 假设在很小的时间段内（t, )加速度是线性变化的 对上式积分两次： 该方法是线性加速度法的扩展（将线性加速度假设延拓至比一个步长更大的时间区间）。其假设： 当 ?＝1 时，为线性加速度法； 当 ?1.37 时，是无条件稳定格式。 该方法是隐式方法。 对于 ，可证明以上假设相当于时间间隔内的线性加速度假设： Newmark方法 假设时刻 的速度和位移演变格式为： 若假设 ，则有 恰好是Newmark方法中假设的速度和位移演变格式 但 时，Newmark方法不是无条件稳定的 Newmark方法 若假设 ，则等价于等平均加速度法（无条件稳定格式）： 恰好是Newmark方法中假设的速度和位移演变格式 Newmark方法是隐式方法（需要求解线性方程组），但是方程组的系数矩阵只需分解一次 将以上位移、速度和加速度表达式代入动力学方程 根据假设的位移、速度与加速度的关系可反推出下一时刻加速度表达式： 得到以下一时刻位移为自变量的线性方程组： 逐步积分法的精度 逐步积分法的精度由空间离散误差和时间离散误差确定。为保障时间离散精度，时间步长的选择应该与关心响应的最小周期相适合，例如一般在最小周期的1/10以下。 逐步积分法的稳定性 无条件稳定性：在任何时间步长下，对于任何初始条件，若解都不会发散，则称算法无条件稳定。 条件稳定性：若对于时间步长小于某个临界值的情形，任意初始条件下的解都不会发散，则称算法为条件稳定。 例如，对于Newmark算法，满足下列条件时，是无条件稳定的 振型叠加法 基本思路：特征值分析＋空间映射（降维） n-DOF方程： 仅取前 m 阶振型 物理自由度 模态坐标 原方程成为： 其中的系数矩阵分别为： 可用杜哈梅(Duhamel)积分求解 其中： 讨论： 振型叠加法将系统位移转换到以固有振型为基向量的子空间中。若m=n,则方程与原方程完全等价。根据问题求解的需要，适当选取m阶低阶振型参与计算（考虑载荷变化的频率、响应的空间复杂程度等确定）。 对于一定形式的阻尼矩阵（如比例阻尼），将得到非耦合的m个单自由度系统的运动方程。 还可采用振型加速算法提高振型收敛速度。 对于一般的阻尼（无法解耦）或非线性问题，振型叠加法应用受到限制。 频率响应分析（模态分解法） 是由外载荷或者基础振动引入的简谐激励（地震波、转子不平衡力等），不妨假设形式为 承受简谐激励的系统的稳态响应方程： u(t)按结构前 q 阶振型展开 其中， yi 代表各振型在响应中所占的分量，是时间的函数。 对于特定的阻尼模型，以上方程可以解耦，进一步减小了计算量。 其稳态解 其中： Guyan缩减法(静态凝聚) 假设从动自由度受主自由度由静力平衡支配： 即 将所有自由度划分