基于概率功能度量的疲勞可靠壽命分析方法

劉勤，孫志禮，錢云鵬，劉英

（1.東北大學(xué)機械工程與自動化學(xué)院，遼寧沈陽110819；2.中國兵器科學(xué)研究院，北京100089）

基于概率功能度量的疲勞可靠壽命分析方法

劉勤1，2，孫志禮1，錢云鵬2，劉英2

（1.東北大學(xué)機械工程與自動化學(xué)院，遼寧沈陽110819；2.中國兵器科學(xué)研究院，北京100089）

在分析疲勞壽命工程常用預(yù)計方法的基礎(chǔ)上，引入功能度量法，建立疲勞可靠壽命模型，提出可靠壽命計算、靈敏度分析方法。該方法將可靠壽命的計算問題轉(zhuǎn)換為一個球面約束的優(yōu)化問題，只需進(jìn)行一次概率功能度量的求解即可獲得可靠壽命；利用計算過程的中間量，導(dǎo)出了可靠壽命對隨機變量均值、標(biāo)準(zhǔn)差以及確定量的靈敏度計算公式。通過對某型車輛扭力軸的可靠壽命分析，證明基于概率功能度量的疲勞可靠壽命預(yù)計方法具有較好的效率和精度，靈敏度結(jié)果有效地反映出各變量對疲勞可靠壽命的影響程度。

兵器科學(xué)與技術(shù)；可靠壽命；疲勞壽命；概率功能度量；靈敏度分析

0 引言

疲勞破壞是機械構(gòu)件最常見的一種失效形式，受許多隨機因素的影響，疲勞壽命分散性較大［1］，因此，在進(jìn)行疲勞壽命預(yù)計時，往往需要同時進(jìn)行可靠性計算。工程中疲勞可靠性分析主要包括兩方面內(nèi)容:一是分析構(gòu)件在給定壽命下的可靠度或失效概率，通常利用應(yīng)力-強度干涉模型、功能函數(shù)法、模擬法等［2-3］；二是預(yù)計構(gòu)件給定可靠度或失效概率下的壽命［1，4］，即可靠壽命，這往往是設(shè)計師更為關(guān)注的問題。

現(xiàn)有的可靠壽命預(yù)計方法，一般是假定壽命服從威布爾分布、對數(shù)正態(tài)分布等，計算分布參數(shù)，利用概率分布函數(shù)確定可靠壽命［5-6］；或者是進(jìn)行多次給定壽命的可靠度計算，反復(fù)迭代，直至結(jié)果逼近給定可靠度值，從而得到可靠壽命值［7］，計算量較大。功能度量法常用于評估可靠性優(yōu)化設(shè)計中的概率約束，即給定概率下的功能函數(shù)值，是一種計算高效、穩(wěn)定的可靠性算法［8-10］。

本文在常規(guī)疲勞壽命預(yù)計方法的基礎(chǔ)上，考慮影響結(jié)構(gòu)疲勞壽命因素的隨機性，利用功能度量法的思路，通過將疲勞壽命表示為隨機變量的函數(shù)，把可靠壽命的預(yù)計問題轉(zhuǎn)化為給定概率下功能函數(shù)值的計算問題，也即為概率功能度量的求解，并在功能度量法的基礎(chǔ)上提出了計算可靠壽命及靈敏度分析方法，為工程應(yīng)用提供一種適用的估算方法。

1 疲勞壽命預(yù)計的工程方法

累積損傷理論是疲勞壽命預(yù)計的基礎(chǔ)，工程上一般采用Miner線性累積損傷理論，疲勞壽命（單位為載荷譜的重復(fù)次數(shù)）表示為

式中:D為臨界損傷因子，較多場合取1；ni為第i級應(yīng)力水平的循環(huán)次數(shù)（或在總循環(huán)次數(shù)中的百分比）；Ni為第i級應(yīng)力水平對應(yīng)的疲勞壽命；k為載荷譜的級數(shù)，當(dāng)k=1時，載荷譜為恒幅譜，直接計算N1即可。各級應(yīng)力水平對應(yīng)的疲勞壽命預(yù)計方法，包括應(yīng)力壽命法、應(yīng)變壽命法、能量法和場強法、裂紋擴(kuò)展壽命法等，下面將簡要介紹常用的3種方法:

