一种用于动力土-结构相互作用分析的混合SEM-MTF-DRM方法

《Soil Dynamics and Earthquake Engineering》：A hybrid SEM-MTF-DRM method for dynamic soil-structure interaction analysis

　　

地震工程中，对动力土-结构相互作用问题进行准确且高效的数值模拟仍是一项挑战，其主要困难在于对土体域的处理。研究人员通过结合谱元法（SEM）用于空间离散、研究人员先前提出的优化多透射公式（MTF）用于人工边界条件、以及域缩减方法（DRM）用于模型内部地震波输入，开发了一种用于动力土-结构相互作用分析的混合SEM-MTF-DRM方法。研究详细阐述了在SEM模型中实现优化MTF边界条件和DRM地震波输入的过程，从而为该集成方法提供了清晰的算法流程。数值测试验证了该方法在宽频率范围和入射角下，对P波和S波入射均具有高精度和鲁棒性。研究还发现，由于DRM将地震波输入与边界解耦，并能在内部域中抵消部分出射波，相较于传统将地震波输入与边界耦合的方法，MTF边界的精度和稳定性均可进一步得到提升。对实现问题的研究表明，当DRM层与MTF边界保持一定距离时，可获得更好的结果。计算效率主要取决于MTF边界的位置（即模型尺寸），而受DRM位置的影响较小。在某些特殊模型，如具有深埋结构或空间分离局部特征的模型中，可采用不同形状或数量的DRM层来简化模拟。

**论文解读：一种用于动力土-结构相互作用分析的混合SEM-MTF-DRM方法**

**研究背景与问题**

动力土-结构相互作用（SSI）或复杂局部场地中地震波传播的数值模拟，是工程抗震设计与地震风险分析的关键工具。与纯结构动力分析不同，土壤或地质构型的参与使得模拟在空间离散方法、人工边界条件（ABC）和地震波输入技术三个方面面临显著挑战。对于空间离散，有限元法（FEM）因其成熟的框架和丰富的商业软件而广泛应用，但谱元法（SEM）凭借其在波形复现和集中质量矩阵自然推导方面的优势，在地震工程领域日益受到关注。在人工边界方面，多透射公式（MTF）作为一种外推型ABC，具有简单、通用且计算成本低的优点，但传统上在FEM中与地震波输入耦合时存在数值稳定性问题（如低频漂移和虚假高频振荡）。尽管将MTF应用于SEM可改善稳定性，但问题仍未完全解决。此外，地震波输入技术是另一个难题：传统方法常通过边界施加波场，但这对复杂ABC（如完全匹配层PML）实现困难，且可能降低边界稳定性。域缩减方法（DRM）作为一种从边界解耦的内部波输入技术，通过在模型内部施加等效节点力，简化了流程，但DRM要求获取输入层的刚度（及质量、阻尼）系数，这在封闭代码中可能受限。因此，开发一种集成SEM高精度、优化MTF的简易性以及DRM解耦优势的混合方法，对于提升SSI模拟的精度、稳定性和效率具有重要意义。

**研究内容与结论**

研究人员提出了一种混合SEM-MTF-DRM方法，用于动力土-结构相互作用分析。该方法将先前提出的优化MTF（即具有多个计算波速的MTF）与SEM和DRM相结合，并详细阐述了优化MTF边界条件和DRM地震波输入在SEM模型中的实施流程，形成了清晰的算法工作流。通过三个基准算例（半圆形峡谷、半圆形山丘和半空间圆形孔洞）的数值测试，验证了该方法在宽频率范围和入射角下对P波和S波入射均具有高精度和鲁棒性。进一步研究发现，由于DRM将地震波输入与边界解耦，并能在内部域中抵消部分出射波，MTF边界的精度和稳定性相较于传统耦合方法得到显著提升。对于实现问题的参数研究表明，当DRM层与MTF边界保持一定距离时效果最佳；计算效率主要取决于MTF边界位置（即模型大小），而受DRM位置影响较小；在某些特殊模型（如深埋结构或空间分离局部特征）中，可采用不同形状或数量的DRM层来简化模拟。该论文发表在《Soil Dynamics and Earthquake Engineering》。

**主要关键方法**

本研究涉及的主要关键技术方法包括：（1）谱元法（SEM）用于波动传播的空间离散，采用基于Gauss-Lobatto-Legendre（GLL）节点的Lagrange插值和相应的高精度数值积分，形成高阶集中质量单元；（2）优化的多透射公式（MTF）用于人工边界条件，通过引入多个计算波速提高宽频范围精度；（3）域缩减方法（DRM）用于模型内部的地震波输入，通过在仅一个单元厚的输入层上施加等效节点力，实现与边界的解耦。所有模拟均为数值试验，不涉及样本队列或试剂操作。

**研究结果**

**1. 混合SEM-MTF-DRM方法的验证**
通过三个基准问题（半圆形峡谷、半圆形山丘和半空间圆形孔洞）验证了混合数值方法的准确性。在第一步中，研究人员利用谱元法（SEM）结合优化多透射公式（MTF）边界高效获取高精度自由场节点位移，用于构建DRM等效节点力；第二步中，采用完整的混合框架进行总场模拟。结果表明，该方法在宽频率范围和不同入射角下（P波和S波）均能获得与参考解高度吻合的响应，验证了其高精度和鲁棒性。

**2. 引入DRM对MTF边界性能的改进**
对比传统将地震波输入耦合在MTF边界上的方法，引入域缩减方法（DRM）后，MTF边界的数值稳定性显著提升。研究发现，DRM通过解耦地震波输入，能够在内部域中抵消部分出射波，从而抑制了传统MTF中常见的低频漂移和虚假高频振荡，使边界精度和稳定性均得到改善。

**3. 若干实现问题的研究**
通过系统参数研究，分析了DRM层位置、MTF边界位置、DRM层形状和多个DRM子层对计算效率和精度的影响。结果表明：当DRM层与MTF边界保持一定距离时，可获得更优的结果；计算效率主要受MTF边界位置（即模型尺寸）控制，DRM位置影响很小；在某些特殊模型（如深埋结构或空间分离局部特征）中，可采用不同形状或数量的DRM层来简化模拟。

**讨论与结论**

研究结论部分指出：通过将先前提出的优化多透射公式、谱元法和域缩减方法相结合，开发了一种用于动力土-结构相互作用分析的混合SEM-MTF-DRM方法。该集成方法兼顾了SEM在波动模拟中的高精度、优化MTF作为人工边界条件的简单通用性，以及DRM地震波输入（从边界解耦）带来的性能提升。数值测试表明，该混合方法在宽频范围和入射角范围内均有效且高精度。相较于传统将地震波输入与边界耦合的方法，DRM的引入显著改善了MTF边界的精度和稳定性，因为DRM可抵消部分出射波。参数研究为实践中DRM和MTF的设置提供了指导：DRM层应与MTF边界保持适当距离以获得最优结果；计算效率主要由MTF边界位置（即模型大小）决定，DRM位置影响可忽略；特殊模型可根据需要调整DRM层的形状或数量。