基于單循環(huán)功能度量法的結(jié)構(gòu)可靠壽命優(yōu)化設(shè)計

劉　勤，孫志禮，涂宏茂，姬廣振

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

?

基于單循環(huán)功能度量法的結(jié)構(gòu)可靠壽命優(yōu)化設(shè)計

劉勤1,2，孫志禮1，涂宏茂1,2，姬廣振2

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

摘要：針對常規(guī)的壽命優(yōu)化設(shè)計，未考慮疲勞、磨損影響系數(shù)等因素的分散性，引入單循環(huán)方法，提出了一種以可靠壽命為目標(biāo)、約束的結(jié)構(gòu)優(yōu)化設(shè)計方法；將可靠壽命迭代求解與優(yōu)化計算迭代同步進(jìn)行，顯著提高計算效率；通過對某型車輛傳動箱的結(jié)構(gòu)可靠壽命優(yōu)化設(shè)計結(jié)果表明該方法具有較好的效果和計算效率。

關(guān)鍵詞：優(yōu)化設(shè)計；可靠壽命；單循環(huán)；功能度量法

本文引用格式：劉勤，孫志禮，涂宏茂，等.基于單循環(huán)功能度量法的結(jié)構(gòu)可靠壽命優(yōu)化設(shè)計[J].兵器裝備工程學(xué)報，2016(6):6-9.

Citationformat:LIUQin,SUNZhi-li,TUHong-mao,etal.StructuralReliableLifeOptimizationDesignBasedonSingle-LoopPerformanceMeasureApproach[J].JournalofOrdnanceEquipmentEngineering,2016(6):6-9.

結(jié)構(gòu)載荷、材料性能及疲勞、磨損影響因素等的不確定性會導(dǎo)致機(jī)械結(jié)構(gòu)壽命具有較大的變異性?？煽繅勖紤]了影響因素的隨機(jī)性，定義為結(jié)構(gòu)給定的可靠度所對應(yīng)的壽命單位數(shù)[1]。常規(guī)的結(jié)構(gòu)優(yōu)化方法由于難以系統(tǒng)考慮不確定性因素的影響，優(yōu)化結(jié)果往往并不可靠[2-4]。以可靠壽命為目標(biāo)的優(yōu)化設(shè)計，在滿足性能、結(jié)構(gòu)布局、費用等條件下，力求使機(jī)械結(jié)構(gòu)壽命最長、耐久度最高，從而有效提升產(chǎn)品的平均故障間隔時間等系統(tǒng)可靠性指標(biāo)。

由于可靠壽命優(yōu)化模型的目標(biāo)、約束考慮了影響因素的隨機(jī)性，其可靠壽命、可靠度計算為一個復(fù)雜的迭代計算過程，因此可靠壽命優(yōu)化屬于概率優(yōu)化問題。在優(yōu)化迭代計算過程中，需要采用轉(zhuǎn)換的策略，按照一定的方式將可靠壽命目標(biāo)、可靠度約束轉(zhuǎn)換為確定性目標(biāo)、約束，從而將概率優(yōu)化轉(zhuǎn)換為常規(guī)確定優(yōu)化，再利用常規(guī)的優(yōu)化算法實現(xiàn)求解。對于概率優(yōu)化求解，主要有三類方法：雙循環(huán)方法，單循環(huán)方法和解耦方法[5-8]。雙循環(huán)方法采用兩個嵌套的優(yōu)化循環(huán)：設(shè)計優(yōu)化循環(huán)(外層)和可靠度分析循環(huán)(內(nèi)層)。單循環(huán)方法[9]是在雙循環(huán)的基礎(chǔ)上改進(jìn)的，內(nèi)層循環(huán)由單次可靠度計算近似代替，實現(xiàn)可靠性分析和優(yōu)化計算的同步收斂，優(yōu)化效率提升顯著。

1　結(jié)構(gòu)可靠壽命優(yōu)化模型

工程中常以可靠壽命作為耐久性指標(biāo)，如汽車底盤首次大修期B10≥10 000km，即要求90%的底盤壽命達(dá)到10 000km，才滿足指標(biāo)要求[10]。

若壽命函數(shù)表示為N(x，d)，式中d為設(shè)計變量向量，x為隨機(jī)變量向量；給定可靠度R時，可靠壽命NR滿足

(1)

式中，P{·}為概率函數(shù)。由此，可靠壽命函數(shù)表示為NR(x，d)。對于結(jié)構(gòu)疲勞壽命的預(yù)計，工程常用的以下幾類方法：應(yīng)力壽命法、應(yīng)變壽命法、斷裂力學(xué)法、損傷力學(xué)法。應(yīng)用上述方法，得到壽命函數(shù)N(x，d)，然后根據(jù)上述概率函數(shù)轉(zhuǎn)換為可靠壽命函數(shù)。

