装备可靠性工程专栏

基于历史数据修正系数的可靠性分配方法

徐建坤1,魏盛兴1,高 剑2,孙 博1,邹天刚2,桂 鹏2

(1.北京航空航天大学 可靠性与系统工程学院,北京 100191;2.中国北方车辆研究所,北京 100072)

摘要:针对传统可靠性分配方法过分依赖历史数据或专家经验,以及对复杂机械系统进行可靠性分配时与实际生产制造过程不匹配两类问题,提出一种基于历史数据修正系数的可靠性分配方法。在改进的方法中,在分配流程上引入了系统部件隶属性矩阵,并在分配方法上提出了考虑新研产品特性的历史数据修正系数。为了证明改进方法的有效性,以某型液力机械综合传动系统为研究对象来进行可靠性指标分配。结果表明,改进的分配方法能够充分考虑历史故障数据和新研产品特性两类分配信息,同时能够有效避免传统方法中“自上而下”的分配流程与生产过程相脱离的情况。

关键词:复杂机械系统;可靠性分配;历史数据;修正系数;隶属性矩阵

0 引言

可靠性是指产品在规定条件下和规定时间内完成规定功能的能力[1]。随着各种机械设备的结构层次和功能划分日趋复杂,执行任务的环境条件愈加严苛,对复杂设备系统的可靠性要求也越来越高。可靠性分配作为可靠性设计的重要组成部分[2],是将可靠性定量指标按要求分配到各分系统、部件/组件、零件/元器件的自上而下的分解过程,目标是综合权衡有限的资源实现整个系统的可靠性目标,确保每一个基本单元的可靠性都满足可靠性设计要求[3-4]

目前,产品设计人员的观念从以往重视产品性能逐渐向注重产品的可靠性维修性方面发展[5-7],因此可靠性分配在现代复杂机械设备的正向研发设计过程中有着至关重要的指导意义,其分配结果直接明确了各组成部分的可靠性要求,并作为可靠性鉴定考核试验与评估的依据,分配结果的准确与否对产品内在可靠性产生直接影响[8]。在实际工程上常采用的可靠性分配方法有等分配法、相似产品法、评分分配法、AGREE分配法等传统分配方法[9-11]。等分配法假设系统各组成单元的可靠性水平相同,分配结果合理性较差。相似产品法利用已有相似系统的可靠性数据对新系统进行可靠性指标分配,但其过分依赖历史故障数据,且未考虑新研系统的技术特性。评分分配法根据专家经验对系统组成单元的技术水平、工作时间、环境条件等因素评分,根据评分结果进行可靠性分配,但该方法受专家主观影响大。AGREE分配法根据系统各组成单元的复杂度、重要度和工作时间进行可靠性分配,但对于复杂机械系统往往难以提前确定各个零部件的工作时间信息。

因此,传统的可靠性分配方法实质上都会过分的依赖历史故障数据或者专家经验的决策而具有应用上的局限性[12],在分配流程上,由于在生产制造过程中并不直接生产分系统,对于分系统划分有部件结构交叉的复杂系统,传统的自上而下的分配过程难以进行统一协调的可靠性指标分配,造成可靠性分配与实际生产过程相脱离的情况。为了克服以上传统分配方法的不足,同时更好的适应实际工程上的需要,本文中综合考虑历史故障数据和新研产品特点,提出一种基于历史数据修正系数的可靠性分配方法,并在分配流程上引入部件隶属性矩阵,为复杂机械产品系统实现科学合理的可靠性分配提供一种新的方法依据。

1 复杂机械系统可靠性分配方法流程

1.1 典型复杂机械系统功能-结构层次划分

现代典型复杂机械系统的功能层次一般划分为:系统级、分系统级、部件级和元件级,在进行可靠性指标分配时需要依次将系统的顶层指标分配到分系统、部件及元件。但复杂机械系统的功能层次划分的界限并不明确,按照功能划分的各个分系统会在结构层次经常会出现交叉,按照传统的自上而下的可靠性分配方法会造成部件级和元件级被重复分配可靠性指标的情况,此时传统方法的由上至下的可靠性分配流程并不适用。

将本文中提出的基于历史数据修正系数的可靠性分配方法适用对象划分为四级,即系统级、分系统级、部件级和元件级,其功能-结构层次如图1所示,其中分系统级按照具体功能划分。

