激波風(fēng)洞高低壓段鋼膜片破裂特性研究1)

聶少軍 汪運鵬,2) 薛曉鵬 姜宗林

*(中國科學(xué)院力學(xué)研究所高溫氣體動力學(xué)國家重點實驗室,北京 100190)

?(中國科學(xué)院大學(xué)工程科學(xué)學(xué)院,北京 100049)

**(中南大學(xué)航空航天學(xué)院,長沙 410083)

引言

激波管是實驗室中能夠產(chǎn)生可控激波并且利用激波壓縮試驗氣體的常用試驗裝置,用來模擬實際工作條件,它能瞬間將試驗氣體升溫、加速和增壓[1].激波風(fēng)洞由一個激波管以及連于其后的噴管、實驗段等風(fēng)洞主要部件組成,包含相連的驅(qū)動段和被驅(qū)動段,驅(qū)動段貯存高壓、高能量的工作氣體,被驅(qū)動段貯存低壓試驗氣體,試驗前兩者通常被膜片隔開.試驗時,上游的驅(qū)動段氣體經(jīng)過爆轟驅(qū)動[2]或者其他驅(qū)動方法[3-4]產(chǎn)生脈沖高壓,巨大的壓差會導(dǎo)致膜片破裂,繼而在下游產(chǎn)生所需的運動激波.在進行激波風(fēng)洞試驗時,膜片不僅需要承受驅(qū)動段的初始充氣壓力并且不發(fā)生破裂,而且還需要保證在點火驅(qū)動時完全打開而且不發(fā)生脫落.膜片能夠瞬間打開并且完全打開是形成激波的關(guān)鍵前提條件,同時也是激波風(fēng)洞試驗中普遍存在的難題.

國外一些激波管專家學(xué)者針對膜片破裂特性進行了相關(guān)研究.在20 世紀(jì)50 年代,一些學(xué)者通過激波管試驗對膜片破裂的規(guī)律進行研究,發(fā)現(xiàn)膜片破裂打開是一個非瞬態(tài)過程[5-6].隨后,有學(xué)者從理論上建立了簡單的膜片破裂的理論模型,通過實驗研究了有效破膜時間與膜片材料的關(guān)系,但是理論模型與實驗獲得的破膜時間存在較大差距[7].隨后,Hickman 等[8]通過實驗研究的方法獲得了膜片開啟過程中的高速陰影照片,觀測到了膜片開啟時存在的擬定常自由射流.范良藻和俞鴻儒[9]進一步研究了膜片材料對有效破膜時間的影響,結(jié)果表明有效破膜時間與激波管直徑和膜片材料密度成正相關(guān)關(guān)系,與膜片壓差成負(fù)相關(guān)關(guān)系.林建民[10]研究了激波風(fēng)洞中不同外形膜片對破膜壓力的影響,比較了金屬成型膜片和平板膜片的承載能力,發(fā)現(xiàn)金屬成型膜片更適合用于爆轟驅(qū)動的激波風(fēng)洞.栗繼偉[11]研究了不同膜片厚度對破膜過程造成的壓力損失的影響,在JFX 氫氧爆轟驅(qū)動激波風(fēng)洞中進行試驗,結(jié)果表明膜片總厚度和有效厚度對破膜壓力損失影響較大,當(dāng)厚度增大時,破膜壓力損失也會相應(yīng)增大.

隨著計算機技術(shù)的發(fā)展,學(xué)者已不再局限于實驗研究方法,針對激波管中膜片破裂過程的數(shù)值模擬方法成為可能.目前已有學(xué)者通過數(shù)值方法模擬了卡魯尼亞大學(xué)的SS304 激波管中的膜片破裂過程,分析了在不同厚度(1 mm,2 mm 或3 mm)的情況下鋁膜片的破裂過程,包括破膜壓力和膜片開啟時間隨膜片厚度的變化規(guī)律,得到的結(jié)論為:當(dāng)膜片的無量綱厚度(即膜片有效厚度與膜片總厚度的比值)分別為0.33 和0.67 時,有效厚度越小,膜片破裂所需的壓力越小,破裂所需的時間越長[12].

