沈佳興， 徐 平， 于英華

(遼寧工程技術(shù)大學(xué) 機(jī)械工程學(xué)院，遼寧 阜新 123000)

煤炭生產(chǎn)行業(yè)具有高危險(xiǎn)性，受到許多不安全因素影響如：瓦斯爆炸、煤塵爆炸、冒頂、透水等。據(jù)不完全統(tǒng)計(jì)，我國每生產(chǎn)100萬噸煤炭約有3.1個(gè)生產(chǎn)人員付出生命[1-2]，可見煤炭生產(chǎn)的安全性是一個(gè)亟需改善的問題。救生艙的概念從2006年首次被提出后就得到廣大國內(nèi)外專家關(guān)注，相關(guān)研究設(shè)計(jì)人員也做了大量研究。但是相關(guān)研究主要集中在采用傳統(tǒng)方法設(shè)計(jì)、優(yōu)化救生艙的結(jié)構(gòu)和參數(shù)，其性能幾乎已經(jīng)達(dá)到極限[3-7]。為使救生艙的性能有根本性的提高，本文提出一種泡沫鋁填充結(jié)構(gòu)救生艙(下文簡稱泡沫鋁救生艙)，充分利用泡沫鋁輕質(zhì)、高比強(qiáng)度、高比剛度、隔熱保溫和抗沖擊吸能的特點(diǎn)以期使救生艙的性能得到根本性的改善和提高。

泡沫鋁是近些年來迅速發(fā)展起來的一種新型結(jié)構(gòu)功能于一體的多孔材料。泡沫鋁的壓縮應(yīng)力—應(yīng)變曲線有一段長長的“平臺”，因此泡沫鋁可以在幾乎恒定的壓力作用下吸收大量的能量，為其成為高效的吸能材料提供了可能；閉孔泡沫鋁由于結(jié)構(gòu)中有密閉的胞體且胞體內(nèi)部有密封的氣體，可以起到很好的隔熱作用，所以閉孔泡沫鋁的導(dǎo)熱率比密實(shí)基體金屬要低8倍～30倍，其可以作為隔熱材料應(yīng)用；泡沫鋁材料的阻尼是密實(shí)基體金屬的5倍～10倍，可以作為減振材料應(yīng)用[8-11]。為此將泡沫鋁作為填充材料應(yīng)用于救生艙的設(shè)計(jì)上，渴望為進(jìn)一步提高救生艙的安全、隔熱、輕質(zhì)和舒適等性能提供新思路和新途徑。

1 救生艙的結(jié)構(gòu)設(shè)計(jì)

選用某型號救生艙為研究原型，并以其中一段艙體為研究對象，其結(jié)構(gòu)如圖1所示，主要結(jié)構(gòu)包括：鋼板蒙皮、法蘭、縱向加強(qiáng)筋及橫向加強(qiáng)筋，其內(nèi)部尺寸為1 550 mm×1 740 mm×1 200 mm，鋼板蒙皮厚為8 mm，加強(qiáng)筋為10#槽鋼。泡沫鋁救生艙的結(jié)構(gòu)如圖2所示，主要結(jié)構(gòu)包括：法蘭、外板、內(nèi)板、泡沫鋁填充體和縱向加強(qiáng)筋及橫向加強(qiáng)筋，泡沫鋁救生艙的內(nèi)部尺寸與原型相同。外板厚為P1、內(nèi)板厚為P2、泡沫鋁厚為P4、縱、橫向加強(qiáng)筋寬分別為P3、P5，具體尺寸如表1所示。

圖1 原型模型Fig.1 Prototype model

圖2 泡沫鋁模型Fig.2 aluminum foam model

代號P1P2P3P4P5數(shù)值/mm3.555505.5

2 救生艙爆炸瞬態(tài)熱分析

根據(jù)分析可知救生艙的結(jié)構(gòu)和受到的載荷以及約束條件均具有對稱性，因此為節(jié)約分析所需內(nèi)存，提高計(jì)算效率。分析時(shí)只分析模型的一半并在剖切的對稱面上施加結(jié)構(gòu)對稱約束來保證分析結(jié)果的準(zhǔn)確性。

