郭建華

(黑龍江工程學(xué)院)

基于ANSYS的發(fā)動機曲軸有限元模態(tài)分析

郭建華

(黑龍江工程學(xué)院)

曲軸的強度對于發(fā)動機整體性能有著非常重要的影響.通過Solidworks軟件對發(fā)動機曲軸進行建模，通過ANSYS有限元分析軟件對曲軸的模態(tài)進行分析，根據(jù)曲軸各階模態(tài)的固有頻率及振型，對曲軸的優(yōu)化設(shè)計提供一定的理論分析基礎(chǔ).

發(fā)動機曲軸;ANSYS;模態(tài)分析;有限元分析

0 引言

曲軸作為發(fā)動機內(nèi)部重要組成零部件之一，在發(fā)動機工作過程中，曲軸與連桿相結(jié)合，使得活塞沖壓汽油噴霧，進而汽油燃燒產(chǎn)生的內(nèi)能轉(zhuǎn)化成傳動軸的轉(zhuǎn)動動能，最終實現(xiàn)動力機構(gòu)的運動.其工作環(huán)境比較惡劣，一直處于一種汽油沖擊內(nèi)能沖擊壓力、傳動軸轉(zhuǎn)動慣性力的周期性載荷作用中，極易出現(xiàn)塑性變形、疲勞失效等情況，因此，保證曲軸的強度對于保證發(fā)動機整體性能，提高發(fā)動機的使用壽命有著非常重要的作用.據(jù)分析，當(dāng)發(fā)動機內(nèi)曲軸的激勵振動頻率與某階固定頻率相近或者相等的時候，曲軸極易產(chǎn)生共振現(xiàn)象，進而產(chǎn)生較大的應(yīng)力，發(fā)生較大的應(yīng)變，對發(fā)動機整體性能產(chǎn)生非常不良的影響，因此，對發(fā)動機曲軸進行模態(tài)分析具有非常重要的意義.

1 有限元分析方法及模態(tài)分析理論

1.1 有限元分析方法

有限元分析方法融合了力學(xué)分析方法與數(shù)值分析方法，從兩個角度去分析解決問題，其分析過程一方面從力學(xué)角度來說，有限元分析方法將連續(xù)體結(jié)構(gòu)通過離散化，將非常復(fù)雜的連續(xù)體結(jié)構(gòu)離散成若干個離散的子單元，子單元與子單元之間通過節(jié)點進行連接;另一方面從數(shù)值分析分析方法角度來說，有限元方法將復(fù)雜的偏微分方程離散化，最終形成簡單的代數(shù)方程進行數(shù)值求解，大大減少了工程數(shù)值計算工作量.

對于具體的連續(xù)體結(jié)構(gòu)的有限元分析過程，首先，對復(fù)雜連續(xù)體通過多種離散方法，離散成有限個子單元，相鄰子單元之間通過節(jié)點進行連接，實現(xiàn)整個復(fù)雜連續(xù)體的離散化;其次，對于每一個單元，在子單元內(nèi)部選擇合適的插值函數(shù)，根據(jù)實際問題需要，對單元差值函數(shù)設(shè)置相應(yīng)的邊界條件，根據(jù)連續(xù)體平衡條件，窮列出一系列以節(jié)點位移為未知量的線性方程組，求解可得出不同節(jié)點的位移;最后，根據(jù)彈性力學(xué)中節(jié)點與單元之間的關(guān)系，運用節(jié)點位移方程，轉(zhuǎn)化計算出各個單元的應(yīng)力和應(yīng)變，然后就可以得知整個連續(xù)體的應(yīng)力和應(yīng)變.

相對于傳統(tǒng)的數(shù)值分析方法，有限元分析方法通過離散化，使得工程數(shù)值計算過程更為簡單，大大減少了計算工作量，也使得計算機計算運行時間大幅度減少.此外，有限元分析方法也使得一些復(fù)雜工程問題、傳統(tǒng)數(shù)值分析方法不能求解的問題得以解決，根據(jù)其分析理論，如果整個連續(xù)體的離散化程度足夠高，得到的計算結(jié)果與實際計算結(jié)果的誤差將會越來越小，其具體誤差在工程實際問題誤差范圍之內(nèi)，因此，有限元分析方法在工程實際問題中的應(yīng)用非常廣泛.

