旗·膜的FSI分析
旗·膜FSI的理论基础
现象的物理
模拟旗飘动现象很难吗？
旗或膜的飘动（flutter）是流体-结构耦合的典型问题。轻量而柔软的结构在流动中发生大变形，因此需要同时处理几何非线性和强耦合效应。
支配方程
膜结构的力学如何描述？
忽略弯曲刚性的膜模型情况下，运动方程为
$T$ 是膜张力，$\Delta p$ 是流体压力差，$\mathbf{n}$ 是膜的法向。考虑弯曲刚性时使用Kirchhoff-Love板理论。
其中 $D = Eh^3 / (12(1-\nu^2))$ 是弯曲刚性。
影响旗稳定性的主要参数是什么？
质量比 $M^* = \rho_s h / (\rho_f L)$ 和无量纲弯曲刚性 $K_B = D / (\rho_f U^2 L^3)$ 是主要参数。$M^*$ 越小（旗越轻）越容易失稳，$K_B$ 越小越容易产生大幅度飘动。
线性稳定性分析可以求得临界流速。无量纲临界流速 $U^*_c$ 是 $M^*$ 和 $K_B$ 的函数，可与Shelley等的理论分析及Connell的实验相对比。
旗飘动的"临界速度"——理论预言的飘动开始时刻
在风中旗开始飘动时，跨越某个临界风速后突然出现振动，这就是"扑振不稳定性的出现"。当风速超过临界值时，作用在旗上的流体力超过结构的阻尼，振幅呈指数增长。从理论上看，单位面积质量ρ_s与流体密度ρ_f的比值（质量比ρ_s/（ρ_f·L），L为旗长）很重要，这个值越小，越容易在低风速下开始扑振。反过来说，纸张或薄膜在水中即使流速很低也会飘动——这个特性被用于"利用水流的柔性膜能量采集器"的最新研究，旗的FSI理论被应用于发电技术。
旗·膜FSI的数值计算手法
数值手法的选择
旗的FSI分析采用什么手法？
由于涉及大变形，IB法或overset mesh比ALE法更合适。
| 手法 | 特点 | 代表代码 |
|---|---|---|
| IB法 | 网格固定。对大变形强 | IBAMR, IB2d |
| ALE-FEM | 界面精度高。需要重新网格划分 | Ansys, COMSOL |
| LBM-IBM | 格子Boltzmann+IB | Palabos |
| SPH-FEM | 粒子法。有自由表面 | OpenFPM |
IB法将膜厚度设为零处理吗？
Peskin型IB法将拉格朗日结构作为1维低的流形处理。2D的旗作为1D纤维，3D的膜作为2D面嵌入。力的传递使用Dirac delta函数的正则化版本。
时间积分和稳定性
膜很轻时时间积分会有问题吗？
质量比小的系统中，显式时间积分的CFL条件变得极其严格。需要半隐式方案（流体隐式、结构显式）或完全隐式方案。
根据Causin型稳定性条件，弱耦合的稳定性取决于 $\rho_s h / (\rho_f \Delta x)$。这个比值越小，弱耦合越容易失稳。
用ALE法计算旗时网格会"融化"——格子变形的极限
旗的FSI计算中使用ALE（任意拉格朗日-欧拉）法时，网格会随旗的大变形而变形。但当旗大幅飘动时，旗的后缘附近的网格会极度扭曲，矩阵接近奇异，导致计算发散——这就是"网格坍塌"。为了避免这种情况，采用网格重新生成（重网格划分）或各向同性方法，但频繁的重网格划分会大幅增加计算成本。近年来沉浸边界法（IBM）受到关注，即使结构大幅变形，流体网格仍保持固定，非常适合旗这样的大变形问题。IBM在计算精度上不如ALE，但在稳定性和计算成本上具有优势——手法的选择是设计人员的功夫所在。
旗·膜FSI的实际应用
模型设置的实践
请教我如何设置2D旗飘动问题的基准测试。
Turek基准（FSI2, FSI3）是广泛使用的标准问题。
| 参数 | FSI2 | FSI3 |
|---|---|---|
| Re | 100 | 200 |
| 结构密度 | 10,000 kg/m³ | 1,000 kg/m³ |
| Young模量 | 1.4 MPa | 5.6 MPa |
| 变位幅度 | O(cm) | O(cm) |
| 特征 | 弱耦合可用 | 需要强耦合 |
网格设置应该注意什么？
旗的前后缘和旗附近的尾迹区域要充分分辨。旗厚度方向至少需要4~6个单元的分辨率。使用ALE法时，考虑大变形会导致网格压缩，网格设计时要预留最大变位1.5倍以上的变形余裕。
膜结构的应用：抛物线太阳能集热器
有什么工程应用例子吗？
太阳光集热器的菲涅尔镜或抛物面槽镜膜会因风荷载而变形，降低集热效率。通过求解FSI来确定满足许可变形量的结构设计。其他应用还有充气结构（可展开式太空结构体、帐篷结构）和能量采集用旗（压电旗）。
"旗的噪音太大睡不着"——薄膜FSI导致的噪声问题
建筑工地的隔离膜或工事围栏在强风下发出"啪啪啪"的巨大噪音，这是薄膜扑振产生的噪声。这种"飘拍声"当扑振频率进入可听范围（20~20,000Hz）时特别明显，成为周边住户的投诉。实务中通过给膜施加一定的张力（预张）来提高固有振动频率，避免共振。将预张从0.5kPa增加到2kPa，扑振开始风速就从7m/s改善到12m/s的实验数据存在。帐篷结构和户外广告幕的设计中也会出现同样问题，建筑业界已积累了薄膜FSI的实用诀窍。
旗·膜FSI的软件比较
工具选择
旗和膜的FSI分析用什么工具最好？
对大变形的支持能力是工具选择的关键。
| 工具 | 大变形支持 | 膜单元 | 特点 |
|---|---|---|---|
| Ansys Fluent + Mechanical | ALE + remeshing | 壳/膜 | System Coupling耦合。工业应用有实绩 |
| STAR-CCM+ | 变形 + overset | 壳（内置FEA） | 自动重网格优秀 |
| COMSOL | ALE | 膜/壳 | 小规模单一耦合 |
| OpenFOAM + preCICE | IB法/ALE | CalculiX/FEniCS | 开源。最适合研究应用 |
| LS-DYNA | ALE/SPH | 膜/壳 | 显式法。含冲击的FSI |
如果在研究中要解Turek基准用哪个好？
preCICE + OpenFOAM + deal.II的组合在研究社区被广泛使用。preCICE的教程包含Turek基准的设置文件，从这里开始是最高效的。
工业应用怎么样？
Ansys的System Coupling在稳定性和易用性的平衡上表现良好。膜单元使用Ansys Mechanical的SHELL181（6自由度壳单元），启用大变形选项（NLGEOM, ON）。
能否用OpenFOAM做旗的FSI？——工具选择的现实
"我想用OpenFOAM计算旗的FSI"这样的咨询经常来自FSI初学者，但实际上门槛很高。OpenFOAM的固体力学（solidDisplacementFoam等）以小变形为前提，不适合旗这样的大变形大挠度。实务中通常采用OpenFOAM的流体求解器与CalculiX（开源FEM求解器）配合，使用preCICE这个开源耦合中间件连接的方法。"OpenFOAM + CalculiX + preCICE"的组合在旗、帆、窗帘等大变形膜FSI计算中被广泛使用，扑振开始风速的预测精度与风洞试验相比在±15%以内。自由工具也能做本格的旗FSI的时代已经来临了。
旗·膜FSI的先进研究
能量收集
听说有从旗的飘动中取出能量的研究。
在柔性旗上贴压电素子（压电旗），将流体能量转换成电能。需要在FSI中加入压电效应方程进行耦合。
$\{D\}$ 是电束密度，$[e]$ 是压电应力常数，$\{S\}$ 是应变，$\{E\}$ 是电场。
走向实用化的课题是什么？
转换效率低（目前不足1%）和疲劳寿命问题。不过作为低功耗物联网传感器的电源有小众应用可能。
织物非织造布的微观FSI
过滤器和织物的空气阻力也用FSI分析吗？
微观上是纤维间流和纤维变形的耦合。使用体积平均Navier-Stokes方程与Darcy定律相结合的多尺度FSI。还被应用于降落伞伞包的透过率评估。
等几何分析（IGA）的应用
膜FSI中使用IGA有什么优点？
使用NURBS基函数，曲面表现精度高，与薄膜壳模型兼容性好。CAD-CAE间无需网格转换，与形状优化循环的整合也容易。Bazilevs研究室的手法是先驱，被应用于风电机叶片和旋转球帆的FSI。
鱼为何游得不累——旗FSI与推进力的秘密
鱼在水中游泳时，尾鳍的运动遵循与旗相同的FSI机制。鱼主动利用肌肉同时巧妙地利用流体的反馈力来最大化能效。MIT的研究团队在2015年发表的分析表明，形似金枪鱼的柔性膜结构"与卡门涡的位相同步振动"，与刚性尾鳍相比可以将能耗降低多达35%。这个发现直接指向水下机器人（仿生UUV）的设计，旗FSI的先端研究已发展成"如何聪明地利用流体力"的工程课题。
旗·膜FSI的故障排除
网格变形导致计算破裂
旗的位移很大时网格压缩会导致计算停止。
这是ALE法的极限。解决方案如下。
1. 增加重网格频率: Ansys Fluent中设置remeshing interval为较小值
2. 使用overset mesh: STAR-CCM+或Ansys Fluent的overlapping grid
3. 改用IB法: 网格固定，对大变形兼容
4. 加强网格光滑: 调整网格平滑的弹簧常数
膜的初始形状问题
如何设置膜在初始状态因自重下垂的形状？
需要进行预应力分析。首先在结构单独上进行自重加初始张力的静力学分析，将变形后的形状作为初始形状引入FSI分析。Abaqus中用IMPORT功能，Ansys中用ICSTATE实现。
膜数值翻折
膜会出现非物理的折叠或自我穿插现象。
是膜单元的自接触判定不足造成的。对策如下
- 添加自接触条件（self-contact）
- 给膜设置微小的弯曲刚性（使用壳单元而不是纯膜单元）
- 添加稳定化项（Baumgarte稳定化等）
LS-DYNA中可用*CONTACT_AUTOMATIC_SINGLE_SURFACE设置自接触。
"旗只向一个方向飘动"——初始条件非对称性引发的计算bug
旗FSI计算中出现"旗只向一个方向振荡，无法再现实际的飘动"症状而被卡住的人很多。原因几乎总是初始条件的对称性过于完美。旗从完全平展（位移为零）开始计算时，流动也完全对称，旗向上还是向下的"决定因素"消失了。理论上对称系统会因微小数值噪声而破坏，但计算机的舍入误差有时不足以产生足够的非对称性。对策是"给旗施加约0.01L（旗长的1%）的初始挠度"。仅此而已就能让计算马上稳定。风洞试验中旗也从"微小的非对称"开始飘动，所以从物理上讲这也是正确的处理。
相关主题
价值
详细
错误