Aabqus或COMSOL实现热流固耦合

分析煤层气开采工况，力场（损伤dp模型)和+渗流场（煤体为双孔介质模型，考虑基质孔隙系统的吸附扩散行为和裂隙系统的达西渗流定律，后期再考虑水气两相渗流)的耦合计算。1 使用abaqus软件，力场本构的umat fortran子程序已有，请帮我实现与渗流场的耦合计算，顺序耦合或者强耦合(uel单元子程序，同时求解位移和孔隙压力)。2 使用comsol软件，两相流分析由软件内置模块实现，请帮我实现耦合计算模块，特别是损伤本构的开发。