不同荷載模式下礦用救生艙受力性能的數(shù)值模擬＊

李冀龍，唐亞男，劉軒銘

（1.哈爾濱工業(yè)大學(xué)結(jié)構(gòu)工程災(zāi)變與控制教育部重點(diǎn)實(shí)驗(yàn)室，黑龍江哈爾濱150090；2.哈爾濱工業(yè)大學(xué)土木工程學(xué)院，黑龍江哈爾濱150090；3.Department of Civil and Environmental Engineering，University of California，Berkeley，CA，the US，94710）

不同荷載模式下礦用救生艙受力性能的數(shù)值模擬＊

李冀龍1，2，3，唐亞男2，劉軒銘2

（1.哈爾濱工業(yè)大學(xué)結(jié)構(gòu)工程災(zāi)變與控制教育部重點(diǎn)實(shí)驗(yàn)室，黑龍江哈爾濱150090；2.哈爾濱工業(yè)大學(xué)土木工程學(xué)院，黑龍江哈爾濱150090；3.Department of Civil and Environmental Engineering，University of California，Berkeley，CA，the US，94710）

針對某煤礦井下救生艙，利用ANSYS和LS－DYNA軟件建立救生艙、瓦斯／瓦斯煤塵和巷道三維有限元模型。分別進(jìn)行了真實(shí)內(nèi)瓦斯／瓦斯煤塵爆炸作用下救生艙流固耦合數(shù)值分析，等效三角波爆炸沖擊荷載作用下救生艙動態(tài)模擬和靜水壓力荷載作用下救生艙的極限承載力和極限變形的數(shù)值計(jì)算，確定救生艙動／靜態(tài)承載能力變化規(guī)律和變形模式，給出瓦斯／瓦斯煤塵爆炸超壓與救生艙承載能力和變形模式的變化規(guī)律。研究結(jié)果表明，等效三角波爆炸沖擊荷載作用下該救生艙整體結(jié)構(gòu)的極限超壓明顯比流固耦合計(jì)算分析結(jié)果大，流固耦合計(jì)算分析結(jié)果更接近于實(shí)驗(yàn)結(jié)果。該救生艙滿足爆炸沖擊波荷載下的變形要求時，其在靜水壓荷載作用下的變形也能滿足。

爆炸力學(xué)；有限元模型；爆炸沖擊；瓦斯煤塵；救生艙；巷道

隨著煤礦產(chǎn)業(yè)日益發(fā)達(dá)，礦井安全越發(fā)引人關(guān)注。每年的礦井事故均造成大量人員傷亡，其中絕大多數(shù)是由于災(zāi)難所形成的有毒煙氣，導(dǎo)致井下人員中毒、窒息而死亡［1］。據(jù)不完全統(tǒng)計(jì)，2001～2010年我國每生產(chǎn)100萬噸煤炭就有3.1個礦工付出生命［2］，礦難給我國造成了不可挽救的人員傷亡和巨大的經(jīng)濟(jì)損失［3］。

隨著對煤礦瓦斯爆炸事故發(fā)生的原因、爆炸傳播的過程及影響因素等方面的認(rèn)識越來越深入，在理論研究方面，D.Bjerketvedt等［4］、A.Grillo等［5］明確了瓦斯爆炸發(fā)生的基本條件，并給出了瓦斯爆炸點(diǎn)火過程的一套完整的基元反應(yīng)表達(dá)式；A.D.Benedetto等［6］、Y.J.Zhu等［7］深入研究了瓦斯爆炸的點(diǎn)火階段和傳播階段，得到了最小點(diǎn)火能量、爆炸參數(shù)、熱力學(xué)參數(shù)等一系列重要的結(jié)果；Y.Inaba等［8］建立了爆轟的經(jīng)典C－J理論，解釋了發(fā)生爆燃和爆轟的條件，廣泛應(yīng)用于實(shí)際工程中。

隨著計(jì)算流體力學(xué)和數(shù)值計(jì)算方法開始應(yīng)用于瓦斯爆炸方面的研究，A.Kobiera等［8］提出了描述管道內(nèi)氣體爆炸過程的數(shù)學(xué)模型，并針對乙炔／空氣的混合氣體爆炸問題進(jìn)行了研究，采用二維Navier－Stokes方程對爆炸過程進(jìn)行了數(shù)值模擬。王志榮［10］、李江濤［11］、胡春明［12］、Jiang Bingyou等［13－14］主要利用Fluent、AutoReaGas等有限元軟件研究了爆炸沖擊波的傳播特性和動力特性。

