基于Adams 的軸承摩擦力矩動力學(xué)仿真與驗證*

黃佳男 王建航 范零峰 雷良育，2

（1.浙江農(nóng)林大學(xué)；2.浙江兆豐機(jī)電股份有限公司）

輪轂軸承主要用于作為汽車車軸處的質(zhì)量承擔(dān)和為輪轂的轉(zhuǎn)動提供精確引導(dǎo)的零部件，既承受徑向載荷又承受軸向載荷，是汽車載重和轉(zhuǎn)動的重要組成部分，關(guān)乎汽車安全和性能。摩擦力矩作為輪轂軸承單元的1 項重要性能指標(biāo)，直接影響到軸承及軸類零件的使用壽命和傳動效率，在輪轂軸承的加工過程中，需要經(jīng)過粗磨、精磨、粗超精及精超精等大量工序來降低輪轂軸承接觸中的摩擦力，摩擦力矩檢測也是輪轂軸承出廠綜合性檢測的一項重要指標(biāo)。因此，在輪轂軸承的分析優(yōu)化環(huán)節(jié)中，要對輪轂軸承進(jìn)行優(yōu)化，首先就需要對摩擦力矩的組成形式及產(chǎn)生方式進(jìn)行分析?；谝陨蠁栴}，文章提出1 種動力學(xué)仿真方法，通過Adams 動力學(xué)仿真對理論研究中的輪轂軸承摩擦力矩進(jìn)行分析，驗證輪轂軸承結(jié)構(gòu)的合理性，經(jīng)過臺架試驗驗證，達(dá)到預(yù)想的效果。

1 重卡輪轂軸承單元的組成

重卡是重型卡車的簡稱，相比家用小型汽車，重型卡車的工作環(huán)境更加惡劣，承載的力更大。因此重卡輪轂軸承各部分間的配合比一般小型汽車輪轂軸承更為精密，在結(jié)構(gòu)上比后者大2 倍以上。某公司研發(fā)生產(chǎn)了一款500 000 km 免維護(hù)的重卡輪轂軸承單元，并實現(xiàn)整個集成化的單元質(zhì)量較以前產(chǎn)品減輕110 kg[1]，由于結(jié)構(gòu)與密封唇形狀的改變，該輪轂軸承的摩擦性能需要重新測定。

該重卡輪轂軸承的主要組成部分由內(nèi)圈、外圈、滾動體、保持架、鎖扣扭、重卡車架(球殼)及螺栓等組成。根據(jù)汽車輪轂軸承的裝配要求和特點，軸承需要能夠同時承受較大的軸向載荷和徑向載荷，因此重卡輪轂軸承單元內(nèi)部裝配雙列圖錐滾子軸承。其工作特點主要為:將重卡輪轂軸承單元裝配在重卡的車橋主軸上，軸承的內(nèi)圈固定在主軸上，由汽車控制車橋主軸旋轉(zhuǎn)帶動重卡球殼旋轉(zhuǎn)，從而帶動外圈和滾子依據(jù)旋轉(zhuǎn)速度旋轉(zhuǎn)，重卡的輪胎裝在螺栓上從而驅(qū)動重卡行走[2]。重卡輪轂軸承的三維圖，如圖1、圖2 所示。

圖1 重卡輪轂軸承三維圖

圖2 重卡輪轂軸承爆炸圖

2 軸承摩擦力矩計算理論

軸承的內(nèi)部充滿了摩擦阻力，克服摩擦阻力使得軸承的內(nèi)外圈產(chǎn)生相對轉(zhuǎn)動的外加的力矩總和即為軸承的摩擦力矩[3]。在測量摩擦力矩時，啟動摩擦力矩和勻速旋轉(zhuǎn)摩擦力矩往往被作為評定軸承性能的重要技術(shù)指標(biāo)，摩擦力矩是體現(xiàn)整個輪轂軸承加工和裝配的綜合性指標(biāo)，軸承零件的制造工藝和裝配精度都會影響摩擦力矩的數(shù)值大小。一般將軸承摩擦力矩產(chǎn)生原因歸于以下4 點[4]。

1）滾動摩擦阻力：由于軸承的材料存在彈性滯后現(xiàn)象，從而會產(chǎn)生純滾動摩擦。2）軸承內(nèi)外圈及滾動體接觸處的滑動摩擦力：內(nèi)外圈及滾動體接觸過程中會產(chǎn)生差動滑動、自轉(zhuǎn)和傾斜等，從而引起的微小滑動，進(jìn)而產(chǎn)生滑動摩擦力。3）滾動體與保持架兜孔間的滑動、保持架與套圈引導(dǎo)表面間的滑動以及滾子端面與套圈擋邊之間的滑動摩擦力：滾子端面與套圈擋邊之間的滑動摩擦對屬于雙列圓錐滾子軸承的重卡輪轂軸承的影響尤其重要。4）攪動阻力引起的粘性摩擦阻力：由于潤滑劑的作用，在軸承運轉(zhuǎn)時產(chǎn)生。

