超高斯隨機(jī)振動(dòng)疲勞加速試驗(yàn)?zāi)Ｐ脱芯?/h1>
2017-05-17 02:35:10陶俊勇
振動(dòng)與沖擊 2017年9期
蔣 瑜， 陶俊勇， 陳 循
(國(guó)防科技大學(xué) 機(jī)電工程與自動(dòng)化學(xué)院裝備綜合保障技術(shù)重點(diǎn)實(shí)驗(yàn)室， 長(zhǎng)沙 410073)
超高斯隨機(jī)振動(dòng)疲勞加速試驗(yàn)?zāi)Ｐ脱芯?/p>
蔣 瑜， 陶俊勇， 陳 循
(國(guó)防科技大學(xué) 機(jī)電工程與自動(dòng)化學(xué)院裝備綜合保障技術(shù)重點(diǎn)實(shí)驗(yàn)室， 長(zhǎng)沙 410073)
通過(guò)理論分析，針對(duì)超高斯隨機(jī)振動(dòng)激勵(lì)建立了相應(yīng)的振動(dòng)疲勞加速試驗(yàn)數(shù)學(xué)模型，給出了模型中未知參數(shù)的具體求解方法，并通過(guò)實(shí)際試驗(yàn)進(jìn)行了驗(yàn)證。該模型十分便于進(jìn)行振動(dòng)疲勞加速試驗(yàn)的定量設(shè)計(jì)，具有很好的工程實(shí)用性，可用于評(píng)估隨機(jī)振動(dòng)環(huán)境下工程結(jié)構(gòu)的疲勞壽命與可靠性。
超高斯； 隨機(jī)振動(dòng)； 振動(dòng)疲勞； 加速模型
振動(dòng)引起的疲勞問(wèn)題作為多個(gè)工程領(lǐng)域廣泛存在的一個(gè)共性問(wèn)題，嚴(yán)重危及重大裝備及結(jié)構(gòu)的可靠性和安全性。如果能夠提前準(zhǔn)確預(yù)測(cè)結(jié)構(gòu)的振動(dòng)疲勞壽命，就能在發(fā)生災(zāi)難性事故之前及時(shí)預(yù)知并采取相應(yīng)的措施，并可為其定壽、延壽提供科學(xué)的依據(jù)，最大限度地發(fā)揮裝備效益。因此，準(zhǔn)確預(yù)測(cè)重大裝備和工程結(jié)構(gòu)在復(fù)雜隨機(jī)動(dòng)態(tài)載荷作用下的振動(dòng)疲勞壽命是提高其可靠性和安全性的關(guān)鍵技術(shù)，被國(guó)家自然科學(xué)基金委員會(huì)列入《機(jī)械工程學(xué)科發(fā)展戰(zhàn)略報(bào)告(2011～2020)》。
當(dāng)前工程中進(jìn)行結(jié)構(gòu)振動(dòng)疲勞壽命分析和預(yù)測(cè)時(shí)，由于缺乏充分的認(rèn)識(shí)和解決相關(guān)問(wèn)題的工具及方法，通常假設(shè)結(jié)構(gòu)承受的隨機(jī)載荷符合平穩(wěn)高斯分布。但是，在實(shí)際環(huán)境中，上述許多隨機(jī)載荷往往呈現(xiàn)比較明顯的非高斯特征，尤其是在工況惡劣或者極端環(huán)境下。近年來(lái)的一些研究結(jié)果表明，隨機(jī)載荷的非高斯特性對(duì)結(jié)構(gòu)振動(dòng)疲勞有著不可忽視的影響，某些情況下會(huì)加速結(jié)構(gòu)的疲勞失效?；诟咚辜僭O(shè)計(jì)算振動(dòng)疲勞損傷往往會(huì)得到偏大的振動(dòng)疲勞壽命估計(jì)結(jié)果，給裝備服役或使用階段埋下巨大的安全隱患。
同時(shí)，隨著結(jié)構(gòu)可靠性水平的提高，結(jié)構(gòu)的振動(dòng)疲勞壽命越來(lái)越長(zhǎng)，為了能夠在實(shí)驗(yàn)室驗(yàn)證其壽命是否達(dá)到要求，加速試驗(yàn)成為必然的選擇。要想通過(guò)加速試驗(yàn)得到結(jié)構(gòu)在實(shí)際服役振動(dòng)環(huán)境下的疲勞壽命，正確的加速試驗(yàn)統(tǒng)計(jì)模型是關(guān)鍵。國(guó)內(nèi)外目前關(guān)于振動(dòng)加速試驗(yàn)?zāi)Ｐ偷难芯坎欢啵绕涫菍?duì)非高斯隨機(jī)振動(dòng)加速試驗(yàn)。王冬梅等[1]對(duì)振動(dòng)加速試驗(yàn)的逆冪律模型進(jìn)行了推導(dǎo)，探討了其適用范圍，指出其適用于窄帶高斯載荷，不適用于寬帶高斯和非高斯載荷。李奇志等[2]提出通過(guò)試驗(yàn)的方式獲得振動(dòng)試驗(yàn)的加速因子，認(rèn)為振動(dòng)加速試驗(yàn)的逆冪律模型對(duì)平穩(wěn)窄帶和寬帶高斯隨機(jī)過(guò)程均是適用的。朱學(xué)旺等[3]應(yīng)用基于窄帶模型的修正方法得到了寬帶隨機(jī)振動(dòng)試驗(yàn)加速因子計(jì)算的通用表達(dá)式，認(rèn)為基于窄帶模型的加速因子表達(dá)式對(duì)于比例載荷的寬帶隨機(jī)振動(dòng)也是適用的，而對(duì)于非比例載荷，則需要應(yīng)用其提出的通用表達(dá)式才可以獲得。Allegri等[4]研究了適用于平穩(wěn)寬帶高斯隨機(jī)振動(dòng)加速試驗(yàn)的逆冪律模型，其主要目的是用來(lái)評(píng)估相對(duì)疲勞損傷，而非用來(lái)精確預(yù)計(jì)壽命。Pothula[5]研究了附著不同阻尼材料的鋁合金梁在高斯隨機(jī)振動(dòng)加速試驗(yàn)中的疲勞壽命差異，探討了結(jié)構(gòu)阻尼對(duì)振動(dòng)試驗(yàn)加速因子的影響。John[6]提出應(yīng)用Fatigue Damage Spectrum(FDS)來(lái)進(jìn)行振動(dòng)加速試驗(yàn)設(shè)計(jì)與評(píng)估，但是只考慮了振動(dòng)激勵(lì)的功率譜密度PSD (Power Spectral Density)，只適用于高斯振動(dòng)加速試驗(yàn)。
總的來(lái)說(shuō)，目前國(guó)內(nèi)外開(kāi)展的振動(dòng)疲勞加速試驗(yàn)研究主要針對(duì)高斯振動(dòng)，尤其是實(shí)際進(jìn)行的非高斯振動(dòng)加速試驗(yàn)很少，需要進(jìn)一步探索適用于非高斯隨機(jī)振動(dòng)的疲勞加速試驗(yàn)定量模型。
1 理論分析
工程上常用偏斜度S和峭度K這兩個(gè)參數(shù)來(lái)描述非高斯隨機(jī)過(guò)程X，定義為
(1)
高斯隨機(jī)過(guò)程的偏斜度值等于0，峭度值等于3；而非高斯隨機(jī)過(guò)程的峭度值肯定不等于3，偏斜度值可以等于零也可以不等于0。偏斜度用來(lái)描述隨機(jī)過(guò)程幅值概率密度曲線偏離對(duì)稱分布的程度，偏斜度值不為零表示服從非對(duì)稱分布。峭度是描述隨機(jī)過(guò)程幅值概率密度曲線拖尾分布特征的參數(shù)，它不僅可用來(lái)區(qū)分高斯和非高斯隨機(jī)過(guò)程，而且還可進(jìn)一步將非高斯隨機(jī)過(guò)程區(qū)分為亞高斯和超高斯隨機(jī)過(guò)程，其中亞高斯隨機(jī)過(guò)程的K<3，超高斯隨機(jī)過(guò)程的K>3。工程中常見(jiàn)的非高斯振動(dòng)信號(hào)往往是具有尖峰分布的對(duì)稱超高斯信號(hào)，因此本文研究的非高斯振動(dòng)疲勞主要針對(duì)超高斯。
1.1 高斯振動(dòng)疲勞加速試驗(yàn)?zāi)Ｐ?/p>
首先，從最基本的描述材料疲勞現(xiàn)象的S～N曲線出發(fā)。通常理想的S～N曲線的數(shù)學(xué)表達(dá)式為
N=cS-b
(2)
式中：S為應(yīng)力幅值；N為引起失效的應(yīng)力循環(huán)次數(shù)；b、c為材料雙對(duì)數(shù)S～N曲線中的常數(shù)，且b的取值范圍一般為6～25。
根據(jù)著名的疲勞損傷累積Miner準(zhǔn)則，不同幅值應(yīng)力共同作用下造成的疲勞損傷為
(3)
式中：ni是幅值為Si的應(yīng)力作用循環(huán)次數(shù)；Ni為試件在幅值為Si的應(yīng)力作用下至疲勞失效的循環(huán)次數(shù)；D為累積疲勞損傷(一般認(rèn)為當(dāng)D=1時(shí)發(fā)生疲勞失效)。
將式(2)代入式(3)可得
(4)
對(duì)于連續(xù)的隨機(jī)應(yīng)力時(shí)間歷程，式(4)可寫成如下積分的形式[7]
(5)
當(dāng)隨機(jī)應(yīng)力響應(yīng)接近平穩(wěn)窄帶高斯分布時(shí)，根據(jù)隨機(jī)過(guò)程理論，應(yīng)力幅值概率密度函數(shù)p(S)將近似服從瑞利(Rayleigh)分布
(6)
式中：σS為應(yīng)力的均方根值(即標(biāo)準(zhǔn)差)。將式(6)代入式(5)進(jìn)行積分可得
(7)
式中，Γ為Gamma函數(shù)。
進(jìn)行振動(dòng)試驗(yàn)的產(chǎn)品或試件一般可近似看成一個(gè)線性系統(tǒng)，而振動(dòng)試驗(yàn)設(shè)備產(chǎn)生的振動(dòng)激勵(lì)可看作該系統(tǒng)的輸入。工程實(shí)踐表明，一般結(jié)構(gòu)件的阻尼比ξ通常遠(yuǎn)小于1，可以取到0.05以下這樣小的值。根據(jù)文獻(xiàn)[8]，在滿足小阻尼線性系統(tǒng)的假設(shè)前提下，σS的近似計(jì)算式為
(8)
式中：f1為試件結(jié)構(gòu)的一階固有頻率；ξ為等效的阻尼比(一般假設(shè)ξ≤0.1)；k為與試件材料相關(guān)的比例常數(shù)；Ga(f1)為輸入振動(dòng)激勵(lì)的加速度功率譜密度在試件固有頻率f1處的量值。
(9)
由于試件結(jié)構(gòu)的傳遞函數(shù)一般都類似一個(gè)窄帶濾波器，因此可以認(rèn)為試件在平穩(wěn)高斯隨機(jī)激勵(lì)作用下(無(wú)論是寬帶還是窄帶)，其應(yīng)力響應(yīng)為平穩(wěn)窄帶高斯分布，可以采用式(9)來(lái)估計(jì)其疲勞累積損傷。
一般認(rèn)為D=1時(shí)結(jié)構(gòu)發(fā)生疲勞失效，根據(jù)式(9)可得高斯隨機(jī)振動(dòng)激勵(lì)下的結(jié)構(gòu)振動(dòng)疲勞壽命TG為
(10)
對(duì)常用工程材料，b的取值范圍一般為4～25。從式(10)可以看出，對(duì)高斯隨機(jī)振動(dòng)激勵(lì)，當(dāng)結(jié)構(gòu)動(dòng)力學(xué)特性參數(shù)如f1、ξ固定時(shí)，結(jié)構(gòu)振動(dòng)疲勞壽命T與高斯隨機(jī)振動(dòng)激勵(lì)的功率譜密度在結(jié)構(gòu)一階固有頻率處的量值大小Ga(f1)成反比關(guān)系。隨著Ga(f1)的增加(或減小)，T按指數(shù)規(guī)律迅速衰減(或增長(zhǎng))。
根據(jù)式(10)進(jìn)一步可得a1、a2兩種不同高斯隨機(jī)加速度激勵(lì)下的結(jié)構(gòu)振動(dòng)疲勞壽命分別為
(11)
(12)
根據(jù)式(11)和式(12)可得
(13)
顯然，式(13)就是傳統(tǒng)文獻(xiàn)中用于描述高斯振動(dòng)加速試驗(yàn)的逆冪律模型。從上述推導(dǎo)過(guò)程來(lái)看，只要工程結(jié)構(gòu)系統(tǒng)的頻響函數(shù)是一個(gè)窄帶濾波器(實(shí)際中大多數(shù)工程結(jié)構(gòu)均滿足這一條件)，并且結(jié)構(gòu)隨機(jī)響應(yīng)服從高斯分布，則就可以運(yùn)用式(13)對(duì)振動(dòng)加速試驗(yàn)進(jìn)行建模。因此，式(13)只適用于描述高斯振動(dòng)疲勞加速試驗(yàn)，下面繼續(xù)討論非高斯振動(dòng)疲勞加速試驗(yàn)?zāi)Ｐ汀?/p>
1.2 超高斯振動(dòng)疲勞加速試驗(yàn)?zāi)Ｐ?/p>
當(dāng)隨機(jī)應(yīng)力響應(yīng)接近平穩(wěn)窄帶非高斯分布時(shí)，可考慮在式(9)的基礎(chǔ)上增加一個(gè)非高斯修正因子λ來(lái)描述應(yīng)力響應(yīng)的峭度值對(duì)振動(dòng)疲勞累積損傷的影響
(14)
非高斯修正因子λ和應(yīng)力響應(yīng)的峭度值Ks直接相關(guān)，可用式(15)來(lái)描述
(15)
式中，參數(shù)α為比例系數(shù)。
顯然，當(dāng)應(yīng)力響應(yīng)為高斯即Ks=3時(shí)，λ=1，式(14)變成式(9)；當(dāng)應(yīng)力響應(yīng)為超高斯即Ks>3時(shí)，λ>1，意味著應(yīng)力響應(yīng)的超高斯特性會(huì)加快疲勞累積損傷進(jìn)程。
下面根據(jù)隨機(jī)振動(dòng)理論分析影響應(yīng)力響應(yīng)峭度值Ks的因素。由于結(jié)構(gòu)的一階模態(tài)對(duì)結(jié)構(gòu)響應(yīng)起決定性作用，因此建立振動(dòng)臺(tái)基礎(chǔ)激勵(lì)作用下的單自由度系統(tǒng)模型進(jìn)行分析，如下圖1所示。
經(jīng)過(guò)推導(dǎo)可得加速度響應(yīng)y與振動(dòng)臺(tái)臺(tái)面的基礎(chǔ)加速度激勵(lì)x之間的傳遞函數(shù)為
(16)
(17)
(18)
式中：ω1=2πf1，f1為試件結(jié)構(gòu)的一階固有頻率；ξ為阻尼比。這兩個(gè)參數(shù)表征了試件結(jié)構(gòu)本身的動(dòng)力學(xué)特性，并且試件結(jié)構(gòu)系統(tǒng)的通頻帶BWH也由這兩個(gè)參數(shù)決定，即
BWH=2ξf1
(19)
圖1 振動(dòng)臺(tái)基礎(chǔ)激勵(lì)作用下的單自由度模型
由于實(shí)際工程結(jié)構(gòu)的阻尼比ξ往往較小，而一階固有頻率f1往往不大，所以大多數(shù)工程結(jié)構(gòu)的通頻帶BWH往往也不大，也就是說(shuō)可以視為一個(gè)窄帶濾波器。
下面先對(duì)響應(yīng)帶寬進(jìn)行分析。設(shè)激勵(lì)x的功率譜密度記作X(f)，響應(yīng)y的功率譜密度記作Y(f)，系統(tǒng)頻響函數(shù)為H(f)。根據(jù)線性系統(tǒng)和隨機(jī)過(guò)程理論有
(20)
由此容易得出響應(yīng)的有效頻譜帶寬
(21)
接下來(lái)討論非高斯激勵(lì)下結(jié)構(gòu)響應(yīng)的幅值分布特性。隨機(jī)過(guò)程x(t)通過(guò)線性系統(tǒng)H(f)后的輸出在時(shí)域可表示為
(22)
式中，h(t)為系統(tǒng)H(f)的沖擊響應(yīng)函數(shù)，上述積分可用極限和形式表示，即
(23)
式中：x(τk)為隨機(jī)變量；Δτk為取樣時(shí)間間隔。根據(jù)中心極限定理，大量統(tǒng)計(jì)獨(dú)立的隨機(jī)變量之和的分布趨于高斯分布。如果輸出隨機(jī)過(guò)程y(t)在任意時(shí)刻t上皆為大量獨(dú)立隨機(jī)變量之和，則y(t)便接近于高斯分布。顯然這里要求兩個(gè)條件：一個(gè)是隨機(jī)變量必須相互獨(dú)立，另一個(gè)是獨(dú)立隨機(jī)變量要累加求和。
現(xiàn)在來(lái)考慮第二個(gè)條件。眾所周知，當(dāng)寬帶隨機(jī)信號(hào)作用于窄帶系統(tǒng)(如圖10所示濾波器頻響結(jié)構(gòu))時(shí)，由于系統(tǒng)有惰性，不能立即對(duì)信號(hào)作出響應(yīng)，它需要一定的建立時(shí)間ts，而ts與系統(tǒng)帶寬BWH成反比，即ts∝1/BWH。這樣，BWH越小，則ts越大，對(duì)信號(hào)響應(yīng)的時(shí)間越長(zhǎng)，對(duì)隨機(jī)輸入各個(gè)取樣(隨機(jī)變量)累積時(shí)間也越長(zhǎng)，于是當(dāng)各個(gè)取樣相互獨(dú)立，且累積時(shí)間又足夠長(zhǎng)，即滿足ts?Δτk時(shí)，則y(t)趨于高斯分布。反之，若非高斯隨機(jī)過(guò)程作用于線性系統(tǒng)，而系統(tǒng)的通頻帶較寬時(shí)，這時(shí)ts較小，若小到ts?Δτk，則輸入隨機(jī)過(guò)程通過(guò)系統(tǒng)后失真很小，于是輸出隨機(jī)過(guò)程y(t)的分布將接近原輸入隨機(jī)過(guò)程x(t)的分布，即為非高斯分布。
綜上所述，當(dāng)τx?Δτk?ts時(shí)，或簡(jiǎn)化為τx?ts時(shí)，系統(tǒng)在非高斯隨機(jī)輸入下，其輸出接近高斯分布。由于τx∝1/BWX及ts∝1/BWH，因此τx?ts也就意味著B(niǎo)WX?BWH。于是上述結(jié)論又可表述為線性系統(tǒng)輸入隨機(jī)過(guò)程的有效頻譜帶寬遠(yuǎn)大于系統(tǒng)帶寬時(shí)，輸出隨機(jī)過(guò)程的幅值分布將趨于高斯分布，而與輸入隨機(jī)過(guò)程是否為高斯分布無(wú)關(guān)；換言之，當(dāng)輸入隨機(jī)過(guò)程的有效頻譜帶寬接近或小于系統(tǒng)帶寬時(shí)，如果輸入是非高斯分布則系統(tǒng)輸出也將為非高斯分布。
根據(jù)上述分析結(jié)論，可采用下式來(lái)描述應(yīng)力響應(yīng)的峭度值Ks和加速度激勵(lì)的峭度Ka、帶寬BWa以及結(jié)構(gòu)本身的帶寬BWH之間的關(guān)系
(24)
式中，參數(shù)β為比例系數(shù)。
綜合式(15)和式(24)可得非高斯修正因子λ的表達(dá)式為
(25)
從式(25)可以看出，非高斯隨機(jī)振動(dòng)激勵(lì)的峭度和帶寬對(duì)結(jié)構(gòu)應(yīng)力響應(yīng)的非高斯特性影響比較明顯，從而對(duì)結(jié)構(gòu)的振動(dòng)疲勞壽命也會(huì)產(chǎn)生明顯影響。
將式(25)代入式(14)，并令D=1，可得非高斯隨機(jī)振動(dòng)激勵(lì)下的結(jié)構(gòu)振動(dòng)疲勞壽命TNG為
(26)
令ε=αβ，則式(26)可進(jìn)一步簡(jiǎn)化為
(27)
有了式(27)，就可以將結(jié)構(gòu)振動(dòng)疲勞壽命與振動(dòng)激勵(lì)的諸多特性以及結(jié)構(gòu)本身動(dòng)力學(xué)特性都緊密地聯(lián)系起來(lái)，用于定量設(shè)計(jì)振動(dòng)加速試驗(yàn)將十分方便。
從式(27)可以看出：一旦確定了試件結(jié)構(gòu)的材料、外形和尺寸，f1、ξ和BWH也隨之確定，可視為已知量；一旦確定了振動(dòng)激勵(lì)，Ga(f1)和BWa也隨之確定，也可以視為已知量。這樣，式(27)中還有3個(gè)未知量待求：b、k1和ε。下面依次探討這3個(gè)未知量的求解方法。
首先討論如何求解參數(shù)b。
對(duì)式(13)兩邊取對(duì)數(shù)，可得
(28)
(29)
這樣就可以通過(guò)進(jìn)行幾組高斯隨機(jī)振動(dòng)試驗(yàn)，得到若干組(X1,Y1)的值，然后進(jìn)行曲線擬合，就可以得到參數(shù)b的估計(jì)值。
接下來(lái)討論如何求解參數(shù)k1。
對(duì)式(10)進(jìn)行變換可得
(30)
Y2=k1X2
(31)
同樣，根據(jù)幾組高斯振動(dòng)疲勞試驗(yàn)的結(jié)果，得到若干組(X2,Y2)的值，然后進(jìn)行曲線擬合，就可以得到參數(shù)k1的估計(jì)值。
最后討論如何求解參數(shù)ε。
對(duì)式(27)進(jìn)行變換可得
(32)
Y3=εX3
(33)
這樣可以通過(guò)進(jìn)行幾組非高斯隨機(jī)振動(dòng)試驗(yàn)，可以得到若干組(X3,Y3)的值，然后進(jìn)行曲線擬合，就可以得到參數(shù)ε的估計(jì)值。
至此，就求出了式(27)中所有的未知參數(shù)，從而就可以預(yù)計(jì)出不同試驗(yàn)條件下的結(jié)構(gòu)振動(dòng)疲勞壽命，并可據(jù)此推斷出試件在某實(shí)際服役振動(dòng)環(huán)境下的疲勞壽命。
2 試驗(yàn)研究
2.1 試驗(yàn)對(duì)象
某型變壓器電子組件在振動(dòng)環(huán)境下其電容管腳容易出現(xiàn)疲勞斷裂，如圖2所示。
(a)
(b)
2.2 試驗(yàn)方案及結(jié)果
根據(jù)前述參數(shù)估計(jì)方法，分別設(shè)計(jì)了一組高斯和一組超高斯振動(dòng)疲勞試驗(yàn)，如表1和表2所示。其中高斯振動(dòng)疲勞試驗(yàn)結(jié)果用于估計(jì)參數(shù)b和k1，超高斯振動(dòng)疲勞試驗(yàn)結(jié)果用于估計(jì)參數(shù)ε。在得到三個(gè)未知參數(shù)的估計(jì)值后，再選擇如表3所示的試驗(yàn)條件進(jìn)行試件的常規(guī)壽命試驗(yàn)，并與采用疲勞加速試驗(yàn)?zāi)Ｐ皖A(yù)測(cè)的結(jié)果進(jìn)行比較。
表1 高斯振動(dòng)疲勞試驗(yàn)
表2 超高斯振動(dòng)疲勞試驗(yàn)
表3 疲勞壽命預(yù)測(cè)驗(yàn)證試驗(yàn)
3 結(jié) 論
本文通過(guò)理論分析，提出了一種能夠適用于超高斯隨機(jī)振動(dòng)加速試驗(yàn)的數(shù)學(xué)模型，并給出了模型中未知參數(shù)的具體求解方法。該模型系統(tǒng)全面地將結(jié)構(gòu)振動(dòng)疲勞壽命與振動(dòng)激勵(lì)的諸多特性以及結(jié)構(gòu)本身動(dòng)力學(xué)特性都緊密地聯(lián)系起來(lái)，用于指導(dǎo)超高斯振動(dòng)疲勞加速試驗(yàn)剖面的定量設(shè)計(jì)將十分方便，具有很好的工程實(shí)用性。
[1] 王冬梅, 謝勁松. 隨機(jī)振動(dòng)試驗(yàn)加速因子的計(jì)算方法[J]. 環(huán)境技術(shù), 2010, 28(2): 47-51.
WANG Dongmei, XIE Jinsong. Acceleration factor calculation method of random vibration test[J]. Environment Technology, 2010, 28(2): 47-51.
[2] 李奇志, 陳國(guó)平, 王明旭, 等. 振動(dòng)加速因子試驗(yàn)方法研究[J]. 振動(dòng)、測(cè)試與診斷, 2013, 33(1): 35-39.
LI Qizhi, CHEN Guoping, WANG Mingxu, et al. Test method of vibration acceleration factor[J]. Journal of Vibration,Measurement&Diagnosis, 2013, 33(1): 35-39.
[3] 朱學(xué)旺, 張思箭, 寧佐貴, 等. 寬帶隨機(jī)振動(dòng)試驗(yàn)條件的加速因子[J]. 環(huán)境技術(shù), 2014(6): 17-20.
ZHU Xuewang, ZHANG Sijian, NING Zuogui, et al. Accelerated testing factor on wideband random vibration[J]. Environment Technology, 2014(6):17-20.
[4] ALLEGRI G, ZHANG X. On the inverse power laws for accelerated random fatigue testing[J]. International Journal of Fatigue, 2008, 30(6): 967-977.
[5] POTHULA A, GUPTA A, KATHAWATE G R. Fatigue failure in random vibration and accelerated testing[J]. Journal of Vibration and Control, 2012, 18(8): 1199-1206.
[6] JOHN V B. Fatigue damage spectrum—a new tool to accelerate vibration testing[J]. Sound and Vibration, 2015, 49(3): 15-17.
[7] CRANDALL S H, MARK W D. Random vibration in mechanical systems[M]. New York: Academic Press, 1963.
[8] KERN D L. Dynamic environmental criteria: NASA-HDBK-7005[S]. [S.l.]: NASA Headquarters, 2001.
Super-Gaussian random vibration fatigue accelerated testing model
JIANG Yu, TAO Junyong, CHEN Xun
(Science and Technology on Integrated Logistics Support Technology, College of Mechatronic Engineering&Automation, National University of Defense Technology, Changsha 410073, China)
Here, a new mathematical model of accelerated vibration fatigue testing was established under super-Gaussian random vibration excitation using theoretical analysis, and detailed solving methods were presented for unknown parameters in the model. The model was verified with the actual tests. It was shown that the model can be effectively used for the practical quantitative design of accelerated vibration fatigue tests, it can be applied to assess the fatigue life and reliability of engineering structures under random vibration environment.
super-Gaussian; random vibration; vibration fatigue; accelerated model
國(guó)家自然科學(xué)基金項(xiàng)目(50905181)
2016-01-08 修改稿收到日期：2016-07-08
蔣瑜 男，博士，副教授，1977年生
O324
A
10.13465/j.cnki.jvs.2017.09.038
猜你喜歡
一半模型童話王國(guó)·奇妙邏輯推理(2024年5期)2024-06-19 16:03:38 振動(dòng)的思考科學(xué)大眾(2023年17期)2023-10-26 07:39:14 《形而上學(xué)》△卷的結(jié)構(gòu)和位置哲學(xué)評(píng)論(2021年2期)2021-08-22 01:53:34 重要模型『一線三等角』中學(xué)生數(shù)理化·七年級(jí)數(shù)學(xué)人教版(2020年10期)2020-11-26 08:24:50 振動(dòng)與頻率天天愛(ài)科學(xué)(2020年6期)2020-09-10 07:22:44 重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布數(shù)學(xué)物理學(xué)報(bào)(2020年2期)2020-06-02 11:29:24 論結(jié)構(gòu)中華詩(shī)詞(2019年7期)2019-11-25 01:43:04 中立型Emden-Fowler微分方程的振動(dòng)性數(shù)學(xué)物理學(xué)報(bào)(2017年6期)2018-01-22 02:26:40 論《日出》的結(jié)構(gòu)影視與戲劇評(píng)論(2016年0期)2016-11-23 05:26:01 3D打印中的模型分割與打包光學(xué)精密工程(2016年6期)2016-11-07 09:07:19
蔣 瑜， 陶俊勇， 陳 循
(國(guó)防科技大學(xué) 機(jī)電工程與自動(dòng)化學(xué)院裝備綜合保障技術(shù)重點(diǎn)實(shí)驗(yàn)室， 長(zhǎng)沙 410073)
超高斯隨機(jī)振動(dòng)疲勞加速試驗(yàn)?zāi)Ｐ脱芯?/p>
蔣 瑜， 陶俊勇， 陳 循
(國(guó)防科技大學(xué) 機(jī)電工程與自動(dòng)化學(xué)院裝備綜合保障技術(shù)重點(diǎn)實(shí)驗(yàn)室， 長(zhǎng)沙 410073)
通過(guò)理論分析，針對(duì)超高斯隨機(jī)振動(dòng)激勵(lì)建立了相應(yīng)的振動(dòng)疲勞加速試驗(yàn)數(shù)學(xué)模型，給出了模型中未知參數(shù)的具體求解方法，并通過(guò)實(shí)際試驗(yàn)進(jìn)行了驗(yàn)證。該模型十分便于進(jìn)行振動(dòng)疲勞加速試驗(yàn)的定量設(shè)計(jì)，具有很好的工程實(shí)用性，可用于評(píng)估隨機(jī)振動(dòng)環(huán)境下工程結(jié)構(gòu)的疲勞壽命與可靠性。
超高斯； 隨機(jī)振動(dòng)； 振動(dòng)疲勞； 加速模型
振動(dòng)引起的疲勞問(wèn)題作為多個(gè)工程領(lǐng)域廣泛存在的一個(gè)共性問(wèn)題，嚴(yán)重危及重大裝備及結(jié)構(gòu)的可靠性和安全性。如果能夠提前準(zhǔn)確預(yù)測(cè)結(jié)構(gòu)的振動(dòng)疲勞壽命，就能在發(fā)生災(zāi)難性事故之前及時(shí)預(yù)知并采取相應(yīng)的措施，并可為其定壽、延壽提供科學(xué)的依據(jù)，最大限度地發(fā)揮裝備效益。因此，準(zhǔn)確預(yù)測(cè)重大裝備和工程結(jié)構(gòu)在復(fù)雜隨機(jī)動(dòng)態(tài)載荷作用下的振動(dòng)疲勞壽命是提高其可靠性和安全性的關(guān)鍵技術(shù)，被國(guó)家自然科學(xué)基金委員會(huì)列入《機(jī)械工程學(xué)科發(fā)展戰(zhàn)略報(bào)告(2011～2020)》。
當(dāng)前工程中進(jìn)行結(jié)構(gòu)振動(dòng)疲勞壽命分析和預(yù)測(cè)時(shí)，由于缺乏充分的認(rèn)識(shí)和解決相關(guān)問(wèn)題的工具及方法，通常假設(shè)結(jié)構(gòu)承受的隨機(jī)載荷符合平穩(wěn)高斯分布。但是，在實(shí)際環(huán)境中，上述許多隨機(jī)載荷往往呈現(xiàn)比較明顯的非高斯特征，尤其是在工況惡劣或者極端環(huán)境下。近年來(lái)的一些研究結(jié)果表明，隨機(jī)載荷的非高斯特性對(duì)結(jié)構(gòu)振動(dòng)疲勞有著不可忽視的影響，某些情況下會(huì)加速結(jié)構(gòu)的疲勞失效?；诟咚辜僭O(shè)計(jì)算振動(dòng)疲勞損傷往往會(huì)得到偏大的振動(dòng)疲勞壽命估計(jì)結(jié)果，給裝備服役或使用階段埋下巨大的安全隱患。
同時(shí)，隨著結(jié)構(gòu)可靠性水平的提高，結(jié)構(gòu)的振動(dòng)疲勞壽命越來(lái)越長(zhǎng)，為了能夠在實(shí)驗(yàn)室驗(yàn)證其壽命是否達(dá)到要求，加速試驗(yàn)成為必然的選擇。要想通過(guò)加速試驗(yàn)得到結(jié)構(gòu)在實(shí)際服役振動(dòng)環(huán)境下的疲勞壽命，正確的加速試驗(yàn)統(tǒng)計(jì)模型是關(guān)鍵。國(guó)內(nèi)外目前關(guān)于振動(dòng)加速試驗(yàn)?zāi)Ｐ偷难芯坎欢啵绕涫菍?duì)非高斯隨機(jī)振動(dòng)加速試驗(yàn)。王冬梅等[1]對(duì)振動(dòng)加速試驗(yàn)的逆冪律模型進(jìn)行了推導(dǎo)，探討了其適用范圍，指出其適用于窄帶高斯載荷，不適用于寬帶高斯和非高斯載荷。李奇志等[2]提出通過(guò)試驗(yàn)的方式獲得振動(dòng)試驗(yàn)的加速因子，認(rèn)為振動(dòng)加速試驗(yàn)的逆冪律模型對(duì)平穩(wěn)窄帶和寬帶高斯隨機(jī)過(guò)程均是適用的。朱學(xué)旺等[3]應(yīng)用基于窄帶模型的修正方法得到了寬帶隨機(jī)振動(dòng)試驗(yàn)加速因子計(jì)算的通用表達(dá)式，認(rèn)為基于窄帶模型的加速因子表達(dá)式對(duì)于比例載荷的寬帶隨機(jī)振動(dòng)也是適用的，而對(duì)于非比例載荷，則需要應(yīng)用其提出的通用表達(dá)式才可以獲得。Allegri等[4]研究了適用于平穩(wěn)寬帶高斯隨機(jī)振動(dòng)加速試驗(yàn)的逆冪律模型，其主要目的是用來(lái)評(píng)估相對(duì)疲勞損傷，而非用來(lái)精確預(yù)計(jì)壽命。Pothula[5]研究了附著不同阻尼材料的鋁合金梁在高斯隨機(jī)振動(dòng)加速試驗(yàn)中的疲勞壽命差異，探討了結(jié)構(gòu)阻尼對(duì)振動(dòng)試驗(yàn)加速因子的影響。John[6]提出應(yīng)用Fatigue Damage Spectrum(FDS)來(lái)進(jìn)行振動(dòng)加速試驗(yàn)設(shè)計(jì)與評(píng)估，但是只考慮了振動(dòng)激勵(lì)的功率譜密度PSD (Power Spectral Density)，只適用于高斯振動(dòng)加速試驗(yàn)。
總的來(lái)說(shuō)，目前國(guó)內(nèi)外開(kāi)展的振動(dòng)疲勞加速試驗(yàn)研究主要針對(duì)高斯振動(dòng)，尤其是實(shí)際進(jìn)行的非高斯振動(dòng)加速試驗(yàn)很少，需要進(jìn)一步探索適用于非高斯隨機(jī)振動(dòng)的疲勞加速試驗(yàn)定量模型。
1 理論分析
工程上常用偏斜度S和峭度K這兩個(gè)參數(shù)來(lái)描述非高斯隨機(jī)過(guò)程X，定義為
(1)
高斯隨機(jī)過(guò)程的偏斜度值等于0，峭度值等于3；而非高斯隨機(jī)過(guò)程的峭度值肯定不等于3，偏斜度值可以等于零也可以不等于0。偏斜度用來(lái)描述隨機(jī)過(guò)程幅值概率密度曲線偏離對(duì)稱分布的程度，偏斜度值不為零表示服從非對(duì)稱分布。峭度是描述隨機(jī)過(guò)程幅值概率密度曲線拖尾分布特征的參數(shù)，它不僅可用來(lái)區(qū)分高斯和非高斯隨機(jī)過(guò)程，而且還可進(jìn)一步將非高斯隨機(jī)過(guò)程區(qū)分為亞高斯和超高斯隨機(jī)過(guò)程，其中亞高斯隨機(jī)過(guò)程的K<3，超高斯隨機(jī)過(guò)程的K>3。工程中常見(jiàn)的非高斯振動(dòng)信號(hào)往往是具有尖峰分布的對(duì)稱超高斯信號(hào)，因此本文研究的非高斯振動(dòng)疲勞主要針對(duì)超高斯。
1.1 高斯振動(dòng)疲勞加速試驗(yàn)?zāi)Ｐ?/p>
首先，從最基本的描述材料疲勞現(xiàn)象的S～N曲線出發(fā)。通常理想的S～N曲線的數(shù)學(xué)表達(dá)式為
N=cS-b
(2)
式中：S為應(yīng)力幅值；N為引起失效的應(yīng)力循環(huán)次數(shù)；b、c為材料雙對(duì)數(shù)S～N曲線中的常數(shù)，且b的取值范圍一般為6～25。
根據(jù)著名的疲勞損傷累積Miner準(zhǔn)則，不同幅值應(yīng)力共同作用下造成的疲勞損傷為
(3)
式中：ni是幅值為Si的應(yīng)力作用循環(huán)次數(shù)；Ni為試件在幅值為Si的應(yīng)力作用下至疲勞失效的循環(huán)次數(shù)；D為累積疲勞損傷(一般認(rèn)為當(dāng)D=1時(shí)發(fā)生疲勞失效)。
將式(2)代入式(3)可得
(4)
對(duì)于連續(xù)的隨機(jī)應(yīng)力時(shí)間歷程，式(4)可寫成如下積分的形式[7]
(5)
當(dāng)隨機(jī)應(yīng)力響應(yīng)接近平穩(wěn)窄帶高斯分布時(shí)，根據(jù)隨機(jī)過(guò)程理論，應(yīng)力幅值概率密度函數(shù)p(S)將近似服從瑞利(Rayleigh)分布
(6)
式中：σS為應(yīng)力的均方根值(即標(biāo)準(zhǔn)差)。將式(6)代入式(5)進(jìn)行積分可得
(7)
式中，Γ為Gamma函數(shù)。
進(jìn)行振動(dòng)試驗(yàn)的產(chǎn)品或試件一般可近似看成一個(gè)線性系統(tǒng)，而振動(dòng)試驗(yàn)設(shè)備產(chǎn)生的振動(dòng)激勵(lì)可看作該系統(tǒng)的輸入。工程實(shí)踐表明，一般結(jié)構(gòu)件的阻尼比ξ通常遠(yuǎn)小于1，可以取到0.05以下這樣小的值。根據(jù)文獻(xiàn)[8]，在滿足小阻尼線性系統(tǒng)的假設(shè)前提下，σS的近似計(jì)算式為
(8)
式中：f1為試件結(jié)構(gòu)的一階固有頻率；ξ為等效的阻尼比(一般假設(shè)ξ≤0.1)；k為與試件材料相關(guān)的比例常數(shù)；Ga(f1)為輸入振動(dòng)激勵(lì)的加速度功率譜密度在試件固有頻率f1處的量值。
(9)
由于試件結(jié)構(gòu)的傳遞函數(shù)一般都類似一個(gè)窄帶濾波器，因此可以認(rèn)為試件在平穩(wěn)高斯隨機(jī)激勵(lì)作用下(無(wú)論是寬帶還是窄帶)，其應(yīng)力響應(yīng)為平穩(wěn)窄帶高斯分布，可以采用式(9)來(lái)估計(jì)其疲勞累積損傷。
一般認(rèn)為D=1時(shí)結(jié)構(gòu)發(fā)生疲勞失效，根據(jù)式(9)可得高斯隨機(jī)振動(dòng)激勵(lì)下的結(jié)構(gòu)振動(dòng)疲勞壽命TG為
(10)
對(duì)常用工程材料，b的取值范圍一般為4～25。從式(10)可以看出，對(duì)高斯隨機(jī)振動(dòng)激勵(lì)，當(dāng)結(jié)構(gòu)動(dòng)力學(xué)特性參數(shù)如f1、ξ固定時(shí)，結(jié)構(gòu)振動(dòng)疲勞壽命T與高斯隨機(jī)振動(dòng)激勵(lì)的功率譜密度在結(jié)構(gòu)一階固有頻率處的量值大小Ga(f1)成反比關(guān)系。隨著Ga(f1)的增加(或減小)，T按指數(shù)規(guī)律迅速衰減(或增長(zhǎng))。
根據(jù)式(10)進(jìn)一步可得a1、a2兩種不同高斯隨機(jī)加速度激勵(lì)下的結(jié)構(gòu)振動(dòng)疲勞壽命分別為
(11)
(12)
根據(jù)式(11)和式(12)可得
(13)
顯然，式(13)就是傳統(tǒng)文獻(xiàn)中用于描述高斯振動(dòng)加速試驗(yàn)的逆冪律模型。從上述推導(dǎo)過(guò)程來(lái)看，只要工程結(jié)構(gòu)系統(tǒng)的頻響函數(shù)是一個(gè)窄帶濾波器(實(shí)際中大多數(shù)工程結(jié)構(gòu)均滿足這一條件)，并且結(jié)構(gòu)隨機(jī)響應(yīng)服從高斯分布，則就可以運(yùn)用式(13)對(duì)振動(dòng)加速試驗(yàn)進(jìn)行建模。因此，式(13)只適用于描述高斯振動(dòng)疲勞加速試驗(yàn)，下面繼續(xù)討論非高斯振動(dòng)疲勞加速試驗(yàn)?zāi)Ｐ汀?/p>
1.2 超高斯振動(dòng)疲勞加速試驗(yàn)?zāi)Ｐ?/p>
當(dāng)隨機(jī)應(yīng)力響應(yīng)接近平穩(wěn)窄帶非高斯分布時(shí)，可考慮在式(9)的基礎(chǔ)上增加一個(gè)非高斯修正因子λ來(lái)描述應(yīng)力響應(yīng)的峭度值對(duì)振動(dòng)疲勞累積損傷的影響
(14)
非高斯修正因子λ和應(yīng)力響應(yīng)的峭度值Ks直接相關(guān)，可用式(15)來(lái)描述
(15)
式中，參數(shù)α為比例系數(shù)。
顯然，當(dāng)應(yīng)力響應(yīng)為高斯即Ks=3時(shí)，λ=1，式(14)變成式(9)；當(dāng)應(yīng)力響應(yīng)為超高斯即Ks>3時(shí)，λ>1，意味著應(yīng)力響應(yīng)的超高斯特性會(huì)加快疲勞累積損傷進(jìn)程。
下面根據(jù)隨機(jī)振動(dòng)理論分析影響應(yīng)力響應(yīng)峭度值Ks的因素。由于結(jié)構(gòu)的一階模態(tài)對(duì)結(jié)構(gòu)響應(yīng)起決定性作用，因此建立振動(dòng)臺(tái)基礎(chǔ)激勵(lì)作用下的單自由度系統(tǒng)模型進(jìn)行分析，如下圖1所示。
經(jīng)過(guò)推導(dǎo)可得加速度響應(yīng)y與振動(dòng)臺(tái)臺(tái)面的基礎(chǔ)加速度激勵(lì)x之間的傳遞函數(shù)為
(16)
(17)
(18)
式中：ω1=2πf1，f1為試件結(jié)構(gòu)的一階固有頻率；ξ為阻尼比。這兩個(gè)參數(shù)表征了試件結(jié)構(gòu)本身的動(dòng)力學(xué)特性，并且試件結(jié)構(gòu)系統(tǒng)的通頻帶BWH也由這兩個(gè)參數(shù)決定，即
BWH=2ξf1
(19)
圖1 振動(dòng)臺(tái)基礎(chǔ)激勵(lì)作用下的單自由度模型
由于實(shí)際工程結(jié)構(gòu)的阻尼比ξ往往較小，而一階固有頻率f1往往不大，所以大多數(shù)工程結(jié)構(gòu)的通頻帶BWH往往也不大，也就是說(shuō)可以視為一個(gè)窄帶濾波器。
下面先對(duì)響應(yīng)帶寬進(jìn)行分析。設(shè)激勵(lì)x的功率譜密度記作X(f)，響應(yīng)y的功率譜密度記作Y(f)，系統(tǒng)頻響函數(shù)為H(f)。根據(jù)線性系統(tǒng)和隨機(jī)過(guò)程理論有
(20)
由此容易得出響應(yīng)的有效頻譜帶寬
(21)
接下來(lái)討論非高斯激勵(lì)下結(jié)構(gòu)響應(yīng)的幅值分布特性。隨機(jī)過(guò)程x(t)通過(guò)線性系統(tǒng)H(f)后的輸出在時(shí)域可表示為
(22)
式中，h(t)為系統(tǒng)H(f)的沖擊響應(yīng)函數(shù)，上述積分可用極限和形式表示，即
(23)
式中：x(τk)為隨機(jī)變量；Δτk為取樣時(shí)間間隔。根據(jù)中心極限定理，大量統(tǒng)計(jì)獨(dú)立的隨機(jī)變量之和的分布趨于高斯分布。如果輸出隨機(jī)過(guò)程y(t)在任意時(shí)刻t上皆為大量獨(dú)立隨機(jī)變量之和，則y(t)便接近于高斯分布。顯然這里要求兩個(gè)條件：一個(gè)是隨機(jī)變量必須相互獨(dú)立，另一個(gè)是獨(dú)立隨機(jī)變量要累加求和。
現(xiàn)在來(lái)考慮第二個(gè)條件。眾所周知，當(dāng)寬帶隨機(jī)信號(hào)作用于窄帶系統(tǒng)(如圖10所示濾波器頻響結(jié)構(gòu))時(shí)，由于系統(tǒng)有惰性，不能立即對(duì)信號(hào)作出響應(yīng)，它需要一定的建立時(shí)間ts，而ts與系統(tǒng)帶寬BWH成反比，即ts∝1/BWH。這樣，BWH越小，則ts越大，對(duì)信號(hào)響應(yīng)的時(shí)間越長(zhǎng)，對(duì)隨機(jī)輸入各個(gè)取樣(隨機(jī)變量)累積時(shí)間也越長(zhǎng)，于是當(dāng)各個(gè)取樣相互獨(dú)立，且累積時(shí)間又足夠長(zhǎng)，即滿足ts?Δτk時(shí)，則y(t)趨于高斯分布。反之，若非高斯隨機(jī)過(guò)程作用于線性系統(tǒng)，而系統(tǒng)的通頻帶較寬時(shí)，這時(shí)ts較小，若小到ts?Δτk，則輸入隨機(jī)過(guò)程通過(guò)系統(tǒng)后失真很小，于是輸出隨機(jī)過(guò)程y(t)的分布將接近原輸入隨機(jī)過(guò)程x(t)的分布，即為非高斯分布。
綜上所述，當(dāng)τx?Δτk?ts時(shí)，或簡(jiǎn)化為τx?ts時(shí)，系統(tǒng)在非高斯隨機(jī)輸入下，其輸出接近高斯分布。由于τx∝1/BWX及ts∝1/BWH，因此τx?ts也就意味著B(niǎo)WX?BWH。于是上述結(jié)論又可表述為線性系統(tǒng)輸入隨機(jī)過(guò)程的有效頻譜帶寬遠(yuǎn)大于系統(tǒng)帶寬時(shí)，輸出隨機(jī)過(guò)程的幅值分布將趨于高斯分布，而與輸入隨機(jī)過(guò)程是否為高斯分布無(wú)關(guān)；換言之，當(dāng)輸入隨機(jī)過(guò)程的有效頻譜帶寬接近或小于系統(tǒng)帶寬時(shí)，如果輸入是非高斯分布則系統(tǒng)輸出也將為非高斯分布。
根據(jù)上述分析結(jié)論，可采用下式來(lái)描述應(yīng)力響應(yīng)的峭度值Ks和加速度激勵(lì)的峭度Ka、帶寬BWa以及結(jié)構(gòu)本身的帶寬BWH之間的關(guān)系
(24)
式中，參數(shù)β為比例系數(shù)。
綜合式(15)和式(24)可得非高斯修正因子λ的表達(dá)式為
(25)
從式(25)可以看出，非高斯隨機(jī)振動(dòng)激勵(lì)的峭度和帶寬對(duì)結(jié)構(gòu)應(yīng)力響應(yīng)的非高斯特性影響比較明顯，從而對(duì)結(jié)構(gòu)的振動(dòng)疲勞壽命也會(huì)產(chǎn)生明顯影響。
將式(25)代入式(14)，并令D=1，可得非高斯隨機(jī)振動(dòng)激勵(lì)下的結(jié)構(gòu)振動(dòng)疲勞壽命TNG為
(26)
令ε=αβ，則式(26)可進(jìn)一步簡(jiǎn)化為
(27)
有了式(27)，就可以將結(jié)構(gòu)振動(dòng)疲勞壽命與振動(dòng)激勵(lì)的諸多特性以及結(jié)構(gòu)本身動(dòng)力學(xué)特性都緊密地聯(lián)系起來(lái)，用于定量設(shè)計(jì)振動(dòng)加速試驗(yàn)將十分方便。
從式(27)可以看出：一旦確定了試件結(jié)構(gòu)的材料、外形和尺寸，f1、ξ和BWH也隨之確定，可視為已知量；一旦確定了振動(dòng)激勵(lì)，Ga(f1)和BWa也隨之確定，也可以視為已知量。這樣，式(27)中還有3個(gè)未知量待求：b、k1和ε。下面依次探討這3個(gè)未知量的求解方法。
首先討論如何求解參數(shù)b。
對(duì)式(13)兩邊取對(duì)數(shù)，可得
(28)
(29)
這樣就可以通過(guò)進(jìn)行幾組高斯隨機(jī)振動(dòng)試驗(yàn)，得到若干組(X1,Y1)的值，然后進(jìn)行曲線擬合，就可以得到參數(shù)b的估計(jì)值。
接下來(lái)討論如何求解參數(shù)k1。
對(duì)式(10)進(jìn)行變換可得
(30)
Y2=k1X2
(31)
同樣，根據(jù)幾組高斯振動(dòng)疲勞試驗(yàn)的結(jié)果，得到若干組(X2,Y2)的值，然后進(jìn)行曲線擬合，就可以得到參數(shù)k1的估計(jì)值。
最后討論如何求解參數(shù)ε。
對(duì)式(27)進(jìn)行變換可得
(32)
Y3=εX3
(33)
這樣可以通過(guò)進(jìn)行幾組非高斯隨機(jī)振動(dòng)試驗(yàn)，可以得到若干組(X3,Y3)的值，然后進(jìn)行曲線擬合，就可以得到參數(shù)ε的估計(jì)值。
至此，就求出了式(27)中所有的未知參數(shù)，從而就可以預(yù)計(jì)出不同試驗(yàn)條件下的結(jié)構(gòu)振動(dòng)疲勞壽命，并可據(jù)此推斷出試件在某實(shí)際服役振動(dòng)環(huán)境下的疲勞壽命。
2 試驗(yàn)研究
2.1 試驗(yàn)對(duì)象
某型變壓器電子組件在振動(dòng)環(huán)境下其電容管腳容易出現(xiàn)疲勞斷裂，如圖2所示。
(a)
(b)
2.2 試驗(yàn)方案及結(jié)果
根據(jù)前述參數(shù)估計(jì)方法，分別設(shè)計(jì)了一組高斯和一組超高斯振動(dòng)疲勞試驗(yàn)，如表1和表2所示。其中高斯振動(dòng)疲勞試驗(yàn)結(jié)果用于估計(jì)參數(shù)b和k1，超高斯振動(dòng)疲勞試驗(yàn)結(jié)果用于估計(jì)參數(shù)ε。在得到三個(gè)未知參數(shù)的估計(jì)值后，再選擇如表3所示的試驗(yàn)條件進(jìn)行試件的常規(guī)壽命試驗(yàn)，并與采用疲勞加速試驗(yàn)?zāi)Ｐ皖A(yù)測(cè)的結(jié)果進(jìn)行比較。
表1 高斯振動(dòng)疲勞試驗(yàn)
表2 超高斯振動(dòng)疲勞試驗(yàn)
表3 疲勞壽命預(yù)測(cè)驗(yàn)證試驗(yàn)
3 結(jié) 論
本文通過(guò)理論分析，提出了一種能夠適用于超高斯隨機(jī)振動(dòng)加速試驗(yàn)的數(shù)學(xué)模型，并給出了模型中未知參數(shù)的具體求解方法。該模型系統(tǒng)全面地將結(jié)構(gòu)振動(dòng)疲勞壽命與振動(dòng)激勵(lì)的諸多特性以及結(jié)構(gòu)本身動(dòng)力學(xué)特性都緊密地聯(lián)系起來(lái)，用于指導(dǎo)超高斯振動(dòng)疲勞加速試驗(yàn)剖面的定量設(shè)計(jì)將十分方便，具有很好的工程實(shí)用性。
[1] 王冬梅, 謝勁松. 隨機(jī)振動(dòng)試驗(yàn)加速因子的計(jì)算方法[J]. 環(huán)境技術(shù), 2010, 28(2): 47-51.
WANG Dongmei, XIE Jinsong. Acceleration factor calculation method of random vibration test[J]. Environment Technology, 2010, 28(2): 47-51.
[2] 李奇志, 陳國(guó)平, 王明旭, 等. 振動(dòng)加速因子試驗(yàn)方法研究[J]. 振動(dòng)、測(cè)試與診斷, 2013, 33(1): 35-39.
LI Qizhi, CHEN Guoping, WANG Mingxu, et al. Test method of vibration acceleration factor[J]. Journal of Vibration,Measurement&Diagnosis, 2013, 33(1): 35-39.
[3] 朱學(xué)旺, 張思箭, 寧佐貴, 等. 寬帶隨機(jī)振動(dòng)試驗(yàn)條件的加速因子[J]. 環(huán)境技術(shù), 2014(6): 17-20.
ZHU Xuewang, ZHANG Sijian, NING Zuogui, et al. Accelerated testing factor on wideband random vibration[J]. Environment Technology, 2014(6):17-20.
[4] ALLEGRI G, ZHANG X. On the inverse power laws for accelerated random fatigue testing[J]. International Journal of Fatigue, 2008, 30(6): 967-977.
[5] POTHULA A, GUPTA A, KATHAWATE G R. Fatigue failure in random vibration and accelerated testing[J]. Journal of Vibration and Control, 2012, 18(8): 1199-1206.
[6] JOHN V B. Fatigue damage spectrum—a new tool to accelerate vibration testing[J]. Sound and Vibration, 2015, 49(3): 15-17.
[7] CRANDALL S H, MARK W D. Random vibration in mechanical systems[M]. New York: Academic Press, 1963.
[8] KERN D L. Dynamic environmental criteria: NASA-HDBK-7005[S]. [S.l.]: NASA Headquarters, 2001.
Super-Gaussian random vibration fatigue accelerated testing model
JIANG Yu, TAO Junyong, CHEN Xun
(Science and Technology on Integrated Logistics Support Technology, College of Mechatronic Engineering&Automation, National University of Defense Technology, Changsha 410073, China)
Here, a new mathematical model of accelerated vibration fatigue testing was established under super-Gaussian random vibration excitation using theoretical analysis, and detailed solving methods were presented for unknown parameters in the model. The model was verified with the actual tests. It was shown that the model can be effectively used for the practical quantitative design of accelerated vibration fatigue tests, it can be applied to assess the fatigue life and reliability of engineering structures under random vibration environment.
super-Gaussian; random vibration; vibration fatigue; accelerated model
國(guó)家自然科學(xué)基金項(xiàng)目(50905181)
2016-01-08 修改稿收到日期：2016-07-08
蔣瑜 男，博士，副教授，1977年生
O324
A
10.13465/j.cnki.jvs.2017.09.038
猜你喜歡
童話王國(guó)·奇妙邏輯推理(2024年5期)2024-06-19 16:03:38
科學(xué)大眾(2023年17期)2023-10-26 07:39:14
哲學(xué)評(píng)論(2021年2期)2021-08-22 01:53:34
中學(xué)生數(shù)理化·七年級(jí)數(shù)學(xué)人教版(2020年10期)2020-11-26 08:24:50
天天愛(ài)科學(xué)(2020年6期)2020-09-10 07:22:44
數(shù)學(xué)物理學(xué)報(bào)(2020年2期)2020-06-02 11:29:24
中華詩(shī)詞(2019年7期)2019-11-25 01:43:04
數(shù)學(xué)物理學(xué)報(bào)(2017年6期)2018-01-22 02:26:40
影視與戲劇評(píng)論(2016年0期)2016-11-23 05:26:01
光學(xué)精密工程(2016年6期)2016-11-07 09:07:19
