单晶高温合金的组织和缺陷的模拟方法与流程

本发明涉及材料技术领域，具体而言，涉及一种单晶高温合金组织和缺陷的模拟方法。

背景技术：

单晶高温合金是航空发动机叶片的首选材料，其变形和损伤机理本质上表现为由微观组织(主要为γ/γ'相)和缺陷(主要为空位和位错)的相互作用，因此，定量描述单晶高温合金组织和缺陷的同时演变，是正确揭示叶片变形和损伤机理的前提和基础。

航空发动机运行时，叶片在高温低应力环境中变形，现有实验技术手段只能观察变形中断后的静态组织和缺陷，无法实现观测高温变形过程中的组织和缺陷，需借助计算模拟手段。而目前没有一种模拟方法能够模拟高温变形过程中真实的γ/γ'相演化，并耦合空位和位错，因此无法实现组织和缺陷的同时演变。

需要说明的是，在上述背景技术部分发明的信息仅用于加强对本发明的背景的理解，因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。

技术实现要素：

本发明的目的在于提供一种单晶高温合金组织和缺陷的模拟方法，解决现有模拟方法不能同时模拟高温变形过程中组织和缺陷的问题。

根据本发明的一个方面，提供一种单晶高温合金的组织和缺陷的模拟方法，其特征在于，包括：

确定模拟物理量：以所述单晶高温合金γ/γ'相中掺杂元素i的摩尔分数、第k个序参量、位错密度、空位密度为模拟物理量；

根据所述模拟物理量，利用如下公式计算所述单晶高温合金的自由能，

式中，f表示自由能，f^chem表示化学自由能，f^el表示弹性能；表示摩尔分数，表示第k个序参量，表示滑移系m的正刃型位错密度，表示滑移系m的负刃型位错密度，表示空位密度；

根据所述自由能，利用模拟方程模拟所述单晶高温合金的组织和缺陷；其中，利用模拟方程模拟所述单晶高温合金的组织和缺陷包括：利用描述原子扩散的模拟方程模拟成分变化，利用描述原子局部有序-无序的转变的模拟方程模拟序参量变化，利用描述位错密度的模拟方程模拟位错变化，利用描述位错攀移的模拟方程模拟空位变化。

在本发明的一种示例性实施方式中，所述化学自由能包括体积自由能和梯度自由能。

在本发明的一种示例性实施方式中，所述计算单晶高温合金的自由能包括：利用如下方程计算所述体积自由能：

式中，f^bulk表示体积自由能，表示γ相的自由能，表示γ'相的自由能，为插值函数。

在本发明的一种示例性实施方式中，所述计算单晶高温合金的自由能包括：利用如下方程计算所述梯度自由能：

式中，f^grad表示梯度自由能，kφ表示梯度自由能系数。

在本发明的一种示例性实施方式中，所述计算单晶高温合金的自由能包括：利用如下公式计算所述弹性能：

式中，ε^el表示弹性应变，σ表示应力，ε^el和σ满足胡克定律。

在本发明的一种示例性实施方式中，所述计算弹性能包括利用如下方程计算所述弹性应变：

ε^el＝ε+ε^inel

式中，ε为总应变，ε^inel为非弹性应变，所述总应变根据公式计算；所述非弹性应变包括由γ/γ'共格界面引起的本征应变、由位错引起的本征应变和由空位引起的本征应变。

在本发明的一种示例性实施方式中，所述所述利用描述原子扩散的模拟方程模拟成分变化包括：利用如下模拟方程模拟成分变化：

在本发明的一种示例性实施方式中，所述所述利用描述原子局部有序-无序的转变的模拟方程模拟序参量变化包括：利用如下模拟方程模拟序参量变化：

在本发明的一种示例性实施方式中，所述所述利用描述位错密度的模拟方程模拟位错变化包括：利用如下模拟方程模拟位错变化：

其中，为位错滑移速度，为位错攀移速度，为位错的增殖和湮灭。

在本发明的一种示例性实施方式中，所述利用描述位错攀移的模拟方程模拟空位变化包括：利用如下模拟方程模拟空位变化：

本发明的模拟方法以介观相场法为框架，在该框架内同时描述单晶高温合金变形时γ/γ’相、空位和位错的演变，在介观尺度内实现了组织和缺陷的共同演变，无需跨尺度，不仅大大降低了建模和数值运算的工作量，而且能够解析组织和缺陷的相互作用机制，进而揭示航空发动机叶片的变形和损伤机理。

应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本发明。

具体实施方式

现在将更全面地描述示例实施方式。然而，示例实施方式能够以多种形式实施，且不应被理解为限于在此阐述的实施方式；相反，提供这些实施方式使得本发明将全面和完整，并将示例实施方式的构思全面地传达给本领域的技术人员。

相关技术中，分子动力学以单个原子为基本模拟单位，可描述极少量空位和位错的演变。离散位错动力学不再分辨单个原子，以单根位错线段为基本模拟单位，可描述少量位错的演变。不同尺度相场法可描述空位、位错和γ/γ'相的演变。不同尺度相场法可描述空位、位错和γ/γ'相的演变。

但分子动力学和离散位错动力学的缺陷主要是：1)空间尺度太小，所能描述的位错数量远小于实验所测实际情况；2)时间尺度太小，无法描述γ/γ'相的演变。水平集方法的主要缺陷有：1)跟踪明锐特性使得计算效率很低；2)无法描述空位和位错。原子尺度相场方法(即晶体相场法)的优劣与分子动力学类似。介观尺度相场法虽能模拟真实的γ/γ'相演化，但还没有可以耦合介观空位和位错模型的相场方法，因而无法实现组织和缺陷的同时演变。

