唐 婕

(西南交通大學(xué)信息科學(xué)與技術(shù)學(xué)院，611756，成都∥助教)

有軌電車在線路建設(shè)之初，需要通過仿真有軌電車運行情況來計算其正線通過能力、折返能力及出入庫能力。我國現(xiàn)代有軌電車主要采用司機目視瞭望和司機手動操作的方式控制電車運行。為有效模擬現(xiàn)代有軌電車的實際運行過程，本文基于Wiedemann模型，根據(jù)有軌電車駕駛員的生理反應(yīng)和心理反應(yīng)分析其駕駛行為，以此來仿真有軌電車的運行過程。

1 有軌電車運行仿真模型的類型

社會車輛運行模型主要有跟馳模型和換道模型。有軌電車具有軌道交通的特性，無法在運行過程中換道超車，而其追蹤模式與社會車輛相似，因此有軌電車運行模型只考慮其跟馳特性。

1.1 社會車輛跟馳模型

社會車輛跟馳模型主要用于模擬分析駕駛員在受外界因素影響時控制車輛運行的微觀行駛數(shù)據(jù)。其主要分為：刺激-反應(yīng)類模型、安全距離模型、人工智能模型和生理-心理類模型。

1.2 有軌電車運行仿真模型的選用

有軌電車在實際運營過程中出現(xiàn)跟馳的情況較少，大多為追蹤前方信號機或站臺等固定停車點的情況。因此，有軌電車不適用研究擁堵情況的刺激-反應(yīng)類模型和安全距離類模型。

生理-心理類模型更貼近實際的駕駛員控車過程，也更能描述日常生活中所見的司機駕駛行為。故本文在典型生理-心理類模型——Wiedemann模型的基礎(chǔ)上研究并建立有軌電車運行仿真模型。

2 有軌電車運行仿真模型建立

2.1 Wiedemann模型

Wiedemann模型根據(jù)前后跟馳車輛的相對速度和相對距離的不同，將車輛運行狀態(tài)劃分為緊急制動區(qū)、跟馳行駛區(qū)、車輛接近區(qū)及自由行駛區(qū)[4]。

2.2 有軌電車運行仿真模型建立

當(dāng)駕駛員在控制有軌電車追蹤前行有軌電車時，有軌電車運行仿真模型按Wiedemann模型也分為4個駕駛員感知區(qū)域：緊急制動區(qū)、跟馳行駛區(qū)、車輛接近區(qū)和自由行駛區(qū)，如圖1所示。圖1中，AX為水平直線，表示列車在靜止時所期望的停車間隔距離；ABX曲線為列車在低速運行時的最小期望跟車運行距離曲線；SDX曲線為最大跟車運行距離曲線；SDV曲線為駕駛員意識到正在接近一輛低速運行列車的臨界曲線；CLDV曲線則表示當(dāng)跟車運行距離較小，駕駛員意識到車間距在減少時的臨界曲線，默認CLDV與SDV曲線重疊；OPDV曲線表示駕駛員注意到前后車之間的距離隨時間而增加的臨界曲線[5]。

圖1 有軌電車運行狀態(tài)感知區(qū)域劃分

當(dāng)有軌電車的運行前方為路口信號機、道岔信號機和車站等固定限制點時，有軌電車駕駛員等效于追蹤一列速度為0且長度為0的列車，此時有軌電車駕駛員模型主要劃分緊急制動區(qū)、列車接近區(qū)及自由行駛區(qū)，如圖2所示。此時SDX、ABX及AX曲線重合。

圖2 追蹤前方為紅燈或車站時模型感知域劃分

當(dāng)有軌電車處于不同的行駛感知區(qū)域時，其運行狀態(tài)不同。由于有軌電車牽引制動特性與社會車輛有所區(qū)別，因此在不同感知區(qū)域中，需要結(jié)合牽引制動特性、駕駛員反應(yīng)、冒進情況及物理動力學(xué)特性等為不同區(qū)域的子計算模型進行建模。

2.2.1 自由行駛區(qū)

當(dāng)有軌電車處于在自由行駛狀態(tài)區(qū)時，司機感知目視前方?jīng)]有限制點或離限制點距離較遠，可根據(jù)區(qū)間限速、列車的牽引-制動特性和舒適度等約束條件來控制電車運行。自由行駛區(qū)子模型計算步驟如下：

步驟一：在線路限速條件下，根據(jù)司機冒險程度，司機控制有軌電車追蹤期望最大運行速度。不同類別司機的冒險程度不同。根據(jù)線路限速和當(dāng)前速度來計算列車最大運行加速度，以便列車運行至司機期望的最大速度，有軌電車在自由行駛區(qū)的最大加速度為：