图1 复杂系统功能-结构层次图

Fig.1 Complex system function-structure hierarchy diagram

1.2 基于隶属性矩阵的可靠性分配流程

典型复杂机械系统的结构分为系统级、分系统级、部件级和元件级,因此其可靠性分配工作一般分3个步骤进行,即将系统顶层的可靠性指标分别分配到分系统、部件和元件,以上3个可靠性分配步骤由于分配对象的不同分配方法也有所差异。由于元件在使用寿命周期中难以得到其历史故障数据,因此本文中提出的基于历史数据修正系数的可靠性分配方法仅探讨将系统基本可靠性指标分配到分系统和部件。

显然,由于各个分系统按照具体功能划分,因此会出现同一个部件属于不同分系统的情况,此时传统方法的由上至下的可靠性分配流程并不适用。为了解决这个问题,提出部件隶属性矩阵的概念,通过隶属性矩阵判断该分系统包括的具体部件。在进行可靠性分配时,首先利用基于历史数据修正系数的可靠性分配方法将系统顶层指标分配到各个部件,再通过隶属性矩阵由部件可靠性指标反推计算得到分系统可靠性指标。具体分配步骤如图2所示。

图2 可靠性分配计算流程图

Fig.2 Flow chart of reliability allocation calculation

1) 确定A系统基本可靠性指标失效率λA

2) 确定各个部件所属分系统,填写部件隶属性矩阵表(如表 1所示),得到隶属性矩阵。其中部件Ci属于该分系统填1,反之填0。隶属性矩阵用M表示:

(1)

表1 部件隶属性矩阵表

Table 1 Component attribute matrix table

B1B2…BnC1b11b12…b1nC2b21b22…b2n︙︙︙…︙Cnbm1bm2…bmn

3) 调研收集历史产品系统各个部件的历史故障数据,历史故障数据集记作I:

I=(I1,I2,…,In)

(2)

4) 根据新研机械系统特性,确定评分特性k个,组织专家根据评分特性给所有部件打分,部件Ci打分结果记作Si:

Si=(Si1,Si2,…,Sik)

(3)

5) 使用基于历史数据修正系数的可靠性分配方法将顶层指标λA分配到各个部件,分配结果记作λci:

λci=(λc1,λc2,…,λcm)

(4)

6) 计算分系统可靠性指标,分配结果记作λBi:

(5)

2 考虑两类信息的可靠性分配方法框架

根据以上分析,针对传统可靠性分配方法的不足,结合典型复杂机械系统的功能-结构特点,提出一种能够综合考虑历史故障数据和新研产品特点两类信息的可靠性分配方法。

2.1 基于相似产品历史数据的可靠性分配

一般来说,新研机械类产品与历史产品具有较高的相似性,因此,历史产品的故障数据对于指导新研产品的可靠性设计具有重要意义,基于相似产品故障数据按比例分配可靠性指标的步骤如下:

1) 定义复杂机械产品的系统总失效率λA;

2) 收集历史产品故障数据,统计各部件故障数,记作I=(I1,I2,…,In);

3) 计算各部件历史故障比例Ki:

(6)

式中,Ki为第i个部件的故障比例; Ii为第i个部件的故障树。

4) 计算新研机械产品各部件的失效率λci1:

λci1=λA×Ki

(7)

2.2 考虑新研产品特性的可靠性分配结果修正

根据3.1节中传统方法进行可靠性分配后,历史未出现故障的部件的失效率为0(即可靠度为1),实际上新研产品的各个部件都会受到技术方案、制造技术成熟度等多方面的影响,其失效率往往大于0(即可靠度小于1)。因此,以相似产品历史故障数据为基础进行的初次分配结果并不准确,需要对其进行修正。

考虑机械产品的具体可靠性特点,如产品复杂度、环境条件、载荷强调与频次、制造技术成熟程度等,利用专家领域知识根据新研产品特点对各部件进行多维度的综合定量评分,根据获得的评分结果对第一部分可靠性分配结果进行二次调整与优化,使得分配结果既能利用相似产品历史测试获得真实故障数据,又能兼顾新研产品的具体特点。其分配步骤如下:

1) 根据新研产品特性,计算各部件新特性评分

