N-第9章 机械结构的动力有限元(li 10.24)

1、第九章机械结构的动力有限元 对于机械结构的动力学分析问题，有限元法也是非常有效的对于机械结构的动力学分析问题，有限元法也是非常有效的数值计算工具。数值计算工具。其主要过程为：首先进行每个单元的特性分析，其主要过程为：首先进行每个单元的特性分析，包括进行单元刚度矩阵、单元质量矩阵的计算、单元阻尼矩阵包括进行单元刚度矩阵、单元质量矩阵的计算、单元阻尼矩阵的计算，再把各个单元的特性矩阵组集起来，组成结构的总刚的计算，再把各个单元的特性矩阵组集起来，组成结构的总刚度矩阵、总质量矩阵和总阻尼矩阵，从而形成整个结构的动力度矩阵、总质量矩阵和总阻尼矩阵，从而形成整个结构的动力学方程学方程，之后再进行求解。，

2、之后再进行求解。本章概述本章概述 本章主要介绍机械结构动力学的有限元方法，然后介绍利用有限元求解机械结构的固有动特性、分析典型工况下的动态响应，并给出了应用举例。MC 机械结构的静力学及动力学有限元方程 (9.1) (9.2)K FMCKF整个结构节点位移向量，常值（静力学），随时间变化（动力学）总刚度矩阵总阻尼矩阵总质量矩阵节点载荷向量KF引入约束条件，求固有特性，求振动响应均为大型带状稀疏矩阵(9.18) 在有限元动力学分析中，单元质量矩阵、单元阻尼矩单元质量矩阵、单元阻尼矩阵阵与单元刚度矩阵都是十分重要的。为了更好地理解单元质量矩阵和单元阻尼矩阵，以下做进一步的解释。 如下式所表达的单元

3、质量矩阵又称为单元协调质量矩单元协调质量矩阵或单元一致质量矩阵阵或单元一致质量矩阵： 单元协调质量矩阵采用了与推导单元刚度矩阵一致的形函数矩阵。TdeeVVmNN (9.19) 计算得到平面三角形常应变单元的协调质量矩阵的具体表达式为20101002010110201001020112101020010102etm (9.20) 对于等参数单元，设其形函数矩阵为，单元协调质量矩阵的表达式为 在有限元动力学计算中，还经常采用单元集中质量矩阵的形式。即假设单元的质量分散在单元的所有节点上，即把每个单元的分布质量按静力学的平行分解原理分配到每个节点上，形成一个阶数等于单元自由度数的主对角线质量矩阵，

4、其它非主角线元素均为0。例如，对于平面三角形常应变单元，其单元集中质量矩阵为111111|d d deT mNNJ (9.21)1000000100000010000001003000010000001etm 对于前面所述的单元阻尼矩阵即单元协调阻尼矩阵，TdeeVvVcNN (9.22) 它取决于假设阻尼力正比于质点运动速度，即通常的粘性阻尼假设。除了这种形式的单元阻尼矩阵之外，在有限元动力学分析中，还经常采用另外一种形式的比例阻尼，即采用单元质量矩阵和单元刚度矩阵的线性组合作为单元阻尼矩阵，如下式 其中，是不依赖于频率的常数，称为比例阻尼系数。比例阻尼又称为瑞利阻尼（Rayleigh da

5、mping）。(9.23)eeecmk (9.30)不考虑阻尼项及激振力项，机械结构的动力学方程为( )( )0ttM+K=它的解可以假设为以下形式0sin()tt (9.31) 其中， 是 n阶向量， 是振动圆频率，t是时间变量，t0是由初始条件确定的时间常数。 将式(9.30)代入式(9.31)，得到如下特征方程（即广义特征值问题)20KM或20KM (9.32)v,d=eig(k,m); % eig求部分特征值v d均为方阵 d对角线即是固有频率 通常把结构动力分析分为结构的瞬态响应分析和结构的谐振响应分析。动力响应问题就是求解动力学方程(9.29)，即在激振力的作用下，求出作为时间或频率的函数的位移。 模态叠加法