aMax_t=(αvPSR-v)

(1)

式中：

aMax_t——列車運行時間為t后的列車允許最大加速度；

α——司機冒險系數(shù)，是指司機能追蹤線路限速的最大比值，即司機期望的最大速度與線路限速的比值(α取值越大表示司機越冒進)；

vPSR——列車所在區(qū)段的限制速度；

v——當(dāng)前列車速度。

步驟二：根據(jù)下一區(qū)段的限速條件和當(dāng)前行駛速度計算提前減速區(qū)。當(dāng)下一區(qū)段的線路限速小于當(dāng)前區(qū)段的限速,且下一區(qū)段限速的起始位置與列車當(dāng)前位置的距離差小于提前減速距離時，司機有必要提前減速，以避免因減速不及時導(dǎo)致的有軌電車超速情況。故需計算提前減速后的最大加速度。將式(1)中的vPSR替換成下一區(qū)段的列車限速vnext_PSR，則提前減速距離S為：

(2)

式中：

adefault——默認減速度。

步驟三：考慮有軌電車運行的牽引特性和制動特性。該特性限制有軌電車在v下能達到的最大加速度或最小減速度，列車運行過程中最大加速度應(yīng)控制在列車牽引-制動特性曲線范圍內(nèi)。

步驟四：考慮舒適度約束條件，要求電車運行平穩(wěn)。司機控制電車達到上述步驟所得的最大加速度時應(yīng)保障乘客的舒適度[6]。

at=a+kt

(3)

式中：

t——列車運行時間常數(shù)，表示列車運行的固定時間；

at——本次列車運行t后列車的加速度；

a——本次列車當(dāng)前的運行加速度；

k——沖擊率，是指加速度相對于時間的變化率；k值越小，旅客的舒適度越高；當(dāng)k小于一定值的時候，加速度變化對人體舒適度的影響很小，可以認為k為一個恒定值。

2.2.2 列車接近區(qū)

當(dāng)有軌電車進入接近行駛狀態(tài)時，司機感知接近前方限制點，于是司機期望將車速從高于前車車速降為低于前車車速。參考Wiedemann模型，列車接近區(qū)子模型為：

(4)

式中：

SABX——列車在低速運行時的最小期望跟車運行距離；

Δv——當(dāng)前列車與前車的速度差；

Δx——當(dāng)前列車與前車的距離；

L——前車的車長；

af——前車的運行加速度。

當(dāng)前方追蹤目標為信號機或者車站時,駕駛員采用的加速度為：

(5)

式中：

Δx1——后車車頭到限制點的距離。

當(dāng)有軌電車進入跟馳行駛狀態(tài)時，司機反應(yīng)到當(dāng)前位置與前車距離較近，期望與前車在安全間隔距離下保持相同的行駛狀態(tài)。以此為目標，則有：

Δx=Ssafe+SABX

(6)

(7)

v0=v+a(t1+t2)+aLt3

(8)

式中：

Ssafe——當(dāng)前列車與前車的安全間隔距離；

vf——前車的速度；

v0——司機反應(yīng)后的列車初速度；

a——本次列車當(dāng)前的的運行加速度；

at——本次列車運行t后的加速度；

2.1 預(yù)防試驗結(jié)果 由表1可知，在預(yù)防試驗中，Ⅰ～Ⅵ組發(fā)病雞的癥狀基本一致,為典型的IBD病變,血清檢查結(jié)果呈陽性；對照組雛雞發(fā)病率明顯高于試驗組；Ⅵ組因未服用中草藥制劑，發(fā)病率高且病情最嚴重,試驗期內(nèi)7只雛雞發(fā)病，發(fā)病率高達17.5%。Ⅳ組雛雞發(fā)病癥狀最輕且發(fā)病率最低，試驗期內(nèi)僅有1只雛雞發(fā)病，發(fā)病率為2.5%，將制得的復(fù)方中草藥制劑按4.0 mL/只服用，預(yù)防效果最好，表明中草藥制劑對雛雞IBD有明顯的預(yù)防效果。

t1——司機反應(yīng)時間；

t2——司機操作時間；

aL——本次列車當(dāng)前的線路阻力加速度；

t3——施加牽引制動所需時間[7]。

將式(7)和(8)代入式(6)，得到：

(9)

2.2.4 緊急制動區(qū)

當(dāng)有軌電車進入緊急制動狀態(tài)時，司機感知本車與前方限制點非常接近，需要控制車輛安全停車以防止冒進信號或發(fā)生碰撞。

司機控制電車追蹤前車時，其追蹤目標為距離前車尾部一定安全距離的停車點，有：

Δx=Ssafe

(10)

(11)

將式(11)和式(8)代入式(10)，得到：