1.2 有限元模態(tài)分析理論

發(fā)動機曲軸屬于一種復(fù)雜無限自由度系統(tǒng)，因此可以離散為有限個自由度的多自由振動系統(tǒng).根據(jù)牛頓第二定律計算出該多自由度振動系統(tǒng)的微分方程:

其中，M，C，K分別為質(zhì)量矩陣、阻尼矩陣和剛度矩陣;u(t)為系統(tǒng)的位移向量，f(t)為系統(tǒng)的激振力向量，為系統(tǒng)的加速度向量，為系統(tǒng)的速度向量，u0為系統(tǒng)初始位移向量，為系統(tǒng)的初速度向量.

對于具有n個自由度的振動系統(tǒng)，需要用n階向量來描述物理參數(shù)模型.

將式(1)兩邊通過拉氏變換可得:

令Z(s)=(K+sC+s2M)得:

其中Z(s)稱為阻抗矩陣，H(s)為傳遞函數(shù)矩陣，將jω代替s進行博氏域處理可得:

設(shè)有一點l，則可得其相應(yīng)表達(dá)式如式(7)所示:

其中qr(ω)為階模態(tài)坐標(biāo)，φlr為測點l的r階模態(tài)振動系數(shù)，N個測點的各階振動系數(shù)組成向量，稱為r階模態(tài)向量.各階模態(tài)向量組成系統(tǒng)的模態(tài)矩陣.各階模態(tài)坐標(biāo)則組成模態(tài)坐標(biāo)向量.

由式(8)(9)(10)可得:

將式(11)代入式(4)可得:

在進行模態(tài)分析時，系統(tǒng)的微分方程可以看做一個多自由度無阻尼系統(tǒng)，對于無阻尼系統(tǒng)的自由振動，式(12)可以寫成一下形式:

則可以得到微分方程的特征方程:

其中ω為系統(tǒng)的固有頻率.

通過求解，就可以得到矩陣的特征值和特征向量，即系統(tǒng)模型的固有頻率和振型.

2 有限元曲軸模型的建立

2.1 曲軸的三維模型建模

由于ANSYS建模的復(fù)雜性和局限性，該文應(yīng)用Solidworks三維建模軟件對曲軸進行三維預(yù)建模，為了降低ANSYS數(shù)值計算計算量，提高計算機計算速度和精度，根據(jù)經(jīng)驗，對Solidworks建模的曲軸結(jié)構(gòu)中的圓角、倒角進行刪除，對曲軸結(jié)構(gòu)中的油道細(xì)小特征進行忽略，通過對曲軸結(jié)構(gòu)進行簡化，最終建模成如圖1所示結(jié)構(gòu)，然后保存為x_t文件格式，以備有限元分析導(dǎo)入.

圖1 曲軸的三維建模

2.2 模型的導(dǎo)入與單元類型的選擇

運用ANSYS 14.5軟件導(dǎo)入曲軸的x_t格式建模文件，改變模型實體單元，添加ANSYS單元種類為Solid Tet 10node 187，該種單元為10節(jié)點四面體單元，就是運用10個節(jié)點來定義子單元體(如圖2所示)，每一個節(jié)點定義3個沿xyz方向平移的自由度，這種單元種類支持塑性、超彈性材料，可以進行材料的大變形、大應(yīng)變、蠕變等分析，而且對于通過其他三維制圖軟件創(chuàng)建的三維模型，該種單元類型的二次位移模式能夠?qū)δＰ瓦M行更好的仿真.

圖2 Solid 187單元體

2.3 網(wǎng)格劃分

該文ANSYS仿真所分析的曲軸的材料選擇為QT800－2，經(jīng)過查詢材料的屬性見表1.

表1 QT800－2材料屬性

網(wǎng)格的劃分先使用smart size確定劃分精度為6，采用四面體網(wǎng)格劃分方法(Tetrahedral)，曲軸網(wǎng)格劃分后的模型如圖3所示，網(wǎng)格劃分后的節(jié)點數(shù)為76518，單元個數(shù)為49119，具體模型結(jié)構(gòu)信息如圖4所示:

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

