煤礦井下作業(yè)經(jīng)常發(fā)生瓦斯、煤塵爆炸，產(chǎn)生大量一氧化碳等有毒氣體，直接造成安全事故[1]。礦用防爆耐壓密閉門（簡(jiǎn)稱防爆門）是規(guī)避事故的重要安全裝備，其抗爆性能已成為煤礦井下生產(chǎn)裝備的研究熱點(diǎn)之一，對(duì)提升我國煤礦安全生產(chǎn)保障水平等具有重要意義[2]。

在抗爆性能研究方面，Katherine A[3]提出了一種分析救生艙的抗爆性能的數(shù)值計(jì)算方法，通過將爆炸的動(dòng)載荷轉(zhuǎn)化為靜載荷的方式，均勻加載到救生艙的艙體上，并進(jìn)行艙體結(jié)構(gòu)的抗爆性能研究。Houf W G等[4]深入研究了瓦斯爆炸后艙體的變形情況，揭示了在瓦斯沖擊波作用下艙體的變形規(guī)律。劉建英等[5]將泡沫鋁夾芯板材料作為艙門，通過加強(qiáng)筋的固定作用強(qiáng)化艙體，通過數(shù)值計(jì)算了移動(dòng)式救生艙的強(qiáng)度要求與安全系數(shù)。宋勝偉等[6]采用ANSYS仿真軟件模擬艙體內(nèi)的自動(dòng)隔離門不同外界環(huán)境下的應(yīng)力應(yīng)變情況。吳廣明等[7]通過LS-DYNA有限元軟件仿真研究了艙室內(nèi)在爆載荷作用下鋪設(shè)凱夫拉防護(hù)材料發(fā)生的動(dòng)態(tài)響應(yīng)規(guī)律，建立了動(dòng)態(tài)響應(yīng)的有限元模型及艦艇艙室結(jié)構(gòu)模型。王文娟等[8]在有限元方法的基礎(chǔ)上通過結(jié)構(gòu)動(dòng)力分析的方式系統(tǒng)研究了煤礦的避難硐室中的密閉門在爆炸載荷作用下的受力及變形情況。

本文針對(duì)防爆門的設(shè)計(jì)和結(jié)構(gòu)模型，采用顯式非線性動(dòng)力分析方法，通過建立有限元仿真模型，數(shù)值分析防爆門抗爆性能響應(yīng)。

1.   防爆門的設(shè)計(jì)要求及結(jié)構(gòu)模型

在實(shí)際應(yīng)用中，防爆門采用一種整體式的應(yīng)用結(jié)構(gòu)，門體為密閉式，防爆門的門體結(jié)構(gòu)如圖1所示。防爆門的門體為矩形結(jié)構(gòu)，寬1120 mm，高2120 mm；防爆門的鉸鏈聯(lián)接板是門體的聯(lián)接零件，長(zhǎng)60 mm，厚225 mm；防爆門的門扇厚約200 mm，并且門扇的邊緣安裝有耐高溫的密封條；防爆門門體的內(nèi)腔尺寸約為1000 mm×2000 mm，防爆門門體的內(nèi)腔四角具有圓形外觀；加強(qiáng)筋是防爆門上關(guān)鍵的固定件，排布在門體的門扇內(nèi)腔中。防爆門的主要聯(lián)接材料是Q345鋼，門體結(jié)構(gòu)主要由Q345鋼的扁條焊接而成，門體中間及門框均采用Q345鋼的扁條聯(lián)接。

防爆門主要零件的材料均采用Q345鋼，如鉸鏈、加強(qiáng)筋、大齒輪、手輪結(jié)構(gòu)、蓋板結(jié)構(gòu)、門板焊件等零件；防爆門傳動(dòng)結(jié)構(gòu)采用Q235鋼來保證部件的剛性及穩(wěn)定性；門體中的把手結(jié)構(gòu)采用ZG45；門體中主要零件的聯(lián)接部件，如螺母、螺栓、銷釘、墊片、鋼蓋、螺母等均使用材料304不銹鋼；通過鉚釘鉚接304不銹鋼表面拉絲板作為外觀門板；防爆門的觀察結(jié)構(gòu)采用高強(qiáng)度的防爆抗壓玻璃；防火材料為MT113隔離抗靜電的膠條材料。建立防爆門的三維模型如圖2所示。在防爆門實(shí)際加工過程中，通過門栓將門框固定起來，門框通過鉸鏈來聯(lián)接，因此在數(shù)值仿真中將防爆門設(shè)定為固定約束[9]。

