大扁平比胎側(cè)解析剛度建模及輪胎動力學(xué)分析

劉志浩 劉釔汛 劉巖 高欽和 馬棟 黃通

摘要：針對重載輪胎大扁平比結(jié)構(gòu)建模問題，從動力學(xué)建模、實驗?zāi)B(tài)分析、結(jié)構(gòu)參數(shù)辨識等方面，基于解析彈性基礎(chǔ)的歐拉梁模型，對重載輪胎的柔性胎體和大扁平比胎側(cè)曲梁的低頻動力學(xué)特性開展研究，建立了考慮充氣預(yù)緊力的歐拉梁胎體模型，利用實驗?zāi)B(tài)方法，探究了不同充氣壓力下的柔性胎體振動特性;考慮胎側(cè)曲梁預(yù)緊力弦效應(yīng)和結(jié)構(gòu)彎曲效應(yīng)，建立了大扁平比胎側(cè)曲梁解析剛度模型;基于模態(tài)測試結(jié)果，進(jìn)行柔性胎體與解析胎側(cè)結(jié)構(gòu)參數(shù)辨識。研究結(jié)果表明：在0～180 Hz頻率范圍內(nèi)，重載輪胎以結(jié)構(gòu)周向彎曲振動為主，可利用基于彈性基礎(chǔ)的柔性梁模型表征;大扁平比胎側(cè)曲梁的解析剛度與胎側(cè)的幾何、結(jié)構(gòu)和充氣壓力參數(shù)直接相關(guān);輪胎充氣壓力影響柔性胎體梁的軸向預(yù)緊力和胎側(cè)的弦剛度，進(jìn)而影響輪胎彎曲振動特性。

關(guān)鍵詞：輪胎動力學(xué);重載子午胎;胎側(cè)曲梁解析剛度模型;彈性基礎(chǔ)柔性梁

中圖分類號：U461.2;U463.341-.3

文獻(xiàn)標(biāo)志碼：A

文章編號：10044523（ 2022）01-015909

DOI： 10.1638 5/j .cnki.issn.10044523.2022.01.017

引 言

重載輪胎[1]作為軍用車輛與地面直接接觸的部件，除去空氣阻力、坡度阻力和加速阻力等少數(shù)作用力，其他作用力都將通過輪胎作用在整車上，因此其力學(xué)特性直接影響車輛動力學(xué)特性。重載輪胎具有氣壓高、阻尼低、花紋粗大、較大扁平比的特點，在承載過程中，胎體、胎側(cè)產(chǎn)生結(jié)構(gòu)柔體變形，同時在滾動過程中，輪胎胎體與路面直接接觸，將路面不平度激勵經(jīng)胎側(cè)、容腔，傳遞至輪輞，其中，胎體一胎側(cè)一輪輞振動傳遞路徑屬于結(jié)構(gòu)振動[2]，胎體一空腔一輪輞振動傳遞路徑屬于結(jié)構(gòu)一空腔耦合振動，則輪輞的傳遞響應(yīng)[3]由結(jié)構(gòu)振動響應(yīng)和結(jié)構(gòu)一空腔耦合振動構(gòu)成。為此，為表征輪胎滾過任意路面時的振動響應(yīng)，需對重載輪胎的結(jié)構(gòu)柔性[4]開展研究，考慮胎體和胎側(cè)柔性建立的基于彈性基礎(chǔ)的柔性胎體模型[5]成為輪胎結(jié)構(gòu)模型的典型代表，其中柔性胎體可模擬輪胎胎體的變形特征，而胎側(cè)單元作為連接柔性胎體和輪輞的部件，可等效為具有剛度特性的彈簧，起著傳遞力和變形的作用。

