（本科）第章 混合型有限元公式课件.pptxVIP

第12章 混合型有限元公式课程主讲人：第12章 混合型有限元公式 12.1 引 言 1 概念描述问题的场方程1）为了描述场问题，如果 u 中任何分量都不能消去，则称公式为 不可约型（irreducible）。2）如果 u 中的某些变量可以消去，则公式为混合型（mixed）。可以构造不同型式的有限元方程。2 二维热传导问题本构方程连续方程边界条件1）将本构方程代入连续方程得不可约型， 为独立变量。2）由连续方程（罚参数）代入本构方程不可约型，q 为独立变量。3）原方程混合型， T和 q 为变量。12.2 热传导问题的混合型公式 1 混合型方程的离散化（1）（2）本构方程等效积分形式连续方程和热流边界条件的等效积分形式单元插值：任意性任意性可导出如下方程讨论：1）混合型式对插值函数连续性要求不同。由于系数矩阵计算中未出 现Nq的导数，插值函数在单元之间可以不连续，而要求NT 是C0连 续的。采用不可约型时，要求C1连续。2）若重点研究q，则混合型给出较高的精度。3）混合型方程的系数矩阵中通常有零对角元素，给标准高斯消去法 解方程带来困难。 称为主变量； 称为约束变量。4）变量增加，方程规模增大。这一困难可以用迭代法求解方程克服。2 混合型方程解的稳定性如果A非奇异此方程有唯一解的条件非奇异可以证明，为使 H 非奇异，要求 的未知量个数 的未知量个数12.3 弹性力学问题的混合型公式 1 以位移为基本变量的有限元公式虚功原理本构关系几何关系位移插值以位移为基本变量的不可约型有限元公式2 位移－应力混合型公式平衡方程和力边界条件的等效积分弱形式B=LNu以位移和应力为基本变量的混合型有限元公式方程要有唯一解，必须满足条件等效Hellinger－Reissner变分原理证明：标量3 位移－应力－应变混合型公式本构方程和几何方程的等效积分以及虚功原理等价于胡－鹫广义变分原理设可得方程其中B=LNu4 解的稳定性要求1）位移－应力混合型2）位移－应力－应变混合型12.4 不可压缩弹性体问题 1 材料的不可压缩性工程中常见材料：橡胶、塑料等为不可压缩弹性体弹性本构关系有限元方程材料不可压：无法求解。1）位移公式无效。位移公式的解出现振荡。2）解决途径： 用罚函数法或Lagrange法将不可压缩条件（体积变形为零）作为附加条件引入泛函。采用以位移和压力为基本变量的混合型公式，能够处理完全不可压缩和接近不可压缩的问题。2 应力和应变偏量平均应力（压力）体积应变压力和体积应变的关系K：体积模量偏应变偏应力偏应力与偏应变之间的关系3 位移-压力混合型公式代入虚功原理必要条件4 位移-压力-体积应变混合型公式12.5 混合型方程简单迭代解法 1 混合型方程的一般形式求解的困难：1）对角元有零元素；2）自由度增加措施：采用迭代法求解可以大大降低求解代价。2 迭代求解的一般过程第 k 次迭代第二个方程的残值收敛加速矩阵12.6 多个弹性子域之间的连接 背景：并行算法的出现区域分解算法求解域不同部分的连接虚功表达式的弱形式交界面上位移连续条件弱形式离散化可以扩展到多个子域的情况也可以建立受拉格朗日乘子约束的变分原理，通过如下泛函的驻值条件得到方程12.7 小结热传导问题的混合型公式混合型方程的离散化解稳定性必要条件弹性力学问题的混合型公式位移-应力混合型公式解稳定性必要条件位移-应力-应变混合型公式解稳定性必要条件不可压缩弹性体问题位移-压力混合公式位移-压力-体积应变混合公式混合型方程简单迭代解法多个弹性子域之间的连接两个子域