隐式与显式动力学的区别-弹性动力学有限元基本解法.pdfVIP

隐式与显式动力学的区别 1.1. 弹性动力学有限元基本解法 结构系统的通用运动学方程为：   M U  CU KU R （1） t 求解该动力学振动响应主要有三类方法：（1）时域法（2 ）频域法（3 ）响应谱法 时域法又可分为：（1）直接积分法，（2 ）模态叠加法。直接积分法又可分为中心差分法（显式），Wilson  （隐式）法以及Newmark （隐 式）法等。 本文介绍中心差分法（显式）与Newmark （隐式）法。 1 中心差分法（显式） 1 2 n n 假定0， ， ，…， 时刻的节点位移，速度与加速度均为已知，现求解 时刻的结构响应。中心差分法对加速度，速度的导数采用 t t t t (t t) 中心差分代替，即为：  1 U (U 2U U ) t t2 tt t tt  1 Ut (Utt Utt ) (2) 2t 将（2 ）式代入（1）式后整理得到 ˆ ˆ M U R （3 ） tt t 式（3 ）中 ˆ 1 1 M M  C t2 2t ˆ 2 1 1 R R  (K  M )U  ( M  C)U t t t2 t t2 2t tt 分别称为有效质量矩阵，有效载荷矢量。 ， ， ， 为结构载荷，质量，阻尼，刚度矩阵。 R M C K 1 求解线性方程组（3 ），即可获得t t 时刻的节点位移向量Utt ，将Utt 代回几何方程与物理方程，可得t t 时刻的单元应力和应变。 ˆ ˆ 中心差分法在求解 瞬时的位移U 时，只需 时刻以前的状态变量 和U ，然后计算出有效质量矩阵 ，有效载荷矢量 ，即可 t t t t U M R