塑性变形与传热耦合分析技术的研究.pdfVIP

维普资讯 、 第 l8卷 第3期 太 原 重 型 机 械 学 院 学 报 Vo!．18 NQ3 1997年 9月 JOURNALOFTAIYUANHEAVYMACHINERYINSTITUTE Sep．1997 塑性变形与传热耦合分析技术的研究 琴 (太太原重型机械学院压力加工 太原030024) ＼ ＼ 形过程中变形与传稿的耦夸分析。 r 量中圈分类剑号T譬G案302翌进竖；嫂 苎竺 投一j 塑性加工中有很多热成形或温成形工艺，即工件的塑性变形是在高温或较高温度下进 行的。此时，温度对金属材料的变形行为影响很大。首先，变形前工件内部的温度分布就可 能不均匀；变形过程中，由于工件表面与环境和模具的热损失以及工件内部变形不均匀产生 的不均匀内热源，使工作内温度梯度较大。这样会造成不同质点材料的屈服应力差别很大， 最终对整个变形过程产生较大影响作用。此外，高温下的塑性变形还会使金属材料产生微 观组织变化，如相变、动态再结晶等等。因此，对于这类工艺问题，必须在进行变形分析的同 时分析温度场的变化，并且考虑二者之间的相互影响作用。 热变形过程中变形与传热的耦合分析早为研究者们所关注，它随着有限元法的发展不 断完善。对于非稳态问题，目前常用耦合迭代法。该方法是将变形传热看作瞬态独立方程， 在两者问迭代求解。耦合迭代法最早由Zienkiewicz[11提出，Rebelo和Kobayashi[2]使之完善 并用于园柱和园环镦粗过程的分析。变形与传热耦合分析的另一种方法是增量区间准静态 耦合技术，作者 】采用基于后差法的技术模拟分析了大型转子JTS工艺试验，著名商品软 件DEFORM[中也采用了该项技术。 塑性成形热传导理论与公式 1．1 有限元公式 金属塑性成形过程的传热属于含内热源的瞬间传热问题，其控制微分方程即能量平衡 收筏 日期：1997—03一毗 本文第一作者；男，1958年5月生．压力加工系翻教授。研究方向为塑性加工模拭与控镧 维普资讯 第 18卷第3期 刘建生等：塑性变形与传热耦合分析技术的研究 227 方程为： (kT．f)．f+ —pd=0 (1) 式中 T和t表示温度和温度率，k是材料热传导系数数，p，c分别为材料密度和比热容， 是内热源率，表示为： 毒=,gi (2) 此时，a是塑性功转变为热能的比例系数，通常a=0．85～0．95，取值与材料有关。 对式(1)采用Galekin法变型，并经有限元离散化后，可得 K T+Ct=Q (3) 式中 K—— 热传导矩阵； C—— 热容矩阵； Q—— 结点热流向量； T，t—— 结点温度和温度率向量 热流向量Q由以下各项组成，表示为 Q：l Hd +l ( 一 )HdS+1^( — )Hds+ ～ J V ～ J ～ J — Jr r s^船(Td—T,)HdSJsqjt-一tdS (4) 其中H是形函数向量， 是工件表面温度， 和 分别是环境和接触面模具温度。式(4) 右端第一项是塑性变形产生的热流；第二、第三项分别是工件 自由表面与环境间辐射和对流 换热的热流；第四项为工件与模具接触面热传导的热流；最后一项是接触表面摩擦产生的热 流，即ql=If』·』“I． 1．2 时间差分格式