段 斌，單穎春，劉獻棟，張紅衛(wèi)，董金松

(1.北京航空航天大學(xué) 交通科學(xué)與工程學(xué)院，北京 100191；2.交通運輸部公路科學(xué)研究院，北京 100088)

0 引 言

半掛汽車列車是物流行業(yè)主要的運輸工具，如在歐盟陸路范圍內(nèi)，74%的貨物運輸是由半掛汽車列車來完成的，并且汽車列車總長度可達18.5 m[1]。西方國家很早就開始進行甩掛運輸研究，1970年 P. M. Leucht[2]對半掛汽車列車彎道制動進行了研究；在仿真分析方面，C. C. Macadam[3]提出了一種最優(yōu)預(yù)瞄控制的駕駛員模型，并進行車輛閉環(huán)駕駛仿真；H. Peng[4]提出一種最優(yōu)預(yù)瞄控制算法；TruckSim工程軟件引入了這一模型，根據(jù)最優(yōu)控制理論建立了車輛仿真方向盤控制模塊和駕駛員閉環(huán)控制模塊[5-6]。

通過性是半掛汽車列車的一個重要性能，通過性的指標(biāo)包括最小轉(zhuǎn)彎半徑、轉(zhuǎn)彎通道寬度、外擺值等[7-9]，它們直接影響著汽車列車行駛的通過性和安全性，牛建峰，等[10]研究了山區(qū)公路緊急避險車道駛?cè)虢茄芯?。組成汽車列車的牽引車和半掛車是兩個牌照，用戶可自由選擇組成列車[11]。外擺值過大會影響半掛汽車列車的轉(zhuǎn)彎機動性能，尤其是在狹小的空間，甚至?xí)秩肫渌嚨溃绊懶熊嚢踩?。GB 1589—2004《道路車輛外廓尺寸、軸荷及質(zhì)量限值》(以下簡稱《限值》)對汽車列車外擺值的限值進行了明確的規(guī)定[8]，準(zhǔn)確計算汽車列車的外擺值對于用戶選擇半掛汽車列車匹配具有重要的參考意義，對汽車生產(chǎn)企業(yè)和汽車行業(yè)的發(fā)展都可以提供重要的技術(shù)支持。

目前國內(nèi)有兩個國家標(biāo)準(zhǔn)對外擺值做出規(guī)定。在《限值》中，規(guī)定汽車和汽車列車(不計具有作業(yè)功能的專用裝置的突出部分)必須能在同一個車輛通道圓內(nèi)通過，車輛通道圓的外圓直徑D1為25.00 m，車輛通道圓的內(nèi)圓直徑D2為10.60 m。汽車和汽車列車由直線行駛過渡到上述圓周運動時，任何部分超出直線行駛時的車輛外側(cè)面垂直面的值(車輛外擺值)T不得大于0.80 m(圖1)。

圖1 GB 1589—2004測量方法示意Fig.1 Measurement in GB 1589—2004

在《限值》和GB/T 12540—2009《汽車最小轉(zhuǎn)彎直徑、最小轉(zhuǎn)彎通道圓直徑和外擺值測量方法》(以下簡稱《方法》)中均規(guī)定了汽車列車外擺值的試驗方法。不同的是，在《限值》中規(guī)定汽車列車起步后由直線行駛過渡到直徑為25 m的圓上，而在《方法》中規(guī)定汽車列車以保持預(yù)定轉(zhuǎn)角狀態(tài)在圓內(nèi)行駛[9]。歐洲標(biāo)準(zhǔn)70/311/EEC《關(guān)于各成員國汽車及其掛車轉(zhuǎn)向機構(gòu)的相關(guān)法律》有不同的規(guī)定[12]。經(jīng)文獻檢索，目前尚未有明確的外擺值數(shù)值計算方法。

