劉詩(shī)佳

（1 北京縱橫機(jī)電科技有限公司，北京 100094；2 中國(guó)鐵道科學(xué)研究院集團(tuán)有限公司 機(jī)車車輛研究所，北京 100081）

高速動(dòng)車組的設(shè)計(jì)和運(yùn)營(yíng)是非常復(fù)雜的系統(tǒng)工程。因此，借助電子計(jì)算機(jī)的高速和精準(zhǔn)運(yùn)算能力，提前分析并提高動(dòng)車組在軌道線路上的運(yùn)行效率，對(duì)提高運(yùn)營(yíng)效果尤為重要。

雖然國(guó)外先于我國(guó)開展高速動(dòng)車組仿真平臺(tái)的研究，例如北美的RAILSIM、歐洲的TRAINSTAR、日本的UTRAS［1］、俄羅斯的EMUTTCS［2］。與此同 時(shí)國(guó)內(nèi)不少機(jī)構(gòu)也展開了符合我國(guó)國(guó)情的高速動(dòng)車組仿真系統(tǒng)的研究開發(fā)。中國(guó)鐵道科學(xué)研究院集團(tuán)有限公司機(jī)車車輛研究所開發(fā)了高速動(dòng)車組牽引制動(dòng)仿真系統(tǒng)［3］，使用動(dòng)車組的動(dòng)力學(xué)模型推算動(dòng)車組的運(yùn)行狀態(tài)［4］，提供了方便直觀的動(dòng)車組控制系統(tǒng)，還繪制了軌道線路縱斷面的細(xì)節(jié)［5-6］。

以上仿真系統(tǒng)大多為離線仿真。為了將仿真系統(tǒng)用于指導(dǎo)司機(jī)實(shí)際操作，就需要對(duì)仿真系統(tǒng)進(jìn)行優(yōu)化，提高其運(yùn)行效率以便及時(shí)為司機(jī)控車給出正確精準(zhǔn)的指導(dǎo)。文中從仿真系統(tǒng)數(shù)據(jù)結(jié)構(gòu)角度出發(fā)，提出優(yōu)化方案，同時(shí)提出高速動(dòng)車組的能耗計(jì)算解決方案，為進(jìn)一步研究和對(duì)比節(jié)能控車策略，提供了幫助。

1 輔助駕駛系統(tǒng)需求

高速動(dòng)車組輔助駕駛系統(tǒng)的目的，是為了給司機(jī)提供控制高速動(dòng)車組行進(jìn)速度的操作指令序列，如在什么時(shí)間或位置，將控制手柄具體放在牽引、制動(dòng)、惰行或定速中的哪個(gè)檔位上。確保在這些操作建議都順利執(zhí)行的情況下，高速動(dòng)車組能按照時(shí)刻表準(zhǔn)時(shí)到達(dá)每個(gè)停車站，同時(shí)所消耗的電能盡量少，并提高旅客的舒適感。這些高速動(dòng)車組操作指令序列，將成為未來(lái)高速動(dòng)車組全自動(dòng)駕駛系統(tǒng)的重要基礎(chǔ)。

1.1 靈活的配置

文中使用的仿真系統(tǒng)在文獻(xiàn)［3］中有所描述，以其中的基本參數(shù)配置為基礎(chǔ)，對(duì)具體數(shù)據(jù)結(jié)構(gòu)、有關(guān)配置界面進(jìn)行了功能補(bǔ)充和性能優(yōu)化。

首先，對(duì)高速動(dòng)車組特性曲線的計(jì)算補(bǔ)充了散點(diǎn)描述方式，當(dāng)沒(méi)有合適的公式來(lái)描述對(duì)應(yīng)曲線時(shí)，我們用“線性平均值”算法基于速度來(lái)推算所需數(shù)值，如圖1 所示。

其次，補(bǔ)充和完善了以下配置及參數(shù)：

（1）時(shí)刻表及模擬乘客上下車后車重變化情況，讓仿真過(guò)程更加貼近實(shí)際運(yùn)營(yíng)情況，其配置如圖2 所示。

圖2 配置時(shí)刻表及模擬停車站上下車情況

（2）分相偏移、車站限速、軌道限速偏移等配置，讓仿真過(guò)程更加貼近司機(jī)在區(qū)間運(yùn)行控車時(shí)的實(shí)際操作，如圖3 所示。

圖3 配置相關(guān)的偏移參數(shù)

（3）模擬司機(jī)加速和減速的操作特性，設(shè)置界面如圖4 所示。通過(guò)在低速階段限制牽引電機(jī)發(fā)揮牽引力的比例，保證規(guī)劃和仿真時(shí)，最大程度上貼近司機(jī)操作高速動(dòng)車組進(jìn)出站時(shí)的實(shí)際情況。

