基于球面波势函数基本解方法的弹性波三维散射与动应力求解

针对弹性波三维散射和动应力集中问题,提出一种基于球面波势函数的基本解方法(SWP-MFS)。方法基于单层位势理论,采用膨胀波和矢量剪切波球面波势函数构造散射波场,根据边界条件建立边界积分方程并配点求解。精度检验表明,该方法具有良好的数值精度及稳定性。以无限空间中三维夹杂体及空洞对平面P、SV波的散射为例,进行方法展示,并揭示了三维夹杂体周围弹性波散射的一些重要规律。结果表明:三维夹杂体随其内部介质刚度降低,位移谱曲线震荡越加强烈。三维球形空洞在P波和SV波水平入射下应力集中规律不同,前者在顶部和底部及附近更明显,后者在纵截面两45°角线附近更明显。与以往的集中力源函数相比,新的波场构造基本解更为简洁易用,为三维弹性波动分析提供了一种新型无网格数值方法。

弹性波三维散射快速多极子间接边界元法求解

边界元方法对于无限域中弹性波散射求解具有独特优势,但求解矩阵的非对称稠密特征极大限制了该方法在大规模实际工程中的应用。为此,基于单层位势理论,结合快速多极子展开技术,通过对球面压缩波和剪切波势函数的泰勒级数展开,建立一种新的快速多极间接边界元方法,以实现大规模弹性波三维散射的精确高效模拟。算例分析表明所提方法能够大幅度降低计算时间和存储量,可在目前普通计算机上快速实现上百万自由度弹性波三维散射问题的快速精确求解。最后以全空间椭球形孔洞群对平面P波、SV波的散射为例,揭示了三维孔洞群周围稳态位移场和应力场的若干分布规律。该文方法对低无量纲频率(ka<5.0)的大规模多体散射问题尤为适合。