含非正态随机变量的机械零件单模可靠度统一模型

含非正态随机变量的机械零件单模可靠度统一模型

0机械零件单模失效可靠度求解理论框架机械产品的可靠性设计不同于电子可靠性设计。在机械系统中，有三个不同形式的相关问题。首先，考虑单件损失机械的随机变量相关性，其次，考虑零件失败模型之间的相关性，第三，机械系统零件之间的相关性。在分析机械系统单件的可靠性时，应充分考虑上述两种相关问题。在当前的研究工作中，重点关注设备零件（单台故障）之间的相关问题。机械零件可靠度计算的基本理论为应力-强度干涉理论(又称为应力-强度相干模型),在以往的文献中,此理论都是在应力和强度完全独立的情况下导出的,但实际上并不是所有的零件可靠性问题中的应力和强度完全独立,例如轴系零件设计变量的荷载、强度等因素之间存在程度不同的相关性引起机械零件广义应力与强度随机性的因素主要有材料、加工、制造、安装、荷载等,如同种材料性能分布的不均匀、制造精度误差.正是由于这些参数的随机性才有了可靠度的概念.要想得到这些因素的实际随机性参数是比较困难的,因此在理论分析中常常假设这些因素服从常见的分布,如正态分布、威布尔分布、指数分布等.现有的大多数文献中假设随机参数服从正态分布,但很多情况下并非如此.这时就需要探求随机变量因素服从其他分布时的可靠度计算公式,探求将非正态分布转化为正态分布的方法,利用正态分布运算的成熟的计算技术来求得零件可靠度.本文针对引起机械零件单模失效的随机变量为不相关正态随机变量、相关正态随机变量、不相关非正态随机变量、相关非正态随机变量4种情况,借鉴结构可靠性求解理论,寻求建立机械零件单模失效可靠度求解模型的统一方法.主要是建立相关非正态随机变量的可靠度求解理论模型.不相关正态随机变量的可靠度理论是求解其他3种情况的基础.最后,通过齿轮传动的实例说明本文所提方法的有效性.1单用户模型的失败可靠性由统一模型构建1.1u3000随机变量y标准面为5.此为求解零件可靠度最简单也是基础的模型,如图1中Ⅰ框所示,计算步骤如下:步骤1正态随机变量正则化为标准正态随机变量.假设影响极限状态方程的随机变量为Y=(Y其中Y步骤2零件单模可靠度的求解.将极限状态方程(1)以正则化后的标准正态随机变量表示为式(3),则式(3)分为两种情况:随机变量Z(1)若极限状态方程(3)为n维标准正态随机变量Z此时可靠度指数可靠度(2)若极限状态方程(3)为非线性组合,则可采用两种方法求解.1一次二阶矩(first-orderu3000second-moment,FOSM)法将式(3)在平均值Z处一阶泰勒展开,仅保留一阶项:则零件单模可靠度可由式(5)、(6)求得.2改进的一次二阶矩(advancedu3000first-ordersecond-moment,AFOSM)法通过求解式(9)得到可靠度指数β,代入式(6)求得可靠度.对于式(9)的求解可以采用Lagrange乘子迭代法导出的公式进行迭代计算,将在2.1中详细介绍.1.2u3000随机变量正态分布相关时零件单模可靠度求解理论假设影响极限状态方程的随机变量为X*=(X*其中步骤1正态随机变量正则化为标准正态随机变量.按照式(2)、(3)的形式将相关的非标准正态随机变量X*=(X*步骤2相关的标准正态随机变量转化为不相关的正态随机变量.将相关的标准正态随机变量Y*=(Y*变换过程如下:设随机变量Y*=(Y*由线性代数理论可知,矩阵C其中A是一个正交矩阵,它的列向量等于矩阵C对角矩阵C步骤3变量正态分布相关时零件单模可靠度求解.经过步骤1、步骤2所述方法变换后,可靠度的求解转化为1.1中讨论的问题,按照1.1中提供的方法就可以得到考虑随机变量正态分布且相关时的零件单模失效下的可靠度.1.3u3000随机变量中含有非正态随机变量且相关时零件单模可靠度求解理论研究发现,自然界中大多数的物理量在经过无数次的统计分析后都近似服从正态分布,因此有关正态分布函数的性质、特征、数字特征等研究也最多,发展了许多可用的结论及大量的可供参考查阅的数据.基于此现行对于由高维极限状态函数表示的可靠度关系式的求解都是以各随机变量服从正态分布为前提下得出的.然而面对多维极限状态函数中含有非正态随机变量的情形,只有将非正态随机变量转化为正态随机变量,才能利用现有成熟的、计算精度高的可靠度理论进行机械零件可靠度的求解,本节将对随机变量中含有非正态随机变量的情形进行研究.假设机械零件单模失效的极限状态方程为其中X=(X步骤1将X中的非正态随机变量转化为正态随机变量.(1)Rosenblatt变换(Rosenblattu3000transformation,RT)Rosenblatt证明:任意非正态随机变量X,理论上都可以通过Rosenblatt变换转化为线性无关的标准正态随机变量Z.但这种变换具有局限性,需要知道n维随机变量X的联合累积分布函数,具体变换形式如下:设n维随机变量X的联合累积分布函数为F按顺序反求式(15)中的方程,可以得到n维标准正态线性无关随机变量Z:式(16)中的条件概率密度函数通过联合概率密度函数按式(17)、(18)求得:同理可通过逆变换按顺序反求式(16)得到X=(X(2)R-F变换设X由式(19)得变换得到正态随机变量X*=(X*步骤2含有非正态随机变量且相关时零件单模可靠度求解.利用Rosenblatt变换方法进行非正态随机变量的正态化时,必须利用随机变量的联合概率密度函数,这一点在实际中较难实现;但经过此种变换后,含有非正态随机变量且相关的随机变量经转换后直接转化为不相关的标准正态随机变量,此时利用1.1中所述的方法就可以快速地求出其可靠度.R-F变换方法简单,较容易在计算机上实现,并不需要联合概率密度函数.将非正态随机变量正则化为标准正态随机变量Y*=(Y2机械零件单模极限状态方程本章将根据第1章介绍的机械零件单模失效可靠度求解统一模型理论探讨可靠度求解的计算机实现算法,机械零件单模极限状态方程为式(14)描述的n维方程,含有非正态随机变量,且各随机变量之间存在相关性.算法中非正态随机变量利用R-F变换方法正则化为正态随机变量;可靠度指数β的求解采用H-L方法,利用Lagrange乘子迭代法导出的公式进行迭代计算.2.1lagrange乘子迭代法Lagrange乘子迭代法详细叙述如下:设求L的极小值得令则式(23)可表示为因此所以将式(28)代入式(26)得将式(29)代入式(21)得采用数值方法可以计算出式中:β则式(24)、(29)、(30)即可组成Lagrange乘子迭代法计算的迭代公式.2.2机械零件单模可靠度计算假设影响零件单模可靠度的非正态随机变量为X=(X步骤2在设计点x=(x步骤4对相关的标准正态随机变量y*=(y步骤7将z步骤8采用步骤7计算得到的β重新计算z重复步骤6~8直到计算结果收敛(以z稳定作为收敛指标).步骤9计算步骤10重复步骤2~9直到计算结果收敛(以步骤11计算步骤12由式(5)求得可靠度指数,式(6)求得机械零件单模失效可靠度.3u3000可靠度分析本文以某1.6MW风力发电机平行级直齿圆柱齿轮传动为研究对象,通过对风机齿轮箱传动系统受力分析得到直齿圆柱齿轮传动各项参数:齿轮材料及热处理为17CrNiMo6低碳合金钢、渗碳淬火,小齿轮转矩T按国家标准规定的方法计算或查线图得到与齿面接触疲劳可靠性相关各参量的均值以上各参数的含义及应力σ(1)假设各随机变量均服从正态分布,其均值和方差如前所述,各随机变量之间相互独立,此种情况下根据文献[21]提供的方法计算得到可靠度指标β=1.374u3000112,R=0.915u3000330u30005.(2)假设各随机变量均服从正态分布,相关随机变量之间的相关性以齿面接触疲劳应力-强度干涉理论中应力σ当相关系数ρ在区间[-1,1]内变化,得到可靠度指数β、可靠度R随ρ变化曲线如图2所示.由图2可知,在相关系数ρ从-1增大到1(即应力、强度相互独立到应力、强度完全相关)时,可靠度指数β从1.063u30001增大到2.400u30008,相应的可靠度R从0.856u3000010增大到0.991u3000802,可靠度R的相对变化率为15.86%,变化较大.此时若忽略应力、强度之间的相关性将会产生两种结果:若实际为正相关,此时利用独立假设理论就会低估系统本身可靠性,通过改变设计来提高其可靠度时设计成本增高,造成不必要的浪费;若实际为负相关,则高估了系统的可靠性,按照此可靠性设计的产品,因过早的失效破坏而达不到设计要求.以简单的应力、强度双变量相关的模型为基础,分析相关系数对不同零件单模可靠度的影响.以应力、强度相互独立时可靠度为基准,求解可靠度随相关系数在[0,1]内的变化曲线,通过对多种具有不同基准可靠度的机械零件ρ-R曲线分析知,基准可靠度越小,考虑与不考虑相关系数时的可靠度差别越大.因此可以得出结论,当在假设变量相互独立情况下求得的可靠度较大时(如0.999),变量之间的相关性可以忽略,以它们相互独立求得可靠度,这样不仅简化了计算,还不会影响机械零件可靠性优化设计的精度.(3)假设各随机变量服从正态分布,通过分析国家标准GBu30003480—1997中各随机参数的计算知某些随机变量之间具有相关性:F根据各随机变量计算式设定相关随机变量相关系数为假设相关随机变量线性变换为非相关变量后的随机变量为z=(z(4)假设随机变量F由以上讨论可知,本文介绍的变量相关的机械零件单模可靠度求解算法是可行的,且具有较高的收敛速度,对于复杂的变量相关的可靠度的求解具有重要的应用价值.另外,从考虑各个变量之间的相关性与将只考虑广义应力、强度之间单相关性求得的可靠度可以看出,差别很大,因此在求解变量相关可靠度时,为使求得的可靠度比较精确,应尽量将各个相关变量的相关性做单独处理,不要将其等效为应力、强度之间的相关性.4u3000结u3000论(1)本文以系统中机械零件的单模失效为研究对象,通过探讨影响单模失效极限状态方程的随机变量的不同特征:是否含有非正态随机变量、随机变量间是否相关,建立了考虑含有非正态随机变量和随机变量具有相关性时的机械零件单模失效下可靠度求解的统一模型,并对统一模型求解的理论进行了详细的论述.(2)基于所建立的机械零件单模失效下的可靠度求解统一模型,提出了高效的计算机求解算法.(3)通过算例说明了随机变量之间的相关性对零件单模失效可靠度的影响,表明在零件单模失效求解中忽略随机变量的相关性将会带来较大的误差.通过分析得到,当基准可靠度较大时(如0.999u30009),变量之间的相关性可以忽略,这样在不影响优化设计精度的前提下简化了可靠度的计算.(4)本文提出的理论为含有非正态随机变量时零件单模失效下可靠度的求解提供了一个很好的方法.实例证明,所提出的计算机实现算法收敛速度快,具有很强的实用性.本文所建