我尝试实现[1]中提出的光谱元素方法1。我模拟了一块铝板，并在板2左上角的中间施加了一个垂直集中力。我应该得到一个带有圆形波峰的波，来自我施加力的点。但显然波前反转了，如下图所示。
（这是特定时刻的 z 位移。）
我的代码一定有问题。可能发生了什么？
组装刚度和质量矩阵时出了什么问题？此链接 包含我的 pdf 格式的 python 源代码和说明波传播的视频。我会先尝试用普通的板元素重写我的代码。
1这里，SEM 只是具有高阶多项式形状函数的有限元方法。其他一些作品使用正弦形状函数。
2我没有滥用对称性，因为我稍后会在盘子上添加一些东西。
更新：我尝试使用普通的板元素，但结果是相同的，因为网格足够细。
编辑 2：我认为我的代码实际上是正确的，因为 Abaqus 模拟还在主波峰的顶部、底部、左侧和右侧显示了这种圆形图案（尽管它们更模糊）。但我仍然对为什么会出现这些圆形图案感到困惑。我想这与集中负载有关，只有当我细化网格（h-refinement）时，这种虚假图案的效果才会减弱。但根据我的观察，p-refinement 没有帮助。