本发明实施方式中提供了一种单晶高温合金组织和缺陷的模拟方法，具备主要优点是：1)空间尺度足够大，可描述与实验所测一致的位错密度；2)时间尺度足够打，单位时间步长达10⁰s，可描述10²h时间内的γ/γ'相的演变；3)以自然形成的扩散界面代替传统明锐界面，无需再跟踪界面坐标，大大提高计算效率；4)适用于高温变形过程中定量描述单晶高温合金组织和缺陷的同时演变。

本发明实施方式的单晶高温合金组织和缺陷的模拟方法，包括：

步骤110，确定模拟物理量：以单晶高温合金γ/γ'相中掺杂元素i的摩尔分数、第k个序参量、位错密度、空位密度为模拟物理量；

步骤210，根据所模拟物理量，利用如下公式计算单晶高温合金的自由能，

式中，f表示自由能，f^chem表示化学自由能，f^el表示弹性能；表示摩尔分数，表示第k个序参量，表示滑移系m的正刃型位错密度，表示滑移系m的负刃型位错密度，表示空位密度；

步骤310，根据自由能，利用模拟方程模拟单晶高温合金的组织和缺陷，包括：利用描述原子扩散的模拟方程模拟成分变化，利用描述原子局部有序-无序的转变的模拟方程模拟序参量变化，利用描述位错密度的模拟方程模拟位错变化，利用描述位错攀移的模拟方程模拟空位变化。

下面对本发明实施方式的单晶高温合金的组织和缺陷的模拟方法进行详细说明：

步骤110中，单晶高温合金γ/γ'相中，γ相为化学无序相，γ'相为化学有序相，例如，在镍基铝合金单晶高温合金的γ/γ'中，γ相为化学无序相，即al原子随机地在ni晶格中占位，γ'相为化学有序相，即al原子在ni晶格中占据特定位置，可能有四种不同的特定占位方式，从而形成四种不同的γ'相变体。因而，γ/γ'相可由掺杂元素i的摩尔分数和第k个序参量描述，序参量表征有序-无序转变。γ相表示为和γ'相的k变体表示为和其中l≠k。和分别为γ相和γ'相的平衡成分。γ相和γ'相的界面(γ/γ'界面)可认为是γ相和γ'相的混合区域，γ'相变体之间的界面(即反相畴界)可认为是γ'相k变体和l变体的混合区域。本发明掺杂元素i指合金中含量较少的元素，如二元合金中，镍基铝合金中的i指代al。