2.1 原型救生艙瞬態(tài)熱性能分析

原型救生艙的材料為Q345R鋼，密度7 860 kg/m3，彈性模量206 GPa，泊松比0.3，屈服強(qiáng)度345 MPa，強(qiáng)度極限510 MPa，線膨脹系數(shù)為1.2×10-5/K，比熱容為460 J/(kg·K)，導(dǎo)熱率為50 W/(m·K)[7,12]。采用六面體劃分網(wǎng)格，網(wǎng)格模型如圖3所示，包含有41 371個(gè)節(jié)點(diǎn)和14 376個(gè)單元。井下瓦斯爆炸時(shí)所形成的高溫是與瓦斯含量及與爆炸點(diǎn)的距離等因素不同而不同的，因此根據(jù)理論計(jì)算救生艙受的高溫較為困難，但根據(jù)國家相關(guān)規(guī)定可知救生艙應(yīng)該能夠承受1 200℃下至少3 s的外界環(huán)境瞬間高溫的能力[13-14]。

圖3 網(wǎng)格模型Fig.3 Mesh model

圖4 原型溫度云圖Fig.4 Prototype temperature

為了提高救生艙的安全性，分析在該高溫環(huán)境下3×2 s(2為安全系數(shù)[7])的熱性能。在救生艙的側(cè)板和頂板及法蘭的外側(cè)施加上述瞬態(tài)溫度載荷，在艙體底面外側(cè)施加25 ℃的地面溫度載荷，在救生艙的內(nèi)表面施加環(huán)境溫度為22 ℃的Stagnant Air-Simplified Case邊界條件。原型救生艙的溫度云圖和溫度曲線如圖4,5所示。

圖5 原型溫度曲線Fig.5 Prototype temperature curve

由圖4可知原型救生艙的最高溫度為138.02 ℃，出現(xiàn)在法蘭底面邊緣處，如圖中MAX標(biāo)記所示，救生艙的蒙皮外側(cè)溫度為90 ℃左右，蒙皮外側(cè)的加強(qiáng)筋溫度為115 ℃左右，蒙皮內(nèi)側(cè)溫度為70 ℃左右。由圖5可知原型救生艙的最高溫度和蒙皮內(nèi)側(cè)溫升略有不同，內(nèi)側(cè)溫度表現(xiàn)為0～1.5 s內(nèi)變化較緩慢，隨后逐步升高直至74.026 ℃，而最高溫度變化趨勢為0～1.5 s內(nèi)變化較快，隨后上升趨勢略有減緩。

2.2 泡沫鋁生艙瞬態(tài)熱性能分析

泡沫鋁救生艙的結(jié)構(gòu)如圖2所述，包括有外板，泡沫鋁填充體及內(nèi)板等結(jié)構(gòu)。文中選用閉孔泡沫鋁，其材料參數(shù)為：密度540 kg/m3，彈性模量12 GPa，泊松比0.33，抗壓強(qiáng)度14 MPa，線膨脹系數(shù)2×10-5/K，比熱容915 J/(kg·K)，導(dǎo)熱率15 W/(m·K)[8]；而模型內(nèi)，外板和加強(qiáng)筋為Q345R鋼，其材料參數(shù)見2.1節(jié)。分析時(shí)模型均采用六面體網(wǎng)格劃分，整體模型包含有82 553個(gè)單元，20 404個(gè)節(jié)點(diǎn)。泡沫鋁救生艙的熱載荷及熱邊界條件與原型救生艙相同。對泡沫鋁救生艙的熱性能求解計(jì)算后其溫度云圖和溫度曲線如圖6，7所示。

圖6 泡沫鋁溫度云圖Fig.6 Aluminum foam temperature

圖7 泡沫鋁溫度曲線Fig.7 Aluminum foam temperature curve

