高锰钢晶界微观应力状态的有限元计算.pdfVIP

维普资讯 第 34卷 第6期 报 vn1．34 No．6 1998年 12月 CHINESEJOURNALOFMECHANICAL ENGINEERING Dee． 1998 陈 达 薛侃时／王 位 王 滨 (上海交通大学材料学院 200030) (机械工业都上海材料研究所) ． 孙大涌 ，r 、 (机械工业部机械科学研究院) 、P 彳 I0 ’ 摘要 通过对高锰钢晶粒内部和晶界偏聚硬壳层区域的有限元应力模拟计算．获得 了不同偏聚程 度条件下的应力状况。计算结果表 明．晶界应力集中程度艟晶界Mr,偏聚浓度和偏聚层厚度的上 升而显著提高，由于晶界偏聚而导致的晶界区域硬壳层对高锰钢的断裂性能是非常有害的。从模 拟计算这一角度证实了晶界应力集中是高锰钢产生脆断的主要塘c因，这与宏微观结构成分和力学 性能试验所得到的结果吻合。 翱调； 法 垄盐苎 叵 铜一 鲤 中圈分粪号：trG144 洳 0 前言 高锰奥氏体钢的低温脆断机理一直是国际上争论的热点，众说纷纭，始终未有公断。上海 材料研究所薛佩时曾提出低温脆断的硬壳层设想_1l。近年来，课题组通过大量试验，进一步 证实了高锰奥氏体钢固溶淬火后，在奥氏体晶界上出现锰的不平衡偏聚，而且弄清了锰的偏聚 规律及合金元紊对锰偏聚的影响口]。试验确认了晶界处锰偏聚层的显微硬度比晶内高。并且 在形变过程中应变强化速率也比晶内高，位错塞积严重。应力集中加剧，在外载荷和内应力共 同作用下，导致材料在低温下晶界脆断 ]。 本文的目的就是从理论计算的角度。考察由于Mn元紊在晶界偏聚引起晶界区域产生固 溶强化及其变形过程的形变强化，从而导致晶界区域应力集中的具体状况。并了解不同Mn偏 聚量和偏聚厚度对应力集中状况的影响，以其对比说明试验所观察到的现象和理论验证所提 出的硬壳层模型。 l 有限元计算模型 1．1 计算模型的选择 为简化应力状态及其相关过程的有限元计算。将偏聚态晶粒设置为一个具有一定厚度外 壳层的立方体。如图1所示。图中 为外拉伸方向，z为模拟晶粒长度，阴影部分为晶界壳层 区域，d为晶界壳层厚度(本文分别取 a为0。5％￡和 2％￡)，T、M、B点所示位置分别处 于晶界壳层的底部、中间和顶部。B点实际上就是晶界壳层与内部晶粒的交界位置。将分别在上 * 罾家 自然科学基金和机械工业技术发展基金资助项 目到祈祷到修改j膏 维普资讯 机 械 工 程 学 报 第 34卷 述三点位置计算记录拉伸应变过程中的轴向应力值。外壳层假设为晶界偏聚硬化层 (图1中 的阴影部分为其横截面)，内部则为基体晶粒 在有限元计 算过程中，沿图1所示外拉伸应力OA的方向逐渐加大应变 量，并计算得出不同拉伸应变量情况下，晶粒内部平均轴向 拉伸应力和晶界壳层区域不同位置(底部、中间和顶部处) 的轴向拉伸应力。计算中假设晶粒为各向同性介质。 虽然用立方体来代表晶粒形态并不与实际情况完全相 符，但由于本文读取数据的位置处在拉伸轴上，即立方体的 中心线上(图1中的B、M 、T点)，应力状态不会受到立方 体边、角区域应力集中(奇异)的影响，这三点处的应力状况 较真实地反映了由于晶界壳层和内部晶粒本身强度差异所 图1有限元计算模型示意图 造成的应力集中行为。从计算获得的整个晶粒的应力分布图也说明了仅在立方体边、角局部 范围存在由于晶粒形态所引起的应力集