abaqus 焊接温度场-组织场-应力场子程序编写问题

小白本人要做焊接焊接方面的模拟，但是对子程序一头雾水，想做的方面大概了解到必须用dflux、usdfld和uexpan三个子程序联用在一起，dflux提供热源；usdfld来求组织分布；uexpan来模拟相变时不同组织转变体积的变化。但是具体怎么来写还是不清楚，各位专家有没有类似的代码或者大概的一个流程能帮到小弟，让我有个思路，万分感谢！！