CN112949020A - 增材钛合金的疲劳寿命的确定方法、系统、设备与介质 - Google Patents
增材钛合金的疲劳寿命的确定方法、系统、设备与介质 Download PDFInfo
- Publication number
- CN112949020A CN112949020A CN202011455131.0A CN202011455131A CN112949020A CN 112949020 A CN112949020 A CN 112949020A CN 202011455131 A CN202011455131 A CN 202011455131A CN 112949020 A CN112949020 A CN 112949020A
- Authority
- CN
- China
- Prior art keywords
- fatigue life
- distribution
- titanium alloy
- bimodal
- reliability
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 85
- 229910001069 Ti alloy Inorganic materials 0.000 title claims abstract description 50
- 239000000654 additive Substances 0.000 title claims abstract description 38
- 230000000996 additive effect Effects 0.000 title claims abstract description 38
- 238000009826 distribution Methods 0.000 claims abstract description 183
- 230000002902 bimodal effect Effects 0.000 claims abstract description 53
- 238000005315 distribution function Methods 0.000 claims abstract description 23
- 238000005070 sampling Methods 0.000 claims abstract 2
- 230000006870 function Effects 0.000 claims description 25
- 238000012360 testing method Methods 0.000 claims description 12
- 238000001514 detection method Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 7
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 241000287196 Asthenes Species 0.000 claims description 2
- 238000004590 computer program Methods 0.000 claims 2
- 238000005457 optimization Methods 0.000 claims 2
- 239000013256 coordination polymer Substances 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 10
- 239000000463 material Substances 0.000 abstract description 5
- 238000004458 analytical method Methods 0.000 description 9
- 150000001875 compounds Chemical class 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000009661 fatigue test Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- TZCXTZWJZNENPQ-UHFFFAOYSA-L barium sulfate Chemical compound [Ba+2].[O-]S([O-])(=O)=O TZCXTZWJZNENPQ-UHFFFAOYSA-L 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及增材疲劳寿命技术领域,具体涉及增材钛合金的疲劳寿命的确定方法、系统、设备与介质。该方法包括:设定疲劳寿命数据服从双峰威布尔分布,构建双峰威布尔分布模型的概率密度函数和分布函数;基于ECM估值法,在求解非线性方程组时采用逐步最大化,得到双峰威布尔分布的显式解耦迭代式,进而求得分布参数π,α1,β1,α2,β2;对数据抽取样本,估计双峰威布尔分布的参数,并计算N*可靠度要求为P的可靠度寿命多次计算,由获得C=95%,P=95%对应疲劳寿命NP,C。采用疲劳寿命数据服从双峰威布尔分布,能够更有效、更准确得对疲劳寿命的分布进行描述;确定指定置信度和可靠度要求的疲劳寿命,结果更精确,提高了高可靠度和高置信度要求下的疲劳寿命估计值。
Description
技术领域
本发明涉及增材疲劳寿命技术领域,具体涉及增材钛合金的疲劳寿命的确定方法、系统、设备与介质。
背景技术
激光沉积成形TA15钛合金在航空飞行器结构上得到应用,但其疲劳失效行为和寿命分布特征与传统锻件不同,迫切需要建立寿命可靠性评价方法。而且DFR是结构细节固有的疲劳性能特征值,是一种对构件质量和耐重复载荷能力的度量。某置信区间和可靠度下的疲劳寿命是计算DFR值的关键中间变量。是否能够精确的估计N90/99.9将直接影响DFR分析的准确性与可靠性。因此,选择合适的分布模型对材料的疲劳寿命分布情况进行描述对于疲劳可靠性评估至关重要。
技术内容
本发明提供了一种增材钛合金的疲劳寿命的确定方法和确定系统,以提高增材钛合金疲劳寿命数据准确描述,以及估算增材钛合金疲劳寿命的准确度。
本发明一个方面,提供一种增材钛合金的疲劳寿命的确定方法,其包括:
以预设应力比,在不同应力水平下获取增材钛合金的疲劳寿命数据;
设定所述疲劳寿命数据服从双峰威布尔分布,构建概率密度函数f(N)和分布函数F(N);
基于ECM估值法,在求解非线性方程组时采用逐步最大化,得到双峰威布尔分布的显式解耦迭代式,进而求得π,α1,β1,α2,β2,其中,π,α1,β1,α2,β2为概率密度函数和分布函数中的分布参数,其中α1,α2为形状参数,β1,β2为范围参数,π为权重参数,其取值范围为0≤π≤1;
本发明另一方面还提供了一种执行所述增材钛合金的疲劳寿命的确定系统,其包括:
获取模块,获取以预设应力比,在不同应力水平下增材钛合金的疲劳寿命数据;
处理模块,设定所述疲劳寿命数据服从双峰威布尔分布,构建概率密度函数f(N)和分布函数F(N);
基于ECM估值法,在求解非线性方程组时采用逐步最大化,得到双峰威布尔分布的显式解耦迭代式,进而求得π,α1,β1,α2,β2,其中,π,α1,β1,α2,β2为概率密度函数和分布函数中的分布参数,其中α1,α2为形状参数,β1,β2为范围参数,π为权重参数,其取值范围为0≤π≤1;
DEDTA15钛合金疲劳寿命分布呈双峰特性的现象,在基于Weibull分布的DFR分析方法的基础上,扩展到双峰威布尔BioWeibull分布,基于BioWeibull分布的分析方法:提出DEDTA15钛合金疲劳寿命服从双峰Weibull分布。基于EM方法,对M步骤进行改进,形成ECM方法,ECM方法原理给出了BioWeibull分布的参数估计方法。验证对比,与Weibull分布相比,BioWeibull分布能够更有效的对疲劳寿命的分布形式进行描述;给出了指定置信度和可靠度要求的疲劳寿命计算方法,通过算例得出采用BioWeibull分布对疲劳寿命进行更加精确的描述,相对与Weibull分布,BioWeibull分布下可靠度和高置信度要求下的疲劳寿命估计值均有所提高。
附图说明
图1(a)为720MPa应力水平下的疲劳寿命频率分布直方图;
图1(b)为760MPa应力水平下的疲劳寿命频率分布直方图;
图1(c)为800MPa应力水平下的疲劳寿命频率分布直方图;
图2(a)为720MPa应力水平下FDH和根据参数估计结果绘制的BioWeibull分布和Weibull分布的PDF曲线;
图2(b)为760MPa应力水平下FDH和根据参数估计结果绘制的BioWeibull分布和Weibull分布的PDF曲线;
图2(c)为800MPa应力水平下FDH和根据参数估计结果绘制的BioWeibull分布和Weibull分布的PDF曲线;
图3(a)为720MPa应力水平下根据参数估计结果绘制的BioWeibull分布和Weibull分布的CDF曲线;
图3(b)为760MPa应力水平下根据参数估计结果绘制的BioWeibull分布和Weibull分布的CDF曲线;
图3(c)为800MPa应力水平下根据参数估计结果绘制的BioWeibull分布和Weibull分布的CDF曲线;
图4(a)为720MPa应力水平下95%置信度,不同可靠度下的威布尔分布(Weibull)和双峰威布尔分布(BioWeibull)下疲劳寿命值;
图4(b)为760MPa应力水平下95%置信度,不同可靠度下的威布尔分布(Weibull)和双峰威布尔分布(BioWeibull)下疲劳寿命值;
图4(c)为800MPa应力水平下95%置信度,不同可靠度下的威布尔分布(Weibull)和双峰威布尔分布(BioWeibull)下疲劳寿命值;
图5(a)为720MPa应力水平下95%可靠度,不同置信度下的威布尔分布(Weibull)和双峰威布尔分布(BioWeibull)下疲劳寿命值;
图5(b)为760MPa应力水平下95%可靠度,不同置信度下的威布尔分布(Weibull)和双峰威布尔分布(BioWeibull)下疲劳寿命值;
图5(c)为800MPa应力水平下95%置可靠度,不同置信度下的威布尔分布(Weibull)和双峰威布尔分布(BioWeibull)下疲劳寿命值;
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实例用于说明本发明,但不用来限制本发明的范围。
本发明一个实施例中提供一种增材钛合金的疲劳寿命的确定方法,其包括:
以预设应力比,在不同应力水平下获取增材钛合金的疲劳寿命数据;
设定所述疲劳寿命数据服从双峰威布尔分布,构建概率密度函数f(N)和分布函数F(N);
基于ECM估值法,在求解非线性方程组时采用逐步最大化,得到双峰威布尔分布的显式解耦迭代式,进而求得π,α1,β1,α2,β2,其中,π,α1,β1,α2,β2为概率密度函数和分布函数中的分布参数,其中α1,α2为形状参数,β1,β2为范围参数,π为权重参数,其取值范围为0≤π≤1;
DEDTA15钛合金疲劳寿命分布呈双峰特性的现象,在基于Weibull分布的DFR分析方法的基础上,扩展到双峰威布尔分布,基于BioWeibull分布的DFR分析方法:
(1)提出DEDTA15钛合金疲劳寿命服从双峰Weibull分布。基于EM方法,对M步骤进行改进,形成ECM方法,ECM方法原理给出了BioWeibull分布的参数估计方法。验证对比,与Weibull分布相比,BioWeibull分布能够更有效的对疲劳寿命的分布形式进行描述。
(2)给出了指定置信度和可靠度要求的疲劳寿命计算方法,通过算例得出采用BioWeibull分布对疲劳寿命进行更加精确的描述,相对与Weibull分布,高BioWeibull分布下可靠度和高置信度要求下的疲劳寿命估计值均有所提高。
(3)给出了指定置信度和可靠度要求的C-P-S-N拟合方法,并基于DFR的定义,给出了基于BioWeibull分布的DFR计算方法。通过BioWeibull分布对数据的分布情况进行更加准确的描述,可以在一定程度上减少因为描述模型不精确导致的对设计许用应力的过度限制,有效的提高材料或结构的使用潜力。
接下来对增材钛合金的疲劳寿命的确定方式进行详细说明:
DFR方法是一种有效的寿命评定方法,DFR值是核心参数,迫切需要建立DFR值测定方法。
DFR值是结构在承受R=0.06的等幅交变载荷时,在95%的可靠度要求和95%的置信度要求下,结构疲劳寿命达到100000循环时的最大名义应力。其本质上结构疲劳强度的固有特征,与载荷环境无关,同时也考虑了和疲劳强度相关的不确定性因素。
以DED钛合金为例,估计该钛合金的疲劳寿命。
步骤110,疲劳寿命测试
进行应力比为R=0.06的三个应力水平下总计61件标准圆棒试件的疲劳,三种应力水平的应力峰值Smax分别为720MPa,760MPa和800MPa,每种应力水平疲劳寿命数据样本容量分别为18,21和22,具体的疲劳寿命数据见表1。
表1
对上表构建三种应力水平下的疲劳寿命频率分布直方图,见图1(a)-图1(c)从图中可以看出三种应力水平下,DED钛合金的疲劳寿命都呈现一定的双峰特性。
步骤120,采用双峰威布尔分布构建概率密度函数和分布函数
通常,采用双峰威布尔分布的对疲劳寿命N进行描述,记为N~W(α,β),其概率密度函数和分布函数见Eq.1和Eq.2。在weibull分布的基础上,采用线性加权的方式构造双峰Weibull分布对呈现双峰分布特性的疲劳寿命数据的分布情况进行描述。
α和β为概率密度函数、分布函数的参数。
假设疲劳寿命N服从双峰Weibull分布,记为N~BW(π,α1,β1,α2,β2),其概率密度函数f(N)和分布函数F(N)分别为Eq.3和Eq.4。
式中,π,α1,β1,α2,β2为分布参数,其中α1,α2为形状参数,β1,β2为范围参数,π为权重参数,其取值范围为,0≤π≤1。f1(N),F1(N),和f2(N),F2(N)分别为两个服从威布尔分布的随机变量的概念密度函数和分布函数。
步骤130,进行参数估计
在本文中以EM算法的核心思想为基础,只是在求解非线性方程组时采用逐步最大化的思想,避免求解复杂的非线性方程组。接下来将从EM算法和EM改进算法(ECM算法)说明参数的估计方法。
EM方法
EM算法是建立在最大似然估计法之上,利用迭代的手段求解复杂分布函数的参数估计值的一种数值求解方法。其核心思想是数据扩张,主要过程是在一步迭代中,分为E步与M步,在E步求对数似然函数在未知参数Ψ期望,在M步将期望极大化。
对于增材钛合金的双峰WEIBULL分布,设N={N1,N2......Nm}为来自于总体NΩ的样本,样本值为N={n1,n2......nm},m为N的样本容量。设分布律中的所有未知参数组成的向量为Ψ=(α1,α2,β1,β2,π)采用极大似然估计法估计未知参数,有似然函数为:
有对数似然函数为:
定义潜在变量Z=(z1,z2,......zm),其中z1,z2,.......zm相互独立,有如下概率分布函数,
且有,当zi=1时,ni~W1(α1,β1),当zi=0时,ni~W2(α2,β2)。
设X=(N,Z),即有xi=(ni,zi)的似然函数为:
对似然函数取对数得:
设在第k+1步迭代中,有估计值ψ(k),由EM算法的E步和M步得到新的估计值ψ(k+1),k为迭代次数,i为第i个样本数。
在E步中,令
容易验证,有
在M步中,寻找期望Q(Ψ;Ψ(k))的最大值,有:
上述方程组即为EM算法求解五个参数的迭代格式。
ECM方法
在求解中,若直接将迭代格式Eq.13交给计算机处理,会遇到诸多问题。一是算法可能会陷入局部最优。二是在求解非线性方程组时由于初值敏感性的原因很有可能会得到错误的结果,在上百步的迭代中,若出现某一次求解错误就需要重新选定当前迭代下的方程组求解初值,使用程序处理起来会非常麻烦,为此,在EM算法的基础上,对M法给出改进——ECM算法。
ECM算法是以EM算法的核心思想为基础,只是在求解非线性方程组时采用逐步最大化的思想,避免求解复杂的非线性方程组。设在第k+1步迭代中,给定求使得Q取极大值。再给定求使得Q取极大值。依次求出完成第K+1步迭代。
进而得到,
由Eq.11可得,
使用一步牛顿迭代法将两个超越方程转化为两个显式的迭代格式,如下:
记(j=1,2)
由
其中g为超越方程函数。
得
由此得到双重二参数威布尔分布的显式解耦的迭代格式。
步骤140,计算置信区间C,可靠度P对应的NP,C
DFR方法作为一种基于PSN曲线和等寿命曲线的疲劳分析方法,已经被广泛地应用于民用飞机结构的疲劳分析和设计。其适用性和便捷性已经被大量的试验和应用所证明。DFR值是结构在承受R=0.06的等幅交变载荷时,在95%的可靠度要求和95%的置信度要求下,结构疲劳寿命达到100000循环时的最大名义应力。其本质上结构疲劳强度的固有特征,与载荷环境无关,同时也考虑了和疲劳强度相关的不确定性因素。
通常认为,在疲劳寿命服从双峰威布尔分布的前提下,对于同种材料,在不同的置信水平和可靠度要求下的P-S-N曲线相互平行。
根据DFR的定义,可以通过进行Rs=0.06,不同应力水平下的成组疲劳试验计算得到不同应力水平下的疲劳寿命分布参数,进而计算得到不同应力水平下,95%的可靠度要求和95%的置信度要求的疲劳寿命N95/95,拟合得到双95要求下的S-N曲线,再根据该曲线可以计算得到DFR。
根据威布尔分布的相关理论,N95/95的计算公式可以按下式计算,
对于双峰威布尔分布来说,由于分布形式较为复杂并且缺乏分布理论的基础,因此需要采用数值方法计算双峰威布尔分布下的N95/95。
指定可靠度P要求下的疲劳寿命记为NP,满足eq.21
P(N≥NP)=P, eq.21
有
对于给定的置信度C,NP,C表示NP置信度要求为C的置信下限,满足eq.23
P(NP,E≥NP,C)=C, eq.23
其中,NP,E表示NP的真实值。
由于BW模型和参数估计方法复杂,因此无法通过解析方法轻松获得。因此采用Bootstrap方法估计BW分布的参数,步骤如下
求取置信水平C和可靠度P对应的疲劳寿命NP,C。当C=95%,P=95%时,计算得到的NP,C即为N95/95。
步骤150,拟合优度检验
采用Anderson–Darling(AD)法对BioWeibull分布和Weibull分布的参数估计结果进行拟合优度检验,其主要步骤如下:
(1)假设样本数据服从X分布,将所检验的样本数据和估计得到分布参数带入Eq.25,计算得到检测统计量An0。
(2)根据Bootstrap的方法,对样本进行再抽样,并估计X分布的分布参数。
(3)根据(2)中估计得到的分布参数,基于原始样本数据计算检测统计量An。
(4)重复(2)(3)步骤5000次,获得5000个检测统计量An的序列。
(5)将序列从小至大排序,获得有序序列其中k表示排序后的序号,则给定置信度c下的分位点为检验临界值,本文中c=95%。若(1)中检测统计量An0大于Aα,则假设分布X不成立,反之则接受分布X为该寿命分布的模型。
本发明还提供了一种执行上述增材钛合金的疲劳寿命的确定系统,其包括:
获取模块,获取以预设应力比,在不同应力水平下增材钛合金的疲劳寿命数据;
处理模块,设定所述疲劳寿命数据服从双峰威布尔分布,构建概率密度函数f(N)和分布函数F(N);
基于ECM估值法,在求解非线性方程组时采用逐步最大化,得到双峰威布尔分布的显式解耦迭代式,进而求得π,α1,β1,α2,β2,其中,π,α1,β1,α2,β2为概率密度函数和分布函数中的分布参数,其中α1,α2为形状参数,β1,β2为范围参数,π为权重参数,其取值范围为0≤π≤1;
接下来以DLD-TA15钛合金为例,获得DLD-TA15钛合金等幅谱下成组疲劳试验结果,如表1,对步骤120、步骤130所给出的BW分布参数估计方法和基于BW分布发展的DFR方法进行验证。同时,为了进一步分析基于BW的DFR方法与传统DFR方法的异同,也采用传统的威布尔分布和基于威布尔分布DFR方法对上述数据进行了处理。
采用步骤130中所给出的ECM方法估计3种应力水平下疲劳寿命的双峰weibull分布参数,ECM方法迭代所用各个参数的初值可以通过图表法初步估计给出。疲劳寿命数据见
表1,估计得到BW分布参数见表2。采用图表法估计得到了3种应力水平下疲劳寿命数据服从单峰weibull分布的分布参数,列入表3。
表2 Estimation Result of Bio-Weibull Distribution Parameters
表3 Estimation Result of Weibull Distribution Parameters
图2(a)-图2(c)为FDH和根据参数估计结果绘制的BioWeibull分布和Weibull分布的PDF曲线。从图中可以看出,与Weibull分布相比,采用BioWeibull分布对算例所用的疲劳寿命数据进行描述,可以更好的反应数据的分布规律,并且能够在一定程度上反应疲劳寿命数据的双峰分布特性。
图3(a)-图3(c)为根据参数估计结果绘制的BioWeibull分布和Weibull分布的CDF曲线,并一并给出了根据秩统计理论求得的疲劳寿命数据对应的累积概率PR。PR可以通过Eq.24计算得到。
式中,n为某一应力水平下疲劳寿命数据的样本数量,i为对全本样本按升序排序后,某个疲劳寿命样本在全部n个样本中的序号。从中可以看出,与Weibull分布相比,BioWeibull能够更好的拟合疲劳寿命数据,更好的反应数据的分布情况。
从图2(a)-图2(c)和图(a)-图3(c)三种,可以定性认为BioWeibull分布可以更好反应类似于算例所用,这一类呈现双峰分布特性的数据的分布情况。为了更进一步的从数学上证明上述结论,基于Eq.25所给的经验函数(Empirical Distribution Function,EDF),采用Anderson–Darling(AD)法对BioWeibull分布和Weibull分布的参数估计结果进行拟合优度检验。
式中,φ(Ni)是概率分布函数在样本Ni处的取值。当检验数据是否服从Weibull分布的时候,φ(N)为Eq.2中W(N),当检验数据是否服从BioWeibull分布时,φ(N)为Eq.4中F(N)。
AD法是一种基于Bootstrap方法的拟合优度检验方法,其主要步骤如下:
(1)假设样本数据服从X分布,将所检验的样本数据和估计得到分布参数带入Eq.25,计算得到检测统计量An0。
(2)根据Bootstrap的方法,对样本进行再抽样,并估计X分布的分布参数。
(3)根据(2)中估计得到的分布参数,基于原始样本数据计算检测统计量An。
(4)重复(2)(3)步骤5000次,获得5000个检测统计量An的序列。
(5)将序列从小至大排序,获得有序序列其中k表示排序后的序号,则给定置信度c下的分位点为检验临界值,本文中c=95%。若(1)中检测统计量An0大于Aα,则假设分布X不成立,反之则接受分布X为该寿命分布的模型。
根据上述方法对三种应力水平下疲劳寿命数据服从BioWeibull分布和Weibull的情况进行拟合优度检验,检验结果和中间边练列入表4。拟合优度检验的结果表明,三种应力水平下疲劳寿命数据均服从双峰威布尔分布,而双峰威布尔分布对疲劳寿命数据的拟合效果不佳。
表4拟合优度检验
确定疲劳寿命N95/95
根据上述方法计算得到三个应力水平下,两种分布的N95/95,见表5。从表中可以看出,在三种应力水平下,采用BioWeibull分布估计得到的N95/95显著大于Weibull分布估计得到的N95/95。
表5 N95/95 calculation results
为了进一步对比较两种分布模型的异同,计算了三个应力水平,两种分布几种典型置信度和可靠度下的寿命,见表6与表7。将表6中数据绘制成图4(a)-图4(c),从图中可以看出,在相同的置信度要求下(95%),采用双峰威布尔分布估计得到的高可靠度下(>=90%)的疲劳寿命高于威布尔分布的估计结果。图5(a)-图5(c)为根据表7绘制的相同可靠度要求下(95%)下,不同置信度要求的疲劳寿命估计结果数据对比图,从图中可以看出,在相同的可靠度要求下,采用双峰威布尔分布估计得到的不同置信度要求下的疲劳寿命均高于威布尔分布的估计结果。
表6 Fatigue life under a confidence level of 0.95and severalreliability levels
表7 Fatigue life under a reliability level of 0.95and severalconfidence levels
另外,本领域内的技术人员应当理解的是,在本发明实施例的申请文件中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括要素的过程、方法、物品或者设备中还存在另外的相同要素。
本发明实施例的说明书中,说明了大量具体细节。然而应当理解的是,本发明实施例的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。类似地,应当理解,为了精简本发明实施例公开并帮助理解各个发明方面中的一个或多个,在上面对本发明实施例的示例性实施例的描述中。
然而,并不应将该公开的方法解释成反映如下意图:即所要求保护的本发明实施例要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说,如权利要求书所反映的那样,发明方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本发明实施例的单独实施例。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种增材钛合金的疲劳寿命的确定方法,其特征在于,
以预设应力比,在不同应力水平下获取增材钛合金的疲劳寿命数据;
设定所述疲劳寿命数据服从双峰威布尔分布,构建概率密度函数f(N)和分布函数F(N);
基于ECM估值法,在求解非线性方程组时采用逐步最大化,得到双峰威布尔分布的显式解耦迭代式,进而求得π,α1,β1,α2,β2,其中,π,α1,β1,α2,β2为概率密度函数和分布函数中的分布参数,其中α1,α2为形状参数,β1,β2为范围参数,π为权重参数,其取值范围为0≤π≤1;
6.根据权利要求1-5任一项所述的增材钛合金的疲劳寿命的确定方法,其特征在于,在所述获得C=95%,P=95%时对应的NP,C之后,所述确定方法还包括:对双峰威布尔分布的参数估计结果进行拟合优化检验。
7.根据权利要求6所述的增材钛合金的疲劳寿命的确定方法,其特征在于,所述拟合优化检验包括:
(2)根据Bootstrap的方法,对所述增材钛合金的疲劳寿命数据进行随机又放回得抽样,并估计X分布的分布参数;
(3)根据(2)中估计得到的分布参数,基于原始抽样数据计算检测统计量An;
(4)重复(2)(3)步骤5000次,获得5000个检测统计量An的序列;
8.一种执行权利要求1-7任一项所述增材钛合金的疲劳寿命的确定系统,其特征在于,包括:
获取模块,获取以预设应力比,在不同应力水平下增材钛合金的疲劳寿命数据;
处理模块,设定所述疲劳寿命数据服从双峰威布尔分布,构建概率密度函数f(N)和分布函数F(N);
基于EM估值法,在求解非线性方程组时采用逐步最大化,得到双峰威布尔分布的显式解耦迭代式,进而求得π,α1,β1,α2,β2,其中,π,α1,β1,α2,β2为概率密度函数和分布函数中的分布参数,其中α1,α2为形状参数,β1,β2为范围参数,π为权重参数,其取值范围为0≤π≤1;
9.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机介质,其特征在于,所述处理器执行计算机程序时实现权利要求1~7中任意一项所述方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实施权利要求1~7中任意一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011455131.0A CN112949020B (zh) | 2020-12-10 | 2020-12-10 | 增材钛合金的疲劳寿命的确定方法、系统、设备与介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011455131.0A CN112949020B (zh) | 2020-12-10 | 2020-12-10 | 增材钛合金的疲劳寿命的确定方法、系统、设备与介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112949020A true CN112949020A (zh) | 2021-06-11 |
CN112949020B CN112949020B (zh) | 2022-08-12 |
Family
ID=76234822
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011455131.0A Active CN112949020B (zh) | 2020-12-10 | 2020-12-10 | 增材钛合金的疲劳寿命的确定方法、系统、设备与介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112949020B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113536490A (zh) * | 2021-08-27 | 2021-10-22 | 西安石油大学 | 一种抽油杆疲劳寿命预测方法及系统 |
CN115186393A (zh) * | 2022-08-16 | 2022-10-14 | 中国航发沈阳发动机研究所 | 一种低样本量的航空发动机结构件c-p-s-n曲线获取方法 |
CN115310260A (zh) * | 2022-07-01 | 2022-11-08 | 南京航空航天大学 | 疲劳寿命分布模型建模方法、系统、装置、计算机可读介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060243055A1 (en) * | 2005-04-28 | 2006-11-02 | Sundermeyer Jeffry N | Systems and methods for determining fatigue life |
WO2011000888A1 (fr) * | 2009-07-01 | 2011-01-06 | European Aeronautic Defence And Space Company Eads France | Procédé d'estimation de la durée de vie d'un circuit électronique intégré décananométrique ( "deep-sub-micron" ) |
CN109213963A (zh) * | 2017-07-03 | 2019-01-15 | 北京航空航天大学 | 一种激光沉积成型钛合金疲劳寿命统计分析方法 |
CN109470549A (zh) * | 2018-09-07 | 2019-03-15 | 北京航空航天大学 | 增材制造材料p-s-n曲线表征方法及其应用 |
CN109829136A (zh) * | 2019-01-22 | 2019-05-31 | 中国人民解放军火箭军工程大学 | 一种带随机跳变的退化设备的剩余寿命预测方法及系统 |
-
2020
- 2020-12-10 CN CN202011455131.0A patent/CN112949020B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060243055A1 (en) * | 2005-04-28 | 2006-11-02 | Sundermeyer Jeffry N | Systems and methods for determining fatigue life |
WO2011000888A1 (fr) * | 2009-07-01 | 2011-01-06 | European Aeronautic Defence And Space Company Eads France | Procédé d'estimation de la durée de vie d'un circuit électronique intégré décananométrique ( "deep-sub-micron" ) |
CN109213963A (zh) * | 2017-07-03 | 2019-01-15 | 北京航空航天大学 | 一种激光沉积成型钛合金疲劳寿命统计分析方法 |
CN109470549A (zh) * | 2018-09-07 | 2019-03-15 | 北京航空航天大学 | 增材制造材料p-s-n曲线表征方法及其应用 |
CN109829136A (zh) * | 2019-01-22 | 2019-05-31 | 中国人民解放军火箭军工程大学 | 一种带随机跳变的退化设备的剩余寿命预测方法及系统 |
Non-Patent Citations (3)
Title |
---|
严永滔等: "具双峰分布特征的全程产量Weibull预测模型研究及应用", 《油气藏评价与开发》 * |
凌丹等: "基于威布尔分布的疲劳剩余寿命可靠性预测方法", 《机械设计》 * |
刘小宁: "基于ECM算法的混合改进威布尔分布的参数估计", 《温州大学学报(自然科学版)》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113536490A (zh) * | 2021-08-27 | 2021-10-22 | 西安石油大学 | 一种抽油杆疲劳寿命预测方法及系统 |
CN113536490B (zh) * | 2021-08-27 | 2022-05-17 | 西安石油大学 | 一种抽油杆疲劳寿命预测方法及系统 |
CN115310260A (zh) * | 2022-07-01 | 2022-11-08 | 南京航空航天大学 | 疲劳寿命分布模型建模方法、系统、装置、计算机可读介质 |
CN115310260B (zh) * | 2022-07-01 | 2023-10-20 | 南京航空航天大学 | 疲劳寿命分布模型建模方法、系统、装置、计算机可读介质 |
CN115186393A (zh) * | 2022-08-16 | 2022-10-14 | 中国航发沈阳发动机研究所 | 一种低样本量的航空发动机结构件c-p-s-n曲线获取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112949020B (zh) | 2022-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112949020B (zh) | 增材钛合金的疲劳寿命的确定方法、系统、设备与介质 | |
Xia et al. | Statistical damage identification of structures with frequency changes | |
CN111339488B (zh) | 基于克里金Kriging的边坡系统失效概率计算方法 | |
CN112464490B (zh) | 增材钛合金的dfr确定方法 | |
CN108132914B (zh) | 有限数量产品性能失效率的高置信度统计推断和确认方法 | |
Caers et al. | Identifying tails, bounds and end-points of random variables | |
WO2014060160A1 (en) | Method and system for probabilistic fatigue crack life estimation | |
CN112949092B (zh) | 一种多阶段定期抽检型产品贮存寿命评估的离散参数回归法 | |
CN113360983B (zh) | 一种边坡可靠度分析与风险评估方法 | |
CN107657121B (zh) | 基于腐蚀级别评定的飞机结构性能预测处理方法及系统 | |
CN112784414A (zh) | 一种多部件整机贮存寿命置信下限评估方法 | |
CN117993340A (zh) | 一种存储芯片的质量评估方法 | |
CN111339487B (zh) | 基于径向基函数rbf的边坡系统失效概率计算方法 | |
CN117332498A (zh) | 基于熵权灰色关联度的车辆轻量化仿真方法、装置及介质 | |
Ruiz et al. | Analysis of correlated multivariate degradation data in accelerated reliability growth | |
CN106771779A (zh) | 电阻器抗辐照能力无损筛选方法及装置 | |
CN112487561B (zh) | 增材钛合金dfr确定方法、确定系统、设备及介质 | |
Vuppala et al. | A new inverse explicit flow curve determination method for compression tests | |
CN115329490A (zh) | 一种航空发动机结构件静力性能最小值确定方法 | |
CN114565232A (zh) | 一种导弹发动机贮存可用度评估方法 | |
CN114492074A (zh) | 一种概率损伤容限评估分析方法 | |
CN111523253A (zh) | 一种基于载荷融合数据确定铁路继电器吸合时间阈值的方法 | |
CN118010953B (zh) | 用于金属结构制造的智能测试方法及系统 | |
Cianetti et al. | How to experimentally monitor the fatigue behaviour of vibrating mechanical systems | |
CN113297773B (zh) | 一种结构损伤识别的不确定性定量分析方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |