多场耦合问题求解策略研究.pdfVIP

多场耦合问题求解策略的研究 宋少云，李世其 （华中科技大学机械学院，武汉，湖北，430074） 摘要：对多场耦合问题从概念、算法和设计平台方面进行了简要总结，阐述了该问题的两种 求解策略：分析的求解策略和综合的求解策略，并分析了其优缺点。在此基础上提出了多场 耦合问题的协同设计策略，该策略使用客户/服务器结构，具有信息集成、过程集成和知识 集成的特点。其核心是使用节点信息作为场间信息交换的媒介，并通过节点载荷的映射与插 值完成不同数值方法之间的载荷转换。 关键词：多场耦合，协同设计，求解策略 1 引言 多场耦合(multiphysics problem)是由两个或两个以上的场通过交互作用而形成的物理现 象[1]，它在客观世界和工程应用中广泛存在。常见的耦合问题有热－结构耦合问题、流固耦 合问题、磁－结构耦合问题、静电－结构耦合问题等。根据耦合的空间属性，耦合问题可分 为域耦合和边界耦合两种类型[2] 。域耦合指在整个或部分域内多场共存，各场间没有边界， 如热应力问题；边界耦合指各场间有明显的边界，场之间通过边界作用来实现耦合，如典型 的流固耦合问题。 基于微元分析法所建立的单场数学模型通常是以时间和空间坐标为自变量的偏微分方 程。数学上求解之主要有四种方法：解析法、半解析法、数值法和实验法，在工程领域广泛 使用的是其中的数值法。数值法又分为两类：有网格法和无网格法。根据网格划分的方式， 有网格法又进一步分为有限元法、边界元法、有限差分法、有限体积法和有限分析法。有限 元法适合处理几何形状和边界条件复杂的固体力学问题；边界元法计算量较小，适合处理无 限域问题；有限差分法发展成熟，在各种规则场中广泛使用；有限体积法[3]能保持体积内的 守恒特性，在流体动力学中得到广泛应用；有限分析法保证局域内连续可微，适合计算函数 导数。无网格法[4]则不需要借助于网格，在一系列离散点上建立近似函数，适于处理特大变 形、奇异性或裂纹动态扩展等问题。上述方法针对不同的对象各有其优缺点，都在发展中且 有扩大其应用领域的趋势，因此很难一概的判断孰优孰劣。 数值方法的实施必须依赖数值仿真软件。在多场耦合方面，国际上用得比较多的软件是 ANSYS 、PHYSICA 、FEMLAB等。ANSYS软件最初为解决固体力学和结构力学问题，在其 更新版本中陆续加入了对流场、声场、热场、电磁场的仿真功能，并构造了一系列多场耦合 的仿真算法。FEMLAB[5]最初是作为一个在交互环境中进行基于等式的多物理场建模而出现 的的仿真工具，1999 年发布第一版，最新版本是FEMLAB3.0 。该软件以建模简单为首要目 的，无缝的集成了数学仿真软件MATLAB ，可求解声场、扩散、电磁场、流体力学、结构 力学以及偏微分方程组的耦合问题。上述二软件均使用有限元法，而PHYSICA[6]则是一个 基于有限体积法的多场求解工具，它提供了四种层次的抽象：模型层、控制层、算法层和工 具层。用户可以方便的与前三层进行交互，从而可以自己编制算法，对于科学研究十分有用。 也有一些研究者针对特殊耦合问题进行研究并编制了专用多场仿真软件[7-9] 。 2 分析的求解策略和综合的求解策略 求解多场耦合问题有两种截然不同的策略。一种是首先对各个场建立其数学模型，然后 使用计算机进行数值求解，以得到系统的局部属性，本文称为分析的求解策略；另一种则通 过对耦合场的适当简化，建立系统的宏观模型并分析此模型，直接得到系统的总体属性，本 文称为综合的求解策略。 2.1 分析的求解策略 针对每个场建立各自的数学模型并联立求解是这类求解策略的主要特点，通常该法建立 [2] 的数学模型是代数微分方程组。对于此耦合方程组有三种求解方法 ：域消除法、集成求解 法（直接法）和分区法（迭代法、间接法）。域消除法使用积分变换或者是模态缩减法以消 除一个或几个域，余下的域则用集成求解法。集成求解法在一个时间步内同时求解所有的方 程，而分区法则在某一个时间内依次求解各个场，并在场之间传递耦合载荷。