一种流体结构耦合计算问题的网格数据交换方法.pdfVIP

第八届全国现代结构工程学术研讨会 一种流体一结构耦合计算问题的网格数据交换方法 于猛’钱若军 (同济大学建筑工程系．上海200092) 提要：在流体和结构耦合分析中，结构常常会产生大的位移，为了保持结构和流体两者网格的相容性，要求在计算中对流体及固 体结构界面的不匹配网格进行数据交换。本文介绍了一种流体一结构网格数据交换的方法。 关键词：流固耦合。cED／CSD。数据交换，插值 一、引言 在流体(空气、水流等)作用下，结构(特别是柔性结构)的位移、速度和加速度响应会很大，对周围流场 会产生较大的影响，而流场的改变会进一步改变作用在结构表面上的流体力大小，从而形成流体与结构的相互耦 合作用，即流体一结构耦合作用(Fluid．StructureInteraction，FSI)。 在流体一结构耦合分析中，流体计算更关注与结构所接触表面周围的流动区域，而结构计算中更关注作用在 结构表面的荷载及荷载对结构内部造成的影响。因此，在进行网格划分时，流体和固体结构的网格密度将会不同， 流体网格密度在物体表面处需要增大，而结构动力学要求固体结构表面网格尽量划分均匀，以便方便求出刚度矩 阵．这就会导致划分后的网格不能完全重合。另外，即使最初的网格划分一致，随着流体一结构相互作用的不断 变化，也不能保证流体一结构界面上的流体网格和固体结构网格完全重合，这就要求在计算中对流体及固体结构 界面的不匹配网格进行数据交换，借此来维持流体和结构网格的相容性。 由此可知，要实现流体一结构耦合计算，如何设计两个网格系统的数据交换界面是很重要的，即寻求一种方 便、高质量的插值方法，将计算结构动力学得到的变网格的位移插值到流体}c9格上，并将流体网格上的荷载插值 到结构网格结点上。 二、一种流体一结构耦合问题数据交换方法 (一)简介 典型的流固耦合闯题可以用以下的图1表示。流固耦合模型主要的任务是通过耦合计算得到流体和固体的响 应。结构模型基于拉格日坐标系统，以位移为未知量，流体域采用欧拉坐标系。然而，对于流固耦合问题必须以 任意拉格朗日坐标系统为依据，因为流固边界存在变形。因此，对于流体计笪蛮詈除诵常的变量(压强，速度等 等)还包括位移。 工业建筑2008增刊 6ll 第八届全国现代结构工程学术研讨会 觯＼ 图2流体一固体结点耦合 流体节点位移根据固体节点位移插值计算。例如，节点2由节点1和节点2插值计算。根据边界节点位移，其它 流体节点位移可以采用一定的方法计算以保持初始网格的质量。 同样，固体节点的流体力也采用边界处流体单元节点的应力插值计算。如图2所示，固体节点≥的流体应力 采用流体节点2和3的内插，固体节点1和3的流体牵引力分别等于流体1和4的应力。 根据公式(1)。流体对固体节点2的力包括在固体节点1。2和3的应力。因为流体应力是压强和剪应力的 · 合力，流体求解变量和固体求解变量由此完全耦合。 ． 因为允许流体和固体域可以分别划分网格，因此在界面上两种网格可以一致。相|一几何的离散对于两种模型 也可以不相同。然而，两种离散边界的距离必须在小的范围。这里定义一个从FSI流体节点到FSl固体节点的相 对距离。 。=懈㈥ ㈤ 这里当为从流体节点到固体离散边界的距离，口为固体边界单元的长度(见图3)·这个相对距离一直检查， 当0享：l，程序给出错误并中止；当o．001(_0≤1程序给出警告信息但程序继续。 流体 几何边界 固体模型 FSl界面