為了研究風(fēng)洞中膜片的破裂特性規(guī)律,一些學(xué)者主要采用傳統(tǒng)風(fēng)洞試驗的方法研究了鋁膜片的厚度對破膜壓力和有效破膜時間的影響,而對鋼膜片的研究比較少,因為傳統(tǒng)風(fēng)洞試驗中鋁膜片的承壓能力已能達(dá)到試驗要求.在目前的激波風(fēng)洞試驗中,鮮有學(xué)者采用數(shù)值模擬的方法研究脈沖型風(fēng)洞中膜片的破裂過程,尤其是膜片的凹槽長度對破裂過程的影響.

因此,本文將采用數(shù)值模擬的方法,利用有限元分析軟件中的顯式動力學(xué)模塊對膜片破裂過程進行瞬態(tài)結(jié)構(gòu)分析,并與復(fù)現(xiàn)高超聲速飛行條件激波風(fēng)洞(以下簡稱“JF-12 復(fù)現(xiàn)風(fēng)洞”)[13]中的試驗結(jié)果進行比對,分析爆轟驅(qū)動的激波風(fēng)洞中不同有效厚度、總厚度和凹槽長度對膜片破裂特性的影響.從工程應(yīng)用的角度來看,探尋激波風(fēng)洞中高低壓段膜片的破裂特性有助于預(yù)測膜片在風(fēng)洞運行中打開時的壓差與破膜時間,節(jié)省試驗成本,縮短試驗時間,提升試驗流場品質(zhì),促進我國激波風(fēng)洞試驗研究的有利發(fā)展.本工作采用有限元數(shù)值模擬的方法,推導(dǎo)了膜片破裂的動力學(xué)方程,研究了膜片破裂特性規(guī)律,以期為JF-12 復(fù)現(xiàn)風(fēng)洞試驗提供參考數(shù)據(jù).

1 膜片破裂的動態(tài)模型建立及研究方法

2012 年中國科學(xué)院力學(xué)研究所成功研制了JF-12 復(fù)現(xiàn)風(fēng)洞,該風(fēng)洞可復(fù)現(xiàn)25～40 km 高空、馬赫數(shù)5～9 的飛行條件的飛行器氣動試驗,有效試驗時間超過了100 ms[14-17].該風(fēng)洞由卸爆段、爆轟驅(qū)動段、激波管段、噴管、試驗段和真空段等部分組成.卸爆段和爆轟驅(qū)動段、爆轟驅(qū)動段和激波管段、激波管段和噴管段均有膜片隔開,爆轟驅(qū)動段和激波管段管道為等截面圓柱形,主膜片是指爆轟驅(qū)動段和激波管段之間的膜片[18-19].常見的爆轟驅(qū)動技術(shù)有正向爆轟和反向爆轟兩種運行模式[20],JF-12 復(fù)現(xiàn)風(fēng)洞采用反向爆轟驅(qū)動技術(shù)[21-22]產(chǎn)生高溫高壓氣體.本文在數(shù)值模擬的過程中,采用階躍載荷壓力模擬試驗中爆轟驅(qū)動產(chǎn)生的壓力.在爆轟驅(qū)動的激波風(fēng)洞試驗中,爆轟驅(qū)動段和激波管段管道的剛度遠(yuǎn)大于膜片的剛度,理論分析時可將管道視作固定端.研究方法是通過有限元軟件對膜片的破裂過程進行動態(tài)模擬,將模擬結(jié)果與JF-12 復(fù)現(xiàn)風(fēng)洞試驗結(jié)果進行對比分析.研究目標(biāo)是通過數(shù)值模擬方法對激波風(fēng)洞中膜片破裂的特性規(guī)律進行研究,為JF-12 復(fù)現(xiàn)風(fēng)洞試驗提供膜片參考數(shù)據(jù).

1.1 膜片幾何模型

激波風(fēng)洞中的主膜片如圖1 所示,膜片可以簡化為一個圓柱體,圓柱體的高即膜片的總厚度,在膜片的中心開有適當(dāng)深度的互相交叉垂直的“十”字形凹槽,凹槽的作用是使膜片順著凹槽順利打開,保證膜片打開的質(zhì)量.采用三維建模軟件繪制膜片模型,如圖2 所示.

圖1 激波風(fēng)洞主膜片示意圖Fig.1 Schematic diagram of main diaphragm of shock tunnel

圖2 主膜片模型Fig.2 Main diaphragm model

膜片有效厚度E是膜片總厚度T與膜片凹槽深度D之差,膜片的凹槽長度是基于爆轟驅(qū)動段和激波管段管道內(nèi)徑設(shè)定的.膜片模型尺寸信息見表1.