1）應(yīng)力壽命法［11］。應(yīng)力壽命法以材料或零件應(yīng)力-壽命（S-N）曲線為基礎(chǔ)，主要用于高周疲勞失效。S-N曲線在雙對數(shù)坐標(biāo)系中通常用兩條直線和5個參數(shù)近似表示，如圖1所示。NS1為第1壽命點的循環(huán)次數(shù)，一般為104，S1為對應(yīng)的應(yīng)力幅值極限；NS2為第2壽命點的循環(huán)次數(shù)，S2為對應(yīng)的應(yīng)力幅值極限，也即疲勞極限σ-1.K2為循環(huán)次數(shù)超過第2壽命點的直線斜率。以σei表示第i級應(yīng)力水平的應(yīng)力幅值，則Ni與σei的函數(shù)關(guān)系為

圖1　材料的S-N曲線示意圖Fig.1 S-N curve of material

當(dāng)應(yīng)力水平的均值不為0時，采用Goodman模型等進(jìn)行平均應(yīng)力修正，并根據(jù)構(gòu)件的應(yīng)力集中情況、尺寸、表面粗糙度、加工方式等，修正材料S-N曲線。

2）應(yīng)變壽命法。應(yīng)變壽命法是以Manson-Coffin應(yīng)變-壽命關(guān)系式為基礎(chǔ)，適應(yīng)于低周疲勞失效。在計算壽命之前，通過諾伯法、有限元方法等計算結(jié)構(gòu)的局部應(yīng)力、應(yīng)變歷程，確定每個循環(huán)的應(yīng)力均值σm和應(yīng)變幅值εa，通過應(yīng)變-壽命關(guān)系，如Morrow修正式［11］，計算壽命N:

式中:σ′f為疲勞強度系數(shù)；Ε為彈性模量；ε′f為疲勞延性系數(shù)；b為疲勞強度指數(shù)；c為疲勞延性指數(shù)。由于（3）式對于N與其他參數(shù)之間的關(guān)系是一個隱式函數(shù)，無法直接寫出表達(dá)式，需要先求解非線性方程，才能得到壽命值。

3）裂紋擴(kuò)展壽命方法。對于含裂紋的構(gòu)件，根據(jù)斷裂力學(xué)理論，由Paris公式，得裂紋擴(kuò)展壽命的計算式［12］為

式中:ac為臨界裂紋尺寸；a0為初始裂紋尺寸；C、m為材料參數(shù)；Δσi為結(jié)構(gòu)件所受的應(yīng)力范圍；Y為裂紋的行位和幾何性質(zhì)修正系數(shù)。ac的計算公式為

式中:K1c為斷裂韌性。由于裂紋長度擴(kuò)展是一個連續(xù)過程，對于變幅應(yīng)力下的壽命估算，以循環(huán)塊計的裂紋擴(kuò)展壽命為

2 疲勞可靠壽命模型

根據(jù)疲勞強度理論，構(gòu)件應(yīng)力、應(yīng)變受載荷、幾何尺寸、材料力學(xué)性能的波動，是一組隨機量［13］；材料的S-N曲線、應(yīng)變-壽命曲線及裂紋擴(kuò)展參數(shù)具有分散性，實際結(jié)構(gòu)疲勞壽命亦具有較大的分散性。在實際應(yīng)用中，應(yīng)統(tǒng)計各參數(shù)的分布類型、均值和標(biāo)準(zhǔn)差，將分散性較大的參數(shù)設(shè)為隨機變量x，分散性較小的參數(shù)不考慮其隨機性，設(shè)為確定性變量v，歸納（1）式～（4）式4個疲勞壽命計算式，則壽命N是載荷、幾何尺寸、材料性能參數(shù)等隨機變量x、確定性變量v的函數(shù)，用N（x，v）統(tǒng)一表示，其中x為各隨機變量組成的向量，v為各確定性變量組成的向量。

若已知壽命的概率分布，由概率密度函數(shù)，即可得到給定可靠度R時的壽命，即可靠壽命NR，如圖2所示。

圖2　可靠壽命與壽命概率密度圖Fig.2 Probability density of life and reliable life

因此，構(gòu)件的疲勞可靠度表示為壽命N（x，v）大于給定壽命NR的概率，即

式中:f（x）為x的聯(lián)合概率密度函數(shù)；FN（·）為壽命N（x，v）的累積分布函數(shù)。此外，使用矩法計算可靠度時，往往將原始隨機變量向量x變換為相互獨立的標(biāo)準(zhǔn)正態(tài)分布向量u，記u=T（x），在獨立標(biāo)準(zhǔn)正態(tài)空間中壽命函數(shù)表示為Nu（u，v），則構(gòu)件的疲勞可靠度也即為

