SPH粒子热边界条件问题？

各位大神，小弟用SPH模拟金属成型热力耦合模拟。图中为建立的一个板，在板的前后、左右、下表面用SPH_SYMMETRY_PLANE约束SPH粒子的位移边界条件（~PS~在论坛及参考书中都说对SPH粒子应该用SPH_SYMMETRY_PLANE约束，而不应该用SPC），模型初始温度为20℃，在分析过程中为什么所有SPH_SYMMETRY_PLANE约束处的SPH粒子温度下降至0℃，应该如何修改？

谢谢