防爆門需要擁有一定的抗爆壓力才能滿足防護(hù)需求，至少能夠承受0.6 MPa的壓力，并且能夠承受瓦斯爆炸產(chǎn)生的三角沖擊波壓力300 ms，其中防爆門的整體結(jié)構(gòu)及關(guān)鍵部位承受的最大三角形沖擊波壓力為1 MPa。

2.   有限元仿真分析模型

防爆門的數(shù)值仿真分析主要考慮爆炸沖擊波作用的外界條件引起門體結(jié)構(gòu)的應(yīng)力破壞及應(yīng)變變形。防爆門門體的位移變形是非線性變化的，因此在數(shù)值仿真中采用ANSYS軟件中的顯式非線性動(dòng)力模塊——LS-DYNA分析模塊來進(jìn)行分析，研究防爆門在爆炸環(huán)境下的變形情況及抗爆性能。

2.1   仿真模型

2.1.1   單元類型的選擇

為了提高防爆門抗爆性能的模擬準(zhǔn)確性與理論計(jì)算求解的正確性，需要根據(jù)防爆門的實(shí)際尺寸進(jìn)行建模，對(duì)防爆門主體結(jié)構(gòu)進(jìn)行網(wǎng)格劃分的同時(shí)，需要合理簡(jiǎn)化門體的一些部件，便于研究防爆門在爆炸沖擊波下的形變情況，這便需要對(duì)仿真的單元類型進(jìn)行合理選擇[11]?？紤]到防爆門的結(jié)構(gòu)及計(jì)算需求，在ANSYS軟件的模塊程序LS-DYNA中，選擇實(shí)體單元與彈簧阻尼單元作為仿真單元類型。在單元類型的選擇中，主要針對(duì)防爆門的門體結(jié)構(gòu)進(jìn)行仿真處理，采用軟件中的SOLID 164實(shí)體單元對(duì)門體結(jié)構(gòu)進(jìn)行整體的網(wǎng)格劃分，采用軟件中的COMBI 165彈簧阻尼單元對(duì)門體的銷釘、螺栓等聯(lián)接處進(jìn)行局部的網(wǎng)格劃分。

2.1.2   仿真網(wǎng)格劃分

防爆門門體的結(jié)構(gòu)特征較為復(fù)雜，需要對(duì)不同的表征參數(shù)進(jìn)行多次實(shí)驗(yàn)研究，數(shù)值仿真也需要更精細(xì)的參數(shù)設(shè)置，本次仿真的網(wǎng)格劃分參數(shù)如下：門體的劃分的最大單元為25 mm，門體上零件聯(lián)接處的細(xì)節(jié)部分劃分的最大單元為15 mm，其比例系數(shù)設(shè)定為默認(rèn)的參數(shù)1。設(shè)置好仿真參數(shù)后，對(duì)門體進(jìn)行網(wǎng)格劃分，劃分的單元共135766個(gè)，如圖3所示。

2.2   材料關(guān)系模型與特性參數(shù)

防爆門在沖擊波作用下會(huì)產(chǎn)生彈塑性變形。對(duì)防爆門材料而言，重點(diǎn)是考慮與穩(wěn)定有關(guān)的彈塑性本構(gòu)模型及材料的應(yīng)變率，有限元分析時(shí)選用LS-DYNA中材料模型Johnson-Cook[12]。

根據(jù)防爆門設(shè)計(jì)中各零部件的材料選擇，防爆門的主要材料為Q345鋼、Q235鋼、ZG45鋼等，材料的特征參數(shù)如表1所示。

2.3   仿真的邊界條件及初始條件

防爆門的結(jié)構(gòu)較為穩(wěn)定，不存在相對(duì)松動(dòng)的零部件，設(shè)定防爆門的初始條件下穩(wěn)定可靠。門體各結(jié)構(gòu)中沒有焊接應(yīng)力存在，焊接部位對(duì)門體結(jié)構(gòu)無影響；防爆門各零件間沒有安裝變形，螺栓等聯(lián)接件的預(yù)應(yīng)力不影響門體結(jié)構(gòu)特征。

2.3.1   邊界條件

設(shè)定防爆門在承受爆炸條件下的沖擊波時(shí)，沖擊波的沖擊壓力直接作用于迎爆面上，門框上的鉸鏈及其他零部件沒有過定義接觸情況的發(fā)生，在模擬條件下，設(shè)定門栓、門框固定，門框邊界無其他壓力作用效果，防爆門可以簡(jiǎn)化成固定約束，在仿真模擬條件下，通過ANSYS軟件對(duì)防爆門的邊界施加固定約束。

