考虑不完全维修的风机齿轮箱优化检修策略(2)

风机厂石家庄风机厂石家庄风机石家庄市风机厂石家庄风机维修石家庄风机销售

1石家庄风机齿轮箱比例强度模型
1.1石家庄风机齿轮箱的维修方式针对石家庄风机齿轮箱这一复杂多部件可修系统而言，在寿命期限内发生故障，不需要对整个设备进行更换，而仅更换或维修部分部件，来维持所预定的功能，故在其整个寿命期内将发生多次故障及维修活动。根据维修类型和维修程度的不同，维修可分为如下几种 [15-16]（1）完全维修。指设备修复如新，如齿轮箱整机更换就属于完全维修。其维修因子：M =0。（2）最小维修。指设备修复如旧，如齿轮箱零部件加固、润滑油更换等属于最小维修。其维修因子 M =1。（3）不完全维修。指设备维修后，功能得以恢复，设备状态介于完全维修和最小维修之间。如齿轮箱零部件更换或维修属于不完全维修。其维修因子 0< M <1。对于石家庄风机齿轮箱常见的故障，如轮齿损坏、轴承磨损等故障类型而言，通常仅对其部分部件进行更换或维修而不进行整机更换。故虽然其功能得以恢复，但其健康状态通常处于“修复如新”与“修复如旧”之间，即不完全维修(0<M<1)。因此，本文针对不完全维修建模研究维修活动与实际情况更加吻合。不完全维修对系统强度函数的影响如图1所示。其中： ( )v t 为齿轮箱的强度函数；mv 为失效临界强度值；1t 、2t 、3t 为齿轮箱故障时刻(假设维修时间相对设备运行时间可以忽略不计)。
1.2 比例强度模型在 PIM 中，每次故障事件是随机的，故障事件通常被视为一个非齐次 Possion 过程，则系统发生失效的强度函数为T0{ ( ) ( ) 1| ( )}( ) lim( )exp( ( ))tP N t t N t Z tv ttv t Z t γ∆ →∞+ ∆ − ≥= =∆⋅(1)其中： t 是齿轮箱的运行时间； ( )v t 为齿轮箱在 t 时刻的强度函数；0 ( )v t 为基本强度函数； ( ) N t 为设备在区间(0, t )内失效次数； ( ) Z t 是系统在时刻t的协变量； γ 为协变量回归参数。比例强度模型包括全参数型和半参数型，当失效规律已知(0 ( )v t 已知)时，为全参数型PIM，能同时估计出基本强度函数的参数和协变量的回归系数，参数估计的准确率相对较高；当失效规律未知(0 ( )v t 未知)时，采用半参数型PIM。参数的各种估计方法，以基本强度函数服从威布尔分布的PIM的估计结果最为准确。 研究表明， 齿轮箱的故障时间服从威布尔分布，其强度函数为(2)其中： β 为形状参数； η 为尺度参数；( ) Z t 为系统在 t 时刻的协变量； γ 为协变量的回归参数。协变量可分为外部协变量和内部协变量。外部协变量包括石家庄风机齿轮箱的运行环境温度、工作负荷等，而内部协变量包括反映齿轮箱故障征兆的振动数据、温度数据等。考虑维修效应对强度函数的影响，既可以将维修效应作为衰减因子得到虚拟寿命过程，也可以将其作为协变量叠加到强度函数中，本文采用后者。
1.3 石家庄风机齿轮箱 PIM 的参数估计利用石家庄风机齿轮箱故障数据对PIM 的参数进行极大似然估计。首先验证故障数据是否符合威布尔分布以确定模型选择的正确性，然后构建似然函数，并在 Matlab 中应用 Newton-Raphson 迭代算法解非线性方程组，计算出估计值，具体步骤如下。第一步：样本分布的检验。采用最常见的威布尔分布检验图，检验图反映了威布尔分布的基本特征，是用于检验样本数据是否符合威布尔分布规律的最直观方法。 若样本服从威布尔分布，则该曲线应近似地拟合为一条斜率大于0的直线， 直线的斜率即为形状参数的估计值， 如图2所示，其故障数据服从威布尔分布。第二步：似然函数的构造。