仿真技術(shù)在滾動阻力測試設(shè)備研發(fā)中的應(yīng)用

張紹國 高 峰 徐國艷 崔 瑩

(北京航空航天大學(xué) 交通科學(xué)與工程學(xué)院，北京100191)

滾動阻力是輪胎滾動時與路面接觸變形產(chǎn)生的遲滯能量損失，以發(fā)熱為主要形式發(fā)散到大氣中.降低輪胎的滾動阻力可以節(jié)約油耗，減少溫室氣體排放.歐盟和美國等國家頒布的輪胎標簽法規(guī)對輪胎滾動阻力提出了分級要求，不滿足滾動阻力認證的輪胎禁止在相關(guān)國家銷售［1］.

輪胎結(jié)構(gòu)設(shè)計中的能量損耗模型主要對輪胎的滾動阻力進行預(yù)測分析，尋求減少滾動阻力的結(jié)構(gòu)設(shè)計方法，但不能滿足輪胎的認證需求［2-3］.滾動阻力測試是高負荷下測量小分力的試驗，早期的輪胎滾動阻力試驗在底盤測功機上通過反拖法進行［4］，室內(nèi)專用的單滾筒大直徑測試設(shè)備，對控制和測試精度要求較為嚴格，高精度的多分量傳感器和合理的定位誤差補償算法直接影響著設(shè)備的產(chǎn)品質(zhì)量.本文基于輪胎滾動阻力專用檢測設(shè)備的結(jié)構(gòu)原理，建立了滾動阻力的動力學(xué)測試模型，模擬設(shè)備各種工況下的測試過程，并利用正交回歸擬合算法對測試數(shù)據(jù)進行誤差修正，最后與試驗數(shù)據(jù)進行對比驗證.

1 測試原理與誤差分析

1.1 測試方法

設(shè)備采用轉(zhuǎn)鼓式結(jié)構(gòu)，GB/T 18861—2002及ISO 28580等滾動阻力試驗標準中規(guī)定滾動阻力測試方法包括測力法和測扭矩法等［5］.由于轉(zhuǎn)鼓的起動扭矩要比穩(wěn)速運轉(zhuǎn)扭矩大得多，且測量必須在無加速扭矩輸入時進行，測試上難于實現(xiàn)，精度很難保證，故現(xiàn)在國外設(shè)備滾動阻力的計算模式多采用測力法，測試原理如圖1所示，通過測量輪胎軸上的反作用力以避免轉(zhuǎn)鼓轉(zhuǎn)動慣性對測量結(jié)果的影響［6］.滾動阻力和滾動阻力系數(shù)換算公式如式(1)和式(2)所示.

式中，F(xiàn)r為輪胎滾動阻力;Ft為輪軸反作用力;L為輪軸中心到轉(zhuǎn)鼓面的距離;R為轉(zhuǎn)鼓半徑.

式中，Cr為輪胎滾動阻力系數(shù);Fy為輪胎負荷.

圖1 滾動阻力測試原理圖

1.2 設(shè)備原理

設(shè)備主要有床身、轉(zhuǎn)鼓、徑向加載裝置和輪軸裝置等組成，如圖2所示.旋轉(zhuǎn)的輪胎與轉(zhuǎn)鼓之間的接觸壓力使輪胎產(chǎn)生變形，并在接觸點的切線方向產(chǎn)生了滾動阻力，通過輪軸的反作用力傳遞給傳感器板.

圖2 試驗設(shè)備結(jié)構(gòu)圖

1.3 誤差分析

設(shè)備的安裝精度直接影響測試數(shù)據(jù)的可靠性.機械定位誤差包括輪軸作用方向與轉(zhuǎn)鼓接觸中心平面的對中誤差;輪軸的中心平面與轉(zhuǎn)鼓的接觸中心平面的垂直誤差(外傾角);輪軸的行進方向與轉(zhuǎn)鼓中心平面方向間的夾角(側(cè)偏角).輪軸與轉(zhuǎn)鼓基準線的對中偏差使輪胎負荷在滾動阻力方向產(chǎn)生偏移分量，影響測試數(shù)據(jù)精度，是設(shè)備的主要誤差源.

設(shè)備的系統(tǒng)誤差還包括傳感器零點漂移以及軸承摩擦等，這些因素產(chǎn)生的縱向力視為附加損失處理，可以通過設(shè)備標定消除，本文重點對設(shè)備對中定位誤差進行討論，數(shù)據(jù)的誤差修正算法按照式(3)進行處理.