式中:g（u）為標(biāo)準(zhǔn)正態(tài)分布向量u的聯(lián)合概率密度函數(shù)。若為壽命的逆分布函數(shù)，則可靠壽命為

式中:β為可靠度指標(biāo)，即R=φ（β）.此時，NR可稱為壽命函數(shù)N（x，v）的概率功能度量，也即為待求的可靠壽命。

3 疲勞可靠壽命預(yù)計

前面給出了疲勞壽命函數(shù)，均為非線性函數(shù)，面對多個隨機變量，壽命的概率密度函數(shù)、逆分布函數(shù)難以直接獲得，一般需要借助近似的方法計算可靠壽命。功能度量法常用于可靠性優(yōu)化設(shè)計中概率約束的評估，經(jīng)過工程實踐發(fā)現(xiàn)，相對于可靠度指標(biāo)法相比，具有更高效、穩(wěn)定和較少依賴于隨機變量的概率分布類型等特點［14-15］。本文將借鑒該方法的思路，預(yù)計疲勞可靠壽命。

為了更好理解如何利用功能度量法求解可靠壽命，先簡要介紹一次2階矩法中可靠度系數(shù)的幾何意義。在獨立標(biāo)準(zhǔn)正態(tài)空間中，可靠度指標(biāo)的意義是原點至極限狀態(tài)曲線（或曲面）的最短距離，如圖3所示。若給定壽命NR，由（7）式知Nu（u，v）= NR為極限狀態(tài)曲線，則通過一次可靠度方法等能夠計算該曲線離原點的距離為β=2，可靠度為φ（2）.

對于可靠壽命的求解，即為給定可靠度值φ（β），確定離原點最短距離為β的等壽命曲線（或曲面）。以兩個隨機變量為例，如圖3所示，Nu（u，v）分別為N1、NR、N2表示一簇等壽命曲線，利用計算可靠度的方法，將會得到對應(yīng)的一系列β值；反之，作半徑為β=2的圓弧，尋找與該圓弧相切的等壽命曲線，切點為u*，即標(biāo)準(zhǔn)正態(tài)空間內(nèi)半徑為β的圓弧面上函數(shù)值最小點，則圖中NR為給定可靠度φ（2）的壽命。因此，在標(biāo)準(zhǔn)正態(tài)空間中求解疲勞可靠壽命，相當(dāng)于求解疲勞可靠度的逆命題。

由此，在獨立標(biāo)準(zhǔn)正態(tài)空間，對于給定可靠度指標(biāo)βt，疲勞可靠壽命的計算表示為如下的數(shù)學(xué)優(yōu)化問題，即確定u*滿足（10）式:

對于這一具有球面約束的優(yōu)化問題，由優(yōu)化問題的Karush-Kuhn-Tucker條件，最優(yōu)點u*滿足:

圖3　標(biāo)準(zhǔn)正態(tài)空間中可靠壽命的幾何意義Fig.3 Illustration of the calculation of reliable life in the standard normal space

對于可靠壽命，給定可靠度一般大于0.5，故βt＞0，代入（11）式。由此得到優(yōu)化計算的迭代公式為

式中:梯度向量ΔNu（uk，v）為壽命函數(shù)在uk點處對各隨機變量的偏導(dǎo)數(shù)，迭代初值可以取標(biāo)準(zhǔn)正態(tài)空間中各隨機變量的均值，當(dāng)‖uk+1‖-‖uk‖小于容許誤差時收斂，此時求得的即為設(shè)計點u*= uk+1，將u*代入壽命函數(shù)就得到可靠壽命，即NR= Nu（u*，v）.

上述可以看出，利用本文方法求解可靠壽命，只需進(jìn)行一次優(yōu)化迭代計算后，就可以直接獲得給定可靠度下的可靠壽命，大幅降低計算量，適合于工程應(yīng)用。

4 疲勞可靠壽命靈敏度分析

4.1可靠壽命對隨機變量的靈敏度分析

在計算可靠壽命的每一步迭代時需要計算的梯度向量ΔNu（u，v），在設(shè)計點u*處可靠壽命NR對隨機變量xi的均值μi偏導(dǎo)數(shù)表示為

式中:

式中:y為隨機變量向量轉(zhuǎn)換為相關(guān)標(biāo)準(zhǔn)正態(tài)空間的向量；L為隨機變量相關(guān)系數(shù)矩陣的Cholesky分解矩陣，即L與相關(guān)系數(shù)矩陣ρ的關(guān)系為ρ=L·LT；σi為隨機變量xi在設(shè)計點u*的當(dāng)量標(biāo)準(zhǔn)差［16-17］，其計算公式為

式中:Φ（·）、φ（·）分別是標(biāo)準(zhǔn)正態(tài)分布函數(shù)和密度函數(shù)，F(xiàn)i（·）、fi（·）分別是隨機變量xi的分布函數(shù)和密度函數(shù)。

可靠壽命對隨機變量的均值靈敏度為

同理，可靠壽命對隨機變量的標(biāo)準(zhǔn)差靈敏度為

4.2可靠壽命對確定性變量的靈敏度分析

利用差分方法，在設(shè)計點u*處，設(shè)定設(shè)計變量vi賦微小的步長Δvi，則可靠壽命對確定性設(shè)計變量的靈敏度為

通過可靠壽命靈敏度分析能夠得到可靠壽命對各設(shè)計變量的靈敏度，均值靈敏度和確定性量靈敏度分別反映了隨機變量均值和確定性量對可靠壽命的影響。靈敏度值正負(fù)可以判斷隨機變量均值和確定性量對可靠壽命的影響趨勢，而靈敏度值大小則反應(yīng)了各設(shè)計變量對可靠壽命的影響程度，對結(jié)構(gòu)設(shè)計改進(jìn)時，著重考慮靈敏度值較大的變量，更容易實現(xiàn)結(jié)構(gòu)壽命設(shè)計。

5 計算實例

某型車輛行動系統(tǒng)扭力軸在工程使用發(fā)生過多起斷裂故障，由于加工工藝、磨損、腐蝕等原因，扭力軸表面容易存在小裂紋，通過對多根扭力軸表面裂紋測量統(tǒng)計，裂紋尺寸a0服從正態(tài)分布，均值為2mm，標(biāo)準(zhǔn)差0.3mm.扭力軸在工作過程中主要承受反復(fù)的扭矩作用，容易產(chǎn)生疲勞斷裂。表1所示為扭矩幅值的塊譜數(shù)據(jù)，總平均扭矩為 8280N·m，當(dāng)量行駛里程為350km.試分析該扭力軸的耐久性，并提出改進(jìn)措施。

表1　扭力軸扭矩譜Tab.1 Torque spectra of a torsion shaft

對于含裂紋的構(gòu)件進(jìn)行壽命預(yù)計，主要采用裂紋擴(kuò)展壽命法，根據(jù)（4）式，需先求解扭力軸的應(yīng)力幅值。由扭力軸的扭矩幅值計算扭力軸光滑圓柱的名義切應(yīng)力幅值

式中:d為扭力軸光滑圓柱部分直徑，服從均值為54mm、標(biāo)準(zhǔn)差為1.6mm的正態(tài)分布。

扭力軸的應(yīng)力集中部位為過渡圓弧至花鍵的結(jié)合處，有效應(yīng)力集中系數(shù)為

式中:R為由光滑圓柱向頭部花鍵過渡的圓弧半徑；kt為扭力軸軸頸的理論扭轉(zhuǎn)應(yīng)力集中系數(shù)，與R/d有關(guān)，查應(yīng)力集中手冊為1.05；σs、σb分別為材料的屈服強度、抗拉強度。

扭力軸材料為高強度合金鋼，抗拉強度σb為1553MPa，抗拉屈服強度σs為1374MPa；斷裂韌度K1c服從對數(shù)正態(tài)分布，均值、標(biāo)準(zhǔn)差分別為經(jīng)試驗統(tǒng)計，裂紋擴(kuò)展參數(shù)C服從均值為2.08×10-11、標(biāo)準(zhǔn)差為0.28×10-11的正態(tài)分布，裂紋擴(kuò)展參數(shù)m服從均值為2.44、標(biāo)準(zhǔn)差為0.05的正態(tài)分布，忽略C、m相關(guān)性。

利用（4）式～（6）式分別預(yù)計給定可靠度為0.5、0.9、0.95、0.99時扭力軸的壽命；同時，為對比本文方法預(yù)計結(jié)果的精度，利用蒙特卡洛模擬法（簡稱模擬法），即通過計算機程序按各參數(shù)隨機分布生成樣本組，依次計算各樣本組對應(yīng)的可靠壽命，為保證模擬法的計算精度，樣本量設(shè)置為10萬次。并且計算本文方法預(yù)計結(jié)果相對于模擬法結(jié)果的誤差，一并列舉于表2.

