2-5有限元法在流体力学中的应用.docxVIP

第五章 有限元法在流体力学中的应用 本章介绍有限元法在求解理想流体在粘性流体运动中的应用。 讨论了绕圆柱体、 翼型和轴对称物体的势流， 分析了求解粘性流动的流函数—涡度法流函数法和速 度—压力法，同时导出粘性不可压流体的虚功原理。 § 1 不可压无粘流动 真实流体是有粘性和可压缩的， 理想不可压流体模型使数学问题简化， 又能较好 地反映许多流动现象。 圆柱绕流 本节详细讨论有限无法的解题步骤。 考虑两平板间的圆柱绕流． 如图 5—1 所示。 为了减小计算工作量，根据流动的对称性可取左上方的 l ／4 流动区域作为计算 区域。 选用流函数方法，则流函数 应满足以下 Laplace 方程和边界条件 203 2 2 x2 y 2 0 (x, y) 0 ( x, y) 2 ( x, y) y ( x, y) 0 ( x, y) n  aec流线 bd 流线 (5-1) ab流线 cd流线 将计算区域划分成 10 个三角形单元。单元序号、总体结点号和局部结点号都按 规律编排．如图 5— 2 所示。 从剖分图上所表示的总体结点号与单元结点号的关系，可以建立联缀表于下 元素序号 1 2 3 4 5 6 7 8 9 10 总体 n1 1 4 4 4 2 2 6 6 5 5 结点 n2 4 5 9 8 6 5 7 10 10 9 号 n3 2 2 5 9 3 6 3 7 8 10 表 5-1 各结点的坐标值可在图 5— 2 上读出。如果要输入计算机运算必须列表。本质边 204 界结点号与该点的流函数值列于下表 边界结点号 1 2 3 4 8 9 10 n 流函数 Φ 2 2 2 1 0 0 0 表 5-2 选用平面线性三角形元素，插值函数为 (3—15)式。对二维 Laplace 方程进行元素分析，得到了单元系数矩阵计算公式 (3—19)和输入向量计算公式 (3—20)。现在对全部元素逐个计算系数矩阵。 例如元素 1，其结点坐标为 x1 =0, y1 =2; x2 =0, y2 =1; x3 =2.5, y3 =2. 由 (3—15)式可得 a1 x3 x2 2.5; a2 x1 x3 2.5 a3 x2 x1 0 , b1 y2 y3 1 ; b2 y3 y1 0 ; b3 y1 y2 1 ; A0 1.25 从 (3—19)式可计算出 K1 1.45 1.25 0.2 K1 1.25 0 对称 0.2 依次可计算出全部子矩阵 0.2 0.2 0 K 2 1.45 1.25 1.25 0.2 0 0.2 K3 1.25 1.25 1.45 205 1.25 1.25 0 K 4 1.45 0.2 0.2 0.5 0.5 0 K 5 0 0.5 0.5 0.5 0 0.5 K 6 0.5 0.5 1 0.5 0.5 0 K 7 1 0.5 0.5 0.5 0 0.5 K8 0.5 0.5 1 0.5 0 0.5 K9 0.5 0.5 1 1 0.5 0.5 K10 0.5 0 0.5 根据联缀表把元素矩阵组合成总体系数矩阵 A= 1.45 0.2 0 1.25 0 2.45 0 0 1.25 1.0 1 0 0 0.5 0.5 2.9 0.4 0 0 1.25 4.9 1 0 0 1.75 0.5 4.0 1 0 0 0.5 对称 2 0 0 0.5 1.45 0.2 0 1.95 0 1.5 206 矩阵中零元素没有一一写出，下三角部分与上三角部分对称。 从(3—20)式计算元素输入向量，由于流函数满足齐次的自然边界条件 qn 0 ,所以输入向量为零，总体输入向量也为零，这样就得了总体有限元 n 方程 . N A B 式中 : N T , , , T B 0,0, ,0 用缩减方程的重新编号修正方法施加边界条件， 本质边界结点的函数值是已知 的。把它们代入方程，修正右端项，再减去相应的方程，整理得 4.9 1 0 1 4 1 0 1 2 解方程得到  5 6 7  2.9 3 1 5 ＝0.845， 6 =1.241， 7 =1.121 这样求出了全部结点上的流函数。 为了求出每个单元形心处的速度， 可以由 单元的流函数近似表达式求导计算。对元素 e 来 说，有 T I y 1 u y x  T I 1 ,a2 , a3 T I y a1 2 A 2A0 3 TI 1 b1 , b2 , b3 I T I x 2 A B 例如单元 1 =2, 2 =1, 3 =3，这样计算得到的速度为 u=1, ＝0。 2 2 二维绕圆柱流动还可以用势函数求解， x2 y2 0 (x, y) 则定解问题可写成 0 边界 cd上 n 0 边界 aec及 bd上 207 n 1 边界 ab上 表示势函数，为了使数值解唯一必须在部分边界