我正在尝试使用Fenics Solid Mechanics Application模拟弹塑性地震波传播。

该应用程序。提供了一些准静态演示来展示立方体/梁/正方形中的弹塑性行为。他们实现的弱形式取自开发者自己写的这篇论文的第 199 页，给出如下：

所有术语都在本文等式下方的一段中定义。该方程与线性弹性获得的弱形式非常相似。必须注意的是，在计算可塑性中，返回映射算法（Simo 和 Hughes 1998）很流行，因此该应用程序使用“最近点投影”返回映射算法以及牛顿法。线性形式“L”的应力 (sigma) 是残余应力（或来自先前迭代的应力）。在 LHS（或双线性形式“a”）中，变量“C”是一致的弹塑性切线模量。其他变量是体力'f'；牵引'h'; 测试函数'v'；和试用功能'u'。

到目前为止，一切都很好。现在，对于弹性线性波传播，实现的弱形式是：

这里，测试函数“w”与等式（6）和（7）中的测试函数“v”相同。唯一的新术语是“intertia term”，其中$/rho$是材料的质量密度和$/ddot{u}$是加速度。

考虑到上述两个方程，可以说我可以使用以下弱形式来实现弹塑性波吗？

这里$C_{ep}$是相同的$C$在等式（6）中。