根據(jù)可靠壽命為約束還是目標(biāo)，結(jié)構(gòu)可靠壽命優(yōu)化模型一般可以分為2種：可靠壽命約束下結(jié)構(gòu)費用、結(jié)構(gòu)質(zhì)量等極小化或性能最優(yōu)化(簡稱可靠壽命為約束的優(yōu)化模型)；在性能、費用、質(zhì)量、強(qiáng)度、剛度、等約束下可靠壽命極大化或(簡稱可靠壽命為目標(biāo)的優(yōu)化模型)。這兩種類型優(yōu)化實質(zhì)上是一致的，是同一個問題的兩種不同表達(dá)方式。

1) 可靠壽命為目標(biāo)的優(yōu)化模型

在一些新的結(jié)構(gòu)設(shè)計時，往往事先規(guī)定了費用、性能等的界限，可靠壽命要求越大越好。這類問題即為，在結(jié)構(gòu)費用、結(jié)構(gòu)質(zhì)量等特性約束下可靠壽命極大化，其模型如下

(2)

2) 可靠壽命為約束的優(yōu)化模型

根據(jù)用戶需求提出的耐久性指標(biāo)，對結(jié)構(gòu)進(jìn)行設(shè)計，以可靠壽命為約束，性能、費用、質(zhì)量等特性極小化，適合于此類設(shè)計，其模型如下：

(3)

式中，f(d)為目標(biāo)函數(shù)，例如結(jié)構(gòu)費用、質(zhì)量、體積等的函數(shù)；Ng為耐久性指標(biāo)，其他同前。

2　基于單循環(huán)功能度量法的結(jié)構(gòu)可靠壽命優(yōu)化求解

1) 結(jié)構(gòu)可靠壽命優(yōu)化求解流程

通過近似功能度量分析，每次進(jìn)行確定性優(yōu)化前將所有可靠壽命約束近似轉(zhuǎn)換為線性約束。針對結(jié)構(gòu)可靠壽命優(yōu)化模型，利用單循環(huán)方法，如圖3所示。

圖1　基于單循環(huán)功能度量法的結(jié)構(gòu)可靠壽命優(yōu)化求解流程

在序列二次規(guī)劃法(SQP)[11]等優(yōu)化算法的每一次優(yōu)化迭代步，增加在設(shè)計點dk處計算可靠壽命目標(biāo)或約束、可靠性約束值，并將其轉(zhuǎn)換成近似函數(shù)，實現(xiàn)優(yōu)化求解與可靠壽命、可靠度求解的同步收斂，從而得到最優(yōu)解和最優(yōu)值(2)可靠壽命目標(biāo)/約束的轉(zhuǎn)換。

在優(yōu)化迭代過程，利用功能度量法[12-13]，構(gòu)造并求解近似規(guī)劃問題前，將可靠壽命目標(biāo)/約束在當(dāng)前設(shè)計點處進(jìn)行泰勒展開，轉(zhuǎn)換成線性約束：

(4)

上式采用功能度量法的迭代公式進(jìn)行一次近似可靠壽命分析，k表示優(yōu)化迭代次數(shù)。若給定可靠度為R=φ(β)，β為可靠度指數(shù)，其中近似可能失效點根據(jù)以下2式計算

式中，▽uNu(dk-1,uk-1)是壽命函數(shù)關(guān)于隨機(jī)向量u的梯度向量。

當(dāng)設(shè)計變量是隨機(jī)變量的統(tǒng)計參數(shù)時，設(shè)計變量梯度向量的對應(yīng)元素由下式計算：

(7)

式中，▽uNu(dk-1,uk-1)可通過功能度量法的迭代公式求解，▽yu為標(biāo)準(zhǔn)正態(tài)變量對原設(shè)計變量概率分布參數(shù)的偏導(dǎo)數(shù)。

當(dāng)設(shè)計變量是確定性參量時，設(shè)計變量梯度向量的對應(yīng)元素由下式計算：

(8)

式中，hi為差分步長，ei為除第i個分量為1，其余分量為零的單位向量。

2) 可靠度約束的轉(zhuǎn)換

利用可靠性優(yōu)化方法對可靠性約束的轉(zhuǎn)換，如對于第i個可靠度約束，令P{gi(d,x)≥0} =Φ(βi)，βi為可靠度系數(shù)；在優(yōu)化迭代過程，該約束在當(dāng)前設(shè)計迭代點轉(zhuǎn)化為線性約束函數(shù)，即：

(9)

(10)

