高璐 高磊 孫禮

摘要：? 為研究汽車(chē)受側(cè)風(fēng)影響時(shí)風(fēng)壓中心線位置對(duì)整車(chē)操縱穩(wěn)定性的影響，在CCM+軟件中計(jì)算不同風(fēng)速和不同流入角下的氣動(dòng)力，獲得氣動(dòng)力與風(fēng)速和車(chē)速的合速度以及流入角的初始關(guān)系曲線。利用MATLAB和Adams/Car聯(lián)合仿真，搭建空氣動(dòng)力學(xué)和車(chē)輛動(dòng)力學(xué)的雙向耦合模型，以風(fēng)壓中心線與質(zhì)心的距離為變量進(jìn)行仿真分析。結(jié)果顯示，隨著風(fēng)壓中心線由質(zhì)心前方移動(dòng)到質(zhì)心后方：在開(kāi)環(huán)仿真工況中，車(chē)輛的側(cè)向位移會(huì)減小，但是如果風(fēng)壓中心線位于車(chē)輛后方超過(guò)一定的距離，車(chē)輛會(huì)在氣動(dòng)力產(chǎn)生的橫擺力矩作用下向另外一側(cè)偏移;閉環(huán)仿真工況中，當(dāng)風(fēng)壓中心線穿過(guò)車(chē)輛質(zhì)心時(shí)，由于側(cè)向力的存在，方向盤(pán)依然需要一個(gè)很小的轉(zhuǎn)角以維持直線行駛，風(fēng)壓中心線需要繼續(xù)后移以保證方向盤(pán)穩(wěn)態(tài)時(shí)的回正。

關(guān)鍵詞：? 氣動(dòng)性; 操縱穩(wěn)定性; 雙向耦合; 風(fēng)壓中心線; 流入角; 質(zhì)心; 側(cè)風(fēng); 聯(lián)合仿真

中圖分類(lèi)號(hào)：? U461.1;TB115.1文獻(xiàn)標(biāo)志碼：? B

收稿日期： 2020-08-04修回日期： 2021-08-31

作者簡(jiǎn)介： 高璐（1986—），女，安徽廬江人，研究方向?yàn)檐?chē)輛動(dòng)力學(xué)，（E-mail） 416431884@qq.comInfluence of aerodynamics on handling stability based on

real time change of inflow angle

GAO Lu， GAO Lei， SUN Li

（CAE Department， Chery Automobile Co.， Ltd.， Wuhu 241000， Anhui， China）

Abstract： To study the influence of the position of the wind pressure center line on the handling and stability of the vehicle under the influence of crosswind， the aerodynamic forces under different wind speeds and inflow angles are calculated in CCM+ software， and then the initial relationship curves of aerodynamic force with wind speed， vehicle speed and inflow angle are obtained. Using the co-simulation of MATLAB and Adams/Car， a bidirectional coupling model of air aerodynamics and vehicle dynamics is established， and the distance between the wind pressure center line and center of mass is taken as the variable for simulation analysis. The results show that， as the wind pressure center line moves from the front of the mass center to the back of the center， the lateral displacement of the vehicle will be reduced in open loop simulation， however， while the wind pressure center line is located behind the vehicle for more than a certain distance， the vehicle will offset to the other side under the yaw moment generated by aerodynamic force. In closed-loop simulation， while the wind pressure center line passes through the vehicle centroid， due to the existence of lateral force， the steering wheel still need a small angle to maintain straight driving， and the wind pressure center line needs to continue to move backward to ensure the alignment of the steering wheel in steady state.

Key words： aerodynamics; handling stability; bidirectional coupling; wind pressure center line; inflow angle; centroid; crosswind; co-simulation

