1、前言：

疲勞分析通常是在時(shí)域進(jìn)行，所有的輸入載荷和輸出應(yīng)力都是基于時(shí)間的信號；時(shí)域疲勞可以通過靜應(yīng)力分析或者模態(tài)瞬態(tài)法進(jìn)行分析；其中模態(tài)瞬態(tài)法一般用于需要考慮共振對疲勞的影響，載荷的加載頻率接近系統(tǒng)的共振頻率；在一些情況下，共振應(yīng)力和輸入載荷卻通過頻域信號來分析，通常用PSD功率譜密度來表達(dá)，基于的PSD頻域疲勞預(yù)測方法比時(shí)域疲勞預(yù)測方法有以下優(yōu)勢。

(1)時(shí)域所得損傷是取自對一段隨機(jī)變化信號的計(jì)數(shù)，因此通過時(shí)域方法獲得的損傷本身就是一個(gè)隨機(jī)變量，無法避免對所得的損傷結(jié)果進(jìn)行統(tǒng)計(jì)推斷。通常，用雨流計(jì)數(shù)法得到的零部件應(yīng)力幅值服從威布爾分布，均值服從正態(tài)分布。這些需要進(jìn)行循環(huán)計(jì)數(shù)，數(shù)據(jù)處理量非常大。而基于PSD的頻域分析方法計(jì)算簡單，不需要循環(huán)計(jì)數(shù)。
(2)隨機(jī)動(dòng)態(tài)應(yīng)力，在時(shí)域內(nèi)需要很長的信號記錄才能準(zhǔn)確地描述隨機(jī)響應(yīng)，同時(shí)處理長的時(shí)域信號非常困難，而得到頻域功率譜應(yīng)力信號則較為方便。
(3)用來進(jìn)行疲勞分析的頻域信號采樣率只要達(dá)到時(shí)域
信號采樣率的1／10就可以得到與用時(shí)域信號預(yù)測同樣精度的結(jié)果，頻域信號的讀取、儲(chǔ)存都比時(shí)域信號方便。

2、隨機(jī)振動(dòng)信號的特征
當(dāng)系統(tǒng)所受到的載荷信號是隨機(jī)不確定的時(shí)候，我們通常采用隨機(jī)振動(dòng)分析的進(jìn)行疲勞分析；假設(shè)所受載荷X(t)在 x和x+dx范圍內(nèi)，在一個(gè)總時(shí)長T的時(shí)間段內(nèi)，載荷出現(xiàn)的概率為fx（x）。



圖1



如果T足夠長，fx（x）可以通過下式表達(dá)：




相對的概率密度函數(shù)PDF可以通過總的時(shí)間段在帶寬X和X+dX段，如圖2所示




圖2



隨機(jī)信號X(t)的均值和均方根值可以表示為：



這里T是總的時(shí)長，當(dāng)ux=0時(shí)，就是隨機(jī)信號X(t)的均方根RMS
Fx（x）服從高斯正態(tài)分布



圖3









圖4

因?yàn)殡S機(jī)振動(dòng)激勵(lì)被假設(shè)為服從高斯正態(tài)分布，各幅值發(fā)生概率為：



圖5

基于這個(gè)特點(diǎn)，在實(shí)際計(jì)算中一般取3 sigma為計(jì)算的上限；高斯正態(tài)分布具有以下重要屬性：如果高級正態(tài)分布激勵(lì)作用在線性系統(tǒng)上，則輸出的激勵(lì)是不同的隨機(jī)過程，但是仍然服從另外一個(gè)高斯正態(tài)分布。
舉個(gè)隨機(jī)信號的處理過程的例子：



圖4



圖5

信號X(t)被分解為幾段，如圖4和圖5所示，對于平穩(wěn)的隨機(jī)過程x(t) ，時(shí)間歷程是非周期的，因此不能用傅里葉級數(shù)表示；而且x(t)是一個(gè)無限長的信號，不能通過傅里葉變換得到該隨機(jī)過程的頻域信息。這個(gè)困難可以通過對該隨機(jī)過程的自相關(guān)函數(shù)Rx(τ) 做傅里葉變換來解決。
由于X(t)是平穩(wěn)隨機(jī)信號其均值和標(biāo)準(zhǔn)差是獨(dú)立于時(shí)間t的，因此



相關(guān)系數(shù)



如果對隨機(jī)過程x(t)的零點(diǎn)進(jìn)行處理，使得該過程的平均值為0，并且假定x(t)不含有周期性分量，那么Rx(τ∞)=0，條件得到滿足。我們就可以得到自相關(guān)函數(shù)的Rx(τ)傅里葉變換和逆變換。



其中函數(shù)Sx(ω)稱為功率譜密度。我們令τ為0，則得到：





圖6

這是Sx(ω)最重要的一個(gè)特性，即功率譜密度曲線下的面積就是平穩(wěn)隨機(jī)過程x(t)的均方值，所以函數(shù)Sx(ω)又叫均方譜密度，其單位是(x的單位)2/(rad/s)，常見的有加速度隨機(jī)激勵(lì) 單位為(mm/s2)2/Hz或G2/Hz；速度隨機(jī)激勵(lì) 單位為(mm/s)2/Hz ；位移隨機(jī)激勵(lì)單位為 (mm)2/Hz。
在上面的推導(dǎo)過程中，圓頻率ω取值是從負(fù)無窮到正無窮，但我們研究振動(dòng)更習(xí)慣用頻率f而不是圓頻率ω，頻率f的取值應(yīng)該是從0到正無窮，單位應(yīng)該是Hz而不是rad/s。所以雙側(cè)譜密度Sx(ω)可以變換為一個(gè)等效的單側(cè)譜密度。







圖7



圖8

圖7中依次為窄帶、寬帶、白噪聲信號，圖8依次為其頻譜圖。

窄帶隨機(jī)過程的時(shí)間歷程類似于振幅和相位隨機(jī)變化的正弦波。根據(jù)窄帶隨機(jī)過程的PSD曲線，我們可以得到它的很多特性，如頻率成分和有效值等，還可以進(jìn)一步得到其峰值分布的信息，即組成這個(gè)過程的一系列正弦波的幅值分布信息。也就是說，我們可以依據(jù)應(yīng)力PSD曲線求得時(shí)間段T內(nèi)的應(yīng)力循環(huán)次數(shù)，以及應(yīng)力幅值在S和S+dS之間的概率Pp(S)dS 。由PSD求得應(yīng)力循環(huán)次數(shù)v和應(yīng)力幅值區(qū)間概率Pp(S)dS的公式推導(dǎo)比較復(fù)雜，建議讀者參考《隨機(jī)振動(dòng)與譜分析概論》一書，本文不再介紹。

對于寬帶隨機(jī)過程，以上述窄帶分析法為基礎(chǔ)進(jìn)行拓展，也可得出計(jì)算疲勞損傷的近似表達(dá)式。常見的寬帶疲勞算法有DirliK算法、Wirsching-Light算法等，其中Dirlik算法的計(jì)算結(jié)果與試驗(yàn)結(jié)果接近，成為基于功率譜密度計(jì)算疲勞失效的首選算法，已被大多數(shù)商用疲勞分析軟件采用。
關(guān)于這部分的詳細(xì)介紹見下一節(jié)：隨機(jī)振動(dòng)疲勞分析——振動(dòng)疲勞的壽命計(jì)算。