式中，Gi(dk-1，uk-1)是功能函數(shù)gi(dk-1，xk-1)由隨機(jī)向量x轉(zhuǎn)換至獨立標(biāo)準(zhǔn)正態(tài)空間向量u后的函數(shù)；k表示優(yōu)化迭代步序號；▽uGi(dk-1，uk-1)是Gi(dk-1，uk-1)對u的梯度向量。下一步迭代的uk由式(11)計算。

(11)

3) 優(yōu)化求解

可靠壽命為目標(biāo)的優(yōu)化模型轉(zhuǎn)換為：

(12)

可靠壽命為約束的優(yōu)化模型轉(zhuǎn)換為：

(13)

轉(zhuǎn)換為常規(guī)的優(yōu)化模型后，利用序列二次規(guī)劃(SQP)等優(yōu)化算法即可求解。

3　應(yīng)用實例

某車輛傳動系統(tǒng)的耐久性指標(biāo)B10為10 000 km，對該箱體進(jìn)行結(jié)構(gòu)優(yōu)化，問題可以描述為：在滿足耐久性指標(biāo)的情況下，通過優(yōu)化箱體結(jié)構(gòu)，減輕箱體質(zhì)量，其一般模型表示為：

式中，Weight( )表示箱體的質(zhì)量函數(shù)，LifeR=0.9()表示箱體可靠壽命函數(shù)。

表1　箱體結(jié)構(gòu)設(shè)計參數(shù)

經(jīng)過有限元計算，獲得箱體最大應(yīng)力為132.3 MPa，在輸入軸區(qū)域，如圖2所示。為確定滿足耐久性指標(biāo)要求時允許的箱體初始裂紋尺寸最大值a0，故選取a0作為箱體設(shè)計參數(shù)，按保守設(shè)計，裂紋位置位于最大應(yīng)力附近。

圖2　箱體局部應(yīng)力分布

利用Promrel軟件，通過拉丁超立方設(shè)計生成100組試驗樣本，將箱體有限元分析參數(shù)化，編制APDL，并利用Paris公式計算裂紋擴(kuò)展壽命。改變樣本參數(shù)值，分別進(jìn)行100次裂紋擴(kuò)展壽命計算，擬合了箱體壽命的二次響應(yīng)面模型，如下：

Life=262.4+1.3×10-2×B0-1.32×B1-

0.83×B2-0.73×H1-0.25×B3-0.15×B4-

1.3×10-5×E-2.2×10-4×B1×B2-

2.5×10-5×B1×B3+2.7×10-4×B0×H1-

8.3×10-4×B0×B3+0.011×B3×B4+

其中，E為材料彈性模量。箱體質(zhì)量W的響應(yīng)面函數(shù)為：

W=4.73×B0+0.115×B1+0.013×B2+

0.039×H1+0.175×B3+0.026×B4+219.4

為驗證本文的方法，分別利用單循環(huán)功能度量法、雙循環(huán)功能度量法進(jìn)行了可靠壽命為目標(biāo)、約束的優(yōu)化求解，其結(jié)果列入表2。

表2　箱體結(jié)構(gòu)耐久性優(yōu)化結(jié)果對比

從表2中可以看出，可靠壽命為目標(biāo)的結(jié)構(gòu)優(yōu)化獲得了更長壽命的方案，可靠壽命為約束的結(jié)構(gòu)優(yōu)化獲得滿足耐久性約束下質(zhì)量最輕的方案；對于可靠壽命優(yōu)化，單循環(huán)功能度量法與雙循環(huán)功能度量法優(yōu)化結(jié)果相近，但單循環(huán)功能度量法的計算效率較雙循環(huán)高14%；以常規(guī)壽命為約束的結(jié)構(gòu)優(yōu)化方案，雖然質(zhì)量更輕，但經(jīng)過驗算后，發(fā)現(xiàn)該點的可靠壽命為8 600 km,不滿足要求。

4　結(jié)論

本文將單循環(huán)功能度量法引入結(jié)構(gòu)壽命優(yōu)化設(shè)計，提出了一種以可靠壽命為目標(biāo)、約束的結(jié)構(gòu)優(yōu)化設(shè)計方法；可靠壽命為目標(biāo)的結(jié)構(gòu)優(yōu)化可獲得更長壽命的方案，可靠壽命為約束的結(jié)構(gòu)優(yōu)化可獲得滿足耐久性約束下質(zhì)量最輕的方案。該方法將可靠壽命的迭代求解，與優(yōu)化計算的迭代同步進(jìn)行，提高計算效率，這對于復(fù)雜結(jié)構(gòu)優(yōu)化設(shè)計具有重要的工程意義。

