材料本构模型参数的获得方法范文

《中国有色金属学报》2015年第十二期

摘要：

为了获得材料在高应变率下的本构模型参数，提出一种基于激光冲击强化实验和有限元模拟相结合的材料本构模型参数获得方法。该方法首先预估材料在高应变率下的本构模型参数，然后进行有限元模拟，以模拟结果与实验结果对比的误差落在一定范围内为优化目标，对预估的本构模型参数进行修正，并最终得到其具体值。基于此方法，得出TiAl合金在高应变率下的Hugoniot弹性极限(HEL)为1.46GPa。

关键词：

钛铝合金；本构模型；高应变率；激光冲击强化；有限元模拟

激光冲击强化是一种新型的表面强化技术，该技术是利用强激光束辐照涂覆在金属表面的吸收层，使其气化并等离子化形成高压冲击波并向材料内部传播，引起材料高应变率(＞1×106s1)动态响应，使材料表层产生塑性变形并形成高残余压应力，从而提高材料表面硬度、抗磨损、抗腐蚀和抗疲劳等力学性能[15]。目前，关于激光冲击强化中冲击波的产生与发展，冲击波作用下材料的高应变率动态响应，材料显微组织变化规律等许多问题的了解仍不是很清楚[6]。因此，许多学者采用实验与有限元模拟相结合的方法来研究激光冲击强化对材料表面性能的影响，但是由于受高应变率的影响，材料的本构模型参数明显不同于准静态情况下材料的本构模型参数(如对于大多数材料，动态屈服强度随应变率的增大而增大)，为有限元模拟带来了困难，因而，确定材料在高应变率下的本构模型参数就显得相当重要。为了研究材料在高应变率下的本构模型，研究者大多使用膨胀环、分离式Hopkinson杆和Taylor实验等方法来获得材料的力学性能数据，然后，通过拟合和外推得到其本构模型参数[79]。HOGGATT等[9]运用爆炸膨胀环技术，对许多工程材料的本构关系进行了测试，赵峰等[10]运用分离式Hopkinson杆对AZ91铸造镁合金在不同应变率下的动态本构关系进行了研究，REN等[11]运用Taylor杆冲击实验对Ti-6Al-4V合金在高应变率加载条件下的动态断裂行为进行了研究。但是，这些实验方法所获得的应变率范围大都在1×102~1×105s1范围内，无法满足激光冲击强化模拟的要求。要想获得大于1×105s1应变率下的材料动态力学数据，就必须采用冲击、炸药爆轰或辐射等方法[8]，这些实验方法不仅非常危险，而且实验装置复杂、成本较高，对材料本构模型参数的确定带来了很大的困难。为了解决这个问题，本文作者在对TiAl合金进行激光冲击强化实验和有限元模拟的基础上，提出了一种获得材料在高应变率下本构模型参数的方法，并利用该方法，得到了TiAl合金在高应变率下的理想弹塑性本构模型参数。

1本构模型参数的获得方法

在本构模型参数获得方法中，首先预估模型参数，然后进行有限元模拟，以模拟结果与实验结果对比的误差落在一定范围内为优化目标，对预估的材料本构模型参数进行修正，最终得到其具体值。在激光冲击强化过程中，高压冲击波会引起材料的高应变率动态响应，因此，材料的本构模型需考虑高应变率的影响，高应变率下材料的本构模型主要有Johnson-cook模型，Zerilli-Armstrong本构模型，理想弹塑性模型(屈服强度为在高应变率下的屈服强度)等。Johnson-cook模型综合考虑了应变、应变率和温度等的影响，可以较好地反映材料的真实本构关系，但是，模型参数复杂，影响因素较多，多通过测量不同应变速率、不同温度条件下的多组应力应变曲线，利用数值拟合的方法确定。DING等[12]使用理想弹塑性模型对35CD450HRC钢进行了模拟研究，模拟结果与实验结果误差在7%左右。花银群等[13]也使用该模型对TC4进行了研究，结果与实验数据对比误差在2%左右，由此可以看出，在利用理想弹塑性模型对激光冲击强化进行模拟时，可以得到较好的结果，并且该模型参数简单，不考虑应变硬化等因素，只需确定材料在高应变率时的Hugoniot弹性极限(HEL)，当材料中应力波压力小于HEL时，材料不发生塑性应变；当应力波压力大于HEL时，将发生永久的塑性变形[6]。