胎側(cè)作為胎體和輪輞的連接和傳力部件，其力學(xué)特征和振動特性影響胎體與輪輞間的傳遞特性，基于彈性基礎(chǔ)的柔性胎體模型將胎側(cè)假設(shè)為線性、沿圓周均勻分布的彈簧，包含胎側(cè)一維徑向剛度的力學(xué)特性模型[6]、具有徑向和切向剛度的胎側(cè)二維剛度彈簧模型等。文獻(xiàn)[8]基于徑向和切向彈簧建立了基于彈性基礎(chǔ)的環(huán)模型，研究了不同胎側(cè)剛度和充氣壓力對輪胎振動特性的影響規(guī)律;Liu等建立了考慮胎側(cè)彎曲剛度的基于彈性連續(xù)基礎(chǔ)的二維胎側(cè)剛度模型[9]，研究了胎體環(huán)與連續(xù)胎側(cè)環(huán)的耦合振動特性。文獻(xiàn)[10]建立了胎側(cè)穩(wěn)態(tài)徑向和切向剛度的解析解。針對任意彈性基礎(chǔ)的柔性環(huán)振動問題，國內(nèi)外學(xué)者開展了研究，該類模型可用于分析具有非均勻胎側(cè)剛度和非均勻質(zhì)量特性輪胎的動態(tài)響應(yīng);Matsubara等研究了基于一般彈性基礎(chǔ)的柔性環(huán)白由振動問題[11]，考慮空間圓周方向的非均勻剛度，獲得了固有頻率和陣型的解析解;Cooley等[12]基于負(fù)載和變形敏感的剛度模型，研究了非線性剛度的彈性基礎(chǔ)對柔性環(huán)動態(tài)特征的影響規(guī)律。普通小扁平比輪胎模型，通常將胎側(cè)等效為線性彈簧[13]，忽略其慣性和非線性剛度特性;而對于大扁平比重載輪胎，其大扁平比胎側(cè)的慣性特性及變形過程中由于結(jié)構(gòu)變化易造成非線性剛度凸顯，需在胎側(cè)剛度建模中將胎側(cè)的拉伸、剪切等結(jié)構(gòu)變形特征考慮在內(nèi)，建立可表征大扁平比胎側(cè)曲梁結(jié)構(gòu)和材料特征的剛度模型。

1 問題描述

重載輪胎在0--180 Hz頻率范圍內(nèi)，表現(xiàn)為胎體的周向彎曲振動，與基于彈性基礎(chǔ)的歐拉梁模型一致。論文以基于彈性基礎(chǔ)的歐拉梁模型為基礎(chǔ)，分別從動力學(xué)建模、實驗?zāi)B(tài)分析和結(jié)構(gòu)參數(shù)辨識等三個方面對模型中的柔性胎體和大扁平比胎側(cè)曲梁的動力學(xué)特性開展研究。基于歐拉梁理論，探究充氣壓力對胎體彎曲振動的影響規(guī)律;建立考慮預(yù)緊力弦效應(yīng)和結(jié)構(gòu)彎曲效應(yīng)的大扁平比胎側(cè)曲梁解析剛度模型，并研究胎側(cè)曲梁非均勻截面特性和幾何、結(jié)構(gòu)參數(shù)對解析剛度的影響規(guī)律。

圖1給出了一維基于彈性基礎(chǔ)的柔性梁輪胎模型的示意圖，包含柔性胎體梁、大扁平比胎側(cè)徑向曲梁和輪輞三部分組成，其中輪輞剛度較大，忽略其柔性，將其等效為剛體質(zhì)量塊，建立考慮胎體柔性和大扁平比胎側(cè)曲梁預(yù)緊力弦效應(yīng)和結(jié)構(gòu)變形效應(yīng)的一維輪胎模型。并開展模態(tài)試驗探究，對比不同充氣壓力作用下的重載輪胎模態(tài)特性，并基于彈性基礎(chǔ)的歐拉梁模型開展輪胎結(jié)構(gòu)參數(shù)辨識方法研究。

2 柔性胎體動力學(xué)建模及試驗

2.1 基于彈性基礎(chǔ)的歐拉梁輪胎模型

假設(shè)：（1）胎體梁截面的中心慣性軸在平面xoy內(nèi)，且外載荷作用在該平面內(nèi);（2）胎體梁在平面內(nèi)作橫向微振動;（3）在低頻振動時可以忽略剪切變形以及截面繞中性軸轉(zhuǎn)動慣量的影響。