針對礦用救生艙在爆炸沖擊波作用下的承載力、動力響應(yīng)等問題，以往的研究中，為了節(jié)省時間，簡化計(jì)算，更多采用三角形或簡單的沖擊波形式對救生艙直接加載進(jìn)行強(qiáng)度分析。李志強(qiáng)等［15］、常德功等［16］、曾一鑫等［17］利用有限元軟件，分析了救生艙在等效三角沖擊波作用下的動力響應(yīng)。該方法雖然計(jì)算簡便，但存在失真問題，會造成救生艙動力響應(yīng)分析的誤差。劉超等［18］、王磊等［19］、王云艷等［20］、白博等［21］結(jié)合流固耦合原理，采用AUTODYN、LS－DYNA等有限元軟件對爆炸沖擊波下救生艙的受力性能進(jìn)行了分析。徐景德等［22］、焦宇等［23］、司榮軍［24］利用實(shí)驗(yàn)巷道進(jìn)行了許多相關(guān)的瓦斯煤塵爆炸研究，并取得了一定的研究成果。使用巷道來進(jìn)行瓦斯煤塵爆炸實(shí)驗(yàn)的耗資比較大，其適用性也較為不足，因此，K.Lebecki等［25］、J.W.Boh等［26］、P.E.Moore等［27］、W.Houf等［28］、浦以康等［29］、董磊［30］開始采用小尺寸管道內(nèi)瓦斯煤塵爆炸實(shí)驗(yàn)來進(jìn)行瓦斯煤塵爆炸方面的研究。

本文中針對某煤礦井下救生艙，利用大型有限元分析軟件建立救生艙、瓦斯／瓦斯煤塵和巷道三維有限元模型。分別討論3種荷載形式作用下的救生艙極限承載力和變形，包括：真實(shí)內(nèi)瓦斯／瓦斯煤塵爆炸作用下救生艙流固耦合數(shù)值分析，等效三角波爆炸沖擊荷載作用下救生艙動態(tài)模擬和靜水壓力荷載作用下救生艙靜態(tài)數(shù)值計(jì)算。通過以上3種形式荷載的多種工況數(shù)值模擬，確定救生艙動／靜態(tài)承載能力變化規(guī)律和變形模式，給出瓦斯／瓦斯煤塵爆炸超壓與救生艙承載能力和變形模式的變化規(guī)律，可為救生艙在大型實(shí)驗(yàn)巷道進(jìn)行足尺實(shí)驗(yàn)提供理論指導(dǎo)。

1 實(shí) 驗(yàn)

1.1 實(shí)驗(yàn)系統(tǒng)

大型實(shí)驗(yàn)巷道如圖1所示［31］，巷道全長896m（其中主巷長710m），主巷斷面積7.2m2，副巷斷面積6.2m2。救生艙長度8m，高度1.72m，寬度1.4m，艙體放置在巷道中橫斷面示意圖如圖2所示。

圖1 大型實(shí)驗(yàn)巷道及斷面示意圖Fig.1 Sketches of system section of large－scale experiment laneway

圖2 艙體在巷道中橫斷面示意圖Fig.2 Sketch of refuge chamber section in large－scale experiment laneway

1.2 實(shí)驗(yàn)過程及結(jié)果

實(shí)驗(yàn)分2個部分進(jìn)行。首先進(jìn)行瓦斯爆炸實(shí)驗(yàn)，救生艙安裝位置距爆源為40m，該位置進(jìn)行2次實(shí)驗(yàn)。接著進(jìn)行瓦斯／煤塵爆炸實(shí)驗(yàn)，救生艙在爆炸荷載作用下，沖擊超壓達(dá)到0.375MPa，其主要變形及破壞情況，如圖3所示。其中，在艙體中部側(cè)面，目測最大相對變形約20cm；在梁柱相交處，出現(xiàn)較大的裂縫，梁柱之間斷開；在柱與加勁肋相交處，出現(xiàn)較大的屈曲，但并未出現(xiàn)裂縫；2節(jié)艙體連接處的法蘭，出現(xiàn)較大的裂縫。