表2　可靠壽命的預(yù)計結(jié)果Tab.2 Predicted results of reliable life

從表2中結(jié)果可以看出:利用本文方法求解的可靠壽命結(jié)果，與模擬法結(jié)果非常接近；但計算量遠(yuǎn)小于模擬法的計算量，對于本例，各次求解的迭代次數(shù)均不超過5次，計算時間約為1s；而10萬次蒙特卡洛模擬計算時間約480s，因此該方法的計算效率高是顯著的。尤其對于調(diào)用CAE軟件計算構(gòu)件應(yīng)力、應(yīng)變，本文的方法求解可靠壽命的優(yōu)勢更加凸顯。

針對表2中第2列結(jié)果，利用一次2階矩法計算給定壽命為4355km的扭力軸疲勞可靠度，結(jié)果為0.897，與0.90相接近，再次驗證了本文方法計算可靠壽命有較高的精度。

對于車輛，通常以B10作為可靠性指標(biāo)，即可靠度為0.90的壽命，為4355km，低于工程要求，需要改進(jìn)設(shè)計。利用（17）式～（19）式分析該可靠壽命對隨機變量均值、標(biāo)準(zhǔn)差的靈敏度，如圖4、圖5分別為扭力軸可靠壽命對主要隨機變量均值、標(biāo)準(zhǔn)差的靈敏度直方圖。

由圖4和圖5可以看出，扭力軸直徑d、裂紋擴(kuò)展參數(shù)C和m、初始裂紋尺寸a0的均值、標(biāo)準(zhǔn)差對可靠壽命的靈敏度值較大，其影響較敏感。在扭力軸改進(jìn)設(shè)計時，適當(dāng)增大扭力軸直徑，選取C、m較小的材料，控制初始裂紋大小，有助于提高扭力軸的可靠壽命。如將扭力軸直徑均值增加2mm，同時加工上下偏差控制在±1.0mm內(nèi)，經(jīng)過重新計算，扭力軸的B10提高至7543km，壽命得到明顯提升。

6 結(jié)論

本文將功能度量法引入疲勞壽命分析，提出了一種可靠壽命及其靈敏度分析方法。該方法將可靠壽命的計算問題轉(zhuǎn)換為一個球面約束的優(yōu)化問題，只需進(jìn)行一次概率功能度量的求解即可獲得可靠壽命，避免了可靠壽命預(yù)計之前需假設(shè)壽命服從某種分布。該方法給出了可靠壽命對隨機變量均值、標(biāo)準(zhǔn)差以及確定量的靈敏度計算公式。

通過對某型車輛扭力軸的可靠壽命分析，并與蒙特卡洛模擬法、一次2階矩法等計算結(jié)果相比較，表明基于概率功能度量的疲勞可靠壽命預(yù)計方法具有較好的效率和精度，靈敏度結(jié)果有效地反映出各隨機變量對疲勞可靠壽命的影響程度。

圖4　可靠壽命對隨機變量均值的靈敏度比較Fig.4 Comparison of sensitivity factors for average value of each random variable

圖5　可靠壽命對隨機變量標(biāo)準(zhǔn)差的靈敏度比較Fig.5 Comparison of sensitivity factors for standard deviation of each random variable

（References）

［1］ 張鈺，王丹，張風(fēng)和，等.一種新的疲勞可靠壽命計算方法［J］.東北大學(xué)學(xué)報:自然科學(xué)版，2000，21（1）:42-45. ZHANG Yu，WANG Dan，ZHANG Feng-he，et al.New model on prediction of reliable fatigue life［J］.Journal of Northeastern University:Natural Science，2000，21（1）:42-45.（in Chinese）

［2］ Rackwitz R.Reliability analysis-a review and some perspectives［J］.Structural Safety，2001，23（4）:365-395.

［3］ Lopez I，Sarigul-Klijin N.A review of uncertainty in flight vehicle structural damage monitoring，diagnosis and control:challenges and opportunities［J］.Progress in Aerospace Sciences，2010，46（7）:247-273.

［4］ Wirsching T Y T，Martin W S.Advanced fatigue reliability analysis［J］.International Journal of Fatigue，1991，13（5）:89-394.