根据具体机械产品确定q个新研产品特性,记作Q=(Q1,Q2,…,Qq),由专家组对产品特性打分,每个特性得分取专家组的打分结果均值,结果记作Xij(j=1,2,…,q-1,q),分值取值范围1~10分,再将6个维度的得分相乘,得到该部件的专家打分结果,分值高代表许用失效率高,可靠性差。记第i个部件新特性评分为Wi,则:

(8)

式中,Wi为第i个部件的新特性评分;Xij为第i个部件第j个新特性评分。

2) 根据新特性专家评分结果,计算各部件失效率λci2

有各部件评分所占比例αi:

(9)

则各部件失效率λci2:

λci2=λA×αi

(10)

3) 根据新研产品特性,确定各部件修正系数

根据相似机械产品的相似程度和历史故障数据的可信程度,结合新研机械产品与可靠性相关的影响因素情况,确定依据历史故障数和新研产品特性评分两者之间的调节系数。定义β为新研产品特性评分调节系数,则历史故障数确定的失效率调节系数为1-β

4) 综合考虑历史故障数据和新研机械产品特点评分结果,计算最终的失效率λci

λci=λci1×(1-β)+λci2×β

(11)

综合以上步骤,基于历史数据修正系数的可靠性分配方法的流程如图3所示。

图3 基于历史数据修正系数的可靠性分配方法流程

Fig.3 Flow of reliability allocation method based on historical data and correction coefficient

3 算例验证

为验证以上方法的有效性,现以某液力机械综合传动系统为研究对象,采用本文中提出的基于历史数据修正系数的可靠性分配方法对该传动系统方案设计阶段的基本可靠性指标进行分配。

3.1 某液力机械综合传动系统基本信息

该传动系统是集机电液为一体的复杂系统,可以提供直驶变速、转向、制动、冷却、驱动、启动等功能。按照系统-分系统-部件对其进行层级划分,则该机械综合传动系统有直驶变速系统、转向系统、中心转向系统、制动系统、冷却系统、驱动系统和启动系统7个分系统,在部件层级有箱体、中间支架、前传动总成、风扇传动总成、离合器总成、油泵组、供油系统、液压操纵系统、操纵电控系统、变矩器支架、行星变速机构、辅助传动、左汇流排、右汇流排、联体泵马达、测试系统、液力减速器、减速器控制阀、液力变矩器和传动总体共20个部件。

3.2 基本可靠性指标与新研产品特性

1) 基本可靠性指标。根据机械综合传动系统执行任务特点,选择平均故障间隔里程(mean kilometers between failures,MKBF)和失效率λA作为其基本可靠性参数,二者之间有如下转换关系:

(12)

式中:λA表示以平均故障间隔里程为参数的失效率; λAt表示以平均故障间隔时间为参数的失效率;MTBF表示平均故障间隔时间。

参考历史产品可靠性指标数据,可初步确定系统级的MKBF指标规定值为4 000 km。根据传统的可靠性分配方法经验、参考近年国内外针对相似产品的可靠性指标分配实例,对MKBF指标采用阶段式的提升方式,即乘以提升系数1.25[13-15],得到可靠性指标目标值,以提高该项指标的裕度。

表2 传动系统基本可靠性指标

Table 2 Basic reliability index of transmission system

基本可靠性指标MKBF/kmλA规定值4 0000.000 25目标值5 0000.000 20

2) 新特性评分。选择能够准确评价机械评分特性是得到准确可靠性分配结果的前提,根据新研机械综合传动系统的具体特点,考虑各评分特性的全面性和相互独立性,确定了评分特性集为:产品复杂程度,环境条件,载荷频次,设计技术成熟程度,制造技术成熟程度,维修难易程度。由专家组对机械综合传动系统根据上述6个评分特性进行评分,评分范围为0~10,评分原则如表3所示。

表3 传动系统新特性评分准则

Table 3 Scoring criteria for new features of transmission system

评分特性评分原则产品复杂程度产品越复杂,评分越高环境条件环境条件越恶劣,评分越高载荷频次载荷频次越多,评分越高设计技术成熟程度成熟程度越低,评分越高制造技术成熟程度成熟程度越低,评分越高维修难易程度维修难度越大,评分越高

3.3 可靠性分配计算及结果

根据图 3中的可靠性分配计算流程,确定机械综合传动系统可靠性分配的顶层指标和相关参数具体值如下:

1) 系统顶层指标:λA=2×10-4,MKBF=5 000 km

2) 部件隶属性矩阵

