CN113221286A - 一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法 - Google Patents
一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法 Download PDFInfo
- Publication number
- CN113221286A CN113221286A CN202110581811.5A CN202110581811A CN113221286A CN 113221286 A CN113221286 A CN 113221286A CN 202110581811 A CN202110581811 A CN 202110581811A CN 113221286 A CN113221286 A CN 113221286A
- Authority
- CN
- China
- Prior art keywords
- stress
- distribution
- threaded connection
- formula
- vibration
- 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
- 238000012360 testing method Methods 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000011156 evaluation Methods 0.000 claims abstract description 24
- 230000007246 mechanism Effects 0.000 claims description 15
- 230000001133 acceleration Effects 0.000 claims description 11
- 238000001276 Kolmogorov–Smirnov test Methods 0.000 claims description 10
- 238000005315 distribution function Methods 0.000 claims description 9
- 230000001419 dependent effect Effects 0.000 claims description 5
- 238000007689 inspection Methods 0.000 claims description 5
- 238000012886 linear function Methods 0.000 claims description 5
- 230000001186 cumulative effect Effects 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000004913 activation Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 claims description 2
- 230000007613 environmental effect Effects 0.000 claims description 2
- 230000008676 import Effects 0.000 claims description 2
- 238000011551 log transformation method Methods 0.000 claims description 2
- 239000002184 metal Substances 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000003746 surface roughness Effects 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提供一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法,其实施步骤如下:一:收集螺纹连接副加速寿命试验数据;二:螺纹连接副寿命分布参数估计;三:建立螺纹连接副的逆幂律-威布尔模型,外推其正常振动应力下螺纹连接副的寿命分布;四:k‑out‑of‑n模型下盘轴螺纹连接器振动可靠性评估。通过以上步骤,本发明给出了盘轴螺纹连接器的振动可靠性评估方法,解决了盘轴螺纹连接器高可靠长寿命,短时间内难以获得足够寿命数据从而无法进行可靠性评估的问题,保证了模型的有效性,计算简单,为人们提供了有效可行技术手段,具有重要的应用价值。
Description
技术领域
本发明提出一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法,它是一种盘轴螺纹连接器振动可靠性评估方法,从盘轴螺纹连接器的组成单元盘轴螺纹连接副入手,首先,基于螺纹连接副的振动加速寿命试验数据,综合运用最小二乘法和K-S检验对螺纹连接副进行多加速振动应力下的寿命分布参数估计;其次,以逆幂律-威布尔模型为基础,在保证失效机理一致(形状参数不变)的前提下,构建尺度参数(特征寿命)与振动应力的逆幂律方程;然后,使用最小二乘法,得到正常振动应力下螺纹连接副寿命分布尺度参数的估计值,从而进行螺纹连接副振动可靠性评估;最后,以盘轴螺纹连接器组成结构为依托,结合k-out-of-n模型,进行盘轴螺纹连接器的振动可靠性评估。所提出的方法可以有效解决正常振动应力下寿命试验数据获取时间长、成本高,即短时间内难以获得足够的寿命数据来进行盘轴螺纹连接器可靠性评估的问题。本专利属于装备剩余寿命预测与可靠性评估等相关技术领域。
背景技术
在航空发动机多级涡轮中,螺纹连接是其较为常见的连接方式,它是保证航空发动机结构完整的关键器件,而螺纹连接一旦失效,发动机乃至飞机的安全性能都会受到严重危害,因此,盘轴螺纹连接器的可靠性非常重要。
在航空发动机盘轴螺纹连接器的全寿命周期中,其所处复杂环境中的多种应力都将影响其可靠性,其中振动应力是对连接器的可靠性影响最大的因素之一。具体来说,螺栓轴向力会在振动应力作用下不断降低。在初期阶段,可能出现被连接件的密封性能降低、结构异响等故障;在后期阶段,螺栓可能发生断裂或脱落,造成重大安全事故。因此,对螺栓连接结构的松动进行分析及预测,具有重要意义。然而,螺栓松动机理非常复杂,其通常为一个非线性过程,同时伴随着接触非线性与材料非线性,在不同工况下,松动机理也有所差别。另一方面,宏观的松动可通过传感器等进行测量,但是试验手段难以监测接触面之间微观的滑动。此外,由于加工精度和表面粗糙度不可能完全相同,导致试验结果具有一定分散性,难以量化螺栓的松动机理,因此,在实际的工程应用中,基于螺栓滑脱的寿命试验数据来分析其松动情况,更具可行性。
通常针对盘轴螺纹连接器振动环境下寿命试验数据,按如下方法进行分析:
记P为试验所用的样本数目,将随机振动下的样本寿命数据按时间先后顺序顺序记录,从而得到的数据:
ti,i=1,2,…,P. (1)
式中:t表示样本寿命,i表示样本序号,P表示试验样本数目;接着,根据式(1)中的试验数据,建立连接器寿命分布模型,进而开展盘轴螺纹连接器可靠性评估。
但盘轴螺纹连接器具有高可靠、长寿命的特征,若在正常振动应力下开展寿命试验,短时间内难以获得足够的寿命数据来进行盘轴螺纹连接器的可靠性评估。考虑到连接器是由多组相同的螺纹连接副组成,它们共同分担连接任务,且盘轴螺纹连接器的失效判据与螺纹连接副的失效紧密相关,本专利提出基于螺纹连接副振动加速寿命试验数据,综合应用最小二乘法、K-S检验(Kolmogorov-Smirnov test)、逆幂律-威布尔模型,先进行螺纹连接副的振动可靠性评估,再基于盘轴螺纹连接器的失效判据,采用k-out-of-n模型(结构为k-out-of-n的系统指系统n个组件中的至少k个组件正常工作,系统正常工作),进行其振动可靠性评估。该专利有助于缩短寿命试验数据收集时间,准确快速地进行盘轴螺纹连接器振动可靠性评估。
发明内容
(1)本发明的目的:本发明针对盘轴螺纹连接器振动寿命试验数据获取时间长且成本高的问题,提供一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法;首先,它基于螺纹连接副的振动加速寿命试验数据,综合运用最小二乘法、K-S检验,得到各振动应力下的螺纹连接副寿命分布;其次,在失效机理一致的前提下,建立逆幂律-威布尔模型,构建特征寿命与振动应力的逆幂律方程;然后,使用最小二乘法外推得到正常振动应力下螺纹连接副的寿命分布;最后,结合k-out-of-n模型,进行盘轴螺纹连接器的振动可靠性评估。
(2)技术方案:基于上述理论和思路,本发明提供一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法,具体的实施步骤如下:
步骤一:收集螺纹连接副加速寿命试验数据
收集不同加速应力水平(振动应力)下,所有样本的失效时间;假设共有r个应力水平,记为{Sl;l=1,2,…,r},每个应力水平下有nl个样本,记录不同应力水平Sl下各样本的寿命数据,并将各个应力水平所对应的样本寿命数据按时间先后顺序进行排列,所得的试验数据集记为
{t(i,l);i=1,2,…,nl;l=1,2,…,r}. (2)
式中:t表示样本寿命,i表示样本序号,l表示应力序号,nl表示应力下样本数,r表示应力水平总数;
步骤二:螺纹连接副寿命分布参数估计
首先,确定螺纹连接副的寿命分布类型;双参数威布尔分布是可靠性工程一种常见的寿命分布,广泛应用于各种寿命试验的数据处理,假设螺纹连接副的寿命服从双参数威布尔分布,则寿命t的概率密度函数表示为
式中,t表示寿命,η表示特征寿命(尺度参数),m表示形状参数;
相应的寿命分布函数为
式中:t表示寿命,η表示特征寿命(尺度参数),m表示形状参数;将该公式进行双对数变换,得到如下线性函数:
式中:t表示寿命,η表示特征寿命(尺度参数),m表示形状参数;
其次,采用最小二乘法估计各振动应力下螺纹连接副的寿命分布参数;根据从小到大排序的寿命数据集{t(i,l);i=1,2,…,nl;l=1,2,…,r},计算中位秩,公式如下:
y(i,l)=Alx(i,l)+Bl. (7)
式中:x(i,l)表示自变量,y(i,l)表示因变量,Al、Bl表示常数,i表示样本序号,l表示应力序号;
然后,使用最小二乘法估计参数Al和Bl,具体公式为
同时,用相关系数法对所得寿命分布进行线性相关性检验:
式中:Dl表示K-S检验统计量,l表示应力序号,nl表示应力下样本数;
若Dl<Dc,则接受威布尔分布模型假设,否则,拒绝该假设;其中,Dc是置信水平为α和样本量为nl下的临界值,可通过查阅D值表获得;
步骤三:建立螺纹连接副的逆幂律-威布尔模型,外推其正常振动应力下螺纹连接副的寿命分布
通过步骤二,可得到r个振动应力水平下的寿命分布函数Fl(t)、特征寿命ηl和形状参数ml,l=1,2,…,r;为得到正常振动应力下螺纹连接副寿命分布,需要在失效机理不变的基础上(具体表示为形状参数m保持不变),建立特征寿命η与振动应力S之间的联系,即逆幂律-威布尔模型;模型具体体现为:产品寿命在统计上相互独立;产品寿命服从双参数威布尔分布;在不同振动应力作用下,形状参数保持不变;产品特征寿命与振动应力之间满足逆幂律方程;
首先,对逆幂律-威布尔模型的合理性进行假设检验;若产品在不同振动应力下的失效机理没有改变,则代表形状参数m相同;因此,对数据集{ml;l=1,2,…,r}进行一致性检验,检验假设H0是否成立:
H0:m1=m2=…=mr. (14)
式中:m表示形状参数,r表示应力水平总数;
由于失效时间T服从双参数威布尔分布W(m,η),则其对数Z=lnT服从极值分布G(μ,σ),分布函数为
H'0:σ1=σ2=…=σr, (16)
式中:m表示形状参数,r表示应力水平总数;
通过极值分布参数的最好线性无偏估计得到威布尔分布参数的估计量:
式中:D(nl,dl,i)和C(nl,dl,i)为无偏系数,该系数不仅与样本大小nl、失效次数dl有关,而且与故障顺序i有关,其值可在可靠性测试表中查找;
由σl的最好线性无偏估计的性质,可知:
式中:C(nl,dl,i)为无偏系数,与样本大小nl、失效次数dl、故障顺序i有关,μl表示极值分布的均值,σl表示极值分布的标准差,t(i,l)表示寿命,l表示应力序号,r表示应力总数;
式中:C(nl,dl,i)为无偏系数,与样本大小nl、失效次数dl、故障顺序i有关,yli表示标准极值分布的顺序统计量,σl表示极值分布的标准差,l表示应力序号,r表示应力总数;
取
作为检验H'0的统计量;
式中:n表示样本大小,d表示失效次数,r表示应力总数,j表示应力序号;
当H'0:σ1=σ2=…=σr成立时,
式中:C(nl,dl,i)为无偏系数,与样本大小nl、失效次数dl、故障顺序i有关,yli表示标准极值分布的顺序统计量,σl表示极值分布的标准差,l表示应力序号,r表示应力总数,j表示应力序号,表示nj个标准极值分布顺序统计量的线性组合;
求出显著性水平α下的临界值Cα,使得:
其次,建立螺纹连接副基于振动应力影响下的逆幂律-威布尔模型;按照温度应力加速模型(阿伦尼斯模型)与振动应力加速模型(逆幂律模型),螺纹连接副寿命t与振动应力S、环境温度T的关系表示为
t=AS-αeΔE/hT, (24)
式中:t表示螺纹连接副寿命,S表示为振动应力,T表示为环境温度,ΔE代表激活能,h代表玻尔兹曼常数,而A和α表示待定系数;在螺纹连接副的加速寿命试验中,环境温度T保持不变,可简化为
t=KS-α, (25)
式中,t表示螺纹连接副寿命,S表示为振动应力,K和α表示待定系数;因为η为特征寿命,η∝t,则
η=βS-α. (26)
进行对数化线性处理得到:
lnη=-αlnS+lnβ, (27)
式中,η表示特征寿命,S表示为振动应力,α和β表示待定系数;导入数据集{(ηl,Sl);l=1,2,…,r},并转化为{(lnSl,lnηl);l=1,2,…,r};同样采用最小二乘法,使用式(8)、式(9)计算估计值与进而求出模型参数估计值
设正常条件下的振动应力为S0,则威布尔分布参数为
最后,得到螺纹连接副正常振动应力下失效时间t的概率密度函数:
步骤四:k-out-of-n模型下盘轴螺纹连接器振动可靠性评估
根据螺纹连接副正常振动应力下寿命的概率密度函数,推出其可靠度函数:
式中,t表示寿命,表示特征寿命(尺度参数)估计值,表示形状参数估计值;而盘轴螺栓连接器由N组螺纹连接副组成,由于螺纹连接副共同承担连接作用,因此连接器系统的可靠性与螺纹连接副的可靠性间的关系可看成一个k-out-of-n模型;假设某个阈值k,若连接器失效判据为:至少k个螺纹连接副失效;则盘轴螺纹连接器的可靠度函数为RS(t):
通过以上步骤,以螺纹连接副为元素,以加速寿命试验数据为基础,以k-out-of-n模型为桥梁,给出了盘轴螺纹连接器的振动可靠性评估方法,解决了盘轴螺纹连接器高可靠长寿命,短时间内难以获得足够寿命数据从而无法进行可靠性评估的问题,结合检验方法,保证了模型的有效性,符合工程实际情况,计算简单,为工程技术人员提供了有效可行技术手段,存在重要的应用价值。
(3)优点:
本发明提出了一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法,其优点如下:
①本发明提出了从盘轴螺纹连接副出发,以加速振动寿命试验数据为基础,以k-out-of-n模型为架构,进行盘轴螺纹连接器振动可靠性评估,解决了盘轴螺纹连接器高可靠长寿命,短时间内难以获得足够寿命数据从而无法进行可靠性评估的问题。
②本发明提出的方法计算简便,容易实现,且更加符合工程实际,方便工程技术人员掌握使用,方法科学,便于应用推广。
附图说明
图1是本发明所述方法的流程图。
具体实施方式
本发明一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法,其流程图如图1所示。
下面以某型航空发动机连接涡轮盘与压气机轴的4组螺纹连接副组成的连接器为例,对本发明做进一步详细说明;
在盘轴螺纹连接器加速寿命试验中,加速因子为振动应力,试验的应力水平所对应的功率谱密度分别设置为0.2g2/Hz,0.4g2/Hz,0.6g2/Hz,0.8g2/Hz和1.0g2/Hz,在每个振动应力下,分别安排一定量的样本(螺纹连接副)进行随机振动,开展加速寿命试验。试验过程中,当螺纹连接副出现滑脱现象,即判定为失效。每出现一次失效,便记录下其失效时间,随后按时间大小排序。
本发明提出一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法,其具体的实施步骤如下:
步骤一:收集螺纹连接副加速寿命试验数据
分别设置0.2g2/Hz,0.4g2/Hz,0.6g2/Hz,0.8g2/Hz和1.0g2/Hz的振动应力水平,通过加速寿命试验试验获得不同应力下排序后的失效时间数据,如表1所示。
表1盘轴螺纹连接器中螺纹连接副的加速寿命试验数据
步骤二:螺纹连接副寿命分布参数估计
首先,根据从小到大排序的寿命数据集{t(i,l);i=1,2,3,4,5;l=1,2,3,4,5},使用式(6)计算对应失效时间下的中位秩:
式中i代表累积失效数,l代表应力序号;通过计算得出中位秩,如表2所示;
表2螺纹连接副加速寿命试验中位秩表
表3各振动应力下参数Al和Bl的估计值
用相关系数法(研究变量之间线性相关程度的方法)对求得的威布尔分布进行线性相关性检验,使用式(10)计算相关系数可以得出线性相关性显著的结论;结果如表4所示,当取显著性水平α=0.05时,不同振动应力下自变量与因变量线性相关;
表4线性相关性检验结果
然后,使用K-S检验验证模型合理性,使用式(11)得出K-S检验统计量Dl,若Dl<Dc,则接受威布尔分布模型(螺纹连接副寿命服从威布尔分布);其中,Dc为临界值,可根据置信水平α和样本量nl,查阅D值表获得;结果如表5所示,当取α=0.05,nl=5时,K-S检验通过,接受威布尔分布模型为适用的模型;
表5 K-S检验结果
表6威布尔分布参数估计值
步骤三:建立螺纹连接副的逆幂律-威布尔模型与外推其正常应力下寿命分布
从步骤二中可看到,根据表1加速寿命试验数据计算出了5个振动应力水平下失效时间的累积分布函数Fl(t)、特征寿命和形状参数的估计值和接下来,需要在保持失效机理不变的基础上,建立起特征寿命η与振动应力S之间的联系;
首先,进行失效机理一致性检验;若产品在不同振动应力下的失效机理没有改变,则代表形状参数m相同;针对数据集{ml;l=1,2,,5},检验假设H0是否成立:
H0:m1=m2=…=m5←—→H0:m1、m2、…、m5不全相等, (36)
式中:m1、m2、…、m5分别表示不同振动应力下的形状参数;
作为检验H0的统计量,根据式(20)得出:
其次,建立逆幂律-威布尔模型;按照温度应力加速模型(阿伦尼斯模型)与振动应力加速模型(逆幂律模型),建立式(24)中螺纹连接副特征寿命与振动应力的逆幂律方程与式(25)中的线性方程;同时根据表6中各振动应力的特征寿命估计值,使用式(8)、式(9)进行最小二乘估计,求出模型参数估计值以及模型公式:
η=14.6875*S-1.044. (38)
设正常条件下的振动应力为S0=0.001g2/Hz,则威布尔分布参数为
最后,得到螺纹连接副正常振动应力下失效时间t的概率密度函数:
步骤四:k-out-of-n模型下盘轴螺纹连接器振动可靠性评估
根据螺纹连接副正常环境下失效时间的概率密度函数,推出其可靠度函数:
因为盘轴螺栓连接器由4组螺纹连接副组成,可建立2-out-of-4:F模型(4部件的模型失效当且仅当有2个及以上的部件失效),得到盘轴螺纹连接器的可靠度函数为
综上,本发明以某型航空发动机连接涡轮盘与压气机轴的4组螺纹连接副组成的连接器为研究对象,针对盘轴螺纹连接器受到振动应力影响失效,而短时间正常应力下的失效数据难以获得所带来的可靠性评估问题,首先,基于螺纹连接副的振动加速寿命试验数据,综合运用最小二乘估计方法、K-S检验对螺纹连接副进行寿命分布参数估计;其次,建立逆幂律-威布尔模型,在保证失效机理一致的前提下构建特征寿命与振动应力的逆幂律方程;然后,使用最小二乘估计方法得到正常振动应力螺纹连接副的威布尔参数估计值以及失效时间概率密度函数;最后,建立k-out-of-n模型,计算出盘轴螺纹连接器的可靠度函数,有效完成了盘轴螺纹连接器的振动可靠性评估。
Claims (1)
1.一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法,其特征在于:具体的实施步骤如下:
步骤一:收集螺纹连接副加速寿命试验数据
收集不同加速应力水平下,所有样本的失效时间;设共有r个应力水平,记为
{Sl;l=1,2,…,r},每个应力水平下有nl个样本,记录不同应力水平Sl下各样本的寿命数据,并将各个应力水平所对应的样本寿命数据按时间先后顺序进行排列,所得的试验数据集记为
{t(i,l);i=1,2,…,nl;l=1,2,…,r} (2)
式中:t表示样本寿命,i表示样本序号,l表示应力序号,nl表示应力下样本数,r表示应力水平总数;
步骤二:螺纹连接副寿命分布参数估计
首先,确定螺纹连接副的寿命分布类型;设螺纹连接副的寿命服从双参数威布尔分布,则寿命t的概率密度函数表示为
式中,t表示寿命,η表示特征寿命,m表示形状参数;
相应的寿命分布函数为
式中:t表示寿命,η表示特征寿命,m表示形状参数;将该公式进行双对数变换,得到如下线性函数:
式中:t表示寿命,η表示特征寿命,m表示形状参数;
其次,采用最小二乘法估计各振动应力下螺纹连接副的寿命分布参数;根据从小到大排序的寿命数据集{t(i,l);i=1,2,…,nl;l=1,2,…,r},计算中位秩,公式如下:
式中:t(i,l)表示寿命,i表示样本序号,l表示应力序号,nl表示应力下样本数;
y(i,l)=Alx(i,l)+Bl (7)
式中:x(i,l)表示自变量,y(i,l)表示因变量,Al、Bl表示常数,i表示样本序号,l表示应力序号;
然后,使用最小二乘法估计参数Al和Bl,具体公式为
同时,用相关系数法对所得寿命分布进行线性相关性检验:
式中:Dl表示K-S检验统计量,l表示应力序号,nl表示应力下样本数;
若Dl<Dc,则接受威布尔分布模型假设,否则,拒绝该假设;其中,Dc是置信水平为α和样本量为nl下的临界值,通过查阅D值表获得;
步骤三:建立螺纹连接副的逆幂律-威布尔模型,外推其正常振动应力下螺纹连接副的寿命分布
通过步骤二,得到r个振动应力水平下的寿命分布函数Fl(t)、特征寿命ηl和形状参数ml,l=1,2,…,r;为得到正常振动应力下螺纹连接副寿命分布,需要在失效机理不变的基础上,建立特征寿命η与振动应力S之间的联系,即逆幂律-威布尔模型;模型具体体现为:产品寿命在统计上相互独立;产品寿命服从双参数威布尔分布;在不同振动应力作用下,形状参数保持不变;产品特征寿命与振动应力之间满足逆幂律方程;
首先,对逆幂律-威布尔模型的合理性进行假设检验;若产品在不同振动应力下的失效机理没有改变,则代表形状参数m相同;因此,对数据集{ml;l=1,2,…,r}进行一致性检验,检验假设H0是否成立:
H0:m1=m2=…=mr (14)
式中:m表示形状参数,r表示应力水平总数;
由于失效时间T服从双参数威布尔分布W(m,η),则其对数Z=lnT服从极值分布G(μ,σ),分布函数为
H'0:σ1=σ2=…=σr, (16)
式中:m表示形状参数,r表示应力水平总数;
通过极值分布参数的最好线性无偏估计得到威布尔分布参数的估计量:
式中:D(nl,dl,i)和C(nl,dl,i)为无偏系数,该系数不仅与样本大小nl、失效次数dl有关,而且与故障顺序i有关,其值可在可靠性测试表中查找;
由σl的最好线性无偏估计的性质,可知:
式中:C(nl,dl,i)为无偏系数,与样本大小nl、失效次数dl、故障顺序i有关,μl表示极值分布的均值,σl表示极值分布的标准差,t(i,l)表示寿命,l表示应力序号,r表示应力总数;
式中:C(nl,dl,i)为无偏系数,与样本大小nl、失效次数dl、故障顺序i有关,yli表示标准极值分布的顺序统计量,σl表示极值分布的标准差,l表示应力序号,r表示应力总数;
取
作为检验H'0的统计量;
式中:n表示样本大小,d表示失效次数,r表示应力总数,j表示应力序号;
当H'0:σ1=σ2=…=σr成立时,
式中:C(nl,dl,i)为无偏系数,与样本大小nl、失效次数dl、故障顺序i有关,yli表示标准极值分布的顺序统计量,σl表示极值分布的标准差,l表示应力序号,r表示应力总数,j表示应力序号,表示nj个标准极值分布顺序统计量的线性组合;
求出显著性水平α下的临界值Cα,使得:
其次,建立螺纹连接副基于振动应力影响下的逆幂律-威布尔模型;按照温度应力加速模型与振动应力加速模型,螺纹连接副寿命t与振动应力S、环境温度T的关系表示为
t=AS-αeΔE/hT, (24)
式中:t表示螺纹连接副寿命,S表示为振动应力,T表示为环境温度,ΔE代表激活能,h代表玻尔兹曼常数,而A和α表示待定系数;在螺纹连接副的加速寿命试验中,环境温度T保持不变,可简化为
t=KS-α, (25)
式中,t表示螺纹连接副寿命,S表示为振动应力,K和α表示待定系数;因为η为特征寿命,η∝t,则
η=βS-α (26)
进行对数化线性处理得到:
lnη=-αlnS+lnβ, (27)
式中,η表示特征寿命,S表示为振动应力,α和β表示待定系数;导入数据集{(ηl,Sl);l=1,2,…,r},并转化为{(lnSl,lnηl);l=1,2,…,r};同样采用最小二乘法,使用式(8)、式(9)计算估计值与进而求出模型参数估计值
设正常条件下的振动应力为S0,则威布尔分布参数为
最后,得到螺纹连接副正常振动应力下失效时间t的概率密度函数:
步骤四:k-out-of-n模型下盘轴螺纹连接器振动可靠性评估
根据螺纹连接副正常振动应力下寿命的概率密度函数,推出其可靠度函数:
式中,t表示寿命,表示特征寿命估计值,表示形状参数估计值;而盘轴螺栓连接器由N组螺纹连接副组成,由于螺纹连接副共同承担连接作用,因此连接器系统的可靠性与螺纹连接副的可靠性间的关系可看成一个k-out-of-n模型;设某个阈值k,若连接器失效判据为:至少k个螺纹连接副失效;则盘轴螺纹连接器的可靠度函数为RS(t):
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110581811.5A CN113221286B (zh) | 2021-05-24 | 2021-05-24 | 一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110581811.5A CN113221286B (zh) | 2021-05-24 | 2021-05-24 | 一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113221286A true CN113221286A (zh) | 2021-08-06 |
CN113221286B CN113221286B (zh) | 2022-08-19 |
Family
ID=77099025
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110581811.5A Active CN113221286B (zh) | 2021-05-24 | 2021-05-24 | 一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113221286B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114970157A (zh) * | 2022-05-26 | 2022-08-30 | 北京航空航天大学 | 电子产品在电压应力作用下的小样本试验寿命预测方法 |
CN115906541A (zh) * | 2023-02-28 | 2023-04-04 | 航天精工股份有限公司 | 基于多竞争失效模式的紧固连接系统可靠性正向设计方法 |
CN118171753A (zh) * | 2024-05-13 | 2024-06-11 | 摩多利智能传动(江苏)有限公司 | 基于振动分析的电机故障智能诊断系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150127308A1 (en) * | 2013-11-06 | 2015-05-07 | Vetco Gray Inc. | Stress Amplification Factor Analysis Methodology for Assessing Fatigue Performance of Threaded Connectors |
CN108027901A (zh) * | 2015-08-13 | 2018-05-11 | 拜耳医药保健有限公司 | 确定其可替换部件的有效寿命的系统和指示其可替换部件的有效寿命的装置及其方法 |
CN112784415A (zh) * | 2021-01-22 | 2021-05-11 | 北京航空航天大学 | 一种定数截尾加速寿命试验机理等同性检验与寿命预测方法 |
-
2021
- 2021-05-24 CN CN202110581811.5A patent/CN113221286B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150127308A1 (en) * | 2013-11-06 | 2015-05-07 | Vetco Gray Inc. | Stress Amplification Factor Analysis Methodology for Assessing Fatigue Performance of Threaded Connectors |
CN108027901A (zh) * | 2015-08-13 | 2018-05-11 | 拜耳医药保健有限公司 | 确定其可替换部件的有效寿命的系统和指示其可替换部件的有效寿命的装置及其方法 |
CN112784415A (zh) * | 2021-01-22 | 2021-05-11 | 北京航空航天大学 | 一种定数截尾加速寿命试验机理等同性检验与寿命预测方法 |
Non-Patent Citations (1)
Title |
---|
巩浩: "《振动条件下螺纹预紧力衰退机理和影响因素研究》", 《机械工程学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114970157A (zh) * | 2022-05-26 | 2022-08-30 | 北京航空航天大学 | 电子产品在电压应力作用下的小样本试验寿命预测方法 |
CN114970157B (zh) * | 2022-05-26 | 2024-05-28 | 北京航空航天大学 | 电子产品在电压应力作用下的小样本试验寿命预测方法 |
CN115906541A (zh) * | 2023-02-28 | 2023-04-04 | 航天精工股份有限公司 | 基于多竞争失效模式的紧固连接系统可靠性正向设计方法 |
CN115906541B (zh) * | 2023-02-28 | 2023-05-05 | 航天精工股份有限公司 | 基于多竞争失效模式的紧固连接系统可靠性正向设计方法 |
CN118171753A (zh) * | 2024-05-13 | 2024-06-11 | 摩多利智能传动(江苏)有限公司 | 基于振动分析的电机故障智能诊断系统 |
Also Published As
Publication number | Publication date |
---|---|
CN113221286B (zh) | 2022-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113221286B (zh) | 一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法 | |
CN108959676B (zh) | 一种考虑有效冲击的退化建模与寿命预测方法 | |
CN102426307B (zh) | 一种基于短期寿命试验数据的快速寿命评价方法 | |
CN107563054B (zh) | 一种基于SWT参数的Weakest-Link方法的涡轮盘概率寿命分析方法 | |
CN112069697A (zh) | 多应力多失效模式相依竞争条件下加速退化试验统计分析方法 | |
CN107145641B (zh) | 叶片振动疲劳概率寿命预测方法 | |
CN112784413B (zh) | 一种zn-40阻尼减振结构剩余贮存寿命评估方法 | |
CN111310314B (zh) | 一种基于人工智能确定机电装置寿命的方法和系统 | |
CN112733332B (zh) | 一种考虑拧紧力矩的螺栓疲劳可靠性评估方法 | |
CN114295491A (zh) | 一种蠕变损伤及变形随时间演化行为的预测方法 | |
CN110414086B (zh) | 一种基于灵敏度的综合应力加速因子计算方法 | |
CN112926698B (zh) | 一种大型旋转装备振动预测与装配评价方法 | |
CN114492074A (zh) | 一种概率损伤容限评估分析方法 | |
CN111307483A (zh) | 一种机械产品长周期疲劳试验数据处理与趋势预判方法 | |
Letot et al. | Dynamic reliability degradation based models and maintenance optimization | |
CN115408755B (zh) | 一种考虑时变效应的组合梁桥动力疲劳可靠性评估方法 | |
CN110895624B (zh) | 基于最大熵谱估计的加速贮存与自然贮存退化数据一致性检验法 | |
CN109490334B (zh) | 一种运用残余应力预测模型的t字型锻件无损测试方法 | |
CN118260564B (zh) | 一种基于贝叶斯算法的零件寿命抽样检验方法及系统 | |
CN112001061A (zh) | 一种弹性非金属材料产品加速寿命预估方法 | |
Leitner et al. | Fatigue damage prediction as a part of technical systems reliability assessment | |
CN113902282B (zh) | 一种面向小样本数据条件的航班运行风险评价方法 | |
CN117262237B (zh) | 考虑装配分散性的飞机座舱盖骨架模拟件疲劳试验方法 | |
Marková et al. | Statistical evaluation of fatigue curve of vibration-loaded structures | |
Bacharoudis et al. | Measurement uncertainty of fatigue properties and its effect on the wind turbine blade reliability level |
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 |