式中，F(xiàn)m為輪軸力測量值;Fd為動態(tài)誤差.

2 三分力傳感器

2.1 傳感器結(jié)構(gòu)

設(shè)備的試驗數(shù)據(jù)通過傳感器測量輪軸力來實現(xiàn)，高精度的三分力傳感器是設(shè)備測試系統(tǒng)核心部件.靈敏度、維間耦合等因素與傳感器的性能密切相關(guān)，基于有限元的傳感器模型具備結(jié)構(gòu)基準坐標明確和對稱性好等優(yōu)點，目前國內(nèi)外已經(jīng)采用有限元方法對傳感器進行結(jié)構(gòu)設(shè)計［7］.

十字型結(jié)構(gòu)傳感器的主要變形部分由4個十字型筋和4個浮動支撐薄片組成，如圖3所示.在力的作用下，與力方向垂直的支撐薄片可以自由活動，該方向的十字型筋沒有變形，而與力作用方向垂直的十字型筋則發(fā)生彎曲變形，從而避免維力之間的干擾.通過分析應(yīng)變區(qū)域彈性體的應(yīng)變分量，將力分解為各自獨立的坐標成分，計算得到彈性體網(wǎng)格區(qū)域的應(yīng)變量.

圖3 傳感器結(jié)構(gòu)和有限元模型

2.2 模擬標定

傳感器滾動阻力方向的標定方式采用砝碼作為力源加載.4個靈敏度完全一致的傳感器以測力板形式與設(shè)備裝配在一起，傳感器分別安裝在臺架和傳感器板之間對應(yīng)的角部位置，命名為D1，D2，D3，D4，如圖4所示.

圖4 傳感器布置方案

靜負荷下，在法蘭參考點沿滾動阻力方向分別施加大小不同的負載，同時采集4個傳感器的應(yīng)變量作為樣本數(shù)據(jù)，按式(4)對數(shù)據(jù)進行和差處理得到總應(yīng)變量D，標定數(shù)據(jù)如表1所示.

表1 滾動阻力方向標定數(shù)據(jù)

依據(jù)最小二乘法對采集的傳感器樣本數(shù)據(jù)進行計量單位轉(zhuǎn)換，應(yīng)變量D和力F的擬合方程為

式中，k，b為回歸系數(shù);計算求出 k=1.27;b=0.23.

3 測試模型

3.1 輪胎有限元模型

輪胎以某成品胎11R22.5為研究對象.輪胎中的橡膠大部分填充碳黑材料，采用Yeoh三次方程，如式(6)所示.

式中，W為橡膠的應(yīng)變能函數(shù);I1為第一應(yīng)變不變量;C10，C20，C30為常數(shù).

結(jié)構(gòu)阻尼能夠很好地描述輪胎滾動時的遲滯特性，為限制測試過程中產(chǎn)生的數(shù)值振蕩，提高計算精度，輪胎模型采用結(jié)構(gòu)阻尼模型［8］.

式中，M為質(zhì)量矩陣;P為外力;K為剛度矩陣常數(shù);C為阻尼矩陣常數(shù).

滾動阻力分析中要考慮熱力耦合現(xiàn)象［9-10］.計算輪胎的溫度場分布時，主要考慮胎體熱源與輪胎外界的熱交換.輪胎的有限單元網(wǎng)格熱能Q由式(8)計算［11］.

式中，σp和εp分別為單元應(yīng)力和應(yīng)變量;ω為輪胎滾動角速度;tanδ為對應(yīng)材料的損耗正切值.

3.2 邊界條件

在輪輞和輪軸裝置之間設(shè)置鉸接屬性，使輪胎在輪軸上具有旋轉(zhuǎn)自由度.輪胎的轉(zhuǎn)動通過轉(zhuǎn)鼓驅(qū)動.穩(wěn)態(tài)條件下，通過對轉(zhuǎn)鼓施加恒定速度，驅(qū)動輪胎以設(shè)定速度做旋轉(zhuǎn)運動.圖5為輪胎在穩(wěn)態(tài)條件下的計算模型.

圖5 滾動阻力測試模型

4 試驗驗證