圖4 曲柄網(wǎng)格劃分信息

3 曲軸的有限元模態(tài)分析

對曲軸進行模態(tài)分析是為了避免在發(fā)動機運轉(zhuǎn)時出現(xiàn)共振現(xiàn)象，因此，在對曲軸進行受到約束時的模態(tài)分析，則需要對曲軸的固有頻率進行計算.因此對于曲軸的模態(tài)分析，通常情況下需要進行自由模態(tài)和約束模態(tài)兩種分析方式，自由模態(tài)分析主要是為了分析出曲軸的固定頻率和振型，進而可以對曲軸的自身材料特點進行分析，在忽視約束條件下對曲軸的性能進行研究，進而可以為曲軸的動力學(xué)分析提供力學(xué)分析依據(jù);約束模態(tài)分析則是針對曲軸在發(fā)動機中所受到的邊界條件限制，通過改變ANSYS中約束條件，進而分析出曲軸的不同模態(tài)變化，約束模態(tài)分析方法更能適用于實際工作環(huán)境.

選擇ANSYS分析方式為Modal，在模態(tài)階數(shù)的選擇方面，由于發(fā)動機自身振動頻率階數(shù)較低，因此，選擇10階模態(tài)分析就可以分析得出曲軸的模態(tài)變化.在模態(tài)提取方法選擇方面，使用Block Lanczos模態(tài)提取方法，因為該種模態(tài)提取方法經(jīng)常應(yīng)用于實體單元或殼單元模型，特別是對于節(jié)點數(shù)在100000以下較大型的模型來說，該種模態(tài)提取方法能夠很好地處理其剛體振型，在模態(tài)分析中的應(yīng)用非常廣泛.最終分析結(jié)果見表2.

表2 曲軸前10階模態(tài)分析頻率

從表2曲軸前10階模態(tài)分析頻率中可以看出，曲軸在受到約束條件下的約束頻率比在自由條件下的自由頻率提高了許多.具體來說，在自由狀態(tài)下，前6階模態(tài)的自由頻率均很小，因此可以忽略，對于第7模態(tài)的自由頻率大小與第1，2模態(tài)的約束頻率相近，第9模態(tài)的自由頻率大小與第4模態(tài)的約束頻率大小相近，因此，在這兩段頻率段就極易產(chǎn)生共振現(xiàn)象，因此需要對這兩段頻率段進行重點分析.

由于約束模態(tài)分析結(jié)果更適用于工程實際，而且曲軸的共振對發(fā)動機的性能影響很大，因此，通過發(fā)動機中曲軸的所受到的約束條件，對曲軸模型兩端進行全約束限制，進而對曲軸進行10階模態(tài)分析，最終得出的10階模態(tài)分析模型如圖5所示.

從圖5曲軸前10階模態(tài)振型圖中可以看出，曲軸在第一階、第四階、第六階、第七階模態(tài)下的振型變化比較大，具體來說，曲軸在第一階模態(tài)下，靠近前端的中間位置曲柄臂和平衡塊部位的變形量較大，靠近前端的中間位置曲柄臂沿著y軸方向的彎曲變形量較大;曲軸在第二階模態(tài)下，靠近前端的中間位置的主軸頸部位位移變化量較大，但是整個曲柄的彎曲振動變化量不大;曲軸在第三階模態(tài)下，靠近前端的曲柄臂及主軸頸的位移變形量大，曲柄整體彎曲振動較小;曲柄在第四階模態(tài)下，靠近輸出端的曲柄臂及平衡塊部位變形量較大，曲柄整體的彎曲振動變化量較大，曲柄整體呈現(xiàn)波形扭曲;曲柄在第六階模態(tài)下，曲柄整體的位移變形量較小，但是曲柄中間平衡塊位置的彎曲振動變化量較大，曲柄整體還出現(xiàn)一定程度的扭轉(zhuǎn);曲柄在第七階模態(tài)下，曲柄整體的位移變形量較小，但是曲柄靠近前端部位的彎曲振動變形量較大;曲柄在第十階模態(tài)下，曲柄整體沿Z軸方向呈現(xiàn)較大程度的扭轉(zhuǎn)變形.