表1 膜片模型尺寸信息(單位:mm)Table 1 Main diaphragm model size information(Unit:mm)

在上述模型的基礎(chǔ)上,利用三維建模軟件建立了不同總厚度T、有效厚度E和凹槽長度L的膜片模型,總厚度有3.5 mm,2.5 mm 和1.5 mm 三種類型,有效厚度有0.5 mm,1 mm 和1.5 mm 三種類型,凹槽長度設(shè)定有400 mm,380 mm 和360 mm 三種類型,具體尺寸如表2 所示.

表2 不同外形的膜片模型尺寸信息(單位:mm)Table 2 Dimensional information of diaphragm models with different shapes(Unit:mm)

1.2 膜片破裂的動力學(xué)模型

針對激波風(fēng)洞試驗,在爆轟壓力的作用下,膜片破裂是一個典型的非線性動力學(xué)問題,膜片破裂時膜片周圍復(fù)雜的流動模式以及膜片內(nèi)部的動態(tài)應(yīng)力-應(yīng)變過程導(dǎo)致該問題十分復(fù)雜,如圖3 所示.為了簡化該問題,便于研究,假設(shè)膜片在沒有任何變形的情況下瞬間破裂;膜片的作用力均勻分布且是開口面積的線性函數(shù);在破裂過程中,膜片彎曲應(yīng)力產(chǎn)生的力矩為常數(shù)[8].

圖3 膜片打開過程示意圖Fig.3 Schematic diagram of diaphragm opening process

根據(jù)剛體定軸轉(zhuǎn)動微分方程,膜片打開角度θ滿足

式中,I是膜片繞軸的轉(zhuǎn)動慣量,M1和M2分別是高低壓力差和膜片彎曲應(yīng)力產(chǎn)生的力矩,表達(dá)式為

式中,P5和P1分別為作用在膜片上游和下游的壓力;A為膜片有效面積;為膜片內(nèi)徑的長度;P4為驅(qū)動段的初始壓力;S是無量綱強度參數(shù),表示膜片彎曲應(yīng)力力矩與膜片最大力矩的比值,根據(jù)假設(shè),在膜片破裂過程中,S為常數(shù).

為了便于分析膜片打開角度θ,令

將式(5)積分一次,得

將式(6)積分一次,即可得到

由式(7)可知,當(dāng)τ=0 時,θ=0°,dθ/dτ=0,膜片處于初始狀態(tài),未發(fā)生破裂;當(dāng)θ=90°時,膜片完全打開完成破裂.

1.3 研究方法

在完成模型建立之后,采用有限元分析軟件中的顯式動力學(xué)對膜片破裂過程進行數(shù)值模擬計算.顯式動力學(xué)是一種基于顯式算法對動力學(xué)方程進行求解的方法,具有較好的穩(wěn)定性,一般不存在收斂性問題[23].顯式動力學(xué)適合用于解決短時間(ms)內(nèi)發(fā)生的事件,如瞬態(tài)、大應(yīng)變、大變形、材料的破壞、材料的完全失效或者伴隨復(fù)雜接觸的結(jié)構(gòu)問題.采用顯式動力學(xué)模塊對膜片進行有限元分析的步驟如下圖4 所示.

圖4 膜片破裂過程數(shù)值計算流程圖Fig.4 Flow chart of numerical calculation for diaphragm rupture process

膜片破裂過程的數(shù)值模擬計算步驟包括計算前處理、計算和計算后處理操作,計算前處理操作主要包括建立膜片幾何模型、設(shè)置膜片材料性質(zhì)、劃分膜片網(wǎng)格以及建立膜片初始邊界條件;計算后處理是對計算結(jié)果進行相關(guān)操作,以輸出變形、應(yīng)變、應(yīng)力等結(jié)果,具體包括查看結(jié)果、輸出結(jié)果、誤差估計等內(nèi)容.

1.3.1 膜片材料性質(zhì)

