1.一种机床运动机构可靠性灵敏度一阶指标分析方法,其特征在于:包括如下步骤:S1:建立机床运动机构可靠性灵敏度一阶指标计算模型RSi:
其中,RSi表示机床可靠性灵敏度一阶指标,R表示机床运动机构的可靠度,Xi表示随机变量,IS表示机床运动机构安全域S的指示函数,E(.)表示数学期望;
S2:将目标随机变量代入所述机床运动机构可靠性灵敏度一阶指标计算模型,获得所述目标随机变量的可靠性灵敏度一阶指标值。
2.根据权利要求1所述的机床运动机构可靠性灵敏度一阶指标分析方法,其特征在于:所述机床运动机构的可靠度R采用如下方法计算,
R=P{Is(X)=1} (2)
其中,X=(x1、x2…xn),X表示随机变量,IS(X)表示机床运动机构的功能函数的安全域S的指示函数,所述IS(X)采用如下方法计算,当随机变量属于所述机床运动机构的功能函数的安全域时,IS(X)的取值为1,当随机变量不属于所述机床运动机构的功能函数的安全域时,IS(X)的取值为0,如下式所示:其中,S表示机床运动机构的功能函数的安全域,所述S采用如下方法计算,S={X:|g(x)|≤ε} (4)其中,X=(x1、x2…xn),X表示随机变量,ε表示误差极限,g(x)表示机床运动机构的功能函数。
3.根据权利要求2所述的机床运动机构可靠性灵敏度一阶指标分析方法,其特征在于:所述机床运动机构的功能函数g(x)采用如下方法计算:
其中,g(x)表示机床运动机构的功能函数,Y(x)表示机床运动机构的运动输出,Y0表示机床的理想运动输出,βc表示输出层的阈值,Vbc表示隐藏层节点与输出层节点间的网络权值,αb表示隐藏层节点的阈值,Wab表示输入层节点与隐藏层节点间的网络权值,f1(.)表示S型非线性函数tansig,f2(.)表示线性函数purelin。
4.一种机床运动机构可靠性灵敏度二阶指标分析方法,其特征在于:S1:建立机床运动机构可靠性灵敏度二阶指标计算模型RSSij:其中,RSSij表示机床可靠性灵敏度二阶指标,Is表示机床运动机构安全域S的指示函数,Vi和Vj均表示单独一个随机变量的方差,R表示运动机构的可靠度,E(.)表示数学期望,其中,Vi和Vj的区别仅在于随机变量的不一样,二者均采用同样的计算模型来进行计算,Vi=V(Is)-E[V(Is|Xi)] (7)
其中,Xi表示随机变量,IS表示机床运动机构安全域S的指示函数,E(.)表示数学期望,V(.)表示方差,V(IS)表示机床运动机构安全域S的指示函数IS的方差;
其中,V(Is)=R(1-R) (8)
其中,R表示机床运动机构的可靠度;
S2:将目标随机变量代入所述机床运动机构可靠性灵敏度二阶指标计算模型,获得所述目标随机变量的可靠性灵敏度二阶指标值。
5.根据权利要求4所述的机床运动机构可靠性灵敏度二阶指标分析方法,其特征在于:机床运动机构的可靠度R采用如下方法计算:
R=P{Is(X)=1} (9)
其中,X=(x1、x2…xn),X表示随机变量,IS(X)表示机床运动机构的功能函数的安全域S的指示函数,所述IS(X)采用如下方法计算,当随机变量属于所述机床运动机构的功能函数的安全域,则IS(X)的取值为1,当随机变量不属于所述机床运动机构的功能函数的安全域,则IS(X)的取值为0,如式(4)所示:其中,S表示机床运动机构的功能函数的安全域,所述S采用如下方法计算,S={X:|g(x)|≤ε} (11)其中,X=(x1、x2…xn),X表示随机变量,ε表示误差极限,g(x)表示机床运动机构的功能函数;
其中,所述机床运动机构的功能函数g(x)采用如下方法计算:
其中,g(x)表示机床运动机构的功能函数,Y(x)表示机床运动机构的运动输出,Y0表示机床的理想运动输出,βc表示输出层的阈值,Vbc表示隐藏层节点与输出层节点间的网络权值,αb表示隐藏层节点的阈值,Wab表示输入层节点与隐藏层节点间的网络权值,f1(.)表示S型非线性函数tansig,f2(.)表示线性函数purelin。
6.一种机床运动机构可靠性灵敏度总指标分析方法,其特征在于:
S1:建立机床运动机构可靠性灵敏度总指标计算模型RSTi:
其中,RSTi表示机床可靠性灵敏度总指标,R表示机床运动机构的可靠度,E(.)表示数学期望,IS表示机床运动机构安全域S的指示函数,X-i表示随机变量X中除去Xi以外的n-1维变量;
S2:将目标随机变量代入所述机床运动机构可靠性灵敏度总指标计算模型,获得所述目标随机变量的可靠性灵敏度总指标。
7.根据权利要求6所述的机床运动机构可靠性灵敏度总指标分析方法,其特征在于:机床运动机构的可靠度R采用如下方法计算:R=P{Is(X)=1} (14)
其中,X=(x1、x2…xn),X表示随机变量,IS(X)表示机床运动机构的功能函数的安全域S的指示函数,所述IS(X)采用如下方法计算,当随机变量属于所述机床运动机构的功能函数的安全域,则IS(X)的取值为1,当随机变量不属于所述机床运动机构的功能函数的安全域,则IS(X)的取值为0,如式(4)所示:其中,S表示机床运动机构的功能函数的安全域,所述S采用如下方法计算,S={X:|g(x)|≤ε} (16)其中,X=(x1、x2…xn),X表示随机变量,ε表示误差极限,g(x)表示机床运动机构的功能函数。
8.根据权利要求7所述的机床运动机构可靠性灵敏度总指标分析方法,其特征在于:所述机床运动机构的功能函数g(x)采用如下方法计算:其中,g(x)表示机床运动机构的功能函数,Y(x)表示机床运动机构的运动输出,Y0表示机床的理想运动输出,βc表示输出层的阈值,Vbc表示隐藏层节点与输出层节点间的网络权值,αb表示隐藏层节点的阈值,Wab表示输入层节点与隐藏层节点间的网络权值,f1(.)表示S型非线性函数tansig,f2(.)表示线性函数purelin。
9.一种机床运动机构的可靠性灵敏度高阶指标分析方法,其特征在于:S1:建立机床运动机构可靠性灵敏度高阶指标计算模型RSHi:
其中,RSHi表示机床可靠性灵敏度高阶指标,Xi表示随机变量,X-i表示变量X中去除Xi以外的n-1维向量,R表示机床运动机构的可靠度,E(.)表示数学期望,V(.)表示方差,Vij所述任意两个变量之间的交互方差;
其中,机床运动机构的可靠度R采用如下方法计算:
R=P{Is(X)=1} (19)
其中,X=(x1、x2…xn),X表示随机变量,IS(X)表示机床运动机构的功能函数的安全域S的指示函数,所述IS(X)采用如下方法计算,当随机变量属于所述机床运动机构的功能函数的安全域,则IS(X)的取值为1,当随机变量不属于所述机床运动机构的功能函数的安全域,则IS(X)的取值为0,如式(4)所示:其中,S表示机床运动机构的功能函数的安全域,所述S采用如下方法计算,S={X:|g(x)|≤ε} (21)其中,X=(x1、x2…xn),X表示随机变量,ε表示误差极限,g(x)表示机床运动机构的功能函数;
所述机床运动机构的功能函数g(x)采用如下方法计算:
其中,g(x)表示机床运动机构的功能函数,Y(x)表示机床运动机构的运动输出,Y0表示机床的理想运动输出,βc表示输出层的阈值,Vbc表示隐藏层节点与输出层节点间的网络权值,αb表示隐藏层节点的阈值,Wab表示输入层节点与隐藏层节点间的网络权值,f1(.)表示S型非线性函数tansig,f2(.)表示线性函数purelin;
S2:将目标随机变量代入所述机床运动机构可靠性灵敏度高阶指标计算模型,获得所述目标随机变量的可靠性灵敏度高阶指标值。
10.根据权利要求9所述的机床运动机构的可靠性灵敏度高阶指标分析方法,其特征在于:所述任意两个变量之间的交互方差Vij采用如下方法计算:Vij=E[E2(IS|Xi,Xj)]-R2-Vi-Vj (23)其中,Xi和Xj表示随机变量,R表示机床运动机构的可靠度,E(.)表示数学期望,Vi和Vj的表示单个变量的方差,Vi和Vj的区别仅在于随机变量的不一样,二者均采用同样的计算模型来进行计算,Vi=V(Is)-E[V(Is|Xi)] (24)
其中,Xi表示随机变量,IS表示机床运动机构安全域S的指示函数,E(.)表示数学期望,V(.)表示方差,V(IS)表示机床运动机构安全域S的指示函数IS的方差。