结构有限元法的数学基础.docVIP

结构有限元法的数学基础 （航空学院 航空结构工程系 王生楠） 一、弹性力学的基本方程 在弹性力学问题里，通常是已知物体的边界、物体的弹性常数、物体所受的面力以及位移边界条件，求解其应力分量、应变分量和位移分量。为了由已知量求出未知量，必须建立这些已知量与未知量之间的关系，以及未知量之间的关系，从而导出求解的方程。 设在卡氏直角坐标系中，体积为的弹性体中任意一点的坐标参数为，该点的位移参数为、应力分量为以及应变分量为。由线弹性力学理论，可以得到如下的用于描述一个弹性静力学小位移变形问题的基本方程式。 （1）力的平衡方程 用张量形式描述： （在内） （1-1） 式中表示体力，表示应力分量对坐标分量的偏导数（以下相同）。 （2）应变位移关系式（几何关系） ， ， ， ， ， 用张量形式描述： （在内） （1-2） 注意：应变张量与工程切应变之间的关系为。 变形协调方程（从应变位移关系式中消去位移分量）： ，或 （3）应力应变关系式、本构方程（物理关系） （1-3） （1-4） 式中为弹性模量系数，为劲度系数，和都具有对称性。 （4）在弹性体的边界上，表面S可划分为两部分：外力已知的边界及位移为已知的边界，即，前者称为面力边界，后者称为位移边界。 在面力边界上， （1-5） 式中为已知边界力，为的边界外法线向量与坐标轴夹角的方向余弦。 在位移边界上， （1-6） 式中为已知边界位移。（1-5）式和（1-6）式统称为“边界条件”。 上述的诸方程共有15个，即3个平衡方程，6个应变位移关系方程，6个物理关系方程。而未知变量也共计15个：6个应力分量，6个应变分量和3个位移分量。因此基本方程是可以求解的。 这里没有考虑变形协调方程，原因是位移已经作为基本未知量。对于任意的单值连续的位移函数，如果设其有三阶的连续导数，则变形协调方程仅仅是几何方程微分的结果，自然地满足，所以位移作为基本未知量时，则不需要考虑变形协调方程。 要使基本方程有确定的解，还要有对应的面力或位移边界条件。 几何方程 几何方程 物理方程 平衡微分方程 位移分量 ui 应变分量 εij 应力分量 σij 体力 F 变形协调方程 位移边界条件 面力边界条件 混合边界条件 已知位移 ui 已知面力Ti 二、弹性力学问题求解的一般性原理 弹性力学问题就是：在给定的边界条件下，求解基本方程，求出未知量。 当然，具体求解弹性力学问题时，并不需要同时求解15个基本未知量，可以而且必须做出必要的简化。根据几何方程和本构方程可见，位移、应力和应变分量之间不是相互独立的。 假如已知位移分量，通过几何方程可以得到应变分量，然后通过物理方程可以得到应力分量。反之，如果已知应力分量，也可通过物理方程得到应变分量，再由几何方程的积分求出位移分量，不过这时的应变分量必须满足一组补充方程，即变形协调方程。 基于上述的理由，为简化求解的难度，选取部分未知量作为基本未知量。 位移解法。若以位移函数作为基本未知量求解，称为位移解法。以位移分量作为基本未知函数（未知量），从基本方程和边界条件中消去应力分量和形变分量（应变分量），导出只含位移分量的方程和相应的边界条件，并由此解出位移分量，然后再由几何方程求出形变分量，由物理方程求出应力分量。此类解法类似于结构力学中的位移法。 通过几何方程将位移函数表达为应变分量，再通过物理方程将其表达为应力分量，代入平衡微分方程即可得到位移解法的基本方程。 应力解法。若以应力函数作为基本未知量，称为应力解法。以应力分量作为基本未知函数（未知量），从基本方程和边界条件中消去位移分量和形变分量（应变分量），导出只含应力分量的方程和相应的边界条件，并由此解出应力分量，然后再由物理方程求出形变分量，由几何方程求出位移分量。此类解法类似于结构力学中的力法。 以应力分量作为基本未知函数求解弹性力学问题时，应力分量必须满足平衡微分方程和面力边界条件。但是仅此还不够，仅仅满足平衡微分方程和面力边界条件的应力分量并不是真正的应力，因为这组应力分量求出的应变分量代入几何方程，将可能得到一组矛盾方程，不可能求出单值连续的位移分量。要使这组方程不矛盾，则要求应力分量不仅满足平衡微分方程和面力边界条件，而且要求应力分量对应的应变分量必须满足变形协调方程。 这个问题也可以从物理上解释，应力分量满足平衡微分方程和面力边界条件，只能保证物体的平衡，但是不能保证物体的连续。只有这组应力分量求出的应变分量满足变形协调方程时，才能保证变形后的物体是