王霽飛，劉龍權(quán)

(上海交通大學 航空航天學院， 上海 200240)

0 引 言

飛行安全是飛機研究中最重視的部分，運輸類民航客機適航標準CCAR第25部803條款要求：對于客座量大于44座的飛行器，必須表明在應急情況下其最大乘坐量的乘員能在90 s內(nèi)從飛機撤離至地面[1]。因此，民航客機在通過適航認證并正式運營之前，必須表明其客艙安全性滿足適航規(guī)章的要求。然而，一次真機試驗需要涉及許多人員、設(shè)備和儀器，費用較高，且具有較高的危險性。以波音公司為例，一次客艙逃生試驗，總計耗費金額達到100萬美元以上[2]。1972-1980年間，F(xiàn)AA的7次全機逃生試驗共涉及人員2 571人，造成其中166人受傷。1981-1991年間，F(xiàn)AA全機逃生試驗共涉及3 761人，受傷人數(shù)達到212人[3]。1991年，在對MD-11寬體機進行應急撤離疏散實驗時，一位女性測試者因受傷而永久性的癱瘓[4]。2006年A380完成逃生測試時一人腳部骨折，多人受傷[5]。與之相比，數(shù)值仿真模擬既安全且可重復性高，又可以模擬各種極限情況，因此飛機客艙逃生數(shù)值模擬可作為試驗的有益補充。

目前國外已長期開展飛機客艙安全逃生相關(guān)仿真研究。據(jù)文獻調(diào)研，F(xiàn)AA早在1962年就已開展乘客應急撤離模型相關(guān)研究[6]；1970年，F(xiàn)AA開發(fā)了GPSS軟件，這是第一個能用來模擬飛機客艙應急撤離仿真的軟件；此后國外陸續(xù)開發(fā)GA、AIREVAC、MACEY、airEXODUS等多個用于民航客機90 s應急撤離的仿真算法和軟件[7]，其中英國格林威治大學開發(fā)出商用軟件airEXODUS能較好地模擬仿真各種應急撤離的場景，但關(guān)于飛機部分其軟件對外保密[8-10]。

相比之下，國內(nèi)對于飛機客艙逃生仿真的研究相對較少。李杰等[11]研究了飛機應急撤離的元胞自動機模型，并在DELMIA平臺實現(xiàn)建模和三維演示；徐進津[12]研究了飛機客艙疏散仿真的人為因素，總結(jié)了客艙和乘客的屬性以及一些人為因素數(shù)據(jù)。上述文獻中均采用元胞自動機模型模擬疏散人群的屬性，把疏散人群視為相同粒子，忽視了個體的差異化，與真實撤離實驗存在一定差別。陳晨、杜紅兵等[13-15]運用元胞自動機模型，并借鑒了智能體模型，開發(fā)仿真軟件CabinEvac，并研究乘客的應急心理狀態(tài)，但其對心理狀態(tài)的考慮并不全面，且未與相同狀態(tài)下的逃生試驗結(jié)果進行對比驗證。

智能體模型把關(guān)注點放在疏散體身上，認為模型對象能在一定的環(huán)境下，具有一定的行為規(guī)則和學習能力，通過感知外界因素調(diào)整自身的狀態(tài)，最終完成目標[16]。這種模型的優(yōu)點在于其能體現(xiàn)個體的差異性，使模擬仿真更接近于真實。

本文采用智能體模型，建立逃生智能體的行為準則，以最優(yōu)路徑算法A*算法為主要行為依據(jù)，根據(jù)擁擠情況選擇在A*算法中適合飛機的估值函數(shù)；結(jié)合個體與個體之間的協(xié)助與阻礙影響，個體與環(huán)境的影響，建立客艙應急撤離的仿真模型；將算法結(jié)果與前人實驗結(jié)果和權(quán)威軟件airEXODUS仿真結(jié)果進行對比。

1 疏散人群的智能體模型設(shè)計

1.1 智能體模型

智能體模型來源于人工智能，指的是模型對象在一定的環(huán)境下，具有一定的行為規(guī)則和學習能力，最終完成目標任務[16]，這是一種自下而上構(gòu)建的社會模型[17]。

基于智能體模型理論，結(jié)合運輸類飛機適航標準，本文中的智能體模型遵循的行為規(guī)則分為三大子規(guī)則，分別為路徑抉擇規(guī)則、移動規(guī)則和交互規(guī)則，如圖1所示。在路徑抉擇時，智能體通過A*算法求取各逃生出口對應的最優(yōu)路徑，然后比較不同出口的最優(yōu)路徑，最終做出選擇；在移動規(guī)則下，疏散體的移動速度受到初始移動速度與其他系數(shù)的影響；在交互規(guī)則下，擁擠、恐慌和互助分別對路徑抉擇中的A*算法估值函數(shù)和影響速度的系數(shù)產(chǎn)生影響。

