CN108470082A - 一种基于指数-逆幂律的航空液压泵加速寿命试验建模方法 - Google Patents
一种基于指数-逆幂律的航空液压泵加速寿命试验建模方法 Download PDFInfo
- Publication number
- CN108470082A CN108470082A CN201810094971.5A CN201810094971A CN108470082A CN 108470082 A CN108470082 A CN 108470082A CN 201810094971 A CN201810094971 A CN 201810094971A CN 108470082 A CN108470082 A CN 108470082A
- Authority
- CN
- China
- Prior art keywords
- rotor
- hydraulic pump
- oil
- failure
- abrasive particles
- 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 48
- 238000000034 method Methods 0.000 title claims description 18
- 238000009826 distribution Methods 0.000 claims abstract description 22
- 239000002245 particle Substances 0.000 claims description 27
- 238000007789 sealing Methods 0.000 claims description 15
- 230000001186 cumulative effect Effects 0.000 claims description 10
- 230000007246 mechanism Effects 0.000 claims description 10
- 238000005299 abrasion Methods 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000001228 spectrum Methods 0.000 claims description 6
- 238000007476 Maximum Likelihood Methods 0.000 claims description 5
- 239000006061 abrasive grain Substances 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000003780 insertion Methods 0.000 claims description 3
- 230000037431 insertion Effects 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 239000007788 liquid Substances 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 abstract description 4
- 238000009825 accumulation Methods 0.000 abstract 1
- 239000003921 oil Substances 0.000 description 21
- 230000035882 stress Effects 0.000 description 9
- 239000010408 film Substances 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 229940030850 avar Drugs 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000005728 strengthening Methods 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006353 environmental stress Effects 0.000 description 1
- 239000010720 hydraulic oil Substances 0.000 description 1
- 230000001050 lubricating effect Effects 0.000 description 1
- 239000010687 lubricating oil Substances 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Rotary Pumps (AREA)
Abstract
本发明提出一种基于指数‑逆幂律的航空液压泵加速寿命试验模型,应用于航空液压泵可靠性领域。该模型通过定量计算揭示出转子配流盘单位时间内的磨损量与压力、转速不是简单的幂函数关系,在此基础上,提出了一种基于失效物理的指数‑逆幂律加速寿命模型。模型采用累积威布尔分布描述压液泵的寿命分布,采用指数‑逆幂律模型描述寿命与载荷的折算关系。仿真和实验结果表明:与传统基于统计经验的逆幂律加速寿命模型相比,本发明提出的指数‑逆幂律加速寿命模型能更深入地反应航空液压泵的失效物理特性,具有更高的估计精度。
Description
技术领域
本发明属于航空液压泵可靠性领域,具体涉及一种基于指数-逆幂律的航空液压泵加速寿命试验建模方法。
背景技术
轴向柱塞泵为飞机液压***提供高压油液以驱动飞控舵面和起落架等***。柱塞泵故障会导致飞行任务无法正常完成,甚至造成机毁人亡的严重事故,因此,对于飞机液压泵寿命和可靠性的要求越来越高。随着液压泵设计寿命的不断提高,对其进行寿命评估的难度和成本也越来越大。因此在工程上往往采用加速寿命试验的方法以压缩试验时间。在加速寿命试验中,通过对液压泵施加比正常工作条件更严酷的载荷,使液压泵更快发生故障,再通过统计方法估计出正常载荷下液压泵的寿命。
加速寿命试验是在不改变产品失效机理的前提下,通过加强应力的办法,加快产品故障、缩短试验时间,在较短的时间内预测出产品在正常应力作用下寿命特征的方法。不改变失效机理是加速寿命试验的前提,加强产品所承受的环境应力或工作应力是进行加速寿命试验的必要手段。加速寿命试验是通过加强应力来缩短试验时间,但如果应力过大,改变了产品的失效机理,则加速寿命试验就失去了意义。如果应力偏小,则会导致试验时间缩短并不明显,加速寿命试验无法得到最佳的效果。如何结合产品的实际工况,确定不改变产品失效机理、且能起到较好的加速作用的加速寿命试验剖面一直是困扰设计人员的难题。
飞机轴向柱塞泵的主要故障模式包括内部元部件的磨损、吸油口压力不足、转子轴承故障等。据统计,内部元部件的磨损是飞机轴向柱塞泵最主要的故障模式。液压泵的工作压力和转速是影响内部元部件磨损的重要因素,较高的转速会直接加快内部元部件的磨损,而较高的工作压力会使内部元部件的润滑效果变差,从而加剧其磨损。因此泵的转速和压力可以作为液压泵加速寿命试验的加速载荷。
理想情况下,液压泵的加速寿命试验应该基于液压泵的失效物理模型。然而,由于轴向柱塞泵失效机理过于复杂,目前还缺乏一种基于失效物理的加速寿命模型。液压泵的失效过程主要由内部元部件的润滑磨损状态决定,加速寿命模型不仅要考虑直接造成液压泵性能退化的因素,还要考虑多种耦合因素的影响。
目前可以检索到国外产品加速寿命试验的参考资料,但大多集中在统计方法的研究,关于加速模型的内容非常少。鉴于国外对我国相关技术采取封闭政策,我们对国外航空液压泵加速寿命试验模型如何建立无从得知,我国对航空液压泵加速寿命试验的研究也刚刚起步,到目前为止我国尚未有适合于工程应用的航空液压泵加速寿命试验模型。
发明内容
本发明基于对液压泵磨损失效过程的分析,提出了一种将失效物理和工程经验相结合的加速寿命模型。在所提出的加速寿命模型中,认为液压泵的基本寿命分布服从累积威布尔分布,而其寿命-载荷折算关系是一种指数-幂律关系。
本发明提出的一种基于指数-逆幂律的航空液压泵加速寿命试验建模方法,具体步骤为:
步骤一、液压泵故障机理及磨损过程分析;
步骤二、转子配流盘间距与泵的工作压力、转速的关系推导;
步骤三、基于指数-逆幂律的加速寿命模型;
步骤四、加速寿命模型误差分析。
本发明的优点和积极效果在于:
(1)在详细分析航空液压泵摩擦磨损失效机理的基础上,通过定量计算揭示出转子配流盘单位时间内的磨损量与压力、转速不是简单的幂函数关系,在此基础上,提出了一种基于失效物理的指数-逆幂律加速寿命模型。
(2)模型采用累积威布尔分布描述压液泵的寿命分布,采用指数-逆幂律模型描述寿命与载荷的折算关系。与传统基于统计经验的逆幂律加速寿命模型相比,本发明提出的指数-逆幂律加速寿命模型能更深入地反应航空液压泵的失效物理特性,具有更高的估计精度。
附图说明
图1是二维菱形磨粒模型;
图2是配流盘上一点处的磨损情况;
图3是磨粒嵌入转子和配流盘的深度;
图4是配流盘几何尺寸;
图5是加速寿命试验装置;
图6是加速寿命试验原理;
图7是额定工况下可靠性曲线和故障概率密度曲线;
图8是本发明的方法流程图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是一种基于指数-逆幂律的航空液压泵加速寿命试验建模方法,流程如图8所示,包括以下几个步骤:
步骤一、液压泵故障机理及磨损过程分析。
柱塞泵的内部有三对主要的摩擦副:转子和配流盘、柱塞和柱塞孔以及斜盘和滑靴。上述三对摩擦副的磨损会造成泄漏量的增加,导致泵的输出流量无法满足要求。理论计算和实验均表明,转子-配流盘磨损造成的泄漏在航空液压泵的泄漏中占主要地位。因此在工程上,通常用泄漏量来监测航空液压泵转子配流盘的磨损状况。
液压泵工作过程中,转子和配流盘的表面并不是直接接触,而是被一层很薄的润滑油膜隔开。然而,液压油并不是纯净的,其中悬浮着一定浓度的硬质颗粒,这些颗粒物会对转子和配流盘的表面造成磨损。将发生故障的配流盘放到电子显微镜下观察,可以在配流盘表面观察到一圈一圈的划痕,这说明配流盘的磨损主要是油液中的硬质颗粒物造成的。
研究表明,颗粒造成的磨损程度主要受颗粒物尺寸和两表面的距离影响。
根据Williams和Hyncica的二维菱形磨粒模型,如图1所示,当D/h≤1/cosβd时,磨粒在油液中自由转动,对转子和配流盘表面都不造成磨损;当D/h>1/cosβd时,磨粒在转子和配流盘两表面之间旋转,达到一个平衡位置,同时磨损两个表面。其中,D为磨粒最大直径,h为两表面间距,βd为磨粒形状角。
如图2所示,考虑配流盘上一点(x,y),油膜厚度为h(x,y),面积为dxdy,讨论(x,y)处油液中所含磨粒对配流盘的磨损情况。
如图3所示,根据几何关系以及转子和配流盘的相对硬度,单个磨粒嵌入转子和配流盘的深度分别为:
其中,ΔA为磨粒嵌入转子表面深度,ΔB为磨粒嵌入配流盘表面深度,H为配流盘和转子表面硬度比。
相应单个磨粒在转子和配流盘上造成的切槽横截面积为:
一定时间内单个磨粒对转子和配流盘造成的磨损体积分别为:
其中,f为磨损系数,ω为液压泵转速,r为点(x,y)对应的半径。
一定时间内(x,y)区域内所有磨粒对转子和配流盘造成的磨损体积分别为:
其中,k(D)为尺寸为D的磨粒对应的浓度。
所以一定时间内转子和配流盘磨损总体积为
其中,Ω表示转子和配流盘上所有需要考虑的区域。
对公式(5)进行整理可得,单位时间内磨粒造成的转子-配流盘接触面体积损失v,满足如下关系式:
这说明单位时间内磨粒造成的转子-配流盘接触面体积损失v与泵的转速ω成正比,与转子配流盘间距h的三次方成正比。
对于加速寿命模型来讲,最关心的是单位时间体积损失与泵的工作压力、转速的关系,所以下一步根据简化的雷诺方程,推导转子-配流盘间距h与泵的转速ω以及泵的工作压力Pd的解析表达式。
步骤二、转子配流盘间距与泵的工作压力、转速的关系推导。
转子配流盘的间距h可以通过柱坐标下简化的雷诺方程求得。在柱坐标下,雷诺方程可以被简化为:
其中,p为油膜动态压力分布,μ为油液动力黏度,ω为液压泵转速,r为径向坐标,θ为方位角坐标。
配流盘的几何尺寸如图4所示,在配流盘上半径为rr的圆周上间隔120°取三个参考点,则配流盘上任意一点(x,y)处的油膜厚度可以表示为这三个固定参考点的函数:
其中,h1 h2 h3分别三个参考点坐标。
转子配流盘的平均距离可以表示为:
柱坐标下,配流盘上任意一点(x,y)可以表示为
将公式(10)代入公式(8),可得
对公式(11)求关于θ的偏导数,可得
将公式(12)代入公式(7),并且等号两边对r进行积分,r为转化的柱坐标,可得:
其中,C1为待定常数,需要根据边界条件求解。
为了简化表达式,令
整理后得
对公式(15)等号两边对r进行积分,可得:
其中,C2为待定常数。
根据排油腔内封油带的边界条件,有
其中,PS为吸油腔压力,Pd为排油腔压力,r1为内密封带外半径,r2为内密封带外半径。
将公式(16)代入边界条件(17),则有
根据上面的方程可以解出待定系数C1和C2:
其中,
对内封油带上的压力分布进行积分,可得内封油带对转子的支撑力
其中,θ2-θ1为配油盘窗口夹角。
类似的,根据排油腔外封油带的边界条件,有
其中,r3为外密封带内半径,r4为外密封带外半径。
则外封油带对转子的支撑力为
其中,待定系数为
其中,
将转子压向配流盘的总压力可以近似计算为
根据转子稳定的条件,总压力应与总支撑力相等,即
Fp=Fsi+Fso (28)
将公式(21),(24),(27)代入公式(28),并对方程进行简化,仅保留以下变量:平均油膜厚度泵工作压力Pd,泵转速ω。其他参数视为常数合并处理,化简后得:
(C01Pd+C02ω+1)M(h0)=ln{N(h0)} (29)
其中,C01,C02是常数,M(h0),N(h0)是关于h0的有理多项式。公式(29)通常情况下无法用解析的方法求解,但是它反应出压力、转速和两表面间距满足某种指数关系,即
h=aPbωcexp{dP+eω} (30)
其中,a,b,c,d,e是待定参数,ω为转子转速,P为液压泵工作压力。
步骤三、基于指数-逆幂律的加速寿命模型。
传统的液压泵加速寿命试验一般采用逆幂律模型作为加速寿命模型(power-lawmodel,简称PL模型):
v=aPbωc (31)
其中,a,b,c是待定参数,ω为转子转速,P为液压泵工作压力。
结合公式(6)和(30)可知,转子-配流盘间距与压力转速不是简单的幂函数关系,故此,区别于传统的逆幂律模型,提出一种新的指数-逆幂律模型(exponential-power-lawmodel,简称EPL模型):
v=aPbωcexp{dP+eω} (32)
其中,a,b,c,d,e是待定参数,ω为转子转速,P为液压泵工作压力。
下面推导基于累积威布尔分布和指数-逆幂律的液压泵加速模型。一个完整的加速寿命模型应包括产品在一定载荷下的基本寿命分布,以及不同载荷间的寿命折算关系。
四个基本假设:
(1)液压泵的寿命在任一载荷S(P,ω)下服从如下Weibull分布:
其中,F(t)为故障概率,t为失效时间,η为特征寿命,m为形状参数。
(2)Weibull分布的形状参数m在任一载荷S(P,ω)下均为常数。
(3)特征寿命η在不同载荷Si(Pi,ωi)和Sj(Pj,ωj)之间可按如下关系折算:
(4)产品剩余寿命只与当前的累积损伤和当前的载荷有关,与累计损伤方式无关。
选择n个样本在变应力载荷谱下进行加速寿命试验,第j个样本经历的载荷历程,其相应的累积故障时间为 其中有n1个产品故障,有n2个产品进行到试验截止仍未故障。
故障样本Zj的累积故障密度为:
对于截尾样本Yj的累积可靠度为:
其中为历史载荷到当前载荷的等效折算时间,τj,0=0,
混合威布尔分布的极大似然函数为:
其中,n1为故障样本个数;n2为截尾样本个数,n=n1+n2。待估参数为
额定载荷S0(P0,ω0)下的可靠性函数为:
额定载荷下的故障概率密度函数为:
平均故障时间(MTTF)为
步骤四、加速寿命模型误差分析;
通常使用Fisher信息矩阵来评价极大似然估计的误差,定义Fisher信息矩阵为:
Fisher信息矩阵中的每一项,分别是对数极大似然函数对于各待估参数二阶偏微分的数学期望。Fisher信息矩阵中关于主对角线对称的项是相等的,所以一共需要计算6+5+4+3+2+1=21个二阶偏微分的数学期望。
借助Fisher信息矩阵,可以进一步计算估计量的误差方差。可靠性函数R0的估计误差方差为
其中
平均故障时间TMTTF的估计误差方差为
其中
本发明的一种基于指数-逆幂律的航空液压泵加速寿命试验模型,应用于航空液压泵可靠性领域。该模型通过定量计算揭示出转子配流盘单位时间内的磨损量与压力、转速不是简单的幂函数关系,在此基础上,提出了一种基于失效物理的指数-逆幂律加速寿命模型。模型采用累积威布尔分布描述压液泵的寿命分布,采用指数-逆幂律模型描述寿命与载荷的折算关系。仿真和实验结果表明:与传统基于统计经验的逆幂律加速寿命模型相比,本发明提出的指数-逆幂律加速寿命模型能更深入地反应航空液压泵的失效物理特性,具有更高的估计精度。
实施例
通过实际加速寿命试验数据验证本发明提出的加速寿命模型。实验对象为某型飞机的电马达驱动泵(Electric Motor Driven Pump,EMP)。加速寿命试验装置如图5所示,试验原理如图6所示。利用压力传感器监测泵的工作压力,利用流量传感器监测泵的壳体回油流量。选择压力和转速作为加速应力,泵的额定压力为21MPa,额定转速为4000r/min,通过负载节流阀调整泵的工作压力,通过调速电机控制泵的转速。泵的失效判据为回油流量超标。
根据摸底试验的结果,不改变失效机理的前提下,可以施加的极限应力为压力28MPa转速6000r/min。根据国军标的规定,飞机液压泵寿命试验应采用变载荷谱,加速寿命试验载荷谱如表1所示,按照表1所示的8个阶段进行循环,试验截尾时间为300小时。共有8个样本,试验的结果是有7台样本故障,1台试验样本截尾,各试验样本具体的失效时间和截尾时间如表2所示。
表1 加速寿命试验载荷谱
表2 加速寿命试验结果
根据表1所示的载荷谱和表2所示的实验结果,进行极大似然估计,参数估计结果如表3所示。
表3 加速寿命试验参数估计结果
根据表3的参数估计结果,计算出额定工况下(压力21MPa,转速为4000r/min)的故障概率密度曲线和可靠性曲线如图7所示。
EPL模型的Fisher信息矩阵为:
PL模型的Fisher信息矩阵为:
平均故障时间TMTTF及其估计误差方差AVar(TMTTF)如表4所示。
表4 平均故障时间及其估计误差方差
从表4可以看出,EPL模型对平均故障时间的估计误差方差AVar(TMTTF)小于PL模型,说明EPL模型具有更高的估计精度。
Claims (1)
1.一种基于指数-逆幂律的航空液压泵加速寿命试验建模方法,包括以下几个步骤:
步骤一、液压泵故障机理及磨损过程分析;
当D/h≤1/cos βd时,磨粒在油液中自由转动,对转子和配流盘表面都不造成磨损;当D/h>1/cos βd时,磨粒在转子和配流盘两表面之间旋转,达到一个平衡位置,同时磨损两个表面;其中,D为磨粒最大直径,h为两表面间距,βd为磨粒形状角;
设配流盘上一点(x,y),油膜厚度为h(x,y),面积为dxdy,单个磨粒嵌入转子和配流盘的深度分别为:
其中,ΔA为磨粒嵌入转子表面深度,ΔB为磨粒嵌入配流盘表面深度,H为配流盘和转子表面硬度比;
相应单个磨粒在转子和配流盘上造成的切槽横截面积为:
一定时间内单个磨粒对转子和配流盘造成的磨损体积分别为:
其中,f为磨损系数,ω为液压泵转速,r为点(x,y)对应的半径;
一定时间内(x,y)区域内所有磨粒对转子和配流盘造成的磨损体积分别为:
其中,k(D)为尺寸为D的磨粒对应的浓度;
所以一定时间内转子和配流盘磨损总体积为
其中,Ω表示转子和配流盘上所有需要考虑的区域;
对公式(5)进行整理可得,单位时间内磨粒造成的转子-配流盘接触面体积损失v,满足如下关系式:
步骤二、转子配流盘间距与泵的工作压力、转速的关系推导;
转子配流盘的间距h通过柱坐标下简化的雷诺方程求得,在柱坐标下,雷诺方程简化为:
其中,p为油膜动态压力分布,μ为油液动力黏度,ω为液压泵转速,r为径向坐标,θ为方位角坐标;
在配流盘上半径为rr的圆周上间隔120°取三个参考点,则配流盘上任意一点(x,y)处的油膜厚度表示为这三个固定参考点的函数:
其中,h1 h2 h3分别三个参考点坐标;
转子配流盘的平均距离表示为:
柱坐标下,配流盘上任意一点(x,y)表示为
将公式(10)代入公式(8),可得
对公式(11)求关于θ的偏导数,可得
将公式(12)代入公式(7),并且等号两边对r进行积分,r为转化的柱坐标,可得:
其中,C1为待定常数;
令
整理后得
对公式(15)等号两边对r进行积分,可得:
其中,C2为待定常数;
根据排油腔内封油带的边界条件,有
其中,PS为吸油腔压力,Pd为排油腔压力,r1为内密封带外半径,r2为内密封带外半径;将公式(16)代入边界条件(17),则有
解出待定系数C1和C2:
其中,
对内封油带上的压力分布进行积分,可得内封油带对转子的支撑力
其中,θ2-θ1为配油盘窗口夹角;
根据排油腔外封油带的边界条件,有
其中,r3为外密封带内半径,r4为外密封带外半径;
则外封油带对转子的支撑力为
其中,待定系数为
其中,
将转子压向配流盘的总压力近似计算为
根据转子稳定的条件,总压力应与总支撑力相等,即
Fp=Fsi+Fso (28)
将公式(21),(24),(27)代入公式(28),并对方程进行简化,仅保留以下变量:平均油膜厚度泵工作压力Pd,泵转速ω;其他参数视为常数合并处理,化简后得:
(C01Pd+C02ω+1)M(h0)=ln{N(h0)} (29)
其中,C01,C02是常数,M(h0),N(h0)是关于h0的有理多项式;公式(29)反应出压力、转速和两表面间距满足:
h=aPbωcexp{dP+eω} (30)
其中,a,b,c,d,e是待定参数,ω为转子转速,P为液压泵工作压力;
步骤三、基于指数-逆幂律的加速寿命模型;
设指数-逆幂律模型为:
v=aPbωcexp{dP+eω} (32)
其中,a,b,c,d,e是待定参数,ω为转子转速,P为液压泵工作压力;
设四个基本假设为:
(1)液压泵的寿命在任一载荷S(P,ω)下服从如下Weibull分布:
其中,F(t)为故障概率,t为失效时间,η为特征寿命,m为形状参数;
(2)Weibull分布的形状参数m在任一载荷S(P,ω)下均为常数;
(3)特征寿命η在不同载荷Si(Pi,ωi)和Sj(Pj,ωj)之间按如下关系折算:
(4)产品剩余寿命只与当前的累积损伤和当前的载荷有关,与累计损伤方式无关;
选择n个样本在变应力载荷谱下进行加速寿命试验,第j个样本经历的载荷历程,其相应的累积故障时间为 其中有n1个产品故障,有n2个产品进行到试验截止仍未故障;
故障样本Zj的累积故障密度为:
对于截尾样本Yj的累积可靠度为:
其中为历史载荷到当前载荷的等效折算时间,τj,0=0,
混合威布尔分布的极大似然函数为:
其中,n1为故障样本个数;n2为截尾样本个数,n=n1+n2;待估参数为额定载荷S0(P0,ω0)下的可靠性函数为:
额定载荷下的故障概率密度函数为:
平均故障时间为:
步骤四、加速寿命模型误差分析;
设Fisher信息矩阵为:
可靠性函数R0的估计误差方差为
其中
平均故障时间TMTTF的估计误差方差为
其中
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810094971.5A CN108470082B (zh) | 2018-01-31 | 2018-01-31 | 一种基于指数-逆幂律的航空液压泵加速寿命试验建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810094971.5A CN108470082B (zh) | 2018-01-31 | 2018-01-31 | 一种基于指数-逆幂律的航空液压泵加速寿命试验建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108470082A true CN108470082A (zh) | 2018-08-31 |
CN108470082B CN108470082B (zh) | 2020-01-10 |
Family
ID=63266278
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810094971.5A Expired - Fee Related CN108470082B (zh) | 2018-01-31 | 2018-01-31 | 一种基于指数-逆幂律的航空液压泵加速寿命试验建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108470082B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109931255A (zh) * | 2019-04-02 | 2019-06-25 | 哈工新欧(岳阳)测控装备有限公司 | 基于泄漏测试与泵壳温度测试的柱塞泵磨损评估***及方法 |
CN109960893A (zh) * | 2019-04-09 | 2019-07-02 | 重庆科技学院 | 一种连续管钻井定向负载参数分布规律实验装置测试方法 |
CN110287507A (zh) * | 2019-03-20 | 2019-09-27 | 北京航空航天大学 | 一种应用于恒压变量液压柱塞泵疲劳寿命分析方法 |
CN112149251A (zh) * | 2020-09-23 | 2020-12-29 | 南京工业大学 | 一种风力发电机主轴轴承模型寿命试验相似准则建立方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101050758A (zh) * | 2007-05-11 | 2007-10-10 | 华中科技大学 | 基于齿轮传动的斜轴式海水柱塞泵 |
CN101581295A (zh) * | 2009-06-22 | 2009-11-18 | 北京航空航天大学 | 一种基于dsp的机载液压泵源故障诊断*** |
CN102629300A (zh) * | 2012-03-15 | 2012-08-08 | 北京航空航天大学 | 一种基于灰色预测模型的步进应力加速退化数据评估方法 |
US20160312552A1 (en) * | 2015-04-27 | 2016-10-27 | Baker Hughes Incorporated | Integrated modeling and monitoring of formation and well performance |
-
2018
- 2018-01-31 CN CN201810094971.5A patent/CN108470082B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101050758A (zh) * | 2007-05-11 | 2007-10-10 | 华中科技大学 | 基于齿轮传动的斜轴式海水柱塞泵 |
CN101581295A (zh) * | 2009-06-22 | 2009-11-18 | 北京航空航天大学 | 一种基于dsp的机载液压泵源故障诊断*** |
CN102629300A (zh) * | 2012-03-15 | 2012-08-08 | 北京航空航天大学 | 一种基于灰色预测模型的步进应力加速退化数据评估方法 |
US20160312552A1 (en) * | 2015-04-27 | 2016-10-27 | Baker Hughes Incorporated | Integrated modeling and monitoring of formation and well performance |
Non-Patent Citations (2)
Title |
---|
葛薇 等: "航空液压泵磨损状况预测", 《北京航空航天大学学报》 * |
黄伯超 等: "基于性能可靠性约束的智能液压泵节能优化", 《北京航空航天大学学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110287507A (zh) * | 2019-03-20 | 2019-09-27 | 北京航空航天大学 | 一种应用于恒压变量液压柱塞泵疲劳寿命分析方法 |
CN109931255A (zh) * | 2019-04-02 | 2019-06-25 | 哈工新欧(岳阳)测控装备有限公司 | 基于泄漏测试与泵壳温度测试的柱塞泵磨损评估***及方法 |
CN109931255B (zh) * | 2019-04-02 | 2023-10-20 | 哈工新欧(岳阳)测控装备有限公司 | 基于泄漏测试与泵壳温度测试的柱塞泵磨损评估***及方法 |
CN109960893A (zh) * | 2019-04-09 | 2019-07-02 | 重庆科技学院 | 一种连续管钻井定向负载参数分布规律实验装置测试方法 |
CN109960893B (zh) * | 2019-04-09 | 2022-04-08 | 重庆科技学院 | 一种连续管钻井定向负载参数分布规律实验装置测试方法 |
CN112149251A (zh) * | 2020-09-23 | 2020-12-29 | 南京工业大学 | 一种风力发电机主轴轴承模型寿命试验相似准则建立方法 |
CN112149251B (zh) * | 2020-09-23 | 2023-07-28 | 南京工业大学 | 一种风力发电机主轴轴承模型寿命试验相似准则建立方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108470082B (zh) | 2020-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108470082B (zh) | 一种基于指数-逆幂律的航空液压泵加速寿命试验建模方法 | |
Manring | Friction forces within the cylinder bores of swash-plate type axial-piston pumps and motors | |
AU2019201038B2 (en) | Hydraulic system for a vehicle and method of using the same | |
CN111766064A (zh) | 一种舰载主轴轴承冲击试验方法 | |
CN110243591B (zh) | 一种变速器润滑油加注量的确定方法 | |
CN110287546B (zh) | 一种高压内啮合齿轮泵多轴疲劳寿命预测方法 | |
CN114018802A (zh) | 一种柱塞泵配流副摩擦系数的测量***及计算方法 | |
Zhang et al. | Effect of piston-slipper assembly mass difference on the cylinder block tilt in a high-speed electro-hydrostatic actuator pump of aircraft | |
CN108133109B (zh) | 一种基于非均匀间隙油膜的滑靴副偏磨磨损预测方法 | |
CN117307474B (zh) | 配流副油膜特性测试试验台及实验*** | |
Chao et al. | Development of an analytical model to estimate the churning losses in high-speed axial piston pumps | |
Zheglova et al. | Numerical and analytical evaluation of service life of the details of axial piston hydraulic machines with complicated configuration under cyclic loading | |
Zhang et al. | Performance degradation analysis of aviation hydraulic piston pump based on mixed wear theory | |
Vodička et al. | Wear behaviour of vanes for a rotary vane expander with various graphite materials under dry sliding conditions | |
CN113833641B (zh) | 机载燃油泵退化试验方案设计与寿命预测方法 | |
CN116257948A (zh) | 一种带液工况下含轴承故障的轴向柱塞泵分析方法及*** | |
Liu et al. | A Mixed Lubrication Model for Predicting the Lubrication Performance Degradation Behavior of Slipper Pair in Early Wear Failure | |
CN115389308A (zh) | 一种轴向柱塞泵配流副摩擦磨损试验机及其工作方法 | |
Arkhipov et al. | Modeling of cyclic life for compressor rotor of gas turbine engine taking into account production deviations | |
CN108256254A (zh) | 一种基于样本拓展的非平衡数据多维度参数估计方法 | |
Mancò et al. | Miniature gerotor pump prototype for automotive applications | |
Itoh et al. | Study on the oil supply system for rotary compressors | |
Han et al. | An abrasive wear model between valve plate and cylinder block in axial piston pumps | |
Nelson et al. | Reliability models for mechanical equipment | |
Huang et al. | Study on wear characteristics of friction pair in high pressure and high-speed axial piston pump |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200110 |
|
CF01 | Termination of patent right due to non-payment of annual fee |