4 結(jié)論

(1)該文利用Solidworks軟件進行發(fā)動機曲柄的建模，通過ANSYS分析軟件對曲柄的模態(tài)進行分析，在分析過程中，為了更為系統(tǒng)的分析曲軸的振動，該文通過分析曲柄在自由狀態(tài)下的自由頻率得出曲柄本身的固有頻率，然后通過曲柄的實際工況，對曲柄添加約束條件，分析得出曲柄的約束頻率，通過對比自由頻率與約束頻率，最終分析得出:在受到約束條件下，曲柄的振動頻率較自由狀態(tài)出現(xiàn)較大的增加.

(2)通過分析實際，對約束狀態(tài)下曲柄的振型進行分析，通過對比分析可以發(fā)現(xiàn):在第一階、第四階、第六階、第七階模態(tài)下，曲柄的平衡塊及主軸頸部位的位移變化量較大，曲柄在第十階模態(tài)下的扭轉(zhuǎn)變形量最大，因此可以判斷曲柄的彎曲共振頻率在1000～1200Hz之間，曲柄的扭轉(zhuǎn)共振頻率在1900Hz以上，因此，在曲柄的設(shè)計制造過程中，應(yīng)該著重關(guān)注這兩段頻率段.

(3)曲柄的受力情況時比較復(fù)雜的，在考慮曲柄受力載荷變化的情況下，需要對曲柄的強度和剛度進行進一步優(yōu)化，使其滿足設(shè)計強度要求，而且，為了保證曲軸的平穩(wěn)，還需要對曲軸的平衡性進行設(shè)計和研究.

［1］王啟興，馬星國，尤小梅，等.基于Pro/E和ANSYS的某型柴油機曲軸模態(tài)分析［J］.機械工程與自動化，2016 (3):53－55..

［2］陳然，曹詠弘，孫華東，等.基于ANSYS的內(nèi)燃機曲軸的模態(tài)分析［J］.煤炭工程，2015，47(4):108－109.

［3］劉昌領(lǐng)，陳建義，李清平，等.基于ANSYS的六缸壓縮機曲軸模態(tài)分析及諧響應(yīng)分析［J］.流體機械，2012，40 (8):17－21.

［4］田俊龍.發(fā)動機曲柄連桿機構(gòu)的扭振及疲勞分析［D］.太原理工大學(xué)，2016.

［5］王康，夷紅志，劉有成，等.一種小型二沖程汽油機曲軸箱有限元動態(tài)特性分析［J］.機械研究與應(yīng)用，2016(1): 26－31.

［6］于海號，鐘佩思，宋慧.基于ADAMS和Workbench的曲軸動力學(xué)仿真和有限元分析［J］.科技展望，2016，26 (22):42－46.

［7］米小珍，郭萬劍，王楓，等.基于模態(tài)應(yīng)力恢復(fù)理論的內(nèi)燃機曲軸疲勞壽命評估［J］.制造業(yè)自動化，2016，38 (10):23－26.

［8］張帥，張繼忠，孫學(xué)良，等.基于內(nèi)燃式壓氣機曲軸系多體動力學(xué)仿真與響應(yīng)分析研究［J］.壓縮機技術(shù)，2016 (4):22－29.

Finite－element Modal Analysis of Engine Crankshafe Based on ANSYS Workbench

Guo Jianhua
(Heilongjiang Institute of Technology)

The strength of the crankshaft has very important impact on the overall performance of the engine.In this paper，the engine crankshaft is modeled by Solidworks Software，and the modal analysis of the crankshaft is analyzed by ANSYS FEM software，including natural frequency and vibration modes of crankshaft modal，and the theoretical basis for the optimum design of crankshaft is provided.

Engine crankshaft;ANSYS;Modal analysis;Finite－element analysis

TB115

:A

:1000－5617(2017)01－0034－05

(責(zé)任編輯:季春陽)

2016－12－15