根據(jù)圖6，7可知泡沫鋁救生艙的最高溫度為141.44 ℃，出現(xiàn)在外板與法蘭連接附近，如圖中MAX標(biāo)記所示。泡沫鋁救生艙的最高溫度比原型救生艙的略有提升，這是因?yàn)榕菽X救生艙的外板厚度與原型蒙皮相比要薄一些，所以在相同的瞬態(tài)熱能量沖擊的條件下其溫度要高。泡沫鋁救生艙的外板外側(cè)溫度為120 ℃左右，內(nèi)板溫度為30 ℃左右，且分布較均勻。比較圖5和圖7可知，兩種救生艙的最高溫度變化趨勢幾乎完全一致，但泡沫鋁救生艙的內(nèi)板內(nèi)側(cè)溫度變化緩慢，最高溫度為35.622 ℃。泡沫鋁救生艙內(nèi)板溫度與原型救生艙的蒙皮內(nèi)板溫度相比減小了51.88%，證明泡沫鋁救生艙具有優(yōu)異的隔熱性能。

3 救生艙爆炸熱-壓力耦合沖擊分析

3.1 救生艙爆炸熱應(yīng)力分析

救生艙受熱后材料會發(fā)生膨脹而產(chǎn)生熱應(yīng)力和相應(yīng)的變形，所以需要分析救生艙在瞬態(tài)熱沖擊載荷下的結(jié)構(gòu)力學(xué)特性。

3.1.1 原型救生艙熱應(yīng)力分析

利用上文分析的熱性能結(jié)果在Workbench軟件中添加結(jié)構(gòu)力學(xué)分析單元，并在艙體底面添加固定約束，為模擬兩段艙體間的連接，在法蘭連接面添加對稱約束。求解后救生艙的熱應(yīng)力和熱變形如圖8，9所示。根據(jù)圖8可知，當(dāng)時(shí)間為3 ms時(shí)救生艙有最大熱應(yīng)力317.9 MPa，出現(xiàn)在法蘭和蒙皮連接處，同時(shí)在加強(qiáng)筋的周圍也出現(xiàn)了比較大的熱應(yīng)力。這是因?yàn)槟Ｐ徒Y(jié)構(gòu)受熱后會膨脹，但是由于蒙皮和加強(qiáng)筋及法蘭的厚度不同(法蘭厚度最大，加強(qiáng)筋厚度居中，蒙皮厚度最小)且法蘭的連接面添有對稱約束。結(jié)構(gòu)受到熱沖擊時(shí)，蒙皮的膨脹趨勢最大，其余兩結(jié)構(gòu)相對較小，蒙皮膨脹受到加強(qiáng)筋和法蘭的擠壓，所以在蒙皮與加強(qiáng)筋和法蘭的聯(lián)接處產(chǎn)生了較大應(yīng)力。根據(jù)圖9可知當(dāng)時(shí)間為1 ms時(shí)救生艙有最大熱變形0.835 16 mm，出現(xiàn)在側(cè)板和頂板聯(lián)接的過渡圓弧處，且變形方向?yàn)橄蛲馔钩?。這是因?yàn)轫敯搴蛡?cè)壁受熱后均有沿著其所在平面方向的膨脹趨勢，但側(cè)板底部受到底面的約束所以只能沿側(cè)板上測膨脹，并在圓弧過渡處兩種膨脹趨勢發(fā)生相互擠壓產(chǎn)生最大變形。

圖8 原型熱應(yīng)力Fig.8 Prototype thermal stress

圖9 原型熱變形Fig.9 Prototype thermal deformation

3.1.2 泡沫鋁救生艙熱應(yīng)力分析

同理根據(jù)泡沫鋁熱性能分析結(jié)果在Workbench軟件中添加結(jié)構(gòu)力學(xué)分析單元，并在泡沫鋁救生艙的相應(yīng)位置添加與原型相同的結(jié)構(gòu)約束。求解結(jié)果如圖10,11所示。根據(jù)圖10可知泡沫鋁救生艙在3 ms時(shí)有最大熱應(yīng)力247.4 MPa，且主要分布在救生艙的側(cè)板下部和頂板中部。根據(jù)比較可知泡沫鋁救生艙最大應(yīng)力比原型降低22.18%。根據(jù)圖11可知泡沫鋁救生艙在1 ms時(shí)有最大熱變形為0.709 97 mm，出現(xiàn)位置與原型相同，且變形方向?yàn)橄蛲馔钩?。，形成熱變形的原因與原型救生艙的原因相同。通過比較可知泡沫鋁救生艙最大變形比原型降低14.99%。

