曹銀萍 李明飛 竇益華 于洋
西安石油大學(xué)
大孔徑、高孔密、深穿透射孔是目前射孔工藝的發(fā)展方向,它將增加射孔彈的裝藥密度、裝藥量與爆轟能量,加劇對套管和射孔槍的損傷。國內(nèi)西部油田和美國 Fort Worth 盆地[1]、Arkoma 盆地[2]以及墨西哥灣[3]等油田都曾遇到過射孔事故。射孔是包含能量、質(zhì)量、速度及其變形轉(zhuǎn)化的復(fù)雜動態(tài)過程,如何查明之間復(fù)雜的轉(zhuǎn)換、轉(zhuǎn)化和相互影響規(guī)律,如何實(shí)現(xiàn)真正意義的爆轟沖擊下液固相變轉(zhuǎn)化和流固耦合,為射孔槍和套管強(qiáng)度安全性分析提供依據(jù)是亟待解決的關(guān)鍵問題。
針對聚能射流對射孔槍的動態(tài)傷害,部分文獻(xiàn)給出了數(shù)值分析結(jié)果,文獻(xiàn)[4]應(yīng)用開孔彈塑性力學(xué)理論,建立了射孔套管單孔應(yīng)力分布模型,推導(dǎo)了應(yīng)力集中系數(shù),得到了整個(gè)圓孔周圍應(yīng)力分布規(guī)律。文獻(xiàn)[5]采用ALE (Arbitrary Lagrange-Euler)法,模擬了射流形成過程及其對混凝土靶板的侵徹,分析了藥型罩錐角對射流形成和侵徹的影響,未考慮對套管、射孔槍毛刺及其應(yīng)力和強(qiáng)度的影響。文獻(xiàn)[6-7]應(yīng)用二維有限元法,分析了爆轟波疊加及干擾對偏轉(zhuǎn)角的影響。文獻(xiàn)[8-9]研究了ALE 算法在相關(guān)領(lǐng)域的應(yīng)用,給出了研究方法和技巧,為進(jìn)一步研究提供了借鑒和參考。文獻(xiàn)[10]研究了非穩(wěn)態(tài)噴射流通過射孔孔眼對套管強(qiáng)度的影響。文獻(xiàn)[11]通過炸藥做功和猛度實(shí)驗(yàn),驗(yàn)證了爆轟波碰撞的聚能效果。文獻(xiàn)[12]采用有限元和自適應(yīng)網(wǎng)格技術(shù),分析了射流和碎片撞擊射孔槍壁的過程,分析了藥型罩上附金屬對脹槍高度的影響。
目前針對爆轟波沖擊、疊加對能量和射流速度及套管強(qiáng)度的影響研究較少。文獻(xiàn)[13]結(jié)合室內(nèi)材料試驗(yàn)及有限元法,分析了承受非均勻載荷套管的力學(xué)性能。文獻(xiàn)[14]應(yīng)用SPH(Smoothed Particle Hydrodynamics)法,建立了不同藥型罩射流模型,考察了不同藥型罩材料形成金屬射流的過程、速度、性質(zhì)等。文獻(xiàn)[15]利用LS-DYNA 軟件,模擬了不同裝藥類型的聚能射流成型及侵徹過程,分析了對均質(zhì)裝甲鋼的侵徹。文獻(xiàn)[16]采用AUTODYNA 軟件,模擬了Y1201 射孔彈壓垮藥型罩和侵徹套管的過程。文獻(xiàn)[17]基于Hamilton 原理,建立一個(gè)包含射流強(qiáng)度、應(yīng)變率效應(yīng)、剪切、表面張力、動力學(xué)黏性和速度梯度等多耦合的金屬射流運(yùn)動方程。文獻(xiàn)[18]應(yīng)用ALE 法,對聚能裝藥在不同起爆方式下的射流形成過程進(jìn)行了三維數(shù)值分析。
綜上所述,目前針對爆轟波沖擊、疊加對射孔槍和套管強(qiáng)度的影響研究較少。常用高速柱狀體代替金屬射流,未模擬爆炸和射流過程。侵徹是具有長度的一段金屬射流,具有連續(xù)和持續(xù)性,簡單用柱狀體代替射流存在差異,影響分析精度。部分研究為二維建模,往往忽略射孔槍和盲孔的存在,而射孔槍盲孔將消耗掉一部分射流能量,并且二維模型無法模擬射孔彈間的互相干擾和影響。還有部分研究直接給出爆轟過程的關(guān)鍵參數(shù),沒有模擬真實(shí)爆轟,參數(shù)與實(shí)際情況匹配性差,同時(shí)也無法分析爆轟瞬間各部分間的能量分配。軍工和火工主要針對武器和礦、煤開采等方面,工況、初始和邊界條件與本次研究的初始高溫和高壓存在很大區(qū)別。為了分析爆轟波疊加對能量和射流速度及套管和射孔槍強(qiáng)度的影響規(guī)律,研究應(yīng)用ALE 方法,建立射孔彈-射孔槍(盲孔)-套管三維有限元模型,實(shí)現(xiàn)爆轟過程、藥型罩固流轉(zhuǎn)化、射流侵徹套管的大變形和流固耦合仿真,提高分析精度,為射孔槍安全優(yōu)化設(shè)計(jì)和爆轟沖擊下套管的力學(xué)性能分析提供思路。
應(yīng)用LS-DYNA 軟件的ALE 方法,建立射孔彈-射孔槍-套管三維有限元模型,具體參數(shù)見表1。模型中炸藥和藥型罩之間,藥型罩、射孔槍、套管與氣體之間的接觸采用Surface to Surf 和Eroding (ESTS),即侵徹接觸類型。由于藥型罩在炸藥作用下形成聚能射流,在侵徹射孔槍和套管的過程中存在大變形和高速流動,采用Lagrange 算法會出現(xiàn)單元畸變的現(xiàn)象,而Euler 算法為了精確捕捉固體材料的變形響應(yīng)需要很精細(xì)的網(wǎng)絡(luò),極大增加了數(shù)值分析的成本。ALE 算法(Arbitrary Lagrange-Euler)則可以很好地解決大變形和流固耦合問題。因此,空氣、炸藥、藥型罩的單元算法為ALE 算法,空腔、套管、射孔槍采用Lagrange 單元算法。
表1 有限元建模關(guān)鍵參數(shù)Table 1 Key parameters of finite element modeling
為了在確保得到準(zhǔn)確模擬結(jié)果的同時(shí)盡量縮短計(jì)算時(shí)間,需要控制單元大小和總數(shù)。為此,細(xì)化了射孔彈及其周圍的模型。首先,在粗網(wǎng)格上采用Oseen 迭代法求解非線性問題;然后,在細(xì)網(wǎng)格上分別并行求解Oseen、Newton、Stokes 線性問題以校正粗網(wǎng)格解,對于空間變量采用有限元離散,時(shí)間變量采用向后Euler 格式離散。所建單枚射孔彈系統(tǒng)有限元模型中共有1 545 678 個(gè)節(jié)點(diǎn),1 465 987 個(gè)單元。如圖1 和圖2 所示為單枚和3 枚射孔彈起爆、壓垮、射流、侵徹射孔槍和套管的全過程。爆轟給藥型罩施加高溫高壓,形成拋擲,在軸線上激烈碰撞,形成細(xì)束金屬射流和密實(shí)完整的杵體,速度梯度使得射流被拉長,最終斷裂成局部杵體。LS-DYNA軟件的ALE 方法是近期分析爆轟和超高速撞擊、侵徹的一個(gè)新方法,文獻(xiàn)[5-7]已經(jīng)初步驗(yàn)證了該方法的可行性,筆者目前在進(jìn)行2 級氫氣炮加速藥型罩彈丸至7 000 m/s 的侵徹套管、水泥環(huán)、圍巖板靶實(shí)驗(yàn),將進(jìn)一步驗(yàn)證該方法的可行性和可靠性。
圖1 單枚射孔彈射孔動態(tài)過程Fig.1 Perforating behavior process of one perforating charge
圖2 3 枚射孔彈射孔動態(tài)過程Fig.2 Perforating behavior process of three perforating charges
基于ALE 聚能射流和流體化的射孔槍盲孔材料動力學(xué)方程,可在Euler 描述下流體動力學(xué)方程基礎(chǔ)上推導(dǎo)得到。ALE 描述與Euler 描述的本質(zhì)區(qū)別是材料時(shí)間導(dǎo)數(shù)項(xiàng)的不一致。
在Euler 描述中,聚能射流和流體化的射孔槍盲孔材料連續(xù)性方程,即質(zhì)量守恒方程為
式中, ρ為流體密度,kg/m3;x為 Euler 坐標(biāo)系下流體的坐標(biāo),m; ?x·為 空間散度;v為流體速度,m/s; ρt[x]為流體坐標(biāo)x固定時(shí),對材料密度求時(shí)間的導(dǎo)數(shù); ?x為空間梯度。
將流體時(shí)間導(dǎo)數(shù)表達(dá)為固定坐標(biāo)下ALE 對時(shí)間的偏導(dǎo)數(shù)和空間梯度的形式為
式中,f為場函數(shù);t[χ]為 流體時(shí)間,s; χ為ALE 坐標(biāo)系下流體的坐標(biāo),m;c為爆炸產(chǎn)物中的聲速,m/s。
應(yīng)用ALE 材料時(shí)間導(dǎo)數(shù)(2)式替換(1)式中的材料時(shí)間導(dǎo)數(shù)項(xiàng),可得ALE 下的連續(xù)性方程為
Euler 描述中,流體化的射孔槍盲孔材料動量守恒方程為
式中,vt[x]為流體坐標(biāo)x固定時(shí),對流體速度求時(shí)間的導(dǎo)數(shù); σ為應(yīng)力張量;b為單位質(zhì)量體力向量。同理,用(2)式替換(4)式中的材料時(shí)間導(dǎo)數(shù)項(xiàng),可得ALE 描述下的動量守恒方程為
Euler 描述中,流體化的射孔槍材料能量守恒方程為
式中,e為單位質(zhì)量內(nèi)能,J/kg;et[x]為流體坐標(biāo)x固定時(shí),對單位質(zhì)量內(nèi)能求時(shí)間的導(dǎo)數(shù);D為變形率;q為單位面積熱流矢量,W/m2;ρs為單位體積內(nèi)的熱源,W/m3。
應(yīng)用(2)式替換(6)式中的材料時(shí)間導(dǎo)數(shù)項(xiàng),可得ALE 描述下的能量守恒方程為
網(wǎng)格運(yùn)動是通過控制網(wǎng)格運(yùn)動的速度v1或控制參考流體運(yùn)動的速度w來實(shí)現(xiàn)。參考材料速度為
式中,v1為網(wǎng)格運(yùn)動速度,m/s;w為流體運(yùn)動速度,m/s;X為Lagrange 坐標(biāo)系下流體的坐標(biāo); χ(X,t)為某時(shí)刻流體坐標(biāo)X固定情況下的流體坐標(biāo); χt[X]為流體坐標(biāo)X固定時(shí),對流體坐標(biāo)求時(shí)間的導(dǎo)數(shù)。
由式(8)、(9),可得到ALE 描述下流固耦合動力學(xué)方程為
所以通過控制w可以達(dá)到控制v1的目的。應(yīng)用LS-DYNA 軟件,結(jié)合ALE 算法的內(nèi)核,即可實(shí)現(xiàn)爆轟過程、藥型罩固流轉(zhuǎn)化、射流高速侵徹射孔槍和套管的大變形和流固耦合仿真。
在數(shù)值分析過程中,記錄了不同時(shí)刻射孔槍與套管的開坑、擴(kuò)孔過程,孔徑變化過程,并分析了孔徑變化的規(guī)律。
如圖3 所示,13 μs 時(shí)金屬射流頭部接觸射孔槍盲孔內(nèi)壁,14 μs 時(shí)射流頭部到達(dá)盲孔外壁。14~16 μs處于開坑階段,孔眼直徑為3.2 mm。射穿射孔槍后,后繼射流持續(xù)侵徹射孔槍,35 μs 時(shí),尾部通過射孔槍外壁。16~35 μs 處于擴(kuò)孔階段,孔眼直徑由4 mm迅速增大到13 mm。因杵體尾部掃過射孔槍內(nèi)壁,孔徑由17.1 mm 增大到24.4 mm,平均孔徑為18.7 mm。
圖3 射孔槍內(nèi)外壁孔徑變化規(guī)律Fig.3 Variation law of the perforation diameter in the inner and outer wall of perforating gun
如圖4 所示,15 μs 時(shí)射流開始接觸套管內(nèi)壁,23 μs 時(shí)射流到達(dá)套管外壁并穿透套管。15~23 μs為開坑階段,孔眼直徑達(dá)到3 mm。后續(xù)射流持續(xù)侵徹套管,30 μs 時(shí)尾部通過套管外壁,侵徹結(jié)束。23~30 μs 為擴(kuò)孔階段,套管外壁孔徑由 3 mm 迅速增到7.8 mm,內(nèi)壁孔徑7.6 mm,平均孔徑7.7 mm。對比圖3 和圖4 可知,杵體對套管孔徑和孔道形態(tài)影響較大,造成內(nèi)壁孔徑增大,形成凹坑。射流穿透射孔槍后能量衰減明顯,使套管較射孔槍孔眼直徑小11 mm。
圖4 套管內(nèi)外壁孔徑變化規(guī)律Fig.4 Variation law of the perforation diameter in inner and outer wall of casing
如圖5(a)所示,13 μs 時(shí)射流頭部接觸射孔槍,射孔槍最大應(yīng)力達(dá)到1.055 GPa;如圖5(b)所示,14 μs時(shí)射流頭部到達(dá)射孔槍外壁,應(yīng)力最大值1.084 GPa;40 μs 時(shí)射流完全侵徹射孔槍,孔眼的周圍有圓環(huán)狀應(yīng)力等值分布,最大應(yīng)力不變。
圖5 單顆射孔彈射孔時(shí)射孔槍應(yīng)力云圖Fig.5 Stress cloud map of the perforating gun during the perforation of one perforating charge
圖6 記錄了射孔槍孔眼近似徑向8 個(gè)點(diǎn)的應(yīng)力值,距孔邊距離22 mm 內(nèi)點(diǎn)的應(yīng)力均超過了材料屈服極限835 MPa,距離30 mm 內(nèi)點(diǎn)的應(yīng)力值也達(dá)到500 MPa,距離58 mm 內(nèi)點(diǎn)的應(yīng)力值接近0。距孔邊距離0~22 mm 范圍內(nèi)射孔槍的材料類似于產(chǎn)生了“冷作硬化”現(xiàn)象,該區(qū)域內(nèi)的材料不再符合原材料的應(yīng)力應(yīng)變規(guī)則,材料強(qiáng)度會略有提高,但韌性降低,材料變脆,沿著卸載曲線變化。
圖6 單顆射孔彈射孔時(shí)射孔槍最大應(yīng)力變化規(guī)律Fig.6 Variation law of the maximum stress of the perforating gun during the perforation of one perforating charge
如圖7(a)所示,13 μs 時(shí)射流侵徹射孔槍,應(yīng)力瞬間增大。如圖7(b)所示,14 μs 時(shí)射流頭部到達(dá)射孔槍外壁,應(yīng)力最大值達(dá)到1 084 MPa;截至此時(shí),應(yīng)力變化與單枚射孔彈時(shí)一致。如圖7(c)所示,21 μs 時(shí)孔眼附近的高應(yīng)力區(qū)持續(xù)擴(kuò)展,還未疊加。如圖7(d)所示,40 μs 時(shí)高應(yīng)力區(qū)范圍疊加、重合,應(yīng)力增大;在三孔相連的長200 mm,寬125 mm 的寬帶范圍內(nèi),其最小應(yīng)力也達(dá)到了774 MPa,接近材料的屈服強(qiáng)度;在高爆轟壓力的作用下,該寬帶會產(chǎn)生向外的鼓脹,即會發(fā)生所謂的“脹槍”。
如圖8 所示,3 枚射孔彈射孔時(shí),孔邊最大應(yīng)力達(dá)到1 076 MPa,與單枚時(shí)基本相符;距離孔邊15 mm以外,單枚射孔的射孔槍應(yīng)力急劇降低;3 枚射孔的射孔槍應(yīng)力增加,在爆轟重疊中心,射孔槍應(yīng)力達(dá)到996 MPa,較單枚時(shí)的878 MPa 提高了118 MPa。
射孔彈爆轟瞬間,高溫高壓高速的金屬射流會侵徹、損傷套管,影響后續(xù)工況射孔段套管的強(qiáng)度安全。
如圖9(a)所示,17 μs 時(shí)射流開始侵徹套管,套管內(nèi)壁與射流接觸區(qū)瞬間產(chǎn)生1.04 GPa 應(yīng)力。如圖9(b)所示,17 μs 時(shí)套管外壁和除接觸點(diǎn)外的大部區(qū)域應(yīng)力為0;繼續(xù)侵徹、開孔,30 μs 時(shí)射穿套管。如圖9(c)所示,30~40 μs 之間屬于射穿套管后的擴(kuò)孔階段。套管外壁受射流影響區(qū)域進(jìn)一步擴(kuò)展,直至如圖9(d)所示,侵徹套管結(jié)束時(shí),距孔邊緣10 mm范圍內(nèi)環(huán)形區(qū)域內(nèi)應(yīng)力可以達(dá)到1.27 GPa,其邊緣應(yīng)力也達(dá)到了965 MPa,超過了TP140 套管的屈服極限,相當(dāng)于產(chǎn)生了“冷作硬化”,射孔結(jié)束后該區(qū)域內(nèi)的材料不再遵循原材料的應(yīng)力應(yīng)變規(guī)則,屈服強(qiáng)度將略有提高,對套管強(qiáng)度的影響有待進(jìn)一步研究。
圖7 3 顆射孔彈射孔時(shí)射孔槍應(yīng)力云圖Fig.7 Stress cloud map of the perforating gun during the perforation of three perforating charges
圖8 單枚和3 枚射孔彈射孔時(shí)套管應(yīng)力變化曲線Fig.8 Variation of the casing stress during the perforation of one perforating charge and three perforating charges
由圖9(d)可知,40 μs 時(shí)套管上仍然存在應(yīng)力超過965 MPa 的危險(xiǎn)區(qū)域,因此記錄并考察了套管孔眼徑向上與孔眼不同距離的12 個(gè)點(diǎn)的應(yīng)力變化情況。如圖10 所示,距孔邊10 mm 內(nèi)所有點(diǎn)的應(yīng)力均超過套管材料的屈服強(qiáng)度;在距離孔眼0~8.3 mm范圍內(nèi),計(jì)算點(diǎn)的最大應(yīng)力降低相對緩慢,但是在8.3 mm 以后,測量點(diǎn)的應(yīng)力迅速下降;距離孔邊越近,其應(yīng)力值落在屈服極限上方的范圍越大。
圖9 單枚射孔彈套管應(yīng)力云圖Fig.9 Stress cloud map of the casing during the perforation of one perforating charge
圖10 套管最大應(yīng)力變化規(guī)律Fig.10 Variation law of the maximum casing stress
如圖11 所示,3 枚射孔彈爆轟不同時(shí)段套管的應(yīng)力云圖。16 μs 時(shí)金屬射流接觸套管內(nèi)壁,較單枚射孔彈提前了 1 μs;22 μs 時(shí)射穿套管,較單枚射孔彈提前了8 μs,此時(shí)孔眼處最大應(yīng)力達(dá)到1.30 GPa,較單枚射孔彈增加了 0.26 Gpa;27 μs 時(shí),應(yīng)力變化區(qū)域的圓環(huán)開始相切,說明此時(shí)應(yīng)力變化范圍邊緣距孔邊距離為21.3 mm。隨后應(yīng)力變化區(qū)域繼續(xù)擴(kuò)展、互相重疊,如圖11(d)所示,重疊區(qū)域的應(yīng)力值較非重疊部分有所提高。
圖11 3 枚射孔彈射孔時(shí)套管應(yīng)力變化云圖Fig.11 Stress cloud map of the casing during the perforation of three perforating charges
為考察兩孔間套管應(yīng)力的變化規(guī)律,針對兩孔間連線上的13 個(gè)點(diǎn)的應(yīng)力值開展研究。如圖12 所示,孔距(兩孔連線中點(diǎn))在12.5 mm 范圍內(nèi),其應(yīng)力值均超過材料屈服強(qiáng)度965 MPa,較單枚射孔彈的范圍擴(kuò)大了2.5 mm 寬度的圓環(huán);孔距在12.5~15.3 mm 之間,套管應(yīng)力隨孔距增加迅速降低。因射孔彈間干擾導(dǎo)致應(yīng)力變化區(qū)重疊,孔距在15.3~21.3 mm 范圍內(nèi),套管應(yīng)力呈增加趨勢,最大應(yīng)力達(dá)965 MPa。因應(yīng)力增加又泄壓,15.3~21.3 mm 的重疊區(qū)域的力學(xué)屬性將發(fā)生變化,強(qiáng)度會略有提高,但韌性降低,材料變脆,將影響套管的強(qiáng)度安全。
圖12 套管相鄰兩孔間應(yīng)力變化曲線Fig.12 Variation curve of the stress between two adjacent perforations in the casing
如圖13 所示單枚和3 枚射孔彈射孔時(shí)套管應(yīng)力變化對比曲線。3 枚射孔彈射孔套管的最大應(yīng)力為1 273 MPa,比單枚時(shí)的1 207 MPa 略高;距離孔邊大于8 mm,單枚射孔彈射孔的套管應(yīng)力衰減速度明顯大于3 枚時(shí)套管應(yīng)力衰減速度;單枚射孔彈射孔套管應(yīng)力持續(xù)減少,3 枚時(shí)套管在距離孔邊超過15.3 mm 以上時(shí),套管應(yīng)力增加。在爆轟重疊區(qū)域中心,套管應(yīng)力為965 MPa,較單枚時(shí)的226 MPa提高了737 MPa。
圖13 單枚和3 枚射孔彈射孔時(shí)套管應(yīng)力變化曲線Fig.13 Variation of the casing stress during the perforation of one perforating charge and three perforating charges
(1)應(yīng)用LS-DYNA 軟件,建立了射孔彈-射孔槍(含盲孔)-套管的三維有限元模型,采用ALE 算子分離算法,實(shí)現(xiàn)了爆轟過程、射流高速侵徹套管、藥型罩固流轉(zhuǎn)化的大變形和流固耦合仿真,分析了射孔槍和套管的強(qiáng)度安全性。
(2)射孔槍內(nèi)外壁孔眼直徑分別為24.4 mm 和13 mm,平均孔徑為18.7 mm。套管內(nèi)外壁孔眼直徑分別為7.6 mm 和7.8 mm,平均孔徑為7.7 mm,相差11 mm,說明穿透射孔槍盲孔過程消耗了大量能量,應(yīng)合理設(shè)計(jì)盲孔壁厚。
(3)射孔槍的高應(yīng)力區(qū)疊加、重合,使得此區(qū)域內(nèi)材料應(yīng)力增大,在射孔槍各孔相連的長250 mm、寬125 mm 的寬帶內(nèi),其最小應(yīng)力也達(dá)到了774 MPa,這一條寬帶應(yīng)力值均超過其屈服強(qiáng)度,在高爆轟壓力的作用下,會產(chǎn)生向外的鼓脹。
(4)爆轟波疊加作用區(qū)域和孔邊以20 mm 為半徑的圓環(huán)應(yīng)力變化區(qū)為射孔爆轟引起的套管應(yīng)力加強(qiáng)區(qū),相當(dāng)于材料產(chǎn)生了“冷作硬化”,射孔結(jié)束后該區(qū)域范圍內(nèi)的材料不再符合原材料的應(yīng)力應(yīng)變曲線規(guī)則,材料力學(xué)屬性會發(fā)生變化,強(qiáng)度會略有提高,但韌性降低,材料變脆,影響套管的強(qiáng)度安全。