筆者基于半掛汽車列車在轉(zhuǎn)彎通道圓內(nèi)轉(zhuǎn)彎過程的作圖法，在《方法》規(guī)定的實驗方法的基礎(chǔ)上，提出了根據(jù)汽車列車的外形尺寸確定外擺值的理論計算方法，推導(dǎo)出具體的計算公式；之后，研究了汽車列車不同的參數(shù)對外擺值的影響關(guān)系；筆者還利用TruckSim軟件對汽車列車進行運動學(xué)建模和仿真分析，將仿真計算得到的外擺值和文中算法得到外擺值進行比較，驗證了文中所提出的計算方法的正確性。實車測試的結(jié)果也驗證了文中算法的正確性。筆者工作對汽車列車的匹配研究提供了重要的參考依據(jù)。

1 半掛汽車列車外擺值計算方法

建立慣性坐標(biāo)系XOY(圖2)，設(shè)半掛車最后外點為A點，后軸中心線中點為B點，牽引銷位置為o點，半掛車縱向中心線與X軸夾角為α，于是A點在坐標(biāo)XOY中可以表示如式(1)、式(2)。

XA=Xo+ cos (α)(-L)-sin(α)M

(1)

YA=Yo+ sin(α)(-L)+ cos (α)M

(2)

式中：L，M分別是A點在半掛車上相對于牽引銷的縱向和橫向的距離。

圖2 汽車列車轉(zhuǎn)向分析Fig.2 Analysis of combination vehicles’ turn

計算原理：將轉(zhuǎn)彎過程分成dc次，每次引銷的位置o的軌跡繞著轉(zhuǎn)彎中心轉(zhuǎn)dd°，根據(jù)《限值》要求，至少在圓周內(nèi)行駛1/2圈(半個圓周)，筆者在計算過程中dc和dd滿足關(guān)系:dd×dc=π。

根據(jù)汽車列車轉(zhuǎn)彎過程的作圖方法可知[13]，轉(zhuǎn)彎時B點在i+1時刻位置的確定方法如下：以i+1時刻牽引銷的位置oi+1點為圓心、半掛車軸距L2為半徑畫一個圓弧，該圓弧與i時刻掛車縱向中心線的交點即近似為i+1時刻點B的位置，進而可知i+1時刻半掛車縱向中心線與X軸的夾角αi+1，再根據(jù)式(1)和式(2)確定i+1時刻A點的位置。按照這種方法持續(xù)進行下去，即可得到A點的軌跡，進一步即可求出外擺值。

基于本計算原理，具體的實現(xiàn)過程如下：設(shè)第i時刻牽引銷o點在坐標(biāo)系XOY中的坐標(biāo)為(Xoi，Yoi)，則第i時刻半掛車后軸線中點的坐標(biāo)為：

XBi=Xoi+ cos (αi)(-L)

(3)

YBi=Yoi+ sin (αi) (-L)

(4)

此時掛車縱向中心線的方程為：

YB(i+1)=(XB(i+1)-Xoi)tan (αi)+Yoi

(5)

在第i+1時刻，牽引銷的坐標(biāo)為(Xo(i+1)，Yo(i+1))，那么以第i+1時刻牽引銷的位置點為圓心、半掛車軸距L2為半徑的圓的方程為：

(6)

式(5)、式(6)兩個方程聯(lián)合求解，即可得到i+1時刻B點的位置坐標(biāo)XB(i+1)，YB(i+1)。連接o(i+1)，B(i+1)，即可得到i+1時刻半掛車縱向軸線；其與X軸夾角用式(7)計算：

(7)

根據(jù)式(1)、式(2)可得到A點的坐標(biāo)為：

XA(i+1)=Xo(i+1)+cos (α(i+1)) (-L) -

sin (α(i+1))M

(8)

YA(i+1)=Xo(i+1)+sin (α(i+1)) (-L) -

cos (α(i+1))M

(9)

重復(fù)以上過程，就可以逐步得到A點在i+2，i+3，…時刻的坐標(biāo)，進而得到A點的軌跡，在此基礎(chǔ)上便可求出外擺值。

在上述的迭代計算過程中，首先要確定初始狀態(tài)的情況，如圖3。

圖3 汽車列車轉(zhuǎn)向初始狀態(tài)Fig.3 Initial state of combination vehicles’ turn