圖4 高速動(dòng)車組加速和減速策略配置

1.2 有限的計(jì)算資源

作為車載系統(tǒng)的輔助模塊，其能夠使用的資源有限。讓輔助駕駛系統(tǒng)每次都能迅速地完成規(guī)劃，并給司機(jī)提供準(zhǔn)確的輔助駕駛信息，要對(duì)其數(shù)據(jù)進(jìn)行優(yōu)化，減少數(shù)據(jù)量。

為實(shí)現(xiàn)快速規(guī)劃，在現(xiàn)有仿真平臺(tái)的基礎(chǔ)上，對(duì)仿真系統(tǒng)的數(shù)據(jù)結(jié)構(gòu)進(jìn)行重點(diǎn)優(yōu)化設(shè)計(jì)［3］。

2 輔助駕駛系統(tǒng)優(yōu)化設(shè)計(jì)

首先，優(yōu)化軌道數(shù)據(jù)的記錄結(jié)構(gòu)，并對(duì)應(yīng)優(yōu)化仿真結(jié)果記錄算法，實(shí)現(xiàn)提高輔助駕駛指令序列規(guī)劃的及時(shí)性、降低數(shù)據(jù)量的目的。然后，設(shè)計(jì)了能耗拓?fù)鋱D功能，通過(guò)描述電能傳輸路徑及各電能傳輸節(jié)點(diǎn)的傳輸效率，計(jì)算該能耗拓?fù)浣Y(jié)構(gòu)所在車型運(yùn)行時(shí)所產(chǎn)生的能耗。

2.1 軌道數(shù)據(jù)結(jié)構(gòu)優(yōu)化

軌道在建設(shè)施工時(shí)通過(guò)“公里標(biāo)+本段軌道長(zhǎng)度”描述。但是，后續(xù)由于種種原因?qū)?shí)際軌道進(jìn)行維護(hù)更新之后，一些軌道分段長(zhǎng)度會(huì)發(fā)生變化，由此產(chǎn)生“長(zhǎng)短鏈”的情況。即有些軌道區(qū)段的長(zhǎng)度，不同于2 個(gè)相鄰公里標(biāo)的差值。例如，某條軌道線路上公里標(biāo)1.1～2.2 之間的軌道區(qū)段，長(zhǎng)度為800 m，而不是公里標(biāo)的差異1.1 km，這種情況就是“短鏈”，反之為“長(zhǎng)鏈”。

“長(zhǎng)短鏈”的存在導(dǎo)致高速動(dòng)車組運(yùn)行的實(shí)際距離不同于軌道上的公里標(biāo)，容易引入誤差，導(dǎo)致程序運(yùn)行混亂。為此提出以下優(yōu)化手段：

首先，用“絕對(duì)位置”替換“公里標(biāo)+本段軌道長(zhǎng)度”記錄軌道數(shù)據(jù)，并定位高速動(dòng)車組在軌道上的位置。絕對(duì)位置，以軌道線路的起點(diǎn)為位置0，軌道上每個(gè)點(diǎn)的位置，以動(dòng)車組從位置0 行駛到該位置時(shí)，實(shí)際行駛的距離標(biāo)注。這個(gè)數(shù)值向軌道終點(diǎn)方向遞增，與軌道是上行還是下行無(wú)關(guān)。

通過(guò)轉(zhuǎn)換，軌道線路上的每個(gè)點(diǎn)都有了唯一的“絕對(duì)位置”，且跟原本的“公里標(biāo)+本段軌道長(zhǎng)度”相互轉(zhuǎn)換、一一對(duì)應(yīng)。“絕對(duì)位置”遞增或遞減的距離，對(duì)應(yīng)高速動(dòng)車組前進(jìn)或后退時(shí)移動(dòng)的實(shí)際距離，這樣可以極大地簡(jiǎn)化動(dòng)車組所在位置的相關(guān)運(yùn)算。

其次，軌道的數(shù)據(jù)由坡度、曲線、橋隧、分相位置4 個(gè)列表組成。坡度記錄軌道每個(gè)區(qū)段的起點(diǎn)公里標(biāo)、長(zhǎng)度、坡度值，以及停車車站的站中點(diǎn)；曲線記錄軌道所有曲線的開始位置、曲線半徑、曲線長(zhǎng)度等；橋隧記錄軌道線路經(jīng)過(guò)的所有橋梁、隧道的公里標(biāo)和長(zhǎng)度；分相列表記錄軌道分相開始和結(jié)束的位置，在分相路段內(nèi)，高速動(dòng)車組失去電力供應(yīng)，被迫惰行。