圖1 智能體模型的移動規(guī)則

1.2 路徑抉擇規(guī)則

在路徑抉擇上，每個疏散智能體均選擇最優(yōu)路徑的方法前進，本文采用在最優(yōu)路徑研究領(lǐng)域被廣泛采用的A*算法。

1.2.1 A*算法原理

A*算法是一種人工智能中典型的路徑啟發(fā)式搜索優(yōu)化算法，通過選用合適的估值函數(shù)F(n)，朝著最可能的值進行求取，最終得到從起點S到終點E的最短且最合適的路線，其原理圖如圖2所示。

圖2 A*算法原理圖

A*算法的步驟：

(1) 設(shè)置開放(open)和關(guān)閉(close)兩個列表，起點定為S，終點定為E。

(2) 從S開始，把S放入open列表中。

(3) 判斷open列表是否為空，若為空,則宣告路徑尋找失敗，退出循環(huán)。

(4) open列表為非空集合，選取open列表中F(n)最小的節(jié)點T，將節(jié)點T從open列表中刪除，放入close列表。

(5) 檢查節(jié)點T是否為終點E，若是終點E，則找到路徑并退出循環(huán)。

(6) 若節(jié)點T不是終點E，則把節(jié)點T周圍的可擴展節(jié)點Ti(上下左右且不出界和非墻節(jié)點)做預處理，計算Ti節(jié)點的F(Ti)，并判斷Ti。若Ti不在open和close列表中，則加入open列表，并定義其父節(jié)點為T。若Ti已經(jīng)在open列表中，則比較當前F(Ti)與原先的F′(Ti)：F(Ti)較小，則說明當前路徑更優(yōu)，并更改父節(jié)點為當前T，用F(Ti)代替F′(Ti)；F′(Ti)較小則說明原路徑更優(yōu)，因此跳過。若Ti已經(jīng)在close列表中，則跳過，處理下一個節(jié)點。

(7) 返回步驟(3)，直到在步驟(3)或步驟(5)時跳出循環(huán)。

1.2.2 A*算法的估值函數(shù)

估值函數(shù)F(n)是A*算法中的一個核心，選擇合適的估值函數(shù)能正確地得到目標至終點的最優(yōu)路徑，A*算法選擇F(n)=G(n)+H(n)作為估值函數(shù)，其中G(n)為從起點S到達一個節(jié)點n所已經(jīng)付出的代價，H(n)為從節(jié)點n到終點E所估計需要付出的代價。

在路徑搜索中曼哈頓距離常被用于估值函數(shù)。曼哈頓距離指的是兩點之間的x軸距離和y軸距離之和[18]，即若存在兩點S(x1，y1)和E(x2，y2),則曼哈頓距離L為

L(S,E)=|x1-x2|+|y1-y2|

(1)

飛機的情況較為復雜，人與人的交互也對估值判斷有較大的影響，因此估值H(n)=L(S,E)+Hinter(n)，其中Hinter(n)為人群交互所產(chǎn)生的估值函數(shù)影響變化。

1.3 交互規(guī)則

人與人之間的交互主要對模型體有兩方面的影響，擁擠和恐慌會對路徑抉擇產(chǎn)生影響，同時恐慌與互相幫助也對人的移動速度產(chǎn)生了影響。

1.3.1 擁擠系數(shù)

擁擠對路徑抉擇具有較大的影響，擁擠的人數(shù)會被折算成對應的等價距離。

定義擁擠距離Hinter(S,E)=μ×people(S,E)，其中people(S,E)指的是從起點到終點兩個頂點所組成的矩形中所包含的人數(shù)，μ為擁擠系數(shù)。

因此，A*算法中估值函數(shù)H(n)=L(S,E)+Hinter(S,E)

1.3.2 恐慌系數(shù)

恐慌主要產(chǎn)生兩方面的影響，一方面適當?shù)目只艜谷司o張，促使疏散體移動速度變快[19]；另一方面過度的緊張會使得疏散人群變得不理智，產(chǎn)生移動停滯或反向移動等現(xiàn)象[20]。

本文采用SIS模型和情緒衰減管理模型[21]模擬疏散人群的恐慌現(xiàn)象。

