寬帶隨機振動試驗條件的加速因子

朱學(xué)旺，張思箭，寧佐貴，劉青林

（中國工程物理研究院總體工程研究所，綿陽 621999）

概述

振動環(huán)境工程研究中，經(jīng)常采用加速方法來進行振動環(huán)境試驗。此時，需要科學(xué)設(shè)計試驗條件的加速因子。GJB 150.16A[1]和美軍標(biāo)MIL-STD-810F[2]推薦的包裝運輸試驗條件和導(dǎo)彈飛行器的振動試驗條件都是加速試驗條件，其中基本包裝運輸試驗條件的時間壓縮比達(dá)到1600km/h，導(dǎo)彈飛行器振動試驗的時間為1h。在這兩個標(biāo)準(zhǔn)中，加速因子的設(shè)計，應(yīng)用了基于疲勞損傷等效的加速公式：

式中， G0、G1；T0、T1分別為加速前后的功率譜密度和試驗時間。

這是依據(jù)線性疲勞損傷累積模型獲得的結(jié)果，其中的冪指數(shù)（1/4）可以依據(jù)材料的疲勞試驗的結(jié)果進行調(diào)整。

韓雪山等利用公式（1）討論了包裝產(chǎn)品的運輸試驗的加速研究[3]；李奇志[4]等則通過疲勞損傷的WL模型[5]分析，認(rèn)為公式（1）不僅對于窄帶隨機的加速是合理的，對于寬帶隨機同樣適用；文獻(xiàn)[6]研究認(rèn)為，公式（1）僅對窄帶隨機是嚴(yán)格成立的，對于寬帶隨機需要對其勞損傷進行專門的分析，才能獲得有效的疲勞損傷等效結(jié)果。

寬帶隨機振動的疲勞損傷依然是一個開放的研究領(lǐng)域。一種觀點認(rèn)為，窄帶模型疲勞損傷分析結(jié)果具有廣泛的適應(yīng)性，可通過對其結(jié)果的適當(dāng)修正以滿足寬帶隨機振動情形[7]，典型修正方法有WL方法、 α0.75方法[8]、TB方法[9][10]等。另一種觀點認(rèn)為，寬帶隨機振動下的疲勞損傷需要建立新的模型，例如Drilik模型[8][11]。

本文研究公式（1）表述的加速因子對于寬帶隨機振動適宜性。寬帶隨機振動的疲勞損傷估計采用窄帶模型修正方法來描述，試驗條件的加速因子分析分別考慮比例載荷與非比例載荷兩種情況。為了描述的方便性，首先給出了窄帶隨機振動試驗加速因子的理論推導(dǎo)過程。

1 窄帶隨機振動試驗條件的加速因子

按照線性疲勞累積損傷Palmgren-Miner假設(shè)，滿足Rayleigh分布的窄帶隨機振動疲勞損傷累積可表述為[4][8]:

容易導(dǎo)出，兩個疲勞損傷累積等效的窄帶隨機過程存在下列關(guān)系式：

當(dāng)這兩個窄帶隨機過程為比例載荷時，即在任何頻率f 處，其功率譜密度（PSD）都有關(guān)系：

式中，k為常數(shù)時，公式（4）可改寫為：

公式（4）與公式（6）為窄帶隨機振動試驗條件的加速因子的兩種表達(dá)形式。實際工程應(yīng)用時要注意兩點：一是公式（6）僅適用于比例載荷情形；二是公式（1）暗示了選用了材料的疲勞參數(shù)b=8，且為比例載荷。

2 寬帶隨機振動試驗條件的加速因子

為了分析公式（4）和公式（6）表達(dá)的加速因子對于寬帶隨機振動是否適用，首先給出寬帶隨機振動的疲勞損傷累積研究中的幾種窄帶模型修正方法，包括：1980年代Wirsching和Light提出的譜寬參數(shù)（也稱譜型不規(guī)則因子）修正方法[5]（WL方法）、2002年和2004年由Benasciutti 和 Tovo提出的 α0.75方法[8]與線性組合方法[9][10]（TB方法）。它們的疲勞損傷累積計算公式分別為（7）、（8）和（9）。

公式（9）中，β可以按照下列關(guān)系式來取值β1或β2：

將寬帶隨機振動疲勞損傷累積的不同描述改寫為統(tǒng)一的格式：

按照寬帶隨機振動疲勞損傷累積等效，有關(guān)系式：

公式（12）為寬帶隨機振動試驗條件加速因子的通用表達(dá)式。當(dāng)公式（13）簡化為公式（4）。容易檢驗，對于公式（5）描述的比例載荷，成立，且可以導(dǎo)出公式（6），此時加速試驗時間可以根據(jù)公式（4）或公式（6）獲得。當(dāng)加速載荷不是比例載荷時，不成立，此時的加速試驗時間需要根據(jù)加速載荷的分布通過公式（13）求得。

3 數(shù)值算例及討論

以圖1所示的寬帶隨機振動及其加速試驗條件為例，計算加速試驗時間。圖中條件a為未加速的原始條件PSD（對應(yīng)的試驗時間1000秒），其對應(yīng)的總均方根值為1.1674g；條件b、條件c和條件d分別為總均方根值加速比均為1.2（即傳統(tǒng)意義上的加速因子為1.4）的三種加速條件PSD，其中條件b為比例放大，條件c和條件d為非比例放大。表1為四種條件的具體參數(shù)值。

表2為根據(jù)公式（13）獲得的加速試驗時間（材料疲勞參數(shù)分別取4，6，8）。為比較方便，表中還給出了窄帶方法分析結(jié)果。窄帶模型修正方法計算時，公式（2）中的v0均用替代。

