周志國 方棟華 方劍烽 喻學(xué)才

(浙江交通職業(yè)技術(shù)學(xué)院汽車學(xué)院1) 杭州 311112) (浙江師范大學(xué)工學(xué)院2) 金華 321004)

基于轉(zhuǎn)動慣量轉(zhuǎn)移的農(nóng)用冷藏物流車動態(tài)穩(wěn)定特性研究及控制仿真*

周志國1)方棟華1)方劍烽1)喻學(xué)才2)

(浙江交通職業(yè)技術(shù)學(xué)院汽車學(xué)院1)杭州 311112) (浙江師范大學(xué)工學(xué)院2)金華 321004)

基于整體運動學(xué)模型，建立了農(nóng)用冷藏物流車模型及其狀態(tài)方程，并對其動態(tài)特性進行數(shù)學(xué)分析，運用整車轉(zhuǎn)動慣量辨識分析算法，對車載貨物重心位置變化引起的整車轉(zhuǎn)動慣量進行辨識評估；基于MATLAB平臺對車輛動態(tài)穩(wěn)定性能進行時域和頻域仿真分析.仿真結(jié)果表明，該方法可以得到車載貨物重心位置變化對車輛動態(tài)穩(wěn)定特性的影響結(jié)果.

現(xiàn)代物流；主動安全；模型分析法；動態(tài)穩(wěn)定特性；計算機仿真

0 引 言

現(xiàn)代高速農(nóng)用物流冷藏物車因為車載貨物重心位置變化容易引起整車轉(zhuǎn)動慣量轉(zhuǎn)移，進而給車輛動態(tài)穩(wěn)定特性造成嚴(yán)重影響.文獻(xiàn)[1-3]提出需要針對冷藏物流運輸車的特殊需求進行新的結(jié)構(gòu)設(shè)計和優(yōu)化.由于箱式物流運輸車車速的提高對車輛的穩(wěn)定控制系統(tǒng)也提出了更高的要求，需要在車輛轉(zhuǎn)向系統(tǒng)上有新的技術(shù)要求[4-5].同時，隨著車載電控系統(tǒng)技術(shù)的發(fā)展，對物流運輸車穩(wěn)定性電控系統(tǒng)進行設(shè)計和開發(fā)也是研究的一個方向[6-7].

但是，車輛運行速度作為影響農(nóng)用冷藏物流運輸車效率的主導(dǎo)因素之一，而冷藏物流車車載貨物的重量和重量分布情況具有很大的不確定性，這種不確定性會引起車輛的重心和轉(zhuǎn)動慣量發(fā)生變化，從而對車輛的動態(tài)穩(wěn)定性產(chǎn)生較大的影響，目前成為又一個新的研究重點和難點.因此，為了降低車載貨物對車輛穩(wěn)定性的影響，文中提出基于轉(zhuǎn)動慣量轉(zhuǎn)移對冷藏物流車動態(tài)穩(wěn)定特性進行研究.

1 車輛模型

選取應(yīng)用較多的4×2箱式冷藏物流運輸車作為研究對象，為了研究車載貨物重心轉(zhuǎn)移對車輛轉(zhuǎn)動慣量的影響，分析轉(zhuǎn)動慣量大小和位置的變化對整車動態(tài)穩(wěn)定性的影響，建立車輛整車運動學(xué)模型.

(1)

式中：mveh為整車質(zhì)量；MZi為整車轉(zhuǎn)動慣量矩；ωr為車聲橫擺角速度；FXi為車輛縱向作用分力；FYi為車輛側(cè)向作用分力；Iz為橫擺轉(zhuǎn)動慣量.

為了便于分析車輛重心位置轉(zhuǎn)移與整車轉(zhuǎn)動慣量變化之間的關(guān)系，以及引起的車輛動態(tài)穩(wěn)定特性的變化，忽略車輛地面切向力對輪胎側(cè)偏特性的影響，忽略空氣動力對車輛的作用，忽略車輛由于載荷變化引起的輪胎動力學(xué)特性的影響以及車輛回正力矩的影響等，對方程(1)中的車輛模型進行簡化得到2自由度單軌車輛模型，見圖1.

圖1 單軌車輛模型

由圖1可知2自由度單軌車輛模型為

(2)

式中：v為車輛側(cè)向速度；u為車輛前進速度；a為前軸到質(zhì)心距離；b為后軸到質(zhì)心距離；δf為前輪轉(zhuǎn)角；kaf為前輪側(cè)偏剛度；kar為后輪側(cè)偏剛度.

動態(tài)系統(tǒng)狀態(tài)方程是對系統(tǒng)進行時域和頻域分析的經(jīng)典分析和描述方式，令

(3)

則方程(2)可轉(zhuǎn)換為狀態(tài)方程，即

Y=CX+DuIn

(4)

式中：

；

2 整車轉(zhuǎn)動慣量計算方法

車輛的轉(zhuǎn)動慣量對整車的操縱穩(wěn)定性、制動性和行駛平順性有很大的影響，對于一般剛體，根據(jù)理論力學(xué)的原理，其對某一軸線XG的轉(zhuǎn)動慣量就是剛體中各個質(zhì)點質(zhì)量mGi與其到該軸的距離rGi的平方之乘積的累積值，其計算式為

