基于无网格法的连铸结晶器内坯壳凝固与热应力的数值模拟

• 作 者: 张雷
• 出 处: 清华大学
• 发表日期: 2005年期
• 关键词: 连铸 凝固 热弹塑性应力 无网格方法 无网格迭代法

无网格方法是一种新的数值计算方法,本文对该方法在连铸结晶器内坯壳凝固与热应力分析中的应用进行了研究。论文研究了径向基函数法和有限点法两种配点型的无网格法在分析实际传热与凝固问题时的实现方法。根据经典的高斯迭代法和松弛迭代法,通过对径向基函数法的求解平衡方程进行变换,提出了迭代求解算法。研究表明:(1)迭代法可以节省内存,适合于求解大规模工程问题。(2)有限点方法的实用性较好;(3)有限点方法在求解二维问题时表现良好的权函数影响域放大系数ac,在求解其它类似问题时也能获得较好的结果。以上结果为无网格方法在连铸结晶器内坯壳凝固与热应力分析中的应用提供了依据。基于无网格有限点方法实现了求解凝固导热问题的算法,并用Onate 的稳定性方法处理了Neumann 边界条件,给出了凝固潜热的求解方程式和实施方案,并且考虑了材料非线性。研究表明:对于Neumann 边界条件的处理,能提高计算精度。用经典的凝固问题验证了算法的正确性。在求解弹性静力学方程的无网格局部迦辽金法的基础上,给出了其求解热弹塑性问题的位移求解平衡方程式;采用增量切线刚度法处理了热弹塑性材料非线性,给出了相应的方程式及实现步骤;通过对分片试验、定向拉伸以及等温膨胀等问题的数值分析,验证了该算法的可靠性。作者开发了连铸结晶器内铸坯热-力耦合无网格分析程序,并对小方坯连铸结晶器内铸坯凝固与热应力进行了分析,结果表明:(1)铸坯内不同部位冷却不均匀造成偏角处塑性应变集中,容易发生变形;(2)气隙的产生降低了坯壳角部热流,使得坯壳厚度薄;(3)气隙产生以及随后的扩张,加剧了坯壳内冷却不均匀,使得偏角处产生较大的应力和塑性应变。另外大方坯实例分析显示结晶器内凝固坯壳厚度的无网格计算结果符合测量结果。论文研究表明,无网格方法适合于连铸结晶器铸坯热-力耦合分析,计算结果符合连铸方坯偏角区域缺陷的形成规律和角部内层应力应变分布特征,同实际相符;分析结果同有限元方法一致。