將以上4 個(gè)列表的每個(gè)分段開始和結(jié)束位置，都轉(zhuǎn)化為“絕對(duì)位置”記錄后，把這些分段的點(diǎn)合并為一個(gè)“軌道描述列表”。這個(gè)列表的每一條都記錄了一段軌道的參數(shù)，包括坡度、車站、曲線、橋隧、分相等參數(shù)，確保在同一條記錄描述的軌道長(zhǎng)度內(nèi)，這些參數(shù)不會(huì)變化。

仿真時(shí)，只要高速動(dòng)車組在“軌道描述列表”同一段軌道內(nèi)運(yùn)行，就不必重復(fù)查詢軌道參數(shù)，且來(lái)自軌道的阻力為恒定值。如果高速動(dòng)車組為恒速運(yùn)行，則其風(fēng)阻也是恒定的，那么對(duì)應(yīng)的手柄控制就是固定的。在該段軌道內(nèi)運(yùn)行狀態(tài)的描述就可以極大地簡(jiǎn)化。

2.2 仿真結(jié)果優(yōu)化

文獻(xiàn)［3］中描述的散點(diǎn)方式記錄高速動(dòng)車組的運(yùn)行曲線。每個(gè)數(shù)據(jù)點(diǎn)記錄了高速動(dòng)車組在軌道上的位置、車速、時(shí)間、手柄位置、加減速相關(guān)的力學(xué)參數(shù)及輪周功率等數(shù)據(jù)。

但是，文獻(xiàn)［3］中使用固定的時(shí)間間隔（如0.1 s）生成記錄，導(dǎo)致2 個(gè)問(wèn)題：其一，當(dāng)2 個(gè)數(shù)據(jù)點(diǎn)之間存在變坡點(diǎn)的時(shí)候，坡度變化引入誤差，且該誤差在不斷地積累；其二，數(shù)據(jù)點(diǎn)過(guò)多，即便高速動(dòng)車組以350 km/h的速度運(yùn)行，每公里也需要至少102個(gè)節(jié)點(diǎn)，以武漢到廣州約1 100 多公里的線路為例，整條線路仿真數(shù)據(jù)點(diǎn)超過(guò)11.3 萬(wàn)個(gè)，如此多的數(shù)據(jù)點(diǎn)顯然會(huì)減慢有關(guān)數(shù)據(jù)的處理速度。

文中提出如下步驟對(duì)數(shù)據(jù)進(jìn)行優(yōu)化。

首先，將仿真結(jié)果中的節(jié)點(diǎn)位置，與“軌道描述列表”中的節(jié)點(diǎn)對(duì)應(yīng)。規(guī)則如下：

（1）如果動(dòng)車組在一個(gè)路段內(nèi)恒速運(yùn)行，則只需一個(gè)對(duì)應(yīng)該路段起點(diǎn)位置的數(shù)據(jù)點(diǎn)記錄有關(guān)狀態(tài)，因?yàn)樵谕粋€(gè)路段內(nèi)高速動(dòng)車組的控制參數(shù)和運(yùn)行狀態(tài)一樣。

（2）如果動(dòng)車組在一個(gè)路段內(nèi)無(wú)法保持恒速運(yùn)行，則2 個(gè)數(shù)據(jù)點(diǎn)之間的時(shí)間差以固定時(shí)間間隔記錄，如0.1 s。但是，如果下一個(gè)數(shù)據(jù)點(diǎn)的位置離開了當(dāng)前路段，則重新定位下一個(gè)數(shù)據(jù)點(diǎn)，將其定位到下個(gè)路段開始的位置。

其次，對(duì)于軌道線路中會(huì)重復(fù)使用的曲線序列，提前進(jìn)行計(jì)算的時(shí)候也遵循以上規(guī)則。這些曲線序列包括：停車車站的進(jìn)站停車曲線、出站發(fā)車曲線、分相惰行區(qū)域的強(qiáng)制惰行曲線、可能超速路段的提前減速曲線。由于加速策略和減速策略在配置中指定，不論后續(xù)的曲線優(yōu)化執(zhí)行多少次，以上這些曲線在最貼近限速駕駛的情況下，都是固定的。仿真系統(tǒng)提前計(jì)算一次，后續(xù)遇到的時(shí)候直接復(fù)制這一段曲線數(shù)據(jù)，即可極大地減少實(shí)時(shí)規(guī)劃時(shí)間。

通過(guò)以上兩步優(yōu)化，可以大量削減描述仿真結(jié)果的“位置—速度曲線”所需的數(shù)據(jù)點(diǎn)總數(shù)，提高計(jì)算效率。