膜片材料是影響膜片破裂的重要因素之一,選擇膜片材料時,應(yīng)保證破膜壓力的穩(wěn)定,即破膜過程中保證破膜壓力不變;同時,膜片應(yīng)當(dāng)具有一定的承壓能力,在初始高壓作用下,膜片不能被剪切破壞,而在階躍載荷作用下,膜片需要順利打開且不發(fā)生脫落.傳統(tǒng)風(fēng)洞試驗中,高壓區(qū)和低壓區(qū)的壓差不大,膜片材料可以選擇純鋁,但是在激波風(fēng)洞試驗中,爆轟壓力可以達(dá)到十幾兆帕甚至幾十兆帕,因此常采用鋼(STEEL 4340)膜片,其具體參數(shù)如表3 所示.

表3 STEEL 4340 材料屬性Table 3 STEEL 4340 material properties

1.3.2 膜片網(wǎng)格劃分

在有限元分析計算中,只對網(wǎng)格單元和節(jié)點進行計算,在計算之前需要對膜片模型劃分網(wǎng)格.增大網(wǎng)格密度可以適當(dāng)提升計算精度,但是不能彌補錯誤的假設(shè)和輸入帶來的錯誤,同時,網(wǎng)格密度過大會導(dǎo)致計算量和CPU 計算時間劇烈延長,而且隨著網(wǎng)格的加密,計算機浮點運算的舍入誤差也會變大.理想的網(wǎng)格密度應(yīng)該首先滿足計算精度的要求,其次,當(dāng)網(wǎng)格細(xì)化到一定程度后,計算結(jié)果的變化可以忽略不計,此時便得到理想的網(wǎng)格密度[24].在經(jīng)過多次嘗試之后,確定了膜片深度為2.5 mm 時網(wǎng)格劃分的一些參數(shù)設(shè)置,如圖5 所示.

圖5 膜片網(wǎng)格Fig.5 Diaphragm mesh

網(wǎng)格總體單元尺寸為3 mm,凹槽面處網(wǎng)格單元尺寸為2.5 mm,網(wǎng)格節(jié)點數(shù)為59 885,網(wǎng)格單元數(shù)為206 433.

1.3.3 初始邊界條件設(shè)定

在完成網(wǎng)格劃分之后,接下來是設(shè)定膜片破裂過程的初始邊界條件,具體包括計算時間設(shè)定、膜片的載荷設(shè)置和固定約束設(shè)置.

有限元計算時間的設(shè)定需要考慮兩個方面的因素:一是必須保證膜片能夠在計算時間內(nèi)完全打開;二是計算時間必須滿足實際情況,不應(yīng)超過激波風(fēng)洞試驗的總時長,綜合考慮將數(shù)值計算時間設(shè)置為1000 μs.

考慮實際情況,膜片受力簡化為壓力載荷作用,如圖6 和圖7 所示.在初始充氣壓力作用下,即0～200 μs 內(nèi),膜片正面受力為2 MPa,背面受力為38 kPa,方向垂直于膜片平面;在t=200 μs 時,產(chǎn)生爆轟壓力,膜片正面壓力驟增至32 MPa,背面壓力不變,即在200～1000 μs 內(nèi),膜片正面壓力為32 MPa,背面壓力為38 kPa,壓力載荷具體計算過程將在第2章中給出.

圖6 高壓區(qū)壓力Fig.6 High pressure zone pressure

圖7 低壓區(qū)壓力Fig.7 Low pressure zone pressure

激波風(fēng)洞管道的剛度遠(yuǎn)大于膜片的剛度,在試驗時膜片被夾膜機夾在管道中,故將膜片外環(huán)的圓環(huán)面和側(cè)面設(shè)置為固定約束,如圖8 所示.

圖8 固定約束設(shè)置Fig.8 Fixed constraint settings

2 計算結(jié)果與驗證分析

2.1 JF-12 復(fù)現(xiàn)風(fēng)洞試驗

膜片破裂的計算前處理操作非常重要,其決定了計算過程的質(zhì)量好壞,從而直接決定了計算結(jié)果的準(zhǔn)確性.為了驗證計算方法的可行性,在JF-12 復(fù)現(xiàn)風(fēng)洞中進行了激波風(fēng)洞試驗.

根據(jù)Chapman 和Jouguet 提出的CJ 模型[25],爆轟波可以視為一個帶有化學(xué)反應(yīng)的并且不斷釋放能量的激波.可燃?xì)怏w高速傳播的爆轟波是帶有放熱化學(xué)反應(yīng)的強間斷面的傳播過程,且化學(xué)反應(yīng)瞬時完成,Rayleigh 線和Hugoniot 曲線的切點代表著自持爆轟的穩(wěn)定終態(tài),即CJ 條件[26-29].