圖3 救生艙重點(diǎn)部位變形圖Fig.3 Deformation pattern refuge chamber’s key part

2 流固耦合數(shù)值分析

2.1 流固耦合數(shù)值分析模型

根據(jù)救生艙以及巷道的實(shí)際形狀和尺寸，建立有限元模型如圖4所示。救生艙放在巷道中，起爆點(diǎn)距離艙體分別40和80m。利用ANSYS／LS－DYNA進(jìn)行數(shù)值計(jì)算，采用流固耦合的計(jì)算方法考慮爆炸空氣沖擊波與結(jié)構(gòu)的相互作用。

圖4 救生艙及巷道有限元模型示意圖Fig.4 Finite element sketch of refuge chamber and laneway

2.2 流固耦合數(shù)值分析結(jié)果

根據(jù)司榮軍［24］的研究結(jié)果結(jié)合本文中所用救生艙的構(gòu)造設(shè)計(jì)，通過調(diào)節(jié)瓦斯的體積分?jǐn)?shù)、煤塵的體積質(zhì)量、起爆點(diǎn)距艙體的距離以及瓦斯空氣混合氣體的體積，從而得到不同的沖擊波超壓。總共分3個工況，如表1所示，表中：d為艙體距爆點(diǎn)距離，φ為瓦斯體積分?jǐn)?shù)，ρ為煤塵體積質(zhì)量，L為混合氣體長度，V為混合氣體體積，p為沖擊超壓，τ為作用時間。工況1：起爆點(diǎn)距艙體40m，瓦斯空氣混合氣體體積為50m3，瓦斯體積分?jǐn)?shù)為4.5%，得到的超壓時程曲線如圖5所示；工況2：起爆點(diǎn)距艙體40m，瓦斯空氣混合氣體體積100m3，瓦斯體積分?jǐn)?shù)分別為3%、4.5%、6%、7%和9.5%，得到的超壓時程曲線如圖6所示；工況3：起爆點(diǎn)距艙體80m，瓦斯空氣混合氣體體積為200m3，瓦斯體積分?jǐn)?shù)為9.5%，煤塵體積質(zhì)量分別為50、80、240、300和360g／m3，得到的超壓時程曲線如圖7所示。

圖5 工況1下沖擊超壓時程曲線Fig.5 History of shock waves overpressure in case 1

圖6 工況2下沖擊超壓時程曲線Fig.6 Histories of shock waves overpressure in case 2

圖7 工況3下沖擊超壓時程曲線Fig.7 Histories of shock waves overpressure in case 3

表1 不同工況下瓦斯／瓦斯煤塵的參數(shù)設(shè)置Table 1 Parameters of gas／gas and grime in different cases

通過上述各種工況的數(shù)值分析，計(jì)算出作用在結(jié)構(gòu)上的沖擊超壓為0.05～0.435MPa，則各工況下救生艙應(yīng)力、應(yīng)變和變形的分布情況如表2所示。

表2 瓦斯／瓦斯煤塵作用下救生艙動態(tài)響應(yīng)Table 2 Dynamic response of refuge chamber under the action of gas／gas and grime

從上述分析結(jié)果可以看出，在模擬真實(shí)巷道內(nèi)瓦斯／瓦斯煤塵爆炸沖擊作用下，該救生艙保證變形不大的條件下，所能承受的最大超壓上限不超過0.05MPa。通過對比0.385 0MPa超壓作用下的模擬結(jié)果和實(shí)驗(yàn)結(jié)果，救生艙均已發(fā)生嚴(yán)重變形，最大相對變形達(dá)到約20cm，并且部分橫梁和縱梁斷裂，如此可以看出，流固耦合數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果擬合的較好，該模擬結(jié)果有一定的可靠性。

3 等效三角波動態(tài)數(shù)值分析

3.1 等效三角波模型