2.3 靈活的能耗計(jì)算

輔助駕駛的重要功能之一是節(jié)能，這需要能耗計(jì)算為基礎(chǔ)。但不同高速動(dòng)車組內(nèi)部對(duì)于電力的使用方式不盡相同，如果針對(duì)每個(gè)車型單獨(dú)編程計(jì)算其能耗，會(huì)極大地增加編程工作量，且靈活性不足。

優(yōu)化設(shè)計(jì)就是增加電能傳輸拓?fù)鋱D功能。通過(guò)描述高速動(dòng)車組內(nèi)部電能傳輸?shù)耐負(fù)浣Y(jié)構(gòu)，靈活準(zhǔn)確地計(jì)算其能耗，如圖5、圖6 所示。

圖5 電能拓?fù)鋱D方案1

圖6 電能拓?fù)鋱D方案2

圖5 中描述了一種能耗拓?fù)浣Y(jié)構(gòu)。每個(gè)方塊為一個(gè)能耗節(jié)點(diǎn)，“節(jié)點(diǎn)類型+編號(hào)”在方塊左上角粗體字標(biāo)注，節(jié)點(diǎn)可設(shè)置不同名稱，效率類型方塊可設(shè)置不同的電能使用效率參數(shù)，如“效率1”（圖中紅色文字標(biāo)注“B，動(dòng)態(tài)效率”）方塊的效率為根據(jù)速度而變化的（其具體變化規(guī)律見紅色文字“D，速度，效率”標(biāo)注處），而“效率2”（圖中紅色文字標(biāo)注“C，固定效率”）方塊的效率為不隨速度變化的固定值?！澳芎慕M件”為直接耗電的模塊，如冷卻風(fēng)機(jī)（圖中紅色文字標(biāo)注“A，直接耗電”）、空調(diào)等；“輪周功率”為高速動(dòng)車組在不同速度下的牽引系統(tǒng)傳遞到驅(qū)動(dòng)輪上發(fā)揮出的功率。整個(gè)能耗拓?fù)浣Y(jié)構(gòu)，用箭頭連接各個(gè)方塊來(lái)描述：電能從供電網(wǎng)，經(jīng)過(guò)各種能量轉(zhuǎn)換節(jié)點(diǎn)，傳送到各個(gè)耗電方塊的過(guò)程。

在圖6中，描述了另一個(gè)型號(hào)的高速動(dòng)車組內(nèi)部能耗拓?fù)鋱D，與圖5 中的差異體現(xiàn)在“能耗組件1”的電力來(lái)源不同，以及“效率1”的效率從隨速度變化改為固定值。

能耗計(jì)算的過(guò)程，就是通過(guò)這個(gè)拓?fù)浣Y(jié)構(gòu)，從“輪周功率”和各“能耗組件”，沿著能耗拓?fù)涞穆肪€倒推，計(jì)算出不同速度下，高速動(dòng)車組從供電網(wǎng)獲得的總功率，再按照時(shí)間積分得到總能耗。計(jì)算過(guò)程中，正數(shù)表示耗電量，負(fù)數(shù)表示再生制動(dòng)的發(fā)電量，可以分別計(jì)算出來(lái)。

高速動(dòng)車組仿真運(yùn)行的每個(gè)數(shù)據(jù)點(diǎn)，都有速度、手柄、牽引力等參數(shù)。根據(jù)速度和牽引力，可以算出每個(gè)數(shù)據(jù)點(diǎn)處高速動(dòng)車組車輪上的功率，即“輪周功率”。輪周功率乘以時(shí)間得到車輪能耗，通過(guò)對(duì)能耗拓?fù)鋱D的逆向推導(dǎo)，可以得到高速動(dòng)車組從供電網(wǎng)上獲取的電能。再對(duì)運(yùn)行曲線上所有節(jié)點(diǎn)的能耗進(jìn)行積分，就可以累計(jì)出高速動(dòng)車組運(yùn)行過(guò)程中消耗的電能，以及再生制動(dòng)生成的電能。

能耗計(jì)算時(shí)，如果不考慮能耗拓?fù)浣Y(jié)構(gòu)，則計(jì)算出的能耗會(huì)遠(yuǎn)小于實(shí)際能耗，因?yàn)槊總€(gè)電能的轉(zhuǎn)換節(jié)點(diǎn)都會(huì)產(chǎn)生額外能耗，而且這個(gè)額外能耗會(huì)根據(jù)速度、運(yùn)行距離的不同而產(chǎn)生顯著變化。

