锆板轧制模拟的显式动力和隐式静力有限元法比较.pdfVIP

第 54卷 第 3期 复 旦 学 报 (自然科学版) Vo1．54No．3 2015年6月 JournalofFudanUniversity(NaturalScience) Jun．2015 文章编号：0427—7104(2015)03—0365—08 锆板轧制模拟的显式动力与隐式静力 有限元法比较 徐宏彬 ，丁淑蓉 (复旦大学 力学与工程科学系，上海 200433) 摘 要 ：针对用作核燃料元件金属包壳的锆合金板建立三维有限元模型，分别运用显式动力与隐式静力有限元 法实现了轧制加工的数值模拟．对两种模拟方法下轧制后工件的几何尺寸、接触面上的轧制压力分布等计算结 果进行了分析与比较．研究结果表明：通过选取合理的质量放大因子，显式算法的计算效率远高于隐式算法；而 隐式静力方法的计算结果比显式动力方法更加平滑；在轧板中部和远离接触表面处 ，两种方法所得轧制后工件 的尺寸吻合得很好；并且两种方法所得稳定轧制阶段的工件表面接触压力分布总体上较为吻合．所做工作可为 核燃料元件轧制数值模拟方法的合理选取提供参考． 关键词：轧制 ；塑性加工；数值模拟；显式动力；隐式静力 中图分类号：TB125 文献标志码：A 板型核燃料元件 由核燃料芯体和包壳组成，通过轧制加工方法制造而成[1]，在研究试验堆 中得到广 泛应用．其加工质量与其堆内使用性能密切相关．锆板被用作制造板型核燃料元件的复合坯件盖板和框 架 ，其亦先由轧制工艺加工而成，再与芯体装配成多材质复合坯件，然后经多道工序的冷、热轧，最终达 到冶金结合．轧制加工可以使其组织性能得到改善，并获得理想的外形尺寸．如何优化轧制工艺参数是提 高板型核燃料元件加工质量所关心的一个重要问题． 随着计算机技术的发展，有限元数值模拟_l3]方法与工艺实验[6]相结合，成为工艺过程优化的新趋 势，并在塑性加工领域得到广泛使用_7]．在轧制模拟的有限元方法中，可按所采用的算法不同而分为两 类[9=：基于牛顿迭代算法的隐式静力学方法和以中心差分算法求解微分方程组的显式动力学方法． 隐式静力学方法的时间积分能真正满足金属轧制的特点，且求解精度高，计算结果无条件稳定，但在 迭代求解非线性刚度矩阵以获得收敛解的过程中，需要耗费大量的计算时间和存储空间_l1D_；而显式动力 学方法总是采用一个对角的或者集 中的质量矩阵求解加速度 ，并用中心差分格式对时间进行离散 ，不必 构造和计算总体刚度矩阵，可大大提高计算效率，但受时间增量稳定性所限，其时间步长往往非常小．针 对所研究的具体问题选择合适的轧制模拟算法，首先要明确不同算法对SL~J模拟结果的影响． 多位学者通过选择显式动力或隐式静力有限元方法I9。。]对金属轧制过程进行计算模拟，而对两种模 拟方法下的计算结果和计算效率的对比研究较少L9 ．有文献[93针对铝板轧制过程，比较了两种方法下 轧辊与轧件接触弧长上的压力分布；有文献E]分别采用两种方法对轧制问题进行二维数值仿真，对计算 效率进行了对比分析．但两种方法对工件轧制后几何尺寸等计算结果的影响有待进一步研究．已有的轧 制模拟研究主要针对铝板和钢板进行，对锆板的相关研究尚未见报道． 本文给出了显式动力和隐式静力两种方法所对应的微分控制方程，指出两种理论模型的主要区别． 针对锆板轧制加工过程建立了三维的有限元模型，分别采用上述两种计算方法，在有限元软件 ABAQus 中实现了对单板SL~J过程的计算模拟．主要比较了两种方法所得到的轧件尺寸变化以及稳定轧制阶段接 触面上的轧制压力分布．此研究可为板型核燃料元件轧制模拟方法的选取提供经验与参考． 收稿 日期 ：2014一O5—28 基金项 目：国家 自然科学基金项 目(11】72068 作者简介：徐宏彬 (1989一一)，男，硕士研究生；丁淑蓉(1971一)，女，副教授 ，通讯联系人，E-mail：dsrl971@163．COnL 复 旦 学 报 (自然科学版) 第54卷 1 控制方程 对于轧制加工这种准静