2.2 基于實驗?zāi)B(tài)的胎體低頻彎曲特性研究

基于實驗?zāi)B(tài)分析方法，探究柔性胎體彎曲振動特性，尤其針對不同充氣壓力下的重載輪胎振動模態(tài)特性開展實驗研究。搭建重載輪胎振動模態(tài)測試系統(tǒng)，針對重載輪胎不同的充氣壓力工況，采用移動力錘的方法，依次激勵沿圓周方向等分的34個胎體點，獲取目標(biāo)傳遞函數(shù)。

為便于討論，定義激勵點A和點B與響應(yīng)點A，如圖4所示。點A與點B間隔90°，其中點A激勵一點A響應(yīng)獲取的傳遞函數(shù)為驅(qū)動點傳遞函數(shù)，點B激勵一點A響應(yīng)獲取的傳遞函數(shù)為跨點傳遞函數(shù)。分別在點A和點B進(jìn)行力錘激勵，將激勵力通過力傳感器，經(jīng)電荷放大器調(diào)理后，傳輸給DEWE-43高速采集器，同時利用PCB振動加速度傳感器，檢測響應(yīng)點A處的振動響應(yīng)，同時傳輸給DEWE-43高速采集器，分別計算驅(qū)動點傳遞函數(shù)和跨點傳遞函數(shù)，并對比分析不同充氣壓力下的輪胎傳遞特性。圖5中分別列舉了重載輪胎在不同充氣壓力，包括：0.3，0.4，0.5，0.6，0.7和0.8 MPa下驅(qū)動點的傳遞函數(shù)和相干函數(shù)，結(jié)果表明：（1）相干函數(shù)除低頻范圍內(nèi)，均在0.9以上，驗證了加速度實驗傳遞函數(shù)的有效性;（2）隨著充氣壓力的增大，輪胎的剛度增大，導(dǎo)致輪胎頻率共振峰朝著頻率增大的方向平移。gzslib202204041748

基于最小二乘復(fù)指數(shù)法14對不同充氣壓力下的輪胎振動模態(tài)和傳遞函數(shù)進(jìn)行辨識，圖6為識別的不同充氣壓力下的輪胎模態(tài)參數(shù)，結(jié)果表明：（1）基于最小二乘復(fù)指數(shù)法的驅(qū)動點傳遞函數(shù)擬合精度低于基于最小二乘復(fù)指數(shù)法的跨點傳遞函數(shù)擬合精度。分析原因為：由于胎體花紋較大阻尼的影響導(dǎo)致驅(qū)動點傳遞函數(shù)衰減嚴(yán)重，造成擬合精度較低，但基于最小二乘復(fù)指數(shù)法的跨點擬合傳遞函數(shù)與實驗傳遞函數(shù)基本一致，驗證了基于最小二乘復(fù)指數(shù)法的傳遞特性分析的有效性;（2）隨著充氣壓力的增大，重載輪胎面內(nèi)各階固有頻率增加，而阻尼比降低;（3）充氣壓力對輪胎振動特性的影響，體現(xiàn)在兩個方面，包括：充氣壓力引起的胎體梁軸向預(yù)緊力和充氣壓力引起的胎側(cè)剛度變化。

3 基于虛功原理的胎側(cè)徑向剛度估計

輪胎胎側(cè)的徑向剛度主要體現(xiàn)在兩個方面：（1）充氣預(yù)緊力作用下的弦剛度。在充氣壓力作用下，在胎側(cè)和胎體的連接點處，存在由于充氣預(yù)緊力產(chǎn)生的拉伸力;（2）胎側(cè)結(jié)構(gòu)變形引起的結(jié)構(gòu)剛度。在變形過程中，胎側(cè)單元產(chǎn)生拉伸、剪切和彎曲變形，造成由于胎側(cè)結(jié)構(gòu)變形產(chǎn)生的結(jié)構(gòu)剛度。

3.1 考慮充氣預(yù)緊力的胎側(cè)弦剛度

圖8為胎側(cè)圓弧的受力分析圖，并建立基于充氣預(yù)緊壓力的胎側(cè)弦剛度解析模型。

