二维平面 UEL单元在Dynamic动态分析中的疑问

1454008070@qq.c

关于RHS，在进行ABAQUS的动态分析时，因为会用到NEW-mark方法因而与常规静态的RHS=AMATRX*U求解方式不同，涉及到该公式见下：
file:///F:\tencentfiles\1454008070\Image\C2C\3Z5X66Z1~KIR616%$`7BAGU.png，那么关于这里的G应该怎么去求解呢？是有其他公式进行替代还是需要利用有限元的方程进行内力矢量的求解。希望广大好友和专家们不吝赐教，另外有需要学习UEL的初学者朋友的可以加入QQ群：168360320，免费获得一些基础的FORTRAN教学资源和UEL算例。非诚勿扰！