等效三角波動態(tài)數(shù)值分析所用的救生艙有限元模型與流固耦合動態(tài)分析中所用的模型一樣，只是荷載種類和加載方式不同。利用等效三角波模型近似模擬瓦斯煤塵爆炸沖擊波荷載，僅考慮救生艙在開放空間下等效三角波作用的動態(tài)響應(yīng)，能夠較快得到結(jié)構(gòu)動態(tài)響應(yīng)結(jié)果。簡化的典型三角斜坡沖擊荷載時程曲線如圖8所示。實(shí)驗(yàn)的加載方式為：對于整個六面體救生艙，除底面外，其他5個面均加載超壓和作用時間均相同的等效三角波，每加載1次，作為1個加載工況，根據(jù)所加載的沖擊超壓的不同，一共設(shè)置了15個不同的工況，對于各工況的分析，如表3所示。

圖8 等效三角波荷載Fig.8 Equivalent triangular wave load

表3 不同等效三角波荷載作用下工況及響應(yīng)Table 3 Loading cases and dynamic responses of structures under different equivalent triangular wave loads

3.2 等效三角波動態(tài)數(shù)值分析過程及結(jié)果

整個結(jié)構(gòu)z方向的最大變形出現(xiàn)在前門門扇處，y方向的最大變形出現(xiàn)在過渡艙頂板處，x方向的最大變形出現(xiàn)在過渡艙側(cè)板處。典型殼單元取過渡艙部分圓弧過渡區(qū)，典型梁單元取過渡艙部分縱梁處。將上述15種等效三角波加載工況的結(jié)構(gòu)響應(yīng)在表3中匯總?？梢钥闯觯摼壬摑M足剛度要求（最大變形小于3cm）彈性響應(yīng)極限沖擊超壓（7ms作用時間）為0.13MPa，隨作用時間的增加，沖擊超壓要降低，保守值為0.05MPa。根據(jù)表3中對于在0.4MPa沖擊超壓作用下艙體狀態(tài)的描述，與實(shí)驗(yàn)結(jié)果進(jìn)行對比，在模擬荷載略大于實(shí)驗(yàn)荷載的情況下，模擬結(jié)果最大變形為8.3cm，實(shí)驗(yàn)結(jié)果最大變形約為20cm，模擬得到的變形結(jié)果明顯偏小。

表3中所示15種工況中超壓荷載大多小于0.7MPa，作用時間大多為7ms。將救生艙體的響應(yīng)結(jié)果按照典型節(jié)點(diǎn)處沖擊超壓所對應(yīng)最大變形、最大等效塑性應(yīng)變、最大等效應(yīng)力分別繪制變化曲線，如圖9～10所示。對圖9進(jìn)行分析：從變形角度來看，救生艙整體3個方向以側(cè)板沿x方向變形最大，若要求最大變形小于3cm，沖擊超壓的極限應(yīng)小于0.13MPa；對圖10進(jìn)行分析，從等效應(yīng)變來看，Q235結(jié)構(gòu)鋼的救生艙骨架（梁單元）和Q345結(jié)構(gòu)鋼的救生艙蒙皮（殼單元）彈性極限應(yīng)變對應(yīng)的超壓分別為0.25和0.5MPa；從等效應(yīng)力來看，救生艙骨架和蒙皮彈性極限應(yīng)力對應(yīng)的超壓分別為0.50和0.27 MPa。梁與蒙皮的強(qiáng)度不匹配，結(jié)構(gòu)極限承載力由強(qiáng)度低的部件控制。

圖9 整體最大變形與峰值超壓關(guān)系Fig.9 Relationship between structural maximum deformation and peak overpressure

圖10 效塑性應(yīng)變和應(yīng)力與峰值超壓關(guān)系Fig.10 Equivalent plastic strains and stress vs.peak pressures

4 靜水壓荷載靜態(tài)響應(yīng)數(shù)值分析

4.1 靜水壓荷載靜態(tài)響應(yīng)分析模型

瓦斯／煤塵爆炸沖擊作用下救生艙的數(shù)值模擬能夠反映結(jié)構(gòu)的動態(tài)響應(yīng)，但是，井下災(zāi)害除了爆炸，還會發(fā)生透水事故，此時救生艙承受靜水壓力作用在設(shè)計(jì)中也需要重點(diǎn)考慮。采用ANSYS有限元分析軟件，針對救生艙進(jìn)行了靜水壓力作用下的數(shù)值模擬計(jì)算。救生艙構(gòu)造形式和尺寸與前面相同，在逃生艙前、后、左、右、上等5個面施加均布靜水壓荷載，分13個工況加載，所對應(yīng)的靜水壓荷載分別為0.050、0.100、0.125、0.150、0.175、0.200、0.225、0.250、0.275、0.300、0.325、0.350和0.375MPa。