［5］ 謝里陽，任俊剛，吳寧祥，等.復(fù)雜結(jié)構(gòu)部件概率疲勞壽命預(yù)測方法與模型［J］.航空學(xué)報，2015，36（8）:2688-2695. XIE Li-yang，REN Jun-gang，WU Ning-xiang，et al.Probabilistic fatigue life prediction method and modeling for complex strctural parts［J］.Acta Aeronautica et Astronautica Sinica，2015，36（8）:2688-2695.（in Chinese）

［6］ Zuo F J，Zhu S P，Gao H Y，et al.Stochastic fatigue life and reliability prediction based on residual strength［J］.Journal of Shanghai Jiaotong University（Science），2015，20（3）:331-337.

［7］ Guan X F，He J J，Rasselkorde E M.Probabilistic fatigue life prediction and structural reliability evaluation of turbine rotors integrating an automated ultrasonic inspection system［J］.Journal of Nondestructive Evaluation，2014，33（1）:51-61.

［8］ Tu J，Choi K K.A new study on reliability based design optimization［J］.Journal of Mechanical Design，1999，121（4）:557-564.

［9］ Ramu P，Qu X，Youn B，et al.Inverse reliability measures and reliability-based design optimization［J］.International Journal of Reliability and Safety，2006，1（1/2）:187-205.

［10］ Li G，Meng Z，Hu H.An adaptive hybrid approach for reliability-based design optimization［J］.Structural&Multidisciplinary Optimization，2015，51（5）:1051-1065.

［11］ Suresh S.Fatigue of materials［M］.Cambridge:Cambridge University Press，1998:148-233.

［12］ 楊濤，王紅巖.裝甲車體疲勞裂紋擴(kuò)展壽命預(yù)測方法研究［J］.兵工學(xué)報，2010，31（2）:129-134. YANG Tao，WANG Hong-yan.Research on fatigue crack growth life prediction technique for the armored hull［J］.Acta Armamentarii，2010，31（2）:129-134.（in Chinese）

［13］ 徐灝.概率疲勞［M］.沈陽:東北大學(xué)出版社，1994:236-314. XU Hao.Probability fatigue［M］.Shenyang:Northeastern University Press，1994:236-314.（in Chinese）

［14］ Younes A，Alaa C.Benchmark study of numerical methods for reliability-based design optimization［J］.Structural and Multidisciplinary Optimization，2010，41（2）:277-294.

［15］ Yi P，Cheng G D，Jiang L.A sequential approximate programming strategy for performance-measure-based probabilistic structural design optimization［J］.Structural Safety，2008，30（2）: 91-109.

［16］ Rahman S，Wei D.Design sensitivity and reliability-based structural optimization by univariate decomposition［J］.Struct Multidisc Optim，2008，35（3）:245-261.

［17］ Au S.Reliability-based design sensitivity by efficient simulation［J］.Computers&Structures，2005，83（14）:1048-1061.

Reliable Fatigue Life Analysis Method Based on Probabilistic Performance Measure

LIU Qin1，2，SUN Zhi-li1，QIAN Yun-peng2，LIU Ying2
（1.School of Mechanical Engineering and Automation，Northeastern University，Shenyang 110819，Liaoning，China；2.Ordnance Science and Research Academy of China，Beijing 100089，China）

Based on the frequently-used engineering methods for fatigue life prediction and the performance measure approach，a reliable fatigue life model is built，and a set of methods for reliable life prediction and sensitivity analysis are proposed.By expressing fatigue life as the function of random variables，the reliable life prediction problem is converted into a spherical constraint optimization problem which can be solved using performance measure approach，and then the reliable life can be gotten directly.Moreover，the sensitivity factors for each random variable and deterministic variable can be calculated by using some intermediate variables in life prediction.The proposed method is demonstrated with a torsion shaft example.The results show that the proposed method has high efficiency and precision，and the different parameters have much different influences on reliable life.

ordnance science and technology；reliable life；fatigue life；probabilistic performance measure；sensitivity analysis

TB114.3

A

1000-1093（2016）08-1530-06

10.3969/j.issn.1000-1093.2016.08.027

2016-01-06

總裝備部預(yù)先研究項目（51319010103）；國防技術(shù)基礎(chǔ)項目（Z092012B001）

劉勤（1981—），男，副研究員，博士研究生。E-mail:qinlow@126.com；孫志禮（1957—），男，教授，博士生導(dǎo)師。E-mail:zhlsun@neu.edu.cn