2.3.2   初始條件

通過防爆門的結(jié)構(gòu)模型可知，門體的整體結(jié)構(gòu)受到爆炸沖擊波后，防爆門接觸沖擊波的外邊界區(qū)域發(fā)生了彈塑性變形，變形的區(qū)域相對(duì)防爆門的整體結(jié)構(gòu)幅度非常小，對(duì)網(wǎng)格劃分的精度無明顯影響。數(shù)值仿真選擇拉格朗日算法，該方法能夠直觀處理物質(zhì)界面并能精確描述材料邊界與物質(zhì)之間的相互關(guān)系[13]。

防爆門在爆炸環(huán)境的作用下可能會(huì)受到?jīng)_擊波的直接破壞，為了減小誤差，在數(shù)值仿真中將防爆門的正面作為迎爆面，并將峰值沖擊波加載到迎爆面上，考慮到爆炸環(huán)境的安全系數(shù)及誤差，本次仿真爆炸沖擊波的最小峰值壓力為0.6 MPa，數(shù)值計(jì)算中將沖擊波的最大峰值壓力為1.0 MPa，沖擊波壓力與時(shí)間的關(guān)系曲線如圖4所示。

3.   防爆門整體結(jié)構(gòu)數(shù)值分析

當(dāng)爆炸沖擊波作用在門體結(jié)構(gòu)上時(shí)，設(shè)定門體結(jié)構(gòu)所受的沖擊壓力持續(xù)300 ms，防爆門的壓力峰值為0.6, 1 MPa，仿真中的計(jì)算時(shí)間為400 ms。通過ANSYS/LS-DYNA進(jìn)行仿真分析，通過數(shù)值仿真分析防爆門結(jié)構(gòu)所受的應(yīng)力、應(yīng)變情況，由仿真圖中的應(yīng)力、應(yīng)變?cè)茍D體現(xiàn)。為了保證防爆門良好的剛性結(jié)構(gòu)，梁柱的最大變形量不能超過10 mm，即所受應(yīng)力小于屈服應(yīng)力。為防止密封失效，有密封要求的聯(lián)接件相對(duì)位移應(yīng)小于1 mm。

3.1   防爆門整體應(yīng)力分析

根據(jù)防爆門整體網(wǎng)格劃分，仿真計(jì)算后提取應(yīng)力結(jié)果，如圖5所示。當(dāng)迎爆面的峰值壓力達(dá)到0.6, 1 MPa時(shí)，防爆門受到的應(yīng)力最大值主要出現(xiàn)在門體的各聯(lián)接位置，計(jì)算時(shí)間為350 ms時(shí)，迎爆面所受的沖擊壓力達(dá)到了339.3 MPa?？梢?，在沖擊波持續(xù)300 ms的時(shí)間內(nèi)，門體結(jié)構(gòu)承受了兩種載荷與沖擊波的壓力。從圖中可知，門體結(jié)構(gòu)所受應(yīng)力沒有超過材料的屈服極限，符合材料的強(qiáng)度需求。

3.2   等效塑性應(yīng)變分析

防爆門門體在兩種外界載荷的作用下，門體加強(qiáng)筋部分第38332單元及第150343單元應(yīng)力達(dá)到了最大值，應(yīng)力的變化曲線如圖6所示。

塑性應(yīng)變最大值出現(xiàn)在最后時(shí)刻，提取該時(shí)刻在波峰為0.6, 1.0 MPa兩種載荷下塑性應(yīng)變?cè)茍D如圖7所示，防爆門門體整體未出現(xiàn)塑性應(yīng)變。

3.3   門體位移情況分析

防爆門的門體位移形變情況如圖8所示，當(dāng)門體在0.6, 1 MPa載荷的作用下，位移形變的峰值出現(xiàn)在防爆門的中間區(qū)域位置，位移變形量峰值分別達(dá)到1.706, 3.461 mm。可見，防爆門在承受波峰為0.6, 1.0 MPa，持續(xù)300 ms三角形沖擊波時(shí)，防爆門的門體的最大位移變形量不超過10 mm。對(duì)其進(jìn)行可靠性分析，防爆門各零件無裂縫及破壞現(xiàn)象產(chǎn)生，符合防爆門的剛性要求，能夠保證防爆門安全的防護(hù)性能。