4.2 靜水壓荷載靜態(tài)響應(yīng)分析過程及結(jié)果

通過對13種工況的計(jì)算后，提取3個方向最大變形對應(yīng)的不同單元：（1）x方向（艙長度方向）的最大變形單元位于正面門的中部；（2）y方向（艙側(cè)面方向）的最大變形單元位于側(cè)面主艙中部；（3）z方向（艙高度方向）的最大變形單元位于頂面緩沖艙中部。得到3個方向上最大變形隨靜水壓荷載變化的曲線如圖11所示；提取3個方向上最大應(yīng)力的不同單元，得到3個方向上最大應(yīng)力隨靜水壓荷載變化的曲線如圖12所示；同樣地，提取3個方向上最大應(yīng)變的不同單元，得到3個方向上最大應(yīng)變隨靜水壓荷載變化的曲線如圖13所示。由圖中可以看出，當(dāng)荷載小于0.17MPa時，危險(xiǎn)點(diǎn)最大變形小于3cm，響應(yīng)的應(yīng)力響應(yīng)小于250MPa，應(yīng)變小于0.2%，結(jié)構(gòu)基本處于彈性狀態(tài)。

對整體結(jié)構(gòu)的應(yīng)力響應(yīng)進(jìn)行分析，根據(jù)荷載作用下的應(yīng)力分布情況，選擇5個代表性的單元（應(yīng)力較大區(qū)域的單元）提取數(shù)據(jù)，典型單元位置如圖14所示。繪制最大應(yīng)力、最大應(yīng)變隨荷載等級變化的曲線，如圖15～16所示。其中，36821為救生艙頂部弧面上最大應(yīng)力對應(yīng)的單元，25522為救生艙頂部與正面交接處最大應(yīng)力對應(yīng)的單元，18825為救生艙側(cè)面緩沖艙與主艙交接處最大應(yīng)力對應(yīng)的單元，17016為救生艙側(cè)面與正面和地面交接處最大應(yīng)力對應(yīng)的單元，15061為救生艙側(cè)面中間部位處最大應(yīng)力對應(yīng)的單元。

圖11 不同單元對應(yīng)的最大變形與靜水壓載荷關(guān)系Fig.11 Relationship between structural maximum deformation of different units and hydrostatic pressure load

通過上述分析可以知道，出現(xiàn)塑性最早的單元一般都在艙體與梁交接處，對應(yīng)的彈性極限靜壓在0.1MPa左右。將上述13種加載工況的結(jié)構(gòu)響應(yīng)匯總，可以得到，該救生艙滿足剛度要求（最大變形小于3cm）彈性響應(yīng)極限靜水壓為0.173MPa，滿足強(qiáng)度要求的彈性響應(yīng)極限靜水壓為0.1MPa。

圖12 不同單元對于的最大應(yīng)力與靜水壓載荷關(guān)系Fig.12 Relationship between structural maximum stress of different units and hydrostatic pressure load

圖13 不同單元對于的最大應(yīng)變與靜水壓載荷關(guān)系Fig.13 Relationship between structural maximum strain of different units and hydrostatic pressure load

圖14 典型單元位置示意圖Fig.14 Position of representative elements

圖15 典型單元最大應(yīng)力和靜水壓載荷關(guān)系Fig.15 Maximum stress varying with hydrostatic pressure load of representative units

圖16 典型單元最大應(yīng)變和靜水壓載荷關(guān)系Fig.16 Structural maximum strain varying with hydrostatic pressure load of representative units

5 結(jié) 論

結(jié)合救生艙爆炸實(shí)驗(yàn)，利用有限元分析軟件ANSYS／LS－DYNA對國產(chǎn)某救生艙艙體進(jìn)行流固耦合、等效三角波和靜水壓荷載作用下數(shù)值模擬分析，通過對救生艙結(jié)構(gòu)在不同荷載作用模式下的承載力和變形極限狀態(tài)的對比分析，可以得到以下結(jié)論：