2有限元模型构建

要想获得正确的本构模型参数，建立合理的有限元模型是关键，在激光冲击强化有限元模拟中，主要涉及到冲击波压力的计算与加载、材料本构模型的选择、网格划分与施加边界条件等问题。由于激光冲击强化过程中冲击波压力的加载时间非常短，在冲击波传播过程中，材料内部会发生各种应力波的反射和相互作用，因而，材料对冲击的响应需要持续相当长的时间才能达到稳定，因此，为得到稳定的残余应力场，有限元分析需要采用显式动态分析与隐式静态分析相结合的方法。显式动态分析采用显式动力有限元算法，分析瞬态冲击下产生的冲击波在材料中的传播以及相互作用过程，得到材料动态响应。然后进行隐式静态分析，释放显式动态分析后的材料内部弹性应变，计算材料内部的平衡状态，得到稳定的残余应力场[14]。

2.1冲击波压力波形计算与加载激光诱导的冲击波在空间和时域的分布将对模拟结果产生极其重要的影响。在空间分布上，由于实验中使用的激光脉冲能量在空间平面内近似均匀分布，因此，可认为激光诱导的冲击波压力在整个光斑内均匀分布。而在冲击波的时域分布上，通常使用FABBRO等[15]提出的模型。该模型认为在水约束状态下，激光产生的冲击波压力大致趋于一个六次多项式分布，冲击波压力脉冲的半峰值脉宽大致为激光脉宽的2~3倍，冲击波的峰值压力如式(2)所示。

2.2材料本构模型参数预估由于TiAl合金在高应变率下的HEL未知，因此，需要首先对其预估，在一维应变下，材料的HEL与动态屈服强度的关系为。因此，只需预估出材料的动态屈服强度，即可求出HEL。根据Campbell和Ferguson的研究，对于大多数材料，动态屈服强度随应变率的增大而增大[8]，图3所示为普通低碳钢屈服强度与应变率的关系[17]。从图3可以看出：在不同应变率下，其屈服强度明显不同，且随着应变率的增大，屈服强度也越来越大。TiAl合金在准静态情况下的屈服强度为600MPa[18]，随着应变率的增大，其屈服强度应该越来越大。因此，在本模拟中，可以首先设置TiAl合金在106s1应变率下的动态屈服强度为700MPa，比准静态情况下高出17%。在模拟中使用的其他材料参数为ρ=3800kg/m3，E=165GPa，υ=0.32。

2.3有限元网格模型及边界条件的确定在有限元模拟中冲击应力波会在模型边界处发生反射，反射波汇聚到中心会对结果造成一定影响，因此，为了防止应力波反射所造成的影响，需要使用半无限大三维实体模型，即在模型的侧面和底面施加无反射边界条件，让应力波透射。模拟时，为了与实验结果有更好的对比，所使用的模型尺寸与实验所用的试件完全一致，即直径为15mm、高为5mm的圆柱体，在圆柱体正中心冲击，由于冲击压力与模型均关于中心对称，为提高计算效率，仅建立1/4模型进行分析计算，并在对称面上施加对称边界条件。最终的有限元模型如图4所示。

3TiAl合金HEL的确定