由防爆門整體位移場(chǎng)分析可知，防爆門的位移變形量最大時(shí)刻門體最大位移值在第13588、第23769兩個(gè)節(jié)點(diǎn)處取得，分別提取兩個(gè)節(jié)點(diǎn)處的位移與時(shí)間關(guān)系曲線即δ-t曲線，如圖9所示。

防爆門門板聯(lián)接處的密封位置設(shè)置有3 mm厚的密封條，當(dāng)防爆門在作用時(shí)間300 ms內(nèi)所受的沖擊壓力峰值為0.6 MPa時(shí)，防爆門顯著點(diǎn)的位移變化如圖10所示?？梢钥闯觯辣T上的聯(lián)接零部件的最大位移變化數(shù)值為0.02 mm，出現(xiàn)在防爆門側(cè)方中部的聯(lián)接位置處。防爆門門板聯(lián)接處的密封位置最大位移未超出1 mm，因此不會(huì)發(fā)生密封失效。

當(dāng)防爆門門體在300 ms的作用時(shí)間內(nèi)所承受的沖擊壓力峰值為1.0 MPa，提取門體各聯(lián)接零部件顯著點(diǎn)的相對(duì)位移數(shù)據(jù)（絕對(duì)值）如表2所示。聯(lián)接件的相對(duì)位移最大值為0.004 mm出現(xiàn)在防爆門頂邊中間位置的聯(lián)接處，未超出1 mm，因此也不會(huì)發(fā)生密封失效。

4.   防爆門內(nèi)腔結(jié)構(gòu)數(shù)值分析

4.1   防爆門內(nèi)腔結(jié)構(gòu)應(yīng)力分析

在0.6, 1.0 MPa兩種波峰載荷下，防爆門門體所受的最大應(yīng)力發(fā)生在門體中間加強(qiáng)筋聯(lián)接區(qū)域。當(dāng)沖擊壓力持續(xù)時(shí)間達(dá)到70 ms時(shí)，最大應(yīng)力達(dá)到169.6 MPa；持續(xù)時(shí)間達(dá)到210 ms時(shí)的最大應(yīng)力為339.3 MPa?？梢姺辣T在承受兩種載荷持續(xù)300 ms的三角形沖擊波壓力下，應(yīng)力未超過材料屈服強(qiáng)度，能滿足強(qiáng)度要求。

4.2   防爆門內(nèi)腔結(jié)構(gòu)位移分析

在0.6, 1.0 MPa兩種波峰載荷下，防爆門腔內(nèi)零部件的結(jié)構(gòu)位移形變?nèi)?/span>圖11所示，門體腔內(nèi)零件結(jié)構(gòu)的中間區(qū)域位置的位移變化達(dá)到了最大值。作用時(shí)間在70 ms時(shí)，位移變化達(dá)到了1.706 mm，作用時(shí)間在210 ms時(shí)，位移變化達(dá)到了3.461 mm。從圖中可知，當(dāng)波峰載荷為0.6 MPa時(shí)，防爆門腔內(nèi)零件結(jié)構(gòu)的位移變形范圍在0~1.706 mm；當(dāng)波峰載荷為1 MPa時(shí)，防爆門腔內(nèi)零件結(jié)構(gòu)的位移變形范圍在0~3.461 mm，門體的內(nèi)腔結(jié)構(gòu)位移變化小于10 mm，沒有出現(xiàn)破壞現(xiàn)象，符合材料的強(qiáng)度需求。當(dāng)爆炸沖擊壓力作用在防爆門的門體結(jié)構(gòu)上時(shí)，門體的位移變形極小，門體的密封性能不受其影響，不影響防爆門安全的使用性能，內(nèi)部結(jié)構(gòu)不會(huì)發(fā)生破壞失效，滿足所需的剛度要求。

5.   防爆門抗爆性能數(shù)值分析結(jié)果

防爆門在抗爆沖擊波最大超壓峰值為0.6, 1.0 MPa 載荷下，數(shù)值仿真分析統(tǒng)計(jì)結(jié)果如表3所示。

根據(jù)表3，沖擊波波峰壓力為0.6 MPa防爆門整體結(jié)構(gòu)所受應(yīng)力的最大值為169.6 MPa，對(duì)應(yīng)的位移變形量的最大值為1.706 mm，門體結(jié)構(gòu)腔內(nèi)聯(lián)接區(qū)域最大變化位移發(fā)生在門體側(cè)邊中間位置，最大變化位移為0.002 mm；沖擊波波峰壓力為1 MPa下防爆門整體結(jié)構(gòu)所受應(yīng)力的最大值為339.3 MPa，對(duì)應(yīng)的位移變形量的最大值為3.461 mm，門體內(nèi)腔結(jié)構(gòu)聯(lián)接部件的最大變化位移發(fā)生在加強(qiáng)筋與門板的聯(lián)接處，最大變化位移為0.004 mm。

