羅龍君 郭悅新

摘?要：曲軸作為發(fā)動(dòng)機(jī)的重要運(yùn)動(dòng)部件，其結(jié)構(gòu)性能直接影響著發(fā)動(dòng)機(jī)的可靠性和使用壽命。以Pro/E建立的直列四缸發(fā)動(dòng)機(jī)曲軸為基礎(chǔ)，利用ANSYS有限元法對(duì)曲軸的靜力和模態(tài)分析。通過改變軸頸圓角半徑大小、曲柄銷長(zhǎng)度、曲柄臂厚度、軸頸重疊度的結(jié)構(gòu)參數(shù)，確定了各參數(shù)對(duì)軸頸圓角最大應(yīng)力的影響。模態(tài)分析結(jié)果表明，曲軸在低階次頻率下以彎曲為主，最大振型位移出現(xiàn)在曲軸兩端軸頸處，在高階次頻率下以扭轉(zhuǎn)為主。分析結(jié)果為曲軸的優(yōu)化設(shè)計(jì)和動(dòng)力學(xué)分析提供了指導(dǎo)。

關(guān)鍵詞：曲軸；有限元；靜力分析；模態(tài)分析

中圖分類號(hào)：TK403；U464133.3

曲軸是發(fā)動(dòng)機(jī)的重要組成部分，尺寸的大小直接影響發(fā)動(dòng)機(jī)的整體尺寸、質(zhì)量、可靠性和使用壽命［1］。曲軸不僅承受著氣缸氣體的作用力、活塞連桿組往復(fù)運(yùn)動(dòng)慣性力，還承受著扭轉(zhuǎn)振動(dòng)與彎曲振動(dòng)。分析曲軸靜態(tài)受力及運(yùn)動(dòng)情況，可以有效避免曲軸出現(xiàn)裂紋、變形、斷裂等情況。

針對(duì)曲軸受力計(jì)算，傳統(tǒng)方法是簡(jiǎn)化為簡(jiǎn)支梁或連續(xù)梁，但由于其結(jié)構(gòu)復(fù)雜，無法完全采用解析法求解。目前，國內(nèi)外對(duì)曲軸模型多采用有限元分析方法，但模型簡(jiǎn)化存在差異，載荷施加也不盡相同［27］。付澤民［5］選取1/4結(jié)構(gòu)模型，分析每一連桿軸頸的受力情況，但缺少主軸頸的相應(yīng)受力。

本文以Pro/E建立的單拐曲軸和整體曲軸的三維實(shí)體模型為基礎(chǔ)，利用ANSYS軟件，分析曲軸的軸頸圓角半徑、曲軸銷長(zhǎng)度、曲柄臂厚度等參數(shù)對(duì)軸頸圓角應(yīng)力影響。同時(shí)，對(duì)曲軸的振動(dòng)特性進(jìn)行研究和分析，得到其自振頻率，從而為其設(shè)計(jì)提供了理論基礎(chǔ)。

1?曲軸模型建立及網(wǎng)格劃分

以直列四缸發(fā)動(dòng)機(jī)曲軸為例，考慮曲軸結(jié)構(gòu)復(fù)雜性，為方便分析計(jì)算，忽略曲軸模型上的油孔、倒角以及軸頸處的過渡圓角，采用Pro/E建模的方法（圖1），對(duì)曲軸模型進(jìn)行了分析，并列出了相應(yīng)的模型材料參數(shù)（表1）。

由于曲軸受力具周期性和對(duì)稱性的特點(diǎn)，同時(shí)參考文獻(xiàn)［3］表明運(yùn)用單拐曲軸有限元模型計(jì)算得出的應(yīng)力值大于整體曲軸模型計(jì)算值，有利于保障曲軸設(shè)計(jì)的可靠性。所以在進(jìn)行靜力分析時(shí)使用單拐曲軸模型、三維模型和網(wǎng)格劃分結(jié)果如圖2所示。其中，網(wǎng)格單元選擇十節(jié)點(diǎn)四面體實(shí)體結(jié)構(gòu)單元Tet?10Node92，采用智能網(wǎng)格劃分。該有限元模型有39179個(gè)節(jié)點(diǎn)、37074個(gè)單元。

2?基于ANSYS的曲軸有限元分析

2.1?載荷與約束條件

