流固耦合方程的建立及耦合数据传递.pdfVIP

第八届全国现代结构工程学术研讨会 流固耦合方程的建立及耦合数据传递 苏波钱若军 (同济大学土木工程学院．上海200092) 提要：流固藕合问题的难点在于其方程的建立及其数据传递过程。本文从基本理论出发，建立了流固耦合有限元方程。从方程出 发。可以对流固耦合数据传递进行分析．从而建立有效的耦合途径．本文为流固耦合的实现提供理论和方法指导． 关键词：流周耦合，数据传递．有限元，边界 一、流体、固体物理域及其边界 对于流固耦合系统·由不同物质形态的质点所构成的空间可分为两部分：流体域Q，和固体域Q．·图1。 ．r： 图1流体域和固体域固不同边界 流体域Q，，边界r，，考虑耦合边界： r，=r；ur；ur；1 式中，r：-rf分别表示流体域的Dir／chletNeumann边界：r：1为流囿耦台边界a 固体域Q，，边界r，，考虑耦合边界： r，=r：ur；ur77 式中，r二．E分别表示固体域的DirichletNeumann边界；r：”为流固耦合边界。 二、流体运动控制方程 流体在Euler坐标系下的动量平衡和连续方程为； 工业建筑2008增刊 第八届全国现代结构工程学术研讨会 II +mV函+ Q7 p ， xl(3) 址一新 、●，●，／ V(x。t)e = 趾 O V(x。t)e口7xl(4) 采用ALE坐标建立动量守恒方程及连续方程¨·2】： xl(5) ．p㈢小心胁卜州 V(i，f)∈Q7 Vim=O Vfi，t)e口7xl(6) { 式中F角标善表不为ALE坐标糸；p为流体的罾厦；f为流体单位体积力：Ⅳg为参考域上网格点的速厦，仃为 柯西应力张量，用下式表示： 。 o=一pI+2雄 (7) 式中； ： (8)、’ 寥=寺(V，+V；口7)，’、 j J ， I为二阶识别张量；∥为流体的粘性系数；P为流体的静水匿力。 将式(7)和式(8)代入(5)、(6)，可展开写为：’ 娶+搴+妻：0 (9)