6.   結(jié)束語

利用顯式非線性動(dòng)力分析方法，通過仿真分析得出：防爆門的整體結(jié)構(gòu)與腔內(nèi)結(jié)構(gòu)所承受的沖擊壓力在材料的屈服強(qiáng)度之內(nèi)，符合材料的抗壓性能的要求；門體門板及加強(qiáng)筋零部件沒有出現(xiàn)塑性變形等情況，最大位移變形量不超過10 mm，沒有出現(xiàn)零件結(jié)構(gòu)失效、斷裂等破壞形式，符合防爆材料的剛度要求；防爆門中各零部件的相對(duì)位移變化小于1 mm，說明防爆門的密閉性能較好；防爆門無失效零部件，能夠保持完整，結(jié)構(gòu)強(qiáng)度和剛度達(dá)到安全要求。在煤礦井下實(shí)際應(yīng)用中防爆門抗爆性能良好，說明防爆門結(jié)構(gòu)模型和仿真數(shù)值分析有效可靠。

參考文獻(xiàn)

[1]羅星娜, 黃平, 錢新明. 礦用救生艙防爆門抗爆結(jié)構(gòu)安全性數(shù)值分析. 高壓物理學(xué)報(bào), 2013,27(6):889doi: 10.11858/gywlxb.2013.06.015

[2]宋維賓, 孫玉寧. 煤礦立風(fēng)井防爆門發(fā)展現(xiàn)狀分析. 中國安全生產(chǎn)科學(xué)技術(shù), 2015(6):108

[3]Katherine A, Margolis, Catherine Y, et al. Underground mine refuge chamber expectations training: program development and evaluation. Saf Sci,2011,49:522

[4]Houf W G, Evans G H, Schefer R W. A study of barrier walls for mitigation of unintended releases of hydrogen. Int J Hydrogen Energy,2011,36(3):2520

[5]Liu J Y, Zhang Q X, Chen Y H. Research of affection by different stiffener to lifesaving cabin property based on computer simulation. J Residuals Sci Tech,2016,13(8):701

[6]宋勝偉, 劉富利, 楊志, 等. 礦井自動(dòng)隔離門強(qiáng)度的數(shù)值模擬分析. 黑龍江科技大學(xué)學(xué)報(bào), 2016,26(1):58doi: 10.3969/j.issn.2095-7262.2016.01.014

[7]吳廣明, 許亞東, 李正國, 等. 艙室復(fù)合材料結(jié)構(gòu)的抗爆性能數(shù)值仿真. 中國艦船研究, 2019,14(4):14

[8]王文娟, 倪永帥, 李丹, 等. 基于有限元方法的煤礦避難硐室防護(hù)密閉門結(jié)構(gòu)動(dòng)力分析. 起重運(yùn)輸機(jī)械, 2017(6):45doi: 10.3969/j.issn.1001-0785.2017.06.013

[9]宋維賓. 煤礦立風(fēng)井防爆門安全防護(hù)理論與技術(shù)研究[學(xué)位論文]. 焦作: 河南理工大學(xué), 2018

[10]陳曉坤, 李海濤, 王秋紅, 等. 瓦斯爆炸載荷作用下礦用柱殼救生艙抗爆性分析及結(jié)構(gòu)優(yōu)化. 爆炸與沖擊, 2018,38(1):124doi: 10.11883/bzycj-2016-0248

[11]解江, 李翰, 周書婷, 等. 爆炸沖擊載荷下航空鋁合金平板動(dòng)態(tài)響應(yīng)數(shù)值分析方法. 應(yīng)用數(shù)學(xué)和力學(xué), 2017,38(4):410

[12]高軒能, 付詩琪. 空間高度對(duì)球面網(wǎng)殼結(jié)構(gòu)內(nèi)爆炸響應(yīng)的影響. 土木建筑與環(huán)境工程, 2017,39(4):107

[13]賈娟娟. 爆炸沖擊荷載作用下大型LNG儲(chǔ)罐動(dòng)力響應(yīng)分析[學(xué)位論文]. 青島: 中國石油大學(xué), 2016

文章來源——金屬世界