SKF 公司根據(jù)這幾個方面的摩擦力矩影響情況推出了目前最常用的軸承摩擦力矩計算公式，非常適用于中等負(fù)載、潤滑良好及中等轉(zhuǎn)速等環(huán)境條件下。該公式如下[5]：

滯后性-滾動摩擦力矩：

式中：Grr—與軸承類型、尺寸、所受載荷相關(guān)的變量；

n—軸承轉(zhuǎn)速，r/min；

v—潤滑油運動粘度，mm2/s。

接觸點的線速度不同-滑動摩擦力矩：

式中：Gs1—與軸承類型、軸承直徑、所受載荷相關(guān)的變量；

μs1—滑動摩擦因數(shù)，與潤滑脂材料相關(guān)。

密封圈-內(nèi)外圈摩擦：

式中：Ks1—軸承類別的常量；

Ks2—軸承類別、密封圈形式相關(guān)的常量；

ds—肩部直徑，mm；

β—軸承類型、密封圈形式相關(guān)常量。

潤滑脂攪拌阻力：

式中：h—油膜厚度，mm；

αoil—粘壓系數(shù)，mm2/s；

S—油膜潤滑系數(shù)；

Dm—軸承節(jié)圓直徑，mm；

aε、ai—內(nèi)外滾道膨脹系數(shù)；

Z—鋼球粒數(shù)。

3 重卡輪轂軸承動力學(xué)仿真分析

3.1 重卡輪轂軸承單元三維建模

本文針對某公司開發(fā)的13T 免維護(hù)重卡輪轂軸承進(jìn)行Solidworks 三維建模，軸承主參數(shù)如表1 所示。

表1 重卡輪轂軸承主參數(shù)

3.2 材料屬性與約束設(shè)置

材料選擇:根據(jù)重卡輪轂軸承實際生產(chǎn)資料，分別對重卡輪轂軸承各部件的材料參數(shù)進(jìn)行設(shè)置，具體材料參數(shù)如表2 所示。

表2 圓錐滾子軸承各部件材料參數(shù)表

約束設(shè)置：重卡在行駛過程中，輪轂軸承內(nèi)圈固定，外圈隨著球殼轉(zhuǎn)動。在內(nèi)圈與大地之間添加平行副，在外圈與球殼之間添加固定副，在球殼與大地之間添加旋轉(zhuǎn)副[6]，分別添加滾子與內(nèi)圈、滾子與外圈、滾子與保持架之間的接觸，具體約束類型與約束數(shù)量如表3所示。

表3 重卡輪轂軸承各零部件約束

3.3 Adams 動力學(xué)仿真分析

在Solidworks 軟件中完成三維建模后，將Solidworks 中的三維模型導(dǎo)入到Adams 動力學(xué)仿真軟件，該操作需要先將Solidworks 裝配體保存為x_t 格式。打開Adams 軟件的File Type 選項，文件類型選擇Parasolid，在file to read 選項框中雙擊，選擇保存的x_t 文件，在Model name 中右鍵，依次選擇MODEL，creat，點擊ok完成導(dǎo)入。為重卡輪轂軸承單元添加模擬真實情況下所存在的接觸，選擇內(nèi)圈固定，外圈旋轉(zhuǎn)，并為整體添加約束副和轉(zhuǎn)速，由于該模型接觸較多，且滾動體轉(zhuǎn)動條件復(fù)雜，可以先在Settings-Solver-Executable 中將Thread Count（線程數(shù)）設(shè)置為8 線程，從而加快求解速度。點擊Simulation 模塊進(jìn)行動力學(xué)仿真，將時長設(shè)為1 s，步數(shù)設(shè)為2 000 步。設(shè)置完成后如圖3 所示[7]。

圖3 ADAMS 動力學(xué)仿真界面

3.3.1 啟動摩擦力矩仿真

為了得到啟動摩擦力矩的仿真值，設(shè)置轉(zhuǎn)速為5 r/min 模仿軸承從靜止開始轉(zhuǎn)動的過程。Adams 動力學(xué)仿真完成后，按F8 進(jìn)入Adams PostProcessor 模塊，其摩擦力矩如圖4 所示，軸承在剛剛啟動經(jīng)過一個過沖量后摩擦力矩最大，隨后慢慢平穩(wěn)得到啟動摩擦力矩。

圖4 啟動摩擦力矩仿真結(jié)果

去除初始過沖量數(shù)據(jù)后，將穩(wěn)定后的數(shù)據(jù)導(dǎo)入Matlab 處理，為了減少誤差，此項步驟重復(fù)操作10 次，取出每次仿真的最大摩擦力矩和啟動摩擦力矩，取平均值。處理得到轉(zhuǎn)速穩(wěn)定部分的平均值為12.17 N·m。

3.3.2 旋轉(zhuǎn)摩擦力矩仿真