曲軸承載來自連桿的作用力，把它轉(zhuǎn)化成扭矩，然后通過曲軸輸出，進(jìn)而帶動(dòng)其他零件工作。由于曲軸的轉(zhuǎn)動(dòng)離心力、周期變化的慣性力、往復(fù)慣性力的影響，使得曲軸在彎曲扭轉(zhuǎn)載荷的作用下受力。

在不考慮油孔壓力峰突然變化的情況下，利用有限寬軸頸油膜的應(yīng)力分布，沿軸頸縱向呈二次拋物線分布，沿軸頸周向呈余弦分布，由于曲軸受力具有周期性和對(duì)稱性，只分析連桿軸頸受到壓力時(shí)的軸頸圓角應(yīng)力大小，為方便載荷施加，將載荷簡(jiǎn)化為沿連桿軸頸上半面的均布載荷，其大小為曲軸最大爆發(fā)力50MPa。

轉(zhuǎn)速慣性載荷可以由曲軸額定轉(zhuǎn)速2500r/min轉(zhuǎn)換為角速度261.7rad/s直接加在曲軸上。另外，再添加重力加速度。

就約束而言，在實(shí)際工作中，曲軸主要受主、縱推力軸承的制約。工作時(shí)，軸頸和滑動(dòng)軸承間的間隙依賴于它們間的油膜來潤(rùn)滑。在負(fù)載的作用下，軸頸上、下軸瓦之間會(huì)出現(xiàn)彎曲變形，而軸頸縱向推力軸承能有效地阻止軸向竄動(dòng)，進(jìn)而確?；钊麢C(jī)構(gòu)的工作狀態(tài)。針對(duì)單拐曲軸的2個(gè)主軸頸的滑動(dòng)軸承支撐，施加主軸頸的表面徑向?qū)ΨQ約束，同時(shí)對(duì)曲軸兩個(gè)主軸頸施加軸向位移約束。

2.2?后處理

由表1曲軸結(jié)構(gòu)參數(shù)得到單拐等效應(yīng)力分布，如圖3所示，位移分布如圖4所示，曲軸軸頸過渡圓角處是應(yīng)力最大部位。

2.3?不同結(jié)構(gòu)參數(shù)的應(yīng)力分析

分別改變曲軸的軸頸圓角半徑、曲柄銷長(zhǎng)度、曲柄臂厚度以及軸頸重疊度后，分析其軸頸最大圓角應(yīng)力的變化，得出這些參數(shù)對(duì)曲軸應(yīng)力集中的影響，如圖5—圖7所示。

上述結(jié)果表明，減小軸頸最大圓角應(yīng)力，可以采取的措施有增大軸頸圓角半徑、減小曲柄銷長(zhǎng)度、增大曲柄臂厚度、增大曲柄重疊度，為曲軸尺寸優(yōu)化提供了依據(jù)。當(dāng)然，上述四個(gè)結(jié)構(gòu)參數(shù)不是孤立存在，比如由于曲軸頂端與連桿軸頸距離有限，連桿軸頸處圓角的增大又受到了很大的限制。

3?整體曲軸的模態(tài)分析

曲軸在干擾力或力矩的作用下會(huì)按激勵(lì)的頻率進(jìn)行受迫振動(dòng)，當(dāng)激勵(lì)的頻率與其固有頻率相同時(shí)，會(huì)產(chǎn)生共振現(xiàn)象，嚴(yán)重時(shí)會(huì)造成曲軸斷裂。因此，通過對(duì)曲軸尤其是轉(zhuǎn)速較高的曲軸進(jìn)行自由模態(tài)分析時(shí)，能夠較準(zhǔn)確地判斷出曲軸的振動(dòng)狀態(tài)，并能迅速地分析曲軸尺寸、平衡塊尺寸與布置、飛輪和曲軸材料等參數(shù)對(duì)曲軸動(dòng)力性能的影響。

曲軸固有振動(dòng)模態(tài)與材料和結(jié)構(gòu)有關(guān)，與外部載荷及約束無關(guān)，因此不施加任何約束。曲軸整體的自由模態(tài)前6階（剛體模態(tài)）都約為零，所以主要分析7—12階模態(tài)。限于篇幅，以7階模態(tài)為例，其振型云圖如圖9所示，其他階數(shù)固有頻率和分析如表2所示。

