[工学]第01章 有限元方法基础.pptVIP

Contents Basic ideals Basic equations of finite element in Mechanics (Omitted intension) 1. Basic ideals 有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 有限元法分析计算的思路和做法可归纳如下：  1） 物体离散化 2） 单元特性分析 3） 单元组集 4） 求解未知节点位移 1） 物体离散化  将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。 　　离散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算进度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。 　　所以有限元中分析的结构已不是原有的物体或结构物，而是同型材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。  2） 单元特性分析 A、选择位移模式 　　在有限单元法中，选择节点位移作为基本知量时称为位移法；选择节点力作为基本知量时称为力法；取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算自动化，所以，在有限单元法中位移法应用范围最广。 　　 2） 单元特性分析  A、选择位移模式 　　当采用位移法时，物体或结构物离散化之后，就可把单元总的一些物理量如位移，应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。通常，有限元法将位移表示为坐标变量的简单函数。这种函数称为位移模式或位移函数，如y=a其中a是待定系数，y是与坐标有关的某种函数。 2） 单元特性分析 B、分析单元的力学性质 　　根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵，这是有限元法的基本步骤之一。 2） 单元特性分析 C、 计算等效节点力 　　物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边传递到另一个单元中去的。因而，这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去，也就是用等效的节点力来代替所有作用在单元上得力。 3） 单元组集  　　利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来，形成整体的有限 元方程，K是整体结构的刚度矩阵；q是节点位移列阵；f是载荷列阵。 4） 求解未知节点位移 　　解有限元方程式得出位移。这里，可以根据方程组的具体特点来选择合适的计算方法。 通过上述分析，可以看出，有限单元法的基本思想是一分一合，分是为了进行单元分析，合则为了对整体结构进行综合分析。 直杆受自重作用的简单拉伸问题 就整个直杆来说，位移函数U(x)是未知的，但对每一单元可以近似地假设一位移函数，它在结点上等于结点位移。此处，假设单元中的位移按线性分布 ，即： 有了位移插值函数，就可以按材料力学公式求出应变和应力用节点位移表示的公式： 外载荷与结点的平衡方程 为第i个结点上承受的外载荷 假定将直杆分割成3个单元，每个单元长为a=L/3，则对结点2，3，4列出的平衡方程为： 解线性代数方程组 联立求解方程组得： 有限单元法解题的一般步骤 结构的离散化 选择位移模式 建立平衡方程 求解节点位移 计算单元中的应力和应变 结构的离散化 将分析的结构物分割成有限个单元体，使相邻的单元体仅在节点处相连接，而以如此单元的结合体去代替原来的结构。 选择位移模式 首先对单元假设一个位移差值函数，或称之为位移模式，得到用节点位移表示单元体内任一点的唯一的关系式 建立平衡方程 可利用最小势能原理建立结构的节点载荷和节点位移之间的关系式，即结构的平衡方程 求解结点位移 将线性代数方程组 代入边界条件后，经解算可求得所有未知的结点位移。 计算单元中的应变和应力 依据求得的结点位移，由 可求得单元中任一点的应变和应力。 平面问题的有限单元法 结构的离散化 用有限元法对结构进行应力分析时，首先要将结构进行离散化。即将一个连续体看成由有限个单元组成的体系。