王春生，張 磊，東海杰

?

基于NESSUS的杜瓦結(jié)構(gòu)可靠性分析

王春生，張 磊，東海杰

（華北光電技術(shù)研究所，北京 100015）

為了提高杜瓦結(jié)構(gòu)的可靠性水平，本文提出了基于NESSUS（Numerical Evaluation of Stochastic Structures Under Stress）軟件的杜瓦結(jié)構(gòu)可靠性分析方法。此方法將杜瓦組件的載荷、材料性能、幾何尺寸等參數(shù)均作為具有某種分布特性的隨機(jī)變量，通過概率分析，可獲得定量的杜瓦結(jié)構(gòu)可靠度值以及各個設(shè)計參數(shù)的概率靈敏度值，并由此確定關(guān)鍵影響因素，為提高杜瓦結(jié)構(gòu)的可靠性提供設(shè)計依據(jù)。本文以杜瓦冷臺支撐結(jié)構(gòu)為例進(jìn)行了結(jié)構(gòu)可靠性分析，分析結(jié)果表明此方法可以顯著提高杜瓦結(jié)構(gòu)的可靠性水平。

結(jié)構(gòu)可靠性；NESSUS；可靠性仿真分析；杜瓦結(jié)構(gòu)

0 引言

杜瓦為紅外焦平面探測器提供良好的低溫工作環(huán)境以及光、機(jī)、電、熱傳輸通道，是紅外焦平面探測器的封裝和保護(hù)裝置，也是紅外焦平面探測器完成工程化、實用化的重要環(huán)節(jié)[1]。因此，杜瓦結(jié)構(gòu)的可靠性對整個紅外焦平面探測器組件的可靠性有著重要的影響。杜瓦中包含十幾個關(guān)鍵零部件，其材料特性及所承受的載荷都存在一定的離散性。例如，在測量材料的物理參數(shù)時必須將實際材料制成一定規(guī)格的試件，施加特定載荷并記錄變形情況，從而計算出物理參數(shù)。但是，通過多次實驗獲得的物理參數(shù)值并不是一個唯一值，而是具備某種分布特性的一組值，表現(xiàn)出一定的分散性。這種分散性將導(dǎo)致設(shè)計結(jié)果的不確定性，帶來一定的風(fēng)險。

目前杜瓦結(jié)構(gòu)的設(shè)計與分析采用的是確定性方法。該方法隱含了下列假設(shè)：①對于設(shè)計參數(shù)可以實現(xiàn)完備的測量；②設(shè)計參數(shù)與結(jié)構(gòu)響應(yīng)間存在確定性關(guān)系；③在給定設(shè)計方案下，其制造、安裝和運行中各種因素可以實現(xiàn)完全控制。顯然在實際工程中，這種假設(shè)是很難實現(xiàn)的。因此，確定性方法的設(shè)計結(jié)果存在風(fēng)險。為此，在確定性設(shè)計中引入一定的安全系數(shù)來降低不確定性風(fēng)險。但安全系數(shù)法只是一種簡單的經(jīng)驗方法，無法準(zhǔn)確定量地評價可靠程度，也不能識別影響組件可靠性的關(guān)鍵參數(shù)或關(guān)鍵失效模式，不利于風(fēng)險控制和優(yōu)化設(shè)計[2-3]。同時，它也無法從理論上解釋按安全系數(shù)法設(shè)計為安全的系統(tǒng)，在使用中卻發(fā)生失效的現(xiàn)象。另外，安全系數(shù)法通常采用最壞情況假設(shè)，會帶來較大的冗余設(shè)計，造成結(jié)構(gòu)體積、重量增加，不適合機(jī)載、星載、彈載等特殊應(yīng)用場合的需求。因此現(xiàn)行的基于安全系數(shù)法的確定性設(shè)計很難滿足紅外焦平面探測器杜瓦結(jié)構(gòu)的可靠性保障要求。隨著紅外探測應(yīng)用領(lǐng)域的不斷拓展，對紅外探測器組件可靠性的要求日益提高，需要開展系統(tǒng)化的杜瓦結(jié)構(gòu)可靠性設(shè)計和仿真分析工作。