按照《方法》，初始狀態(tài)下汽車列車后軸中線與坐標(biāo)系XOY的X軸重合，前輪處于合適轉(zhuǎn)向位置以保證在轉(zhuǎn)向過程中牽引車左前最外側(cè)正好沿著直徑為25 m的圓弧，在轉(zhuǎn)彎過程中，轉(zhuǎn)向輪的轉(zhuǎn)向角保持不變。初始狀態(tài)下，坐標(biāo)系XOY的O點與牽引銷o點的連線與X軸的夾角為β0，β0可通過汽車列車的外廓尺寸計算得出。牽引銷o點到轉(zhuǎn)向中心O點的距離S也可計算得出，且在列車轉(zhuǎn)彎過程中保持不變。算法流程如圖4。

圖4 算法流程Fig.4 Algorithm flow chart

2 外擺值實例計算分析

2.1 利用筆者提出的算法進行計算

以某牽引車和某半掛車進行匹配，在中國人民解放軍總裝備部定遠試驗場進行了實際測試。筆者對實驗車輛使用文中的計算方法進行分析。算例中應(yīng)用的具體數(shù)據(jù)如表1。

表1 計算數(shù)據(jù)Table 1 Parameters of combination vehicles /m

根據(jù)以上理論的分析和具體數(shù)據(jù)，進行編程計算得到汽車列車的外擺值為0.275 1 m，這個計算結(jié)果是符合《限值》對外擺值的限值要求的(注：此計算結(jié)果是按照每次0.1°迭代計算得到的)。

文中計算方法是根據(jù)作圖法推導(dǎo)出來的，按照不同的迭代步長計算得到的結(jié)果不一樣。由表2可以看出，隨著迭代步長的減小，計算得出的外擺值也在不斷減小?？上攵?，當(dāng)?shù)介L無限小的時候，得到的計算結(jié)果就可以和真實解無限接近。筆者采用0.1°迭代步長。

表2不同迭代步長下的外擺值

Table2Outvaluesofdifferentiterationsteplengths

迭代步長/(°)外擺值/m迭代步長/(°)外擺值/m300．589530．3279200．518910．2914100．44870．10．275150．3614

2.2 外擺值隨參數(shù)變化的規(guī)律研究

在歐美發(fā)達國家，以牽引車拖帶掛車組成的半掛汽車列車的運輸量占到總運輸量的近80%。目前我國牽引車與掛車比達到4∶5[14]，但與世界發(fā)達國家1∶2.5[15]以上的標(biāo)準(zhǔn)仍有較大差距。牽引車和半掛車的參數(shù)設(shè)置對半掛汽車列車的通過性都有影響。

為了分析前置距變化對外擺值的影響，對前置距E從0.8 m變化到1.5 m(步長為0.1 m)分別進行外擺值計算,得到外擺值隨前置距E的變化曲線。同樣，對牽引車的軸距L1、半掛車軸距L2取不同的數(shù)進行計算并得到相應(yīng)的變化曲線。外擺值變化曲線見圖5。

圖5 外擺值變化曲線Fig.5 Curve of out value change

由圖5可以看出，汽車列車的外擺值隨著前置距的增大，呈現(xiàn)增大的趨勢；隨著牽引車軸距的增大，呈現(xiàn)單調(diào)增大的趨勢；而隨著半掛車軸距L2的增大，汽車列車的外擺值呈現(xiàn)減小的趨勢；隨半掛車后輪中心距B2的增加也呈現(xiàn)減小的趨勢。前置距E、牽引車軸距L1、半掛車后輪中心距B2等對外擺值影響較小，而半掛車軸距L2的影響相對較大。同時，外擺值在參數(shù)變化范圍內(nèi)均遠小于0.8 m的國標(biāo)要求。

3 外擺值計算方法的驗證

3.1 使用TruckSim工程軟件進行仿真驗證

將半掛汽車列車實際數(shù)據(jù)在TruckSim仿真軟件中輸入進行建模，根據(jù)《限值》給定的汽車列車運行軌跡，并在模型中安裝測試半掛汽車列車最后外點的位移傳感器，可根據(jù)其位移曲線計算外擺值的數(shù)值。根據(jù)TruckSim的仿真動畫和曲線也可以明顯地看出汽車列車在轉(zhuǎn)彎時的運動趨勢。圖6為TruckSim仿真運行的動畫截圖。圖7為TruckSim仿真運行計算得到的汽車列車最后外點的軌跡。