3) 各部件历史故障数据:

I=(9,2,3,2,2,16,5,9,7,5,13,1,5,3,9,10,3,5,9,19)

4) 各部件专家打分结果

5) 调节系数β=0.3。

计算得到各分系统和各部件的可靠性分配结果分别如表4、表5所示。

表4 分系统可靠性分配结果

Table 4 Subsystem reliability allocation results

序号分系统分系统可靠性指标失效率/km-1MKBF/km1直驶变速系统1.73E-045.78E+032转向系统1.76E-045.67E+033中心转向系统1.42E-047.04E+034制动系统8.46E-051.18E+045冷却系统1.07E-051.19E+046驱动系统8.40E-051.19E+047启动系统1.14E-048.80E+03

表5 部件可靠性分配结果

Table 5 Component reliability allocation results

序号部件部件可靠性指标失效率/km-1MKBF/km1箱体部件1.59E-056.29E+042中间支架4.50E-062.22E+053前传动总成5.44E-061.84E+054风扇传动总成3.48E-062.87E+055离合器总成3.24E-063.09E+056油泵组2.66E-053.76E+047供油系统6.55E-061.53E+058液压操纵系统9.74E-061.03E+059操纵电控系统7.67E-061.30E+0510变矩器支架6.78E-061.47E+0511行星变速机构1.58E-056.31E+0412辅助传动1.54E-066.47E+05

续表(表5)

序号部件部件可靠性指标失效率/km-1MKBF/km13左汇流排1.33E-057.51E+0414右汇流排3.75E-062.67E+0515联体泵马达1.14E-058.78E+0416测试系统1.22E-058.22E+0417液力减速器6.35E-061.58E+0518减速器控制阀1.06E-059.40E+0419液力变矩器1.17E-058.58E+0420传动总体2.35E-054.26E+04

4 结论

可靠性分配的过程本质上是一个全面权衡多种因素的工程决策过程,针对功能划分界限不清晰的复杂机械系统,传统的可靠性分配方法难以得出科学合理的结果。本文中给出的基于历史数据修正系数的可靠性分配方法,在计算方法上充分考虑了历史故障数据和新研产品特性两类信息,在分配流程上引入了部件隶属性矩阵,通过先计算部件可靠性指标再反推分系统可靠性指标的方法也更加贴合复杂机械系统的实际研制生产过程,为产品的可靠性设计提供了参考。

参考文献:

[1] KIM K O,YANG Y,ZUO M J.A new reliability allocation weight for reducing the occurrence of severe failure effects[J].Reliability Engineering &System Safety,2013,117(3):81-88.

[2] WANG Z,HAN Y,LI W,et al.An information fusion reliability allocation method considering limited sample size circumstances[J].Quality and Reliability Engineering International,2021,37(5):1863-1874.

[3] WEI L H,YAO C J,WANG H L.Research on reliability allocation method of rv reducer system[C]//IOP Conference Series:Earth and Environmental Science.IOP Publishing,2019,237(5):052046-052054.

[4] CATELANI M,CIANI L,GUIDI G,et al.Reliability Allocation:an iterative approach for complex systems[C]//2021 IEEE International Symposium on Systems Engineering (ISSE).IEEE,2021:1-6.

[5] 张义民.机械可靠性设计的内涵与递进[J].机械工程学报,2010(14):22-29.

ZHANG Yiming.Connotation and development of mechanical reliability-based design[J].Journal of Mechanical Engineering,2010(14):22-29.

[6] 严宇飞,李方义.基于蒙特卡洛重要抽样法的结构时变可靠性分析[J].重庆理工大学学报(自然科学),2021,35(6):104-112.

YAN Yufei,LI Fangyi.Time-variant reliability analysis of structures based on important sampling Monte Carlo method[J].Journal of Chongqing University of Technology(Natural Science),2021,35(6):104-112.

[7] 顾士林.机械工程的可靠性优化设计分析[J].内燃机与配件,2019(24):216-217.

GU Shilin.Reliability optimization design analysis of mechanical engineering[J].Internal Combustion Engine &Parts,2019(24):216-217.

[8] 张义民,黄贤振.机械产品研发的可靠性规范[J].中国机械工程,2010(23):2773-2785.

ZHANG Yiming,HUAng xianzhen.reliability guiding principles for research and development of mechanical products[J].China Mechanical Engineering,2010(23):2773-2785.