近年來國內(nèi)外相關(guān)機(jī)構(gòu)針對紅外探測器杜瓦結(jié)構(gòu)進(jìn)行了可靠性評價與加速壽命實驗等方面的研究，文獻(xiàn)[4]中介紹法國SOFRADIR公司的杜瓦真空壽命已達(dá)到20年，文獻(xiàn)[5]分析了SOFRADIR公司的杜瓦在振動、機(jī)械沖擊以及溫度沖擊下的失效模式。國內(nèi)近年來通過加速壽命實驗對杜瓦真空壽命進(jìn)行了研究，文獻(xiàn)[6-8]均采用阿列尼烏斯模型對紅外杜瓦的真空壽命進(jìn)行了研究。但是迄今為止，在設(shè)計和仿真分析領(lǐng)域，杜瓦結(jié)構(gòu)設(shè)計仍然采用的是確定性方法，無法滿足工程實際的需求。本文結(jié)合杜瓦結(jié)構(gòu)的特點和可靠性要求，提出了基于結(jié)構(gòu)可靠性仿真軟件NESSUS（Numerical Evaluation of Stochastic Structures Under Stress）的杜瓦結(jié)構(gòu)可靠性仿真分析方法，并以冷臺面支撐結(jié)構(gòu)為例進(jìn)行了可靠度計算和設(shè)計參數(shù)靈敏度分析，分析結(jié)果對杜瓦結(jié)構(gòu)的可靠性設(shè)計具有重要的指導(dǎo)意義。

1 基于NESSUS的結(jié)構(gòu)可靠性仿真分析方法

1.1 結(jié)構(gòu)可靠性定量分析

工程實際中，杜瓦結(jié)構(gòu)的強(qiáng)度和應(yīng)力都是呈某種分布的隨機(jī)變量，因此，以廣義應(yīng)力-強(qiáng)度干涉模型作為可靠性仿真分析的理論基礎(chǔ)。圖1表示出了應(yīng)力-強(qiáng)度干涉模型的具體形式，將應(yīng)力和強(qiáng)度的分布曲線畫在同一坐標(biāo)系中，橫坐標(biāo)代表應(yīng)力和強(qiáng)度的數(shù)值，縱坐標(biāo)代表二者的概率密度。

圖1中應(yīng)力和強(qiáng)度相互交疊的區(qū)域即為可能發(fā)生失效的區(qū)域，在此區(qū)域內(nèi)強(qiáng)度可能小于應(yīng)力。結(jié)構(gòu)的可靠度就是強(qiáng)度大于應(yīng)力的概率，由式(1)所示：

＝(＜)＝[(－)＞0] (1)

式中：為結(jié)構(gòu)可靠度；為廣義應(yīng)力；為廣義強(qiáng)度。

理論上，結(jié)構(gòu)可靠度可由公式(2)計算得出：

但在工程實際中式(2)是很難求得解析解的。結(jié)構(gòu)可靠性分析就是通過對影響可靠性的各種性能指標(biāo)進(jìn)行概率計算，得到定量的可靠性值，并以此為依據(jù)對結(jié)構(gòu)設(shè)計進(jìn)行改進(jìn)，滿足可靠性要求。同時，在量化分析的基礎(chǔ)上得到較小的零件尺寸、體積和重量，從而在研發(fā)過程中平衡設(shè)計保守和設(shè)計不足的矛盾，解決安全系數(shù)不能綜合量化可靠性的問題。

1.2 靈敏度分析

在概率分析中，靈敏度表示概率隨分布參數(shù)（如均值和方差）的改變而變化的程度，由下式定義：

＝?/?(3)

式中：為結(jié)構(gòu)的失效概率；為隨機(jī)變量的分布參數(shù)。

靈敏度分析可定量地給出影響結(jié)構(gòu)可靠性的主要因素，可用于指導(dǎo)結(jié)構(gòu)設(shè)計改進(jìn)和結(jié)構(gòu)故障預(yù)測[9]。若某參數(shù)對可靠性影響較大，則應(yīng)采取相應(yīng)的優(yōu)化設(shè)計和工藝控制手段，降低其不確定性。反之，若某參數(shù)對結(jié)構(gòu)可靠性的影響不顯著，則可在下一步的可靠性改進(jìn)設(shè)計中將其作為確定量值處理，以減少隨機(jī)變量的數(shù)目，提高計算效率。

1.3 基于NESSUS的可靠性仿真分析

