本发明涉及增材材料疲劳寿命确定技术领域，具体涉及增材钛合金dfr的确定方法、确定系统、设备及介质，尤其适用于军机。
背景技术：
dfr是结构细节固有的疲劳性能特征值，是一种对构件质量和耐重复载荷能力的度量。激光沉积成形ta15钛合金在航空飞行器结构上得到应用，但其疲劳失效行为和寿命分布特征与传统锻件不同，迫切需要建立寿命可靠性评价方法。在传统的军机dfr分析过程中，n90/99.9可以根据对数正态分布的理论估计得到。通常认为，在疲劳寿命服从对数正态分布的假设下，不同置信水平和可靠度要求下的s-n曲线在对数坐标系下相互平行。因此，dfr值可以通过在特定的应力水平和应力比下的一组成组疲劳试验结果结合标准s-n曲线参数计算得到。但是，对于通过特殊工艺，比如增材制造工艺制备的材料，当疲劳寿命分布发生改变时，需要使用新的分布模型对其疲劳寿命进行描述，此时，dfr值的计算方法也需要进行相应的改进。
技术实现要素：
本发明提供的增材钛合金dfr的确定方法和系统，采用双峰对数正态分布对增材钛合金的疲劳寿命进行描述，并减少了模型不准确造成的过度限制。
本发明一个方面提供了一种增材钛合金dfr的确定方法，其包括：
以预设应力比，在不同应力水平下获取增材钛合金的疲劳寿命数据；
采用双峰对数正态分布模型
利用参数估计em方法，并引入潜在变量z＝(z1,z2,...,zm)，其中z1,z2,...,zm相互独立，m为指定应力水平下成组疲劳试验获得的样本数量，对所述增材钛合金的疲劳寿命数据进行处理，获取双峰对数正态分布的分布参数α，μ1，σ1，μ2，σ2；
获得指定置信区间c要求下的可靠性寿命nc/p，其满足p(np≥nc/p)＝c，np为可靠度p要求下的疲劳寿命；
对增材钛合金疲劳寿命数据拟合得到dfr所要求的置信水平与可靠度下的s-n曲线参数，通过m90/99.9lgsmax lgn90/99.9＝c90/99.9描述s-n曲线，smax为应力水平的应力峰值，m90/99.9和c90/99.9为90％置信水平和99.9可靠度要求下的s-n曲线的曲线参数；
利用
本发明另一方面还提供了一种增材钛合金dfr的确定系统，包括：
获取模块，用于以预设应力比，在不同应力水平下获取增材钛合金的疲劳寿命数据；
处理模块，用于采用双峰对数正态分布模型
利用
本文提供的增材钛合金dfr确定方法，基于双峰对数分布模型，发展了适用于增材钛合金的军机dfr值计算方法，采用双峰对数正态分布描述增材钛合金疲劳寿命分布，依据em算法，建立了分布参数估计方法；并给出了指定置信度和可靠度要求的疲劳寿命计算方法。与对数正态分布相比，双峰对数正态分布模型更加精确的描述了增材钛合金的疲劳寿命分布，高可靠度和高置信度要求下的疲劳寿命估计值明显提高。另一方面，本发明给出了基于双峰对数正态分布的增材钛合金军机dfr值估计方法，该方法可以在一定程度上减少因为描述模型不精确导致的对设计许用应力的过度限制，提高增材钛合金的使用潜力。
附图说明
图1(a)为720mpa应力水平下的概率密度函数曲线；
图1(b)为760mpa应力水平下的概率密度函数曲线；
图1(c)为800mpa应力水平下的概率密度函数曲线。
图2(a)为720mpa应力水平下的90％置信区间不同可靠度下的疲劳寿命；
图2(b)为760mpa应力水平下90％置信区间不同可靠度下的疲劳寿命；
图2(c)为800mpa应力水平下90％置信区间不同可靠度下的疲劳寿命。
图3(a)为720mpa应力水平下的可靠度99.9％不同置信区间下的疲劳寿命；
图3(b)为760mpa应力水平下可靠度99.9％不同置信区间下的疲劳寿命；
图3(c)为800mpa应力水平下可靠度99.9％不同置信区间下的疲劳寿命；
图4为置信区间c＝90％,可靠度p＝99.9％下双峰对数正态分布和对数正态分布两种不同分布模型的s-n曲线。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实例用于说明本发明，但不用来限制本发明的范围。
本发明一个实施例中提供一种增材钛合金的dfr确定方法，包括：
以预设应力比，在不同应力水平下获取增材钛合金的疲劳寿命数据；
采用双峰对数正态分布模型
利用参数估计em方法，并引入潜在变量z＝(z1,z2,...,zm)，其中z1,z2,...,zm相互独立，m为指定应力水平下成组疲劳试验获得的样本数量，对所述增材钛合金的疲劳寿命数据进行处理，获取双峰对数正态分布的分布参数α，μ1，σ1，μ2，σ2；
获得指定置信区间c要求下的可靠性寿命nc/p，其满足p(np≥nc/p)＝c，np为可靠度p要求下的疲劳寿命；
对增材钛合金疲劳寿命数据拟合得到dfr所要求的置信水平与可靠度下的s-n曲线参数，通过m90/99.9lgsmax lgn90/99.9＝c90/99.9描述s-n曲线，smax为应力水平的应力峰值，m90/99.9和c90/99.9为90％置信水平和99.9可靠度要求下的s-n曲线的曲线参数；
利用
采用双峰对数正态分布描述增材钛合金疲劳寿命分布，依据em算法，建立了分布参数估计方法；并给出了指定置信度和可靠度要求的疲劳寿命计算方法。与对数正态分布相比，双峰对数正态分布模型更加精确的描述了增材钛合金的疲劳寿命分布，高可靠度和高置信度要求下的疲劳寿命估计值明显提高。另一方面，给出了基于双峰对数正态分布的增材钛合金军机dfr值估计方法，该方法可以在一定程度上减少因为描述模型不精确导致的对设计许用应力的过度限制，提高增材钛合金的使用潜力。
常规来说，对于军用飞机结构，dfr代表应力比r＝0.1时，在90％置信水平和99.9％可靠度要求下，疲劳寿命达到50000此循环对应的最大应力。根据上述定义，dfr值可以通过90％置信水平和99.9％可靠度要求下测试计算得到的r＝0.1的s-n曲线计算得到。
n90/99.9是计算dfr值的关键中间变量。是否能够精确的估计n90/99.9将直接影响dfr分析的准确性与可靠性。因此，选择合适的分布模型对材料的疲劳寿命分布情况进行描述对于疲劳可靠性评估至关重要。在军用飞机结构疲劳评估中，通常采用对数正态分布模型来描述材料的疲劳寿命分布情况，记为
式中，疲劳寿命用随机变量n表示，μlg为对数期望，σlg为对数寿命标准差。这些参数可以采用极大似然估计法(maximumlikelihoodestimation,mle)估计得到。在军用飞机dfr的过程中，通常认为σlg是一个仅与材料有关的常数。
在传统的军机dfr分析过程中，n90/99.9可以根据对数正态分布的理论估计得到。通常认为，在疲劳寿命服从对数正态分布的假设下，不同置信水平和可靠度要求下的s-n曲线在对数坐标系下相互平行。因此，dfr值可以通过在特定的应力水平和应力比下的一组成组疲劳试验结果结合标准s-n曲线参数计算得到。但是对于增材钛合金来说，探索采用双峰对数正态分布对疲劳寿命进行描述。dfr值的计算方法也需要进行相应的改进。
接下来对增材钛合金的dfr确定方式进行详细说明：
以ded钛合金为例，估计该钛合金的疲劳寿命。
步骤110，获取疲劳寿命测试数据
试验在instron8801-100kn上进行，试验频率为10hz。三种应力水平的应力峰值(smax)分别为720mpa、760mpa和800mpa，共计获得61件有效数据，见表1。三种应力水平下均出现了混合失效模式，所有试件可以根据裂纹源的类型分为两类，即内部缺陷起裂(si)和表面起裂(ss)。
表1疲劳寿命数据
将对数疲劳寿命绘制成对数水平频率分布直方图(frequencydistributionhistogram,fdf)。从图中可以看出，在对数坐标下，ded-ta15钛合金的疲劳寿命分布具有典型的双峰分布特征。
步骤120，构建双峰对数正态分布的概率密度函数和累积分布函数
假设疲劳寿命n服从双峰对数正态分布模型，用
上面的pdf描述了cdf的变化趋势，即曲线的斜率。
步骤130，采用em方法进行参数估计
令y＝lgn，则y服从双峰正态分布，记为
引入潜在变量z＝(z1,z2,...,zm)，其中z1,z2,...,zm相互独立，m为指定应力水平下成组疲劳试验获得的样本数量，有，
式中，i＝1,2,...,m。z满足如(10)所示概率分布函数。
令x＝(y,z)，则有如式(11)所示似然函数，
式中ψ＝(α，μ1，σ1，μ2，σ2)，对式(12)取对数并去掉与所估计参数无关项，有：
令
容易验证，其中有，
在m步中，解下列方程
解得：
式(16)即为双峰对数正态分布参数估计的迭代格式，k为迭代次数，i为第i个样本。对于表1中数据进行统计分析，根据式(14)和式(16)估计得到的双峰对数正态分布的分布参数见表2。
表2双峰对数正态分布参数估计结果
步骤140，基于双峰对数正态分布的dfr方法
首先，估计n90/99.9
可靠度p要求下的疲劳寿命np满足式(17)
p(n≥np)＝p(17)
且有，
对于指定置信水平c要求下的可靠性寿命nc/p满足式
p(np≥nc/p)＝c(19)
由于双峰对数正态分布的形式较为复杂，无法方便地通过解析方法获得nc/p。本文采用bootstrap方法进行nc/p的估计。过程如下：
①从原始疲劳寿命样本集{n1,n2,…,nm}中有放回随机抽样，得到一个新的疲劳寿命样本集
②对于
③根据式(17)，计算指定可靠度要求下的疲劳寿命
④重复若干次(本文取5000次)①到③过程，获得若干个
当置信水平c＝0.9和可靠度p＝0.999，上文nc/p即为所求dfr计算过程中的关键中间量，n90/99.9。根据上述方法，估计得到的三个应力水平的n90/99.9如表3所示。如图2(a)-图2(c)为90％置信区间不同可靠度下的疲劳寿命，图3(a)-图3(c)为99.9％可靠度不同置信水平下的疲劳寿命。
表3双峰对数正态分布下n90/99.9估计结果
接下来描述s-n曲线：
在计算dfr值之前，通过进行几组不同应力水平下的多组疲劳试验测试结果拟合得到dfr计算所要求的置信水平与可靠度下的s-n曲线参数。通常，在军机dfr分析中，要求置信水平c＝0.9，可靠度p＝0.999。
在疲劳寿命基于双峰对数正态分布的假设下，s-n曲线仍然可以采用basquin方程进行描述。其对数形式如式(20)所示。
m90/99.9lgsmax lgn90/99.9＝c90/99.9(20)
式中，m90/99.9和c90/99.9为90％置信水平和99.9可靠度要求下的s-n曲线的曲线参数，可以采用最小二乘法拟合得到，如式(21)所示。
其中，smax,i表示第i组成组疲劳试验的应力水平，n90/99.9,i表示估计得到的该应力水平下，90％置信水平、99.9％可靠度要求下的疲劳寿命，l为第i个应力水平。
接着，计算dfr值
根据上述方法计算得到s-n曲线，可以计算得到90％置信水平、99.9％可靠度要求下，指定疲劳寿命nd对应的疲劳应力峰值
在军机dfr分析中，dfr值代表应力比r＝0.1时，在90％置信水平和99.9％可靠度要求下，疲劳寿命达到50000次循环对应的最大应力。因此还需要根据等寿命曲线，将nd＝50000次循环对应的
式中，σm0为等寿命曲线与横轴交点的横坐标，本文取σm0＝620mpa。根据上述方法，计算得到的90％置信水平和99.9％可靠度要求下的s-n曲线参数，nd＝50000次循环对应的
表4基于双峰对数正态分布的s-n曲线参数和dfr值
另外，本领域内的技术人员应当理解的是，在本发明实施例的申请文件中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、物品或者设备中还存在另外的相同要素。
本发明实施例的说明书中，说明了大量具体细节。然而应当理解的是，本发明实施例的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。类似地，应当理解，为了精简本发明实施例公开并帮助理解各个发明方面中的一个或多个，在上面对本发明实施例的示例性实施例的描述中。
然而，并不应将该公开的方法解释成反映如下意图：即所要求保护的本发明实施例要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说，如权利要求书所反映的那样，发明方面在于少于前面公开的单个实施例的所有特征。因此，遵循具体实施方式的权利要求书由此明确地并入该具体实施方式，其中每个权利要求本身都作为本发明实施例的单独实施例。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
1.一种增材钛合金dfr确定方法，其特征在于，包括：
以预设应力比，在不同应力水平下获取增材钛合金的疲劳寿命数据；
采用双峰对数正态分布模型
利用参数估计em方法，并引入潜在变量z＝(z1,z2,...,zm)，其中z1,z2,...,zm相互独立，m为指定应力水平下成组疲劳试验获得的样本数量，对所述增材钛合金的疲劳寿命数据进行处理，获取双峰对数正态分布的分布参数α，μ1，σ1，μ2，σ2；
获得指定置信区间c要求下的可靠性寿命nc/p，其满足p(np≥nc/p)＝c，np为可靠度p要求下的疲劳寿命；
对增材钛合金疲劳寿命数据拟合得到dfr所要求的置信水平与可靠度下的s-n曲线参数，通过m90/99.9lgsmax lgn90/99.9＝c90/99.9描述s-n曲线，smax为应力水平的应力峰值，m90/99.9和c90/99.9为90％置信水平和99.9可靠度要求下的s-n曲线的曲线参数；
利用
2.根据权利要求1所述的增材钛合金dfr的确定方法，其特征在于，所述双峰对数正态分布模型的概率密度函数和累积分布函数分别为：令y＝lgn，则y服从双峰正态分布，
n为疲劳寿命，μlg为对数期望，σlg为对数寿命标准差。
3.根据权利要求1所述的增材钛合金dfr的确定方法，其特征在于，所述获取双峰对数正态分布的分布参数的方法包括：
利用
对所述增材钛合金的疲劳寿命数据进行处理，获取双峰对数正态分布的分布参数α，μ1，σ1，μ2，σ2，k为迭代次数，i为第i个样本。
4.根据权利要求1所述的增材钛合金dfr的确定方法，其特征在于，所述nc/p的计算方法包括：
从原始疲劳寿命样本集{n1,n2,…,nm}中有放回随机抽样，得到一个新的疲劳寿命样本集
对于
根据p(n≥np)＝p，计算指定可靠度要求下的疲劳寿命
重复多次以上步骤，获得多个对应的
5.根据权利要求1所述的增材钛合金dfr的确定方法，其特征在于，样本抽取的方式为：随机又放回得抽取。
6.根据权利要求1-5任一项所述的增材钛合金dfr的确定方法，其特征在于，所述m90/99.9和所述c90/99.9的获取方法为：
7.一种权利要求1-6任一项所述增材钛合金dfr确定方法的确定系统，其特征在于，包括：
获取模块，用于以预设应力比，在不同应力水平下获取增材钛合金的疲劳寿命数据；
处理模块，采用双峰对数正态分布模型
8.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机介质，其特征在于，所述处理器执行计算机程序时实现权利要求1～6中任意一项所述方法的步骤。
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实施权利要求1～6中任意一项所述方法的步骤。
技术总结