[发明专利]一种基于代表体积单元的Ni3Al基合金本构模型建立方法

说明书

技术领域

本发明属于航空材料技术领域，具体涉及一种基于代表体积单元的Ni₃Al基合金本构模型建立方法。

背景技术

在研究Ni₃Al基合金内部γ′/γ结构时，研究者通常都是以单胞的形式进行研究，但实际上材料内部强化相γ′相的分布并不是完全有序的，而采用单胞进行分析就忽略了强化相分布的不均匀性对材料性能带来的影响。图4给出了单胞与RVE以及微观结构的关系，其中图4(a)为拉伸试验的的试件，图4(b)为通过扫描电子显微镜(SEM)拍摄的IC10合金微观结构图，在微观结构中选择一个包含足够多强化相的区域，这个包含了足够多微观信息的区域就是代表体积单元(RVE)(如图4(c)所示)，代表性体积单元中包括大量的单胞(unit cell)，单胞与强化相的尺度相同(如图4(d)所示)。从图4中可以看出，强化相的大小和形状都不是完全相同的，而且强化相也不是完全有序的排列在基体相中，其内部结构充满了不确定性，虽然单胞在一定程度上体现了材料微观结构的特性，但它更多描述的是材料的理想平均微结构特征，与材料的实际微结构相去甚远，不能包含材料的微观不确定性，难以有效描述材料内部的力学性能。代表体积单元由大量的强化相和基体相组成，其尺度比单胞的尺度更高，包含了更多的微观信息，为考虑材料微观信息的不确定性，可以基于强化相的大小和形状特征参量的统计分布规律来确定各强化相的大小和形状分布，对代表性体积单元(RVE)进行分析，可以更有效的描述材料的力学性能。

发明内容

发明目的：为了克服现有技术中存在的不足，本发明提供一种基于代表体积单元(RVE)的Ni₃Al基合金本构模型建立方法，能够提高了材料在受载后力学响应的预测精度。

技术方案：为实现上述目的，本发明采用的技术方案为：

一种基于代表体积单元的Ni₃Al基合金本构模型建立方法，包括如下步骤：

1)根据材料的微观结构构造生成材料的微细观代表体积单元RVE，在RVE中强化相和基体相均采用粘塑性模型；

2)基于晶体塑性理论，建立材料的本构模型并通过试验获取模型的参数；

3)建立基于RVE的Ni₃Al基合金的本构模型；

4)模型验证：模型建立后，对Ni₃Al基合金IC10在600℃下的拉伸性能进行模拟验证。

进一步的，步骤2)中，所述参数包括：

2-1)对Ni₃Al单晶进行拉伸试验得到强化相的弹性常数C₁₁、C₁₂和C₄₄；

2-2)对纯Ni进行纳米压痕试验得到基体相的弹性常数C₁₁、C₁₂和C₄₄；

2-3)对Ni₃Al基合金IC10进行不同应变载荷下的拉伸试验，对拉伸后的试件进行切片，通过透射电子显微镜观测材料内部的位错分布；并用图像处理软件统计不同应变载荷下材料的位错密度，通过最小二乘法对不同应变下的位错进行拟合，得到ρ_s参考饱和位错密度，ρ₀初始位错密度以及拟合参数θ的值。

进一步的，步骤3)建立基于RVE的Ni₃Al基合金的本构模型的具体方法为：根据RVE中各个强化相的几何特征为强化相和基体通道选定材料，选定材料后先计算各个组分的位错密度，根据位错密度求出滑移系分切应力门槛值，然后带入相应的流动准则计算滑移系剪切率和求出应力张量，最后求出RVE的等效应力。

进一步的，步骤3)建立基于RVE的Ni₃Al基合金的本构模型的具体方法为：

3-1)根据晶体塑性理论得到材料应力率和应变率以及各滑移系剪切率之间的关系：

其中为中间构形的Jaumann应力率，EMT为瞬时弹性模量张量，D、P和B为晶体塑性理论的中间量，为滑移系剪切率；

3-2)位错密度的演化过程为：

其中上式右边的两项分别为随机位错和几何位错，ρ_s为参考饱和位错密度，ρ₀为初始位错密度，b为伯格斯矢量，d表示强化相的大小，h^(αβ)为位错关联因子，k_s和θ为拟合系数；