【2017年整理】板材热轧热力耦合有限元模拟.docVIP

板材热轧热力耦合有限元模拟* 刘　才　崔振山(燕山大学　秦皇岛　066004) 摘要　建立了三维热力耦合问题弹塑性有限变形有限元方法，并对板材热轧过程进行了计算机模拟。材料流动应力模型中考虑了应变历史、应变速率和温度的影响，导出了与其相应的本构关系矩阵。应用该方法可以给出轧件变形过程中诸如温度场、应变场和应变速率场等各种热力结果。算例表明，该方法具有较好的精度。叙词：热力耦合　板材热轧　有限元模拟前言　　板材热轧过程一般是三维大变形热弹塑性问题，变形过程中塑性功转化为热，变形体与工具和周围介质间存在热交换，从而在变形体内产生较大的温度梯度，而温度、应变、应变速率和流动应力相互影响，并进一步影响到变形后金属的微观组织及其机械性质，因而热轧过程是高度非线性的热力耦合问题。研究热变形过程的温度场、应变场、应变速率场的分布与变化，是控制金属微观组织和力学性能的基础。　　基于修正的Lagrange描述(U.L法)的三维大变形弹塑性有限元法已经成功地模拟了冷轧板过程［1，2］，本文在此基础上进一步考虑温度和变形的耦合过程，构造了板材热轧变形过程的耦合联立求解有限元方程。在变形分析中，金属流动应力取为应变历史、应变速率和温度的函数，因而在应力—应变关系矩阵中包含了应变速率和温度的影响。在温度分析中，考虑了塑性变形功和摩擦功转化的热量，并考虑了轧件和轧辊及周围介质间的各种热量传递的边界条件。文后给出的算例证明这种方法在模拟热轧过程时具有很好的精度。 1　变形场有限元模型　　采用U.L法描述，基于有限变形的右伸长张量，可构造线性的共旋应变(称LCR应变) 式中　R——转动张量分量　u——位移分量　　a——参照构形坐标　　在Δt时间内，第二类Piola-Kirchhoff(简称PK2)应力增量为 式中　ΔσJ——Jaumann应力增量　　　　σ——Euler应力　　　　Δe——无限小Euler应变在应变步长不很大时，可近似认为(δim为Kronecker符号)，故ΔSij=δimδTjn(Δσmn)　　在金属热变形过程中，流动应力为等效应变、等效应变速率ε和温度T的函数， 对于服从Mises屈服条件的材料，取其等效应力为塑性位势 基于塑性势流动理论，可得到 式中　Depijkl——流动应力只受应变历史影响时的弹塑性本构矩阵元素 　　　　　α——线膨胀系数dσvij，dσTij——流动应力受应变速率和温度影响时产生的表观应力增量　　对于有限变形情况，将式(7)应力、应变取为Jaumann应力增量和LCR应变增量形式，并代入式(3)和式(2)得到 采用U.L法求解时，将式(8)和式(1)代入增率型虚功原理［3］，并经等体积变形约束［4］修正后，得到有限元平衡方程 (K0+Kσ+KV)Δd=ΔR　　　　　　　　　　　　　　　　　(10) 式中　K0，Kσ，KV——有限变形增量刚度矩阵、初应力刚度矩阵和等体积变形修正矩阵 2　温度场有限元模型　　金属塑性变形温度场是一个具有内热源的不稳定传热问题，其控制方程和边界条件是 式中　ρ,c,k——材料常数　　r——内热源　　　　q——热流速率　　　　　h——换热系数　　　σ，η——Stefan-Boltsmann常数和物体综合黑度　　　Ts，Tamb——工件表面温度和环境温度热量分配系数、接触表面摩擦力和相对运动速度　　在热变形过程中，内热源为塑性变形功转化的热量，而摩擦生热应分配在变形金属和工具上。对于变形金属，取分配系数m=Tt/(Tt+Tb)，其中Tt、Tb分别为工具和工件的温度。　　应用Galerkin法，得到考虑以上边界条件的温度场有限元方程为 式中　K，H，C——传热系数矩阵　　Q—热流密度矢量再对式(16)时间离散，即可求解。 3　求解方法　　交替迭代式(10)和式(16)即得到热力耦合问题的解。在求解式(10)时，采用了割线增量法，应力计算采用了中点正交法则(Mean-Normal Method)［4］和静水压力的间接计算法。对于轧件进入辊缝的过程，采用了摩擦元拖动模型［5］。 板材热轧热力耦合有限元模拟* 刘　才　崔振山(燕山大学　秦皇岛　066004) 摘要　建立了三维热力耦合问题弹塑性有限变形有限元方法，并对板材热轧过程进行了计算机模拟。材料流动应力模型中考虑了应变历史、应变速率和温度的影响，导出了与其相应的本构关系矩阵。应用该方法可以给出轧件变形过程中诸如温度场、应变场和应变速率场等各种热力结果。算例表明，该方法具有较好的精度。叙词：热力耦合　板材热轧　有限元模拟中图分类号：TG335 0　前言　　板材热轧过程一般是三维大变形热弹塑性问题，变形过程中塑性功转化为热，变形体与工具和周围介质间存在热交换，从而在变形体内产生较大的温度梯度