(5)

在計算汽車的轉(zhuǎn)動慣量時，理論上可以定義車輛是一個連續(xù)的剛體，但是，汽車本身又是一個不均勻的、空間分布比較復(fù)雜的結(jié)構(gòu)體，因此，很難采用上述理論方法積分得到車輛的轉(zhuǎn)動慣量.但是，對應(yīng)設(shè)計好的車輛，其轉(zhuǎn)動慣量也隨之確定了，因此，工業(yè)上通常采用實驗測試法和經(jīng)驗公式法來確定車輛的轉(zhuǎn)動慣量.

這里重點討論經(jīng)驗公式法[8-9]，這是一種建立在大量實驗數(shù)據(jù)基礎(chǔ)之上計算整車轉(zhuǎn)動慣量的方法，方法如下.

(6)

IXG=KmL2

(7)

(8)

式中：K和K1為經(jīng)驗系數(shù).

式(6)只有在很特殊的情況下才可以使用，針對冷藏物流運輸車的載荷分布情況，見圖2，這里選取式(7)作為計算基礎(chǔ)經(jīng)驗公式.

圖2 農(nóng)用冷藏物流車載荷分布示意圖

由圖2可知，由于對于自重引起的車輛轉(zhuǎn)動慣量和質(zhì)心位置在運行過程中變化比較小，可以采用已測試數(shù)值，而對于車載貨物的轉(zhuǎn)動慣量和質(zhì)心位置則會隨著實際情況的變化而發(fā)生較大變化，這也是引起整車轉(zhuǎn)動慣量和質(zhì)心位置變化的重要因素，因此，將這2部分載重量進行分離計算，則有

(9)

(10)

則可得整車轉(zhuǎn)動慣量變化值為

(11)

整車質(zhì)心位置距離前軸變化值為

Δa=a-a′

(12)

3 車輛動態(tài)穩(wěn)定性仿真分析

本節(jié)對車載貨物引起的整車轉(zhuǎn)動慣量和質(zhì)心位置變化對其穩(wěn)定性的影響進行仿真分析，參數(shù)見表1.

表1 仿真分析所用參數(shù)

3.1 前輪轉(zhuǎn)角6°時，不同車速下動態(tài)響應(yīng)仿真分析

在整車轉(zhuǎn)動慣量變化為±10%，前輪轉(zhuǎn)角6°時，分別針對10,20,30 m/s不同車速工況時對整車動態(tài)穩(wěn)定性進行仿真分析，仿真分析結(jié)果見圖3～8.

圖3 時域車輛動態(tài)響應(yīng)圖(10 m/s)

圖4 頻域車輛動態(tài)響應(yīng)圖(10 m/s)

圖5 時域車輛動態(tài)響應(yīng)圖(20 m/s)

圖6 頻域車輛動態(tài)響應(yīng)圖(20 m/s)

圖7 時域車輛動態(tài)響應(yīng)圖(30 m/s)

圖8 頻域車輛動態(tài)響應(yīng)圖(30 m/s)

從仿真結(jié)果可知，對于時域結(jié)果，低轉(zhuǎn)動慣量對輸入的響應(yīng)在回到穩(wěn)態(tài)值之前出現(xiàn)了超調(diào)，反映了其系統(tǒng)的低阻尼特性,而當(dāng)轉(zhuǎn)動慣量增大后,車輛響應(yīng)沒有超調(diào)現(xiàn)象，表現(xiàn)為高阻尼特性.

對于頻域響應(yīng)，轉(zhuǎn)動慣量減小則表現(xiàn)出之后較小，響應(yīng)較快，增大則相反.

對于車速響應(yīng)，在低速時轉(zhuǎn)動慣量增減對車輛動態(tài)性能影響不大，當(dāng)車速提高時，較大的轉(zhuǎn)動慣量對車輛動態(tài)響應(yīng)的影響增大，反之則減小.

3.2 前輪轉(zhuǎn)角15°時，不同車速下動態(tài)響應(yīng)仿真分析

在整車轉(zhuǎn)動慣量變化為±10%，前輪轉(zhuǎn)角增加為15°時，分別針10，20，30 m/s不同車速工況時對整車動態(tài)穩(wěn)定性進行仿真分析，仿真分析結(jié)果分別見圖9～14.

圖9 時域車輛動態(tài)響應(yīng)圖(10 m/s)

圖10 頻域車輛動態(tài)響應(yīng)圖(10 m/s)

圖11 時域車輛動態(tài)響應(yīng)圖(20 m/s)

圖12 頻域車輛動態(tài)響應(yīng)圖(20 m/s)

圖13 時域車輛動態(tài)響應(yīng)圖(30 m/s)

圖14 頻域車輛動態(tài)響應(yīng)圖(30 m/s)

