abaqus UEL 二次开发

本人研究非饱和土方向，需要编写一个自由度为位移、孔隙水压力和孔隙气压力的UEL单元，采用隐式动力分析步，自己也编写了UEL的大致框架，控制方程的形式为我上传的图片，我有相关的理论论文