3.1激光冲击强化实验及结果实验时所用的试件照片如图5所示，材料为TiAl合金，名义成分为Ti-45.5Al-2Cr-2Nb-0.15B。在由中科院沈阳自动化研究所搭建的型号为SIA-LSP-1的激光冲击强化设备上进行，所用的高能脉冲激光器型号为Extra20，脉冲能量为9J，激光波长为1064nm，脉宽为20ns，采用直径为3mm圆形脉冲在中心冲击一次。在激光冲击之前，先将试件在HY4050型豪克能应力消除设备(频率为40kHz，输出振幅为50μm)上处理5min，以消除表面残余应力，再用酒精清洗并用氮气吹干，然后，在试件表面粘贴一层100um厚的黑胶带作为吸收层，并用厚度为2mm的流动水膜作为约束层。然后，在试件中心进行单次激光冲击，冲击之后使用Proto-LXRD型X射线应力测试仪测量激光冲击区域表层与深度方向的残余应力，测量方法采用倾斜固定ψ法，其衍射条件列于表1，深度方向残余应力采用0.5mol/L的Na2SO4+H2SO4溶液进行逐层电解腐蚀测量，X射线辐照区域直径2mm，即测量的是激光冲击区域中心2mm内的平均应力[19]。由于模型与冲击压力的对称性，x方向与y方向残余压力基本相等，分布规律也一致，因此，以下讨论中均使用y方向的残余应力进行表示。测量结果如图6所示，数据误差为±20MPa，从图6可以看出，残余压应力最大值出现在次表层，最大值为350MPa，之后，随着深度的增加，残余压应力逐渐降低，在0.65mm左右减小为0。

3.2误差范围的确定为了获得正确的材料本构模型参数，必须给定合理的误差范围，若误差范围太大，则得到的结果不可靠；误差太小，由于压力模型及实验测量结果等因素的误差，可能根本得不到符合要求的结果。在激光冲击强化有限元模拟中，产生误差的来源很多，当高功率的激光脉冲辐照金属表面后，涂覆在金属表面的吸收层气化并等离子化形成等离子体，该等离子体继续吸收激光能量，并受到约束层的约束形成冲击波，作用在材料表面，因此，该过程是一个非常复杂的物理过程，要对其进行真实模拟是非常困难的。在有限元模拟中，将该过程简化为压力脉冲直接作用在材料表面，从而该过程就引入了误差，具体表现在压力峰值大小、空间分布和时域分布等方面。另外，采用的材料本构模型是理想弹塑性模型，该模型本身就是对真实本构模型的一种简化，因此，该模型存在原理上的误差。在有限元求解时，网格划分的大小，所采用的材料阻尼等也会对结果产生影响；最后，本实验中采用逐层剥离法测量深度方向的残余应力，由于在剥离部分材料后，材料的表面状态发生了改变，会引起残余应力的重新分布，因此，也会引入适量的误差。由于各个误差产生的复杂性，对每个误差进行精确估计是非常困难的，且各个误差之间的相互影响规律也很难确定，因此，若通过计算来确定总体误差范围是非常困难的。而从其他学者的有限元模拟结果与实验结果对比分析中可以发现，当采用相似的实验条件与模拟条件时，残余压应力最大值误差与残余压应力影响深度误差总是在一定的范围内变动，具有一定规律，因此，可以这样认为，该范围即是综合了各种误差之后总的误差范围，可以对该误差范围加以分析，从而得到本例中所需确定的误差范围。图7所示为从文献[1213,2023]中得到的不同材料激光冲击强化实验结果与模拟结果对比的误差，其中残余压应力最大值误差(Errorofthemaximumcompressiveresidualstress，EMCRS)，残余压应力影响深度误差(Errorofaffecteddepthofcompressiveresidualstress，EADCRS)的负值表示模拟值比实验值小，正值表示比实验值大。从图7可以看出，残余压应力最大值与实验值之间的误差主要集中在±10%之内，并且正负分布比较均匀，因此，本文作者将残余压应力最大值误差范围设定为±10%。残余压应力深度误差变化范围比较大，从为37.5%~14%，但误差仍集中在20%~14%之间，因此，将残余压应力影响深度误差范围设定为20%~14%。

