一种基于最大熵有序加权平均的数控机床可靠性分配方法与流程

本发明涉及高档数控机床与基础制造装备的设计领域，具体是涉及数控机床的可靠性分配方法

背景技术

可靠性分配是确定产品可靠性和竞争力时需要考虑的最重要的因素之一，也是可用于提高系统可靠性的重要分析工具。近年来，由于用于高科技工业过程的工程系统的复杂性增加，对系统可靠性的关注有所增加

可靠性分配是为了把系统的可靠性指标按照一定的准则分配给系统各组成单元而进行的工作。其目的就是将整个系统的可靠性指标转换为每一个分系统或单元的可靠性指标，使之协调一致。它是一个由整体到局部，自上到下的分解过程。

目前中常用的可靠性分配方法主要有等分配法、考虑重要度和复杂度的分配法、比例组合法和评分分配法，这些传统的可靠性分配方法没有考虑到可靠性分配的影响因素，或者考虑得不够全面。在针对数控机床这种集机、电、液、光、气、光为一体，且任务异常复杂的系统进行可靠性分配时，上述传统方法很明显不能达到我们的目的。

技术实现要素：

针对现有机床可靠性分配方法的不足，本发明提出一种基于最大熵有序加权平均的数控机床可靠性分配方法，有效地解决了传统目标可行性技术的因子i、s、p、e的测量尺度的等权重问题。

本发明的特征在于通过在目标可行性技术的基础上引入了条件参数α(提供了有关子系统的有价值的信息，可以用于更好地协助设计者做出正确的可靠性分配决策)，并结合通过平均算子导出多项式方程确定的最大熵下的最优加权向量，从而确定最大熵下数控机床各子系统的最优ispe值，从而能够在数控机床子系统合理分配的可靠性水平内准确有效地分配可靠性等级，控制合理支持成本，降低机床维护成本。

一种基于最大熵有序加权平均的数控机床可靠性分配方法，主要包括以下步骤：

s1列出数控机床系统和子系统的结构；

s2定义数控机床可靠性和任务执行时间；

s3为每个子系统分别确定复杂性i、技术发展水平s、工作时间p、环境e的标度，子系统分配因子作为i、s、p和e的数值等级的函数来计算。

公式(1)：λs＝-ln(r)/t

公式(2)：ω'k＝r'ik×ω1+r'sk×ω2+r'pk×ω3+r'ek×ω4

公式(3)：

公式(4)：

公式(5)：

s4计算该数控机床每一个子系统的分配等级r’,然后根据对机床可靠性的要求，确定该目标机床条件参数α，α＝0.5,0.6,0.7,0.8,0.9,1.0，α的值从低到高表示设计人员对目标机床可靠性的乐观程度，α＝1用于表示决策者最乐观的情况，当决策者面临适度评估时，α＝0.5。当可用信息在可靠性分配设计阶段不准确，不完整或不确定时，条件参数尤其有用，根据确定的条件参数然后确定最大熵下的权重向量，ω1、ω2、ω3、ω4，不同条件参数下的权重向量数值可参见表1，然后利用公式(2)求得第k个子系统的总体等级ω’k

s5利用公式(3)和公式(5)求得第k个子系统的复杂度c’k

s6根据数控机床系统规范来计算系统故障率，假设r是系统可靠性，t是任务持续时间，系统故障率λs计算如下：λs＝-ln(r)/t

s7利用公式(4)计算分配的子系统故障率

s8通过第七步求得的各子系统故障率，互相比较，可以确定目标数控机床系统的最优可靠性分配等级。

其中ω'k指的是第k个子系统的评级，r'ik、r'sk、r'pk、r'ek指的是第k个子系统的ispe值，c'k指的是第k个子系统的复杂度，ω'指的是整个系统的评级，为分配给第k个子系统的故障率,λs指的是整个系统的故障率，r为整个系统可靠性。

表1：在最大熵下的最优权重向量(n＝4)

与现有技术相比较，本发明具有如下技术效果。

在数控机床可靠性分配领域，提高可靠性的可靠性分配方法为之不多且效率低下，在目标可行性技术的基础上引入了条件参数，并结合通过平均算子导出多项式方程确定的最大熵下的最优加权向量，从而确定最大熵下数控机床各子系统的最优ispe值，从而能够在数控机床子系统合理分配的可靠性水平内准确有效地分配可靠性等级，控制合理支持成本，降低机床维护成本，大大提高了机床的可靠性。

附图说明

图1数控机床系统的结构图。

图2是本发明方法的实施流程图。

具体实施方式

本发明以某数控机床为例，对该发明能否达到最优分配来进行验证。

具体包括如下步骤:

步骤一：列出数控机床系统及子系统的结构，如图1所示

步骤二：基于设计要求和系统运行环境，数控机床系统的系统可靠性设置为0.995214，任务时间设为1000个小时，由公式(1)可计算该数控机床系统的故障率，即

λs＝-ln(r)/t＝0.0000047975＝47.975×10^-7

步骤三：为数控机床的每个子系统分别确定i、s、p、e的标度(ispe值为1～10内的一个数)，然后利用专家评判的方法分别确定数控机床各个子系统的ispe值，如表2所示

表2:数控机床子系统的ispe值

步骤四：假设目标数控机床条件参数α＝0.5，由图1可知其在最大熵下的权重向量分别为ω1＝0.25ω2＝0.25ω3＝0.25ω4＝0.25，以传动系统为例，其ispe值分别为7、5、9、2则根据公式(2)可求得其系统的评级值：ω'k＝r'ik×ω1+r'sk×ω2+r'pk×ω3+r'ek×ω4＝7×0.25+5×0.25+9×0.25+2×0.25＝5.75以此方法可依次求得剩余系统的评级值，由公式(5)可知，将上述求得各子系统的评级值累加可得该数控机床的总的评级值ω'，当α为其他值时，计算过程类似，就不再赘述了，各结果可参见表3。

表3：各子系统在条件参数不同时的评级值

步骤五：根据公式(3)以及步骤四求得的各子系统的评分以及系统总的评分，可求得各子系统的复杂度，以传动系统为例，当α＝0.5时，其他系统结果参见表4。

表4：各子系统在条件参数不同时的复杂度

步骤六：然后可以根据公式(5)来计算传动系统的故障率，其系统故障率(每10⁷小时)，其余系统结果可参见表5。

表5：各子系统在条件参数不同时的故障率(每10⁷小时)

结论：

(1)就α＝0.5而言，由表4可知，进给系统以及数控系统的值最高，自然，其可靠性分配等级在α＝0.5时应最高

(2)在考虑其他α值各子系统的值并比较后，可知，数控系统的值综合来说是最高的，所以数控系统应当在可靠性分配时最先考虑，其综合可靠性分配等级应该是最高的，其次是进给系统。