基于谱元法的弹性波传播数值模拟研究报告.docVIP

基于谱元法的弹性波传播数值模拟研究报告

一、谱元法的核心原理与技术框架

谱元法（SpectralElementMethod,SEM）是20世纪80年代末由Patera提出的一种高精度数值计算方法，它融合了有限元法的几何灵活性和谱方法的高精度特性，在弹性波传播模拟领域展现出独特优势。其核心思想是将计算区域划分为若干非重叠的单元，在每个单元内部采用高次多项式（通常为拉格朗日插值多项式）近似场变量，通过变分原理建立离散化的控制方程。

在弹性波传播问题中，控制方程基于牛顿第二定律和胡克定律推导而来。对于各向同性弹性介质，位移场(\mathbf{u}(\mathbf{x},t))满足的运动方程为：[\rho\frac{\partial^2\mathbf{u}}{\partialt^2}=\nabla\cdot\boldsymbol{\sigma}+\mathbf{f}]其中(\rho)为介质密度，(\boldsymbol{\sigma})为应力张量，(\mathbf{f})为体积力向量。胡克定律将应力与应变联系起来：[\boldsymbol{\sigma}=\lambda(\nabla\cdot\mathbf{u})\mathbf{I}+2\mu\boldsymbol{\varepsilon}]其中(\lambda