曾桓濤,余 志,張 輝
(中山大學 智能交通研究中心 廣東省智能交通系統(tǒng)重點實驗室,廣東 廣州510006)
在我國,交通安全問題日益突出,公路交通中重大、特大事故屢有發(fā)生,造成交通事故的各種因素中,人的因素占的比例最大(占55% ~90%)。因此,采集、記錄和分析車輛行駛過程中的行駛姿態(tài)或駕駛員的操作過程,對于事故原因分析、事故責任認定以及駕駛員行為研究(如反應速度、操縱習慣等)都具有重要意義[1~3]。
目前,獲取駕駛員對方向盤的操縱過程,最直接的方法是采用方向盤轉(zhuǎn)角傳感器,這種傳感器造價較高,目前已應用于配備有主動控制系統(tǒng)(如ESP,VSC 等)高檔轎車中,但數(shù)據(jù)主要用于控制系統(tǒng)內(nèi)部處理,一般不易獲取,很難被用于監(jiān)控分析[4~6]。而作為目前車輛行駛監(jiān)控的主要裝置,行駛記錄儀(黑匣子)已經(jīng)得到了廣泛的應用,特別在美國、歐共體、巴西等國家和地區(qū)都在法律上強制規(guī)定特定車型安裝行駛記錄儀[7,8],在我國受到的重視也日益廣泛。隨著慣性傳感器的快速發(fā)展,加速度傳感器和角速度傳感器等慣性測量單元(inertial measurement unit,IMU)被引入到一些新型的行駛記錄儀中,其采集記錄得到的慣性數(shù)據(jù)為采集和提取駕駛員的操縱行為提供了途徑,為行駛記錄儀實現(xiàn)了新的功能。
在已知汽車模型、汽車的運動狀態(tài)(汽車響應,包括速度、加速度等)的基礎(chǔ)上,反求駕駛員的操縱輸入,屬于汽車操縱動力學“逆問題”,到目前為止,關(guān)于汽車操縱逆動力學問題的研究還處于初期階段。文獻[9,10]根據(jù)動力學模型,利用逆系統(tǒng)理論提出了汽車方向盤轉(zhuǎn)角識別的方法,該方法需要確定一些汽車參數(shù),但如車輪側(cè)偏剛度、整車繞鉛垂軸轉(zhuǎn)動慣量等參數(shù)通常難以獲得。本文基于運動學模型,行駛記錄儀的行駛速度和三軸加速度信息,提出了汽車轉(zhuǎn)向角的估計方法,并證明了其可行性,避免了對汽車參數(shù)的過度依賴,為駕駛行為采集和監(jiān)控提供了有效手段。
四輪汽車可以簡化為一個H 形車架的單軌模型,用前軸中點轉(zhuǎn)向角描述前輪轉(zhuǎn)向角如圖1 所示,由正弦定理可得[11]
ρ/sin[90° -(δ-αf)]=l/sin[δ-(αf-αr)].
圖1 前輪轉(zhuǎn)向汽車的轉(zhuǎn)向幾何模型Fig 1 Steering geometry model of front-wheel-steer vehicle
根據(jù)小角度近似可得
其中,δ 為前輪轉(zhuǎn)角,l 為汽車軸距,ρ 為轉(zhuǎn)彎曲率半徑,αf,αr分別為前后輪胎側(cè)偏角。
車輛曲率半徑
式中 vx,vy分別為汽車質(zhì)心縱向和側(cè)向車速,β 為質(zhì)心側(cè)偏角,r 為橫擺角速度。
側(cè)向加速度
將方程(2),(3)代入方程(1)得出
當汽車前后輪側(cè)偏差異較小時,可得
該模型中,l 為可知參數(shù),ay和vx可通過傳感器測得。因此,可由該模型實現(xiàn)轉(zhuǎn)角的近似估計。
該轉(zhuǎn)向角估計模型的推導基于二自由度單軌模型,并假設(shè)側(cè)偏角足夠小。實際上,由于汽車彈性部件的影響,模型的精度受側(cè)偏與側(cè)傾兩方面因素制約。
根據(jù)輪胎彈性力學,側(cè)向力表示為
式中 k 為輪胎彈性剛量,α 為側(cè)偏角。
汽車整體側(cè)向受力為
式中 m 為車質(zhì)量。
前后輪受力分別為
式中 a 為質(zhì)心到前軸距離,b 為質(zhì)心到后軸的距離。
聯(lián)立方程式(6),(8),(9)即可得到
式中 K 被定義為不足轉(zhuǎn)向梯度。根據(jù)定義,當K =0,即αf- αr= 0 時,汽車的轉(zhuǎn)向特性為中性轉(zhuǎn)向[12]。此時,式(5)所得的估計模型不受汽車側(cè)偏特性影響。但實際上,汽車一般會存在一定的轉(zhuǎn)向不足或轉(zhuǎn)向過度的情況。對于小型轎車,輪胎的側(cè)偏剛度kfkr值一般在28 000 ~80 000 N/rad 范圍內(nèi)[13],側(cè)偏剛度越大,側(cè)偏的影響越小。根據(jù)車輛質(zhì)心的位置,a/(a+b)=a/l 一般在[0.4 ~0.6]之間。對于前重車輛(a/l =0.4),αf- αr> 0 ,此時,汽車為不足轉(zhuǎn)向,輪胎側(cè)偏會導致轉(zhuǎn)向角估計過小;對于后重車輛(a/l =0.6),則相反。
由于汽車車廂是懸架結(jié)構(gòu),在轉(zhuǎn)彎的時候產(chǎn)生橫向力,會導致廂體發(fā)生側(cè)傾[15]。而行駛記錄儀一般固定于廂體底部,所以,行駛加速度的測量會受到汽車側(cè)傾的影響。
側(cè)向加速度的測量值為
式中 φ 為側(cè)傾角。
由于側(cè)傾角φ 一般較小,根據(jù)泰勒展開,方程(13)可近似為
而側(cè)傾角 φ =kφay,kφ為汽車平均側(cè)傾角增益,根據(jù)統(tǒng)計,轎車平均側(cè)傾角增益值大約為 7.00°/gn[13](約 為0.12 rad/gn)。
側(cè)向加速度的真實值為
聯(lián)立方程式(4),(5),(10),(14),可得
其中,γ 為相對測量誤差,δest為轉(zhuǎn)角的估計值,δ 為實際的理論值。
假設(shè)汽車質(zhì)量m=1 500 kg,軸距為l=2.5 m,為前重汽車((b-a)/l=0.2)。前后輪使用相同輪胎,側(cè)偏剛度kf=kr=28 000 N/rad,以使側(cè)偏的影響最大化。kφ取 7.00°/gn,代入式(10),式(15),則相對誤差與車速的關(guān)系為如圖2所示。
圖2 相對測量誤差與車速的關(guān)系Fig 2 Relation between relative measurement error and speed
由圖3 可知,當車速較低時(v <27 km/h),δest> δ,側(cè)傾的影響為估計誤差的主要因素,導致轉(zhuǎn)角的估計值偏大,相對誤差在10%以內(nèi);當車速較高時(27 km/h <v <35 km/h),δest<δ,側(cè)偏的影響為主要因素,導致轉(zhuǎn)角估計值偏小,相對誤差控制在7%以內(nèi)。根據(jù)方程式(12),(15)可知,當輪胎的側(cè)偏剛度越大時,估計誤差將越小,重現(xiàn)的轉(zhuǎn)向角可信度越高。
使用Mechanical Simulation Corporation 公司基于密歇根大學公路交通研究所(UMTRJ)多年的試驗與研究開發(fā)的整車動力學仿真軟件 CarSim 建立仿真模型[14,15]。CarSim 可以通過參數(shù)化建模建立仿真汽車模型、道路模型、傳感器類型和布局等,以駕駛員的角度設(shè)計控制策略,包括人車路閉環(huán)控制,以轉(zhuǎn)角和車速為控制變量的開環(huán)控制等,可以模擬穩(wěn)態(tài)行駛、或極限工況等各種場景,主要用來預測和仿真汽車整車的操縱穩(wěn)定性、制動性、平順性、動力性和經(jīng)濟性。
實驗中,選取 C 級車型 C-Class Hatchback,其中,車質(zhì)量為1 274 kg,軸距2.578 m,于后軸上方布置加速度傳感器組,使汽車以50 km/h 的目標時速進行雙變線(躲避障礙物)的開環(huán)控制,如圖3 所示。并以100 Hz 采樣,記錄前輪轉(zhuǎn)角、側(cè)向加速度、行駛速度等信息。
圖3 CarSim 仿真場景Fig 3 CarSim simulation scene
將CarSim 仿真輸出汽車后軸處側(cè)向加速度ay和行駛速度v,代入方程(5)。將求得的轉(zhuǎn)向角與實際轉(zhuǎn)向角比較,如圖4,實線為實際仿真輸出轉(zhuǎn)向角,虛線為通過加速度間接求得的轉(zhuǎn)向角,兩條曲線基本吻合,在轉(zhuǎn)彎的過程(第3 s ~第12 s),兩者平均誤差為0.002 5°。由于估計方法假設(shè)前后輪側(cè)偏差異較小,當在側(cè)向加速度較大時,側(cè)偏特性影響較為嚴重,從圖中可見,重現(xiàn)的轉(zhuǎn)角與實際值之間存在一定偏差,據(jù)統(tǒng)計,駐點處,兩者最大絕對誤差為0.18°,最大相對誤差為5.7%。
仿真實驗中,車質(zhì)量m 為1 274 kg,軸距為2.578 m,輪胎側(cè)偏剛度為kf=kr=48 700 N/rad,質(zhì)心到前軸距離a =1.016 m,將參數(shù)代入式(15),則相對測量誤差與車速的關(guān)系為如圖5 所示。
圖5 仿真實驗相對測量誤差與車速的關(guān)系Fig 5 Relation between relative measurement error and speed
可見,當輪胎側(cè)偏剛度較大時,側(cè)偏特性影響較輕微,即使在高速操縱(50 km/h)的情形下,重現(xiàn)的轉(zhuǎn)角仍有較高的可信度(相對誤差<6%)。
實驗采用一部前輪轉(zhuǎn)向的SUV 汽車,車質(zhì)量約1.4 t,軸距2.510 m,將配置有GPS、加速度傳感器等傳感器的行駛記錄儀固定于車廂后軸上方位置,其中,GPS 采用GStar GS-87 模塊,數(shù)據(jù)采樣率為 1 Hz,加速度傳感器采用ADXR300 模塊,數(shù)據(jù)采樣率1 000 Hz。在車輛運動過程中同步記錄傳感器數(shù)據(jù),以一段急速右轉(zhuǎn)過程為例,側(cè)向加速度數(shù)據(jù)經(jīng)濾波降噪后,與從GPS 數(shù)據(jù)中獲取的汽車行駛速度v 共同代入測量方程(5),得到前輪轉(zhuǎn)角的估計值,如圖6所示。
圖6 前輪轉(zhuǎn)角的估計值隨時間的變化Fig 6 Estimated value change of steering angle of front wheel with time
在汽車實際的行駛過程中,由于前輪轉(zhuǎn)角一般難以直接獲取,所以,考慮采用仿真還原驗證,以證明估計結(jié)果的可靠性。根據(jù)實驗汽車的參數(shù),在CarSim 中建立相近的汽車模型,將前輪轉(zhuǎn)角的估算結(jié)果和記錄的車速作為仿真操縱模型的輸入,進行開環(huán)控制,還原汽車的行駛過程。
圖7 所示為實驗汽車在如上所述右轉(zhuǎn)過程的行駛軌跡。其中,實線表示GPS 采集的軌跡,虛線為將前輪轉(zhuǎn)角估計值輸入CarSim 仿真還原得到的行駛軌跡。在轉(zhuǎn)彎過程中,平均時速為30 km/h(約8.33 m/s),兩條軌跡最大誤差為4.3 m。實驗證明:估計得出的汽車轉(zhuǎn)向角能夠一定程度上反映真實的操縱過程。
圖7 行駛軌跡Fig 7 Trajectory
根據(jù)汽車運動學模型,提出了基于汽車行駛速度和加速度等慣性數(shù)據(jù)的轉(zhuǎn)向角計算方法。測試結(jié)果的精度主要受汽車側(cè)偏和側(cè)傾2 種現(xiàn)象影響,經(jīng)分析可知,在一般的行駛工況下(時速35 km/h),2 種因素的綜合影響在10%以內(nèi),并且對于輪胎剛度較大的車輛(如重型車、貨車等),精度更高。通過仿真和真車等實驗,證明通過運動學模型解算汽車轉(zhuǎn)向角基本能夠還原駕駛員的操縱過程。這種方法的優(yōu)點在于計算簡便,無需涉及汽車不確定參數(shù),而且無需對汽車原有結(jié)構(gòu)和部件進行改造,適用于汽車后裝信息采集系統(tǒng)的應用(如行駛記錄儀),可以為事故分析、駕駛操縱分析、汽車狀態(tài)估計等相關(guān)應用提供了簡便可靠的途徑。參考文獻:
[1] 郭孜政.駕駛行為險態(tài)辨識理論與方法[D].成都:西南交通大學,2009.
[2] Brown I D.Drivers’margins of safety considered as a focus for research on error[J].Ergonomics,1990,33(10 - 11):1307 -1314.
[3] 鄒鐵方,余 志,蔡 銘,等.人車碰撞事故再現(xiàn)仿真結(jié)果不確定性問題研究[J].公路交通科技,2010,27(7):110 -3.
[4] 楊 財,周艷霞.方向盤轉(zhuǎn)角傳感器研究進展[J].傳感器與微系統(tǒng),2007,26(11):1 -4.
[5] Noguchi A,Yamawaki K,Yamamoto T,et al.Development of a steering angle and torque sensor of contact-type[J].Furukawa Review,2004,25(1):36 - 41.
[6] Gao J,Petovello M,Cannon M E.Integration of steering angle sensor with global positioning system and micro-electro-mechanical systems inertial measurement unit for vehicular positioning[J].Journal of Intelligent Transportation Systems,2008,12(4):159 -167.
[7] 陸文昌.汽車行駛記錄儀的現(xiàn)狀與發(fā)展趨勢[J].汽車電器,2005,5(1):1 -2.
[8] Zhao W M.New generation of vehicle traveling data recorder based on data collection system chip[J].Auto Electric Parts,2008,12(1):57 - 61.
[9] 趙又群,尹 浩,張麗霞,等.汽車操縱逆動力學的現(xiàn)狀與發(fā)展[J].中國機械工程,2005,16(1):77 -82.
[10] 張麗霞,潘福全.基于逆系統(tǒng)方法的汽車方向盤轉(zhuǎn)角識別[J].農(nóng)業(yè)機械學報,2009,40(8):40 -43.
[11] Manfred M H W.汽車動力學[M].北京:清華大學出版社,2009.
[12] Gillespic T D.車輛動力學基礎(chǔ)[M].北京:清華大學出版社,2008.
[13] 余志生.汽車理論[M].北京:機械工業(yè)出版社,2009.
[14] 李志魁.基于CarSim 的整車動力學建模與操縱穩(wěn)定性仿真分析[D].長春:吉林大學,2007.
[15] Kinjawadekar T,Dixit N,Heydinger G J,et al.Vehicle dynamics modeling and validation of the 2003 Ford Expedition with ESC using CarSim[R].SAE Technical Paper,2009.