流固耦合问题的研究一直是众多科研人员的研究热点,商业软件公司也一直想把流固耦合的最新求解技术融入到产品中,随着计算技术的发展和流固耦合理论的不断成熟,使流固耦合的求解更加趋于明朗。在过去几十年里,计算流体力学(CFD)和计算结构力学(CSM)都取得了长足的发展,CFD技术已发展到全三维非定常粘性数值模拟阶段,非线性有限元分析法已相当成熟。流固耦合理论也从早期的线性化已经发展到了非线性阶段,对于强迫响应问题,Sayma A. I 等人早已提出了非线性的流固耦合数值模型,即便理论和技术方面都有所突破,但要真正完全求解流固耦合问题还相当困难。本文对流固耦合的分析侧重于耦合作用仅仅发生在流体和固体的交界面上,其耦合效应是通过在方程中引入两相耦合面边界条件的平衡及协调关系来实现的,而对两相域难以明显分开、部分或全部重叠的物理现象不在本文探讨范围内。
平板在流场中振动是一个经典的流固耦合算例。如图1所示,在12m ×0.06m ×6m 的流场控制域中垂直放置一个1m ×0.4m ×0.06m 的平板,流场控制域四周封闭,前后表面设置为对称边界条件,平板下端全约束,在平板和流场控制域相接面为流固耦合面。在平板左侧施加大小为100 Pa 平面压力使得平板向右侧弯曲,0.5 s 后该作用力释放,平板左右震荡试图回到垂直的位置。由于流体的阻尼的存在,使得平板振幅趋于减缓,最终平板将回到垂直位置。已知平板的材料参数为密度2550 kg / m3 ,弹性模量:2.5e06 Pa ;泊松比:0.35;假定流体密度为1 kg / m3 ,动力粘性系数为0.5 Pa ⋅ s 。平板和流体域网格用六面体单元划分,得到平板网格数为3000个,流体控制域单元为13.4万个。设置计算时间为25s(time,25.0),固定时间步长为0.1s。