[发明专利]基于混合硬化模型的微粒流冲击金属材料力学行为预测方法

说明书

本发明公开了一种基于混合硬化模型的微粒流冲击金属材料力学行为预测方法，包括以下步骤：(1)构建微粒位置在空间随机分布的微粒流冲击金属材料的有限元模型；(2)建立综合考虑应变强化、应变率强化、相变强化以及随动硬化效应的材料率相关混合硬化弹塑性本构模型；(3)编写ABAQUS‑VUMAT用户动态材料子程序，运用应力补偿更新算法实现提出的本构模型；(4)将子程序嵌入到ABAQUS中，对微粒流冲击金属构件的过程进行计算，进一步预测金属材料的力学行为。本发明可以用于模拟在空间随机分布的任意数量的微粒、对任意形式的金属结构进行冲击的力学行为预测。

技术领域

本发明涉及金属表面强化处理领域，尤其涉及一种基于混合硬化模型的微粒流冲击金属材料力学行为预测方法。

背景技术

微粒流冲击技术是一项先进的金属表面处理技术，相比喷丸处理技术，其微粒直径更小，而冲击速度却更高。微粒流冲击金属构件的过程，诱导表层组织发生塑性变形，引起晶粒细化、并且引入有益的残余压应力，因此，能够提高金属构件表面强度和硬度、抗疲劳性能以及抗应力腐蚀性能。由于实验研究耗时长、成本高，目前有限元法成为力学行为预测的主要研究手段。

在实际的工程中，大量微粒喷出后是随机撞击到待处理的金属表面上的。之前研究中的数值模型主要有单个微粒模型、多个微粒模型等形式，其冲击位置是固定的，忽略了真实过程中冲击凹坑位置的随机性。近几年，有学者[1]采用有限元计算软件ABAQUS提供的python 语言开发一种弹丸在空间位置随机分布的喷丸模型，可以更加真实地模拟喷丸冲击过程。为金属材料在微粒流冲击下力学行为的有限元预测提高很好的借鉴。

准确描述材料的本构关系是预测金属构件在微粒流冲击下力学行为的前提。在目前的研究中材料的本构关系多采用Johnson-Cook模型[2]，该模型反映了材料的应变硬化和应变率硬化。但是在大量微粒反复冲击金属构件的过程中，金属表面承受非线性弹塑性循环载荷的作用，会产生包申格效应。随动硬化模型能够考虑包申格效应的存在，即认为材料在拉伸方向强化，而在压缩方向同等地软化。而且，微粒流冲击金属构件引起表层组织发生塑性变形，塑性应变诱导金属组织发生相变，因此还需要考虑材料的相变硬化的影响。Sanjurjo[3]等提出的本构方程考虑了应变率硬化和随动硬化效应，但是没有考虑相变硬化。凌祥等[4]的研究中关注了马氏体相变强化的影响，然而屈服强度的提高仅依据经验公式进行推算，缺乏可靠性。综上所述，需要提出一种综合考虑应变硬化、应变率硬化、相变硬化以及随动硬化的材料本构模型。

参考文献：

[1]李源,雷丽萍,曾攀.弹丸束喷丸有限元模型数值模拟及试验研究[J].机械工程学报, 2011,47(22):43-48.

[2]李雁淮,王飞,吕坚,等.单丸粒喷丸模型和多丸粒喷丸模型的有限元模拟[J].西安交通大学学报,2007,41(3):348-352.

[3]Sanjurjo P.,et al.Influence of the target material constitutivemodel on the numerical simulation of a shot peening process[J].Surface andCoatings Technology,2014,258:822-831.

[4]凌祥,彭薇薇,倪红芳.喷丸三维残余应力场的有限元模拟[J].机械工程学报,2006, 42(8):182-189.

发明内容

本发明的主要目的在于克服现有技术中的不足，提供一种基于混合硬化模型的微粒流冲击金属材料力学行为预测方法。

为解决上述技术问题，本发明的解决方案是：一种基于混合硬化模型的微粒流冲击金属材料力学行为预测方法，包括以下步骤：

(1)构建微粒位置在空间随机分布的微粒流冲击金属材料的有限元模型；