圖10 泡沫鋁熱應(yīng)力Fig.10 Aluminum foam thermal stress

圖11 泡沫鋁熱變形Fig.11 Aluminum foam thermal deformation

3.2 救生艙爆炸熱-壓力耦合沖擊分析

煤礦井下瓦斯爆炸瞬間會產(chǎn)生大量的高溫、高壓氣體迅速向周圍膨脹形成高溫，高壓的超壓沖擊波，因此非常有必要研究救生艙在熱沖擊和壓力沖擊耦合情況下的力學(xué)特性。

瓦斯爆炸形成的壓力沖擊波主要包括壓力激增階段t1、正壓階段t2、負(fù)壓階段t3，如圖12所示，壓力沖擊破壞主要由正壓階段引起。由于壓力沖擊波在傳播過程中遇到結(jié)構(gòu)表面時(shí)會形成反射，反射波超壓大小又與入射波超壓峰值、入射角、結(jié)構(gòu)剛度、面積大小等因素密切相關(guān)，因此計(jì)算結(jié)構(gòu)實(shí)際所受壓力較為困難。根據(jù)國家相關(guān)規(guī)定[13-14]可知井下瓦斯爆炸沖擊波曲線可簡化為線性下降的三角形沖擊波壓力衰減曲線如圖12中點(diǎn)劃線所示，t4為等效超壓作用時(shí)間。根據(jù)規(guī)定要求可知救生艙的抗沖擊壓力不低于0.3 MPa，考慮到實(shí)際工況的復(fù)雜性分析時(shí)取超壓峰值ΔP1為0.3 MPa×2(2為安全系數(shù))；t1取2 ms；t4應(yīng)不小于300 ms，取300 ms[13-14]。

圖12 爆炸沖擊波衰減曲線Fig.12 Explosion shock wave attenuation curve

Q345為應(yīng)變率敏感材料，所以救生艙在遭受爆炸形成的壓力沖擊載荷時(shí)必須考慮材料應(yīng)變率效應(yīng)的影響，材料的應(yīng)變率效應(yīng)通過Cowper-Symond模型來考慮，即：

(1)

泡沫鋁的動態(tài)性能是隨著泡沫鋁基體材料、相對密度、孔穴結(jié)構(gòu)(開孔，閉孔)、孔穴尺寸及孔穴形狀等表征的不同，而表現(xiàn)出不同的特征。根據(jù)文獻(xiàn)[12]可知泡沫鋁的抗沖擊特性受應(yīng)變率的影響目前還存在爭議。但文獻(xiàn)[12]中DESHPANDE研究的泡沫鋁種類、性能及其研究的應(yīng)變率變化范圍與本文中泡沫鋁的載荷情況較為吻合，所以作者采納DESHPANDE教授的研究結(jié)論，即認(rèn)為本文中的泡沫鋁材料是應(yīng)變率無關(guān)的。

3.2.1 原型救生艙熱及壓力耦合沖擊分析

結(jié)合3.1.1節(jié)中分析內(nèi)容，救生艙的約束和熱載荷不變，在救生艙的側(cè)壁和頂板及法蘭的外表面施加上述線性下降的瞬態(tài)三角形壓力沖擊載荷，設(shè)置壓力峰值為0.6 MPa，且t1=2 ms，t4=300 ms。設(shè)置總計(jì)算時(shí)間為0.302 s，迭代步長為0.5 ms。計(jì)算結(jié)果如圖13,14所示。

