基于鞍点估计的机械零部件可靠性灵敏度分析摘要对机械结构来说，可靠性指标一般随材料特性、几何参数、工作环境等不确定性因素变化而减弱，所以结构的可靠度、灵敏度就显得尤为重要，对机械零部件可靠性灵敏度的分析也是亟不可待。

本文利用鞍点估计技术可以无限逼近非正态变量空间中线性极限状态函数概率分布的特点，能有效解决统计资料或实验数据较少而难以确定设计变量的分布规律的问题。

将可靠性设计理论、灵敏度分析技术与鞍点逼近理论相结合，以前面可靠性数学模型为基础，系统地推导了基于鞍点估计的可靠性灵敏度公式，讨论了基于鞍点估计法的机械零部件可靠性灵敏度计算问题，为进一步分析机械零部件的可靠性稳健设计奠定了理论基础。

关键词：不确定性鞍点灵敏度可靠性第一章绪论1.1机械可靠性设计理论研究进展很早以来人们就广泛采用“可靠性”这一概念来定性评价产品的质量问题，这只是靠人们的经验评定产品可靠还是不可靠，并没有一个量的标准来衡量；从基于概率论的随机可靠性到基于模糊理论的模糊可靠性再到非概率可靠性以及最近提出的结构系统概率-模糊-非概率混合可靠性，表明定量衡量产品质量问题的理论方法从产生到现在已有了长足的发展；对于复杂结构的复杂参数由单纯的概率非概率可靠性分析方法发展到可靠性灵敏度分析的各种分析方法，使得这一理论日续丰富和完善，并深入渗透到各个学科和领域。

可靠性当今已成为产品效能的决定因素之一，作为一个与国民经济和国防科技密切相关的科学，未来的科技发展中也必将得到广泛的研究和应用。

20世纪初期把概率论及数理统计学应用于结构安全度分析，已标志着结构可靠性理论研究的初步开始。

20世纪40年代以来，机械可靠性设计理论有了长足的发展，目前为止己经相当成熟，尤其是许多国家幵始研究在结构设计规范中的应用，使机械可靠性设计理论的应用进入一个新的时期。

1.2机械可靠性设计理论研究现状在实际工程中，不确定因素的存在在所难免，可靠性分析与这些不确定性紧密相关。

在过去的几十年中，概率论在各种工业系统的可靠性评估方面获得了巨大的成功，概率可靠性方法成为处理不确定性的最为普遍的方法。

但随着科学技术的发展，人们逐渐认识到，除了随机性以外，在工程中还存在着另一类重要信息——模糊信息。

所以传统的可靠性方法就是用概率论和模糊理论处理不确定性，但概率可靠性和模糊可靠性模型都需要用较多的数据去定义参数的概率分布或隶属函数，a计算量较大。

近年来的有关研究[1]表明，概率可靠性对概率模型参数很敏感,概率数据的小误差可导致结构可靠性计算出现较大误差[2]，说明在没有足够的数据信息描述概率模型时，在主观的假设下概率可靠性计算的结果是不可靠的[3]。

模糊性和随机性是不同的两类不确定性，其产牛机理和物理意义均有一定差异，在机械和结构系统的分析和设计中，由于各种因素的影响，常不可避免地同时存在随机的和模糊的不确定性[4]。

此时，必须同时考虑随机性和模糊性。

对此问题，常用的方法是依据Zadeh提出的模糊概率计算公式[5]，综合考虑功能状态变量的随机分布和模糊隶属情况，给出一确定的失效概率或可靠度值。

Ben-Haim首先提出了不用概率定义即非概率的可靠性概念。

对于不确定信息的描述，不采用随机变量，不用极限状态函数和概率密度函数，而采用凸集合模型描述，经过分析可得到输出（或响应）的变化范围，将此变化范围与要求的变化范围比较即可得到安全程度的度量指标。

非概率可靠性提出后引起了很多人关注。

郭书样等[6]基于区间分析，提出一种非概率可靠性度量指标，来衡量不确定参数为区间变量时系统的安全性，并将其用于结构优化，但对于其它凸模型情况，区间算术计算结果趋于保守。

