魏效玲,顧林賓
(河北工程大學(xué)機(jī)械與裝備工程學(xué)院,河北 邯鄲 056038)
曲軸是柴油機(jī)核心機(jī)件[1]。它承受周期變化的氣體壓力、連桿機(jī)構(gòu)產(chǎn)生往復(fù)和旋轉(zhuǎn)慣性力,這些力的作用引起曲軸上產(chǎn)生交變載荷,從而使曲軸發(fā)生疲勞破壞,這種破壞惡劣時會造成曲軸突然斷裂[2]。曲軸失效的大部分原因是圓角應(yīng)力集中引起疲勞裂紋、長時間工作導(dǎo)致的復(fù)合疲勞破壞等。在機(jī)器運(yùn)行狀態(tài)下,有著接近80%的破壞是彎曲疲勞產(chǎn)生的。
國內(nèi)外在零部件再制造的相關(guān)領(lǐng)域進(jìn)行許多的研究,文獻(xiàn)[3]提出的線性連續(xù)損傷疲勞模型。文獻(xiàn)[4]提出了基于層次分析法的評價模型,并建立再制造指標(biāo)體系和定量方法求出再制造相對加工指標(biāo)。但針對柴油機(jī)曲軸的剩余壽命預(yù)測和可再制造性評估的內(nèi)容尤其不足。對廢舊的零部件的再制造可行性分析,對節(jié)能減材具有重要作用。研究和預(yù)測曲軸的疲勞壽命,分析疲勞壽命與再制造要求之間的定性關(guān)系,對曲軸再制造的可行性分析具有重要作用。
有限元方法是模擬曲軸在工作下應(yīng)力狀態(tài)的常用方法,也是構(gòu)建曲軸疲勞壽命預(yù)測模型的重要途徑。以494柴油機(jī)曲軸為具體研究對象,柴油機(jī)主要性能參數(shù)為:額定轉(zhuǎn)速為3200 r/min,標(biāo)定功率 95.5 kW,最大扭矩350 N.m,最大爆發(fā)壓力15 MPa;柴油機(jī)曲軸的主要結(jié)構(gòu)參數(shù)為:主軸頸直徑和長度分別為60 mm和26 mm,曲柄銷直徑和長度分別為46 mm和22 mm,圓角半徑為5 mm,曲柄厚度為18 mm,曲軸的三維造型如圖1所示,曲軸的性能參數(shù)如表1所示。

圖1 曲軸三維實(shí)體模型

表1 材料機(jī)械性能及相關(guān)參數(shù)
應(yīng)用三維建模軟件對曲軸進(jìn)行實(shí)體建模,將根據(jù)曲軸結(jié)構(gòu)參數(shù)建立的三維實(shí)體模型轉(zhuǎn)化成IGES格式并導(dǎo)入ANSYS中。在ANSYS中的Material Modeis功能選項(xiàng)對該模型進(jìn)行材料曲軸模型的參數(shù)設(shè)置,對其進(jìn)行合適的體切割、面選定,然后施加盡量逼近實(shí)際情況的邊界條件。選擇ANSYS分析中單元類型8節(jié)點(diǎn)六面體單元,選用網(wǎng)格大小為3 mm對曲軸進(jìn)行網(wǎng)格劃分。網(wǎng)格劃分結(jié)果為:對主軸頸和過渡圓角等應(yīng)力集中區(qū)域進(jìn)行網(wǎng)格細(xì)化,共劃分568 634個節(jié)點(diǎn),172 973個單元;利用有限元分析軟件ANSYS中Mesh功,得到曲軸的網(wǎng)格化模型圖如圖2所示,其工作狀態(tài)下的應(yīng)力云圖如圖3所示。

圖2 曲軸網(wǎng)格劃分圖

圖3 曲軸變形云圖
由ANSYS有限元軟件分析得到該曲軸的最大應(yīng)力值為11.67 MPa,處在曲柄到軸頸的過渡圓角處。
大量研究表明裂紋擴(kuò)展率和應(yīng)力兩者之間存在正相關(guān)關(guān)系[5]。材料的裂紋擴(kuò)展速率如圖4所示。