由以上圖表還可以看出曲軸的7—12階振型中，最低頻率為第7階振型286.365Hz，隨著階次上升，其頻率也相應(yīng)增大，振型從彎曲過渡到扭轉(zhuǎn)。曲軸在低階頻率時(shí)，其主要的振動(dòng)模式是彎曲模態(tài)，最大的振動(dòng)位移發(fā)生在曲軸的兩端軸頸，而在高階頻率時(shí)，主要是扭轉(zhuǎn)模態(tài)。

該曲軸的工作轉(zhuǎn)速范圍為1500～3500r/min，由此可以計(jì)算出其激振頻率范圍為50～116.7Hz，從模態(tài)分析可以看出，工作轉(zhuǎn)速對(duì)應(yīng)的頻率比7階的最小頻率286.365低，所以曲柄發(fā)生共振的概率很低。

4?結(jié)論

文中對(duì)直列四缸發(fā)動(dòng)機(jī)曲軸采用Pro/E建立三維模型，利用ANSYS軟件對(duì)單拐曲軸的靜力學(xué)以及整體曲軸的模態(tài)進(jìn)行了有限元分析，獲得如下結(jié)論：

（1）曲軸的主要失效部位為軸頸過渡圓角處，可以采取增大軸頸圓角半徑、曲柄臂厚度、曲柄重疊度或減小曲柄銷長(zhǎng)度的方式來避免強(qiáng)度失效。

（2）曲軸振動(dòng)的主要形式是彎曲和扭轉(zhuǎn)，隨著頻率的增高，曲軸兩端的承載部件會(huì)受到很大的交變載荷，因此必須對(duì)其兩端軸承的剛度和強(qiáng)度進(jìn)行加強(qiáng)并校核。

（3）曲柄臂與主軸頸、曲柄臂與連桿軸頸的交界部位是曲軸振動(dòng)的危險(xiǎn)區(qū)，所以曲軸的設(shè)計(jì)要充分考慮曲柄臂的設(shè)計(jì)參數(shù)以及曲柄臂和曲拐相連處的圓角尺寸。

參考文獻(xiàn)：

［1］楊佳敏，李瑞雪.基于ANSYS的汽油機(jī)曲軸結(jié)構(gòu)分析［J］.農(nóng)機(jī)使用與維修，2019（10）：1417.

［2］余佳奎，李舜酩，李想，等.基于ANSYS的發(fā)動(dòng)機(jī)曲軸有限元靜力與模態(tài)分析［J］.河南科技，2020（23）：3641.

［3］吳軍良，趙清，徐毅煜，等.某柴油機(jī)曲軸三維有限元分析與優(yōu)化［J］.南通航運(yùn)職業(yè)技術(shù)學(xué)院學(xué)報(bào)，2020，19（02）：3944.

［4］施佳裕，王忠，殷文元，等.基于ANSYS的曲軸應(yīng)力及變形敏感度分析［J］.機(jī)械設(shè)計(jì)與制造，2020（06）：231234.

［5］付澤民，李延平，常勇，等.ANSYS環(huán)境中柴油機(jī)曲軸靜動(dòng)特性的有限元分析［J］.柴油機(jī)，2006（01）：3438.

［6］徐兆華，崔志琴，張騰.基于ANSYS的6300柴油機(jī)曲軸的模態(tài)分析［J］.煤礦機(jī)械，2012，33（2）：102103.

［7］付貴，郭湘川.基于ANSYS?Workbench的某活塞發(fā)動(dòng)機(jī)曲軸有限元模態(tài)分析［J］.科技創(chuàng)新與應(yīng)用，2018（06）：2021.

基金項(xiàng)目：2022年華中科技大學(xué)實(shí)驗(yàn)技術(shù)項(xiàng)目《基于工程實(shí)踐的勞動(dòng)教育實(shí)驗(yàn)教學(xué)開發(fā)》

作者簡(jiǎn)介：羅龍君（1988—?），男，湖北荊門人，碩士，工程師，研究方向：工程實(shí)踐創(chuàng)新教育。

*通訊作者：郭悅新（2001—?），女，河北唐山人，本科，助理工程師，研究方向：機(jī)械設(shè)計(jì)制造及自動(dòng)化。