comsol模拟胶水固化对两个物体造成的变形，求教

参考论文《复合材料层合板固化过程多物理场耦合数值模拟》，用到固体传热、域微分方程和微分代数方程、固体力学三个物理场。仅固体传热与域微分方程耦合时，能够求解出固化度分布与温度分布，以及随时间的变化。但三场耦合时，求解进度卡在4%，然后报不收敛。按照网上的教程，调初始值、加阶跃函数、改求解设置，出现过的问题包括但不限于“找不到一致的初始值”“除零错误”“试图计算负数的非整数次幂”“最后一个时步不收敛”“非线性瞬态求解器不收敛”。求教各位大佬，还能怎么调整？