CN104992011B - 一种基于fmea分析的数控机床可靠性综合分配方法 - Google Patents
一种基于fmea分析的数控机床可靠性综合分配方法 Download PDFInfo
- Publication number
- CN104992011B CN104992011B CN201510358432.4A CN201510358432A CN104992011B CN 104992011 B CN104992011 B CN 104992011B CN 201510358432 A CN201510358432 A CN 201510358432A CN 104992011 B CN104992011 B CN 104992011B
- Authority
- CN
- China
- Prior art keywords
- subsystem
- failure
- reliability
- severity
- distribution
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000011058 failure modes and effects analysis Methods 0.000 title claims abstract description 40
- 238000004458 analytical method Methods 0.000 title claims abstract description 27
- 230000015572 biosynthetic process Effects 0.000 title 1
- 238000003786 synthesis reaction Methods 0.000 title 1
- 238000012546 transfer Methods 0.000 claims abstract description 13
- 238000011156 evaluation Methods 0.000 claims description 45
- 239000013598 vector Substances 0.000 claims description 40
- 239000011159 matrix material Substances 0.000 claims description 25
- 238000005192 partition Methods 0.000 claims description 25
- 238000006243 chemical reaction Methods 0.000 claims description 18
- 238000007476 Maximum Likelihood Methods 0.000 claims description 4
- 238000013077 scoring method Methods 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 5
- 241000039077 Copula Species 0.000 abstract description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 230000010354 integration Effects 0.000 abstract description 2
- 238000007796 conventional method Methods 0.000 abstract 1
- 230000001105 regulatory effect Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 4
- 125000006850 spacer group Chemical group 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000007789 sealing Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 210000004907 gland Anatomy 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Landscapes
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明一种基于FMEA分析的数控机床可靠性综合分配方法,属于机械系统可靠性分析领域,本发明基于FMEA三阶转换函数及Copula函数进行系统可靠性的综合分配,解决现有的机床可靠性分配方法的考虑因素单一及只能调节高严重度或失效频率的可靠性区分度等问题;该方法的转换函数形式简单,易于理解及操作,同时涵盖了传统指数转换函数的优势,并拓展了FMEA分析结果低区间及中区间的可靠性区分度调节能力;该方法同时考虑了多种分配指标的影响及故障相关性的作用,具有集成性好,便于工程应用的特点,通过该方法分析得到的子系统可靠性与传统方法相比更低,能够更好地降低制造及加工成本要求,具有十分重要的工程价值。
Description
技术领域
本发明属于机械系统可靠性分析领域,具体涉及一种基于FMEA分析的数控机床可靠性综合分配方法。
背景技术
机械可靠性分配是根据一定的原则和方法,将整机系统可靠性指标自上而下逐级分配到下属各级产品,最终达到整机系统可靠性指标的过程,机械系统的可靠性在一定程度上取决于设计初期可靠性指标分配过程的合理性;
数控车床作为制造业的主要设备,其可靠性程度对制造业的发展起到关键作用;由于数控车床结构复杂且故障模式繁多,其设计初期的可靠性分配结果常出现不合理或可信度不足的现象,进而影响机床的使用效率及寿命;因此,讨论如何合理地对数控车床进行分配或再分配十分重要;
数控机床可以看作是一种典型的串联系统,针对串联系统的可靠性分配,目前主要有平均分配法、基于FMEA的直接分配法、基于FMEA的指数转换函数分配法、综合分配法、基于模糊理论的可靠性分配方法等;然而,不同的理论均有其各自的优势与不足;例如,相比于基于FMEA的平均分配法,基于FMEA的指数转换函数分配法虽然增强了高严重度及高失效频率的可靠性区分度,但是在工程实际中,大部分子系统的严重度及失效频率均分布在中区间内(FMEA评价值为5-8),该方法并不能很好地适用。
发明内容
针对现有技术的不足,本发明提出一种基于FMEA分析的数控机床可靠性综合分配方法,以达到根据工程实际需要调节可靠性分配所关注的严重度区间,调节子系统可靠性的区分度,在子系统分配较低可靠度的同时使整个系统达到要求的可靠度,以达到降低加工制造成本的目的。
一种基于FMEA分析的数控机床可靠性综合分配方法,包括以下步骤:
步骤1、设置目标数控机床的目标可靠度,并确定目标数控机床所包括的子系统,及每个子系统的可靠性评价指标和失效模式;
步骤2、判断可靠性评价指标中是否包含子系统相关失效严重度,即是否考需要考虑故障相关性的影响,若是,则执行步骤3,否则执行步骤4;
步骤3、根据目标数控机床子系统的失效概率,获得子系统相关失效严重度;
步骤4、结合实际情况,采用FMEA分析方法确定每个子系统中不同失效模式的FMEA评价值,包括不同失效模式的失效严重度和不同失效模式的失效频率,根据上述失效严重度和失效频率,进一步获得所有子系统的相对失效严重度平均值和相对失效频率平均值;
步骤5、将FMEA分析方法中的严重程度和失效程度进行分区,包括低、中、高三区;
步骤6、判断每个子系统的相对失效严重度平均值和相对失效频率平均值所属分区,并设定该分区内子系统的可靠性分区度;
步骤7、根据不同分区所设定的子系统的可靠性分区度和三阶转换函数,将每个子系统的严重度和失效率有线性关系转换为非线性关系;
步骤8、根据转换后获得的每个子系统失效严重度和失效频率,及其他可靠性评价指标,获得每个子系统在不同可靠性评价指标下的分配系数和不同可靠性评价指标的分配权重,进而获得子系统分配向量;
步骤9、根据获得的子系统分配向量,获得最终每个子系统的可靠度,具体如下:
当可靠性评价指标中不包括子系统相关失效严重度时,具体如下:
步骤9-1、将步骤1所设定的目标可靠度转换为目标失效率,根据子系统分配向量和目标失效率,获得每个子系统的失效率;
步骤9-2、将每个子系统的失效率转换为可靠度,即获得最终每个子系统的可靠度;
当可靠性评价指标中包括子系统相关失效严重度时,具体如下:
步骤9-a、根据子系统分配向量,获得子系统失效率比值;
步骤9-b、根据步骤1所设定的目标可靠度、子系统的相关系数、子系统的故障相关关系和子系统失效率比值,获得每个子系统的失效率;
步骤9-c、将每个子系统的失效率转换为可靠度,即获得最终每个子系统的可靠度;
步骤10、根据所获得的每个子系统的可靠度,设计目标机床。
步骤3所述的获得子系统相关失效严重度,具体步骤如下:
步骤3-1、统计目标数控机床子系统的失效概率;
步骤3-2、采用极大似然估计法获得任意两个子系统的相关系数;
步骤3-3、建立目标数控机床子系统的相关系数矩阵;
步骤3-4、根据相关系数矩阵计算获得子系统相关失效严重度。
步骤4所述的根据上述失效严重度和失效频率,进一步获得所有子系统的相对失效严重度平均值和相对失效频率平均值,具体为:在每个子系统不同失效模式的失效严重度中选择最大值作为相对失效严重度,再获得所有子系统的相对失效严重度平均值;在每个子系统不同失效模式的失效频率中选择最大值作为相对失效频率,再获得所有子系统的相对失效频率平均值。
步骤5所述的低区取值范围为1~3级,中区取值范围为4~7级,高区取值范围为8~10级。
步骤7所述的根据不同分区所设定的子系统的可靠性分区度和三阶转换函数,将每个子系统的严重度和失效率有线性关系转换为非线性关系,具体如下:
三阶转换函数公式如下:
其中,E表示FMEA严重程度或失效程度,ξ表示相对失效严重度平均值或相对失效频率平均值,a0、c0和c1均为待定系数;
三阶转换函数与可靠性分区度关系式如下:
其中,ΛS表示可靠性分区度;
将公式(2)代入至公式(1)中求解a0、c0和c1系数的取值,即获得三阶转换函数即完成每个子系统的严重度和失效率从有线性关系向非线性关系的转换。
步骤8所述的获得子系统分配向量,具体步骤如下:
步骤8-1、将转换后获得的每个子系统失效严重度和失效频率作为两个可靠性评价指标,并根据其他可靠性评价指标,建立每个可靠性评价指标情况下,子系统两两之间的失效率分配相对值矩阵;
步骤8-2、根据失效率分配相对值矩阵,进一步获得每个子系统在不同可靠性评价指标下的相对分配系数矩阵;
步骤8-3、采用专家打分的方法对每个可靠性评价指标的权重进行打分,获得分配指标权重向量;
步骤8-4、将分配指标权重向量与相对分配系数矩阵相乘,获得子系统分配向量。
步骤9-b所述的根据步骤1所设定的目标可靠度、子系统的相关系数、子系统的故障相关关系和子系统失效率比值,获得每个子系统的失效率,具体公式如下:
其中,Rs *(t)表示设定的目标可靠度,m表示具有相关性的子系统个数,s表示相互独立的子系统个数,n表示子系统的个数,Ri(t)表示相关第i子系统的可靠度,λi表示第i个子系统失效率比值,ai表示分配向量的第i个元素,a1表示分配向量的第1个元素,t表示时间,表示不相关第子系统的可靠度,表示具有故障相关性的子系统序数;表示具有故障独立的子系统序数;θ表示相关系数。
本发明优点:
本发明提出一种基于FMEA分析的数控机床可靠性综合分配方法,基于FMEA三阶转换函数及Copula函数进行系统可靠性的综合分配,解决现有的机床可靠性分配方法的考虑因素单一及只能调节高严重度(或失效频率)的可靠性区分度等问题;该方法的转换函数形式简单,易于理解及操作,同时涵盖了传统指数转换函数的优势,并拓展了FMEA分析结果低区间及中区间的可靠性区分度调节能力;该方法同时考虑了多种分配指标的影响及故障相关性的作用,具有集成性好,便于工程应用的特点,通过该方法分析得到的子系统可靠性与传统方法相比更低,能够更好地降低制造及加工成本要求,具有十分重要的工程价值。
附图说明
图1为本发明一种实施例的系统组成图;
图2为本发明一种实施例的基于FMEA分析的数控机床可靠性综合分配方法流程图;
图3为本发明一种实施例的不考虑故障相关性时的操作流程图;
图4为本发明一种实施例的考虑故障相关性时的操作流程图;
图5为本发明一种实施例的数控车床主轴分系统的组成图。
具体实施方式
下面结合附图对本发明一种实施例做进一步说明。
本发明实施例,以HTC2050i型号机床为例,如图1所示,其主要涉及子系统为床身、XYZ轴、主轴系统、尾座、刀架系统;其中,主轴系统可细分为总支承、检测组件、前位支承、后位支承、前密封组件、后密封组件、传动组件、制动组件;其中,总支承包括主轴箱,检测组建包括光栅尺侧头、光栅尺、垫板、电缆,前位支承包括角接触球轴承、双列圆柱滚子轴承、圆螺母、隔套、压盖、套、套筒、前隔套,后位支承包括隔套、双列圆柱滚子轴承、封闭环、法兰盘,前密封组件包括密封盖、迷宫盖、O形橡胶密封圈、压盖、隔套,后密封组件包括密封环、法兰盘、防护罩,传动组件包括轴、同步带轮、锁紧螺母、键,制动组件包括刹车盘、液压直动夹钳、支架。
本发明实施例中,基于FMEA分析的数控机床可靠性综合分配方法,方法流程图如图2所示,本发明包括考需要考虑故障相关性的影响与不考需要考虑故障相关性的影响两种情况,图3和图4分别展示了两种情况的分配过程;
基于FMEA分析的数控机床可靠性综合分配方法包括以下步骤:
步骤1、设置目标数控机床的目标可靠度,并确定目标数控机床所包括的子系统,及每个子系统的可靠性评价指标和失效模式;
对数控机床系统进行分析,分别得到机床的子系统数目及名称,各子系统的失效模式及可靠性分配的评价指标,分别对应图3的分支②,③和①;
本发明实施例,以HTC2050i型号机床的故障相关性的机床主轴系统为例,如图5所示,本实施例中选取的7项分配指标,分别为:失效频率(k=1),失效严重度(k=2),复杂度(k=3),加工难度(k=4),工作环境(k=5),成本(k=6)及相关失效严重度(k=7);k表示分配指标;
本发明实施例,不同失效模式下的FMEA分析结果如表1所示:
表1 数控车床主轴系统FMEA分析
步骤2、判断可靠性评价指标中是否包含子系统相关失效严重度,即是否考需要考虑故障相关性的影响,若是,则执行步骤3,否则执行步骤4;
步骤3、根据目标数控机床子系统的失效概率,获得子系统相关失效严重度,具体步骤如下:
步骤3-1、统计目标数控机床子系统的失效概率;
步骤3-2、采用极大似然估计法获得任意两个子系统的相关系数;
本发明实施例中,如图4所示,考虑分支⑥,根据现场数据统计得到各个子系统同一故障模式下的失效率
步骤3-3、建立目标数控机床子系统的相关系数矩阵;
本发明实施例中,利用极大似然估计计算失效模式下任意两个子系统i,j的相关性系数
其中,ui表示子系统i的故障率,uj表示子系统j的故障率,选用GumbelCopula函数:
根据下式得到子系统相关系数:
计算任意两个子系统i,j的失效相关系数τij:
建立相关系数矩阵为:
本发明实施例中,考虑故障相关性时,对机床的监测故障信息进行统计,得到相关性矩阵为:
步骤3-4、根据相关系数矩阵计算获得子系统相关失效严重度;
计算相关失效严重度:
本发明实施例中,各子系统的相关严重度向量:
τi={2.47 1.86 1.25 1.98 2.16 1.73 1.48 1.17};
步骤4、结合实际情况,采用FMEA分析方法确定每个子系统中不同失效模式的FMEA评价值,包括不同失效模式的失效严重度和不同失效模式的失效频率,根据上述失效严重度和失效频率,进一步获得所有子系统的相对失效严重度平均值和相对失效频率平均值;
在每个子系统不同失效模式的失效严重度中选择最大值作为相对失效严重度或相对失效频率;
其中,Si为第i个子系统的相对严重度或失效频率;为第i个子系统第1个,第2个,......个失效模式的严重度值或失效频率;
再获得所有子系统的相对失效严重度平均值或相对失效频率平均值;
其中,ξ表示相对失效严重度平均值或相对失效频率平均值,n表示子系统的个数;
步骤5、将FMEA分析方法中的严重程度和失效程度进行分区,包括1~3级低区、4~7级中区、8~10级高区;
本发明实施例中,在低区内具体包括:无(1级)、很低(2级)、低(3级);中区具体包括:中等偏低(4级)、中等(5级)、中等偏高(6级),高区具体包括:高(7级)、很高(8级)、极高且预警(9级)、极高无预警(10级);
步骤6、判断每个子系统的相对失效严重度平均值和相对失效频率平均值所属分区,并设定该分区内子系统的可靠性分区度;
步骤7、根据不同分区所设定的子系统的可靠性分区度和三阶转换函数,将每个子系统的严重度和失效率有线性关系转换为非线性关系;
三阶转换函数的导数形式:
其中,待定转换函数;
三阶转换函数公式如下:
其中,E表示FMEA严重程度或失效程度E=1,2,...,10,ξ表示相对失效严重度平均值或相对失效频率平均值,a0、c0和c1均为待定系数;
三阶转换函数与可靠性分区度关系式如下:
其中,ΔS表示可靠性分区度;
将公式(2)代入至公式(1)中求解a0、c0和c1系数的取值,即获得三阶转换函数即完成每个子系统的严重度和失效率从有线性关系向非线性关系的转换;
步骤8、根据转换后获得的每个子系统失效严重度和失效频率,及其他可靠性评价指标,获得每个子系统在不同可靠性评价指标下的分配系数和不同可靠性评价指标的分配权重,进而获得子系统分配向量;
所述的获得子系统分配向量,具体步骤如下:
步骤8-1、将转换后获得的每个子系统失效严重度和失效频率作为两个可靠性评价指标,并根据其他可靠性评价指标,建立每个可靠性评价指标情况下,子系统两两之间的失效率分配相对值矩阵B(k);
本发明实施例中,失效率分配相对值矩阵B(k)如下:
其中,k表示分配指标,如失效频率为k=1,严重度为k=2等,为第k个分配指标下,第i个子系统与第j个子系统的失效率分配相对值,i,j=1,2,...n;
对于失效频率为k=1与严重度为k=2,有
其中,Δi表示子系统i的相对失效频率或失效严重度;Δj表示子系统j的相对失效频率或失效严重度;
本发明实施例中,考虑故障相关性时,k=3至k=6的分配指标值如表2所示;
表2 子系统k=3至k=6指标值
k=7的分配指标值可根据τi={2.47 1.86 1.25 1.98 2.16 1.73 1.48 1.17}计算得到,其中,k=3,4,...,7时,有
其中,Δi,Δj分别为子系统i,j的复杂度,加工难度,工作环境,成本或相关失效严重度;
步骤8-2、根据失效率分配相对值矩阵,进一步获得每个子系统在不同可靠性评价指标下的相对分配系数矩阵Γ;
本发明实施例中,系数矩阵B(k)的每一列代表与该列序号对应的子系统与其它所有子系统的失效率分配相对值;取第i列值的平均数,代表第i个子系统在第k个分配指标下的相对分配系数,写为:
考虑系统共有m个分配指标,根据式(9),可得到相对分配系数矩阵Γ:
本发明实施例中,不考虑故障相关性时,获得相对分配系数矩阵Γ为:
本发明实施例中,考虑故障相关性时,获得相对分配系数矩阵Γ为:
步骤8-3、采用专家打分的方法对每个可靠性评价指标的权重进行打分,获得分配指标权重向量W;
本发明实施例中,如图3所示,引入分支④,对各分配指标的权重进行评分,并得到分配指标权重向量;权重向量可通过专家打分方法获得:
W={w1,w2,…,wm} (17)
其中,
其中,wkr表示第r个专家对第k个分配指标的重要程度的评价分数,N表示参与评分的专家人数,w′j表示第j个分配指标的重要程度的评价分数。
本发明实施例中,不考虑故障相关性时,分配指标权重向量为:
W={0.22 0.11 0.15 0.18 0.16 0.18};
本发明实施例中,考虑故障相关性时,分配指标权重向量为:
W={0.22 0.11 0.15 0.18 0.16 0.18 0.11};
步骤8-4、将分配指标权重向量W与相对分配系数矩阵Γ相乘,获得子系统分配向量A;
A={a1,a2,…,an}=W·Γ (19)
其中,a1,a2,…,an为分配向量中的元素。
本发明实施例中,不考虑故障相关性时,子系统分配向量A为:
A={0.482 0.540 0.559 0.551 0.422 0.379 0.586 0.480};
本发明实施例中,考虑故障相关性时,子系统分配向量A为:
A={0.528 0.593 0.623 0.602 0.471 0.434 0.645 0.545}
步骤9、根据获得的子系统分配向量,获得最终每个子系统的可靠度,具体如下:
当可靠性评价指标中不包括子系统相关失效严重度时,具体如下:
步骤9-1、将步骤1所设定的目标可靠度转换为目标失效率,根据子系统分配向量和目标失效率,获得每个子系统的失效率;
本发明实施例中,如图3所示,考虑分支⑤,引入目标失效率λs *,计算子系统失效率公式如下:
本发明实施例中,不考虑故障相关性时,子系统的失效率为:
λ1∶λ2∶…∶λ8=1∶1.12∶1.16∶1.143∶0.876∶0.786∶1.216∶0.996;
其中,λi表示第i个子系统的失效率,ai表示分配向量中的第i个元素,aj表示分配向量中的第j个元素,λs *表示目标可靠度转换为目标失效率;
步骤9-2、将每个子系统的失效率转换为可靠度,即获得最终每个子系统的可靠度;
根据可靠性计算公式:
其中,Ri表示第i个子系统的可靠度;
将公式(21)代入式(20)的计算结果,得到各子系统的可靠性分配结果;
当可靠性评价指标中包括子系统相关失效严重度时,具体如下:
步骤9-a、根据子系统分配向量,获得子系统失效率比值;
本发明实施例中,考虑故障相关性时,子系统的失效率为:
λ1∶λ2∶…∶λ8=1∶1.124∶1.181∶1.141∶0.893∶0.823∶1.223∶1.033;
步骤9-b、根据步骤1所设定的目标可靠度、子系统的相关系数、子系统的故障相关关系和子系统失效率比值,获得每个子系统的失效率;
本发明实施例中,如图4所示,引入分支⑦,所述的根据步骤1所设定的目标可靠度、子系统的相关系数、子系统的故障相关关系和子系统失效率比值,获得每个子系统的失效率,具体公式如下:
其中,Rs *(t)表示设定的目标可靠度,m表示具有相关性的子系统个数,s表示相互独立的子系统个数,n表示子系统的个数,Ri(t)表示相关第i子系统的可靠度,λi表示第i个子系统失效率比值,ai表示分配向量的第i个元素,a1表示分配向量的第1个元素,t表示时间,表示不相关第子系统的可靠度,表示具有故障相关性的子系统序数;表示具有故障独立的子系统序数;θ表示相关系数;
本发明实施例中,根据主轴故障频率分析,可知8个子系统的故障模式均直接或间接相关,近似取θ-1=0.3,可知:
步骤9-c、将每个子系统的失效率转换为可靠度,即获得最终每个子系统的可靠度;
t=500h时,子系统的可靠度为:
R1=0.9947,R2=0.9940,R3=0.9937,R4=0.9939
R5=0.9952,R6=0.9956,R7=0.9935,R8=0.9945
步骤10、根据所获得的每个子系统的可靠度,设计目标机床。
Claims (5)
1.一种基于FMEA分析的数控机床可靠性综合分配方法,其特征在于,包括以下步骤:
步骤1、设置目标数控机床的目标可靠度,并确定目标数控机床所包括的子系统,及每个子系统的可靠性评价指标和失效模式;
步骤2、判断可靠性评价指标中是否包含子系统相关失效严重度,即是否考需要考虑故障相关性的影响,若是,则执行步骤3,否则执行步骤4;
步骤3、根据目标数控机床子系统的失效概率,获得子系统相关失效严重度;
步骤3-1、统计目标数控机床子系统的失效概率;
步骤3-2、采用极大似然估计法获得任意两个子系统的相关系数;
步骤3-3、建立目标数控机床子系统的相关系数矩阵;
步骤3-4、根据相关系数矩阵计算获得子系统相关失效严重度;
步骤4、结合实际情况,采用FMEA分析方法确定每个子系统中不同失效模式的FMEA评价值,包括不同失效模式的失效严重度和不同失效模式的失效频率,根据上述失效严重度和失效频率,进一步获得所有子系统的相对失效严重度平均值和相对失效频率平均值;步骤5、将FMEA分析方法中的严重程度和失效程度进行分区,包括低、中、高三区;
步骤6、判断每个子系统的相对失效严重度平均值和相对失效频率平均值所属分区,并设定该分区内子系统的可靠性分区度;
步骤7、根据不同分区所设定的子系统的可靠性分区度和三阶转换函数,将每个子系统的严重度和失效率有线性关系转换为非线性关系;
步骤8、根据转换后获得的每个子系统失效严重度和失效频率,及其他可靠性评价指标,获得每个子系统在不同可靠性评价指标下的分配系数和不同可靠性评价指标的分配权重,进而获得子系统分配向量;
步骤9、根据获得的子系统分配向量,获得最终每个子系统的可靠度,具体如下:
当可靠性评价指标中不包括子系统相关失效严重度时,具体如下:
步骤9-1、将步骤1所设定的目标可靠度转换为目标失效率,根据子系统分配向量和目标失效率,获得每个子系统的失效率;
步骤9-2、将每个子系统的失效率转换为可靠度,即获得最终每个子系统的可靠度;
当可靠性评价指标中包括子系统相关失效严重度时,具体如下:
步骤9-a、根据子系统分配向量,获得子系统失效率比值;
步骤9-b、根据步骤1所设定的目标可靠度、子系统的相关系数、子系统的故障相关关系和子系统失效率比值,获得每个子系统的失效率;
步骤9-b所述的根据步骤1所设定的目标可靠度、子系统的相关系数、子系统的故障相关关系和子系统失效率比值,获得每个子系统的失效率,具体公式如下:
其中,Rs *(t)表示设定的目标可靠度,m表示具有相关性的子系统个数,s表示相互独立的子系统个数,n表示子系统的个数,Ri(t)表示相关第i子系统的可靠度,λi表示第i个子系统失效率比值,ai表示分配向量的第i个元素,a1表示分配向量的第1个元素,t表示时间,表示不相关第子系统的可靠度,表示具有故障相关性的子系统序数;表示具有故障独立的子系统序数;θ表示相关系数;
步骤9-c、将每个子系统的失效率转换为可靠度,即获得最终每个子系统的可靠度;
步骤10、根据所获得的每个子系统的可靠度,设计目标机床。
2.根据权利要求1所述的基于FMEA分析的数控机床可靠性综合分配方法,其特征在于,步骤4所述的根据上述失效严重度和失效频率,进一步获得所有子系统的相对失效严重度平均值和相对失效频率平均值,具体为:在每个子系统不同失效模式的失效严重度中选择最大值作为相对失效严重度,再获得所有子系统的相对失效严重度平均值;在每个子系统不同失效模式的失效频率中选择最大值作为相对失效频率,再获得所有子系统的相对失效频率平均值。
3.根据权利要求1所述的基于FMEA分析的数控机床可靠性综合分配方法,其特征在于,步骤5所述的低区取值范围为1~3级,中区取值范围为4~7级,高区取值范围为8~10级。
4.根据权利要求1所述的基于FMEA分析的数控机床可靠性综合分配方法,其特征在于,步骤7所述的根据不同分区所设定的子系统的可靠性分区度和三阶转换函数,将每个子系统的严重度和失效率有线性关系转换为非线性关系,具体如下:
三阶转换函数公式如下:
其中,E表示FMEA严重程度或失效程度,ξ表示相对失效严重度平均值或相对失效频率平均值,a0、c0和c1均为待定系数;
三阶转换函数与可靠性分区度关系式如下:
其中,Λs表示可靠性分区度;
将公式(2)代入至公式(1)中求解a0、c0和c1系数的取值,即获得三阶转换函数即完成每个子系统的严重度和失效率从有线性关系向非线性关系的转换。
5.根据权利要求1所述的基于FMEA分析的数控机床可靠性综合分配方法,其特征在于,步骤8所述的获得子系统分配向量,具体步骤如下:
步骤8-1、将转换后获得的每个子系统失效严重度和失效频率作为两个可靠性评价指标,并根据其他可靠性评价指标,建立每个可靠性评价指标情况下,子系统两两之间的失效率分配相对值矩阵;
步骤8-2、根据失效率分配相对值矩阵,进一步获得每个子系统在不同可靠性评价指标下的相对分配系数矩阵;
步骤8-3、采用专家打分的方法对每个可靠性评价指标的权重进行打分,获得分配指标权重向量;
步骤8-4、将分配指标权重向量与相对分配系数矩阵相乘,获得子系统分配向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510358432.4A CN104992011B (zh) | 2015-06-24 | 2015-06-24 | 一种基于fmea分析的数控机床可靠性综合分配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510358432.4A CN104992011B (zh) | 2015-06-24 | 2015-06-24 | 一种基于fmea分析的数控机床可靠性综合分配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104992011A CN104992011A (zh) | 2015-10-21 |
CN104992011B true CN104992011B (zh) | 2018-08-28 |
Family
ID=54303825
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510358432.4A Active CN104992011B (zh) | 2015-06-24 | 2015-06-24 | 一种基于fmea分析的数控机床可靠性综合分配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104992011B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105700474A (zh) * | 2016-01-18 | 2016-06-22 | 长春工业大学 | 一种Copula函数可靠性动态影响度分析方法 |
CN106202812A (zh) * | 2016-07-26 | 2016-12-07 | 东北大学 | 一种基于子系统相关性模型的数控车床可靠性分配方法 |
CN107038321B (zh) * | 2017-05-24 | 2020-05-22 | 重庆大学 | 基于元动作单元的任务可靠性预计分析方法 |
CN107479494B (zh) * | 2017-07-28 | 2019-05-28 | 东北大学 | 一种基于数控机床子系统相关性的可靠性分配方法 |
US10814448B2 (en) | 2018-08-28 | 2020-10-27 | Dalian University Of Technology | Comprehensive performance evaluation method for CNC machine tools based on improved pull-off grade method |
CN110309550B (zh) * | 2019-06-10 | 2021-06-22 | 北京交通大学 | 一种基于势能场与网络效率的高速列车系统可靠性分析方法 |
CN110322148B (zh) * | 2019-07-05 | 2022-04-05 | 江苏徐工工程机械研究院有限公司 | 一种多层级任务可靠性分配方法及系统 |
CN110619179B (zh) * | 2019-09-20 | 2022-02-11 | 北京航空航天大学 | 一种机械密封可靠性分析方法及系统 |
CN112528551B (zh) * | 2020-10-21 | 2023-09-01 | 广东石油化工学院 | 一种双伺服刀架的频率可靠性评估方法 |
CN112668861B (zh) * | 2020-12-23 | 2022-10-25 | 吉林大学 | 一种基于二元语义和pfowa算子的数控机床可靠性分配方法 |
CN114564853B (zh) * | 2022-04-29 | 2022-07-26 | 希维科技(广州)有限公司 | 一种基于fmea数据的评价报告的生成方法及电子设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064340A (zh) * | 2011-10-21 | 2013-04-24 | 沈阳高精数控技术有限公司 | 一种面向数控机床的故障预测方法 |
CN104595283A (zh) * | 2015-01-24 | 2015-05-06 | 吉林大学 | 数控机床液压系统可靠性试验装置及试验方法 |
-
2015
- 2015-06-24 CN CN201510358432.4A patent/CN104992011B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103064340A (zh) * | 2011-10-21 | 2013-04-24 | 沈阳高精数控技术有限公司 | 一种面向数控机床的故障预测方法 |
CN104595283A (zh) * | 2015-01-24 | 2015-05-06 | 吉林大学 | 数控机床液压系统可靠性试验装置及试验方法 |
Non-Patent Citations (2)
Title |
---|
"Comprehensive Reliability Allocation Method for CNC Lathes Based on Cubic Transformed Functions of Failure Mode and Effects Analysis";YANG Zhou等;《CHINESE JOURNAL OF MECHANICAL ENGINEERING》;20150331;第28卷(第2期);第315-324页 * |
基于Copula理论的数控装备故障相关性;张英芝;《吉林大学学报工学版》;20111130;第41卷(第6期);第1636-1640页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104992011A (zh) | 2015-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104992011B (zh) | 一种基于fmea分析的数控机床可靠性综合分配方法 | |
CN109685340B (zh) | 一种配电设备健康状态评估方法及系统 | |
CN110303380B (zh) | 一种数控机床刀具剩余寿命预测方法 | |
Chang et al. | Applying the concept of exponential approach to enhance the assessment capability of FMEA | |
US20140067327A1 (en) | Similarity curve-based equipment fault early detection and operation optimization methodology and system | |
CN104318482A (zh) | 一套智能配电网综合评估体系和方法 | |
Dehning et al. | Factors influencing the energy intensity of automotive manufacturing plants | |
CN102848266B (zh) | 一种机床主轴精度预测方法 | |
CN109615121B (zh) | 基于数据驱动支持向量机的高速列车车轴温度预测方法 | |
Jiang et al. | Applying grey forecasting to predicting the operating energy performance of air cooled water chillers | |
CN102736546A (zh) | 一种流程工业复杂机电系统的状态监测装置及方法 | |
CN102880808A (zh) | 一种雷达结构维修性综合评价方法 | |
CN106202812A (zh) | 一种基于子系统相关性模型的数控车床可靠性分配方法 | |
CN109085813A (zh) | 故障的诊断方法和装置 | |
Amoah et al. | FDI, inflation, exchange rate and growth in Ghana: Evidence from causality and cointegrated analysis | |
CN104991515A (zh) | 数控机床全生命周期重要性测度分析方法 | |
CN104062904B (zh) | 基于cca-pls的大化工过程分布式建模方法 | |
CN110414086B (zh) | 一种基于灵敏度的综合应力加速因子计算方法 | |
CN106204324A (zh) | 一种确定电厂复杂设备关键监测参数及各参数权重分配的方法 | |
Liu et al. | A novel GREY‒ASMAA model for reliability growth evaluation in the large civil aircraft test flight phase | |
CN112116213A (zh) | 基于fahp的数控机床直线进给系统可靠性分配方法 | |
Jiang et al. | Construction and simulation of failure distribution model for cycloidal gears grinding machine | |
Hu et al. | A new and practical reliability allocation method for a complex system of NC turrets | |
CN109190968B (zh) | 空管运行效率等级评估方法及其装置 | |
CN108985564B (zh) | 一种基于qfd及二元语义的fmea方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |