The governing equations are discredited by using Galerkin method, and the matrix for elastoplastic constitutive model is also deduced.

• 采用加权残值法对耦合方程组进行有限元离散，并推导相应的弹塑性矩阵。