秦 瑤,縱文斌,張偉偉,屈宗長(zhǎng),馬曉紅
(1.西安交通大學(xué)能源與動(dòng)力工程學(xué)院,陜西西安710049;2.吉林油田勘察設(shè)計(jì)院,吉林 松原138000)
基于雨流計(jì)數(shù)法的壓縮機(jī)曲軸疲勞壽命分析
秦 瑤1,縱文斌1,張偉偉1,屈宗長(zhǎng)1,馬曉紅2
(1.西安交通大學(xué)能源與動(dòng)力工程學(xué)院,陜西西安710049;2.吉林油田勘察設(shè)計(jì)院,吉林 松原138000)
以4M12-2.5/250型往復(fù)壓縮機(jī)的曲軸為研究對(duì)象,研究其疲勞失效的過(guò)程,估算其疲勞壽命。對(duì)曲軸進(jìn)行多柔體動(dòng)力學(xué)仿真后得到了曲軸的動(dòng)態(tài)載荷和其危險(xiǎn)節(jié)點(diǎn)的應(yīng)力-時(shí)間歷程。采用雨流計(jì)數(shù)法對(duì)節(jié)點(diǎn)的應(yīng)力-時(shí)間歷程進(jìn)行統(tǒng)計(jì)處理,將其處理成一系列的全循環(huán)或半循環(huán)載荷,編制程序載荷譜,結(jié)合疲勞累積損傷模型和材料的S-N曲線(xiàn),計(jì)算其危險(xiǎn)節(jié)點(diǎn)的疲勞壽命,為曲軸的結(jié)構(gòu)設(shè)計(jì)與優(yōu)化提供一定的理論依據(jù),對(duì)提高曲軸的可靠性有一定的意義。
雨流計(jì)數(shù)法;載荷譜;曲軸疲勞失效
曲軸作為壓縮機(jī)中主要運(yùn)動(dòng)部件之一,其性能好壞在很大程度上決定了壓縮機(jī)的安全性、可靠性和使用壽命。由于曲軸結(jié)構(gòu)比較復(fù)雜,主軸線(xiàn)方向的界面形狀變化比較大,故應(yīng)力分布在曲軸整體結(jié)構(gòu)上具有不均勻性,不同位置的應(yīng)力相差比較大,且存在著嚴(yán)重的應(yīng)力集中。在外界載荷周期性的作用下,曲軸容易產(chǎn)生疲勞破壞。而載荷是隨機(jī)變化的,因此要合理地預(yù)測(cè)零部件的疲勞壽命,必須對(duì)載荷時(shí)間歷程或局部應(yīng)力、應(yīng)變歷程進(jìn)行統(tǒng)計(jì)處理,從而進(jìn)行損傷計(jì)算,完成壽命估算[1]。
目前,常用計(jì)數(shù)法對(duì)于隨機(jī)載荷進(jìn)行統(tǒng)計(jì)處理,包括單參數(shù)計(jì)數(shù)法,范圍對(duì)法,稱(chēng)對(duì)法及雨流計(jì)數(shù)法等。相比于其它方法,雨流計(jì)數(shù)法計(jì)數(shù)過(guò)程與材料真實(shí)應(yīng)力-應(yīng)變特性相符,有堅(jiān)實(shí)的理論基礎(chǔ),并能準(zhǔn)確地計(jì)算結(jié)果。近年來(lái),雨流計(jì)數(shù)法與現(xiàn)代計(jì)算機(jī)技術(shù)相結(jié)合,利用其計(jì)算機(jī)程序快速準(zhǔn)確計(jì)算的特點(diǎn),實(shí)現(xiàn)了其計(jì)算過(guò)程的程序化,因此得到了極大的推廣[2]。本文采用雨流計(jì)數(shù)法對(duì)節(jié)點(diǎn)的應(yīng)力-時(shí)間歷程進(jìn)行統(tǒng)計(jì)處理,編制程序載荷譜,結(jié)合疲勞累積損傷模型和材料的S-N曲線(xiàn),計(jì)算其危險(xiǎn)節(jié)點(diǎn)的疲勞壽命。
雨流計(jì)數(shù)法又為“塔頂法”,是在1968年Mat-sishi和Endo等人在研究材料應(yīng)力-應(yīng)變時(shí)提出來(lái)的,該方法認(rèn)為塑性的存在是疲勞損傷的必要條件。取一段應(yīng)力-時(shí)間歷程曲線(xiàn)時(shí)間為橫坐標(biāo),應(yīng)力為縱坐標(biāo),將該曲線(xiàn)順時(shí)針旋轉(zhuǎn)90°,數(shù)據(jù)點(diǎn)之間的連線(xiàn)就如一些列不同層的屋面,而數(shù)據(jù)記錄點(diǎn)就如一些列的雨點(diǎn),統(tǒng)計(jì)的過(guò)程為:雨點(diǎn)依次順著不同層的屋面內(nèi)側(cè)向下流,直到最底層。其計(jì)算規(guī)則為:
(1)找出應(yīng)力-時(shí)間歷程中的最高峰和最低谷,比較兩者絕對(duì)值的大小,在絕對(duì)值大的點(diǎn)打斷應(yīng)力-時(shí)間歷程,重新對(duì)接,以絕對(duì)值大的點(diǎn)為起點(diǎn),得到新的應(yīng)力-時(shí)間歷程;
(2)雨流按順序依次從應(yīng)力-時(shí)間歷程中的峰值(谷值)的內(nèi)側(cè)向下流,當(dāng)遇到比起始點(diǎn)更大的峰值(更低的谷值)時(shí),便停止流動(dòng);
(3)當(dāng)雨流遇到上面流下的雨流時(shí),則停止流動(dòng);
(4)根據(jù)上面3條規(guī)則,統(tǒng)計(jì)出所有的全循環(huán)或者半循環(huán),記錄并計(jì)算的各個(gè)循環(huán)的幅值和均值。雨流計(jì)數(shù)法示意圖如圖1所示。
將雨流計(jì)數(shù)法的計(jì)數(shù)原理轉(zhuǎn)化為計(jì)算機(jī)程序語(yǔ)言,具體的實(shí)施步驟為:
(1)數(shù)據(jù)壓縮。數(shù)據(jù)壓縮分為相鄰等值數(shù)據(jù)點(diǎn)的剔除和局部峰谷值的提取。為了更真實(shí)的表述某節(jié)點(diǎn)的應(yīng)力-時(shí)間歷程,Ansys施加邊界條件時(shí),盡可能的設(shè)置多的載荷步和子步,本文Ansys分析過(guò)程中的載荷步為2858步,但是過(guò)多的載荷步會(huì)導(dǎo)致運(yùn)用雨流計(jì)數(shù)法處理時(shí),產(chǎn)生很多的無(wú)效數(shù)據(jù)(如圖2所示),在一段實(shí)際的應(yīng)力-時(shí)間里程中,相鄰數(shù)據(jù)點(diǎn)等值或者很多的數(shù)據(jù)點(diǎn)并非局部的峰值、谷值點(diǎn),例如1、2點(diǎn)等。因此,要對(duì)數(shù)據(jù)點(diǎn)進(jìn)行數(shù)據(jù)壓縮。
首先,剔除相鄰等值數(shù)據(jù)點(diǎn)。設(shè)Ansys中輸入的數(shù)組為E(i),i為數(shù)組中元素的編號(hào),若將數(shù)組導(dǎo)入程序過(guò)程中,出現(xiàn)E(i)=E(i+1)則舍去E(i+1),繼續(xù)判斷下一個(gè)數(shù)據(jù)點(diǎn)是否輸入,若E(i)≠E(i+1),則繼續(xù)輸入下一個(gè)數(shù)據(jù)點(diǎn)。其次,峰谷值的提取,對(duì)于剔除相鄰等值數(shù)據(jù)點(diǎn)之后得到新的數(shù)組,并非所有的數(shù)據(jù)點(diǎn)為峰值點(diǎn)或者谷值點(diǎn),對(duì)于峰谷點(diǎn)的判斷為依次相連3個(gè)數(shù)據(jù)點(diǎn)之間的差值的乘積,如果為負(fù)值,則中間點(diǎn)為峰谷點(diǎn);若為正值,則中間點(diǎn)非峰谷值,繼續(xù)進(jìn)行下一組數(shù)據(jù)點(diǎn)的判斷,提取峰谷值組成新的數(shù)組,作為雨流計(jì)數(shù)法處理的數(shù)組。
(2)循環(huán)數(shù)的提取。通常分為一次雨流計(jì)數(shù)法、打斷峰值(谷值)點(diǎn)重新對(duì)接組成新的應(yīng)力-時(shí)間歷程、二次雨流計(jì)數(shù)法。對(duì)于數(shù)組進(jìn)行數(shù)據(jù)壓縮過(guò)之后,對(duì)數(shù)據(jù)點(diǎn)實(shí)施雨流計(jì)數(shù)法提取循環(huán),雨流計(jì)數(shù)法規(guī)則在計(jì)算機(jī)程序中實(shí)現(xiàn)可采用“四點(diǎn)法”,即符合計(jì)數(shù)規(guī)則又易于在程序中實(shí)現(xiàn)。如圖3所示:在應(yīng)力-時(shí)間歷程中能提取循環(huán)的只有這兩種情況,其中a、b、c分別為其兩端點(diǎn)值差的絕對(duì)值,若b≤a且b≤c,則提取一個(gè)循環(huán)圖中三角形部分,并且剔除E(i)和E(i+1),保留E(i-1)和E(i+2),E(i-1)下標(biāo)值不變,而E(i+2)的下標(biāo)號(hào)從i+2開(kāi)始整體提前2個(gè)變?yōu)镋(i),即E(i)=E(i+1),否則保留E(i-1),并從下標(biāo)i開(kāi)始進(jìn)行下一組4個(gè)數(shù)據(jù)點(diǎn)的“四點(diǎn)法”的比較。
第一次雨流計(jì)數(shù)法提取循環(huán)數(shù)并記錄之后,剩下的數(shù)據(jù)點(diǎn)構(gòu)成的波形為發(fā)散-收斂型,這時(shí)利用雨流計(jì)數(shù)法無(wú)法在剩下的數(shù)據(jù)點(diǎn)中提取完整的循環(huán),故需將剩下的波形進(jìn)行重新的對(duì)接。在最大值(最小值)處將波打斷為左右兩段,將左段波的起點(diǎn)對(duì)接在右段波的末點(diǎn),組成新的波形。若右段波的末點(diǎn)的數(shù)值與左段波的起點(diǎn)的數(shù)值不等,則在左段波的起點(diǎn)之前補(bǔ)充一個(gè)與右段波的末點(diǎn)數(shù)值相等的點(diǎn)作為新的起點(diǎn)與右段波對(duì)接。對(duì)接之后的波形為收斂發(fā)散型,對(duì)其再次使用雨流計(jì)數(shù)法進(jìn)行處理,再次提取循環(huán)數(shù)直到剩下最后3個(gè)數(shù)據(jù)點(diǎn)為止,最后3個(gè)數(shù)據(jù)點(diǎn)組成一個(gè)循環(huán)。至此,整個(gè)雨流計(jì)數(shù)法程序完成。
程序譜又稱(chēng)塊譜即按照一定的程序施加載荷,在疲勞過(guò)程中,若塊譜中的塊較少時(shí),塊譜的不同編排順序?qū)?gòu)件的疲勞壽命影響較大,但若塊較多時(shí),編排順序?qū)?gòu)件的疲勞壽命幾乎沒(méi)有影響。故將雨流計(jì)數(shù)法得到的結(jié)果編制為程序載荷譜,具體的編制步驟為:
(1)對(duì)雨流計(jì)數(shù)法得到的結(jié)果進(jìn)行統(tǒng)計(jì)得到應(yīng)力幅的直方圖;
(2)判斷載荷譜所服從的分布以及分布參數(shù);
(3)對(duì)載荷譜進(jìn)行擴(kuò)充到106次,若在106中有出現(xiàn)一次的載荷,則取之為最大載荷幅;
(4)對(duì)載荷幅值分等級(jí)編制。
圖4為文獻(xiàn)[3]中對(duì)曲軸進(jìn)行多柔體動(dòng)力仿真及瞬態(tài)響應(yīng)得到的,曲軸在曲柄銷(xiāo)合力最大的時(shí)刻的應(yīng)力分布云圖。
由圖4可觀察到,一級(jí)、二級(jí)、三級(jí)、四級(jí)曲柄銷(xiāo)在所受全力最大時(shí)刻應(yīng)力最大的節(jié)點(diǎn)(危險(xiǎn)點(diǎn))分別為45031、45141、143050、141589。通過(guò)時(shí)間歷程處理后得到的各個(gè)危險(xiǎn)節(jié)點(diǎn)等效應(yīng)力隨時(shí)間變化的曲線(xiàn),如圖5所示。
對(duì)節(jié)點(diǎn)141589在一個(gè)模擬周期的應(yīng)力-時(shí)間歷程進(jìn)行雨流計(jì)數(shù)法處理,將處理得到的結(jié)果80MPa分為一個(gè)區(qū)間,統(tǒng)計(jì)的結(jié)果如表1所示。
用Matlab中hist函數(shù)得到應(yīng)力幅直方圖如圖6所示。
應(yīng)力幅一般服從威布爾分布,威布爾分布的概率密度函數(shù)和分布函數(shù)分別為
式中 m——形狀參數(shù)
η——尺度參數(shù)
γ——位置參數(shù)
運(yùn)用Matlab中wblplot函數(shù)可判斷應(yīng)力幅密度分布是否為威布爾分布,執(zhí)行函數(shù)wblplot(x)若得到的為一條直線(xiàn)即判斷服從威布爾分布,若為一條曲線(xiàn)則判斷不服從威布爾分布。執(zhí)行函數(shù)wblplot(x)結(jié)果如圖7所示。
獲得的虛線(xiàn)為一條直線(xiàn),證明應(yīng)力幅符合威布爾分布,這時(shí)γ=0,運(yùn)用Matlab中wblfit函數(shù)得到m=1.3478,η=112.5112,將數(shù)據(jù)代入到公式(1)和(2)中得
統(tǒng)計(jì)各個(gè)應(yīng)力水平下的循環(huán)作用次數(shù),并將其擴(kuò)充為1×106次循環(huán)中出現(xiàn)一次最大應(yīng)力的累積頻次范圍。最大的應(yīng)力幅值取1×106次循環(huán)中出現(xiàn)一次的數(shù)值;由概率分布函數(shù)可得到最大幅值的應(yīng)力xmax,即
求解函數(shù)可得到 xamax=γ+13.8161/mη
將其中的參數(shù)代入結(jié)果中可得等到最大幅值為789.4093 MPa。
應(yīng)力幅值等級(jí)一般取8級(jí),各級(jí)應(yīng)力幅值依次取最大應(yīng)力幅值的1倍、0.95倍、0.85倍、0.725倍、0.575倍、0.425倍、0.275倍、0.125倍。故此節(jié)點(diǎn)得到的各級(jí)應(yīng)力幅依次為:789.4093 MPa、749.9388 MPa、670.9979 MPa、572.3218 MPa、453.9104 MPa、335.4990 MPa、217.0876 MPa、98.6762 MPa。
運(yùn)用材料的韌性耗散理論應(yīng)用在疲勞損傷過(guò)程中,推導(dǎo)出了多級(jí)載荷循環(huán)作用下的疲勞累積損傷模型,得到其任意應(yīng)力幅加載下的累積損傷為[3]
式中 i=2,3,…,當(dāng)i=1時(shí),T1=1
σai——第i級(jí)應(yīng)力幅
Nfi——應(yīng)力幅σai對(duì)應(yīng)的疲勞壽命
Ni——第i級(jí)應(yīng)力幅循環(huán)次數(shù)
公式(7)中σf′,b可以根據(jù)文獻(xiàn)[11]提供的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合,另y=σai,x=Nfi,則有y=σf′(2x)b,對(duì)方程兩邊同時(shí)取對(duì)數(shù)得
上式說(shuō)明應(yīng)力幅和應(yīng)力幅對(duì)應(yīng)的疲勞壽命在雙對(duì)數(shù)坐標(biāo)系下呈線(xiàn)性變化關(guān)系,A、b為常數(shù),其值可通過(guò)雙對(duì)數(shù)坐標(biāo)系下的最小二乘法擬合得到。
運(yùn)用Matlab中的ispoly函數(shù)擬合實(shí)驗(yàn)數(shù)據(jù)得到A、b的數(shù)值,并將其取平均值,最后得到, σf′=936.9 MPa,故有,將其表達(dá)式轉(zhuǎn)化為
根據(jù)公式(9)計(jì)算所得的8級(jí)應(yīng)力幅值及對(duì)應(yīng)的循環(huán)周次,應(yīng)用材料韌性耗散的疲勞累積損傷模型和材料的S-N曲線(xiàn),可以得到多級(jí)應(yīng)力水平下的疲勞損傷,計(jì)算結(jié)果如表2所示。
故總的損傷量D為
當(dāng)累積損傷量達(dá)到臨界值時(shí),曲軸發(fā)生失效斷裂,查閱相關(guān)文獻(xiàn)[3],臨界損傷量取1,故循環(huán)次數(shù)N為
N=1/D=1.38420820×106
在此工況下,節(jié)點(diǎn)141589可以達(dá)到的循環(huán)次數(shù)為1.38420820×106次。
按照同樣的處理方法可以得到此工況下其他節(jié)點(diǎn)的疲勞累積損傷量以及循環(huán)數(shù),具體如下:
節(jié)點(diǎn)45031一個(gè)循環(huán)周期內(nèi)的疲勞累積損傷量為4.65891×10-8,可以達(dá)到的循環(huán)次數(shù)為2.14642687×107次;
節(jié)點(diǎn)45141一個(gè)循環(huán)周期內(nèi)的疲勞累積損傷量為7.66292×10-8,可以達(dá)到的循環(huán)次數(shù)為1.30498476×107次;
節(jié)點(diǎn)143050一個(gè)循環(huán)周期內(nèi)的疲勞累積損傷量為5.03973×10-7,可以達(dá)到的循環(huán)次數(shù)為1.98423250×106次。
從計(jì)算4個(gè)節(jié)點(diǎn)的損傷量和循環(huán)次數(shù)可以看出,循環(huán)次數(shù)從高到低的節(jié)點(diǎn)依次為45031、45141、143050、141589,4級(jí)曲柄銷(xiāo)節(jié)點(diǎn)處的損傷量最大,循環(huán)次數(shù)最少,最易發(fā)生疲勞失效,對(duì)于成品的曲軸,在實(shí)際的工作過(guò)程中,時(shí)時(shí)檢測(cè)四級(jí)曲柄銷(xiāo)與曲柄連接處的損傷情況,關(guān)注疲勞裂紋的產(chǎn)生和擴(kuò)展情況,可以有效的預(yù)防和避免曲軸失效突然斷裂的情況發(fā)生,避免事故的發(fā)生以及財(cái)產(chǎn)的損失;對(duì)于在設(shè)計(jì)過(guò)程中的曲軸,4級(jí)曲柄銷(xiāo)與曲柄連接處是曲軸結(jié)構(gòu)設(shè)計(jì)和優(yōu)化的重點(diǎn),例如可以增加過(guò)渡圓角的半徑、加強(qiáng)曲軸表現(xiàn)質(zhì)量系數(shù)等措施,改善曲軸的惡劣的受力情況,提高曲軸的安全性和壽命,提高曲軸的性能。
本文在現(xiàn)有的研究基礎(chǔ)上,采用雨流計(jì)數(shù)法對(duì)曲軸上危險(xiǎn)節(jié)點(diǎn)在一個(gè)周期的應(yīng)力-時(shí)間歷程進(jìn)行處理。將應(yīng)力-時(shí)間歷程經(jīng)雨流計(jì)數(shù)法處理得到的幅、均值二維矩陣進(jìn)行統(tǒng)計(jì)處理,并編制程序載荷譜,結(jié)合材料韌性耗散的累積損傷模型及材料的S-N曲線(xiàn),估算了曲軸上不同的危險(xiǎn)節(jié)點(diǎn)的損傷量以及可循環(huán)次數(shù)。對(duì)比4個(gè)節(jié)點(diǎn)處的計(jì)算結(jié)果,節(jié)點(diǎn)141589的損傷量最大,可循環(huán)次數(shù)最少,最易發(fā)生疲勞破壞,故4級(jí)曲柄銷(xiāo)與曲柄連接處是曲軸結(jié)構(gòu)設(shè)計(jì)和優(yōu)化的重點(diǎn)。
本文的研究成果對(duì)實(shí)際工作中有效的預(yù)防和避免曲軸失效突然斷裂的情況發(fā)生,避免事故的發(fā)生以及財(cái)產(chǎn)的損失有重大意義;對(duì)于在設(shè)計(jì)過(guò)程中的曲軸,4級(jí)曲柄銷(xiāo)與曲柄連接處是曲軸結(jié)構(gòu)設(shè)計(jì)和優(yōu)化的重點(diǎn),例如可以增加過(guò)渡圓角的半徑、加強(qiáng)曲軸表現(xiàn)質(zhì)量系數(shù)等措施,改善曲軸惡劣的受力情況,提高曲軸的安全性和壽命,提高曲軸的性能。
[1] 金德新.改進(jìn)的雨流計(jì)數(shù)法應(yīng)用于隨機(jī)載荷下的壽命預(yù)測(cè)[J].鞍鋼技術(shù),2000,(5).
[2] 賈軍利,崔新維.雨流計(jì)數(shù)法及其載荷譜的編制[J].現(xiàn)代商貿(mào)工業(yè),2008,20(6):326-326.
[3] 秦瑤,縱文斌.4M12型壓縮機(jī)曲軸多柔體動(dòng)力學(xué)仿真分析[J].壓縮機(jī)技術(shù),2015,(5).
[4] 姚衛(wèi)星.結(jié)構(gòu)疲勞壽命分析[M].北京:國(guó)防工業(yè)出版社,2004.
[5] 尚德廣,王德俊.多軸疲勞強(qiáng)度[M].北京:科學(xué)出版社,2007.
[6] 武秀根,鄭百林,楊青,等.疲勞累積損傷理在曲軸疲勞分析中的應(yīng)用[J].同濟(jì)大學(xué)學(xué)報(bào),2008,36(5):655-658.
[7] 鄭葉紅.艦船柴油機(jī)曲軸全壽命周期疲勞主動(dòng)可靠性研究[D].江蘇科技大學(xué),2009.
[8] 董樂(lè)義,羅俊,程禮.雨流計(jì)數(shù)法及其在程序中的具體實(shí)現(xiàn)[J].航空計(jì)測(cè)技術(shù),2004,24(3):38-40.
[9] 胡明敏.基于韌性耗散模型的損傷定量分析方法[J].力學(xué)季刊,2000,21(3):387-391.
[10]張茶花,陳繼傳,段巍.基于雨流計(jì)數(shù)法的風(fēng)力機(jī)葉片疲勞載荷統(tǒng)計(jì)分析[J].科教導(dǎo)刊(上旬刊),2013,02:190-192.
[11]葉篤毅,王德俊.疲勞損傷過(guò)程中45號(hào)鋼剩余力學(xué)性能的變化特征[J].材料研究學(xué)報(bào),1996,10(4):357-361.
The Compressor Crankshaft Fatigue Life Analysis Based on Rain Flow Count Method
QIN Yao1,ZONG Wen-bing1,ZHANG Wei-wei1,QU Zong-chang1,MA Xiao-hong2
(1.School of Energy and Engineering Xi′an Jiaotong University,Xi′an 710049,China;2.Survey and Design Institute of Jilin Oil Field,Songyuan 138000,China)
This paper carried out research works on the fatigue failure of the crankshaft for 4M12-2.5/250 reciprocating compressor. The dynamic load of crankshaft and the actual stress-time after the multi-body dynamics simulation can be got.The rain flow counting method was applied to deal with the stress-time history which could be treated as a series of full cycle and half cycle.Programmed load spectrum combined with fatigue cumulative damage model and the S-N curve of material was applied to calculate the fatigue life of dangerous nodes.It had provided certain theoretical basis for the structural design and the optimization of the crankshaft,and of some significance to improve the reliability of the crankshaft.
rain flow counting method;load spectrum;fatigue failure of the crankshaft
TH457
A
1006-2971(2015)06-0001-07
秦瑤(1990-),女,陜西人,西安交通大學(xué)在讀研究生。
2015-01-20
國(guó)家重大科技專(zhuān)項(xiàng) (2011ZX05054)