通過對某型車輛傳動箱的結(jié)構(gòu)可靠壽命優(yōu)化設(shè)計，并與雙循環(huán)功能度量法、常規(guī)壽命優(yōu)化設(shè)計等計算結(jié)果相比較，表明基于單循環(huán)功能度量法的結(jié)構(gòu)可靠壽命優(yōu)化設(shè)計方法具有較好的效果和計算效率。相對于雙循環(huán)功能度量法優(yōu)化方法，單循環(huán)功能度量法的結(jié)果相近，但計算效率較雙循環(huán)高14%；以常規(guī)壽命為約束的結(jié)構(gòu)優(yōu)化方案，雖然質(zhì)量更輕，但經(jīng)過驗算后，發(fā)現(xiàn)該點的可靠壽命不滿足要求。

參考文獻(xiàn)：

[1]張銀龍，申兆祥，卞士川，等．裝備可靠性、耐久性與壽命之間的關(guān)系[J]．四川兵工學(xué)報，2013(8)：76-79.

[2]RACKWITZ R.Reliability analysis-a review and some perspectives[J].Structural Safety，2001，23(4):365-395.

[3]GUAN X F，HE J J，RASSELKORDE El M.Probabilistic Fatigue Life Prediction and Structural Reliability Evaluation of Turbine Rotors Integrating an Automated Ultrasonic Inspection System[J].J.Nondestruct Eval，2014，33:51-61.

[4]謝里陽，任俊剛，吳寧祥,等.復(fù)雜結(jié)構(gòu)部件概率疲勞壽命預(yù)測方法與模型[J].航空學(xué)報，2015，36(8):2688-2695.

[5]RAMUP，QUX，YOUNB，etal.Inversereliabilitymeasuresandreliability-baseddesignoptimization[J].IntJReliabilityandSafety，2006，1(1/2):187-205.

[6]LIG，MENGZ，HUH.Anadaptivehybridapproachforreliability-baseddesignoptimization[J].StructMultidiscOptim(2015) 51:1051-1065.

[7]YOUNESA，ALAAC.Benchmarkstudyofnumericalmethodsforreliability-baseddesignoptimization[J].StructuralandMultidisciplinaryOptimization，2010，41(2):277-294.

[8]YIP，CHENGGD，JiangL.Asequentialapproximateprogrammingstrategyforperformance-measure-basedprobabilisticstructuraldesignoptimization[J].StructuralSafety，2008,30(2):91-109.

[9]LIMJ，LEEB.Asemi-single-loopmethodusingapproximationofmostprobablepointforreliability-baseddesignoptimization[J].StructMultidiscOptim，2015(11):1-13.

[10]劉勤，錢云鵬，姬廣振，等．隨機(jī)載荷作用下扭力軸耐久性優(yōu)化設(shè)計[J]．兵工學(xué)報，2015(5)：933-937.

[11]袁亞湘，孫文瑜．最優(yōu)化理論與方法[M]．北京：科學(xué)出版社，2005.

[12]RAHMANS，WEID.Designsensitivityandreliability-basedstructuraloptimizationbyunivariatedecomposition[J].StructMultidiscOptim，2008,35(3):245-261.

[13]LIL，LIUJH,LIUSH.AnefficientstrategyformultidisciplinaryreliabilitydesignandoptimizationbasedonCSSOandPMAinSORAframework[J].StructMultidiscOptim，2014(49):239-252.

[14]趙克剛，李健，姚偉浩．車輛傳動系統(tǒng)固有特性的識別[J]．重慶理工大學(xué)學(xué)報(自然科學(xué))，2014(4)：5-9．

(責(zé)任編輯周江川)

doi:【裝備理論與裝備技術(shù)】10.11809/scbgxb2016.06.002

收稿日期：2016-01-17；修回日期：2016-02-20

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

作者簡介：劉勤(1981—)，博士，副研究員，主要從事機(jī)械可靠性設(shè)計分析研究。

中圖分類號：TH122

文獻(xiàn)標(biāo)識碼：A

文章編號：2096-2304(2016)06-0006-05

StructuralReliableLifeOptimizationDesignBasedonSingle-LoopPerformanceMeasureApproach

LIUQin1,2,SUNZhi-li1,TUHong-mao1,2,JIGuang-zhen2

(1.SchoolofMechanicalEngineeringandAutomation,NortheasternUniversity,Shenyang110819,China; 2.OrdnanceScienceandResearchAcademyofChina,Beijing100089,China)

Abstract:The traditional structure life optimization doesn’t consider randomness of factors for fatigue, wear, etc. A structure optimization approach, which set reliabe life as object and took constraint structure optimization as the design method, was proposed by using single-loop method. The new sequential approximate programming realizes concurrent convergence of both design optimization and reliabe life calculation, and the present algorithm is more efficient. Finally, the optimization design of transmission case in vehicle was employed. The result demonstrates the good effectiveness of the presented method.

Key words:optimization design; reliable life; single-loop; performance measure approach