單胎多點測試數(shù)據(jù)重復(fù)性精度決定著設(shè)備的測試能力.輪軸在大負荷作用下產(chǎn)生變形，加載方向與轉(zhuǎn)鼓中心線產(chǎn)生一定偏離，負荷對輪軸力Ft產(chǎn)生橫向干擾，同時由于控制精度等因素設(shè)備產(chǎn)生了系統(tǒng)累計誤差Fd.試驗數(shù)據(jù)在美國MTS公司轉(zhuǎn)鼓式輪胎滾動阻力試驗機上進行，如圖6所示.輪軸與轉(zhuǎn)鼓軸中心基準位置標定偏差為0.04°，為保證設(shè)備測試精度，需對測試數(shù)據(jù)進行補償.

圖6 滾動阻力試驗

4.1 修正模型

測力法通常采用輪胎正反轉(zhuǎn)，用二者輪軸力的平均值作為補償結(jié)果，但正反轉(zhuǎn)測到的輪軸力差異較大，補償精度較低.有限元模型能夠模擬設(shè)備各種理想狀態(tài)下的工作特性，從而建立合理的誤差補償模型.圖7為額定工況下輪胎基準位置和偏差位置測試數(shù)據(jù)對比曲線，在偏差狀態(tài)下，負荷對滾動阻力方向測得的輪軸力產(chǎn)生干擾，因而偏差位置的測試數(shù)據(jù)幅值波動較大，與基準位置的數(shù)據(jù)產(chǎn)生一定差異.

圖7 滾動阻力應(yīng)變波形

一次正交回歸試驗可以建立試驗指標與多個試驗因素之間的一次回歸關(guān)系，具有在每個試驗點獲得最大有用信息的特點，使統(tǒng)計分析具有較優(yōu)性質(zhì).

采用L8(27)正交表進行試驗，對設(shè)備各控制因素進行上、下水平編碼，編碼如下:

式中，v為試驗速度;p為氣壓.

按照不同工況下進行水平試驗，測試數(shù)據(jù)如表2所示.經(jīng)過編碼后建立各控制參數(shù)和測試結(jié)果的一元正交回歸方程如下:

其中，F(xiàn)'d為在仿真模型中得到的誤差值.

表2 正交試驗數(shù)據(jù)

由表2得:a0=0.435，a1=0.359，a2=0.029，a3=-0.056.

方程的回歸系數(shù)不受試驗因素量綱和數(shù)值的影響，回歸系數(shù)反映該因素影響的大小.代入回歸方程，化簡后得到指數(shù)方程(13)，公式表明負荷和胎壓對測試數(shù)據(jù)誤差值影響權(quán)重較大.

4.2 趨勢性試驗

圖8和圖9為額定速度下滾動阻力與負載和充氣壓力的變化關(guān)系，仿真結(jié)果和試驗數(shù)據(jù)的變化趨勢一致.輪胎滾動阻力與負載成正比關(guān)系，負載增加使輪胎變形變大，輪胎的內(nèi)摩擦滾動阻力也隨之增大.

圖8 滾動阻力與負載變化關(guān)系

滾動阻力隨氣壓變大而減小，主要是氣壓升高使輪胎的胎體簾線張緊，剛度增大，相同負載下輪胎的下沉量減小，導(dǎo)致滾動阻力減小［12］.

圖9 滾動阻力與充氣壓力變化關(guān)系

4.3 重復(fù)性試驗

ISO—28580中規(guī)定轎車和載重輪胎滾動阻力測試標準偏差不能超過0.2 N·kN-1，滿足此條件設(shè)備具備對輪胎進行質(zhì)量分級的能力.由于精確的輪胎有限元模型構(gòu)建困難，并且為提高計算效率，對測試模型進行合理簡化，將轉(zhuǎn)鼓、輪輞作為剛體部件處理，仿真結(jié)果比試驗數(shù)據(jù)高出20%左右.仿真結(jié)果和實驗數(shù)據(jù)之間的關(guān)聯(lián)模型如式(14)所示:

η受負載脈動、速度波動和閉氣工況等因素影響，由于變化范圍較小，體現(xiàn)在誤差數(shù)據(jù)上量值有限，為簡化計算，取平均值0.8.

式(13)不僅對定位誤差進行修正，還能夠消除設(shè)備由于控制精度等因素對測試數(shù)據(jù)產(chǎn)生的擾動.在MTS滾動阻力試驗設(shè)備上，對輪胎進行6次相同工況下的閉氣試驗，結(jié)果如表3所示.修正后滾動阻力系數(shù)試驗精度為0.0679，設(shè)備滿足趨勢性試驗標準要求.

表3 重復(fù)性試驗數(shù)據(jù)