恐慌情緒會對路徑?jīng)Q策產(chǎn)生影響，人群分為兩類人：易感者(Susceptibles)與感染者(Infectives)。易感者為情緒接收者，不向外傳遞情緒；感染者為情緒傳遞者，同時也能接收情緒。易感者與感染者可互相轉(zhuǎn)變，其轉(zhuǎn)變零界點為情緒閾值η。當情緒值達到ηc時，其行為呈現(xiàn)非理智狀態(tài)，將出現(xiàn)停滯或隨機移動的行為現(xiàn)象。

初始時刻的恐慌系數(shù)k通過正態(tài)分布確定[21]，如式(2)所示。

kt=0～N(μk,σk)

(2)

式中：k為恐慌系數(shù)；μk為正態(tài)分布的均值；σk為正態(tài)分布的標準差。

由于恐慌程度與年齡性別有關(guān)，建立不同年齡段和性別的初始恐慌因子，如表1所示。

表1 不同年齡、性別人群的初始恐慌因子

某一時刻某一人物恐慌系數(shù)k的計算公式，如式(3)所示。

(3)

一般情緒傳遞半徑為一格到兩格[21]，本文選用兩格作為情緒傳遞半徑。

(4)

疏散人群的恐慌會隨時間衰減，并且會離出口相接近時大幅度衰減。因此情緒衰減速率kpdecay的計算公式，如式(5)所示。

(5)

式中：α為疏散體在靠近出口處的情緒衰減速率，取值0.8；β為疏散體正常情況下的情緒衰減速率，取值0.1；tp為疏散體收到情緒傳遞的時間；d為疏散體離最近出口的曼哈頓距離；Rdoor為出口相近處的臨界曼哈頓距離。

當恐慌系數(shù)kp<ηc時，疏散體呈現(xiàn)較為理智狀態(tài)，此時kmove=λkp，其中λ為激勵調(diào)節(jié)因子；當恐慌系數(shù)kp≥ηc，疏散體呈現(xiàn)不理智狀態(tài)，呈現(xiàn)停滯。

1.3.3 乘客之間的互相協(xié)助因子

互相幫助可以使得疏散者的移動速度變快。根據(jù)AASK的飛機事故數(shù)據(jù)庫調(diào)查顯示，49.5%的乘客會在飛行過程中與同伴一起同行，其中65%為家庭關(guān)系，16%為朋友關(guān)系，8%為同事關(guān)系，11%為其他關(guān)系,不同的關(guān)系決定了幫助的程度[22-23]。在同伴關(guān)系中，9%的乘客在真實飛行過程中，提供了他人幫助，其中男性占據(jù)68%，女性占據(jù)36%。

因此協(xié)助系數(shù)s：

(6)

1.4 移動規(guī)則

1.4.1 乘客的速度

人員最大的移動速度取決于年齡、性別，是否攜帶嬰兒等其他因素影響，因此采用公式Vk=random(Vkmax,Vkmin)×λ，其中Vkmax為人員移動速度最大值，Vkmin為人員移動速度最小值，λ為速度影響系數(shù)。

人員最大移動速度采用GALEA整理的民機人員應急撤離速度分布統(tǒng)計數(shù)據(jù)，結(jié)合Fruin實驗統(tǒng)計[20]，數(shù)據(jù)是對于目前應急撤離實驗錄像影音中人群的速度統(tǒng)計得到，如表2所示。

表2 不同年齡性別人群應急撤離速度統(tǒng)計[24]

小于2歲的嬰兒被視為無行為能力者，且依附于其他人員且弱化該人員逃生速度，本文采用減少被依附人員30%的移動速度。

1.4.2 速度影響系數(shù)

速度影響系數(shù)λ主要受到恐慌系數(shù)k、互助系數(shù)s和障礙物的影響。

在客艙中障礙物分為兩種，可翻越式與不可翻越式。其中可翻越式主要為乘客的行李、座位處的把手，不可翻越式主要為座位、機上廚房與洗手間。根據(jù)《運輸類飛機適航標準》附錄J[1]，在適航逃生時，會有乘客人數(shù)一半的行李、毯子、枕頭等其他類似物品放在過道或出口造成輕微障礙。不同障礙物的速度影響系數(shù)如表3所示，得到速度影響系數(shù)λ=kmove×s×λelse。

表3 不同障礙物的速度影響系數(shù)

2 飛機客艙模型

2.1 機型選擇B737-700機型

本文選用B737-700120座級機型進行數(shù)值仿真。該機型客艙分布圖如圖3所示，左右兩側(cè)呈現(xiàn)對稱狀態(tài)，每側(cè)各有三個應急出口，共120個座位，為B737-700全普艙型機型。

