飞机锁机构部件磨损退化与功能退化竞争失效分析方法与流程

本发明涉及飞机锁机构部件磨损退化与功能退化竞争失效分析方法，属于飞机锁机构部件可靠性分析
技术领域：
。
背景技术：
：飞机锁机构的可靠性水平直接影响飞机特定机构的工作状态，进而影响整机的飞行安全。由于磨损、腐蚀、老化等原因，工程机构在使用过程中通常存在着部件的磨损退化，即便是正常的工作环境，这些导致磨损的原因也无法避免。与此同时，机构在设计之初都有其特定的功能，机构功能是否可完成直接决定了机构的工作状态是否正常，而机构功能的完成又由机构部件共同决定，机构部件的磨损退化同时也会引发机构功能的退化。对于一个确定的机构而言，机构功能量是机构各个部件的函数，该函数由机构的传递关系和机构组成特性共同决定，且机构部件的磨损退化量和功能量的退化是相关的，如何正确建立退化模型并正确描述这种相关性是准确分析机构退化竞争失效的基础。现有的有关竞争退化的论文主要针对处于外界冲击失效和自身退化失效情形下的系统进行建模，外界冲击失效以遵循特定速率的泊松过程描述，退化过程则通常使用退化路径模型或随机过程来描述。对于相关性的随机变量，基本的方法是建立变量的联合概率分布函数，而由于建立联合概率分布函数存在的计算量的困难，更常使用的是利用copula函数建立描述两个随机变量之间的相关性，并据此建立联合概率分布函数以进行可靠性分析计算。目前针对飞机锁机构的失效分析方法中，对处于部件磨损退化和功能退化下的飞机锁机构进行可靠度分析的时候，需要对两种失效模式分别建模并考虑相关性。技术实现要素：本发明目的是为了解决现有飞机锁机构的失效分析方法中，需要对部件磨损退化和功能退化分别建模再考虑相关性，分析结果可靠性低的问题，提供了一种飞机锁机构部件磨损退化与功能退化竞争失效分析方法。本发明所述飞机锁机构部件磨损退化与功能退化竞争失效分析方法，它包括以下步骤：步骤一：根据飞机锁机构，定义功能量，再确定影响功能量实现的部件，并在部件中确定所有磨损退化分部件；步骤二：利用非线性漂移布朗运动对每个磨损退化分部件建立退化模型，确定不同时间点各个磨损退化分部件的随机特性；步骤三：引入代理模型，建立各个磨损退化分部件的随机特性与锁机构功能量之间的传递关系，进而得到锁机构功能量的随机分布特性；步骤四：通过fgmcopula函数建立所有磨损退化分部件的磨损退化量分布与锁机构功能量的随机分布之间的联合分布函数；步骤五：根据联合分布函数分别计算考虑竞争情况下的不同失效模式的竞争失效概率，并计算得到锁机构的整体可靠度。本发明的优点：本发明涉及飞机锁机构失效分析和可靠性计算，特别是锁机构处在部件磨损退化失效以及锁机构功能失效两种失效模式下，且两种失效模式是相互竞争关系。本发明用非线性漂移布朗运动对部件磨损退化失效进行建模，并用代理模型建立各个磨损退化分部件与功能量之间的函数传递关系，以此寻找功能量的退化规律。fgmcopula函数用来建立部件磨损退化量和功能退化量之间的联合概率密度函数，并基于此计算不同失效模式下的竞争失效概率和锁机构的整体可靠度。通过实验对比，本发明方法的计算结果更接近实际工程情况。与此同时，与蒙特卡洛计算所得到的结果相比，本方法所得到的计算结果具有足够的计算精度，可保证在使用过程的准确性。附图说明图1是本发明所述飞机锁机构部件磨损退化与功能退化竞争失效分析方法的流程图；图2是具体实施方式二中飞机舱门上位锁机构打开状态示意图；图3是飞机舱门上位锁机构锁闭状态示意图；图4是相关系数|α|＝0.8时的fgmcopula函数分布图；图5是铰链a磨损退化失效概率曲线图；图6是飞机舱门上位锁机构功能失效概率曲线图；图7是锁机构的整体可靠度曲线图。具体实施方式具体实施方式一：下面结合图1说明本实施方式，本实施方式所述飞机锁机构部件磨损退化与功能退化竞争失效分析方法，它包括以下步骤：步骤一：根据飞机锁机构，定义功能量，再确定影响功能量实现的部件，并在部件中确定所有磨损退化分部件；步骤二：利用非线性漂移布朗运动对每个磨损退化分部件建立退化模型，确定不同时间点各个磨损退化分部件的随机特性；步骤三：引入代理模型，建立各个磨损退化分部件的随机特性与锁机构功能量之间的传递关系，进而得到锁机构功能量的随机分布特性；步骤四：通过fgmcopula函数建立所有磨损退化分部件的磨损退化量分布与锁机构功能量的随机分布之间的联合分布函数；步骤五：根据联合分布函数分别计算考虑竞争情况下的不同失效模式的竞争失效概率，并计算得到锁机构的整体可靠度。具体实施方式二：下面结合图1至图7说明本实施方式，本实施方式对实施方式一作进一步说明，步骤一中，所述飞机锁机构为飞机舱门上位锁机构，飞机舱门上位锁机构由摇臂1、摇臂中连杆2、摇臂下连杆3、锁钩4和锁环5组成，定义锁钩4的钩末端顶端到锁环5中心的距离作为运动功能量，摇臂中连杆2的一端与摇臂1中连接点连接，确定摇臂中连杆2另一端与固定部件连接通过的铰链a为磨损退化分部件；图2和图3中对锁机构的工作过程进行了简单示意，飞机舱门上位锁机构的功能是确保锁钩4可以在两个摇臂连杆等构件的带动下锁住锁环5，以完成其锁定功能；但由于制造误差、装备误差、运动副间隙等原因，锁闭状态的时候锁钩不一定能够在正确的位置钩住锁环，使运动精度达不到进而影响锁机构功能的实现。根据锁的功能描述和运动原理，锁的功能是确保锁钩可以正确地钩住锁环。为了描述锁的功能，定义锁钩的顶端到锁环中心的距离为运动功能量，如图3局部放大图中虚线段所示。锁机构的功能由组成锁的诸多组件共同决定，因此这些组件的几何长度都会影响锁机构功能。在锁的实际使用过程中发现，铰链a出现了比较明显的磨损，选定铰链a作为磨损退化分部件，磨损量用磨损导致铰链间隙的变化衡量。同时，铰链a的磨损导致相应部件的几何长度变化，进一步影响功能量的大小。当功能量与理想情况下的功能量误差超过允许值，锁机构即发生功能失效，当磨损退化部件的磨损量超过允许值，锁机构发生部件磨损退化失效。步骤二中，利用非线性漂移布朗运动对铰链a建立退化模型：非线性漂移布朗运动的公式为：x1(t)＝σb(t)+μtq+d0，式中x1(t)为退化量，b(t)为标准布朗运动，σ为第一参数,μ为第二参数，d0为退化量的初值，tq为时间项，q为指数项；根据飞机舱门上位锁机构的设计要求，铰链a的磨损量不能超过既定的阈值xth，故铰链a失效域为{x1(t)＞xth}，其可靠度表示为：退化模型中相关参数为表1所示：表1σ(mm)μ(mm)qd0(mm)xth(mm)0.10.00111.150.050.6为了正确描述铰链a的退化过程，采用非线性漂移布朗运动建模。典型的漂移布朗运动公式为x(t)＝σb(t)+μt+d0，其中x(t)表示退化量，b(t)为标准布朗运动，σ,μ为参数，d0表示退化量的初值。典型漂移布朗运动公式中的时间项t是线性的，无法更准确的表述退化过程中的速度变化，为了更好的解决此问题，步骤二中在非线性漂移的布朗运动中应用了指数形式的时间项。将典型漂移布朗运动公式中的t转变为tq，当指数项q的值大于1时，整个退化过程是加速的。当新的时间项引入后，根据概率论和数理统计理论，退化量x1(t)在t时刻遵循均值为d0+μtq，方差为σ2t正态分布，即x1(t)～n(d0+μtq,σ2t)。步骤三中，功能量由锁机构各个部件共同完成，因此所述运动功能量为各个部件几何长度的函数，用lpm表示该锁机构的运动功能量，lpm＝g(lcb,lab,lbd,ldf,lef)，定义摇臂1的上连接点为c点，中连接点为b点，下连接点为d点，铰链a的连接点为a点，锁钩4的旋转连接点为e点，状态连接点为f点，a点与b点之间通过摇臂中连杆2连接，f点和d点之间通过摇臂下连杆3连接，lcb为c点与b点之间的杆件长度，lab为a点与b点之间的杆件长度、lbd为b点与d点之间的杆件长度、ldf为d点与f点之间的杆件长度、lef为e点与f点之间的杆件长度；摇臂中连杆2的分布特性是随着时间变化的，根据有效长度理论，锁机构功能量的函数进一步表示为：lpm＝g(lcb,lab(t),lbd,ldf,lef)；当锁机构功能量与理想状态下功能量的差值超过允许值的时候锁机构发生功能失效，该功能失效的失效域为{lpm(t)-lpm理想＞δ}，式中δ为锁机构功能量与理想状态下功能量的差值允许值；其可靠度为：rps(t)＝p{lpm(t)-lpm理想＜δ}＝p(lpm(t)＜lpm理想+δ)；通过代理模型建立锁机构运动功能量分布和各个分部件几何长度分布的传递关系，选用一次响应面模型，各个分部件的分布性质如表2所示：表2步骤三中，由于铰链a是磨损铰链，磨损会影响到与之相连杆件的几何长度，因此杆件ab即摇臂中连杆2的分布特性是随着时间变化的。对于同一批机构，各个构件的几何长度是符合某一分布的变量，因此可以通过代理模型建立锁机构功能量分布和各个构件几何长度分布的传递关系。在这里选用一次响应面模型。各个部件的分布性质如表2所示。步骤四中，针对磨损退化量和锁机构功能退化量通过fgmcopula函数建立联合分布函数：fgmcopula函数的表达式为：c(u1,u2)＝u1u2(1+α(1-u1)(1-u2))，其中u1：0≤u1≤1，为铰链a的磨损退化量的边缘分布，且由于退化的存在该分布随时间发生变化；u2：0≤u2≤1，为锁机构功能量的边缘分布；α：|α|≤1，为u1与u2的相关系数；对于飞机舱门上位锁机构，联合分布函数为：通过前三个步骤，可以知晓各个部件和锁机构功能量在任意时刻t下的分布，针对磨损退化量和锁机构功能退化量通过fgmcopula函数建立其联合分布函数，并据此进行进一步的分析。步骤五中，对于磨损退化失效，其竞争失效概率f1为：对于锁机构功能失效，其竞争失效概率f2为：锁机构的整体可靠度r锁机构为：步骤五中计算得到的竞争失效概率结果与不考虑竞争情况下的失效概率和利用蒙特卡洛方法计算得到的实效概率结果进行对比，可以进一步印证本发明方法的准确性和实用性，如图5及图6所示。任一种失效的出现都会导致系统的失效，利用得到的联合概率密度函数，计算锁系统的整体可靠度，并得到可靠度随时间的变化规律。计算结果如图7所示。通过对各个曲线图的分析可知，本发明所提的方法可以有效的针对处于部件磨损失效和功能失效下的锁机构进行失效分析和可靠度计算，且计算结果更接近实际工程情况。与此同时，与蒙特卡洛计算所得到的结果相比，本方法所得到的计算结果具有足够的计算精度，可保证在使用过程的准确性。当前第1页12