增材制造技术中基于时域有限元方法的声波—弹性波耦合(一):理论

摘要

基于时域有限元方法,为增材制造技术中常涉及的声波—弹性波耦合问题构建理论框架。在三维直角坐标系下引入流体介质所满足的关于压力场的时域声波方程和固体介质所满足的关于位移场的时域弹性波方程。将时域波动方程转换到频域进行复坐标系拉伸,再对拉伸后的频域波动方程进行波场分裂。复拉伸因子是表征声波/弹性波衰减的物理量,在声波方程/弹性波方程中体现为近乎于复介声常数/复介弹常数的反映特性,与介声常数/介弹常数、声导率/弹性导率、声波/弹性波传播速度等参数构成了环环相扣的有机整体,但表现形式与电磁波理论中复介电常数的情形不同。提出一种通过引入过渡变量求解卷积的处理策略,将分裂后的频域波动方程转化回时域,得到带有完全匹配层(PML)的时域波动方程,其中声波方程和弹性波方程中最大未知量个数分别仅为7个和27个。推导带有PML的时域波动方程在有限元方法实现中所需的等效积分弱形式,进而基于声波和弹性波在流体和固体区域交界面的双向转换机制,导出带有PML的声波—弹性波耦合方程的等效积分弱形式。为节约PML计算成本,根据波动衰减方位数量的不同,对角点采取分门别类的处理方法,使得通常情况下带有PML模型的计算资源消耗相比模型不带有PML时仅增加1倍左右。波动方程及其PML内在机理的深入探究和计算效率的深度优化,为增材制造技术中涉及声波—弹性波耦合问题的应用提供了前提和保障,而带有PML的声波—弹性波耦合方程的等效积分弱形式,则是促进基于时域有限元方法的声波—弹性波耦合在增材制造技术中得以应用的理论精髓所在。