邱志平[7]等指出了Ben-Haim鲁棒可靠性准则即响应凸集合与失效凸集合为不相交关系的错误，提出了结构的安全与失效的关系应该对应于凸集合间的偏序关系，这仅仅是基于一种方式对不确定量进行的描述，而计算必须建立在先前的实验数据上，以至该模型也有一定的局限性.曹鸿钧[8]等在区间可靠性指标的基础上，提出一衡量超椭球模型与区间模型并存情况下的非概率可靠性指标由于非概率可靠性指标是用一个极小极大模型定义的，虽然可以采用增广设计变量的方法转化为常规极值问题，但以该指标为设计约束的优化问题求解时仍然十分困难。

在结构的分析和设计中，需要合理地定量处理一些影响其性能的不确定性，虽然概率理论在不确定性的处理及可靠性分析方面得到了成功的应用，但不确定性并不等于随机性。

不确定性的模拟既可以是概率的，也可以是非概率的，同一问题中可能同时含有概率变量和非概率变量。

因此，非概率方法及其与概率方法的混合模型[9]的研究也有着重要的理论和实际意义。

非概率可靠性并不是完全否认概率可靠性，它是概率可靠性有益的补充，很多研究都在试图将两种方法结合起来对系统性能进行评估。

基于结构可靠性分析中的概率可靠性模型和非概率集合可靠性模型，王军等[10]等提出一种新的结构可靠性分析的概率-非概率混合模型，该模型首先将功能函数进行非概率可靠性分析，后将标准化区间变量空间所有区域的可靠度惊喜求和计算，从而给出结果的可靠度。

为了有效地处理结构系统的混合可靠性问题，基于模糊随机可靠性模型及非概率集合可靠模型，尼早等[11]建立了结构概率—模糊-非概率混合可靠性模型，该混合可靠性模型为分析和设计决策提供更全面、更真实的有用信息。

可靠性设计的精确性和先进性是建立在应力、强度、寿命等数据的真实性、精确性的基础上的，重视试验数据的收集和分析，对设计新产品时有很重大的参考价值。

1.3.机械零件可靠性设计理论的发展趋势对当前机械产品而言，如何提高设计质量、完善设计理论、改进设计技术、缩短设计周期是最重要的，而这些都与可靠性有着密切的联系。

可靠性技术己深入到机械零部件结构设计、强度设计以及失效分析中，机械零件可靠性理论研究工作已经成为机械工程中的研究热点，目前大量论文和专著，已证实了结构系统可靠性分析和计算方法相当成熟，就目前的发展趋势看如下几方面应是工程机械结构可靠性理论研究的热点：(1)可靠性优化设计(2)可靠性灵敏度设计(3)可靠性稳健设计(4)可靠性试验。

1.4本文基本思路和主要研究内容可靠性设计理论的应用已深入到机械零部件的选材、结构设计、强度设计、失效概率分析以及产品的创新设计中。

产品的设计变量和预设计参数向量对机械产品的可靠性起重要作用，而在实际工程中这些变量的不确定性是客观存在的，这些不确定性有可能导致机械零部件的性能指标有较大的波动，甚至失效，因此确定设计参数对机械零部件可靠性影响程度十分必要，即可靠性灵敏度设计目前，机械零部件可靠性灵敏度研究发展较为迅速，文献首次提出机械零部件可靠性灵敏度分析的基本概念，研究了可靠度对随机参数的敏感性。

随后国内外部分学者又提出了许多可靠性灵敏度分析方法[12][13][14],张义民吕振宙等对机械可靠性灵敏度进行了较为系统的研究。

当前可靠性灵敏度分析的方法主要有矩方法和Monte-Carlo数值模拟法等，矩方法是基于设计变量均值的可靠性灵敏度分析，计算速度较快，但只考虑了设计变量的前两阶矩，影响计算的精度。

Monte-Carlo数值模拟法，计算量较大，尤其是针对小失效概率问题，以至在实际工程问题中很难应用。