[9] 程强,孙东洋,刘志峰,等.基于模糊数和AHP的数控机床可靠性分配方法[J].组合机床与自动化加工技术,2020(2):34-37,42.

CHENG Qiang,SUN Dongyang,LIU Zhifeng,et al.Reliability assignment method of cnc machine tool based on fuzzy number and AHP[J].Modular Machine Tool &Automatic Manufacturing Technique,2020(2):34-37,42.

[10] 李博远,胡丽琴,陈珊琦,等.基于故障树和层次分析的可靠性分配方法[J].安全与环境工程,2015,22(1):117-120.

LI Boyuan,HU Liqin,CHEN Shanqi,et al.Reliability assignment method based on fault tree and hierarchy analysis[J].Safety and Environmental Engineering,2015,22(1):117-120.

[11] 石闯,陈永冰.基于评分法的某新型船舶操舵控制系统基本可靠性分配[J].舰船电子工程,2021,41(3):115-118.

SHI Chuang,CHEN Yongbing.Basic reliability assignment of a new ship steering control system based on scoring method[J].Ship Electronic Engineering,2021,41(3):115-118.

[12] 程力旻,徐巍,张文俊.基于免疫遗传优化的可靠性分配方法[J].舰船电子工程,2022,42(2):152-156.

CHEN Limin,XU Wei,Zhang Wenjun.Reliability assignment method based on immunogenetics optimization[J].Ship Electronic Engineering,2022,42(2):152-156.

[13] 侯同刚.空空导弹可靠性指标论证方法探讨[J].航空标准化与质量,2010(1):40-42.

HOU Tonggang.Discussion on air-to-air missile reliability index demonstration method[J].Aeronautic Standardization &Quality,2010(1):40-42.

[14] 李根成.生产定型可靠性鉴定试验检验下限研究[J].航空兵器,2019,26(2):81-85.

LI Gencheng.Research on the lower limit of inspection of production type reliability appraisal test[J].Aero Weaponry,2019,26(2):81-85.

[15] 李源.可靠性目标值的检验时机研究[J].航空兵器,2021,28(6):100-103.

LI Yuan.Research on the testing time of reliability target value[J].Aero Weaponry,2021,28(6):100-103.

Reliability allocation method based on correction coefficient and historical data

XU Jiankun1,WEI Shengxing1,GAO Jian2,SUN Bo1,ZOU Tiangang2,GUI Peng2

(1.School of Reliability and Systems Engineering,Beihang University,Beijing 100191,China; 2.China North Vehicle Research Institute,Beijing 100072,China)

AbstractIn order to solve the two problems of traditional reliability allocation methods that excessively rely on historical data or expert experience,and that the reliability allocation of complex mechanical systems does not match the actual manufacturing process,a reliability allocation method based on historical data correction coefficient is proposed.In the improved method,the subordinate attribute matrix of system components is introduced in the allocation process,and the historical data correction coefficient considering the characteristics of newly developed products is proposed in the allocation method.To demonstrate the effectiveness of the improvement method,a certain type of hydraulic mechanical comprehensive transmission system is taken as the research object to allocate reliability indicators.The results show that the improved allocation method can fully consider the historical fault data and the characteristics of newly developed products,and can effectively avoid the separation of the “top-down” allocation process from the production process in the traditional method.

Key wordscomplex mechanical system; reliability allocation; historical data; correction coefficient; membership matrix

本文引用格式:徐建坤,魏盛兴,高剑,等.基于历史数据修正系数的可靠性分配方法[J].兵器装备工程学报,2023,44(9):172-177.

Citation format:XU Jiankun,WEI Shengxing,GAO Jian,et al.Reliability allocation method based on correction coefficient and historical data[J].Journal of Ordnance Equipment Engineering,2023,44(9):172-177.

中图分类号:N945.17

文献标识码:A

文章编号:2096-2304(2023)09-0172-06

收稿日期:2022-10-28;

修回日期:2022-12-21

作者简介:徐建坤(1997—),男,硕士研究生,E-mail:ZY2114211@.buaa.edu.cn。

通信作者:孙博(1979—),男,博士,副研究员,博士生导师,E-mail:sunbo@.buaa.edu.cn。

doi:10.11809/bqzbgcxb2023.09.022

科学编辑 黄景德 博士 (吉林大学教授、博导)

责任编辑 唐定国