对流-固耦合问题进行拓扑优化的新方法

发布时间：2024-04-29 05:31:51 丨 浏览次数：141

工程师在开发涉及流-固耦合 (FSI) 的系统（例如，桥梁、飞机和生物医学应用）时，可能会希望在设计早期就使用拓扑优化。然而，许多拓扑优化方法仅适用于学术范围内的 FSI 问题。本文，研究人员开发了一种新的方法来帮助解决这个问题……

流固耦合 （FSI） 描述了流体动力学规律和结构力学规律之间多物理场的相互作用。在设计各种系统（例如汽车、飞机和航天器）时，考虑 FSI 尤为重要。如果忽视这些影响，可能会导致灾难。

1940 年 11 月，在华盛顿州塔科马市发生了一起悲剧，Tacoma Narrows 大桥在开通仅四个月后就倒塌了。它在被摧毁时是世界上第三长的吊桥，连接着 Tacoma 和 Kitsap Peninsula 岛。

刚开始，许多人认为该桥梁倒塌是由于共振，但后来我们了解到是由于气动弹性颤振导致的。这种大规模的结构设计失败证明了研究 FSI 的重要性，尤其是在设计大跨度桥梁时。（Tacoma Narrows 大桥后来被重建成一对双悬索桥）。

FSI 在血液流动建模中也发挥着重要作用。例如，生物医学工程师可以通过模拟 FSI 来准确预测人类心脏瓣膜的行为，并解决各种心脏健康问题。

拓扑优化是优化的三个主要类别之一（另外两个是参数优化和形状优化）。

拓扑优化是一个非常强大的工具，因为它可以帮助我们为给定目标和约束条件的工程结构生成新的拓扑。但是，它的用途通常仅限于单物理场设计。

为了解决这个问题，来自圣保罗大学（University of S?o Paulo）和加州大学圣地亚哥分校（University of California San Diego）的研究人员着手开发了一种新方法，将拓扑优化方法应用在结构的 FSI 分析中。

在圣保罗大学，研究人员正在努力推进天然气的可持续利用。由于他们的大部分工作都涉及 FSI 问题，因此圣保罗大学多物理场系统优化实验室开发了一种新方法，即通过 LiveLink? for MATLAB? 接口产品将拓扑优化方法——二元结构拓扑优化法（Topology Optimization of Binary Structures，TOBS）集成到有限元软件包 COMSOL Multiphysics? 软件中。

二元结构拓扑优化方法最近由 Renato Picelli（曾是英国卡迪夫大学的博士后助理）和来自加州大学圣地亚哥分校的 Raghavendra Sivapuram 独立开发。该方法使用了顺序近似问题、敏感性过滤、整数线性规划，以及在处理与设计与多物理场问题时非常有用的二进制设计变量。圣保罗大学的多物理场系统优化实验室团队开发了一个基于几何修整（geometry trimming，GT）的框架，创建了用于流体结构设计的 TOBS-GT 方法。

Shahin Ranjbarzadeh 是圣保罗大学的博士后，同时也是圣保罗大学多物理场系统优化实验室团队成员，他表示：“一开始，二元结构拓扑优化方法只是应用于结构的拓扑优化。通过开发 TOBS-GT方法，我们将这一想法应用在涉及层流、湍流和非牛顿流的 FSI 问题。”这项工作代表了在 FSI 问题中考虑整数线性规划的第一次尝试。

将二元结构拓扑优化方法与 LiveLink? for MATLAB? 结合使用的主要好处是结果中没有灰度，因为二元结构拓扑优化方法基于二元设计变量 (0,1)。“二元结构拓扑优化方法基于 0 和 1，所以在拓扑优化结果中没有任何中间值，” Ranjbarzadeh 解释道。此外，使用 TOBS-GT 方法，我们可以根据物理场要求调整网格。

为了让您更好地了解所提方法的特点，研究人员分享了他们 2020 年发表在 Springer Nature 上的论文中使用的示例，这是他们首次发表该方法的文献。

“当我们写第一篇论文时，我们正在寻找一个真正的应用。最终，研究人员选择展示该方法在 FSI 载荷下降低线弹性系统柔度的能力，所有这些载荷都受到体积限制。”Ranjbarzadeh 说道。所有的示例显示了仿真结果的平稳收敛。

研究论文中的一些示例是基于气动装置的。

Ranjbarzadeh 说：“我认为这种方法将有助于预测气动设备的故障，例如压缩机或旋转机械中的故障。”

除了已发表论文中展示的结果之外，他们还对 FSI 拓扑优化进行了 3D 模拟，如下示例所示。

文中介绍的方法也可用于对其他的多物理场耦合执行拓扑优化。

目前，研究人员正在开展一个新项目，以使用相同的方法开发流体拓扑优化。“我们正在做的事情非常简单。该程序分为两部分，” Ranjbarzadeh 说。首先，他们正在使用 COMSOL Multiphysics 解决 CFD 问题。通过定义几何形状、网格划分并计算灵敏度分析接下来，他们将灵敏度提供给 MATLAB? 中的优化点网格（TOBS 方法是在 MATLAB? 环境中编写的），从这里开始，该团队开始计算和分析拓扑优化。重复这个过程，直到优化问题得到解决。