NESSUS是美國西南研究院的結(jié)構(gòu)可靠性分析軟件，它本身具有非常豐富的概率算法，并且通過調(diào)用有限元分析軟件來共同完成結(jié)構(gòu)的可靠性計算。該軟件可定量評估零件、系統(tǒng)的可靠度，識別重要隨機(jī)變量，為結(jié)構(gòu)可靠性設(shè)計提供依據(jù)，并由此開發(fā)出比傳統(tǒng)確定性設(shè)計方法更經(jīng)濟(jì)、更有效、更可靠的產(chǎn)品?；贜ESSUS的結(jié)構(gòu)可靠性仿真分析流程如圖2所示[10]，分析是自上而下逐步進(jìn)行的，首先進(jìn)行問題描述和隨機(jī)變量定義，然后是模型定義、分析類型選擇、執(zhí)行分析，最后顯示分析結(jié)果。

從可靠性的角度分析，紅外探測器組件杜瓦結(jié)構(gòu)的可靠性工程具有如下特點：

1）受力條件復(fù)雜

杜瓦結(jié)構(gòu)件除受到靜載荷、振動載荷外，還受到極大的熱應(yīng)力作用，載荷的復(fù)雜性使得可靠性仿真分析中的極限狀態(tài)方程無法求得解析解，必須借助有限元模型進(jìn)行數(shù)值求解。而概率分析中需對有限元模型中的相關(guān)參數(shù)進(jìn)行擾動處理，并多次進(jìn)行計算。該過程極為復(fù)雜，并且計算量極大，必須采用專業(yè)的概率分析軟件進(jìn)行。

圖2 基于NESSUS的結(jié)構(gòu)可靠性仿真分析流程

2）存在關(guān)鍵結(jié)構(gòu)件

杜瓦內(nèi)存在多個關(guān)鍵的結(jié)構(gòu)件，各關(guān)鍵結(jié)構(gòu)件都有其自身的工作特點和可靠性指標(biāo)，對關(guān)鍵件進(jìn)行定量的可靠性仿真分析是提高杜瓦整體可靠性水平的重要手段。同時，通過靈敏度分析可以明確對可靠性影響最大的關(guān)鍵參數(shù)，找出最有效的改進(jìn)方法。例如，冷臺面支撐結(jié)構(gòu)作為紅外探測器的承載平臺，是杜瓦結(jié)構(gòu)的關(guān)鍵件，需對其進(jìn)行定量的可靠性分析。該支撐結(jié)構(gòu)為薄壁、細(xì)長的懸臂結(jié)構(gòu)，其材料特性、外載荷以及零件尺寸的變化均影響其可靠性，可以通過靈敏度分析找出影響最大的參數(shù)，并加以改進(jìn)。

以下，本文將以杜瓦的冷臺面支撐結(jié)構(gòu)為例，詳細(xì)論述基于NESSUS的杜瓦結(jié)構(gòu)可靠性仿真分析方法。

2 建立有限元模型

如前所述，杜瓦結(jié)構(gòu)的極限狀態(tài)方程難以求得解析解，必須借助有限元模型進(jìn)行數(shù)值求解。有限元模型在可靠性仿真分析中所處的地位如圖3所示。

本文采用ANSYS軟件進(jìn)行有限元建模與求解。NESSUS通過ANSYS的參數(shù)化功能對ANSYS的模型參數(shù)進(jìn)行隨機(jī)參數(shù)映射，然后，由NESSUS修改指定輸入文件中相應(yīng)的概率統(tǒng)計值，并遞交給ANSYS進(jìn)行求解。最后，NESSUS通過解析指定的輸出文件，得到響應(yīng)值。

冷臺面支撐結(jié)構(gòu)的有限元分析模型如圖3所示，約束條件為施加在法蘭端面的固定約束，載荷為隨機(jī)振動載荷，模型的關(guān)鍵參數(shù)為支撐柱的內(nèi)徑、外徑、長度以及材料的彈性模量和泊松比。

圖3 結(jié)構(gòu)可靠性分析解決方案示意圖

3 建立可靠性仿真分析模型

可靠性仿真分析是進(jìn)行可靠性改進(jìn)設(shè)計的先決條件，只有通過重要度、靈敏度分析后才能確定改進(jìn)的方向，并采取針對性的改進(jìn)措施，例如優(yōu)化設(shè)計參數(shù)、加強(qiáng)工藝過程控制、提高加工一致性等，進(jìn)而實現(xiàn)可靠性的提升。要對杜瓦結(jié)構(gòu)進(jìn)行可靠性分析，首先要明確其功能、主要載荷以及與之相對應(yīng)的主要故障模式，并且要找到影響結(jié)構(gòu)性能與故障的主要設(shè)計參數(shù)。

3.1 極限狀態(tài)函數(shù)定義

