有限元动力学问题有限单元法.pptVIP

第9章力有限元法洪

力学第1引言第2量矩和阻尼矩第3直接分法第4振型叠加法第5解的定性第6大型特征的解法第7减系自由度的方法第8小2

第1有限元力学方程的建立u力学中最常遇到的是构力学，它有两研究象。一是在运状下工作的机械或构，例如，高速旋的机，往复运的内燃机，以及高速运行的行器，如何保它运行的平性及构的安全性是极重要的研究。另一是承受力荷作用的工程构，例如建于地面的高建筑和厂房，核站的安全壳和交器，些构的破裂、覆和坍塌等破坏事故的生，将人民的生命造成巨大失。正确分析和构，在理和上都是具有重要意的。u力学研究的另一重要域是波在介中的播。3

三性力学的基本方程是：平衡方程几何方程（1.1）（在V域内）（在V域内）（1.2）（在V域内）（1.3）物理方程界条件（在S域内）（1.4）（在S域内）（1.5）uσ初始条件（1.6）4

力学方程建立：在荷作用下，于任一瞬，元点生虚位移，元内也生相的虚位移和虚。元内生的虚能:元除受荷外，有加速度和速度引起的性力和阻尼力，其中ρ材料密度，v是性阻尼系数。外力所做的虚功：式中，{P}、{P}、{P}分作用于元上的体力、面力vsc和集中力；V元面；A元面。

由于且形函数坐x、y、z的函数，与无关，因此有根据虚位移原理，有代入整理，可得元运方程

式中分称元的度矩、量矩和阻尼矩，它就是决定元性能的特性矩。称元点荷列，它是作用在元上的体力、面力和集中力向元点移置的果。在分析和静力分析中，元的度矩是相同的，外部荷的移置原理也一。

力学有限元分析基本步如下：（1）区域的离散化（2）构造插函数由于只空域行离散，所以元内位移u,v,w的插分表示:（1.7）其中8

（3）形成系的求解方程（1.8）其中分是系的点加速度向量和点速度向量，M,C,K和Q(t)分是系的量、阻尼、度和点荷向量。9

如果忽略阻尼的影响，运方程化如果上式的右端零，上式一步化（1.9）（1.10）是系的自有振方程，又称力特性方程。（4）求解运方程（5）算构的和力10

从以上步可以看出，和静力分析相比，在力分析中，由于性力和阻尼力出在平衡方程中，因此引入了量矩和阻尼矩，最后得到求解方程不是代数方程，而是常微分方程。其它的算步和静力分析是完全相同的。关于二常微分方程的解法有两：直接分法和振型叠加法。直接分法是直接运方程分。而振型叠加法是首先求解一无阻尼的自由振方程，然后用解得的特征向量，即固有振型运方程式行。力分析的算工作量很大，因此提高效率，省算工作量的数方案和方法是力分析研究工作中的重要成部分。目前两种普遍用的减自由度的方法是减法和力子构法。11

第2量矩和阻尼矩一、量矩和集中量矩元量矩称量矩。集中量矩假定元的量集中在点上，得到的量矩是角矩。以下分体元和构元行。1.体元介两种常用方法（1）第一种方法（2.1）每一行的主元其中，n是元的点数。式的力学意是：Mele素等于Me12中行所有元素之和，而非主元素零。

第2量矩和阻尼矩（1）第二种方法（2.2）式的力学意是：Mle每一行的主元素等于Me以放因子a，而非主元素零。中行主元素乘13

第2量矩和阻尼矩例1算平面力（）元的量Me矩和集中量矩Mle。元采用3点三角形元。（1）量矩位移插函数是（2.3）其中I是2×2位矩。（2.4）M,Ce,Ke和Q分是元的量、阻尼、度和荷矩。ee14

第2量矩和阻尼矩算得元的量矩（2.5）其中，W＝ρtA是元的量，t是元的厚度。15

第2量矩和阻尼矩（2）集中量矩按第一种方法算，得到集中量矩（2.6）此式的力学意是：在元的每个点上集中1/3的量。16

第2量矩和阻尼矩按第二种方法算，得到集中量矩与第一种方法果一。注：于8点矩形元，两种方法得到的集中量矩不同。在分析中，更多的是推荐用第二种方法来算集中量矩。2.构元2点典梁元、量矩和集中量矩如下所示：（1）量矩位移插函数是（2.7）其中17

第2量矩和阻尼矩算得元的量矩（2.8）（2.9）其中，l是元度，W＝ρlA是元的量，A是截面面。（2）集中量矩此式的力学意是在每个点上集中1/2的元量。18

第2量矩和阻尼矩需要指出，然量矩M在理上是正定的，但通常需要在算中行精确分才能保此性。如果算中采用低的分，M可能是奇异的，将使后的力分析生困，因此在Me的分次予注意。19

第2量矩和阻尼矩二、振型阻尼矩它是假定阻尼力正比于点运速度的果，通常均将介阻尼化种情况。元矩比例于元量矩。在以后的中，将知道系的固有振型于M和K是具有正交性的，因此固有振型于M和K的阻尼矩C也是具有正交性的。所以种阻尼矩称比例阻尼或振型阻尼。在分析中要精确地决定阻尼矩是相当困的，通常允将构的阻尼矩化M和K的性合。种振型阻尼称Rayleigh阻尼。20

求解方法求解运方程直接分法模叠加法式分式分完整矩法减矩法完整矩法减