從上述仿真結(jié)果可知，當(dāng)轉(zhuǎn)角增大到15°時，車輛轉(zhuǎn)動慣量的增減趨勢和在6°時基本一致，但車速增加后，尤其是在高速階段，車輛橫擺角速度增加比較明顯，車輛趨于失穩(wěn).

3.2 根軌跡分析

為了在最大車速范圍內(nèi)對車輛轉(zhuǎn)動慣量和質(zhì)心位置變化對車輛性能的影響進行全面的分析，對整車轉(zhuǎn)動慣量分別變化為標(biāo)準(zhǔn)轉(zhuǎn)動慣量的±10%和標(biāo)準(zhǔn)轉(zhuǎn)動慣量3種情況下最大車速范圍內(nèi)的系統(tǒng)根軌跡圖進行分析，仿真結(jié)果見圖15.

圖15 車輛系統(tǒng)根軌跡應(yīng)圖(5～30 m/s)

從上述仿真結(jié)果可知，在車速較小時特征值全落在實軸上，隨車速的不斷加大才逐漸接近虛軸，可見該車型在設(shè)計速度范圍內(nèi)是穩(wěn)定的，但是，隨著轉(zhuǎn)動慣量的變化，車輛的根軌跡距離虛軸較近，造成車輛的穩(wěn)定性儲備變小.

4 結(jié) 束 語

結(jié)合現(xiàn)代農(nóng)用冷藏物流運輸車的特點，針對其載貨重物質(zhì)量和質(zhì)心位置變化對物流運輸車的動態(tài)穩(wěn)定性的影響的問題進行了研究，建立了車輛動力學(xué)模型，并對關(guān)鍵變量進行分析，計算得到貨重物質(zhì)量和質(zhì)心位置變化對物流運輸車的動態(tài)穩(wěn)定性的數(shù)學(xué)對應(yīng)關(guān)系，并且通過仿真分析，驗證了分析的有效性，為改善和提高農(nóng)用冷藏物流運輸車的穩(wěn)定性提供了一個有效的分析方法，對優(yōu)化車輛設(shè)計，提高其安全性和效率具有重要的參考價值.

[1]王秀景,魯力群,冷藏物流車運輸車匹配設(shè)計及分析[J].物流技術(shù)與應(yīng)用，2016(5):55-58.

[2]陳琦峰.汽車物流跨界運輸車研究與設(shè)計[J].物流技術(shù)與應(yīng)用，2010(5):156-160.

[3]付智慧，趙長春.廂式運輸車新技術(shù)的應(yīng)用[J].專用汽車，2003(3):67-72.

[4]黎法豪.汽車物流智能運輸車控制系統(tǒng)的開發(fā)與實現(xiàn)[J].高新技術(shù)，2015(6):98-105.

[5]李達(dá)，張銳,楊偉明，等.農(nóng)產(chǎn)品物流運輸車遠(yuǎn)程定位監(jiān)測系統(tǒng)[J].保鮮與加工，2012,12(3):20-23.

[6]LI L. A novel vehicle dynamics stability control algorithm based on the hierarchical strategy with constrain of nonlinear tyre forces[J]. Vehicle System Dynamics,2015,53(8):1093-1116.

[7]OLEKSOWICZ S A, BURNHAM K, SOUTHGATE A, et al. Regenerative braking strategies, vehicle safety and stability control systems: critical use-case proposals[J]. Vehicle System Dynamics,2013(1):1-16.

[8]高云凱,馮海星,方劍光，等.基于質(zhì)量線法的駕駛室慣性參數(shù)識別試驗研究[J].振動與沖擊,2013,32(6):189-193.

[9]MAHMOUDIAN N, GEISBERT J, WOOLSEY C. Approximate analytical turning conditions for underwater gliders: implications for motion control and path planning[J]. IEEE Journal of Oceanic Engineering,2010，35(1):131-143.

Study on Dynamic Response and Computer Simulation for Refrigerated Logistics Vehicles Based on the Moment of Inertia Transfer

ZHOU Zhiguo1)FANG Donghua1)FANG Jianfeng1)YU Xuecai2)

(AutomobileSchoolofZhejiangInstituteofCommunications,Hangzhou311112,China)1)(SchoolofEngineering,ZhejiangNormalUniversity,Jinhua321004,China)2)

Based on the kinematic model of vehicle, the model and state equations of modern high-speed refrigerated logistic vehicle are established to analyze the dynamic stability characteristics of the vehicle. By the identification method of rotational inertia of vehicle, identification evaluation is conducted for the rotational inertia induced by the variation of the center of gravity of the goods in the vehicle. By using MATLAB, both time- and frequency-domain analyses of the vehicle’s dynamic stability characteristics are performed. The simulation results show that this method can calculate the effects of variation of center of gravity on the dynamic stability of the refrigerated logistics vehicles.

modern high-speed logistics; active safety; model analysis; dynamic stability; computer simulation

2016-10-17

*浙江省自然科學(xué)基金項目資助(LY13E080010)

TP391.9

10.3963/j.issn.2095-3844.2016.06.007

周志國(1978—)：男，碩士，副教授，主要研究領(lǐng)域為汽車智能檢測與診斷技術(shù)、車輛控制與仿真