極限狀態(tài)指的是應(yīng)力與強(qiáng)度相等時的臨界狀態(tài)，當(dāng)超過此狀態(tài)時結(jié)構(gòu)就會發(fā)生失效。以數(shù)學(xué)公式表達(dá)為極限狀態(tài)函數(shù)：

＝－(3)

式中：為極限狀態(tài)函數(shù)；為廣義應(yīng)力；為廣義強(qiáng)度。

在本例中，冷臺面支持結(jié)構(gòu)承受隨機(jī)振動載荷的作用，并產(chǎn)生變形，由此引起探測器芯片相對于光軸發(fā)生位置變化。當(dāng)探測器的位置變化達(dá)到一定程度時將影響成像質(zhì)量，尤其是當(dāng)系統(tǒng)對MTF（Modulation Transfer Function，調(diào)制傳遞函數(shù)）要求比較高時，過大的位置變化將成為一個主要的失效模式。因此，本例中以探測器芯片位置變化小于7mm為極限狀態(tài)函數(shù)，如式(4)所示：

＝0.007－＝disp(i,o,x, Pr) (4)

式中：為極限狀態(tài)函數(shù)；為探測器芯片位置變化量；i為支撐柱內(nèi)徑；o為支撐柱外徑；為支撐柱長度；x為彈性模量；Pr為泊松比。

3.2 隨機(jī)變量定義

進(jìn)行可靠性仿真分析首先要確定設(shè)計變量及其分布類型。本文將冷臺面支撐結(jié)構(gòu)的材料特性（包括彈性模量x、泊松比Pr）以及幾何尺寸（包括內(nèi)徑i、外徑o、長度）均作為隨機(jī)變量，其分布類型及參數(shù)如表1所示。

3.3 仿真分析結(jié)果

在NESSUS軟件中選擇選擇概率分析方法為AMV＋方法，分析計算后得到概率累積分布函數(shù)和參數(shù)靈敏度分別如圖4、圖5所示。

由分析結(jié)果可見，在隨機(jī)振動作用下，冷臺面支撐結(jié)構(gòu)位移失效的可靠度為82%（置信度為95%），可見該結(jié)構(gòu)在該工況下無法滿足可靠性要求，需對設(shè)計方案加以改進(jìn)。

概率靈敏度分析表明彈性模量x是影響冷臺面支撐結(jié)構(gòu)可靠性的主要因素。因此，選用更高規(guī)格的材料，并嚴(yán)格控制材料的采購、標(biāo)識和使用環(huán)節(jié)，降低材料特性的離散程度。同時，由靈敏度分析可見，幾何參數(shù)的不確定性對可靠度的影響較小，在進(jìn)一步的分析中可將其視為確定性量，以簡化概率分析模型。

3.4 設(shè)計改進(jìn)

如表2所示，忽略幾何參數(shù)的影響，將材料彈性模量x的標(biāo)準(zhǔn)差縮小，重新進(jìn)行概率分析。

重新計算的累積分布結(jié)果如圖6所示。

由分析結(jié)果可見，僅改進(jìn)一個設(shè)計參數(shù)（彈性模量x的標(biāo)準(zhǔn)差），便可將可靠度提升到99%，由此可充分體現(xiàn)出基于NESSUS的可靠性仿真分析方法的優(yōu)勢。

4 結(jié)論

本文結(jié)合紅外探測器杜瓦結(jié)構(gòu)的特點和可靠性要求，提出了基于NESSUS的杜瓦結(jié)構(gòu)可靠性仿真分析方法，并以冷臺面支撐結(jié)構(gòu)為例進(jìn)行了可靠性仿真分析。該方法將載荷、材料性能、幾何尺寸等設(shè)計參數(shù)均視為具有某種概率分布的統(tǒng)計量，應(yīng)用概率統(tǒng)計方法與強(qiáng)度理論求出零部件的可靠度，從而更加真實地反映可靠性水平，有效解決了傳統(tǒng)確定性設(shè)計不能量化可靠度的問題。采用該方法還可以通過重要度、靈敏度分析識別影響可靠度的主要因素和主要失效模式，進(jìn)而有針對性地加以改進(jìn)，有效提高杜瓦結(jié)構(gòu)的固有可靠性水平。

表1 原設(shè)計隨機(jī)變量表

圖4 累積分布函數(shù)

圖5 隨機(jī)變量靈敏度

表2 改進(jìn)設(shè)計隨機(jī)變量表

圖6 改進(jìn)后的累積分布函數(shù)

