我有一个 2D 方形网格和一个 3D 梁形网格，我想在其中传播地震波。我正在尝试使用开源 FEM 代码（fenics）来模拟它们。我将顶部表面留空，因为这是我们会看到响应（地面运动）的地方。除了底部的一到零狄利克雷条件外，我还设置了其他边界。这些将导致撞击它们的波浪反射回来。通常，这些会吸收边界条件或 PML，但现在我让它们反射。我相信输入地震激发将设置为底部边界（2D / 3D）上的狄利克雷条件。在文献中，Ricker 小波脉冲被认为适合表示地震输入。

我的问题是：

1. 输入源是作为底部边界上的狄利克雷条件给出的，还是作为网格内靠近底部边界的某处的点源里克小波给出？

2. 如果是点源ricker小波；我将如何在我想要的地方使用 c++ 或 python 在 fem 网格中实现它？有没有人有一段代码会有所帮助？是否添加为源术语 '$f$' 在里面 '$f∗v∗dx$' 波动方程弱形式的一部分（其中$v$是测试函数）还是将其投影为网​​格上的初始条件？

非常感谢您的指导！