超塑变形中硬化规律在多尺度数值模拟中的应用_韩晓宁.pdf

超塑变形中硬化规律在多尺度数值模拟中的应用韩晓宁1，杨俊宙2，吴建军2，李金山1(1.西北工业大学 材料学院，陕西 西安 710072；2.西北工业大学 机电学院，陕西 西安 710072)摘要：为了提高TC4超塑成形多尺度仿真数值模拟的有效性，本文针对高温下材料的硬化规律进行了研究。基于材料流变行为的特点，提出了应变硬化和应变软化的表征模型，并将其应用于超塑成形多尺度数值模拟过程中。研究结果表明，本文所采用的应变补偿模型，可用来描述TC4钛合金在超塑变形中的流变行为。此外，对比晶粒尺度、空洞体积分数等微观尺度变量与多尺度物理本构模型的结果可知，所提出的方案对热变形中微观组织演变的数值模拟是有效的。关键词：超塑变形；数值模拟；应变补偿；本构模型；钛合金中图分类号：TG146.2+3文献标识码：A文章编号：1000-8365(2023)01-0075-05Application of the Hardening Law in the Multiscale Simulation forSuperplastic DeformationHAN Xiaoning1,YANG Junzhou2,WU Jianjun2,LI Jinshan1(1.School of Materials Science and Engineering,Northwestern Polytechnical University,Xian 710072,China;2.School ofMechanical Engineering,Northwestern Polytechnical University,Xian 710072,China)Abstract：To improve the precision of multiscale numerical simulations for the superplastic formation of TC4 alloys,thehardening law of metal materials at high temperatures was investigated.Based on the flow behavior of high-temperaturetensile tests,a simple model considering strain hardening and softening was proposed.Then,the established model wasapplied to the multiscale numerical simulation of superplastic forming.The results show that the strain compensation modelcan be employed to describe the flow behavior of TC4 titanium alloy during superplastic deformation.In addition,comparing the simulated results of microscale variables such as grain size and cavity volume fraction with the calculatedresults of a multiscale physically based constitutive model,it can be concluded that the proposed method could be effectivefor the numerical simulation of flow stresses and microstructure evolution during hot deformation.Key words：superplastic deformation;numerical simulation;strain compensation;constitutive model;titanium alloy收稿日期:2022-12-02基金项目:国家重点研发计划(2021YFB3702603)；国家自然科学基金(52174377)作者简介:韩晓宁,1981年生，博士，研究员.研究方向：钛合金超塑成形技术研究.E-mail：引用格式:韩晓宁，杨俊宙，吴建军，等.超塑变形中硬化规律在多尺度数值模拟中的应用J.铸造技术，2023,44(1):75-79.HANXN,YANGJ Z,WUJ J,et al.Application of the hardening law in the multiscale simulation for superplastic deformationJ.FoundryTechnology,2023,44(1):75-79.TC4钛合金由于其良好的机械性能，已经是航空、医疗、汽车等领域中的重要材料1-3。超塑性是指材料在特定温度、应变率区间内可以达到很大伸长量而不发生断裂行为的性能4。超塑成形技术常用于生产具有复杂几何结构的零件，常见于三层板、四层板和点阵结构中。有限元数值模拟是超塑成形加载路径设计中的一种有效手段5，根据虚拟试验结果寻找不同温度、应变率下的变形规律，可以低成本地得到超塑成形的工艺窗口。因此，应用可靠的有限元手段来描述超塑变形过程是成形工艺研究的基础，同时显著提高数值模拟的有效性对描述该过程也非常重要，这也是本文关注的重点。目前，针对超塑成形过程中力学行为和变形规律的数值模拟已被广泛研究。其中，大量文献聚焦于材料成形性与成形工艺之间影响规律的研究，尤其是针对复杂结构中成形极限的判断方法6。然而，对于超塑变形中其复杂的变形机理所导致的结构件成形窗口判断的复杂问题却少有研究。在实际生产中，目前较为常用的数值模拟方法主要有2种：稳定性原则方法和成形极限图(forming limit diagrams,FLD)。Hart的稳定性准则是一种常见的准则，而Ghosh经过研究后发现Hart并不具有普适性7。此外，由于不DOI：10.16410/j.issn1000-8365.2023.2344铸造技术FOUNDRY TECHNOLOGYVol.44 No.01Jan.2023试验研究Experimental Research75同的稳定性准则通常是在特定的本构模型的基础上建立的，导致超塑变形模型应变率影响的复杂性又为该方法带来了新的挑战。因此，在研究塑性变形成形极限时，FLD曲线受到了更多的青睐，该曲线是基于Hill失稳准则和Swift失稳准则而建立的8，综合考虑了超塑变形应变率的影响。Krhn等9将传统的FLD技术扩展到超塑成形工艺中，并利用超塑成 形 极 限 图(superplastic forming limit diagrams,SPFLD)研究了胀形超塑性过程的塑性稳定性。此外，Alabort等10针对三层空心结构采用SPFLD进行成形极限分析，进一步验证了该方法的有效性。同时，在针对三层结构成形过程的研究中，还充分利用多尺度数值模拟技术，以考虑成形过程中晶粒尺寸与空洞体积分数的演变，即在仿真成形过程中不仅考虑了宏观力学性能与结构之间的关系，还考虑了其微观组织的演变规律。因此，这种方案能更加全面地探究结构件成形工艺，本课题组也对这一方法进行了探索11。在数值模拟过程中，材料的硬化规律采用CREEP蠕变子程序进行表征5，虽然该模型已经被广泛应用于超塑成形数值模拟过程12，但其最大的弊端就是只能捕捉材料稳定变化时的流变行为。Du等13在对该结构进行超塑成形处理时，发现材料在变形过程中其晶粒发生了明显的生长现 象，变形后材料的平均晶粒尺寸从5 m增加到10 m，这本质上是应变影响的结果。Alabort等10在处理这一问题时将参数进行了拟合，这种方法思路比较简单，但是存在求解过程中计算量过大等问题，且容易引起应力表征出现硬化规律表征时的波动情况。因此，有必要建立一种针对多尺度数值模拟中硬化规律的表征方法，进而提高仿真的有效性。基于以上分析和讨论，本文根据材料的流变行为对其硬化规律进行表征，随后将其应用于多尺度数值模拟中，对材料的超塑成形过程进行仿真，进而验证该方法的有效性。1超塑变形力学性能表征1.1高温拉伸实验本研究采用的实验材料为1.0 mm厚的TC4钛 合金板材。高温拉伸试验在岛津电子万能试验机(AG-Xplns100KN)上进行，该 设 备 可 实 现 室 温1100 区间内的拉伸试验，最大拉应力为100kN，拉伸速率为0.001800.000 mm/min。加热炉有3个加热区，分别由3个热电偶测量温度，并由闭环温控系统实现控温，保证试件在试验过程中保持所需的温度，以满足实验设计的要求。本试验温度为920，初始应变率采用10-4/s、10-3/s和10-2/s 3种应变速率，该变形条件下TC4材料具有良好的超塑性能14。试样以20/s升温到设定温度，随后保温15 min，以确保材料组织和试样温度的均匀性，然后按设定温度和应变率进行高温拉伸试验。拉伸后的流变应力曲线如图1所示。由图1可知，TC4合金在热变形初期，当初始应变率为10-2/s时，合金的流变应力随其变形程度的增大而快速上升，在达到其峰值后迅速下降，而后下降趋势减缓，直至实验结束。当初始应变率为10-4/s和10-3/s时，合金流变应力随着变形程度的增加而显著上升，在达到最大值后逐渐趋于稳定直至实验结束，表明TC4合金有一定的应变强化效应。此外，图1所示为在相同温度下，当应变速率从10-4/s升至10-2/s过程中，真应变-真应力曲线呈上升趋势，表明TC4合金有一定的应变率强化效应，可知TC4合金为明显应变速率敏感材料。1.2蠕变子程序硬化关系建模对于超塑变形过程而言，材料的应力应变关系通常采用下式进行描述：=k?mp(1)式中，为流动应力；?p为应变速率；k为材料系数；m为应变速率敏感性指数，一般在0.30.9之间。有限元数值模拟是基于CREEP蠕变子程序进行的，因此首先要确定模型中参数k和m的值。对式(1)两边求对数可得式(2)，如下所示：ln=lnk+mln?p(2)式中，为真应力；?p为真应变率。m和lnk分别为ln-ln?p线性拟合的斜率和截距，进一步计算可得不同应变下m和k的值，如图2中散点所示。由图2可知，在不考虑应变对材料力学性能影响的情况下，采取二者的平均值进行计算，其结果分别为m=0.226，k=112.339。根据文献15中的研究结果，对图1 TC4合金高温拉伸试验流变曲线Fig.1 The flow behavior of the TC4 alloy by high-temperaturetensile testsVol.44 No.01Jan.2023FOUNDRY TECHNOLOGY76型涉及的本构模型采用用户定义材料，用户材料选取蠕变子程序。由于需要模型结果中应力试验值作为状态变量进行输出，因此，材料属性中需设置非独立变量为7，分别对应的是应变率、应变、应力、晶粒尺寸、位错密度、再结晶体积分数和空洞体积分数。本研究中，由于标距段内变形是均匀的，因此采用单元3038的计算结果为代表，计算每个增量步中流变应力和微观内变量的结果。其中，内变量的计算结果采用之前研究的结果17。在有限元数值模拟时，对成形应变速率进行控制时，其本质是依赖于对目标应变率进行控制加载。在一个增量步中，新的载荷Pnew是根据控制变形区内所有积分点或单位集中等效应图4有限元软件中的高温拉伸仿真模型16Fig.4 Simulation model of the high temperature tensile test16图2模型中参数拟合结果：(a)k值，(b)m值Fig.2 The parameter fitting results of the constitutive model:(a