圖6 TruckSim仿真運行的動畫截圖Fig.6 Screenshot of animation

圖7 汽車列車最后外點的軌跡Fig.7 Track of last and outermost point of combination vehicle

筆者計算汽車列車外擺值是依據(jù)軌跡側(cè)向坐標(biāo)數(shù)據(jù)，在圖7中也可以看出軌跡側(cè)向坐標(biāo)從0時刻起有一個先下降再上升的過程，從而可以依據(jù)曲線的坐標(biāo)計算外擺值。

在TruckSim中，將仿真得到的數(shù)據(jù)(圖7)導(dǎo)出，并將汽車列車最后外點外擺最大時刻對應(yīng)的坐標(biāo)值與0時刻對應(yīng)的坐標(biāo)值相減，即可得到汽車列車的外擺值。將表1的實驗車輛數(shù)據(jù)輸入模型，通過仿真得到外擺值的計算結(jié)果為0.307 361 m，通過文中算法計算得到的外擺值是0.275 1 m，誤差較小，算法正確。

通過改變汽車列車某參數(shù)，可以得到改變參數(shù)后的一系列外擺值數(shù)據(jù)(表3～表5)，由此可以畫圖觀察外擺值隨該參數(shù)的變化趨勢(圖8)。筆者針對汽車列車外擺值E、牽引車軸距L1、半掛車軸距L2的變化進行仿真，得到了外擺值隨這些參數(shù)變化的曲線，并與通過筆者提出的算法計算得到的變化曲線進行比較，驗證算法的正確性。

表3 不同的前置距E對應(yīng)的外擺值Table 3 Out values of different fifth wheel leads /m

表4 不同的牽引車軸距L1對應(yīng)的外擺值Table 4 Out values of different tractor wheel bases /m

表5 不同的半掛車軸距L2對應(yīng)的外擺值Table 5 Out values of different semitrailer axle bases /m

圖8 外擺值變化曲線對比Fig.8 Comparison of out value change curves

由圖8可以看出，通過計算和仿真得到的汽車列車外擺值隨前置距E、牽引車軸距L1、半掛車軸距L2的變化趨勢是相同的，且文中方法所得結(jié)果與TruckSim仿真結(jié)果非常相近，可以看出筆者提出的外擺值算法得到了很好的驗證。

3.2 實車測試驗證

利用南京智真電子科技有限公司開發(fā)的基于視覺的汽車列車行駛軌跡測試儀，將文中進行計算和仿真的半掛汽車列車在中國人民解放軍總裝備部定遠試驗場進行了實際測試。

本方法先在地面上繪制5條標(biāo)識線，同時隔一定的距離繪制一條垂直于標(biāo)識線的橫向線，如圖9。通過觸發(fā)使攝像機開始采集形成視頻，并通過圖像上橫向線與標(biāo)識線的交點，求出圖像中點在地面上的實際位置。

圖9 汽車列車行駛軌跡測量鋪設(shè)參照圖像Fig.9 Reference for measuring traveling track

按照《限值》規(guī)定的測試方法，利用該測試儀對半掛汽車列車的原始狀態(tài)和最后外側(cè)點的軌跡進行記錄，通過計算記錄的數(shù)據(jù)得到外擺值，如圖10。定遠試驗場進行的半掛汽車列車外擺值測試結(jié)果是0.285 m，通過文中算法計算得到的外擺值是0.275 1 m，誤差為3.5%，誤差很小，驗證了本算法的正確性。

圖10 汽車列車行駛軌跡測量試驗Fig.10 Test of measuring traveling track

4 結(jié) 論

筆者提出了根據(jù)汽車列車的外形尺寸確定外擺值的計算方法，針對具體的實驗車輛，利用該方法進行了外擺值的計算與分析，且通過TruckSim工程軟件仿真和實車測試對該方法進行了驗證，表明筆者所提計算方法結(jié)果準(zhǔn)確、公式簡潔、便于理解和工程實用，可為半掛汽車列車的通過性研究提供指導(dǎo)。經(jīng)過筆者的計算、分析，可得出以下結(jié)論：

