采取共旋应变的三维热弹塑性有限变形有限元法ξ.pdfVIP

第 15 卷第 2 期 计 算 力 学 学 报 V o l. 15 N o. 2 1998 年 5 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS M ay 1998 采用共旋应变的三维热弹塑性 有限变形有限元法 崔振山　　刘　才 (燕山大学, 秦皇岛, 066004) 摘　要　本文采用线性化共旋应变张量和增率型虚功原理, 建立了有限变形热力耦合 弹塑性有限元法。在该方法中, 材料的流动应力取为应变总量、应变速率和温度的函数, 推导了包含这种函数关系的本构矩阵。另外在温度场分析中, 考虑了塑性功和摩擦功转 化的热量。文后给出的算例表明该方法可以很好地模拟热加工过程。 关键词　热力耦合; 弹塑性有限元; 有限变形; 数值模拟 分类号　 3436; 242 21 O O 1　引　　言 塑性成型领域中的绝大多数问题都是三维有限变形热弹塑性问题, 变形过程中金属的温 度、应变、应变速率和流动应力相互作用, 并影响到金属的硬化和软化机理, 因而热成型过程是 高度非线性的热力耦合问题, 三维热力耦合模拟是塑性成型研究的重要发展方向。 用于金属塑性成型数值模拟的主要方法是刚塑性有限元法和弹塑性有限元法。刚塑性有 限元法采用了对能量泛函变分而直接得到速度场的方法, 避开了几何非线性问题, 使迭代速度 [ 1 ] 加快, 例如 Kobayash i 等人的工作 。而弹塑性有限元法的特点是能够按照加载路径得到塑性 区的发展情况, 以及应力、应变、应变速率等的分布规律, 而且能处理卸载问题, 求出变形体的 残余应力。因此从获得的计算信息角度要优越于刚塑性有限元法。然而, 如果弹塑性有限元法 是基于无限小应变构造的, 则增量步长要比较小, 而如果是基于 Green 应变构造的, 则在每步 [ 2 ] 中又要经过较多次数的迭代, 这些均在计算时间上造成了困难。 等人 基于右伸长张 P illinger ( ) 量构造了一种线性的与转动无关的应变 共旋应变, 也称L CR 应变 , 较好地利用了前一步的 计算信息, 可以加大每步的步长, 较好地解决了计算时间问题。 本文采用共旋应变和增率型虚功原理, 构造了热力耦合三维弹塑性大变形过程的 有 U L 限元法。在变形分析中, 金属流动应力取为应变历史、应变速率和温度的函数, 推导了包含应变 速率和温度影响的应力应变关系矩阵。在温度分析中, 考虑了塑性变形功和摩擦功转化的热 量, 以及变形体与周围介质间的各种热量传递边界条件。文后给出一个模拟热轧的算例。 2　共旋应变和本构关系 国家教委博士点专项基金和机械部教育司科研基金资助项目 　收稿日期: 修改稿收到日期: 　崔振山: 男, 1963 年生, 硕士, 副教授 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 2 期　　　　　　　　崔振山等: 采用共旋应变的三维热弹塑性有限变形有限元法