关于vuel自定义单元如何传热问题

使用abaqus模拟传热过程，图中为两个C3D8T单元，中间叉号位置为vuel子程序写的自定义单元，这三个单元是在一个part中，如何才能将下方高温单元的热量通过vuel自定义单元传递到上方的低温单元呢？

目前查到的资料是像位移自由度对应的是节点力，那么相应的 温度自由度对应的是节点热流，所以需要在vuel中的rhs中更新节点热流，这个节点热流该如何计算呢？以及这种方法是否正确？