一种动车组子系统关键部件运用可靠性评估方法
一种动车组子系统关键部件运用可靠性评估方法
【技术领域】
[0001] 本发明涉及轨道车辆安全领域，尤其是涉及一种动车组子系统关键部件运用可靠 性评估方法。
【背景技术】
[0002] 动车组是极端服役条件下的具有先进功能的复杂技术系统。近年来，一方面我国 动车组上线运用数量日益增多，另一方面，国内在动车组运用与维修方面也面临着先期经 验匮乏的客观事实，评估动车组子系统关键部件运用可靠性水平，一方面可以为保障车辆 的运行安全提供依据，另一方面还可以为列车的修程修制优化提供支撑，从而实现列车运 行过程中"安全性"与"经济性"的有机结合。为了切实保证动车组运行安全和经济可用， 迫切需要一套基于上线运用现场数据的动车组子系统关键部件运用可靠性评估方法。然而 目前还没有有效的可靠性评估方法。
【发明内容】
[0003] 本发明的目的就是为了克服上述现有技术存在的空白而提供一种有效的动车组 子系统关键部件运用可靠性评估方法。
[0004] 本发明的目的可以通过以下技术方案来实现：
[0005] 一种动车组关键部件运行可靠性评估方法，包括以下步骤：
[0006] (1)构建标准化动车组子系统结构树，该动车组子系统结构树的最底层为动车组 运用现场最小不可拆分部件，每棵树所有层次的部件名称和数量相同；
[0007] (2)根据动车组实际上线运用情况，按投入运用时间划分车组批次；
[0008] (3)记录同一批次动车组的标准化部件现场发生故障的时间点Xi，i= 1，2,3. ..，r，r为故障个数；
[0009] (4)根据动车组现场动态维修记录的离散时间与运行里程信息的关系，计算对应 的动车组运行里程Xpi= 1，2, 3. . .，r;
[0010] (5)计算Xi的均值0、方差s、二阶矩y2、三阶矩y3、四阶矩y4、偏度Cs、峰度(； 和对数化样本的偏度Cs'、峰度C/ ;
[0011] (6)计算部件故障样本的平均失效率，得到失效率表；
[0012] (7)选择部件寿命分布模型，确定部件可靠度函数估计方法；
[0013] (8)确定单个部件平均故障间隔时间和可靠度函数；
[0014] (9)确定车组级部件平均故障间隔时间和可靠度函数；
[0015] (10)利用可靠度函数对各个部件进行可靠性评估。
[0016] 所述的步骤（2)中，从第一批动车投入运行开始，同一批次的所有车组投入运用 时间相差不超过一年。
[0017] 所述的步骤（4)中，采用分段线性插值的方式来获得Xi。
[0018] 所述的步骤（5)具体包括以下步骤：
[0019] (5-1)根据下列公式计算Xi的均值0、方差s、二阶矩y2、三阶矩y3、四阶矩y4、 偏度Cs、峰度C6:
[0025] Cs=y3/03
[0026] ce= y4/〇4；
[0027] (5-2)将Xi作对数变换，重复步骤（5-1)得到对数化样本的偏度Cs'、峰度Ce'。
[0028] 所述的步骤（6)具体包括以下步骤：
[0029] (6-1)将动车组总运行时间分成k个时间区间Ax，计算每个时间区间长度：
[0030] Ax= (La-Sm)/k
[0031] 其中，k=l+3.31gr，LaSXi的最大值，S"^Xi的最小值；
[0032] (6-2)计算部件故障样本的平均失效率RAX,.):
[0034] 其中，为第i个时间区间Axi内的失效率频数，i= 1，2, 3. . .，k，r为进 入第i个时间区间时的样本数，1^4=r-rh，:^指进入第i个时间区间的累积失效率。
[0035] 所述的步骤（7)具体包括以下步骤：
[0036] (7-1)判断部件故障数据个数，若大于10,则进入步骤（7-2);否则返回步骤 (7-1)；
[0037] (7-2)若| 0 -s| < ( 0 +s)/5,则部件寿命服从指数分布，否则进入步骤（7-3);
[0038] (7-3)若|Cs |〈0. 5且|Ce_3 |〈0. 5,则部件寿命服从正态分布，否则进入步骤 (7-4)；
[0039] (7-4)若|Cs' |〈0. 5且|C/ -3 |〈0. 5,则部件寿命服从对数正态分布，否则进入步 骤（7-5);
[0040] (7-5)根据步骤（6)中得到的平均失效率表，计算AAi:
[0042] 若
对于i= 1，2, ...k均成立，则部件寿命服从指数分布，否则进 入步骤（7-6);
[0043] (7-6)对步骤（7-5)中AAp记AApO的个数为a，AA^ 〇个数为b，若a/ b多3/4,则部件寿命服从威布尔分布，否则部件寿命分布无规律，可靠性评估采用非参数 方法。
[0044]所述的步骤⑶具体为：
[0045] (8-1)若部件寿命服从指数分布，则将动车组部件寿命数据作为无替换定时截尾 情形，即有n件产品投入使用，到规定的时间X(l进行数据收集，依照时间的先后记录截止到 规定时间的失效时间，对应的里程数据为XiS X2彡...彡X y
[0046] 根据定时截尾样本数据，获得该样本的可靠度似然函数L( 0 ):
[0048]其中：
为总运行时间，对L( 0 )取对数并求导，求解似然 方程，得到0和A的极大似然点估计为：
[0050] (8-2)若部件寿命服从威布尔分布，则密度函数为：
[0052] 根据定时截尾样本数据，得到该样本的可靠度似然函数：
[0054] 使用参数估计迭代数值求解算法求解以下方程组：
[0056] 得到可靠度函数为：
[0058] (8-3)若部件寿命服从对数正态分布，则密度函数为：
[0060] 定时截尾时可靠度似然函数为：
[0062] 设ZQ=(InxQ-y) / 〇，标准正态分布函数①（_ZQ) = 1-①（ZQ)，并且记①（ZQ)为 标准正态分布密度函数，则似然方程为：
[0064] 使用参数估计近似数值求解算法求解上述方程组，即可得到参数y、〇的极大似 然估计，从而得到可靠度函数；
[0065] (8-4)若部件寿命服从正态分布，正态分布的求解为：将对数正态分布密度函数 f(x)中的lnx替换为X，其余步骤与对数正态分布的求解方法相同；
[0066] (8-5)若部件寿命分布无规律，令收集到的部件数据为XpXyXi. . .Xy当Xi是故障 数据时，令Si=l;当&是右截尾数据时，令S,=〇,将数据记为（Xi，6山1= 1，2，...，1*， 将这些七按从小到大排列，为xx2彡..?彡x ^
[0067] 可靠度函数的乘积限估计为：
[0069] 平均寿命估计如下：
[0071] 本发明具有以下优点：
[0072] 1)针对我国动车组运用现场实际情况，将动车组标准化并建立子系统结构树，有 助于建立标准化数据分析模型，解决动车组上线运用可靠性分析缺乏的问题。
[0073] 2)通过对故障数据分布类型进行判别，根据分布类型计算可靠度函数，提高可靠 性估计的准确度。
[0074] 3)本发明的基于现场数据的动车组关键部件运行可靠性评估方法也适用于其他 轨道车辆部件的运行可靠性分析。
【附图说明】
[0075] 图1为本发明的流程图。
【具体实施方式】
[0076] 下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案 为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于 下述的实施例。
[0077] 如图1所示，本实施例的方法实现具体包括以下步骤：
[0078] S01，构建动车组子系统结构树。
[0079] 动车组子系统结构树的最低层为动车组运用现场最小不可拆分部件，每棵树所有 层次的部件名称和数量相同。
[0080] S02,划分车组批次。
[0081] 根据动车组实际上线运用情况，可将车组按投入运用时间分成不同批次。如，不同 批次车组投入运用时间间隔一般在2到3年以上，故投入运用时间间隔不超过一年的车组 列可认为同一批次。
[0082] S03,故障时间点记录。
[0083] 记录同一批次动车组的标准化部件现场发生故障的时间点ivi= 1，2, 3. ..r，r 为故障个数。
[0084] S04,故障运行时间计算。
[0085] 根据动车组现场动态维修记录的离散时间_运行里程信息，采用分段线性插值的 方式计算各个时间点Xi上对应的动车组运行里程xpi= 1，2, 3. . .r。
[0086] S05,统计数据参数计算。
[0087] 将动车组总运行时间分成k个时间区间Ax，计算Xi的均值0、方差s、二阶矩y2、 三阶矩U2、四阶矩U4、偏度Cs、峰度(；和对数化样本的偏度Cs'、峰度C/，具体包括以下步 骤：
[0088] S0501，将动车组总运行时间分成k个时间区间Ax，计算每个时间区间长度
[0089] Ax= (La-Sm)/k
[0090]其中，k= 1+3. 31gr，LaSx1最大值，SmSx1最小值。
[0091] S0502,根据下列公式计算Xi的均值0、方差s、二阶矩y2、三阶矩y3、四阶矩y4、 偏度Cs、峰度C6:
[0097] Cs=y3/03
[0098] ce=y4/〇4。
[0099] S0503,将Xi作对数变换，重复步骤（5-2)得到对数化样本的偏度Cs'、峰度Ce'。
[0100] S06,计算失效率表。
[0101] 计算部件故障样本的平均失效率ICAxj，得到失效率表：
[0103] 其中，Ari为第i个时间区间Ari内的失效率频数，i= 1，2,3. 为进入 第i个时间区间时的样本数，1^4=r-rh，:^指进入第i个时间区间的累积失效率。
[0104] S07,寿命分布模型选择。
[0105] 选择部件寿命分布模型，确定部件可靠度函数估计方法，具体包括以下步骤：
[0106] S0701，判断部件故障数据个数，若大于10,则进入步骤S0702。
[0107] S0702,若| 0 -s| < ( 0 +s)/5,则部件寿命服从指数分布，否则进入步骤S0703。
[0108] S0703,若|CS|〈0. 5且|Ce-3|〈0. 5,则 部件寿命服从正态分布，否则进入步骤 S0704。
[0109] S0704,若|CS' |〈0. 5且|Ce' -3|〈0. 5,则部件寿命服从对数正态分布，否则进入步 骤S0705。
[0110] S0705,根据步骤（6)中得到的平均失效率表，计算AAi:
[0112] 若
对于i= 1，2，...k均成立，则部件寿命服从指数分布，否则进 入步骤S0706。
[0113]S0706,对步骤S0705中A记A0的个数为a，AAi<0个数为b，若a/ b多3/4,则部件寿命服从威布尔分布，否则部件寿命分布无规律，可靠性评估采用非参数 方法。
[0114]S08,可靠度函数计算。具体包括以下步骤：
[0115] S0801，若部件寿命服从指数分布，则将动车组部件寿命数据作为无替换定时截尾 情形，即有n件产品投入使用，到规定的里程X(l进行数据收集，依照里程的先后记录截止到 规定里程的失效里程得Xl<x2彡...彡x^
[0116] 根据定时截尾样本数据，写出该样本的可靠度似然函数L(0):
[0118]其中，
，为总运行时间，对L( 0 )取对数并求导，求解似然 方程，得到0和A的极大似然点估计为：
[0120] S0802,若部件寿命服从威布尔分布，则密度函数为：
[0122] 威布尔分布可以利用概率值很容易地推断出它的分布参数，n为比例参数，m为 形状参数。代入里程数据\，根据定时截尾样本数据，得到该样本的可靠度似然函数：
[0124]其中，L是似然函数，S,是似然函数参数，使用用参数估计迭代数值求解算法求解 以下方程组：
[0126] 得到可靠度函数为：
[0128] 其中，参数估计迭代数值求解算法具体步骤为：
[0129] S0802a，记：
[0131] S0802b，选择初值mQ= 1。
[0133] S0802d，重复步骤S0802c，直到 |mk+1-mk| < 0. 0001，得到m=mk+1。
[0134] S0802e，根据下式计算n:
[0136] S0803,若部件寿命服从对数正态分布，则密度函数为：
[0138] 定时截尾时可靠度似然函数为：
[0140] 设ZQ=(InxQ-li)/〇，标准正态分布函数①（_ZQ) = 1-①（Z。)，并且记①（ZQ)为 标准正态分布密度函数，则似然方程为：
[0142] 使用参数估计近似数值求解算法求解上述方程组，即可得到参数y、〇的极大似 然估计，从而得到可靠度函数。
[0143] 其中，参数估计近似数值求解算法具体包括以下步骤：
[0144] S0803a，选择初值yQ，〇0。
[0147]S〇8〇3b，计算
[0149] 其中U(zs) = [V(zs)_zs]。
[0150] S0803c，重复步骤（8-3b)，直到 |yk+1-yk | < 0? 0001 且
，得 到最终的y，〇极大似然估计。
[0151] S0804,若部件寿命服从正态分布，正态分布的求解为：将对数正态分布密度函数 f(t)中的lnx替换为X，其余步骤与对数正态分布的求解方法相同。
[0152] S0805,若部件寿命分布无规律，令收集到的部件数据为Xpx2,Xy. .Xp当Xi是故 障数据时，令\=1;当、是右截尾数据时，令S,= 〇,这样数据可记为（Xi，\)，1 = 1，2, . . .，r，将这些\按从小到大排列，对应的里程数据为x#x2彡...彡x^
[0153] 可靠度函数的乘积限估计为：
[0155] 平均寿命估计如下：
[0157] 在获得上述单个部件平均故障间隔时间和可靠度函数后，可根据动车组编组情 况，计算车组级部件平均故障间隔时间和可靠度。
[0158] 本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发 明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的 普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各 种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。
【主权项】
1. 一种动车组关键部件运行可靠性评估方法，其特征在于，包括以下步骤： (1) 构建标准化动车组子系统结构树，该动车组子系统结构树的最底层为动车组运用 现场最小不可拆分部件，每棵树所有层次的部件名称和数量相同； (2) 根据动车组实际上线运用情况，按投入运用时间划分车组批次，得到有η件产品投 入使用； (3) 记录同一批次动车组的标准化部件现场发生故障的时间点h，i = 1，2, 3. ..，r，r 为故障个数； (4) 根据动车组现场动态维修记录的离散时间与运行里程信息的关系，计算各个时间 点h上对应的动车组运行里程X i，i = 1，2, 3. . .，r ; (5) 计算Xi的均值Θ、方差S、二阶矩μ 2、三阶矩μ3、四阶矩μ4、偏度cs、峰度(；和对 数化样本的偏度Cs'、峰度C/ ; (6) 计算部件故障样本的平均失效率，得到失效率表； (7) 选择部件寿命分布模型，确定部件可靠度函数估计方法； (8) 确定单个部件平均故障间隔时间和可靠度函数； (9) 确定车组级部件平均故障间隔时间和可靠度函数； (10) 利用可靠度函数对各个部件进行可靠性评估。2. 根据权利要求1所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤（2)中，从第一批动车投入运行开始，同一批次的所有车组投入运用时间相差不 超过一年。3. 根据权利要求1所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤（4)中，采用分段线性插值的方式来获得4. 根据权利要求1所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤（5)具体包括以下步骤： (5-1)根据下列公式计算Xi的均值Θ、方差S、二阶矩μ 2、三阶矩μ 3、四阶矩μ 4、偏度 cs、峰度(；：Ce= y 4/s ; (5-2)将Xi作对数变换，重复步骤（5-1)得到对数化样本的偏度C s'、峰度(V。5. 根据权利要求1所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤（6)具体包括以下步骤： (6-1)将动车组总运行里程分成k个里程区间Λ X，计算每个里程区间长度： Δ X = (La-Sm) /k 其中，k = 1+3. 31gn，LaS X i的最大值，S X i的最小值； (6-2)计算部件故障样本的平均失效率无(Δχ,):其中，Λ A为第i个里程区间Λ X i内的失效率频数，i = 1，2, 3. ..，k，r 为进入第 i个时间区间时的样本数，rs，H= r-r η，:^指进入第i个时间区间的累积失效率。6. 根据权利要求1所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤（7)具体包括以下步骤： (7-1)判断部件故障数据个数，若大于10,则进入步骤（7-2);否则返回步骤（7-1); (7-2) I若I Θ-S I〈（ Θ+s)/5,则部件寿命服从指数分布，否则进入步骤（7-3); (7-3)若|CS|〈0.5且|Ce-3 I〈0.5,则部件寿命服从正态分布，否则进入步骤（7-4); (7-4)若|CS' I〈0.5且|C/-3 I〈0.5,则部件寿命服从对数正态分布，否则进入步骤 (7-5)； (7-5)根据步骤（6)中得到的平均失效率表，计算Λ λ i:若夂对于i = 1，2，... k均成立，则部件寿命服从指数分布，否则进入步 骤（7-6); (7-6)对步骤（7-5)中Λ λ i，记Λ λ，〇的个数为a，Λ λ'〇个数为b，若a/b彡3/4， 则部件寿命服从威布尔分布，否则部件寿命分布无规律，可靠性评估采用非参数方法。7. 根据权利要求6所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤⑶具体为： (8-1)若部件寿命服从指数分布，则将动车组部件寿命数据作为无替换定时截尾情形， 即有η件产品投入使用，到规定的时间Xtl进行数据收集，依照时间的先后记录截止到规定 时间的失效时间获得对应的走行里程数据，其中： X1^ X2^... ^ X X 0? r<n 根据定时截尾样本数据，获得该样本的可靠度似然函数L( Θ ):其中，，为总运行里程，对L(0)取对数并求导，求解似然方 程，得到Θ和λ的极大似然点估计为：(8-2)若部件寿命服从威布尔分布，则密度函数为：故障样本走行里程为Xi，其中截尾时间Χ(ι，引入似然函数参数δ i，根据定时截尾样本 数据，得到该样本的可靠度似然函数：使用参数估计迭代数值求解算法求解以下方程组：得到可靠度函数为：(8-3)若部件寿命服从对数正态分布，则密度函数为：定时截尾时可靠度似然函数为：设Ztl= (Inxtl-μ)/s，标准正态分布函数Φ (-?) = l-φ (Ztl),并且记Φ (Ztl)为标准正 态分布密度函数，则似然方程为：使用参数估计近似数值求解算法求解上述方程组，即可得到参数μ、S的极大似然估 计，从而得到可靠度函数； (8-4)若部件寿命服从正态分布，正态分布的求解为：将对数正态分布密度函数f(t) 中的Inx替换为X，其余步骤与对数正态分布的求解方法相同； (8-5)若部件寿命分布无规律，令收集到的部件数据为Xl，x2, Xi. 当Xi是故障数据 时，令δ i= 1 ;当X肩右截尾数据时，令δ F 〇,将数据记为（X D δ山i = 1，2, . . .，r，将 这些\按从小到大排列，得X # X 2彡..·彡X y 可靠度函数的乘积限估计为：平均寿命估计如下：
【专利摘要】本发明涉及轨道车辆安全领域，具体涉及一种动车组子系统关键部件运用可靠性评估方法，主要包括：建立动车组子系统结构树并进行现场数据预处理，获得故障间隔运行里程数据；计算所得样本的均值、方差、二阶矩、三阶矩、四阶矩、偏度Cs、峰度Ce、平均失效率以及对数化样本的偏度Cs’、峰度Ce’；确定部件寿命分布和参数估计算法；计算平均故障间隔时间和可靠度。本发明有效处理动车组实际运行数据并对动车组部件进行运行可靠性估计，解决了动车组上线运用可靠性分析缺乏的问题，同时，本发明的基于现场数据的动车组关键部件运行可靠性评估方法也适用于其他轨道车辆部件的运行可靠性分析。
【IPC分类】G06F19/00
【公开号】CN104899423
【申请号】CN201510225640
【发明人】罗敏, 马哲, 陆正刚, 李道军, 李 浩
【申请人】同济大学
【公开日】2015年9月9日
【申请日】2015年5月6日
【技术领域】
[0001] 本发明涉及轨道车辆安全领域，尤其是涉及一种动车组子系统关键部件运用可靠 性评估方法。
【背景技术】
[0002] 动车组是极端服役条件下的具有先进功能的复杂技术系统。近年来，一方面我国 动车组上线运用数量日益增多，另一方面，国内在动车组运用与维修方面也面临着先期经 验匮乏的客观事实，评估动车组子系统关键部件运用可靠性水平，一方面可以为保障车辆 的运行安全提供依据，另一方面还可以为列车的修程修制优化提供支撑，从而实现列车运 行过程中"安全性"与"经济性"的有机结合。为了切实保证动车组运行安全和经济可用， 迫切需要一套基于上线运用现场数据的动车组子系统关键部件运用可靠性评估方法。然而 目前还没有有效的可靠性评估方法。
【发明内容】
[0003] 本发明的目的就是为了克服上述现有技术存在的空白而提供一种有效的动车组 子系统关键部件运用可靠性评估方法。
[0004] 本发明的目的可以通过以下技术方案来实现：
[0005] 一种动车组关键部件运行可靠性评估方法，包括以下步骤：
[0006] (1)构建标准化动车组子系统结构树，该动车组子系统结构树的最底层为动车组 运用现场最小不可拆分部件，每棵树所有层次的部件名称和数量相同；
[0007] (2)根据动车组实际上线运用情况，按投入运用时间划分车组批次；
[0008] (3)记录同一批次动车组的标准化部件现场发生故障的时间点Xi，i= 1，2,3. ..，r，r为故障个数；
[0009] (4)根据动车组现场动态维修记录的离散时间与运行里程信息的关系，计算对应 的动车组运行里程Xpi= 1，2, 3. . .，r;
[0010] (5)计算Xi的均值0、方差s、二阶矩y2、三阶矩y3、四阶矩y4、偏度Cs、峰度(； 和对数化样本的偏度Cs'、峰度C/ ;
[0011] (6)计算部件故障样本的平均失效率，得到失效率表；
[0012] (7)选择部件寿命分布模型，确定部件可靠度函数估计方法；
[0013] (8)确定单个部件平均故障间隔时间和可靠度函数；
[0014] (9)确定车组级部件平均故障间隔时间和可靠度函数；
[0015] (10)利用可靠度函数对各个部件进行可靠性评估。
[0016] 所述的步骤（2)中，从第一批动车投入运行开始，同一批次的所有车组投入运用 时间相差不超过一年。
[0017] 所述的步骤（4)中，采用分段线性插值的方式来获得Xi。
[0018] 所述的步骤（5)具体包括以下步骤：
[0019] (5-1)根据下列公式计算Xi的均值0、方差s、二阶矩y2、三阶矩y3、四阶矩y4、 偏度Cs、峰度C6:
[0025] Cs=y3/03
[0026] ce= y4/〇4；
[0027] (5-2)将Xi作对数变换，重复步骤（5-1)得到对数化样本的偏度Cs'、峰度Ce'。
[0028] 所述的步骤（6)具体包括以下步骤：
[0029] (6-1)将动车组总运行时间分成k个时间区间Ax，计算每个时间区间长度：
[0030] Ax= (La-Sm)/k
[0031] 其中，k=l+3.31gr，LaSXi的最大值，S"^Xi的最小值；
[0032] (6-2)计算部件故障样本的平均失效率RAX,.):
[0034] 其中，为第i个时间区间Axi内的失效率频数，i= 1，2, 3. . .，k，r为进 入第i个时间区间时的样本数，1^4=r-rh，:^指进入第i个时间区间的累积失效率。
[0035] 所述的步骤（7)具体包括以下步骤：
[0036] (7-1)判断部件故障数据个数，若大于10,则进入步骤（7-2);否则返回步骤 (7-1)；
[0037] (7-2)若| 0 -s| < ( 0 +s)/5,则部件寿命服从指数分布，否则进入步骤（7-3);
[0038] (7-3)若|Cs |〈0. 5且|Ce_3 |〈0. 5,则部件寿命服从正态分布，否则进入步骤 (7-4)；
[0039] (7-4)若|Cs' |〈0. 5且|C/ -3 |〈0. 5,则部件寿命服从对数正态分布，否则进入步 骤（7-5);
[0040] (7-5)根据步骤（6)中得到的平均失效率表，计算AAi:
[0042] 若
对于i= 1，2, ...k均成立，则部件寿命服从指数分布，否则进 入步骤（7-6);
[0043] (7-6)对步骤（7-5)中AAp记AApO的个数为a，AA^ 〇个数为b，若a/ b多3/4,则部件寿命服从威布尔分布，否则部件寿命分布无规律，可靠性评估采用非参数 方法。
[0044]所述的步骤⑶具体为：
[0045] (8-1)若部件寿命服从指数分布，则将动车组部件寿命数据作为无替换定时截尾 情形，即有n件产品投入使用，到规定的时间X(l进行数据收集，依照时间的先后记录截止到 规定时间的失效时间，对应的里程数据为XiS X2彡...彡X y
[0046] 根据定时截尾样本数据，获得该样本的可靠度似然函数L( 0 ):
[0048]其中：
为总运行时间，对L( 0 )取对数并求导，求解似然 方程，得到0和A的极大似然点估计为：
[0050] (8-2)若部件寿命服从威布尔分布，则密度函数为：
[0052] 根据定时截尾样本数据，得到该样本的可靠度似然函数：
[0054] 使用参数估计迭代数值求解算法求解以下方程组：
[0056] 得到可靠度函数为：
[0058] (8-3)若部件寿命服从对数正态分布，则密度函数为：
[0060] 定时截尾时可靠度似然函数为：
[0062] 设ZQ=(InxQ-y) / 〇，标准正态分布函数①（_ZQ) = 1-①（ZQ)，并且记①（ZQ)为 标准正态分布密度函数，则似然方程为：
[0064] 使用参数估计近似数值求解算法求解上述方程组，即可得到参数y、〇的极大似 然估计，从而得到可靠度函数；
[0065] (8-4)若部件寿命服从正态分布，正态分布的求解为：将对数正态分布密度函数 f(x)中的lnx替换为X，其余步骤与对数正态分布的求解方法相同；
[0066] (8-5)若部件寿命分布无规律，令收集到的部件数据为XpXyXi. . .Xy当Xi是故障 数据时，令Si=l;当&是右截尾数据时，令S,=〇,将数据记为（Xi，6山1= 1，2，...，1*， 将这些七按从小到大排列，为xx2彡..?彡x ^
[0067] 可靠度函数的乘积限估计为：
[0069] 平均寿命估计如下：
[0071] 本发明具有以下优点：
[0072] 1)针对我国动车组运用现场实际情况，将动车组标准化并建立子系统结构树，有 助于建立标准化数据分析模型，解决动车组上线运用可靠性分析缺乏的问题。
[0073] 2)通过对故障数据分布类型进行判别，根据分布类型计算可靠度函数，提高可靠 性估计的准确度。
[0074] 3)本发明的基于现场数据的动车组关键部件运行可靠性评估方法也适用于其他 轨道车辆部件的运行可靠性分析。
【附图说明】
[0075] 图1为本发明的流程图。
【具体实施方式】
[0076] 下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案 为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于 下述的实施例。
[0077] 如图1所示，本实施例的方法实现具体包括以下步骤：
[0078] S01，构建动车组子系统结构树。
[0079] 动车组子系统结构树的最低层为动车组运用现场最小不可拆分部件，每棵树所有 层次的部件名称和数量相同。
[0080] S02,划分车组批次。
[0081] 根据动车组实际上线运用情况，可将车组按投入运用时间分成不同批次。如，不同 批次车组投入运用时间间隔一般在2到3年以上，故投入运用时间间隔不超过一年的车组 列可认为同一批次。
[0082] S03,故障时间点记录。
[0083] 记录同一批次动车组的标准化部件现场发生故障的时间点ivi= 1，2, 3. ..r，r 为故障个数。
[0084] S04,故障运行时间计算。
[0085] 根据动车组现场动态维修记录的离散时间_运行里程信息，采用分段线性插值的 方式计算各个时间点Xi上对应的动车组运行里程xpi= 1，2, 3. . .r。
[0086] S05,统计数据参数计算。
[0087] 将动车组总运行时间分成k个时间区间Ax，计算Xi的均值0、方差s、二阶矩y2、 三阶矩U2、四阶矩U4、偏度Cs、峰度(；和对数化样本的偏度Cs'、峰度C/，具体包括以下步 骤：
[0088] S0501，将动车组总运行时间分成k个时间区间Ax，计算每个时间区间长度
[0089] Ax= (La-Sm)/k
[0090]其中，k= 1+3. 31gr，LaSx1最大值，SmSx1最小值。
[0091] S0502,根据下列公式计算Xi的均值0、方差s、二阶矩y2、三阶矩y3、四阶矩y4、 偏度Cs、峰度C6:
[0097] Cs=y3/03
[0098] ce=y4/〇4。
[0099] S0503,将Xi作对数变换，重复步骤（5-2)得到对数化样本的偏度Cs'、峰度Ce'。
[0100] S06,计算失效率表。
[0101] 计算部件故障样本的平均失效率ICAxj，得到失效率表：
[0103] 其中，Ari为第i个时间区间Ari内的失效率频数，i= 1，2,3. 为进入 第i个时间区间时的样本数，1^4=r-rh，:^指进入第i个时间区间的累积失效率。
[0104] S07,寿命分布模型选择。
[0105] 选择部件寿命分布模型，确定部件可靠度函数估计方法，具体包括以下步骤：
[0106] S0701，判断部件故障数据个数，若大于10,则进入步骤S0702。
[0107] S0702,若| 0 -s| < ( 0 +s)/5,则部件寿命服从指数分布，否则进入步骤S0703。
[0108] S0703,若|CS|〈0. 5且|Ce-3|〈0. 5,则 部件寿命服从正态分布，否则进入步骤 S0704。
[0109] S0704,若|CS' |〈0. 5且|Ce' -3|〈0. 5,则部件寿命服从对数正态分布，否则进入步 骤S0705。
[0110] S0705,根据步骤（6)中得到的平均失效率表，计算AAi:
[0112] 若
对于i= 1，2，...k均成立，则部件寿命服从指数分布，否则进 入步骤S0706。
[0113]S0706,对步骤S0705中A记A0的个数为a，AAi<0个数为b，若a/ b多3/4,则部件寿命服从威布尔分布，否则部件寿命分布无规律，可靠性评估采用非参数 方法。
[0114]S08,可靠度函数计算。具体包括以下步骤：
[0115] S0801，若部件寿命服从指数分布，则将动车组部件寿命数据作为无替换定时截尾 情形，即有n件产品投入使用，到规定的里程X(l进行数据收集，依照里程的先后记录截止到 规定里程的失效里程得Xl<x2彡...彡x^
[0116] 根据定时截尾样本数据，写出该样本的可靠度似然函数L(0):
[0118]其中，
，为总运行时间，对L( 0 )取对数并求导，求解似然 方程，得到0和A的极大似然点估计为：
[0120] S0802,若部件寿命服从威布尔分布，则密度函数为：
[0122] 威布尔分布可以利用概率值很容易地推断出它的分布参数，n为比例参数，m为 形状参数。代入里程数据\，根据定时截尾样本数据，得到该样本的可靠度似然函数：
[0124]其中，L是似然函数，S,是似然函数参数，使用用参数估计迭代数值求解算法求解 以下方程组：
[0126] 得到可靠度函数为：
[0128] 其中，参数估计迭代数值求解算法具体步骤为：
[0129] S0802a，记：
[0131] S0802b，选择初值mQ= 1。
[0133] S0802d，重复步骤S0802c，直到 |mk+1-mk| < 0. 0001，得到m=mk+1。
[0134] S0802e，根据下式计算n:
[0136] S0803,若部件寿命服从对数正态分布，则密度函数为：
[0138] 定时截尾时可靠度似然函数为：
[0140] 设ZQ=(InxQ-li)/〇，标准正态分布函数①（_ZQ) = 1-①（Z。)，并且记①（ZQ)为 标准正态分布密度函数，则似然方程为：
[0142] 使用参数估计近似数值求解算法求解上述方程组，即可得到参数y、〇的极大似 然估计，从而得到可靠度函数。
[0143] 其中，参数估计近似数值求解算法具体包括以下步骤：
[0144] S0803a，选择初值yQ，〇0。
[0147]S〇8〇3b，计算
[0149] 其中U(zs) = [V(zs)_zs]。
[0150] S0803c，重复步骤（8-3b)，直到 |yk+1-yk | < 0? 0001 且
，得 到最终的y，〇极大似然估计。
[0151] S0804,若部件寿命服从正态分布，正态分布的求解为：将对数正态分布密度函数 f(t)中的lnx替换为X，其余步骤与对数正态分布的求解方法相同。
[0152] S0805,若部件寿命分布无规律，令收集到的部件数据为Xpx2,Xy. .Xp当Xi是故 障数据时，令\=1;当、是右截尾数据时，令S,= 〇,这样数据可记为（Xi，\)，1 = 1，2, . . .，r，将这些\按从小到大排列，对应的里程数据为x#x2彡...彡x^
[0153] 可靠度函数的乘积限估计为：
[0155] 平均寿命估计如下：
[0157] 在获得上述单个部件平均故障间隔时间和可靠度函数后，可根据动车组编组情 况，计算车组级部件平均故障间隔时间和可靠度。
[0158] 本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发 明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的 普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各 种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。
【主权项】
1. 一种动车组关键部件运行可靠性评估方法，其特征在于，包括以下步骤： (1) 构建标准化动车组子系统结构树，该动车组子系统结构树的最底层为动车组运用 现场最小不可拆分部件，每棵树所有层次的部件名称和数量相同； (2) 根据动车组实际上线运用情况，按投入运用时间划分车组批次，得到有η件产品投 入使用； (3) 记录同一批次动车组的标准化部件现场发生故障的时间点h，i = 1，2, 3. ..，r，r 为故障个数； (4) 根据动车组现场动态维修记录的离散时间与运行里程信息的关系，计算各个时间 点h上对应的动车组运行里程X i，i = 1，2, 3. . .，r ; (5) 计算Xi的均值Θ、方差S、二阶矩μ 2、三阶矩μ3、四阶矩μ4、偏度cs、峰度(；和对 数化样本的偏度Cs'、峰度C/ ; (6) 计算部件故障样本的平均失效率，得到失效率表； (7) 选择部件寿命分布模型，确定部件可靠度函数估计方法； (8) 确定单个部件平均故障间隔时间和可靠度函数； (9) 确定车组级部件平均故障间隔时间和可靠度函数； (10) 利用可靠度函数对各个部件进行可靠性评估。2. 根据权利要求1所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤（2)中，从第一批动车投入运行开始，同一批次的所有车组投入运用时间相差不 超过一年。3. 根据权利要求1所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤（4)中，采用分段线性插值的方式来获得4. 根据权利要求1所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤（5)具体包括以下步骤： (5-1)根据下列公式计算Xi的均值Θ、方差S、二阶矩μ 2、三阶矩μ 3、四阶矩μ 4、偏度 cs、峰度(；：Ce= y 4/s ; (5-2)将Xi作对数变换，重复步骤（5-1)得到对数化样本的偏度C s'、峰度(V。5. 根据权利要求1所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤（6)具体包括以下步骤： (6-1)将动车组总运行里程分成k个里程区间Λ X，计算每个里程区间长度： Δ X = (La-Sm) /k 其中，k = 1+3. 31gn，LaS X i的最大值，S X i的最小值； (6-2)计算部件故障样本的平均失效率无(Δχ,):其中，Λ A为第i个里程区间Λ X i内的失效率频数，i = 1，2, 3. ..，k，r 为进入第 i个时间区间时的样本数，rs，H= r-r η，:^指进入第i个时间区间的累积失效率。6. 根据权利要求1所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤（7)具体包括以下步骤： (7-1)判断部件故障数据个数，若大于10,则进入步骤（7-2);否则返回步骤（7-1); (7-2) I若I Θ-S I〈（ Θ+s)/5,则部件寿命服从指数分布，否则进入步骤（7-3); (7-3)若|CS|〈0.5且|Ce-3 I〈0.5,则部件寿命服从正态分布，否则进入步骤（7-4); (7-4)若|CS' I〈0.5且|C/-3 I〈0.5,则部件寿命服从对数正态分布，否则进入步骤 (7-5)； (7-5)根据步骤（6)中得到的平均失效率表，计算Λ λ i:若夂对于i = 1，2，... k均成立，则部件寿命服从指数分布，否则进入步 骤（7-6); (7-6)对步骤（7-5)中Λ λ i，记Λ λ，〇的个数为a，Λ λ'〇个数为b，若a/b彡3/4， 则部件寿命服从威布尔分布，否则部件寿命分布无规律，可靠性评估采用非参数方法。7. 根据权利要求6所述的一种动车组关键部件运行可靠性评估方法，其特征在于，所 述的步骤⑶具体为： (8-1)若部件寿命服从指数分布，则将动车组部件寿命数据作为无替换定时截尾情形， 即有η件产品投入使用，到规定的时间Xtl进行数据收集，依照时间的先后记录截止到规定 时间的失效时间获得对应的走行里程数据，其中： X1^ X2^... ^ X X 0? r<n 根据定时截尾样本数据，获得该样本的可靠度似然函数L( Θ ):其中，，为总运行里程，对L(0)取对数并求导，求解似然方 程，得到Θ和λ的极大似然点估计为：(8-2)若部件寿命服从威布尔分布，则密度函数为：故障样本走行里程为Xi，其中截尾时间Χ(ι，引入似然函数参数δ i，根据定时截尾样本 数据，得到该样本的可靠度似然函数：使用参数估计迭代数值求解算法求解以下方程组：得到可靠度函数为：(8-3)若部件寿命服从对数正态分布，则密度函数为：定时截尾时可靠度似然函数为：设Ztl= (Inxtl-μ)/s，标准正态分布函数Φ (-?) = l-φ (Ztl),并且记Φ (Ztl)为标准正 态分布密度函数，则似然方程为：使用参数估计近似数值求解算法求解上述方程组，即可得到参数μ、S的极大似然估 计，从而得到可靠度函数； (8-4)若部件寿命服从正态分布，正态分布的求解为：将对数正态分布密度函数f(t) 中的Inx替换为X，其余步骤与对数正态分布的求解方法相同； (8-5)若部件寿命分布无规律，令收集到的部件数据为Xl，x2, Xi. 当Xi是故障数据 时，令δ i= 1 ;当X肩右截尾数据时，令δ F 〇,将数据记为（X D δ山i = 1，2, . . .，r，将 这些\按从小到大排列，得X # X 2彡..·彡X y 可靠度函数的乘积限估计为：平均寿命估计如下：
【专利摘要】本发明涉及轨道车辆安全领域，具体涉及一种动车组子系统关键部件运用可靠性评估方法，主要包括：建立动车组子系统结构树并进行现场数据预处理，获得故障间隔运行里程数据；计算所得样本的均值、方差、二阶矩、三阶矩、四阶矩、偏度Cs、峰度Ce、平均失效率以及对数化样本的偏度Cs’、峰度Ce’；确定部件寿命分布和参数估计算法；计算平均故障间隔时间和可靠度。本发明有效处理动车组实际运行数据并对动车组部件进行运行可靠性估计，解决了动车组上线运用可靠性分析缺乏的问题，同时，本发明的基于现场数据的动车组关键部件运行可靠性评估方法也适用于其他轨道车辆部件的运行可靠性分析。
【IPC分类】G06F19/00
【公开号】CN104899423
【申请号】CN201510225640
【发明人】罗敏, 马哲, 陆正刚, 李道军, 李 浩
【申请人】同济大学
【公开日】2015年9月9日
【申请日】2015年5月6日
最新回复(0)