0引言汽車(chē)技術(shù)的發(fā)展、新材料和新工藝的應(yīng)用，以及輕量化水平的提升，導(dǎo)致整車(chē)輪荷降低。為滿足顧客的個(gè)性化追求，汽車(chē)造型不得不在氣動(dòng)性方面作一些犧牲。這些改變均對(duì)整車(chē)操控產(chǎn)生不利影響。為研究氣動(dòng)性對(duì)整車(chē)操控的影響，實(shí)現(xiàn)在不影響整體造型風(fēng)格的前提下，提出有利于車(chē)輛穩(wěn)定性的建議，部分學(xué)者嘗試結(jié)合空氣動(dòng)力學(xué)和車(chē)輛動(dòng)力學(xué)研究汽車(chē)的側(cè)風(fēng)穩(wěn)定性，結(jié)果表明：側(cè)風(fēng)的大小、方向和波形等變化，均會(huì)對(duì)車(chē)輛的操縱穩(wěn)定性產(chǎn)生較大影響[1-14]。但是，這些研究只考慮風(fēng)載荷變化對(duì)汽車(chē)姿態(tài)的影響，實(shí)際上，汽車(chē)行駛姿態(tài)的變化也會(huì)反過(guò)來(lái)影響汽車(chē)側(cè)風(fēng)氣動(dòng)特性，且風(fēng)載荷距離車(chē)輛質(zhì)心的位置對(duì)車(chē)輛的操縱穩(wěn)定性有重要影響。本文著重于在空氣動(dòng)力學(xué)和車(chē)輛動(dòng)力學(xué)雙向耦合、相互影響下研究側(cè)風(fēng)對(duì)車(chē)輛穩(wěn)定性的影響。同時(shí)，通過(guò)提取風(fēng)壓中心線，分析風(fēng)壓中心線與車(chē)輛質(zhì)心的位置關(guān)系，改變氣動(dòng)力在Adams/Car模型中的加載位置，模擬風(fēng)壓中心線的移動(dòng)，對(duì)風(fēng)壓中心線與車(chē)輛質(zhì)心的距離進(jìn)行變量分析，為汽車(chē)造型開(kāi)發(fā)提供更加準(zhǔn)確的量化建議。

1雙向耦合聯(lián)合仿真模型首先，在CCM+軟件中獲得氣動(dòng)力和力矩與車(chē)速、流入角之間的關(guān)系曲線，再將這些曲線作為二維查表，導(dǎo)入MATLAB與Adams/Car聯(lián)合仿真模型中，通過(guò)Adamc/Car輸出車(chē)速和車(chē)輛在風(fēng)場(chǎng)中的姿態(tài)，計(jì)算出風(fēng)速與車(chē)速的合速度和流入角，再通過(guò)查表讀取氣動(dòng)力和力矩，并施加到Adams/Car模型中，形成閉環(huán)仿真。

1.1風(fēng)場(chǎng)的定義風(fēng)場(chǎng)主要模擬城市隔離帶中間縫隙帶來(lái)的側(cè)風(fēng)井噴效應(yīng)，其定義示意見(jiàn)圖1。車(chē)輛沿右側(cè)車(chē)道向前行駛，同時(shí)車(chē)身右側(cè)受側(cè)風(fēng)影響，具體風(fēng)速設(shè)置見(jiàn)表1[15]。

1.2氣動(dòng)力和氣動(dòng)力矩原始數(shù)據(jù)獲取通常，車(chē)輛在側(cè)風(fēng)中直線行駛時(shí)，流入角小于30°，當(dāng)風(fēng)速為10 m/s、車(chē)速為27.7 m/s（100 km/h）時(shí)，兩者的合速度小于145 km/h。在CCM+軟件中，流入角τ取0°～30°，車(chē)速取100～145 km/h進(jìn)行計(jì)算，獲得氣動(dòng)力和氣動(dòng)力矩與流入角以及合速度的原始關(guān)系曲線，見(jiàn)圖2～7。將該曲線作為MATLAB和Adams/Car聯(lián)合仿真模型中氣動(dòng)力和力矩二維查表的原始曲線。

1.3實(shí)時(shí)流入角的計(jì)算

1.3.1坐標(biāo)系定義取車(chē)輛局部坐標(biāo)系與Adams/Car坐標(biāo)系相同，x軸指向車(chē)尾，y軸指向車(chē)身右側(cè)，z軸垂直向上。風(fēng)場(chǎng)絕對(duì)坐標(biāo)系為：x軸由北指向南，y軸由西指向東，z軸垂直于地面向上。以上2個(gè)坐標(biāo)系的繞坐標(biāo)軸旋轉(zhuǎn)運(yùn)動(dòng)均符合右手定則。