3.2 考慮曲梁變形特征的胎側(cè)結(jié)構(gòu)剛度

輪胎胎側(cè)除充氣預(yù)緊力的弦效應(yīng)外，其在外力作用下，會產(chǎn)生拉伸力、剪切力和彎曲力矩，因此需將胎側(cè)結(jié)構(gòu)變形考慮在內(nèi)，如圖9所示;同時在力和力矩的作用下產(chǎn)生相應(yīng)的結(jié)構(gòu)變形，因此為準(zhǔn)確表征大扁平比胎側(cè)曲梁的剛度特性，提出胎側(cè)曲梁解析剛度建模方法。基于虛功原理，探究胎側(cè)單元在外力作用下產(chǎn)生的彎曲位移ub、拉伸位移ust和剪切位移ush。

4 柔性胎體與胎側(cè)結(jié)構(gòu)參數(shù)辨識

4.1 柔性胎體結(jié)構(gòu)參數(shù)辨識

基于彈性基礎(chǔ)的一維歐拉梁模型，建立了式（8）所示的重載輪胎模型，輪胎的振動特性與輪胎的結(jié)構(gòu)和幾何參數(shù)息息相關(guān)，輪胎的參數(shù)中，除幾何參數(shù)便于測量外，其結(jié)構(gòu)參數(shù)，包括胎體彎曲剛度EI、胎體密度ρA和胎側(cè)徑向剛度kr，無法實現(xiàn)直接測量。結(jié)合2.2節(jié)中獲取的不同充氣壓力下的實驗?zāi)B(tài)參數(shù)，建立關(guān)于結(jié)構(gòu)參數(shù)的重載輪胎模態(tài)共振頻率解析解，利用待定系數(shù)法，求取重載輪胎未知結(jié)構(gòu)參數(shù)。

胎體柔體動力學(xué)模型為偏微分方程組，利用模態(tài)疊加法，將偏微分方程組轉(zhuǎn)化為空間和時間的常微分方程進(jìn)行求解，推導(dǎo)出輪胎各階固有頻率與輪胎結(jié)構(gòu)參數(shù)間的關(guān)系。令：

式（19）為重載輪胎關(guān)于輪胎結(jié)構(gòu)、幾何和模態(tài)階數(shù)相關(guān)的模態(tài)共振頻率的解析形式，分別利用待定系數(shù)法，結(jié)合2.2節(jié)中獲取的不同充氣壓力下的輪胎實驗?zāi)B(tài)參數(shù)，對不同充氣壓力所對應(yīng)的式（19）中的結(jié)構(gòu)參數(shù)進(jìn)行識別，計算結(jié)果如圖10所示。

結(jié)果表明：（1）圖lO（a）為不同充氣壓力下的胎體彎曲剛度計算值，表明胎體彎曲剛度不受輪胎充氣壓力的影響，其平均胎體彎曲剛度為133.69 N/m;（2）圖10（b）為不同充氣壓力下的胎體密度計算值，表明胎體密度不受輪胎充氣壓力的影響，其平均胎體密度為32.50 kg/m2; （3）圖10（c）為不同充氣壓力下的胎側(cè)徑向剛度計算值，計算結(jié)果表明胎側(cè)徑向剛度近似與充氣壓力呈線性關(guān)系，與胎側(cè)解析剛度公式結(jié)論一致;（4）不同充氣壓力下的輪胎胎體彎曲剛度和胎體密度基本為定值，而胎體徑向剛度則近似于充氣壓力呈線性關(guān)系。

4.2考慮弦剛度和結(jié)構(gòu)剛度的大扁平比胎側(cè)曲梁

解析剛度模型

為計算胎側(cè)曲梁沿徑向的彈性模量，將圖10（c）中計算的胎側(cè)在0.8 MPa充氣壓力下的胎側(cè)剛度值代人式（17），通過計算Es= 1.3×10 7。

為驗證計算結(jié)果的準(zhǔn)確性，利用計算的胎側(cè)彈性模量參數(shù)，分析不同充氣壓力下的胎側(cè)解析剛度，與圖10（c）進(jìn)行對比，結(jié)果如圖11所示。