[1] 王小坤, 朱三根, 龔海梅. 長線列碲鎘汞紅外焦平面微型杜瓦的研究[J]. 激光與紅外, 2006, 36(11): 1047-1050.

WANG Xiaokun, ZHU Sangen, GONG Haimei. The study of the micro dewar for long linear HgCdTe IRFPA[J]., 2006, 36(11): 1047-1050.

[2] Moon Byung, Kang Gyung J U, Kang Beom Soo. Dynamic and reliability analysis of stochastic structure system using probabilistic finite element method[J]., 2004, 18(1): 125-135.

[3] 陳孝添. 在型號設(shè)計中應(yīng)用結(jié)構(gòu)可靠性設(shè)計方法[J]. 強(qiáng)度與環(huán)境, 2003, 30(2): 43-48.

CHEN Xiaotian. Applying the design of structural reliability to structure engineering[J]., 2003, 30(2): 43-48.

[4] Brenière X, Tribolet P. IR detectors design and approach for tactical applications with high reliability without maintenance[C]//, 2008,6940: 69400H-1-69400H-13.

[5] Molina M, Brenière X, Tribolet P. IR detector dewar and assemblies for stringent environmental conditions[C]//, 2007, 6542: 65422N-1- 65422N-11.

[6] 黃燕, 吳全信. 阿列尼烏斯模型在紅外微型金屬杜瓦的真空壽命實驗研究中的應(yīng)用[J]. 低溫與超導(dǎo), 2005, 33(2): 69-72.

HUANG Yan, WU Quanxin. The application of arrhenius model in the experimental study on the vacuum life of infrared miniature metal dewar[J]., 2005, 33(2): 69-72.

[7] 林日東, 劉偉, 王冠, 等. 紅外焦平面探測器杜瓦組件真空壽命分析[J]. 激光與紅外, 2011, 41(7): 779-783.

LIN Ridong, LIU Wei, WANG Guan, et al. Vacuum life analyses of infrared detector ＆ dewar assembly[J]., 2011, 41(7): 779-783.

[8] 張亞妮, 朱三根, 龔海梅. 微型杜瓦組件真空加速壽命實驗的研究[J].紅外, 2008, 29(3): 11-15.

ZHANG Yani, ZHU Sangen, GONG Haimei. Research on vacuum life of micro metal dewar[J]., 2008, 29(3): 11-15.

[9] 周紅, 劉永壽, 高宗戰(zhàn), 等. 某型導(dǎo)彈艙段連接結(jié)構(gòu)強(qiáng)度可靠性靈敏度分析[J]. 固體火箭技術(shù), 2011, 34(6): 768-771.

ZHOU Hong, LIU Yongshou, GAO Zongzhan, et al. Analysis on strength reliability and parameter sensitivity of a missile cabin structure[J]., 2011, 34(6): 768-771.

[10] 周新建, 萬正平, 周長國. 基于NESSUS的汽車前橋可靠性靈敏度分析[J]. 機(jī)械設(shè)計, 2009, 26(1): 66-68.

ZHOU Xinjian, WAN Zhengping, ZHOU Changguo. Sensitivity analysis on reliability of automobile front axle based on NESSUS[J]., 2009, 26(1): 66-68.

NESSUS-Based Reliability Analysis for Dewar Structure

WANG Chunsheng，ZHANG Lei，DONG Haijie

(,100015,)

To improve the reliability level of a dewar structure, in this paper, a NESSUS(numerical evaluation of stochastic structures under stress)based reliability analysis method for a dewar structure has been proposed. In this method, the design parameters, such as the loads, material characters, and geometry parameters, are treated as random variables with specified distributions. Through probabilistic analysis, this method can obtain the quantitative reliability index of the dewar structure and the sensitivity factors of every design parameter can be calculated. According to the analysis results, the key parameters of dewar structure reliability can be recognized which can be used to improve the dewar design. As an example, the reliability of a cold platform was analyzed in this paper. The results indicate that this method would dramatically improve the reliability of a dewar structure.

structure reliability，NESSUS，reliability simulation，dewar structure

TN215

A

1001-8891(2017)10-0903-05

2016-09-20；

2016-10-20.

王春生（1974-），男，黑龍江延壽人，華北光電技術(shù)研究所高級工程師，博士，研究方向為紅外探測器杜瓦結(jié)構(gòu)設(shè)計、有限元分析及可靠性分析。E-mail：w_chunsheng2003@163.com。