1.3.2風(fēng)速和車(chē)速的夾角車(chē)輛在風(fēng)場(chǎng)中受到的風(fēng)速與車(chē)速的關(guān)系示意見(jiàn)圖8。根據(jù)車(chē)輛在風(fēng)場(chǎng)絕對(duì)坐標(biāo)系中繞z軸的旋轉(zhuǎn)角度θ和汽車(chē)的質(zhì)心側(cè)偏角β計(jì)算風(fēng)速與車(chē)速的夾角α，具體計(jì)算公式[15-16]見(jiàn)表2。α用于計(jì)算風(fēng)速和車(chē)速的合速度。

1.3.3風(fēng)速和車(chē)速合速度計(jì)算根據(jù)風(fēng)速和車(chē)速的夾角α以及車(chē)速和風(fēng)速的大小，通過(guò)余弦定理計(jì)算風(fēng)速和車(chē)速的合速度，作為聯(lián)合仿真模型中氣動(dòng)力和力矩二維查表的輸入?yún)?shù)之一。V2合=V2風(fēng)+V2車(chē)-2V風(fēng)V車(chē)cos（π-α）????? （1）式中：V車(chē)為車(chē)速;V風(fēng)為風(fēng)速;V合為風(fēng)速和車(chē)速的合速度;α為風(fēng)速與車(chē)速的夾角。

1.3.4實(shí)時(shí)流入角的計(jì)算當(dāng)車(chē)輛位于第一、四象限時(shí)，風(fēng)速與車(chē)速合速度的流入角τ等于合速度與車(chē)速的夾角減去車(chē)輛的質(zhì)心側(cè)偏角β;當(dāng)車(chē)輛位于第二、三象限時(shí)，τ等于合速度與車(chē)速的夾角加上質(zhì)心側(cè)偏角β。通過(guò)余弦定理計(jì)算合速度與車(chē)速的夾角δ，即V2風(fēng)= V2合+V2車(chē)-2V合V車(chē)cos δ???????? （2）將流入角τ作為聯(lián)合仿真模型中氣動(dòng)力和力矩二維查表的另外一個(gè)輸入?yún)?shù)，其計(jì)算公式見(jiàn)表3。

2單向模型說(shuō)明在分析側(cè)風(fēng)影響時(shí)，單向模型假設(shè)側(cè)風(fēng)始終垂直于汽車(chē)縱向，忽略行駛中汽車(chē)質(zhì)心側(cè)偏角的存在，且車(chē)速方向?yàn)檠仄?chē)縱向，同時(shí)忽略汽車(chē)受側(cè)風(fēng)影響后發(fā)生橫擺而造成流入角的變化。

3仿真工況的設(shè)計(jì)

3.1角階躍工況單向模型和雙向耦合模型的差異在角階躍工況下，分析單向模型和雙向耦合模型的差異，見(jiàn)圖9～11。隨著時(shí)間推移，單向模型和雙向耦合模型中的橫擺角、質(zhì)心側(cè)偏角和側(cè)向加速度結(jié)果差異越來(lái)越大，并且有隨時(shí)間增加而持續(xù)增大的趨勢(shì)。具體模型結(jié)果對(duì)比見(jiàn)表4。

3.2基于雙向耦合模型分析風(fēng)壓中心線和車(chē)輛質(zhì)心相對(duì)位置對(duì)操縱穩(wěn)定性的影響3.2.1開(kāi)環(huán)仿真時(shí)直線行駛當(dāng)行駛中的車(chē)輛受側(cè)風(fēng)影響時(shí)，評(píng)價(jià)無(wú)駕駛員干預(yù)情況下車(chē)輛的偏航程度，設(shè)定車(chē)速為100 km/h。為分析風(fēng)壓中心線位置對(duì)操縱穩(wěn)定性的影響，在CCM+軟件中提取風(fēng)壓中心線，發(fā)現(xiàn)車(chē)輛質(zhì)心位于風(fēng)壓線后方553 mm、上方95 mm處。想要得到具有準(zhǔn)確風(fēng)壓中心線位置的新造型非常困難，因此通過(guò)調(diào)整氣動(dòng)力在Adams/Car模型中的加載點(diǎn)模擬風(fēng)壓中線位置的變化?；陂_(kāi)環(huán)仿真和閉環(huán)仿真2個(gè)工況，評(píng)估風(fēng)壓中心線位置對(duì)操縱穩(wěn)定性的影響。在開(kāi)環(huán)仿真中，車(chē)輛行駛時(shí)風(fēng)壓中心線到車(chē)輛質(zhì)心距離對(duì)操縱穩(wěn)定性影響的分析結(jié)果見(jiàn)圖12。

