论文部分内容阅读
结构与流体相互作用的声固耦合问题大量存在于汽车舱室、空间飞行器以及建筑声学等领域。本文研究并给出了一种基于非结构化网格的时域有限体积法,并将其应用到声固耦合问题中,本文主要研究内容如下:1、声固耦合是两种介质的耦合,因此处理多区域网格是研究声固耦合问题的前提。本文解决了前处理程序中非结构化多区域网格的读入问题。本课题前处理程序调用CGNS-API中级程序库读取CGNS格式网格文件网格信息。程序通过实数HASH表法可以查找到相邻区域中的对应点,再由对应节点组成内部连接面。三维非结构多区域网格读入问题的解决具有很重要的意义。使用本课题开发的前处理程序可以实现流体域和固体域的耦合分析,可以保证压力、速度、流量等在流固交界面上的连续。2、本文给出了一种时域非结构化网格有限体积法解声固耦合问题。运用FORTRAN 90语言编程,模拟流场中的声传播及固体中的弹性波传播,处理流体-固体耦合交界面处边界条件。在算法实现的过程中,本文将声波动方程和位移、应力混合变量弹性波方程各自作为声场、结构的控制方程,采用基于非结构网格的数值求解方法,时间上均采用二阶中心差分格式进行离散,空间上均应用有限体积法进行离散,推导出数值模拟方法。在解固体中的弹性波传播问题时,本文基于力平衡方程和胡克定律的相互关系求解位移,采用欧拉显式解法解方程中各项。位移矢量存储在单元节点上,应变张量和应力张量都存储在单元中心上。与其他方法相比较,有限体积法有着良好的守恒性,并且适用于不规则边界形状。3、本文模拟了分别加载应力驱动和声压驱动时二维流体-固体耦合情况下波的传播。本文分别模拟了二维固体中弹性波传播情况,二维静止流体与固体耦合时波的传播情况,水与空气耦合时波的传播特性。这些算例的计算证明了时域非结构化网格有限体积法解声固耦合问题是可行的。