馬尚禎
(河南能源化工集團鶴煤公司,河南鶴壁 458030)
煤礦風(fēng)井防爆門是隨著煤礦通風(fēng)設(shè)備的出現(xiàn)而出現(xiàn)的,當井下發(fā)生瓦斯爆炸時,防爆門裝設(shè)在井口,能夠快速開啟,泄放爆炸沖擊波,從而對主通風(fēng)機免遭破壞起到保護作用。煤礦風(fēng)井防爆門是一種可以防止瓦斯、煤塵爆炸時毀壞主要通風(fēng)機的安全設(shè)備,且在主要通風(fēng)機停運時打開,起到防止井下硐室及主要回風(fēng)道瓦斯積聚的作用[1-2]。然而,當發(fā)生瓦斯爆炸時,防爆門在強大的爆炸壓力沖擊波作用下,容易發(fā)生變形甚至損壞而無法關(guān)閉,致使反風(fēng)時風(fēng)流短路,井下有毒氣體排出困難,影響井下救援工作[3]。因此,開展煤礦風(fēng)井防爆門的研究是煤礦安全生產(chǎn)、提升煤礦安全生產(chǎn)保障水平的迫切需求。
國內(nèi)外對防爆門的研究主要集中在爆破波的發(fā)展與形成上,對其在空氣中的傳播和衰減研究內(nèi)容相對較少,已有學(xué)者針對防爆門的密封結(jié)構(gòu)、冬季易結(jié)冰、備用防爆門快速復(fù)位及鎖扣等問題進行了相關(guān)研究[4-6],但考慮實際爆炸沖擊波作用下的防爆門動態(tài)響應(yīng)過程進行分析仍然較少,與此同時,通過數(shù)值模擬分析防爆門強度及結(jié)構(gòu)優(yōu)化也不常見。因此,本文在針對某煤礦風(fēng)井進行防爆門結(jié)構(gòu)設(shè)計的基礎(chǔ)上,利用FLUENT 軟件對瓦斯爆炸壓力沖擊波作用下防爆門動態(tài)響應(yīng)過程進行數(shù)值模擬,得出瓦斯爆炸沖擊波在不同時刻作用于防爆門上的壓力分布及其運動規(guī)律,同時運用ANSYS軟件對防爆門進行動態(tài)受力分析。在此基礎(chǔ)上,對防爆門現(xiàn)有結(jié)構(gòu)進行了優(yōu)化和強度校核。
目前,煤礦風(fēng)井防爆門結(jié)構(gòu)形式主要有傘式、花瓣式、蝴蝶式等。在礦井災(zāi)變時期,瓦斯爆炸或瓦斯與煤塵復(fù)合爆炸當量不可估計,因此防爆門抗沖擊強度無法準確確定,容易導(dǎo)致上述幾種結(jié)構(gòu)受到爆炸壓力作用時沖擊波氣流在防爆門中心底部積聚,致使其發(fā)生變形破壞[7]。本文設(shè)計的風(fēng)井防爆門是一種雙開門式結(jié)構(gòu),兩扇門可以繞各自的轉(zhuǎn)動軸自由轉(zhuǎn)動,轉(zhuǎn)動角度范圍為0~90°,其作用是在風(fēng)井需要密封時,將井口罩蓋,當瓦斯氣流發(fā)生二次爆炸或多次爆炸時,此防爆門能自動打開,并且在爆炸過后可自動復(fù)位。
雙開式防爆門總體結(jié)構(gòu)如圖1所示。風(fēng)井井筒直徑為4.8 m,方形密封池尺寸為7.5 m ×7.5 m ×4.1 m;單扇防爆門尺寸為2.735 m × 5.2m,厚度為3 mm,門框為50 角鋼,單扇防爆門質(zhì)量為480 kg,配重質(zhì)量為775 kg。為了井筒示意圖和模擬結(jié)果簡潔,配重在圖1 中未畫出,在防爆門運動計算過程中加以考慮。
圖1 防爆門總體結(jié)構(gòu)設(shè)計
為了分析防爆門在瓦斯爆炸沖擊波作用下的安全可靠性,采用FLUENT軟件對瓦斯爆炸壓力場進行數(shù)值模擬,得出防爆門上的壓力分布及防爆門運動情況[8]。
根據(jù)防爆門結(jié)構(gòu)建立1∶1 物理模型,劃分計算網(wǎng)格,如圖2 所示。為提高計算精度,靠近壁面處的計算網(wǎng)格進行局部加密,圖2 中ac和bc為雙開式防爆門,分別設(shè)為可繞點a和點b轉(zhuǎn)動的剛體,其余均為固定壁面。防爆門轉(zhuǎn)動角加速度根據(jù)牛頓第二定律確定,即與所受壓力和重力(包括活動門和配重)的力矩代數(shù)和成正比。井筒和方形密封池充滿化學(xué)當量比的CH4和空氣預(yù)混氣體。初始時刻,在井筒底部利用局部高溫(2400 ℃)方法點火引爆。流動為非穩(wěn)態(tài)湍流流動,湍流模型采用k-ε 湍流方程,壁面采用標準壁面函數(shù),燃燒反應(yīng)采用適用于湍流燃燒的EBU 渦擴散模型,并采用收斂較好的PISO算法進行迭代求解,時間步長設(shè)為1 ×10-5s。
圖2 煤礦風(fēng)井防爆門物理模型及計算網(wǎng)格
防爆門動態(tài)開啟過程如圖3 所示。圖3(a)~(b)分別為爆炸壓力沖擊波作用下防爆門開啟至25°和90°時的壓力場。可以看出,防爆門開啟初期,由于防爆門開啟角度小、泄爆口較小,使得氣流積聚壓力持續(xù)升高,防爆門所受的最高壓力由26 kPa上升至782 kPa,活動門上平均壓力由23 kPa 上升至467 kPa。當防爆門開啟角度增大時,所積聚壓力得到足夠釋放,壓力開始逐漸下降。
圖3 防爆門動態(tài)開啟過程
圖4~5 分別顯示了開啟角度為25°和90°時防爆門上的壓力分布。當開啟角度為25°時,防爆門上壓力分布不均勻,靠近轉(zhuǎn)動軸位置的壓力最高,泄爆口附近壓力最低,說明防爆門泄爆作用十分顯著,但由于泄爆口較小,使得防爆門上的最高壓力仍然較高。當開啟角度為90°時,防爆門上的壓力較為均勻,且已下降至220 kPa左右。
圖4 防爆門開啟25°時的壓力分布
圖5 防爆門開啟90°時的壓力分布
采用ANSYS(9.0 版本)結(jié)構(gòu)分析軟件,對單扇活動門在動態(tài)非線性壓力作用下的應(yīng)力及變形進行分析計算,并通過校核材料強度,判斷變形類型屬于彈性變形還是塑性變形,從而檢驗活動門的安全可靠性。
單扇防爆門尺寸為2.735 m × 5.2 m,門框為50 角鋼,門框中間部分為米字型布置,門框與門板焊接連接,門板厚度為3 mm。圖6 所示為防爆門局部結(jié)構(gòu)。在應(yīng)力分析過程中,門框選用Beam188 梁結(jié)構(gòu),防爆門及門框50 角鋼局部如圖6所示。門板選用Shell63 薄板結(jié)構(gòu)。定義防爆門彈性模量為2.06× 1011Pa,泊松比為0.3,密度為7800 kg/ m3。施加運動自由度約束條件:左門框(軸承附件)進行自由度約束,Ux=Uy=Uz=Rx=Rz=0,僅保留y方向的轉(zhuǎn)動自由度。施加載荷包括壓力載荷和重力載荷。其中,壓力是隨時間變化的非線性動態(tài)載荷,如表1 所示。在應(yīng)力計算中,將流場計算中的24 ms 時刻作為初始時刻,建立活動門上所受平均壓力與時間的變化關(guān)系。相比于最高壓力,平均壓力更能反映對活動門整體的沖擊作用。選擇大位移瞬態(tài)求解器進行迭代求解。求解控制時間與防爆門總開啟時間對應(yīng),總時間為0.03 s,時間步長為0.001 s,每一時間步的載荷與壓力分析結(jié)果對應(yīng)。
表1 不同時刻防爆門所受壓力
圖6 防爆門局部結(jié)構(gòu)
圖7 所示為防爆門的應(yīng)力云圖??梢钥闯?,其最大應(yīng)力為230 MPa,接近鋼的許用應(yīng)力245 MPa;最大變形量為38.9 mm,主要為彈性變形。由于防爆門是受到?jīng)_擊波強大的扭力作用,因此應(yīng)力集中在防爆門的左側(cè)主要表現(xiàn)為彎曲應(yīng)力。
圖7 防爆門應(yīng)力分布
作為防爆門的保護對象,煤礦主通風(fēng)機的抗爆性能應(yīng)作為防爆門優(yōu)化設(shè)計的重要依據(jù)。防爆門的防爆性能過高屬于過度保護,而保護性能過低,則起不到保護作用,因此,對防爆門的結(jié)構(gòu)進行了優(yōu)化設(shè)計[9-10]。單扇門在受力運動時除了滿足強度要求之外,還需要滿足剛度要求,這主要是因為,如果門的剛度不夠,會造成門受力發(fā)生彎曲變形,使得門在再次關(guān)閉時出現(xiàn)問題,所以必須提高防爆門的抗彎強度。從上面受力分析可知,防爆門在壓力波沖擊作用下所受最大應(yīng)力主要集中在左側(cè),即靠近轉(zhuǎn)動軸側(cè)附近,因此可在防爆門左側(cè)加裝兩根加強筋,其幾何結(jié)構(gòu)和優(yōu)化后應(yīng)力分布結(jié)果如圖8 所示。可以看出,結(jié)構(gòu)優(yōu)化后的最大應(yīng)力下降為189 MPa,低于鋼的許用應(yīng)力245 MPa,最大變形量為32 mm,未發(fā)生塑性變形,說明防爆門結(jié)構(gòu)優(yōu)化方案合理,防爆門整體可靠性符合要求,優(yōu)化后效果較好,大大提高了煤礦防爆門的安全性能。
圖8 結(jié)構(gòu)優(yōu)化后防爆門應(yīng)力分布
本文針對煤礦風(fēng)井防爆門在煤礦發(fā)生爆炸時,在強大的爆炸壓力沖擊波作用下,易發(fā)生變形甚至損壞的問題,采用了調(diào)查分析、數(shù)值模擬、理論分析與計算相結(jié)合的方法,通過數(shù)值模擬分析防爆門的強度及結(jié)構(gòu),在此基礎(chǔ)上,對防爆門現(xiàn)有結(jié)構(gòu)進行了優(yōu)化和強度校核[11]。
(1)在防爆門結(jié)構(gòu)設(shè)計基礎(chǔ)上,利用FLUENT 軟件對瓦斯爆炸流場進行數(shù)值模擬,獲得了不同時刻防爆門上的壓力分布及其動態(tài)響應(yīng)規(guī)律。
(2)將不同時刻壓力載荷導(dǎo)入ANSYS 受力分析軟件,進而對防爆門進行動態(tài)受力分析,得到了防爆門應(yīng)力分布,在此基礎(chǔ)上進行了結(jié)構(gòu)優(yōu)化,材料強度校核結(jié)果顯示優(yōu)化效果較好。