沒(méi)有以上能耗拓?fù)鋱D功能的仿真軟件，為了準(zhǔn)確計(jì)算不同車型的能耗，需要根據(jù)每個(gè)車型的能耗拓?fù)鋱D進(jìn)行單獨(dú)編程，這會(huì)大大增加軟件的規(guī)模和出錯(cuò)的概率。

3 仿真結(jié)果

以某線路及某高速動(dòng)車組車型為例，按照時(shí)刻表的出站和到站時(shí)間為標(biāo)準(zhǔn)，進(jìn)行仿真運(yùn)行。將對(duì)比數(shù)據(jù)點(diǎn)個(gè)數(shù)、精準(zhǔn)停車、能耗計(jì)算這幾個(gè)方面。

高速動(dòng)車組運(yùn)行曲線記錄的數(shù)據(jù)點(diǎn)個(gè)數(shù)優(yōu)化前后對(duì)比，如圖7、圖8 所示。在這2 個(gè)圖中，左側(cè)為速度軸，下部為軌道線路參數(shù)，包括坡度、曲線，上部為圖注，中部為限速、分相、車站、位置—速度曲線。其中，“位置—速度曲線”以散點(diǎn)表示。圖中的路段為某站到公里標(biāo)1 273.6（這也是一個(gè)變坡點(diǎn)）這段區(qū)間內(nèi)的“位置—速度曲線”。

圖7 優(yōu)化前“位置—速度曲線”上的數(shù)據(jù)點(diǎn)分布

圖8 優(yōu)化后“位置—速度曲線”上的數(shù)據(jù)點(diǎn)分布

優(yōu)化前的“位置—速度曲線”數(shù)據(jù)點(diǎn)體現(xiàn)在圖7中，優(yōu)化后的“位置—速度曲線”數(shù)據(jù)點(diǎn)體現(xiàn)在圖8 中。可見，在某站到公里標(biāo)1 273.6處，優(yōu)化前的數(shù)據(jù)點(diǎn)上百個(gè)，而優(yōu)化后的數(shù)據(jù)點(diǎn)只有1 個(gè)。此外，在某站左側(cè)的數(shù)據(jù)點(diǎn)對(duì)比中，優(yōu)化后的數(shù)據(jù)點(diǎn)比優(yōu)化前同樣大幅度減少。

精準(zhǔn)停車方面，優(yōu)化前后的“位置—速度曲線”對(duì)比如圖9、圖10 所示，其中用顏色標(biāo)注了高速動(dòng)車組的運(yùn)行狀態(tài)，粉色為制動(dòng)，藍(lán)色為牽引，綠色為惰行。

圖9 優(yōu)化前“位置—速度曲線”的停車誤差

圖10 優(yōu)化后“位置—速度曲線”的停車誤差

優(yōu)化前的停車位置如圖9 所示，停車存在大約300 m的誤差。優(yōu)化后的停車位置如圖10 所示，停車誤差幾乎不可見。

能耗計(jì)算時(shí)，使用圖5 和圖6的拓?fù)浣Y(jié)構(gòu)會(huì)產(chǎn)生不同的總能耗。使用圖10 中的“位置—速度曲線”數(shù)據(jù)，按圖5的電能拓?fù)浞桨赣?jì)算的總能耗為17 244.86 kW?h，而按圖6的電能拓?fù)浞桨赣?jì)算的總能耗則為20 522.68 kW?h。

4 結(jié)論

通過(guò)以上仿真對(duì)比分析，可以看到優(yōu)化后的高速動(dòng)車組輔助駕駛系統(tǒng)實(shí)現(xiàn)了既定目標(biāo)。

通過(guò)軌道數(shù)據(jù)記錄的優(yōu)化設(shè)計(jì)，并對(duì)應(yīng)調(diào)整算法，確保仿真結(jié)果數(shù)據(jù)點(diǎn)大幅度減少；由于“位置—速度曲線”的數(shù)據(jù)點(diǎn)能對(duì)應(yīng)上“軌道描述列表”中的數(shù)據(jù)點(diǎn)，停車位置更加精準(zhǔn)。

通過(guò)能耗拓?fù)鋱D功能，不必根據(jù)每個(gè)車型重新編碼，而只需根據(jù)車型構(gòu)建能耗拓?fù)鋱D，從而極大地減少了編碼工作量。這樣，不但可以靈活準(zhǔn)確地計(jì)算不同控車策略產(chǎn)生的能耗，還能方便地對(duì)比不同車型產(chǎn)生的能耗差異。

實(shí)現(xiàn)了這些優(yōu)化之后，后續(xù)可以進(jìn)行基于能耗優(yōu)化的控車策略研究。