圖3 B737-700客艙模型圖

2.2 應急設(shè)備使用情況

2.2.1 安全門使用情況

按照《運輸類飛機適航標準》25.807條[1]，許可的最大乘客數(shù)取決于機身上每一側(cè)的應急出口類型和數(shù)量，可容納最大乘客數(shù)規(guī)定如表4所示。

表4 應急出口類型與最大乘客數(shù)對應表

C型出口需在10 s內(nèi)完全展開，其他型出口需要在6 s內(nèi)完全展開[1]。結(jié)合文獻[25]可得C型、Ⅰ型與Ⅱ型門大致開門時間如表5所示。

表5 各類型門開門時間

在適航審定過程中，局方要求使用的安全門是開啟一半的門，根據(jù)中軸線對稱的每兩扇門使用其中一扇作為可用出口，用來測試能否通過90 s的應急撤離試驗。由于在仿真計算過程中，左門開啟與右門開啟影響甚微，因此全部使用一側(cè)作為開口。

2.2.2 逃生滑梯的使用情況

查閱文獻[26]，得到在正常離地高度3.5 m，滑梯展開與地面呈35°時，乘客與滑梯之間的摩擦系數(shù)為0.4，單人下滑時間為1.91 s。B737-700艙門離地高度約為2.7 m。因此t單人≈1.91/3.5×2.7，為1.47 s。按照CTSO-C69c標準，一段滑梯上應能同時承載3人[27]，因此當t=1.47/3=0.49 s時，下一位疏散體即可使用滑梯。

2.3 應急撤離時間

應急撤離時長由以下兩部分組成：

(1)T1——乘務員抵達相應安全門，若對應大門可用，則打開門，滑梯充氣的時間。

(2)T2——乘客開始準備撤離，在開門過程中排隊，開完門后依次乘坐應急逃生滑梯離開客艙的時間。

總撤離時長T=T1+T2。

3 模型有效性驗證

為了驗證數(shù)值仿真模擬的有效性，本文采用模擬B737-700機型，將本文數(shù)值模擬結(jié)果與airEXODUS軟件仿真結(jié)果[24]和前人實驗結(jié)果[28]進行對比。模擬客艙內(nèi)含有60名志愿者，包含38名青年男性和22名青年女性。客艙段選用了前十排座椅(總計60個座位)和兩個前排I型門。在真實實驗中一共重復四次，airEXODUS軟件仿真模擬12次，本文采用相同情況進行模擬20次結(jié)果，取最大與最小的模擬時間與實驗結(jié)果進行對比圖如圖4～圖5所示。

圖4 應急撤離數(shù)值模擬結(jié)果、前人實驗結(jié)果與airEXODUS軟件模擬結(jié)果對比圖

圖5 應急撤離數(shù)值模擬20次結(jié)果頻數(shù)分布圖

從圖4～圖5可以看出：本文的數(shù)值模擬結(jié)果與真實實驗結(jié)果和軟件airEXODUS仿真結(jié)果吻合(airEXOUDS仿真軟件為航空器撤離權(quán)威學者E.R.Galea帶領(lǐng)團隊開發(fā)[8-9，22-24])，在15～30 s，本文仿真結(jié)果較真實撤離有更寬的覆蓋范圍，這是因為真實實驗一共由四次組成，因此并沒有全面覆蓋應急撤離的真實時間范圍；同時，平均撤離時間較為接近，在25～40 s部分實線(本文模擬結(jié)果)包在長條點狀線(前人實驗結(jié)果)以外，較為合理。因此本文的仿真結(jié)果有效。

55名乘客疏散時間的數(shù)值模擬結(jié)果與前人實驗結(jié)果和airEXODUS仿真結(jié)果對比如表6所示,可以看出：本文算法模擬的乘客疏散時間與前人實驗結(jié)果和airEXODUS仿真結(jié)果接近。

表6 55名乘客疏散時間的數(shù)值模擬結(jié)果和

4 結(jié) 論

(1) 將適用于飛機疏散的A*算法仿真結(jié)果與前人實驗結(jié)果和國外權(quán)威軟件仿真結(jié)果進行對比，對比結(jié)果表明本文仿真算法是有效的。

(2) 本算法與傳統(tǒng)應急疏散撤離中使用的元胞自動機模型相比，考慮了疏散體之間的差異性，是一種從疏散體本身出發(fā)的算法。利用本算法，可以使仿真模擬更接近真實疏散情況。但本算法的求解效率還有待進一步研究與改善。