隨著風(fēng)壓中心線后移，車(chē)輛的側(cè)向位移在一定范圍內(nèi)逐漸減小，但是當(dāng)風(fēng)壓中心線位于車(chē)輛質(zhì)心后方很遠(yuǎn)處時(shí)，車(chē)輛會(huì)在氣動(dòng)力矩的作用下向另一側(cè)偏移。

3.2.2閉環(huán)仿真時(shí)直線行駛當(dāng)行駛中的車(chē)輛受側(cè)風(fēng)影響時(shí)，評(píng)價(jià)駕駛員通過(guò)修正方向盤(pán)維持直線行駛的難易程度。閉環(huán)仿真中車(chē)輛行駛時(shí)風(fēng)壓中心線位置對(duì)操縱穩(wěn)定性的影響結(jié)果見(jiàn)圖13。當(dāng)風(fēng)壓中心線穿過(guò)車(chē)輛質(zhì)心，即氣動(dòng)力加載點(diǎn)后移553 mm、上移95 mm時(shí)，雖然氣動(dòng)力產(chǎn)生的橫擺力矩幾乎為0，但此時(shí)側(cè)向力仍然存在，故穩(wěn)態(tài)時(shí)方向盤(pán)仍需要較小的轉(zhuǎn)角以維持直線行駛，不能回正。隨著風(fēng)壓中心線繼續(xù)后移，當(dāng)氣動(dòng)力產(chǎn)生的橫擺力矩可以平衡側(cè)向力時(shí)，方向盤(pán)穩(wěn)態(tài)角度幾乎為0，即方向盤(pán)可以回正。此時(shí)，若風(fēng)壓中心線繼續(xù)后移，氣動(dòng)力產(chǎn)生的橫擺力矩會(huì)繼續(xù)增大，導(dǎo)致方向盤(pán)需向另一個(gè)方向給定一個(gè)角度，才能維持直線行駛，

4結(jié)論提出基于流入角實(shí)時(shí)變化的空氣動(dòng)力學(xué)和車(chē)輛動(dòng)力學(xué)雙向耦合分析方法，該方法可以更準(zhǔn)確地反應(yīng)氣動(dòng)力對(duì)操縱穩(wěn)定性的影響程度。以風(fēng)壓中心線與車(chē)輛質(zhì)心的相對(duì)位置為變量，分析開(kāi)環(huán)仿真和閉環(huán)仿真時(shí)氣動(dòng)力對(duì)操縱穩(wěn)定性的影響，得到以下結(jié)論。（1）隨著風(fēng)壓中心線由質(zhì)心前方移到質(zhì)心后方，開(kāi)環(huán)仿真時(shí)，車(chē)輛行駛的側(cè)向位移會(huì)減小，但是如果風(fēng)壓中心線位于車(chē)輛后方超過(guò)一定的距離，車(chē)輛會(huì)在氣動(dòng)力產(chǎn)生的橫擺力矩作用下向另外一側(cè)偏移。（2）閉環(huán)仿真時(shí)，若風(fēng)壓中心線穿過(guò)車(chē)輛質(zhì)心，由于側(cè)向力的存在，方向盤(pán)依然需要一個(gè)很小的轉(zhuǎn)角維持直線行駛，風(fēng)壓中心線需要繼續(xù)后移以保證方向盤(pán)穩(wěn)態(tài)時(shí)的回正，本文模型中風(fēng)壓中心線約位于質(zhì)心后方100 mm處可實(shí)現(xiàn)方向盤(pán)的穩(wěn)態(tài)回正。文中未考慮車(chē)輛發(fā)生小角度俯仰、側(cè)傾時(shí)氣動(dòng)力的變化，而事實(shí)上車(chē)輛的俯仰和側(cè)傾也會(huì)導(dǎo)致氣動(dòng)力和力矩變化，可以借用相似的手段進(jìn)行綜合分析。參考文獻(xiàn)：