（1）流固耦合分析結(jié)果表明，在保證該救生艙變形不大的情況下，所能承受的動態(tài)沖擊超壓上限為0.05MPa，與實(shí)驗(yàn)結(jié)果最為接近，表明該種方法的分析結(jié)果較為可靠。采用簡化的等效三角波荷載作用下，整體結(jié)構(gòu)沖擊超壓最小值為0.13MPa，與流固耦合分析結(jié)果相比較，該超壓明顯偏大。

（2）雖然流固耦合動態(tài)分析方法能夠比較真實(shí)地模擬救生艙在井下巷道受瓦斯爆炸沖擊作用下的響應(yīng)特性，但是其建模和分析相對比較復(fù)雜、耗時。通過調(diào)整各面等效三角波荷載的沖擊超壓、到達(dá)時間和作用時長，能夠很好的模擬救生艙在爆炸沖擊波作用下的動態(tài)響應(yīng)，這種方法會更加高效。

（3）對于救生艙結(jié)構(gòu)受靜水壓荷載，其滿足剛度要求的彈性響應(yīng)極限靜水壓值要比動態(tài)響應(yīng)分析中的沖擊超壓值大，因此，如果能保證救生艙結(jié)構(gòu)在爆炸沖擊波荷載下的變形要求，則其在靜水壓荷載作用下的變形也能滿足要求。

（4）不同荷載模式作用下，救生艙失效部位主要有中部艙體位置、艙門迎爆面位置以及艙體與梁交接位置，設(shè)計(jì)時應(yīng)對這些部位予以加強(qiáng)。從變形、等效應(yīng)變、等效應(yīng)力角度來看，救生艙的失效狀態(tài)是不同的，在保證梁與蒙皮強(qiáng)度匹配的條件下，根據(jù)不同要求，選擇合適的加強(qiáng)方案。

［1］楊俊玲，馬躍征，楊魯偉，等.煤礦救生艙防爆性能分析［J］.煤炭學(xué)報(bào)，2013，38（增1）：159－163.Yang Junling，Ma Yuezheng，Yang Luwei，et al.Explosion－proof performance analysis of mine refuge chamber［J］.Journal of China Coal Society，2013，38（Suppl 1）：159－163.

［2］姚文亮.中國礦難的政府規(guī)制研究［D］.西安：西北大學(xué)，2011.

［3］梅瑞斌，李長生，蔡般，等.爆炸沖擊下救生艙抗爆能力有限元分析［J］.東北大學(xué)學(xué)報(bào)（自然科學(xué)版），2013，34（1）：85－94.Mei Ruibin，Li Changsheng，Cai Ban，et al.FEM analysis of anti－deformation capability for coal mine refuge chamber suffered to gas explosion［J］.Journal of Northeastern University（Natural Science），2013，34（1）：85－94.

［4］Bjerketvedt D，Bakke J R，Wingerden K V.Gas explosion handbook［J］.Journal of Hazardous Materials，1997，52（1）：1－150.

［5］Grillo A，Slack M W.Shock tube study of ignition delay times in methane oxygen nitrogen argon mixtures［J］.Combustion ＆Flame，1976，27（3）：377－381.

［6］Benedetto A D，Sarli V D，Salzano E，et al.Explosion behavior of CH4／O2／N2／CO2，and H2／O2／N2／CO2，mixtures［J］.International Journal of Hydrogen Energy，2009，34（16）：6970－6978.

［7］Zhu Y J，Chao J，Lee J H S.An experimental investigation of the propagation mechanism of critical deflagration waves that lead to the onset of detonation［J］.Proceedings of the Combustion Institute，2007，31（2）：2455－2462.

［8］Inaba Y，Nishihara T，Groethe M A，et al.Study on explosion characteristics of natural gas and methane in semiopen space for the HTTR hydrogen production system［J］.Nuclear Engineering ＆Design，2004，232（1）：111－119.

［9］Kobiera A，Kindracki J，Zydak P，et al.A new phenomenological model of gas explosion based on characteristics of flame surface［J］.Journal of Loss Prevention in the Process Industries，2007，20（3）：271－280.