圖4 不同應(yīng)力比下裂紋擴(kuò)展速率圖
文獻(xiàn)[6]可知,曲軸裂紋尺寸與平均應(yīng)力之間的表達(dá)式為:
σm=a·Sc+c
(1)
式中,a和c是待求的常數(shù),Sc為疲勞缺陷,單位為μm。
采用軟件MATLAB中的CF工具對模型參數(shù)進(jìn)行估計,取置信度為0.95,可以得到應(yīng)力與裂紋的定量表達(dá)式為:
σm=5.708×105Sc+679.3
(2)
根據(jù)上式可得出疲勞裂紋尺寸與平均應(yīng)力之間的關(guān)系并估算不同裂紋尺寸曲軸在工作狀態(tài)下的應(yīng)力,進(jìn)一步研究曲軸的疲勞壽命預(yù)測。
裂紋閉合效應(yīng)可以阻礙裂紋擴(kuò)大,進(jìn)而達(dá)到延長其剩余壽命的作用[7]。CHABOCHE在此基礎(chǔ)上,提出一種新的非線性連續(xù)疲勞及損傷模型:
dD=f(σmax,σm,D)dN
(3)
式中,D、N、σm分別為損傷變量,載荷作用次數(shù)和平均應(yīng)力。
CHABOCH提出式(4)來表示裂紋與疲勞載荷作用次數(shù)之間的關(guān)系[8],表達(dá)式如下:
(4)
其中,β,α為試驗(yàn)常數(shù),σa為屈服疲勞極限,M0為平均力矩,單位為(N·m)。
該模型可作為疲勞壽命評估的基本模型,仍需根據(jù)理論分析與實(shí)際工況的差異進(jìn)行修正,使之用于曲軸的疲勞壽命模型。
在經(jīng)典疲勞極限理論的基礎(chǔ)上,結(jié)合曲軸的實(shí)際工況,ELBER提出有效應(yīng)力強(qiáng)度因子范圍的概念[9]:
kΔeff=kmax-kop
(5)
其中,kΔeff,Kmax,Kop分別為有效應(yīng)力強(qiáng)度因子,最大應(yīng)力強(qiáng)度因子和張開應(yīng)力強(qiáng)度因子。
各種材料的裂紋擴(kuò)展有效因子表達(dá)式如表2所示。其中球墨鑄鐵的裂紋擴(kuò)展有效因子為:

ξ=0.75+0.3R+0.25R2(-1≤R≤1) (6)
考慮裂紋閉合效應(yīng)影響后,對式(4)修正,其表達(dá)式為:
dD=f(σmax,σm,ξ,D)dN
(7)
式中,Dξ為考慮裂紋閉合效應(yīng)的損傷,對式(4)數(shù)次運(yùn)算,得到裂紋閉合效應(yīng)下的非線性疲勞累積損傷表達(dá)式:
(8)
通過比較,這兩種疲勞壽命模型之間存在一定的誤差,誤差在可接受的公差范圍內(nèi),應(yīng)用該模型確定曲軸的再制造臨界閾值和用于曲軸再制造性可行性評價指標(biāo)。
在臨界狀態(tài)下,再制造的最低要求是曲軸的疲勞壽命等于預(yù)期使用壽命[10]。因此,臨界狀態(tài)下的疲勞缺陷尺寸是保證疲勞壽命滿足再制造最低要求的最大值,該疲勞缺陷定義為曲軸再制造臨界閾值A(chǔ)c:曲軸再制造臨界閾值建模為:
(9)

應(yīng)用再制造臨界閾值定性地表示疲勞性能或缺陷狀態(tài)是否滿足再制造要求。
為定量評價曲軸的再制造性能,豐富現(xiàn)有的再制造性能評價體系,提出疲勞缺陷和疲勞壽命指數(shù)[11],得出再制造性評價方法。如圖5所示。

圖5 曲軸的再制造性評價方法圖
基于疲勞壽命的再制造性評價方法的過程可以概括為:以曲軸的疲勞壽命為基礎(chǔ),根據(jù)理論壽命與預(yù)期壽命的關(guān)系確定再制造臨界閾值,定量表達(dá)再制造的基本要求能否滿足。通過對兩個可再制造性評價指標(biāo)進(jìn)行評價,可得曲軸的可再制造性評價函數(shù)。
再制造參數(shù):

(10)
疲勞缺陷指數(shù):

(11)
疲勞壽命指數(shù):

(12)
再制造性評價函數(shù)為:
RD=f(RC,μAC,μN(yùn))=RC·μAC·μN(yùn)
(13)
①當(dāng)疲勞缺陷不滿足再制造要求(Sc>Ac)時,再制造性評價函數(shù)RD=0;
②當(dāng)疲勞缺陷滿足再制造要求(Sc ③當(dāng)曲軸疲勞壽命=預(yù)期壽命時,再制造性評價函數(shù)最小值RD0=1; 以前面提到的曲軸為例,其機(jī)械性能及其參數(shù)見表1;再制造要求曲軸的預(yù)期壽命的循環(huán)次數(shù)設(shè)定為Ne=1×109,利用以上方法計算,得出曲軸實(shí)際疲勞缺陷為Sc=20.44 μm;根據(jù)曲軸再制造臨界閾值的定義和建立曲軸疲勞壽命模型,以預(yù)期壽命Ne為基礎(chǔ),計算出曲軸再制造臨界閾值為Ac=32.9 μm。曲軸的實(shí)際疲勞缺陷Sc=20 μm,此時Sc 通過對曲軸的疲勞裂紋模型進(jìn)行有限元分析,模擬工作中曲軸的應(yīng)力狀態(tài)。由試驗(yàn)結(jié)果,提出平均應(yīng)力σm與疲勞裂紋Sc之間的定量關(guān)系式,對其進(jìn)行有限元分析,建立疲勞壽命預(yù)測模型。其次,提出再制造臨界閾值A(chǔ)c和再制造性能參數(shù)RC的定義,并建立再制造性能評價參數(shù)的定量表達(dá)式。最后,提出曲軸疲勞壽命的再制造性評價方法,建立基于疲勞壽命的曲軸再制造性評價函數(shù),為曲軸的再制造可行性提供依據(jù),為再制造零部件提供科學(xué)評價。3.3 曲軸再制造可行性分析
4 結(jié)論