中国西北工业大学&美国北卡罗来州大学MSEA:基于晶体塑性有限元法的FCC金属绝热剪切局部数值模拟
【前言】
剪切带是塑性应变面心立方(FCC)金属中最常见和公认的微结构特征之一。在重轧制金属板和普通应变压缩材料中,它以集中塑性流动带状区域的形式出现。在这种情况下,剪切带作为一种替代变形模式来适应塑性应变,这是导致应力集中局部化的原因。剪切带现象在近年来受到了研究人员极大关注,这不仅是因为它们代表了许多金属在塑性应变时微观结构的特征,还因为它们对带内微观结构和再结晶织构的演变产生了重大影响。在具有不同堆垛层错能(SFE)的FCC金属中,剪切带可分别为“铜型”和“黄铜型”剪切带。
研究人员目前已经对FCC金属中的剪切局部进行了一些实验和数值研究。然而,许多细节仍不清楚,需要进一步阐述。绝热剪切带(ASBs)与上述离散剪切带有很大不同,因为前者塑性变形过程中的热效应不可忽略。由于热软化导致的材料强度损失超过了由于应变硬化或应变速率硬化导致的强度增加。最终,ASBs会发展成为热塑性不稳定的结果。
事实上,在FCC金属中,ASBs很少在动态单轴压缩下被发现。研究人员通过设计具有特殊几何形状的样品(帽形或截头圆锥体样本)来进行观测,当样本处于机械载荷下时,某些局部区域会引入集中的大剪切应变。但是这种特殊的几何形状并没有抓住ASB材料的“内在”表现。研究人员在ASBs内部对这些特殊设计的样品进行微观观察,以研究局部变形过程中的微观结构演变。近年来,研究人员在BCC金属中发现晶粒尺寸的减小能够有效地促进在ASBs的形成。然而,对于FCC金属,直到最近,即使当材料的粒度细化到纳米级(小于100 nm )研究人员仍旧没有在高速压缩下观察到ASBs。因此,在FCC金属中,形成ASBs的实验证据的缺乏阻碍了对这些金属中ASBs的深入研究。直到最近,在298 K和77 K的动态单轴压缩后,在具有FCC晶体结构的纳米结构(NS) Cu-Al合金中研究人员发现了共轭的ASBs。高应变速率下的各向异性行为和纹理测量表明,严重塑性变形(SPD)期间形成的预先存在的纹理应该是试样在高速压缩下加载时剪切定位的原因。晶粒细化导致的强度增强可以提供足够的变形能量,但是塑性加工产生的热量会导致最终热塑性的不稳定。
【成果简介】
近日,来自中国西北工业大学的索涛教授和美国北卡罗来州大学的Qiuming Wei教授(共同通讯)联合在Materials Science and Engineering: A上发表文章,题为:“Numerical simulations of adiabatic shear localization in textured FCC metal based on crystal plasticity finite element method”。在这项工作中,作者提出了具有不同初始纹理的多晶模型中动态单轴压缩和简单剪切的晶体塑性有限元模拟( CPFEM )。目的是研究在高应变率载荷下,FCC金属在严重塑性变形或再结晶后经常观察到的典型织构对ASBs形成的影响。材料的响应由弹粘塑性连续滑移本构关系描述,其中也考虑了滑移系统的电阻对温度演化的依赖。模拟结果表明,在高速压缩下,不同的纹理导致剪切带的不同取向和剪切定位发生时的不同临界应变。研究人员发现高速简单剪切载荷有助于ASB的形成。即使在无纹理,即随机纹理中,也观察到模型明显的剪切局部化。然而,滑移系统激活对纹理的依赖导致ASB引发的临界剪切应变水平变化巨大。此外,作者还讨论了一些因素,如应变率、温升和材料强度对FCC织构化金属中ASBs形成的影响来更深入地解释形成过程。作者还将实验观察和数值模拟相结合,阐述了催化裂化材料中很少报道ASBs形成的原因。
【图文导读】
图1. FCC 格子的示意图与一些 {111} 滑动面实体图
FCC 格子的示意图与一些 {111} 滑动面实体图。
图2. 关于轴d在坐标系的旋转从XYZ到X’Y’Z’的特点是三旋转角δ, ψ和ω
关于轴d在坐标系的旋转从XYZ到X’Y’Z’的特点是三旋转角δ, ψ和ω。
图3. 加载边界条件
(a)单轴压缩,(b)简单剪切,(c)加载历史,以及(d)用于CPFEM模拟的样本的有限元表示的加载边界条件。
图4. 本文研究的FCC金属中的典型纹理
除了(a)随机纹理,(001)极图还显示了其他四种类型的纹理: (b)立方体{100} <001>、(c)高斯{110} <100>、(d)黄铜{110} <112>和( e )铜{112} <111>;
(f) 晶体单元中滑移系统的示意图,优选取向为高斯、黄铜和铜纹理;
图5. FCC单晶模型的初步结果显示,当所有具有相同取向的元素(Goss (110) [001])在三个正交方向上加载时,在动态加载(104 s−1)下,工程应变高达50 %,变形模式不同
(a)在X方向上均匀变形;
(b)和(c) Y或Z方向的局部剪切;
(d)单轴动态压缩下单晶和多晶模型的真实应力应变响应;
图6. FCC金属中五种典型织构的应力-应变曲线
FCC金属中五种典型织构在(a)单轴动态压缩(应变率~104 s−1)和(b)动态简单剪切(应变率~ 104 s−1)下的应力-应变曲线。
图7. 五种不同类型纹理的局部形变
(a1)-(a6) 五种不同类型的纹理在动态单轴压缩下的变形模式(ε̇=104 s−1) : 1。随机2. 立方体;3. 高斯;4. 铜;5和6。黄铜纹理, 为了清楚地显示变形行为, 在图中也给出了黄铜纹理模型的位移分布a5;
(b) 实验观察了具有强高斯纹理的纳米结构 FCC 合金的局部形变。
图8. 五种纹理动态加载的绝热剪切
(a1)-(a5) 在简单剪切下的五种纹理动态加载的绝热剪切定位(γ̇=104s−1) : 1。随机2. 立方体;3. 高斯;4. 黄铜;5. 铜质地。(b) 实验结果表明无氧铜铜的绝热剪切局部化的形成。
图9. 不同载荷下的多晶样品
在(a)准静态载荷(0.1 s−1)和(b)动态载荷(104 s-1)下压缩具有高达50 %工程应变的Goss织构的多晶样品。
图10. 模型的应变与应变率
(a) 具有铜织构的多晶模型在5×102 s−1至2×104 s−1的不同应变速率下的简单剪切变形;
(b) 局部应变与剪切应变率。
图11. 局部区域元素记录的机械响应(真实应力-应变曲线)和绝热升温
(a)在动态单轴压缩下;(b)在动态简单剪切下。
图12. 不同温度对模型的影响
(a)在高应变率(104 s-1)下,不同加载温度下具有Goss纹理的压缩模型的机械响应。插图显示了变形过程中硬化率的变化;
(b)不同加载温度下动态加载后(工程应变为0.5 )的升温分布: ( b1 ) 298 K,( b2 ) 573 K
图13. 集中变形研究
(a)具有铜纹理的压缩模型的机械响应,其中针对滑移系统给出不同的临界分解剪应力(CRSS)值,以模拟不同强度的不同材料;
(b) 显示了动态加载后(工程应变为0.5时)的升温分布,以了解强度或塑性应变能对变形集中的影响。
【总结】
作者建立了一个动态热-机械耦合晶体塑性模型来研究织构化FCC多晶材料中ASB的形成。本文考虑了五种典型的纹理类型,并采用旋转对称的错位高斯分布来生成更实用的微观纹理。为了与实验结果进行比较,作者模拟了高应变率压缩和简单剪切载荷条件,模拟结果与文献中的实验结果吻合良好。
文献链接:Numerical simulations of adiabatic shear localization in textured FCC metal based on crystal plasticity finite element method, (Materials Science and Engineering: A, 2018, DOI: doi.org/10.1016/j.msea.2018.08.105).
本文由材料人编辑部金属组Z. Chen供稿,材料牛编辑整理。
欢迎大家到材料人宣传科技成果并对文献进行深入解读,投稿邮箱tougao@cailiaoren.com。
投稿以及内容合作可加编辑微信:cailiaokefu。
文章评论(0)