對比結(jié)果表明：（1）利用胎側(cè)彈性模量Es=1.3×10 7計算的不同充氣壓力下的胎側(cè)解析剛度與圖10（c）識別的胎側(cè)剛度擬合較好，驗證了胎側(cè)彈性模量參數(shù)和胎側(cè)解析剛度模型的準(zhǔn)確性;（2）充氣壓力影響胎側(cè)弦剛度，而不影響胎側(cè)結(jié)構(gòu)剛度;（3）充氣壓力低于0.44 MPa時，胎側(cè)結(jié)構(gòu)剛度大于胎側(cè)弦剛度，胎側(cè)結(jié)構(gòu)剛度為胎側(cè)解析剛度的主體，隨著充氣壓力的增大，胎側(cè)弦剛度大于胎側(cè)結(jié)構(gòu)剛度，成為胎側(cè)解析剛度的主體。

5 討 論

式（17）給出了胎側(cè)幾何、結(jié)構(gòu)參數(shù)構(gòu)成的胎側(cè)解析剛度表達(dá)式，為分析不同幾何和結(jié)構(gòu)參數(shù)對胎側(cè)剛度的影響規(guī)律，開展胎側(cè)解析剛度影響分析，主要包括胎側(cè)的結(jié)構(gòu)和幾何參數(shù)。

圖12描述了胎側(cè)結(jié)構(gòu)和幾何參數(shù)對胎側(cè)解析剛度的影響規(guī)律，包括胎側(cè)徑向彈性模量Es（圖12（a））、胎側(cè)截面厚度hs（圖12（b））、胎側(cè)曲梁弧長ls（圖12（c））和胎側(cè)曲梁弧度夾角θ（圖12（d））。結(jié)果表明：

（1）胎側(cè)徑向彈性模量影響胎側(cè)結(jié)構(gòu)剛度，成近似正向線性關(guān)系，彈性模量低于23.5 MPa時，胎側(cè)結(jié)構(gòu)剛度小于胎側(cè)弦剛度，隨著彈性模量的增大，胎側(cè)結(jié)構(gòu)剛度增大，彈性模量大于23.5 MPa時，胎側(cè)結(jié)構(gòu)剛度大于胎側(cè)弦剛度;

（2）胎側(cè)截面厚度影響胎側(cè)結(jié)構(gòu)剛度，成近似正向線性關(guān)系;

（3）胎側(cè)曲梁弧長影響胎側(cè)結(jié)構(gòu)剛度，隨著曲梁弧長的增加，胎側(cè)結(jié)構(gòu)剛度降低，成非線性變化趨勢，曲梁弧長大于0.225 m時，胎側(cè)弦剛度大于胎側(cè)結(jié)構(gòu)剛度，隨著曲梁弧長的減小，胎側(cè)結(jié)構(gòu)剛度增大，曲梁弧長小于0.225 m時，胎側(cè)結(jié)構(gòu)剛度大于胎側(cè)弦剛度;

（4）胎側(cè)曲梁夾角同時影響胎側(cè)弦剛度和胎側(cè)結(jié)構(gòu)剛度，隨著胎側(cè)曲梁夾角的增大，胎側(cè)弦剛度和結(jié)構(gòu)剛度降低，成非線性關(guān)系;

（5）為增大胎側(cè)徑向剛度，需增大胎側(cè)徑向方向的彈性模量，增加胎側(cè)曲梁厚度，減小胎側(cè)曲梁弧長和曲梁夾角。gzslib202204041748

6 結(jié) 論

