《基于不同时间步长时域非结构有限体积法模拟声-弹性耦合问题 NUMERICAL SIMULATION OF ACOUSTIC-ELASTIC INTERACTION BASED ON MULTI-TIME-STEP UNSTRUCTURED FINITE VOLUME TIME DOMAIN METHOD》.pdfVIP

第34卷第2期 固体力学学报 V01．34NO．2 CHINESE OFSOLIDMECHANICS 2013年4月 JOURNAL 2013 April 基于不同时间步长时域非结构 有限体积法模拟声一弹性耦合问题+ 宣领宽 张文平 明平剑“ 李 川 龚京风 (哈尔滨工程大学动力与能源工程学院，哈尔滨，150001) 摘 要介绍一种改进的时域非结构有限体积法(FVM)，并将其应用于声一弹性耦合问题．在流体与固体介质 中分别求解声波动方程与弹性波方程，根据交界面上的力平衡与质点振速连续条件考虑二者的相互作用．同时考 虑双线性四边形单元的线性变化项及常数项，并结合常应变三角形单元处理混合网格问题．分别对三角形单元和 四边形单元进行色散分析，给出声波动方程的稳定性条件．在不同介质中采用不同时间步长，提高计算效率．求解 弹性波问题、声一弹性耦合问题，结果表明，改进后的方法求解声一弹性耦合问题是有效和准确的，并且具有良好的数 值稳定性． 关键词声一弹性耦合，有限体积法，不同时间步长，时域，非结构化网格 0 引言 耦合问题，并且采用多时间步长提高计算效率，但相 波在流体与固体介质中的传播是非常普遍的现 对于FVM来说，FEM需要求解单元刚度矩阵会导 象，目前，丰要有两类方法模拟其传播过程．第一类 致计算量增加，其计算效率较低[11|． 方法是在流体与固体中均采用弹性波方程描述波的 流体区域与固体区域采用不同数值方法求解 传播，在流体中不存在剪切波，这种方法不需要对边 时，由于不同区域之间数据结构不同，数值传递会影 界特殊处理，因此比较简便．文献[1—5]均是采用这 响计算精度，甚至导致数值发散．因此在不同区域采 种方法求解声一弹性耦合问题．但Choi等[51发现这用同一种数值方法求解是有必要的．目前求解固体 种方法会在流体中声源附近产生明显的噪声，同时， 应力问题最有效的方法是FEM，而求解声学问题常 在流体中用弹性波方程模拟声波会增加计算量及内 常采用边界元法(BEM)与FEM．另一方面，FVM 存需求．第二类方法在流体中采用声波动方程而在 曾被应用于求解一系列的弹性波问题[12。16I，这为求 固体中采用弹性波方程分别描述波的传播，二者之 解声一弹性耦合提供一种可选择的途径．Baliga等u列 间的耦合通过衔接条件或者区域重叠方法实现．目 提出一种基于常应变三角形单元的控制容积有限元 前，第二类方法得到了更为广泛的应用．Komatitsch法(CVFEM)用于求解对流扩散问题，这种方法本 等[63采用谱元法模拟声一弹性耦合问题，采用力平衡 质上来说是交错的非结构FVM。Liu等口8]改进这种 与法向速度连续作为交界面耦合条件，但这种方法 方法的计算效率并将其应用于地震波的求解． 会导致计算量增加口]．Soares等[81考虑声一声、弹性一Zhang[71引入这种方法用于模拟地震波在流体与固 弹性与声一弹性耦合问题，其中弹性波基于有限元法 体中传播．文献E7，18]都引入了四边形单元，但双线 (FEM)求解，声波基于有限差分法(FDM)或FEM性的四边形单元均被处理为常应变单元，这样不仅 求解，两者的耦合通过计算域的重叠实现，但不同的 不能发挥出四边形单元的优势，而且会引进新的误 区域采用不同的离散方法，数值传递会影响计算精 差甚至导致错误．本文考虑四边形单元的双线性特 度．Lee等[93提出一种基于单元中心的时域FDM求性，改进这一缺陷． 解声一弹性耦合问题，但FDM对于不规则区域的适 声