5 結(jié)論

針對輪胎滾動阻力測試大負荷小分量的特點，通過建立設(shè)備的有限元測試模型，分析設(shè)備在各種狀態(tài)下的測試特性.計算結(jié)果表明，以測力板形式安裝的傳感器布局方案能夠?qū)崿F(xiàn)對目標測試量的有效采集，在相同的控制參數(shù)下，仿真結(jié)果與試驗數(shù)據(jù)表征的趨勢特征一致.在設(shè)備現(xiàn)有工況的基礎(chǔ)上，利用一次正交回歸方程建立的誤差補償算法綜合考慮設(shè)備機械定位及控制精度等因素對測試數(shù)據(jù)產(chǎn)生的擾動，能夠減少系統(tǒng)的累計誤差，提高現(xiàn)有設(shè)備的測試精度.該結(jié)論為高性能輪胎滾動阻力測試設(shè)備的研制提供理論基礎(chǔ)和數(shù)據(jù)處理方法.

References)

［1］ Jerome B，Jason B.Reducing tire rolling resistance to save fuel and lower emissions［J］.SAE International Journal of Passenger Cars-Mechanical Systems，2009，1(1):9-17

［2］ Jun A，Chen J F，Michael Rotter J，et al.Assessment of rolling resistance models in discrete element simulations［J］.Powder Technology，2011，206(3):269-282

［3］ Shida Z，Koishi M，Kogure T，et al.A rolling resistance simulation of tires using static finite elements analysis［J］.Tire Science and Technology，1999，27(2):84-105

［4］ 王建強，何鳳江，張立斌，等.汽車動力性檢測模型的建立［J］.中國公路學(xué)報，2001:14(3):109-112 Wang Jianqiang，He Fengjiang，Zhang Libin，et al.A model for testing automobile dynamic property［J］.China Journal of Highway and Transport，2001，14(3):109-112(in Chinese)

［5］任禮行，劉青，張艾謙，等.輪胎滾動阻力測量與分析［J］.汽車工程，2000，22(5):316-319 Yam Laihang，Liu Qing，Zhang Aiqian，et al.Measurement and analysis of tire rolling resistance ［J］.Automotive Engineering，2000，22(5):316-319(in Chinese)

［6］韓秀芝.子午線輪胎穩(wěn)態(tài)滾動阻力及水滑特性的研究［D］.北京:北京化工大學(xué)機電工程學(xué)院，2009:7-10 Han Xiuzhi.Research on steady-state rolling resistance and hydroplaning of radial tire［D］.Beijing:School of Mechanical and Electrical Engineering，Beijing University of Chemical Technology，2009:7-10(in Chinese)

［7］王國泰，易秀容，王理麗.六維力傳感器發(fā)展中的幾個問題［J］.機器人，1997，19(6):447-478 Wang Guotai，Yi Xiurong，Wang Lili.Problems related to the development ofsix axes force/torque sensor［J］.Robot，1997，19(6):474-478(in Chinese)

［8］管迪華，曾祥生，范成建.輪胎動態(tài)模型的阻尼和對滾動阻力及動態(tài)響應(yīng)影響的分析［J］.汽車工程，2006，28(7):643-646 Guan Dihua，Zeng Xiangsheng，F(xiàn)an Chengjian.Measurement and analysis of tire rolling resistance［J］.Automotive Engineering，2006，28(7):643-646(in Chinese)

［9］ Lin Y J，Hwang S J.Temperature prediction of rolling tires by computer simulation［J］.Mathematics and Computers in Simulation，2004，67(3):235-249

［10］ Yin H S，Hu Y S，Zhang H，et al.Truck tire thermal-mechanical FEA and DMA with application to endurance evaluation［J］.Tire Science and Technology，2006，34(4):220-236

［11］龔科家，危銀濤，葉進雄.基于熱力學(xué)有限元分析的輪胎滾動阻力仿真［J］.中國機械工程，2009，20(5):626-629 Gong Kejia，Wei Yintao，Ye Jinxiong.Simulation techniques of tire rolling resistance based on thermal-mechanical FEA［J］.China Mechanical Engineering，2009，20(5):626-629(in Chinese)

［12］莊繼德.汽車輪胎學(xué)［M］.北京:北京理工大學(xué)出版社，1995:180-183 Zhuang Jide.Vehicle tire study［M］.Beijing:Beijing Institute of Technology Press，1995:180-183(in Chinese)