[1]李堯. 汽車(chē)側(cè)風(fēng)氣動(dòng)特性及操穩(wěn)性影響的研究[D]. 成都： 西南交通大學(xué)， 2018.

[2]劉立寧， 石碩， 郗學(xué)峰， 等. 隔離帶對(duì)側(cè)風(fēng)條件下超車(chē)過(guò)程車(chē)輛氣動(dòng)力的影響[J]. 汽車(chē)工程， 2018， 40（4）： 437-442. DOI： 10.19562/j.chinasae.qcgc.2018.04.010.

[3]劉立寧. 側(cè)風(fēng)環(huán)境中不同運(yùn)動(dòng)狀態(tài)下廂式貨車(chē)瞬態(tài)氣動(dòng)特性研究[D]. 濟(jì)南： 山東大學(xué)， 2019， 57-108.

[4]趙叢琳， 馮悅， 鄭華. 軸荷對(duì)整車(chē)操穩(wěn)性能影響的仿真分析[J]. 計(jì)算機(jī)輔助工程， 2013， 22（6）： 7-11.

[5]余鳳梅， 朱暉， 楊志剛. 客車(chē)在不同側(cè)風(fēng)下的受力情況及安全性評(píng)價(jià)[J]. 計(jì)算機(jī)輔助工程， 2018， 27（1）： 7-14. DOI： 10.13340/j.cae.2018.01.002.

[6]呂佳琪. 側(cè)風(fēng)環(huán)境下的乘用車(chē)操縱穩(wěn)定性研究[D]. 長(zhǎng)春： 吉林大學(xué)， 2015.

[7]張喆. 側(cè)風(fēng)作用下重型載貨汽車(chē)氣動(dòng)特性與穩(wěn)定性的聯(lián)合研究[D]. 沈陽(yáng)： 吉林大學(xué)， 2012.

[8]吳萍. 側(cè)風(fēng)作用下轎車(chē)的氣動(dòng)特性的研究[D]. 廣州 ：華南理工大學(xué)， 2011.

[9]陳曉東. 大跨橋梁側(cè)風(fēng)行車(chē)安全分析[D]. 上海： 同濟(jì)大學(xué)， 2007.

[10]張輝. 基于主動(dòng)四輪轉(zhuǎn)向系統(tǒng)的高速汽車(chē)側(cè)風(fēng)穩(wěn)定性控制研究[D]. 杭州： 浙江大學(xué)， 2014.

[11]張琨. 側(cè)風(fēng)對(duì)公路行車(chē)安全影響的仿真分析與對(duì)策研究[D]. 西安： 長(zhǎng)安大學(xué)， 2015.

[12]陳寧. 側(cè)風(fēng)作用下橋上汽車(chē)行車(chē)安全性及防風(fēng)措施研究[D]. 成都： 西南交通大學(xué)， 2010.

[13]陳哲. 基于四輪轉(zhuǎn)向和差動(dòng)制動(dòng)聯(lián)合控制系統(tǒng)的高速汽車(chē)側(cè)風(fēng)穩(wěn)定性控制研究[D]. 杭州： 浙江大學(xué)， 2015.

[14]江路銘. 側(cè)風(fēng)下隊(duì)列行駛車(chē)輛的氣動(dòng)特性研究[D]. 長(zhǎng)春： 吉林大學(xué)， 2016.

[15]傅立敏. 汽車(chē)設(shè)計(jì)與空氣動(dòng)力學(xué)[M]. 北京： 機(jī)械工業(yè)出版社， 2010： 243-260.

[16]曼弗雷德， 亨寧. 汽車(chē)動(dòng)力學(xué)[M]. 5版. 陳蔭三， 余強(qiáng)， 譯. 北京： 清華大學(xué)出版社， 2019： 42-547.（編輯武曉英）