［10］王志榮.受限空間氣體爆炸傳播及其動力學(xué)過程研究［D］.南京：南京工業(yè)大學(xué)，2005.

［11］李江濤.微型圓管結(jié)構(gòu)對火焰?zhèn)鞑ヒ种谱饔玫臄?shù)值模擬［D］.大連：大連理工大學(xué)，2006.

［12］胡春明.預(yù)混火焰在平板狹縫中傳播與淬熄的研究［D］.大連：大連理工大學(xué)，2006.

［13］Jiang Bingyou，Lin Baiquan，Shi Shulei，et al.Numerical simulation on the influences of initial temperature and initial pressure on attenuation characteristics and safety distance of gas explosion［J］.Combustion Science and Technology，2012，184（2）：135－150.

［14］江丙友，林柏泉，朱傳杰，等.瓦斯爆炸沖擊波在并聯(lián)巷道中傳播特性的數(shù)值模擬［J］.燃燒科學(xué)與技術(shù)，2011，17（3）：250－254.Jiang Bingyou，Lin Baiquan，Zhu Chuanjie，et al.Numerical simulation on shock wave propagation characteristics of gas explosion in parallel roadway［J］.Journal of Combustion Science ＆Technology，2011，17（3）：250－254.

［15］李志強(qiáng)，白博，謝青海，等.沖擊載荷下礦用移動式救生艙動態(tài)響應(yīng)的數(shù)值模擬［J］.振動與沖擊，2013，32（16）：146－151.Li Zhiqiang，Bai Bo，Xie Qinghai，et al.Numerical simulation on the dynamic response of movable mine rescue chamber subjected to impact loading［J］.Journal of Vibration ＆Shock，2013，32（16）：146－151.

［16］常德功，王吉利，李國星.基于LS－DYNA的礦用救生艙殼結(jié)構(gòu)爆炸沖擊分析［J］.礦山機(jī)械，2013，41（11）：130－134.Chang Degong，Wang Jili，Li Guoxing.Explosion impact analysis on shell of mine refuge chamber based on LSDYNA［J］.Mining ＆Processing Equipment，2013，41（11）：130－134.

［17］曾一鑫，白春華，李建平，等.巷道救生艙抗沖擊數(shù)值仿真［J］.煤炭學(xué)報(bào)，2012，37（10）：1705－1708.Ceng Yixin，Bai Chunhua，Li Jianping，et al.Numerical simulation of rescue cabin under blasting in the tunnel［J］.Journal of the China Coal Society，2012，37（10）：1705－1708.

［18］劉超，鮑久圣.救生艙抗爆炸沖擊載荷數(shù)值模擬［J］.金屬礦山，2014，43（1）：141－144.Liu Chao，Bao Jiusheng.Numerical Simulation of explosion impact load on refuge chamber［J］.Metal Mine，2014，43（1）：141－144.

［19］王磊，范世平，龔建宇，等.礦井KJYF－96／8型救生艙抗爆有限元分析［J］.煤炭科學(xué)技術(shù)，2014，42（3）：69－72.Wang Lei，F(xiàn)an Shiping，Gong Jianyu，et al.Finite element analysis of antiknock performance for KJYF－96／8refuge chamber［J］.Coal Science ＆Technology，2014，42（3）：69－72.

［20］王云艷，何寧.基于AUTODYN、LS－DYNA的救生艙抗爆仿真試驗(yàn)［J］.安全與環(huán)境學(xué)報(bào)，2014，14（2）：9－12.Wang Yunyan，He Ning.Simulation test of rescue chamber under the explosive load based on the AUTODYN and LS－DYNA model［J］.Journal of Safety ＆Environment，2014，14（2）：9－12.

［21］白博，李志強(qiáng)，王志華，等.礦用救生艙抗爆炸TNT當(dāng)量與沖擊載荷研究［J］.科學(xué)技術(shù)與工程，2014，14（17）：1－5.Bai Bo，Li Zhiqiang，Wang Zhihua，et al.Research on the TNT equivalent and impact load in anti－explosion of mine rescue capsule［J］.Science Technology ＆Engineering，2014，14（17）：1－5.