根據(jù)圖13可知原型救生艙在熱和壓力沖擊的共同作用下在9.5 ms時(shí)在救生艙側(cè)板的第三、第四條水平加強(qiáng)筋處有最大應(yīng)力385.31 MPa，該最大應(yīng)力已經(jīng)超過材料的屈服強(qiáng)度，發(fā)生塑性破壞；根據(jù)圖14可知在3 ms時(shí)在救生艙的側(cè)板中部有最大變形3.706 1 mm，且變形方向?yàn)橄騼?nèi)凹陷。模型形成與單純熱沖擊相反的變形的原因是因?yàn)槟Ｐ褪艿酱怪迸擉w表面向內(nèi)的爆炸壓力沖擊，模型有向內(nèi)側(cè)變形的趨勢，且在3ms時(shí)原型艙體的最高溫度約為25.24 ℃其熱變形趨勢較小，而該時(shí)刻的壓力沖擊導(dǎo)致的變形趨勢要遠(yuǎn)大于熱變形，兩者耦合作用下模型產(chǎn)生內(nèi)凹的變形。

圖13 原型耦合應(yīng)力Fig.13 Prototype coupling stress

圖14 原型耦合變形Fig.14 Prototype couplingdeformation

3.2.2 泡沫鋁救生艙熱及壓力耦合沖擊分析

同理在3.1.2節(jié)分析的基礎(chǔ)上在泡沫鋁救生艙相應(yīng)位置施加與原型分析時(shí)相同的壓力沖擊載荷，并設(shè)置求解總時(shí)間為0.302 s，迭代步長為0.5 ms，其余設(shè)置不變。結(jié)果如圖15，16所示。

圖15 泡沫鋁耦合應(yīng)力Fig.15 Aluminum foam coupling stress

圖16 泡沫鋁耦合變形Fig.16 Aluminum foam coupling deformation

根據(jù)圖15可知泡沫鋁救生艙在6 ms時(shí)在側(cè)板的中下部的鋼板有最大應(yīng)力268.97 MPa，內(nèi)部的泡沫鋁應(yīng)力為10.44 MPa，均未超過材料的屈服強(qiáng)度，通過與原型比較可知最大應(yīng)力減小30.19%，證明泡沫鋁救生艙的強(qiáng)度得到大幅度提升。根據(jù)圖16可知泡沫鋁救生艙在3 ms時(shí)有最大變形2.938 1 mm，出現(xiàn)位置與方向與原型救生艙耦合沖擊下的分析結(jié)果相同，且其變形的原理與原型救生艙的相同。但比較兩種結(jié)構(gòu)救生艙的耦合沖擊下的變形可知泡沫鋁救生艙的最大變形比原型減小20.72%，證明泡沫鋁救生艙的剛度得到大幅度提升。通過強(qiáng)度和剛度對比證明泡沫鋁救生艙的安全性有了顯著改善。

分別提取原型救生艙和泡沫鋁救生艙在熱和壓力耦合沖擊作用下的變形和應(yīng)力時(shí)間歷程曲線如圖17，18所示。由圖17可知兩救生艙在相同的載荷下泡沫鋁救生艙的變形波動時(shí)間明顯縮短。圖18中點(diǎn)A、B為壓力和熱沖擊影響大小的分界點(diǎn)，點(diǎn)A、B之前壓力沖擊起到主導(dǎo)作用，之后熱沖擊起到主導(dǎo)作用。

圖17 耦合變形Fig.17 Coupling deformation

圖18 耦合應(yīng)力Fig.18 Coupling stress

4 結(jié) 論

分別對兩種救生艙作熱-壓力耦合沖擊仿真分析，并將結(jié)果對比分析得到如下結(jié)論：

(1) 在相同的熱沖擊載荷下，泡沫鋁內(nèi)板最高溫度比原型救生艙蒙皮內(nèi)側(cè)最高溫度降低51.88%。證明泡沫鋁救生艙具有優(yōu)異的隔熱性能。

(2) 在相同的熱沖擊載荷下，泡沫鋁救生艙的熱應(yīng)力和熱變形比原型分別減小22.18%和14.99%。

(3) 在相同的熱及壓力耦合沖擊載荷下，泡沫鋁救生艙的應(yīng)力和變形分別降低30.19%和20.72%，證明泡沫鋁救生艙的安全性得到顯著提升。