趨勢一，對于比例載荷，不同的方法獲得的加速試驗時間相等，即對于比例載荷，公式（4）或（6）提供的加速因子對于寬帶隨機振動的加速試驗設(shè)計是有效的；

趨勢二，對于非比例載荷的寬帶隨機振動加速，從公式（4）無法導(dǎo)出公式（6），其加速試驗時間會隨加速條件譜型的不同而改變，相同的能量加速因子下，其加速試驗時間可以大于比例載荷的加速試驗時間，也可以小于比例載荷的加速試驗時間；

趨勢三，能量加速因子相同時，如果低頻不加速而高頻加速，其加速試驗時間明顯小于按比例載荷獲得的加速試驗時間，反之，如果高頻不加速而僅在低頻加速，其加速試驗時間則比按比例載荷求得的加速試驗時間要長。這是符合物理規(guī)律的，因為相同的時間內(nèi)，高頻載荷的總作用次數(shù)多，造成的疲勞損傷累積相對較大，而低頻載荷的總作用次數(shù)少，造成的疲勞損傷累積相對較小。

圖1 四種寬帶隨機振動試驗條件的PSD 曲線

趨勢四，如果不考慮載荷的譜型，直接按照公式（4）來設(shè)計加速試驗時間，則當(dāng)PSD類似條件c（以高頻PSD加速為主）時，會產(chǎn)生過試驗；而當(dāng)PSD類似條件d（以低頻PSD加速為主）時，會產(chǎn)生欠試驗。

4 結(jié)論

文中導(dǎo)出了寬帶隨機振動試驗加速因子的通用表達(dá)式，并設(shè)計了加速因子相同的三種加速試驗條件PSD，應(yīng)用寬帶隨機振動的疲勞損傷累積研究中的幾種窄帶模型修正方法，分別計算得出了相應(yīng)的加速試驗時間。結(jié)果顯示，對于比例載荷的加速試驗，傳統(tǒng)的方法（公式（4）和公式（6））適用于寬帶隨機振動的加速試驗設(shè)計；而對于非比例型載荷（文中僅以分段比例載荷為例），按照傳統(tǒng)方法設(shè)計的加速試驗會因為加速試驗條件PSD的不同而出現(xiàn)過試驗或欠試驗，當(dāng)加速試驗條件PSD以高頻分量加速為主時，多為過試驗，當(dāng)加速試驗條件PSD以低頻分量加速為主時，多為欠試驗。

文中加速因子的定義采用了公式（4）和（13）的形式而沒有采用公式（1）和公式（6），這是因為討論非比例載荷加速試驗條件的需要。實際上，對于非比例載荷的加速試驗而言，公式（1）與公式（6）所定義的加速因子并非常數(shù)而是與頻率相關(guān)的變量。

表1 四種寬帶隨機振動試驗條件的PSD 值

表2 加速試驗時間（材料疲勞參數(shù)取b=4,6,8）

文中的討論僅局限于疲勞損傷的一類頻域描述方法——窄帶模型修正方法，當(dāng)采用其它方法描述疲勞損傷時的加速因子分析可類似處理，盡管具體結(jié)果的數(shù)值可能有差異，但是趨勢會相同：相對于比例載荷加速試驗時間而言，以高頻分量加速為主的加速試驗時間更短，以低頻分量加速為主時則加速試驗時間更長。

[1]GJB 150.16A－2009,軍用裝備實驗室環(huán)境試驗方法第16 部分振動試驗[S].

[2]MIL-STD-810F, Environmental engineering considerations and laboratory tests, METHOD 514.5VIBRATION[S].

[3]韓雪山,牛淑梅,何淵井.運輸包裝隨機振動加速試驗探討[J].中國包裝, 2011, 11: 51 ～54.

[4]李奇志,陳國平,王明旭,等.振動加速因子試驗方法研究[J].振動、測試與診斷, 2013, 33(1): 35～39 .

[5]Wirsching PH, Light CL. Fatigue under wide band random stresses[J]. ASCE, J Struct Div 1980; 106(7):1593–607.

[6]王冬梅,謝勁松.隨機振動試驗加速因子的計算方法[J].環(huán)境技術(shù),2010, 28(2): 47-51.

[7]朱學(xué)旺,張思箭,寧佐貴等. 3σ 準(zhǔn)則的應(yīng)用對疲勞損傷累積分析的影響[C].中國宇航學(xué)會結(jié)構(gòu)強度與環(huán)境工程專委會、中國航天第十專業(yè)信息網(wǎng)、五院科技委環(huán)境工程與可靠性專業(yè)組2014年度學(xué)術(shù)交流會論文集, 航天器環(huán)境工程, 2014 : 571-576.

[8]Matjaz Mr?nik, Janko Slavic, Miha Boltezar, Frequency-domain methods for a vibration-fatigue-life estimation Application to real data[J]. International Journal of Fatigue, 2013, 47: 8-17.

[9]Tovo R. Cycle distribution and fatigue damage under broadband random loading[J]. Int J Fatigue 2002, 24(11): 1137–47.

[10]Benasciutti D, Tovo R. Spectral methods for lifetime prediction under wideband stationary random processes[J]. Int J Fatigue,2005, 27(8): 867–77.

[11]Dirlik T. Application of computers in fatigue analysis[D]. Ph.D.hesis. The University of Warwick, 1985.

[12]張翼,楊晨,羅楊陽.隨機振動載荷下導(dǎo)彈吊掛疲勞壽命分析[J].機械科學(xué)與技術(shù), 2013, 32(11): 1675-1679.