1)半掛汽車列車外擺值隨前置距E、牽引車軸距L1、半掛車軸距L2、半掛車后輪中心距B2的變化趨勢是，外擺值隨前置距和牽引車軸距的增大而增大，隨半掛車軸距、半掛車后輪中心距的增大而減小。

2)前置距和半掛車后輪距對半掛汽車列車外擺值的影響很小，牽引車軸距和半掛車軸距對半掛汽車列車外擺值的影響較大。

[1] European Commission.Road Transport Policy[R].Brussels:Commission of the European Communities,2006:2-3.

[2] Leucht P M.The directional dynamics of the commercial tractor-semitrailer vehicle during braking[J].TRB of the national academies,1970,71(3):242-254.

[3] Macadam C C.Application of an optimal preview control for simulation of closed-loop automobile driving[J].IEEE Transactions on Systems,Man and Cybernetics,1981,11(6):393-399 .

[4] Peng H.Optimal Preview Control for Vehicle Lateral Guidance[R].California,U.S.:Institute of Transportation Studies University of California at Berkeley,1991.

[5] Carsim.The vehiclesim steer controller[EB/OL].(2012-10-25) www.carsim.com.

[6] Dunn A L,Heydinger G,Rizzoni G,et al.In-depth analysis of the influence of high torque brakes on the jackknife stability of heavy trucks[J].Society of Automotive Engineers,2003,112(2):255-272 .

[7] 郭正康.現(xiàn)代汽車列車設(shè)計與使用[M] .北京:北京理工大學(xué)出版社,2006:114-116.

Guo Zhengkang.Design and Application of Modern Vehicle Combinations[M].Beijing:Beijing Institute of Technology Press,2006:114-116.

[8] GB 1589—2004 道路車輛外廓尺寸、軸荷及質(zhì)量限值[S].北京:中國標(biāo)準(zhǔn)出版社,2004.

GB 1589—2004 Limits of Dimensions,Axle Load and Masses for Road Vehicles[S].Beijing:China Standards Press,2004.

[9] GB 12540—2009 汽車最小轉(zhuǎn)彎直徑、最小轉(zhuǎn)彎通道圓直徑和外擺值測量方法[S].北京:中國標(biāo)準(zhǔn)出版社,2009.

GB 12540—2009 Minimum Turning Circle Diameter,Minimum Turning Clearance Circle Diameter and Out Value Test Method for Motor Vehicles[S].Beijing:China Standards Press,2009.

[10] 牛建峰,王俊驊,董憲元.山區(qū)公路緊急避險車道駛?cè)虢茄芯縖J].重慶交通大學(xué)學(xué)報:自然科學(xué)版,2010,29(4) :604-608.

Niu Jianfeng,Wang Junhua,Dong Xianyuan.Research on the angle to truck escape ramps of mountain roads[J].Journal of Chongqing Jiaotong University:Natural Sience,2010,29(4):604-608.

[11] 王增權(quán),趙希望.論汽車列車的合理匹配 [J].商用汽車,2010(3):82-84.

Wang Zengquan,Zhao Xiwang.A thought and consideration over reasonable and ideal connection of tractor and trailer[J].Commercial Vehicles,2010(3):82-84.

[12] 70/311/EEC On the Approximation of the Laws of the Member States Relating to the Steering Equipment for Motor Vehicles and Their Trailers[S].[S.l.]:EEC Commission,1999.

[13] 何光里.汽車運用工程師手冊[M].北京:人民交通出版社,1997:112-113.

He Guangli.Automobile Application Engineer’s Handbook[M].Beijing:China Communications Press,1997:112-113.

[14] 張紅衛(wèi),王維.貨運車輛安全研究與應(yīng)用[J].物流技術(shù)與應(yīng)用:貨運車輛,2009(5):44-47.

Zhang Hongwei,Wang Wei.Freight vehicle safety research and application[J].Logistics Technology and Application:Freight Vehicles,2009(5):44-47.

[15] 羅茂堂.淺議甩掛運輸行業(yè)的科學(xué)發(fā)展[J].科技創(chuàng)新導(dǎo)報,2008(32):216-217.

Luo Maotang.Discussion on the scientific development of hang dump the transportation industry[J].Science and Technology Innovation Herald,2008(32):216-217.