在點火起爆后,驅(qū)動段的壓力劇烈增大,考慮到實際情況,為了簡化膜片破裂的問題,此時假設(shè)驅(qū)動段的壓力瞬間達(dá)到CJ 爆轟的壓力狀態(tài),即作用在主膜片正面的壓力突變?yōu)镃J 壓力,在該壓力條件作用下,膜片瞬間被打開.后續(xù)工作將考慮采用非定常爆轟壓力條件進行更為深入的破膜特性研究.

2.2 計算結(jié)果驗證

本節(jié)模擬的是爆轟驅(qū)動段與激波管段中的主膜片破裂情況,膜片外徑為540 mm,內(nèi)徑為400 mm,總厚度為3.5 mm,有效厚度為2.3 mm,中心部分開有“十”字形凹槽,凹槽長度為400 mm,寬度為3 mm.在JF-12 復(fù)現(xiàn)風(fēng)洞中進行試驗時,高壓段爆轟初始壓力為2 MPa,低壓段充氣壓力為38 kPa,根據(jù)CJ 爆轟理論,爆轟驅(qū)動后的高壓區(qū)壓力為32 MPa,低壓區(qū)壓力為38 kPa 保持不變,試驗結(jié)果和數(shù)值計算結(jié)果分別如圖9 所示.

圖9 膜片計算結(jié)果驗證(左圖為JF-12 復(fù)現(xiàn)風(fēng)洞試驗結(jié)果,右圖為數(shù)值計算結(jié)果)Fig.9 Verification of diaphragm calculation results(The left is the wind tunnel test results reproduced by JF-12,and the right is the numerical calculation result)

將數(shù)值模擬計算的結(jié)果與JF-12 復(fù)現(xiàn)風(fēng)洞中試驗中的膜片結(jié)果進行定性比較,可以看出膜片的破裂程度基本一致,說明了數(shù)值計算結(jié)果比較準(zhǔn)確,保證了計算結(jié)果的可靠性.因此,可以用數(shù)值計算結(jié)果模擬激波管試驗的結(jié)果,且模擬效果比較理想.

3 膜片破裂過程數(shù)值分析

本節(jié)針對不同總厚度、有效厚度和凹槽長度的鋼膜片,保持膜片寬度為3 mm 不變,通過數(shù)值模擬方法研究了對膜片破膜壓力和有效破膜時間的影響.破膜壓力ΔP指膜片破裂時高、低壓區(qū)的壓力差,有效破膜時間τ0為膜片從開始破裂到完全打開的時間差,總破膜時間t為施加壓力起至完全打開的時間差.

圖10 是鋼膜片計算過程示意圖,膜片總厚度為3.5 mm,有效厚度為2.3 mm,凹槽長度為400 mm.膜片在第544.45 μs 時發(fā)生破裂,在第1100.00 μs 時完成破裂,破裂過程持續(xù)了555.55 μs.膜片首先在凹槽中心處發(fā)生破裂,產(chǎn)生一個小開口,在相同作用力下,凹槽中心處的結(jié)構(gòu)最易破裂.隨著時間推移,開口順著凹槽擴展,膜片凹槽逐漸破裂,最后帶動整個膜片呈現(xiàn)花瓣狀的向外撕裂.膜片破裂過程持續(xù)了一段時間,而非瞬態(tài)破裂,這與實際情況相符,有效破膜時間大約是總破膜時間的一半,比較接近文獻[30] 提供的結(jié)果.

圖10 膜片計算結(jié)果示意圖Fig.10 Schematic diagram of diaphragm calculation results

3.1 膜片有效厚度對破膜影響

膜片有效厚度E通過改變膜片凹槽處的深度來影響膜片的破裂特性,在試驗時,應(yīng)當(dāng)設(shè)置合適的膜片有效厚度保證試驗的可行性.表4 展示了鋼膜片在不同膜片有效厚度時的計算結(jié)果,膜片總厚度為3.5 mm,膜片凹槽長度為360 mm.為了便于觀察分析,將表4 中的計算結(jié)果繪制成坐標(biāo)曲線圖,橫坐標(biāo)為膜片有效厚度,縱坐標(biāo)分別為有效破膜壓力和有效破膜時間,如圖11 和圖12 所示.

表4 不同膜片有效厚度計算結(jié)果Table 4 Calculation results of effective thickness of different diaphragm