3.3模拟结果与实验结果的对比在进行有限元模拟时，首先设定TiAl合金的动态屈服强度为700MPa进行模拟。由于实验中X射线应力测试仪的辐照区域直径为2mm，测量的是中心2mm内的平均应力，所以模拟结果也必须以2mm内的平均应力为基准进行对比分析。因此，在模拟冲击的中心区域选取如图8所示直径为2mm的圆，在其上选取均匀分布的9个点作为测量点，对其残余压应力求平均值，作为有限元模拟的结果。最终得到的残余压应力分布如图9(a)所示，其中FEA表示有限元模拟结果，Experiment表示实验结果。模拟得到的残余压应力最大值为383MPa，比实验值350MPa大9.4%；在误差范围内，但残余压应力影响深度达到0.95mm，比实验值大46%，严重超出了误差范围，并且在深度方向上模拟得到的值均比实验值大，只有个别数据能够落在实验所测数据的误差范围之内，说明预估的材料动态屈服强度比实际值小，使材料过早进入屈服状态，塑性变形加大，造成残余压应力变大。改变材料的动态屈服强度为850MPa，再次进行模拟，其结果如图9(b)所示，残余压应力最大值为300MPa，影响深度为0.68mm，与实验值的误差分别为14.2%和4.6%，残余压应力最大值超出误差范围，并且压应力在深度方向上均比实验值小，且大部分数据都落在实验数据的误差范围之外，说明所选的动态屈服强度大于实际值，使塑性变形比实际值要小，造成残余压应力减小。因此，TiAl合金的动态屈服强度应在700~850MPa之间，取为中间值775MPa进行模拟，所得结果如图9(c)所示，残余压应力最大值为326MPa，比实验值小6.8%；影响深度为0.7mm，比实验值大7.7%，均在误差范围内。并且在0.15~0.7mm内，模拟值与实验值吻合的非常好，说明TiAl合金的动态屈服强度为775MPa是可信的，由此得到的HEL为1.46GPa。但实验所得的残余压应力最大值出现在次表层，而模拟所得的最大值出现在表层，残余压应力最大值出现的位置不相符，造成这种现象的原因主要有两方面：一方面，由于有限元模拟仅仅考虑了激光冲击强化的力学效应，而忽略了其热效应。在实验中，吸收层在激光辐照下会迅速气化并等离子化，在该过程中，等离子化的气体温度高达数千摄氏度，若吸收层太薄，则等离子化的气体会对试件表面产生热影响，使材料表面温度快速升高，当冲击完毕后，冷却水会快速带走大部分热量，使材料快速冷却收缩并在材料表层形成拉应力，这将抵消掉材料表面部分残余压应力，使表层残余压应力减小；另一方面，当冲击波压力为2倍的HEL时，材料表层塑性变形将达到饱和，表层残余压应力达到最大值，继续增大冲击波压力，表层残余压应力将会减小，残余压应力最大值将会向材料内部移动，造成材料表层残余压应力的减小[6]。由于模拟中使用的冲击波峰值压力是2.65GPa，小于2倍的HEL，因此，可能是热效应造成了实验中表层残余压应力的减小，但具体原因仍需后续深入研究。

3.4本构模型参数获得的方法在3.3节中得到TiAl合金的HEL为1.46GPa，由于目前尚无TiAl合金在高应变率下的材料本构模型参数，因此，无法对该结果进行直接验证，但从图9(c)可以看出，模拟结果与实验结果在趋势与数值上都是非常吻合的，证明该结果是可信的，该参数获得方法是有效的。由于该方法对材料没有特殊要求，只要能获得激光冲击强化实验数据，且选定合理的误差范围，即可对任何具有单一参数的本构模型进行很好的确定。但是对于材料在高应变率下的其他本构模型，如Johnson-cook模型，本构关系如式(4)所示。

4结论