计算机械零部件的失效概率或可靠度的前提是，必须知道极限状态函数的概率密度函数或联合概率密度函数，但是实际工程通常为小样本情况，统计资料或者试验数据往往缺少，所以很难确定基本随机参数的分布规律，特别是分布函数、概率密度函数、实际分布函数尾部与概率密度函数尾部拟合不一致的情况，对结构的可靠性或者失效概率分析的精度会有较大影响[15]。

在小样本情况下，如何能在不降低计算效率的情况下改善可靠性及可靠性灵敏度分析方法的精度，是本文研究的重点和难点。

本文利用鞍点估计技术可以无限逼近非正态变量空间中线性极限状态函数概率分布的特点，获得了外载荷作用下机械零部件结构响应的概率密度函数和分布函数。

在此基础上，针对基本参数存在不确定性的可靠性分析问题，将鞍点估计理论与可靠性灵敏度分析方法相结合，系统地推导了基于鞍点估计的可靠性灵敏度公式，研究了机械零部件可靠性灵敏度分析方法。

并将其引入到新型旋转式立体车库载车台主梁结构的可靠性灵敏度分析中，通过与改进一次二阶矩显式迭代方法计算所得的结果做比较，显示出本方法计算结果的准确度高、计算速率较快优点。

第二章 可靠性灵敏度简述我们当前应用的可靠性灵敏度设计是在可靠性设计的基础上进行机械结构、机构设计参数的灵敏度分析，来确定基本设计变量的改变对机械结构、机构可靠度的影响程度，可靠性灵敏度分析可以定量化反映各基本随机参数对机械结构、机构失效的影响程度，即敏感性。

目前求得灵敏度的常用方法是采用可靠度对设计变量均值和方差的偏导数。

若可靠度对某个设计变量求得的灵敏度数值大于0，则说明该设计变量均值的增加会使机械结构/机构趋于更加可靠；若对某设计变量求得灵敏度数值小于0，则说明该设计变量均值的增加会使机械结构/机构趋于更加不可靠。

同样对于方差的灵敏度亦如此。

可靠度对设计参数均值的灵敏度表达式为X R dX dR g g T ∂∂∂∂∂∂=μμββ （2-1）可靠度对设计参数方差的灵敏度表达式为)()(d X Var R X dVar R g g ∂∂∂∂∂∂=σσββ （2-2）式中：)(βϕβ=∂∂R ，g σμβ1g =∂∂，⎥⎦⎤⎢⎣⎡∂∂∂∂∂∂=∂∂n X g X g X g 21g σβ， 2g g g σμσβ=∂∂，⎥⎦⎤⎢⎣⎡∂∂⊗∂∂=∂∂X g X g X Var g g σσ21)( 在极限状态函数为线性或非线性程度不强的情况下，上述可靠性灵敏度计算公式是实用的。

但在非线性程度较强的情况下，采用随机摄动技术或差分法计算来近似代替微分时，基本随机参数的变化量取得过小或过大都直接影响灵敏度分析的精度，甚至会出现错误。

因此对上述可靠性灵敏度计算公式进行修正。

在对设计参数向量求偏导时，通常认为极限状态函数的方差和设计参数的均值是完全独立的，所以有0g=∂∂T X σ （2-3）事实上，在极限状态函数方程为线性时，则式（2-3)成立，若极限状态函数方程为非线性即为关于X 的二阶或高阶方程时，则式（2-3)就不成立。

所以修正后的可靠性灵敏度计算公式应为Tg g T g g T X R X R dX R ∂∂∂∂∂∂+∂∂∂∂∂∂=μμββσσββd （2-4） 在这里要说明的是方差的灵敏度计算公式无需进行修正，因为极限状态函数的均值和设计参数向量的方差是完全相互独立的，而且中不包含Var(X)项，所以有0)(g=∂∂X Var μ （2-5） 所以)()()()(d X Var R X Var R X Var R X dVar R g g g g g g ∂∂∂∂∂∂=∂∂∂∂∂∂+∂∂∂∂∂∂=σσββμμββσσββ不变。