有限单元法几个专题.pptVIP

6.1 动力学问题的有限元法 6.2 非线性问题的有限元法 6.1 动力学问题的有限元法 动力学问题在国民经济和科学技术的发展中有着广泛的应用领域。最经常遇到的是结构动力学问题有两类： 一类是在运动状态下工作的机械或结构： 另一类是承受动力载荷作用的工程结构： 这些结构的破裂、倾覆和垮塌将给人民的生命财产造成巨大的损失。 如何保证它们运行的平稳性及结构的安全性，正确分析和设计这类结构，在理论和实际上也都是具有意义的课题。 高速电机、汽轮机、离心压缩机，内燃机、冲压机床，高速车辆、飞行器等，它们承受着本身惯性及与周围介质或结构相互作用的动力载荷。 承受强风、水流、地震以及波浪等各种动力载荷的作用的高层建筑、桥梁、核电站的安全壳海洋石油平台等。 动力学研究的另一重要领域是波在介质中的传播问题。 6.1.1 弹性动力学问题的基本方程： 平衡方程(or运动方程)： 几何方程： 物理方程： 边界条件： 初始条件： 力边界： 位移边界： 式中ρ是质量密度，μ是阻尼系数。 平衡方程中出现惯性力和阻尼力这是弹性动力学和静力学相区别的基本特点之一。 位移、应变、应力也是时间的函数。也正因为如此，动力学问题的定解条件中还应包括初始条件。 6.1 动力学问题的有限元法 6.1.2 有限元法求解动力问题的基本步骤： 6.1 动力学问题的有限元法 在动力分析中引入了时间坐标。在有限元分析中一般采用部分离散的方法，即只对空间域进行离散，这样一来，此步骤和静力分析相同。 2、构造插值函数 1、连续区域的离散化 此时，结点位移是时间的函数。 6.1.2 有限元法求解动力问题的基本步骤： 6.1 动力学问题的有限元法 3、形成系统的求解方程 平衡方程及力的边界条件的等效积分形式为： 分部积分，并代入物理方程， 再将离散后的位移代入上式，注意到结点位移变分的任意性 (请同学们课后自己证明) 最终得到系统的求解方程： 如果忽略阻尼的影响，可简化为: 系统的求解方程为： [M]是系统的质量矩阵、 [C ]是系统的阻尼矩阵、 [K]是系统的刚度矩阵、 {Q(t)}是系统的结点载荷向量， 分别由各自的单元矩阵和向量集成。 6.1.2 有限元法求解动力问题的基本步骤： 6.1 动力学问题的有限元法 4、求解运动方程 运动方程的求解方法是本章着重讨论的内容。 5、计算结构的应变和应力 6.1.2 有限元法求解动力问题的基本步骤： 6.1 动力学问题的有限元法 从以上各步骤可以看出： B、最后得到的求解方程不是代数方程组，而是常微分方程组。其他计算步骤和静力分析是完全相同的。 A、和静力分析相比较，在动力分析中由于惯性力和阻尼力出现在平衡方程中，因此引入了质量矩阵和阻尼矩阵。 6.1.3 单元质量矩阵 6.1 动力学问题的有限元法 一致(协调)质量矩阵： 质量矩阵的导出和导出刚度矩阵所根据的原理及所采用位移插值函数是一致的，同时质量分布也是按照实际分布情况考虑的。 在有限元法中还经常采用所谓集中〔或团聚)质量矩阵。它假定单元的质量集中在结点上，这样得到的质量矩阵是对角线矩阵。 集中质量矩阵： 以3结点三角形单元为例： 其中W =ρtA是单元的质量, t是单元的厚度。 6.1.4 单元阻尼矩阵 6.1 动力学问题的有限元法 1、基于和协调质量矩阵的同样理由称为协调阻尼矩阵。它是假定阻尼力正比于质点运动速度的结果，通常均将介质阻尼简化为这种情况。这时单元阻尼矩阵比例于单元质量矩阵。 2、除此而外，还有比例于应变速度的阻尼，例如由于材料内摩擦引起的结构阻尼通常可简化为这种情况，这时单元阻尼矩阵可表示成： 此单元阻尼矩阵比例于单元刚度矩阵。 6.1.5 直接积分法： 6.1 动力学问题的有限元法 常微分方程组的解法，原则上可利用求解常微分方程组的常用方法(例如Runge-K utter方法)求解。 但是在有限元动力分析中，因为矩阵阶数很高，用这些常用算法一般是不经济的，所以只对少数有效的方法有兴趣。这些方法可以分为两类:直接积分法和振型叠加法。 直接积分是指在积分运动方程之前不进行方程形式的变换，而直接进行逐步数值积分。 通常的直接积分法是基于两个概念： 一、将在求解域0tT内的任何时刻t 都应满足运动方程的要求，代之以仅在一定条件下近似地满足运动方程，例如可以仅在相隔Ot的离散的时间点满足运动方程。 二、在一定数目的Δt区域内，假设位移、速度、加速度的函数形式。 在以下的讨论中，假定时间t =0的位移、速度、加速度已知。并假定时间求解域0-T被等分为n个时间间隔Δt=T/n，在讨论具体算法时，假定0, Δt, 2Δt, … t时刻的解已经求得，目的在