圖11 膜片有效厚度對破膜壓力的影響Fig.11 Effect of effective thickness of diaphragm on diaphragm rupture pressure

圖12 膜片有效厚度對有效破膜時間的影響Fig.12 Effect of effective thickness of diaphragm on effective diaphragm rupture time

圖11 和圖12 展示了鋼膜片在膜片總厚度和凹槽長度一定時,其破裂情況隨膜片有效厚度的變化情況.隨著膜片有效厚度的增大,膜片的有效破膜壓力增加,有效破膜時間減少.其原因是當(dāng)膜片有效厚度增加時,即膜片凹槽深度減小,此時需要更大的破膜壓力將膜片凹槽處的材料擊破,隨后再帶動整個膜片破裂.當(dāng)膜片凹槽中心破裂時,隨著破膜壓力的增加,此時會加快整個膜片破裂的過程,膜片持續(xù)打開的時間減少,即膜片的有效破膜時間會隨著膜片有效厚度的增大而縮短.

3.2 膜片總厚度對破膜影響

在膜片凹槽打開后,膜片的總厚度T通過改變瓣膜的厚度來影響破膜壓力,膜片的總厚度對激波管試驗也會產(chǎn)生比較大的影響.試驗中,膜片的總厚度過大會導(dǎo)致膜片在短時間內(nèi)不易破裂,而過薄的膜片無法承受巨大的爆轟脈沖壓力,容易脫落對流場產(chǎn)生不利的影響.因此,膜片的總厚度選擇十分重要,應(yīng)該選擇合適的膜片厚度進行試驗.考慮實際情況,在計算過程中,保持膜片無量綱厚度E/T(即膜片有效厚度與膜片總厚度的比值)不變,改變膜片總厚度進行研究.表5 展示了無量綱厚度為0.43 和0.57時鋼膜片的計算結(jié)果,膜片凹槽長度為400 mm.同樣,將表5 中的計算結(jié)果繪制成坐標(biāo)曲線圖,如圖13和圖14 所示.

表5 不同膜片總厚度計算結(jié)果Table 5 Calculation results of different diaphragm total thickness

圖13 和圖14 展示的是當(dāng)無量綱厚度為0.43 和0.57 時,膜片總厚度對破膜壓力和有效破膜時間的影響.發(fā)現(xiàn)有效破膜時間與膜片總厚度呈現(xiàn)正相關(guān)線性關(guān)系,當(dāng)膜片總厚度越大時,其膜片有效厚度也越大,需要更大的破膜壓力將膜片凹槽處的材料擊破,同時需要更大的破膜壓力將膜片瓣膜完全打開,因此,破膜壓力隨著膜片總厚度的增大而增大.由于膜片總厚度和膜片有效厚度的增大,整個膜片的打開過程變得緩慢,所以膜片有效時間隨著膜片總厚度呈現(xiàn)正相關(guān)線性關(guān)系.

圖13 膜片總厚度對破膜壓力的影響Fig.13 Effect of total thickness of diaphragm on diaphragm rupture pressure

圖14 膜片總厚度對有效破膜時間的影響Fig.14 Effect of total thickness of diaphragm on effective diaphragm rupture time

在圖13 和圖14 中,曲線變化趨勢一致,當(dāng)膜片總厚度一定時,無量綱厚度變大,此時的膜片破膜壓力會增加,而有效破膜時間會減小.造成這一現(xiàn)象的原因與膜片有效厚度影響規(guī)律一致,即膜片有效厚度的增大會導(dǎo)致破膜壓力的增大和有效破膜時間的縮短.

3.3 膜片凹槽長度對破膜影響

在目前關(guān)于激波管膜片的研究中,尚未有關(guān)于膜片凹槽長度對膜片破裂的影響規(guī)律.凹槽長度L通過改變膜片凹槽的打開時間來影響膜片的破膜壓力和有效破膜時間.表6 展示了鋼膜片凹槽長度對膜片破裂特性的影響,膜片總厚度為3.5 mm,無量綱厚度分別為0.43 和0.57.將表6 中的計算結(jié)果進行處理,繪制成坐標(biāo)曲線圖,如圖15 和圖16 所示.

圖15 膜片凹槽長度對破膜壓力的影響Fig.15 Effect of diaphragm groove length on diaphragm rupture pressure

