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

有限元与数值方法第九讲 第五章 弹性动力学问题的有限元法;5.1 概述：弹性动力学问题的有限元方法;弹性动力学问题的有限元提法;虚功原理：;考虑具有n个节点的单元：;应变矩阵;单元质量阵;动力学有限单元列式：;结构动力有限元分析（三维实体）步骤;阻尼矩阵;结构动力有限元分析步骤;5.2. 质量矩阵和阻尼矩阵;阻尼的概念;阻尼的概念;粘滞阻尼 粘滞阻尼假定阻尼力与速度成正比，无论对简谐振动还是非简谐振动得到的振动方程均是线性方程，不仅求解方便，而且能够方便地表达阻尼对频率、共振等的影响，是应用最为广泛的阻尼模型，通过将阻尼系数与结构体系的质量、刚度相联系，可以方便地构造出具体的阻尼系数，是目前最常用的阻尼表达方法。 ;滞回阻尼(复阻尼) 滞回阻尼假定应力应变间存在着相位差，从而振动一周有耗能发生。其特点是可以得到不随频率改变的振型阻尼比。 滞回阻尼将导致复数形式的刚度(所以这种阻尼又称复阻尼)，这对于一般时程分析而言，计算将比较复杂，耗时耗力，因而复阻尼实际应用并不多。;假设：;5.3 几种结构动力分析问题：自由振动分析;固有模态分析; 模态可以相差一个常数（模态向量因而不是唯一确定的）;振型矩阵(模态矩阵)：;逆迭代法（最低阶） 正迭代法（最高阶） （迁移）子空间迭代法 广义Jacobi法 Rayleigh商迭代算法 Lanczos算法 ……;1.逆迭代法;逆迭代法步骤：;2. Raileigh商方法;3. 子空间迭代法;5.4 动力响应问题;单自由度系统的振动方程的求解;直接积分法;;2. 基本的直接积分格式：线性加速度法;2. 基本的直接积分格式：线性加速度法;Wilson-?法;;但 时，Newmark方法不是无条件稳定的;Newmark方法是隐式方法（需要求解线性方程组），但是方程组的系数矩阵只需分解一次;逐步积分法的精度与稳定性;振型叠加法;原方程成为：;可用杜哈梅(Duhamel)积分求解;讨论： 振型叠加法将系统位移转换到以固有振型为基向量的子空间中。若m=n,则方程与原方程完全等价。根据问题求解的需要，适当选取m阶低阶振型参与计算（考虑载荷变化的频率、响应的空间复杂程度等确定）。 对于一定形式的阻尼矩阵（如比??阻尼），将得到非耦合的m个单自由度系统的运动方程。 还可采用振型加速算法提高振型收敛速度。 对于一般的阻尼（无法解耦）或非线性问题，振型叠加法应用受到限制。;频率响应分析（模态分解法）;则原振动方程就可转化如下： 其中： 原问题近似转化为求解q 阶方程组（降阶）;减缩自由度的方法;例如，对于无阻尼自由振动问题，缩减自由度后的方程成为：;，边界自由度 ，;(2)约束模态（依次释放界面上的各自由度并令其取单位位移所得到的静位移）：;通常略去高阶固定主模态，将 n 个物理自由度用 k+o 个模态表示成：;;耦合场问题－流固耦合;流体域控制微分方程及定解条件;其变分为零给出;将给出;流固耦合系统的离散形式;对于不可压缩流体