1)提出了一种基于激光冲击强化实验和有限元模拟相结合的材料本构模型参数获得方法，该方法能对材料在高应变率下具有单一参数的本构模型进行很好的确定，对需要确定多参数的本构模型有一定局限性。2)建立激光冲击强化有限元模拟的分析模型，并提出材料本构模型参数的预估方法。3)提出本构模型构建方法中误差范围的确定办法，得出了TiAl合金在高应变率下的HEL为1.46GPa。

参考文献

[1]PEYREP,FABBROR,MERRIENP,LIEURADEHP.Lasershockprocessingofaluminiumalloys.Applicationtohighcyclefatiguebehaviour[J].MaterialsScience&Engineering,1996,210(1):102113.

[2]乔红超,高宇,赵吉宾,陆莹,赵亦翔.激光冲击强化技术的研究进展[J].中国有色金属学报,2015,25(7):17441755.QIAOHong-chao,GAOYu,ZHAOJi-bin,LUYing,ZHAOYi-xiang.Researchprocessoflaserpeeningtechnology[J].TheChineseJournalofNonferrousMetals,2015,25(7):17441755.

[3]HUYong-xiang,LIKang-mei,QIChen-jie,YAOZQ,GRANDHIRV.Sizeeffectonindentationdepthofoxygen-freehighpuritycopperinducedbylasershockprocessing[J].TransactionsofNonferrousMetalsSocietyofChina,2012,22(S2):s573s578.

[4]张青来,鲍士喜,王荣,钱阳,张永康,李兴成.激光冲击强化对AZ31和AZ91镁合金表面形貌和电化学腐蚀性能的影响[J].中国有色金属学报,2014,24(10):24652473.ZHANGQing-lai,BAOShi-xi,WANGRong,QIANYang,ZHANGYong-kang,LIXing-cheng.EffectoflasershockprocessingonsurfacemorphologyandelectrochemicalcorrosionresistanceofAZ31andAZ91alloys[J].TheChineseJournalofNonferrousMetals,2014,24(10):24652473.

[5]ZHANGXing-quan,CHENLiu-san,YUXiao-liu,ZUOLi-sheng,ZHOUYu.Effectoflasershockprocessingonfatiguelifeoffastenerhole[J].TransactionsofNonferrousMetalsSocietyofChina,2014,24(4):969974.

[6]李应红.激光冲击强化理论与技术[M].北京:科学出版社,2013:20,133.LIYing-hong.Lasershocktheoryandtechnology[M].Beijing:SciencePress,2013:20,133.

[7]FABBROR,PEYREP,BERTHEL.Physicsandapplicationsoflaser-shockprocessing[J].JournalofLaserApplications,1998,10(6):265.

[8]MEYERSMA.材料的动力学行为[M].张庆明,刘彦,黄风雷,吕中杰,译.北京:国防工业出版社,2006:225.MEYERSMA.Dynamicbehaviorofmaterials[M].ZHANGQing-ming,LIUYan,HUANGFeng-lei,LÜZhong-jie,transl.Beijing:NationalDefenceIndustryPress,2006:225.

[9]HOGGATTCR,RECHTRF.Stress-straindataobtainedathighratesusinganexpandingring[J].ExperimentalMechanics,1969,9(10):441448.

[10]赵峰,李玉龙,索涛,黄卫东,刘建睿.高应变率下铸造镁合金AZ91的动态压缩性能及破坏机理[J].中国有色金属学报,2009,19(7):11631168.ZHAOFeng,LIYu-long,SUOTao,HUANGWei-dong,LIUJian-rui.DynamiccompressivebehavioranddamagemechanismofcastmagnesiumalloyAZ91[J].TheChineseJournalofNonferrousMetals,2009,19(7):11631168.

[11]RENYu,TANCheng-wen,ZHANGJing,WANGFu-chi.DynamicfractureofTi-6Al-4ValloyinTaylorimpacttest[J].TransactionsofNonferrousMetalsSocietyofChina,2011,21(2):223235.