圖16 膜片凹槽長度對有效破膜時間的影響Fig.16 Effect of diaphragm groove length on effective diaphragm rupture time

在圖15 和圖16 中,當(dāng)膜片總厚度和無量綱厚度一定時,膜片的破膜壓力與膜片凹槽長度呈現(xiàn)負(fù)相關(guān)線性關(guān)系,膜片有效破膜時間與凹槽長度呈現(xiàn)正相關(guān)線性關(guān)系.膜片凹槽長度增大時,增大的凹槽處有效厚度減小,相應(yīng)所需要的破膜壓力也減小,所以破膜壓力會相應(yīng)減小.但是膜片的總厚度保持不變,破膜壓力的減小會減緩膜片瓣膜打開的過程,所以相對應(yīng)的有效破膜時間增加.

3.4 膜片尺寸優(yōu)化問題

在激波風(fēng)洞中進行試驗時,當(dāng)膜片選擇確定后,根據(jù)膜片破裂的特性規(guī)律可以有效預(yù)測膜片的破膜壓力和有效破膜時間.同時,當(dāng)激波風(fēng)洞中的破膜壓力一定時,研究激波管中的膜片破裂特性有助于風(fēng)洞試驗中膜片的選取.在激波風(fēng)洞中,我們希望膜片在破裂時能夠順利打開,而且打開的過程越快越好,理想情況下我們認(rèn)為膜片瞬間完全打開,因此膜片總破膜時間越短越好.同時,膜片破裂的過程中不能出現(xiàn)脫落的情況,保持破膜壓力不變,不能影響激波風(fēng)洞氣流品質(zhì).針對JF-12 復(fù)現(xiàn)風(fēng)洞,根據(jù)膜片破裂的特性提出了一種設(shè)計方案,適合采用鋼作為膜片材料,如表7 所示.

表7 JF-12 復(fù)現(xiàn)風(fēng)洞膜片優(yōu)化方案Table 7 Diaphragm optimization scheme of JF-12 reproduction wind tunnel

膜片外形因素在一定程度上影響膜片的破裂情況,針對具體風(fēng)洞試驗,需要根據(jù)實際情況判斷破膜壓力的大小.在破膜壓力較大時,需要選取相應(yīng)較厚而且凹槽長度較小的膜片.一般總破膜時間較短,因此常選用總厚度適中,有效厚度較小,凹槽長度較小的鋼膜片進行激波管試驗.

4 結(jié)論

在激波風(fēng)洞試驗前,膜片將激波管分為兩段,將高低壓氣體隔開,創(chuàng)造相互獨立的氣體環(huán)境.在試驗時,膜片在脈沖壓力的作用下破裂完全打開,在被驅(qū)動段產(chǎn)生運動激波,創(chuàng)造所需的來流條件.傳統(tǒng)的實驗研究鮮有學(xué)者研究爆轟驅(qū)動的激波風(fēng)洞中膜片的破裂特性.本文采用有限元分析軟件模擬計算了膜片破裂過程,計算結(jié)果與試驗結(jié)果吻合得比較理想.基于有限元數(shù)值模擬技術(shù)研究了膜片外形參數(shù)(膜片有效厚度、總厚度和凹槽長度) 對于膜片破裂過程的影響,并給出了一種激波管用膜片的優(yōu)化設(shè)計方案.

根據(jù)數(shù)值計算結(jié)果,膜片首先從凹槽中心處開始破裂,然后開口順著凹槽破裂,隨后帶動整個膜片瓣膜破裂.在鋼膜片總厚度一定時,破膜壓力與膜片有效厚度呈現(xiàn)正相關(guān)線性關(guān)系,有效破膜時間與膜片有效厚度近似呈現(xiàn)負(fù)相關(guān)線性關(guān)系.在鋼膜片無量綱厚度一定時,破膜壓力和有效破膜時間均與膜片總厚度呈現(xiàn)正相關(guān)線性關(guān)系.破膜壓力與膜片凹槽長度成負(fù)相關(guān)線性關(guān)系,有效破膜時間與膜片凹槽長度成正相關(guān)線性關(guān)系.

研究激波風(fēng)洞中主膜片的破裂特性將有助于預(yù)測破膜壓力和有效破膜時間,便于激波風(fēng)洞中參數(shù)設(shè)置,為激波風(fēng)洞試驗提供計算依據(jù),具有現(xiàn)實意義和工程實踐意義.