(12)

當(dāng)追蹤前方限制點為信號機或者車站等固定干擾點時，司機的反應(yīng)時間和操作時間很小，可忽略不計，則：

(13)

3 模型參數(shù)校正及仿真結(jié)果驗證

對有軌電車駕駛員模型進行仿真程序設(shè)計，得到有軌電車運行仿真系統(tǒng)結(jié)構(gòu)如圖3所示。

圖3 有軌電車運行仿真系統(tǒng)結(jié)構(gòu)

3.1 模型參數(shù)校正

Wiedemann模型原本是描述社會交通車輛運行過程的模型。為使模型更貼合有軌電車的運行場景，需要對Wiedemann模型進行參數(shù)校正。

模型參數(shù)可利用有軌電車在行駛過程中受牽引制動能力、舒適度和嚴格的速度限制等限制條件的約束來進行校正。約束條件可表示為：

(14)

式中：

vn——列車在運行了n個t時的速度；

vmaxn——列車啊運行了n個t時的最大允許運行速度；

an——列車在運行了n個t時的加速度；

an-1——列車運行了(n-1)個t時的加速度；

kn——列車在運行了n個t時的沖擊率；

abraking_n——有軌電車的允許制動加速度；

amotoring_n——有軌電車的允許牽引加速度。

由于自由行駛區(qū)子模型是根據(jù)約束條件建立的，所以有軌電車需要校正的感知區(qū)域為車輛接近區(qū)和緊急制動區(qū)，即與ABX和AX曲線相鄰的區(qū)域，故需要校正的參數(shù)主要為與SAX(列車在靜止時所期望的停車間隔距離)相關(guān)的加法因子和乘法因子。校正過程如下：

1) 獲取多組深圳龍華新區(qū)現(xiàn)代有軌電車線(以下簡為“龍華線”)大布頭站到河?xùn)|站區(qū)間的運行數(shù)據(jù)，并對每組數(shù)據(jù)的運行時間進行正態(tài)分布擬合。擬合得出，實際運行時間的期望值為65.468 8 s，置信度為99.9%的置信區(qū)間為[62.460 5,68.477]。

2) 將與SAX相關(guān)的加法因子和乘法因子數(shù)值代入至仿真程序生成運行結(jié)果。檢測該結(jié)果是否滿足式(14)的約束條件，且其區(qū)間運行時間是否在實際運行時間的置信區(qū)間范圍內(nèi)，同時對比仿真與實際的運行曲線是否相符。若滿足以上條件則表明模型參數(shù)校正完成，否則繼續(xù)代入其他數(shù)值進行仿真檢測。通過反復(fù)代入檢測得出，當(dāng)與SAX相關(guān)的加法因子為3.6vn+10，乘法因子為1.25時，模型滿足約束條件，且當(dāng)α=0.8時駕駛員在不同反應(yīng)時間下生成的仿真運行時間基本落在的置信區(qū)間范圍內(nèi)。具體仿真結(jié)果如表2所示。

表1 龍華線大布頭站—河?xùn)|站α=0.8時仿真結(jié)果

仿真列車運行曲線與現(xiàn)場運行曲線對比如圖4所示。

由圖4可見，仿真運行時間符合實際情況，仿真所得的站間運行總時長為65.627 s。在站間7 482.32 m處，有軌電車運行速度為14.944 3 m/s，行駛狀態(tài)由自由行駛區(qū)進入接近行駛區(qū)；在7 655.39 m處電車速度為2.042 07 m/s，電車行駛狀態(tài)進入緊急制動區(qū)，并準備到站停車。

圖4 龍華線大布頭站—河?xùn)|站仿真與現(xiàn)場數(shù)據(jù)對比

3.2 仿真結(jié)果驗證

驗證線路選定龍華線河西站至觀城站線路。利用校正后的模型仿真計算河西站至觀城站的有軌電車速度-距離圖，仿真結(jié)果與實際運行結(jié)果對比如圖5所示。

由圖5可見，仿真結(jié)果基本符合現(xiàn)場實際運行情況，證明本模型經(jīng)校正后可用于其他線路的有軌電車運行仿真。

圖5 龍華線河西站—觀域站仿真與現(xiàn)場數(shù)據(jù)對比

4 結(jié)語

基于Wiedemann模型，應(yīng)用龍華線數(shù)據(jù)進行仿真。仿真結(jié)果表明：本文所提出的有軌電車運行仿真模型能有效地模擬有軌電車的運行過程。在有軌電車項目投標階段，使用針對有軌電車司機駕駛行為和有軌電車牽引制動特性的模型能為分析有軌電車的通過能力及評價有軌電車信號系統(tǒng)運營能力提供依據(jù)。