介观尺度不辨别单根位错，而是以位错密度代表位错的存在。以和分别表示滑移系m的正刃型和负刃型位错密度。同理，介观尺度不辨别单个空位，而是以空位密度代表空位的存在。空位无需区分滑移系，可用表示。

在本示例性实施方式中，步骤210中，自由能f中的化学自由能f^chem包括体积自由能f^bulk和梯度自由能f^grad。体积自由能表示离散原子状态相对组合形成均匀平衡相状态的能量，梯度自由能表示相界面处由于非均匀成分/结构状态相对相内均匀成分/结构状态的能量。

在本示例性实施方式中，体积自由能f^bulk利用如下方程计算：

式中，为γ相的自由能，为γ'相的自由能，为插值函数，以把γ/γ'界面作为γ和γ'相的混合；表示体积自由能对γ/γ'界面能和反相畴界能的贡献，其确切数值可通过ω系数调控。

在本示例性实施方式中，梯度自由能f^grad利用如下方程计算：

f^grad是为了进一步补偿γ/γ'界面能和反相畴界能，其数值通过梯度自由能系数kφ调控。

在本示例性实施方式中，弹性能f^el利用如下方程计算：

其中，ε^el表示弹性应变，σ表示应力，ε^el和σ遵循胡克定律

在本示例性实施方式中，弹性应变ε^el＝ε+ε^inel，ε为总应变，ε^inel为非弹性应变，其中，在小变形范围内，总应变ε与位移u的关系为ε^inel包含：1)由γ/γ'共格界面引起的本征应变2)由位错引起的本征应变ε^dis＝∑mη^mp^m，3)由空位引起的本征应变ε^vac＝∑mξ^mq^m；其中，b^m为滑移系m的柏氏矢量，n^m为滑移系m的法向量，p^m为滑移m的投影张量，q^m为滑移系m的柏氏矢量张量。

由于弹性波的传播速度远大于组织和缺陷的演变速度，因此，可以假设力平衡瞬间完成，便可算得每一步的位移，应变和应力。

在本示例性实施方式中，步骤310中，成分的演变本质上是原子的扩散，因此，模拟方程可以为allen-cahn方程，该方程即为描述原子扩散的介观唯象公式：

在本示例性实施方式中，步骤310中，序参量的演变本质上为原子的局部有序-无序转变，因此，模拟方程为cahn-hilliard方程，该方程即为描述有序-无序转变的介观唯象公式：

在本示例性实施方式中，步骤310中，位错密度的演变实质上为位错的滑移、攀移、增殖、湮灭等，位错滑移和滑移本质上为位错线的平动，数学上用平移项描述，位错增殖/湮灭本质上是位错线的直接增加/减少，数学上用源项描述，因此，位错的模拟方程为。：

其中，为位错滑移速度，为位错攀移速度，为peach-koehler力，为位错的增殖和湮灭。

在本示例性实施方式中，步骤310中，空位的演变本质上与位错的攀移相关联，因此，这一过程的介观公式总结即为空位的模拟方程：

本发明物理基础坚实可靠，相场法所用的演变判断“体系自发地向自由能降低的方向演变”是热力学第二定律在恒压或恒容情况下的推论。在介观尺度构建包含组织和缺陷的物理参量并构建由这些物理参量描述的体系自由能，便可在同一尺度内(介观尺度)实现了组织和缺陷的共同演变，无需跨尺度，不仅大大降低了建模和数值运算的工作量，而且能够解析组织和缺陷的相互作用机制，进而揭示航空发动机叶片的变形和损伤机理。

用语“一个”、“一”、“该”、“所述”和“至少一个”用以表示存在一个或多个要素/组成部分/等；用语“包括”和“具有”用以表示开放式的包括在内的意思并且是指除了列出的要素/组成部分/等之外还可存在另外的要素/组成部分/等。

本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本发明的真正范围和精神由所附的权利要求指出。