針對重載輪胎大扁平比結(jié)構(gòu)建模問題，基于解析彈性基礎(chǔ)的歐拉梁模型，分別從動力學(xué)建模、實驗?zāi)B(tài)分析、結(jié)構(gòu)參數(shù)辨識等三個方面對重載輪胎的柔性胎體和大扁平比胎側(cè)曲梁的低頻動力學(xué)特性開展研究。研究結(jié)果表明：（1）在0--180 Hz頻率范圍內(nèi)，重載輪胎以結(jié)構(gòu)周向彎曲振動為主，與歐拉梁模型相一致，可利用基于彈性基礎(chǔ)的柔性梁模型表征;（2）大扁平比胎側(cè)曲梁的解析剛度與胎側(cè)的幾何、結(jié)構(gòu)和充氣壓力參數(shù)直接相關(guān)，可用于快速指導(dǎo)輪胎胎側(cè)的結(jié)構(gòu)設(shè)計和振動特性分析;（3）輪胎充氣壓力影響柔性胎體梁的軸向預(yù)緊力和胎側(cè)的弦剛度，進(jìn)而影響輪胎彎曲振動特性，充氣壓力越高，輪胎的諧波彎曲振動模態(tài)頻率增加，而模態(tài)阻尼降低。

參考文獻(xiàn)：

[1]劉志浩，高欽和，劉準(zhǔn)，等.重載輪胎面內(nèi)剛?cè)狁詈辖＜皞鬟f特性分析 [J] .兵工學(xué)報 ， 2018 ， 39 （ 2 ） ： 224233 .

LIU Zhi-hao， GAO Qin-he， LIU Zhun， et al. In-planerigidelastic coupling dynamic modeling and vibration response prediction of heavy duty radial tire [ J] . Acta Armamentarii， 2018 ， 39（ 2） ： 224-233 .

[2]Zhu T. Li B. Implementation of an extended model formultiaxle articulated vehicle with nonlinear tire model[Jl. SAE Technical Paper Series. 2017. DOI：10.4271/2017-01-0436.

[3]Maleska Markus . Petry Frank， Fehr Duc ， et al. Longitudinal hydroplaning performance of passenger car tires[ J] . Vehicle System Dynamics ， 2021 ， 59（ 3） ： 415432.

[4]Wei C. Olatunbosun O A. The effects of tyre materialand structure properties on relaxation length using finiteelement method [Jl. Materials& Design， 2016， 102：1420 .

[5]趙又群，一種非充氣機械彈性安全車輪的研究進(jìn)展[J].機械工程學(xué)報 ，2019， 55（ 24） ： 105-116.

ZHAO Youqun. Research progress on kind of non-pneu-matic mechanical elastic safety wheel[ J] . Journal of Mechanical Engineering ， 2019 ， 55（ 24） ： 105-116.

[6]TomiatiNicolo ， ColomboAlessandro . MagnaniGianan-tonio. A nonlinear model of bicycle shimmy [ J] . VehicleSystem Dynamics ， 2019 ， 57（ 3） ： 315335.

[7]Pacejka H B. Tire and Vehicle Dynamics L M ] . Elsevier，2005.

[8]Yu Xudong， Huang Haibo， Zhang Tao. A theoreticalthreedimensional ring based model for tire high-orderbending vibration [Jl. Joumal of Sound and Vibration.2019 .459 ：114820.

[9]Liu Z H， （Jao Q H. Analytical investigation on tire dynamics by rigidelastic coupled tire model with nonlinearsidewall stiffness [Jl. Journal of the Brazilian Society ofMechanical Sciences and Engineering， 2018， 40 （2） ：1-14 .

[10]Yamashita H. Matsutani Y . Sugiyama H. Longitudinaltire dynamics model for transient braking analysis ： ANCFLu（Jre tire model [ J] . Journal of Computational andNonlinear Dynamics ， 2015 ， 10（ 3） ： 031003 .

[11]Matsubara M， Tsujiuchi N， Koizumi T. et al. Naturalfrequency analysis of tire vibration using a thin cylindrical shell model [ J] . SAE Technical Paper Series ， 2015 ，2015 . DOI ：

10 .42 71/201501-2198 .

[12]Cooley C G ， Parker R G. Limitations of an inextensiblemodel for the vibration of high-speed rotating elasticrings with attached spacefixed discrete stiffnesses LJl .European Journal of Mechanics A/Solids. 2015 （ 54） ：187-197 .

[13]Lee J， Wang S， Pluymers B， et al. A modified com-plex modal testing technique for a rotating tire with aflexible ring modeI[J] . Mechanical Systems and SignalProcessing ， 2015 ， 60 ： 604618 .