［22］徐景德，徐勝利，楊庚宇.礦井瓦斯爆炸傳播的試驗(yàn)研究［J］.煤炭科學(xué)技術(shù)，2004，33（7）：55－57.Xu Jingde，Xu Shengli，Yang Gengyu.Experimental study on mine gas explosion diffusion［J］.Coal Science ＆Technology，2004，33（7）：55－57.

［23］焦宇，周心權(quán)，段玉龍，等.瓦斯爆炸煙流濃度和溫度的擴(kuò)散規(guī)律［J］.煤炭學(xué)報(bào)，2011，36（2）：293－297.Jiao Yu，Zhou Xinquan，Duan Yulong，et al.Concentration and temperature diffusion laws of smoke plume after methane explosions［J］.Journal of the China Coal Society，2011，36（2）：293－297.

［24］司榮軍.礦井瓦斯煤塵爆炸傳播規(guī)律研究［D］.濟(jì)南：山東科技大學(xué)，2007.

［25］Lebecki K，Sliz J，Cybulski K，et al.Efficiency of triggered barriers in dust explosion suppression in galleries［J］.Journal of Loss Prevention in the Process Industries，2001，14（6）：489－494.

［26］Boh J W，Louca L A，Choo Y S.Numerical assessment of explosion resistant profiled barriers［J］.Marine Structures，2004，17（2）：139－160.

［27］Moore P E，Spring D J.Design of Explosion Isolation Barriers［J］.Process Safety ＆Environmental Protection，2005，83（2）：161－170.

［28］Houf W，Schefer R，Evans G，et al.Evaluation of barrier walls for mitigation of unintended releases of hydrogen［J］.International Journal of Hydrogen Energy，2011，36（36）：2520－2529.

［29］浦以康，胡山.粉塵火焰加速現(xiàn)象的實(shí)驗(yàn)研究［J］.爆炸與沖擊，1995，15（2）：97－106.Pu Yikang，Hu Shan.Experimental studies on the phenomena of dust flame acceleration［J］.Explosion and Shock Waves，1995，15（2）：97－106.

［30］董磊.障礙物體形對爆炸波傳播的影響［D］.天津：天津大學(xué)，2008.

［31］樊小濤.礦用救生艙抗爆性能試驗(yàn)研究［J］.礦業(yè)安全與環(huán)保，2010，37（3）：25－30.Fan Xiaotao.Experiment study on the explosion－proof performance of mine refuge chamber［J］.Mining Safety and Environmental Protection，2010，37（3）：25－30.

Simulation analysis of mine refuge chamber performance in different loading modes

Li Jilong1，2，3，Tang Ya’nan2，Liu Xuanming2
（1.Key Lab of Structures Dynamic Behavior and Control，Harbin Institute of Technology，Harbin 150090，Heilongjiang，China；2.School of Civil Engineering，Harbin Institute of Technology，Harbin 150090，Heilongjiang，China；3.Department of Civil and Environmental Engineering，University of California，Berkeley，CA，the US，94710）

In this work，by using the ANSYS and LS－DYNA softwares，we established the three－dimensional finite element model of the refuge chamber，the gas／gas and grime and the laneway.The ultimate strength and deformation of the refuge chamber were simulated under the action of three kinds of loads，mainly including the fluid－solid interaction numerical analysis of the real gas／gas and grime explosion effect，the equivalent triangle wave explosion numerical dynamic simulation and the hydrostatic pressure load calculations.Through the above three forms of numerical simulations of loading，we determined the dynamic／static load capacities of the refuge chamber as well as the deformation rules with the given gas／gas and coal dust explosion overpressure variation.The simulation results indicate that the limit overpressure of the refuge chamber under the equivalent triangle wave explosion is bigger than the fluid－solid interaction numerical analysis results，which approaches more closely to the experimental results.If the refuge chamber can meet the demand of transformation under explosion loading，it can also satisfy the transformation requirement under hydrostatic pressure load.

mechanics of explosion；finite elements models；blast impact；gas and grim；refuge chamber；laneway

O383；TD77.4國標(biāo)學(xué)科代碼：13035

A

10.11883／1001－1455（2017）01－0140－10

（責(zé)任編輯 王易難）

2015－05－04；

2015－10－27

李冀龍（1968— ），男，博士，副教授，leekeelung＠126.com。