[12]DINGK,YEL.Lasershockpeeningperformanceandprocesssimulation[M].NewYork:CrcPress,2006:89,109.

[13]花银群,蔡峥嵘,陈瑞芳,姜辉,吉光.TC4钛合金激光搭接冲击强化的实验和数值模拟[J].激光技术,2010,34(5):632635,639.HUAYin-qun,CANZheng-rong,CHENRui-fang,JIANGHui,JIGuang.ExperimentandnumericalsimulationofoverlappinglasershockprocessinginTC4titaniumalloy[J].LaserTechnology,2010,34(5):632635,639.

[14]胡永祥.激光冲击处理工艺过程数值建模与冲击效应研究[D].上海:上海交通大学,2008.HUYong-xiang.Researchonthenumericalsimulationandimpacteffectsoflasershockprocessing[D].Shanghai:ShanghaiJiaoTongUniversity,2008.

[15]FABBROR,FOURNIERJ,BALLARDP,DEVAUXD,VIRMONTJ.Physicalstudyoflaser-producedplasmainconfinedgeometry[J].JournalofAppliedPhysics,1990,68(2):775.

[16]戴峰泽.基于纳秒激光冲击波效应的金属表面形貌与性能研究[D].镇江:江苏大学,2014.DAIFeng-ze.Mechanismresearchofnano-secondlaserinducedshockwaveonthesurfacetopographyandperformanceofmetals[D].Zhenjiang:JiangsuUniversity,2014.

[17]柳永宁,朱金华,周惠久.普通低碳钢的强度与温度、应变速率的关系[J].材料科学进展,1990,4(4):285290.LIUYong-ning,ZHUJin-hua,ZHOUHui-jiu.Strengthofalowcarbonsteelatdifferenttemperaturesandstrainrates[J].MaterialsScienceProcess,1990,4(4):285290.

[18]李兴华,杨绍利.钛铝合金制备技术现状及新进展[J].材料导报,2011,25(7):94100.LIXing-hua,YANGShao-li.TheprogressofpreparationandprocessingtechnologyofTiAlalloy[J].MaterialsReview,2011,25(7):94100.

[19]乔红超,赵亦翔,赵吉宾,陆莹.激光冲击强化对TiAl合金组织和性能的影响[J].光学精密工程,2014,22(7):17661773.QIAOHong-chao,ZHAOYi-xiang,ZHAOJi-bin,LUYing.EffectoflaserpeeningonmicrostructuresandpropertiesofTiAlalloy[J].OpticsandPrecisionEngineering,2014,22(7):17661773.

[20]VOOTHALURUR,LIUCR,CHENGGJ.Finiteelementanalysisofthevariationinresidualstressdistributioninlasershockpeeningofsteels[J].JournalofManufacturingScienceandEngineering,2012,134(6):30203023.

[21]胡永祥,姚振强,胡俊.激光冲击强化残余应力场的数值仿真分析[J].中国激光,2006,33(6):846851.HUYong-xiang,YAOZhen-qiang,HUJun.Numericalsimulationofresidualstressfieldforlasershockprocessing[J].ChineseJournalofLasers,2006,33(6):846851.

[22]郭乃国.激光冲击处理40Cr钢及其残余应力场数值模拟[D].镇江:江苏大学,2007.GUONai-guo.Lasershockprocessingof40CrsteelandnumericalsimulationofresidualstressfieldinducedbyLSP[D].Zhenjiang:JiangsuUniversity,2007.

[23]武敬伟.激光冲击Fe-Ni合金残余应力场的数值模拟[D].镇江:江苏大学,2007.WUJing-wei.NumericalsimulationofresidualstressfieldinducedbylasershockprocessingonFe-Nialloy[D].Zhenjiang:JiangsuUniversity,2007.

作者：陈松林 赵吉宾 乔红超 杨灏 单位：中国科学院 沈阳自动化研究所 中国科学院大学