旋轉(zhuǎn)摩擦力矩也稱動態(tài)摩擦力矩，即高轉(zhuǎn)速穩(wěn)定狀態(tài)下的摩擦力矩，為了得到啟動摩擦力矩的仿真值，設(shè)置轉(zhuǎn)速為600 r/min，模仿軸承勻速轉(zhuǎn)動的過程[7]。Adams 動力學(xué)仿真完成后，按F8 進(jìn)入Adams PostProcessor 模塊，其摩擦力矩如圖5 所示.

圖5 旋轉(zhuǎn)摩擦力矩仿真結(jié)果

與啟動摩擦力矩處理方法相同，旋轉(zhuǎn)摩擦力矩數(shù)據(jù)平均值為7.96 N·m，根據(jù)團(tuán)體標(biāo)準(zhǔn)《商用車輪轂軸承單元》，該檢測裝置對于13T 重卡輪轂軸承旋轉(zhuǎn)摩擦力矩的合格標(biāo)準(zhǔn)為4～8 N·m，啟動力矩在10～15 N·m 之間，本次動力學(xué)仿真的啟動摩擦力矩與旋轉(zhuǎn)摩擦力矩仿真數(shù)據(jù)都在該范圍內(nèi)，初步驗證了試驗對象在按照標(biāo)準(zhǔn)的模擬下摩擦力矩的合理性，為了進(jìn)一步驗證Adams 動力學(xué)分析對重卡輪轂軸承摩擦力矩仿真的準(zhǔn)確性，將進(jìn)行臺架試驗分析。

4 測試設(shè)備

本試驗采用的重卡輪轂軸承摩擦力矩檢測裝置如圖6 所示，圖中：1 為機(jī)架；2為頂升裝置；3為試驗工裝；4為驅(qū)動電機(jī)；5為氣液增壓缸；6為升降氣缸；7為扭矩傳感器；8為下壓頭。該裝置主要由扭矩傳感器、氣液增壓缸、頂升部件和定位部件組成。其中，定位部件由預(yù)前擋氣缸組件、后阻擋氣缸組件及RFID 模塊等組成，用于對托盤的定位與管理；頂升裝置部件由下頂升組件、下滑臺橫移擋塊組件及電機(jī)旋轉(zhuǎn)模塊等組成，用于旋轉(zhuǎn)軸承內(nèi)圈及定位產(chǎn)品高度；上滑臺升降部件由上滑臺組件、扭矩檢測模塊等組成，用于撥動外圈并檢測輪轂的轉(zhuǎn)動力矩。該設(shè)備的實體圖如圖7 所示。

圖6 摩擦力矩在線檢測裝置結(jié)構(gòu)圖

圖7 重卡輪轂軸承摩擦力矩檢測裝置實體圖

5 重卡輪轂軸承摩擦力矩試驗結(jié)果

為了保證測量的精確性，降低誤差影響，本次試驗選取100 個試驗工裝，每個試驗工裝分別測試2 次，共得到200 個數(shù)據(jù)[8]。從輪轂軸承單元靜止、啟動至設(shè)定轉(zhuǎn)速為止。試驗機(jī)在每1 次檢測完摩擦力矩之后，會把啟動摩擦力矩和旋轉(zhuǎn)摩擦力矩記錄在人機(jī)盒中，通過拷貝可以得到完整的EXCEL 數(shù)據(jù)文件，提取重卡輪轂軸承摩擦力矩曲線數(shù)據(jù)，將曲線數(shù)據(jù)導(dǎo)入MATLAB 中進(jìn)行繪圖[9]?？梢缘玫剑簡恿氐钠骄禐?1.81 N·m，勻速旋轉(zhuǎn)力矩的平均值為7.43 N·m。折線圖如圖8 所示。

圖8 重卡輪轂軸承單元摩擦力矩試驗數(shù)據(jù)

6 結(jié)語

利用Adams 動力學(xué)仿真軟件對重卡輪轂軸承進(jìn)行動力學(xué)仿真分析，得到13T 重卡輪轂軸承啟動摩擦力矩以及旋轉(zhuǎn)摩擦力矩的仿真結(jié)果，仿真結(jié)果均符合最新的行業(yè)標(biāo)準(zhǔn)，隨后通過重卡輪轂軸承摩擦力矩試驗機(jī)檢測其真實的摩擦力矩數(shù)據(jù)并與仿真結(jié)果比較，結(jié)果顯示仿真與試驗的數(shù)據(jù)相近，由此肯定了本次針對重卡輪轂軸承的動力學(xué)仿真可靠性，且可以得出：該重卡輪轂軸承的啟動摩擦力矩在10～13 N·m 之間，勻速轉(zhuǎn)動摩擦力矩在6～9 N·m 之間。該重卡輪轂軸仿真分析結(jié)果真實可靠，且與真實臺架試驗結(jié)果基本一致，可以應(yīng)用于實際工程分析中，指導(dǎo)生產(chǎn)實踐。