• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    空氣不耦合裝藥孔壁初始沖擊壓力的計(jì)算

    2020-07-23 00:36:48張恒根王衛(wèi)華王永強(qiáng)
    工程爆破 2020年3期
    關(guān)鍵詞:孔壁炮孔裝藥

    張恒根,王衛(wèi)華,王永強(qiáng)

    (中南大學(xué)資源與安全工程學(xué)院,長(zhǎng)沙 410083)

    在確定不耦合裝藥系數(shù)等爆破參數(shù)時(shí),孔壁初始沖擊壓力是重要的依據(jù)之一,對(duì)它的研究吸引了眾多學(xué)者的注意。在以往的研究中,有些學(xué)者主張利用計(jì)算爆轟產(chǎn)物壓力的方法得到不耦合裝藥時(shí)的孔壁壓力,如費(fèi)鴻祿等[1]以爆轟氣體動(dòng)力學(xué)及巖石斷裂理論為基礎(chǔ),推導(dǎo)了裝藥不耦合系數(shù)的計(jì)算公式。宗琦等[2]將爆生氣體膨脹充滿炮孔時(shí)的準(zhǔn)靜壓力作為孔壁初始?jí)毫?,得到了空氣墊層裝藥結(jié)構(gòu)軸向不耦合系數(shù)的計(jì)算方法。也有學(xué)者主張利用計(jì)算沖擊波壓力的方法得到不耦合裝藥時(shí)的孔壁壓力,如朱紅兵等[3]利用沖擊波理論分析了空氣間隔裝藥炮孔內(nèi)壓力波的傳播過(guò)程,得到了合理空氣層比例的確定方法。陳明等[4]推導(dǎo)了爆炸沖擊波與彈性壁碰撞后壓力增大倍數(shù)的理論解,提出了輪廓爆破孔壁壓力峰值的計(jì)算模型。杜俊林等[5]利用沖擊波的基本關(guān)系式計(jì)算了空氣不耦合裝藥時(shí)的孔壁壓力。還有學(xué)者用模擬及實(shí)驗(yàn)等其他方法對(duì)孔壁壓力進(jìn)行了研究,如王志亮等[6]對(duì)水不耦合裝藥爆破進(jìn)行了數(shù)值模擬,分析了孔壁壓力與不耦合系數(shù)之間的關(guān)系。V.R.Feldgun等[7]對(duì)爆破荷載的變化過(guò)程進(jìn)行了數(shù)值模擬,并與實(shí)驗(yàn)結(jié)果進(jìn)行了分析對(duì)比。凌偉明[8]采用錳銅壓阻傳感器進(jìn)行了孔壁壓力測(cè)量的實(shí)驗(yàn),驗(yàn)證了錳銅壓阻傳感器用于測(cè)試孔壁壓力的可行性。劉云川等[9]按照能量等效原則,給出了徑向不耦合裝藥條件下炮孔初始?jí)毫Φ挠?jì)算公式。樓曉明等[10]引入St-arfield迭加法,分析得到了上部空氣間隔裝藥和中部空氣間隔裝藥條件下孔壁初始沖擊壓力的計(jì)算公式。眾多其他學(xué)者[11-14]也研究了孔壁初始沖擊壓力問(wèn)題,提出了許多計(jì)算孔壁壓力的方法。

    目前常用的計(jì)算不耦合裝藥孔壁初始沖擊壓力的方法是,先計(jì)算爆轟產(chǎn)物在炮孔內(nèi)做等熵膨脹后作用于孔壁的準(zhǔn)靜壓力,然后將準(zhǔn)靜壓力乘以壓力增大系數(shù)n,近似得到了孔壁巖石所受到的初始?jí)毫?。n一般取8~11。爆轟產(chǎn)物在炮孔內(nèi)的膨脹可以按照一階段或兩階段等熵膨脹計(jì)算[15]。當(dāng)不耦合系數(shù)比較小時(shí),通常按照一階段等熵膨脹計(jì)算:

    (1)

    當(dāng)不耦合系數(shù)比較大時(shí),按照兩階段等熵膨脹計(jì)算:

    (2)

    式(1~2)中:pt為透射入孔壁的初始沖擊壓力;ρw為炸藥的密度;D為爆轟波波陣面的速度;R為炮孔半徑;r為藥卷半徑;pk為兩階段交接點(diǎn)處的臨界壓力,中等威力的硝銨類炸藥pk一般取2.0×108Pa。

    上述常用公式在理論上并不嚴(yán)密,其計(jì)算結(jié)果與實(shí)測(cè)值往往相差較大。以往對(duì)于不耦合裝藥孔壁初始沖擊壓力計(jì)算方法的研究始終圍繞著不耦合裝藥的破巖機(jī)理,即由沖擊波破巖和爆生氣體破巖的理論展開(kāi)。由于確定爆破近區(qū)空氣沖擊波衰減和爆生氣體膨脹規(guī)律的困難性,以及爆轟波、空氣沖擊波與孔壁之間相互作用的復(fù)雜性和測(cè)試系統(tǒng)的局限性等原因,致使空氣不耦合裝藥條件下孔壁初始沖擊壓力較難確定。以往利用計(jì)算爆轟產(chǎn)物壓力的方法計(jì)算孔壁壓力時(shí),往往事先假定忽略炮孔內(nèi)空氣層的存在,進(jìn)而忽略空氣沖擊波對(duì)孔壁的影響,這與事實(shí)不符,因?yàn)槔糜?jì)算沖擊波壓力的方法確定孔壁壓力時(shí),往往將孔壁面作彈性面簡(jiǎn)化處理,且較少考慮空氣沖擊波傳播過(guò)程中的衰減。作者在前人研究的基礎(chǔ)上,通過(guò)分析不耦合裝藥爆轟后沖擊波作用于孔壁的物理過(guò)程,建立了空氣沖擊波在炮孔內(nèi)衰減傳播模型,利用該模型計(jì)算了空氣沖擊波對(duì)孔壁的入射壓力,再利用空氣沖擊波正入射時(shí)的界面處連續(xù)條件,獲得了孔壁初始沖擊壓力的計(jì)算公式,并將計(jì)算結(jié)果與常用公式的計(jì)算值、模擬值和實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比分析。

    1 物理過(guò)程

    為了比較合理地確定孔壁初始?jí)毫?,首先要清楚爆轟后沖擊波作用于孔壁的物理過(guò)程??諝獠获詈涎b藥爆轟后,空氣沖擊波隨時(shí)間變化作用于孔壁的過(guò)程如圖1所示。

    注:(1)區(qū)為未爆轟的藥卷;(2)區(qū)為炮孔內(nèi)的空氣層;(3)區(qū)為爆轟后還未膨脹的爆轟產(chǎn)物;(4)區(qū)為受膨脹波擾動(dòng)后開(kāi)始膨脹的爆轟產(chǎn)物;(5)區(qū)為空氣沖擊波過(guò)后被壓縮的空氣。圖1 空氣沖擊波隨時(shí)間變化作用于孔壁的過(guò)程Fig.1 The process of air shock wave acting on the hole wall with time

    當(dāng)柱狀藥卷(1)孔底起爆時(shí),爆轟波以半球面波的形式自藥卷底部向孔口方向傳播。爆轟波的波陣面在藥卷壁處以較小的傾斜角度β猛烈地撞擊藥卷周圍(2)區(qū)內(nèi)的空氣,并在接觸空氣的瞬間產(chǎn)生反射和透射。這個(gè)傾角β即為爆轟波對(duì)空氣介質(zhì)的入射角,由于入射角很小,可以近似地認(rèn)為爆轟波以正入射的方式接觸藥卷周圍的空氣。由于爆轟波后(3)區(qū)內(nèi)爆轟產(chǎn)物的沖擊阻抗遠(yuǎn)大于(2)區(qū)內(nèi)空氣介質(zhì)的沖擊阻抗,因此爆轟波在空氣介質(zhì)處的反射波為一膨脹波,反射的膨脹波可以理想化為經(jīng)由(3)區(qū)內(nèi)還未膨脹的爆轟產(chǎn)物向藥卷中心軸線上傳播的半徑不斷縮小的柱面波,受膨脹波的擾動(dòng),(3)區(qū)內(nèi)未膨脹的爆轟產(chǎn)物即開(kāi)始膨脹。膨脹波的波陣面最終匯聚于藥卷中心的軸線上,并且在軸線上產(chǎn)生一個(gè)在(4)區(qū)內(nèi)已經(jīng)膨脹的爆轟產(chǎn)物中向外傳播的新反射波。爆轟波透射形成的空氣沖擊波經(jīng)(2)區(qū)內(nèi)的空氣層向孔壁傳播,由于壓縮(5)區(qū)內(nèi)的空氣散失能量致使空氣沖擊波在空氣層中傳播時(shí),其峰值壓力和波陣面的速度會(huì)逐漸衰減。最后,衰減后的空氣沖擊波會(huì)以正入射的方式作用于孔壁巖石上,產(chǎn)生反射和透射。由于巖石的沖擊阻抗遠(yuǎn)大于(2)區(qū)內(nèi)空氣介質(zhì)的沖擊阻抗,因此空氣沖擊波在孔壁處的反射波仍為沖擊波,反射的沖擊波同樣可以理想化為經(jīng)由(4)區(qū)內(nèi)已經(jīng)膨脹的爆轟產(chǎn)物向炮孔中心軸線上傳播的半徑不斷減小的柱面波,并在其波陣面匯聚于炮孔中心的軸線上后向外傳播一個(gè)新的反射波。由于這個(gè)新的反射波在爆轟產(chǎn)物中傳播時(shí),氣態(tài)爆轟產(chǎn)物已經(jīng)迅速膨脹,充滿了整個(gè)炮孔,故可以認(rèn)為在爆轟氣體接觸孔壁之前,爆轟擾動(dòng)空氣產(chǎn)生的空氣沖擊波在孔壁的透射沖擊波參數(shù)即為孔壁初始沖擊參數(shù)。

    2 衰減規(guī)律

    爆轟擾動(dòng)空氣形成的初始空氣沖擊波以柱面波的形式向孔壁傳播,由于柱面半徑的增大,其波陣面單位面積上的能量將不斷減少,被壓縮后的空氣由于熵值的增大,會(huì)將部分動(dòng)能不可逆的轉(zhuǎn)化為熱能,這些能量的散失致使空氣沖擊波在到達(dá)孔壁時(shí),其峰值壓力和波陣面的速度已迅速下降,另外隨著沖擊波寬度的增加,對(duì)空氣的正壓作用時(shí)間也不斷加長(zhǎng)。

    空氣沖擊波在炮孔空氣層中衰減傳播的計(jì)算模型如圖2所示,其中Da、pa、μa、ρa(bǔ)分別表示空氣沖擊波在炮孔內(nèi)傳播距離為x時(shí)波陣面的速度、壓力、質(zhì)點(diǎn)移動(dòng)速度和密度;p0、ρ0、μ0分別表示空氣沖擊波前未受擾動(dòng)的空氣的壓力、密度和質(zhì)點(diǎn)移動(dòng)速度。未受擾動(dòng)的空氣處于靜止?fàn)顟B(tài),μ0=0,ρ0=1.2 kg/m3。

    圖2 空氣沖擊波在炮孔空氣層內(nèi)衰減傳播模型Fig.2 Attenuation propagation model of air shock wave in the air layer of blasthole

    假定當(dāng)空氣沖擊波在炮孔內(nèi)的傳播距離為x時(shí),途經(jīng)的空氣受其擾動(dòng),質(zhì)量全部集中于波陣面附近厚度為Δr的薄層之中。由于空氣沖擊波波陣面的壓力pa很強(qiáng),波前靜止空氣的壓力p0可以忽略,于是空氣沖擊波滿足強(qiáng)沖擊波的基本關(guān)系式[16]如下:

    (3)

    (4)

    (5)

    式中:k為絕熱指數(shù),一般取1.3。

    (6)

    (7)

    式中:pw為爆生氣體的平均爆轟壓力,pw=ρwD2/2(γ+1);Vw為裝藥體積,Vw=πr2l,l為炮孔長(zhǎng)度,裝藥系數(shù)為1;Vk為爆生氣體壓力為pk時(shí)的體積;px為爆生氣體膨脹到距藥卷中心x處時(shí)的平均壓力;Vx為爆生氣體壓力為px時(shí)的體積,Vx=πx2l。

    由界面條件可知,在薄層與爆生氣體的交界處,界面兩側(cè)的壓力相等,則

    (8)

    令pc等于波陣面壓強(qiáng)pa的θ倍,對(duì)薄層內(nèi)的空氣應(yīng)用牛頓第二定律,得

    (9)

    式中:m為被空氣沖擊波所席卷的空氣的質(zhì)量,m=π(x2-r2)lρ0。

    將式(3)和式(4)代入式(9)化簡(jiǎn)得

    (10)

    可以將上式寫(xiě)為

    (11)

    對(duì)兩邊積分得Da=A(x2-r2)θ-1,其中A為積分常數(shù)。

    (12)

    單位質(zhì)量的理想氣體,其內(nèi)能為

    e=pν/(k-1)

    (13)

    式中:p和ν分別表示氣體的壓力和比容。

    假設(shè)空氣薄層內(nèi)壓力呈線性變化,則薄層內(nèi)空氣的內(nèi)能ET為

    (14)

    將式(4)和式(5)代入式(14)并替換掉m和Da可得

    (15)

    爆炸傳遞給沖擊波的總能量EO為

    (16)

    在對(duì)爆轟擾動(dòng)空氣形成的初始沖擊波參數(shù)的計(jì)算過(guò)程中可以發(fā)現(xiàn),空氣沖擊波的傳播速度要大于爆轟產(chǎn)物的噴流速度[15],故可以認(rèn)為空氣沖擊波對(duì)空氣的裹挾要快于爆轟產(chǎn)物膨脹對(duì)空氣的壓縮,進(jìn)而認(rèn)為爆轟產(chǎn)物在膨脹到孔壁前為自由膨脹,對(duì)炮孔內(nèi)的空氣做功為零。爆轟傳遞給沖擊波的總能量EO為一定值,即EO為與x無(wú)關(guān)的常數(shù),故

    2θ-1=0

    (17)

    得θ=1/2。

    借助于式(8)求得的pc值,可以計(jì)算出當(dāng)空氣沖擊波傳播到孔壁時(shí),即當(dāng)x=R時(shí),空氣沖擊波對(duì)孔壁的入射壓力pa為

    (18)

    利用強(qiáng)沖擊波的基本關(guān)系式和pa值可得,當(dāng)空氣沖擊波傳播到孔壁時(shí),入射空氣沖擊波波陣面的速度Da和質(zhì)點(diǎn)移動(dòng)速度μa分別為

    (19)

    (20)

    將式(16)變換后代入Da=A(x2-r2)θ-1,并化簡(jiǎn)得

    (21)

    將式(21)代入式(3)和式(4)可得

    (22)

    (23)

    由式(21)~(23)得到空氣沖擊波參數(shù)隨傳播距離的衰減關(guān)系(見(jiàn)圖3~圖4)。

    圖3 空氣沖擊波波陣面速度和質(zhì)點(diǎn)移動(dòng)速度隨傳播距離的衰減Fig.3 Attenuation of wave front velocity and particle moving velocity of air shock wave with propagation distance

    圖4 空氣沖擊波波陣面壓力隨傳播距離的衰減Fig.4 Attenuation of wave front pressure of air shock wave with propagation distance

    從圖3和圖4可以看出,空氣沖擊波在炮孔內(nèi)空氣層中傳播時(shí),波陣面的速度、壓力和質(zhì)點(diǎn)移動(dòng)速度會(huì)隨傳播距離的增加呈指數(shù)型衰減,波陣面壓力的衰減速度要比波陣面速度和波陣面質(zhì)點(diǎn)移動(dòng)速度衰減的快。空氣沖擊波參數(shù)的大小與爆炸傳遞給沖擊波的總能量和藥卷長(zhǎng)度之間存在線性關(guān)系,爆炸傳遞給沖擊波的總能量越強(qiáng),則沖擊波的參數(shù)越大??諝鉀_擊波的衰減規(guī)律如下:

    (24)

    3 初始參數(shù)計(jì)算

    空氣沖擊波傳播到孔壁時(shí)即產(chǎn)生反射和透射,在求解透射波的初始參數(shù)時(shí),以往的研究中多將孔壁面作為彈性面處理,認(rèn)為透射波是直接以應(yīng)力波的形式在孔壁巖石中傳播的,進(jìn)而應(yīng)用彈性理論求解孔壁的初始沖擊壓力[16]。實(shí)際上,空氣沖擊波正入射孔壁之后,在已經(jīng)膨脹的爆轟產(chǎn)物中傳播的反射波和在孔壁巖石中傳播的透射波均為沖擊波,且反射沖擊波和透射沖擊波波陣面的壓力會(huì)因巖石波阻抗的增大而大于入射空氣沖擊波波陣面的壓力??諝鉀_擊波在孔壁處的壓力變化如圖5所示。

    注:pa,μa分別表示空氣沖擊波傳播到孔壁時(shí)波陣面的壓力和質(zhì)點(diǎn)速度;pt,μt分別表示在孔壁巖石中傳播的透射沖擊波波陣面的壓力和質(zhì)點(diǎn)速度;μs表示反射沖擊波波陣面的質(zhì)點(diǎn)速度,其與空氣沖擊波的傳播方向相反;pr,μr,ρr分別表示透射沖擊波前未受擾動(dòng)的巖石的壓力、質(zhì)點(diǎn)速度和密度。圖5 空氣沖擊波在孔壁處的透、反射Fig.5 The transmission and reflection of the air shock wave at the hole wall

    未受擾動(dòng)的巖石處于靜止?fàn)顟B(tài),μr=0。空氣與巖石交界面兩側(cè)質(zhì)點(diǎn)的移動(dòng)速度是相等的,因此在交界面上有

    μt=μa-μs

    (25)

    對(duì)反射沖擊波應(yīng)用沖擊波后質(zhì)點(diǎn)速度方程可得

    (26)

    式中:ca,ct分別為反射沖擊波波前和波后介質(zhì)的比容。

    反射沖擊波的沖擊絕熱線可以寫(xiě)為

    (27)

    將式(26)和式(27)代入式(25)可得

    (28)

    將式(18)和式(20)以及ca=(k-1)/(k+1)ρ0代入式(28)可得

    (29)

    對(duì)在巖石中傳播的透射沖擊波可用沖擊波后的質(zhì)點(diǎn)速度方程計(jì)算,由于pt?pr,方程中的pr可以忽略,得

    (30)

    式中:ρt為受透射沖擊波擾動(dòng)后巖石的密度。

    巖石的狀態(tài)方程為[17]

    (31)

    式中:Cr為未受擾動(dòng)的巖石中的縱波波速;B為常數(shù)。

    當(dāng)炸藥的密度ρw,爆速D,炮孔半徑R和藥卷半徑r已知時(shí),由式(18)和式(29)~式(31)即可求得孔壁初始沖擊壓力pt。

    4 結(jié)果對(duì)比與討論

    文獻(xiàn)[8]采用錳銅壓阻傳感器實(shí)驗(yàn)測(cè)量了水泥砂漿試件在空氣不耦合裝藥條件下的孔壁壓力,文獻(xiàn)[12]利用ANSYS/LS-DYNA軟件對(duì)文獻(xiàn)[8]中的實(shí)驗(yàn)?zāi)P瓦M(jìn)行了數(shù)值模擬。采用文獻(xiàn)[8]中的巖石和炸藥參數(shù)以及不耦合系數(shù),炸藥密度ρw=1.0 g/cm3,爆速D=3 200 m/s,水泥砂漿試件未受擾動(dòng)時(shí)的密度ρr=2.2 g/cm3,縱波波速Cr=3 200 m/s,壓力增大系數(shù)n取8,炮孔直徑40 mm,分別代入2個(gè)常用公式和本文推導(dǎo)的方程組中計(jì)算孔壁初始?jí)毫?,將?jì)算結(jié)果與文獻(xiàn)[8]的實(shí)測(cè)值及文獻(xiàn)[12]的模擬值共同繪制于圖6中。

    圖6 本文計(jì)算結(jié)果與模擬值、實(shí)測(cè)值的比較Fig.6 Comparison of calculated results with simulated and measured values

    從圖6可以看出,用本文推導(dǎo)方程組計(jì)算的孔壁初始?jí)毫Ρ瘸S霉降挠?jì)算值和模擬值更接近于實(shí)測(cè)值。當(dāng)不耦合系數(shù)為1.14時(shí),方程組計(jì)算的孔壁初始?jí)毫?,與實(shí)測(cè)值誤差為8%;當(dāng)不耦合系數(shù)為1.6時(shí),誤差為21%,均小于模擬值和常用公式計(jì)算值相對(duì)于實(shí)測(cè)值的誤差。方程組計(jì)算值與實(shí)測(cè)值隨不耦合系數(shù)變化的趨勢(shì)也較為一致。數(shù)值模擬的結(jié)果與實(shí)測(cè)值相比,可靠度為80%左右[12]。

    常用公式的計(jì)算結(jié)果比實(shí)測(cè)值小很多。常用公式依據(jù)沖擊波在炮孔壁的透、反射致使孔壁入射壓力增強(qiáng)的觀點(diǎn)得出,計(jì)算時(shí)將爆生氣體的準(zhǔn)靜態(tài)壓力作為沖擊波對(duì)孔壁的入射壓力,理論上并不嚴(yán)密。其中n的取值直接影響到孔壁壓力的計(jì)算值。n的取值與裝藥結(jié)構(gòu)、炸藥特性和巖體介質(zhì)等多種因素有關(guān),且應(yīng)當(dāng)隨不耦合系數(shù)的變化而變化,因此想要比較合理的確定n值是比較困難的[4]。

    如果在方程組的推導(dǎo)過(guò)程中將爆生氣體在炮孔內(nèi)的膨脹按照一階段的等熵膨脹計(jì)算,則此時(shí)用方程組計(jì)算的孔壁入射壓力值即為爆轟波波陣面的壓力,可得到耦合裝藥時(shí)孔壁壓力的計(jì)算方程組,其與文獻(xiàn)[11]中推導(dǎo)的耦合裝藥時(shí)巖體中爆炸沖擊波參數(shù)的計(jì)算公式一致。

    文獻(xiàn)[18]利用彈性波動(dòng)理論得出了耦合裝藥時(shí)透射波初始?jí)毫Φ挠?jì)算公式為

    (32)

    將文獻(xiàn)[8]中的巖石與炸藥參數(shù)代入式(32)計(jì)算得耦合裝藥時(shí)的孔壁壓力為3.52 GPa,與實(shí)驗(yàn)值和模擬值相差較大。方程組計(jì)算的耦合裝藥時(shí)的孔壁壓力為8.14 GPa,與實(shí)驗(yàn)值誤差僅為4.5%,遠(yuǎn)小于模擬值和式(32)計(jì)算值較實(shí)驗(yàn)值的誤差。

    筆者在求解薄層內(nèi)側(cè)壓力時(shí),假設(shè)爆生氣體內(nèi)壓力均勻分布,這樣會(huì)導(dǎo)致當(dāng)不耦合系數(shù)較小時(shí),計(jì)算的空氣沖擊波對(duì)孔壁的入射壓力值偏大,從而使孔壁壓力計(jì)算值偏大。其原因是當(dāng)不耦合系數(shù)較小時(shí),炮孔內(nèi)孔隙較小,爆生氣體的膨脹空間有限,將爆生氣體內(nèi)壓力做均壓處理會(huì)使計(jì)算的入射沖擊波壓力增大,導(dǎo)致計(jì)算的孔壁壓力值偏大。

    5 結(jié)論

    1)空氣不耦合裝藥時(shí),爆轟擾動(dòng)形成的空氣沖擊波在炮孔內(nèi)的空氣層中傳播時(shí),波陣面的速度、壓力和質(zhì)點(diǎn)移動(dòng)速度會(huì)隨距離的增加呈指數(shù)型衰減。

    2)空氣不耦合裝藥時(shí)孔壁初始沖擊壓力會(huì)隨裝藥不耦合系數(shù)的增大先迅速減小后趨于穩(wěn)定值。

    3)與經(jīng)驗(yàn)公式計(jì)算值和模擬值相比,本文推導(dǎo)方程組的計(jì)算值較實(shí)驗(yàn)值誤差更小。如果將爆生氣體的膨脹作一階段等熵膨脹考慮時(shí),該方程組也能較好的計(jì)算耦合裝藥時(shí)的孔壁壓力。

    猜你喜歡
    孔壁炮孔裝藥
    基于智能巖性識(shí)別的炮孔裝藥量計(jì)算*
    爆破(2022年4期)2022-12-17 02:45:26
    火炸藥學(xué)報(bào)(2022年3期)2022-07-04 07:31:00
    孔內(nèi)爆炸填塞效應(yīng)的數(shù)值模擬與分析
    基于礦用鉆孔成像儀的孔壁模擬實(shí)驗(yàn)裝置研制
    砂土地區(qū)鉆孔樁孔壁穩(wěn)定性與評(píng)價(jià)方法研究
    中外公路(2020年4期)2020-09-14 08:24:00
    某發(fā)射裝藥結(jié)構(gòu)改進(jìn)設(shè)計(jì)的新思路
    溶(土)洞地區(qū)沖孔樁孔壁聲波探測(cè)施工技術(shù)研究
    深孔不耦合裝藥爆破技術(shù)卸壓效果驗(yàn)證
    3Dmine 在雅滿蘇井下礦扇形中深孔爆破炮孔設(shè)計(jì)中的應(yīng)用
    新疆鋼鐵(2015年3期)2015-11-08 01:59:57
    一種凸多邊形的爆破網(wǎng)絡(luò)等時(shí)線算法研究
    另类精品久久| 久久久精品免费免费高清| av天堂久久9| 大话2 男鬼变身卡| av线在线观看网站| 亚洲国产最新在线播放| 一本色道久久久久久精品综合| 亚洲av电影在线进入| 亚洲精品国产一区二区精华液| 日本av手机在线免费观看| 欧美激情极品国产一区二区三区| √禁漫天堂资源中文www| 亚洲少妇的诱惑av| 成人影院久久| 街头女战士在线观看网站| 999久久久国产精品视频| 毛片一级片免费看久久久久| 在线观看三级黄色| 欧美日韩成人在线一区二区| 亚洲av日韩精品久久久久久密 | 十八禁人妻一区二区| 国产欧美亚洲国产| 精品视频人人做人人爽| 如何舔出高潮| 曰老女人黄片| 国产免费视频播放在线视频| 国产精品一区二区在线不卡| 亚洲国产av新网站| 国产女主播在线喷水免费视频网站| 尾随美女入室| 久久国产精品男人的天堂亚洲| 97精品久久久久久久久久精品| 9191精品国产免费久久| 中国国产av一级| 天堂8中文在线网| 黑人欧美特级aaaaaa片| 国产欧美日韩综合在线一区二区| 中文字幕另类日韩欧美亚洲嫩草| 在线观看免费视频网站a站| 老司机靠b影院| 欧美日韩av久久| 亚洲婷婷狠狠爱综合网| 国产高清国产精品国产三级| 在线观看免费午夜福利视频| 中文天堂在线官网| 国产欧美亚洲国产| 欧美黑人精品巨大| 免费黄网站久久成人精品| 老司机亚洲免费影院| 久久人人爽av亚洲精品天堂| 亚洲美女视频黄频| 黄频高清免费视频| 国产成人一区二区在线| 亚洲伊人色综图| 侵犯人妻中文字幕一二三四区| 中文字幕最新亚洲高清| 久久久久精品久久久久真实原创| 母亲3免费完整高清在线观看| 19禁男女啪啪无遮挡网站| 国产日韩欧美在线精品| 国产日韩欧美视频二区| 国语对白做爰xxxⅹ性视频网站| 激情五月婷婷亚洲| 精品一区二区免费观看| 交换朋友夫妻互换小说| 免费在线观看完整版高清| 亚洲国产毛片av蜜桃av| 一区二区三区乱码不卡18| 少妇人妻精品综合一区二区| 欧美 亚洲 国产 日韩一| 91老司机精品| 女人精品久久久久毛片| a级毛片在线看网站| 国产一区二区激情短视频 | 永久免费av网站大全| 国产探花极品一区二区| 一个人免费看片子| 建设人人有责人人尽责人人享有的| 一区二区三区四区激情视频| 大片电影免费在线观看免费| 99精品久久久久人妻精品| 中文字幕人妻丝袜制服| 国产精品二区激情视频| 亚洲国产看品久久| 十八禁高潮呻吟视频| 亚洲三区欧美一区| 色婷婷久久久亚洲欧美| 高清av免费在线| 最黄视频免费看| 欧美日韩综合久久久久久| 久久精品熟女亚洲av麻豆精品| a 毛片基地| 午夜日韩欧美国产| 色吧在线观看| 亚洲精品日韩在线中文字幕| 岛国毛片在线播放| 午夜影院在线不卡| 丝袜在线中文字幕| 丝袜脚勾引网站| 亚洲四区av| 国产成人a∨麻豆精品| 国产精品一区二区在线不卡| 亚洲国产av影院在线观看| 纵有疾风起免费观看全集完整版| 日本黄色日本黄色录像| 国产成人午夜福利电影在线观看| 亚洲图色成人| 黄色视频在线播放观看不卡| 亚洲四区av| 精品卡一卡二卡四卡免费| 久久久久人妻精品一区果冻| 一级片'在线观看视频| 午夜老司机福利片| 大香蕉久久成人网| 51午夜福利影视在线观看| 欧美激情极品国产一区二区三区| 亚洲成av片中文字幕在线观看| 丝袜喷水一区| 国产爽快片一区二区三区| 狂野欧美激情性bbbbbb| 亚洲四区av| 最近手机中文字幕大全| 欧美成人午夜精品| 精品少妇一区二区三区视频日本电影 | 久久亚洲国产成人精品v| 新久久久久国产一级毛片| 纯流量卡能插随身wifi吗| 久久久久国产精品人妻一区二区| 一本大道久久a久久精品| 99九九在线精品视频| 午夜免费鲁丝| 我的亚洲天堂| 丰满乱子伦码专区| 高清av免费在线| 一本—道久久a久久精品蜜桃钙片| 日韩制服骚丝袜av| 又大又黄又爽视频免费| 国产精品一国产av| 国产成人精品久久二区二区91 | 男女之事视频高清在线观看 | 纯流量卡能插随身wifi吗| 女的被弄到高潮叫床怎么办| 热re99久久精品国产66热6| 日韩 亚洲 欧美在线| 亚洲国产精品成人久久小说| 一二三四在线观看免费中文在| 国产 精品1| 99热全是精品| 丝袜在线中文字幕| 国产欧美日韩一区二区三区在线| 国产在线一区二区三区精| 十八禁高潮呻吟视频| 国产精品蜜桃在线观看| 国产亚洲av高清不卡| 亚洲成色77777| 精品酒店卫生间| 人人妻人人澡人人爽人人夜夜| 美女高潮到喷水免费观看| 岛国毛片在线播放| 亚洲精品美女久久久久99蜜臀 | 亚洲五月色婷婷综合| 最近2019中文字幕mv第一页| 欧美日韩精品网址| 韩国高清视频一区二区三区| 老司机靠b影院| 亚洲av在线观看美女高潮| www.自偷自拍.com| 少妇精品久久久久久久| 无遮挡黄片免费观看| 日韩一区二区视频免费看| 国产一卡二卡三卡精品 | 国产女主播在线喷水免费视频网站| 日韩电影二区| 亚洲成人免费av在线播放| 一区二区三区激情视频| 黄网站色视频无遮挡免费观看| 麻豆精品久久久久久蜜桃| 晚上一个人看的免费电影| 黄网站色视频无遮挡免费观看| 国产无遮挡羞羞视频在线观看| 另类精品久久| 亚洲国产欧美日韩在线播放| 999精品在线视频| 欧美 日韩 精品 国产| 999久久久国产精品视频| 国产欧美日韩一区二区三区在线| 嫩草影视91久久| 国产日韩欧美亚洲二区| 午夜影院在线不卡| 久久婷婷青草| 精品视频人人做人人爽| 国产成人午夜福利电影在线观看| 欧美人与善性xxx| 丝袜美腿诱惑在线| 欧美xxⅹ黑人| 亚洲精品aⅴ在线观看| 电影成人av| 国产精品久久久久久人妻精品电影 | 久久影院123| 青春草亚洲视频在线观看| 午夜免费鲁丝| 国产成人免费无遮挡视频| 99九九在线精品视频| 国产精品久久久久成人av| 国产精品免费视频内射| 欧美成人午夜精品| 三上悠亚av全集在线观看| 国产精品蜜桃在线观看| 在线亚洲精品国产二区图片欧美| 嫩草影视91久久| 久久久久久久国产电影| 十八禁高潮呻吟视频| 两个人免费观看高清视频| 精品免费久久久久久久清纯 | 啦啦啦在线免费观看视频4| 久久久精品94久久精品| 欧美人与性动交α欧美精品济南到| 一二三四在线观看免费中文在| 久久亚洲国产成人精品v| 七月丁香在线播放| 国产欧美日韩一区二区三区在线| 亚洲婷婷狠狠爱综合网| 在线观看免费高清a一片| 久久久久精品性色| 欧美日韩亚洲高清精品| 国产野战对白在线观看| h视频一区二区三区| 国产精品.久久久| 国产精品 国内视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品乱久久久久久| 国产片内射在线| 亚洲国产欧美网| 亚洲精品第二区| 亚洲av成人不卡在线观看播放网 | 欧美激情高清一区二区三区 | 美女国产高潮福利片在线看| 少妇 在线观看| 久久久久久久国产电影| 免费少妇av软件| 人妻一区二区av| 久久女婷五月综合色啪小说| 欧美精品人与动牲交sv欧美| 99精品久久久久人妻精品| 一本色道久久久久久精品综合| 自拍欧美九色日韩亚洲蝌蚪91| 中国国产av一级| 国产av国产精品国产| 国产日韩欧美视频二区| 999久久久国产精品视频| tube8黄色片| 亚洲欧美成人精品一区二区| 亚洲精品美女久久久久99蜜臀 | 精品国产一区二区三区四区第35| 亚洲国产精品999| 国产免费现黄频在线看| 国产精品 国内视频| 国产熟女午夜一区二区三区| 中国三级夫妇交换| 一级毛片电影观看| 黄片播放在线免费| 我要看黄色一级片免费的| 一二三四在线观看免费中文在| 免费黄网站久久成人精品| 另类亚洲欧美激情| 90打野战视频偷拍视频| 超色免费av| 女人精品久久久久毛片| 精品第一国产精品| 2021少妇久久久久久久久久久| 一级黄片播放器| 国产亚洲精品第一综合不卡| av视频免费观看在线观看| 女人久久www免费人成看片| 亚洲欧美精品自产自拍| 成人手机av| 欧美日韩亚洲国产一区二区在线观看 | 久久久精品国产亚洲av高清涩受| 免费久久久久久久精品成人欧美视频| 91精品国产国语对白视频| 色婷婷久久久亚洲欧美| 青春草国产在线视频| 久久精品国产a三级三级三级| 久久精品aⅴ一区二区三区四区| 日韩熟女老妇一区二区性免费视频| 最近最新中文字幕大全免费视频 | 丝袜美腿诱惑在线| 色精品久久人妻99蜜桃| 日本av免费视频播放| 王馨瑶露胸无遮挡在线观看| 午夜日韩欧美国产| 久久久久网色| 亚洲国产欧美日韩在线播放| 午夜免费观看性视频| 精品国产露脸久久av麻豆| 免费高清在线观看视频在线观看| 亚洲一级一片aⅴ在线观看| 国产精品女同一区二区软件| 精品久久久精品久久久| 天堂俺去俺来也www色官网| 亚洲av在线观看美女高潮| 在线观看免费视频网站a站| 久久久国产精品麻豆| 亚洲激情五月婷婷啪啪| 亚洲av福利一区| 国产成人a∨麻豆精品| 精品国产超薄肉色丝袜足j| 18禁国产床啪视频网站| 中文天堂在线官网| 国产成人91sexporn| 欧美中文综合在线视频| 日日爽夜夜爽网站| 亚洲欧美激情在线| 国产精品香港三级国产av潘金莲 | 男人添女人高潮全过程视频| 国产一区二区三区综合在线观看| 肉色欧美久久久久久久蜜桃| 国产精品秋霞免费鲁丝片| 久久午夜综合久久蜜桃| 下体分泌物呈黄色| 亚洲国产毛片av蜜桃av| 国产男女超爽视频在线观看| 男人操女人黄网站| 久久ye,这里只有精品| 国产精品久久久久久精品电影小说| 久久久久久久国产电影| 久久亚洲国产成人精品v| 国产成人精品无人区| 亚洲成国产人片在线观看| 毛片一级片免费看久久久久| 飞空精品影院首页| 精品国产一区二区三区四区第35| 成人午夜精彩视频在线观看| 国产成人精品福利久久| 国产在线免费精品| 午夜av观看不卡| 在线 av 中文字幕| 性色av一级| 精品视频人人做人人爽| 亚洲国产精品一区二区三区在线| 亚洲免费av在线视频| 99九九在线精品视频| 精品亚洲乱码少妇综合久久| 在线天堂中文资源库| 亚洲色图 男人天堂 中文字幕| 91老司机精品| 韩国高清视频一区二区三区| 成年动漫av网址| 一区二区三区激情视频| 国产成人欧美| 国产不卡av网站在线观看| 亚洲国产欧美一区二区综合| 色综合欧美亚洲国产小说| 欧美国产精品va在线观看不卡| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕另类日韩欧美亚洲嫩草| 国产麻豆69| 久久精品国产综合久久久| 亚洲精品在线美女| 99国产精品免费福利视频| av在线观看视频网站免费| 亚洲欧洲精品一区二区精品久久久 | 国产精品免费大片| www日本在线高清视频| 嫩草影视91久久| 国产在线视频一区二区| 99热国产这里只有精品6| 中文字幕av电影在线播放| 久久久久久人人人人人| 免费高清在线观看日韩| 国产国语露脸激情在线看| 啦啦啦在线免费观看视频4| 久久精品久久久久久久性| 不卡视频在线观看欧美| bbb黄色大片| 日韩欧美精品免费久久| 亚洲第一av免费看| 久久天躁狠狠躁夜夜2o2o | 一级片'在线观看视频| 水蜜桃什么品种好| 黄片小视频在线播放| 欧美日韩福利视频一区二区| 黄色毛片三级朝国网站| videos熟女内射| 欧美xxⅹ黑人| 亚洲欧洲日产国产| 另类精品久久| 女人久久www免费人成看片| 操美女的视频在线观看| 国产精品熟女久久久久浪| 久久婷婷青草| 久久精品亚洲av国产电影网| 我的亚洲天堂| 亚洲伊人色综图| 国产av精品麻豆| 啦啦啦在线观看免费高清www| 97精品久久久久久久久久精品| 免费在线观看完整版高清| 免费观看a级毛片全部| 丝袜美足系列| 涩涩av久久男人的天堂| 亚洲av欧美aⅴ国产| 欧美久久黑人一区二区| 久久久久久久国产电影| 国产免费现黄频在线看| 最近最新中文字幕免费大全7| 国产精品熟女久久久久浪| 人妻一区二区av| 69精品国产乱码久久久| 久久 成人 亚洲| 丁香六月天网| 日韩人妻精品一区2区三区| 亚洲国产精品一区三区| 女人被躁到高潮嗷嗷叫费观| 人人妻人人澡人人爽人人夜夜| av片东京热男人的天堂| av女优亚洲男人天堂| 亚洲人成网站在线观看播放| 亚洲色图 男人天堂 中文字幕| 美女福利国产在线| 午夜影院在线不卡| 亚洲一区二区三区欧美精品| 啦啦啦视频在线资源免费观看| 精品亚洲成国产av| 老鸭窝网址在线观看| 欧美日韩视频高清一区二区三区二| 欧美日韩亚洲高清精品| 免费观看人在逋| 夫妻午夜视频| 飞空精品影院首页| 成人亚洲欧美一区二区av| 熟女少妇亚洲综合色aaa.| 韩国高清视频一区二区三区| 亚洲精品美女久久久久99蜜臀 | 中文字幕人妻熟女乱码| av天堂久久9| 在线亚洲精品国产二区图片欧美| 国产男女超爽视频在线观看| 国产一区亚洲一区在线观看| 波野结衣二区三区在线| 男女边摸边吃奶| 91精品国产国语对白视频| 少妇精品久久久久久久| 亚洲激情五月婷婷啪啪| 久久青草综合色| 天天影视国产精品| 香蕉国产在线看| 亚洲欧美清纯卡通| 免费在线观看完整版高清| 少妇人妻久久综合中文| 精品少妇一区二区三区视频日本电影 | 嫩草影视91久久| 亚洲七黄色美女视频| 欧美久久黑人一区二区| 国产日韩欧美亚洲二区| 色精品久久人妻99蜜桃| av国产精品久久久久影院| 日本91视频免费播放| 香蕉国产在线看| 婷婷色综合www| 亚洲,欧美精品.| 日日撸夜夜添| 亚洲熟女毛片儿| 国产精品av久久久久免费| 成人18禁高潮啪啪吃奶动态图| 人人妻,人人澡人人爽秒播 | 亚洲精品国产av蜜桃| 天天操日日干夜夜撸| 久久久久精品国产欧美久久久 | 亚洲精品aⅴ在线观看| 丁香六月天网| 成年av动漫网址| 97精品久久久久久久久久精品| 免费观看性生交大片5| 麻豆乱淫一区二区| 亚洲中文av在线| 欧美日韩亚洲国产一区二区在线观看 | 丰满饥渴人妻一区二区三| 国产有黄有色有爽视频| 综合色丁香网| 久久狼人影院| 日韩一卡2卡3卡4卡2021年| 水蜜桃什么品种好| 国产精品一区二区在线观看99| 国产片内射在线| 无遮挡黄片免费观看| av免费观看日本| 狂野欧美激情性bbbbbb| 欧美日韩亚洲综合一区二区三区_| 国产成人av激情在线播放| 成年动漫av网址| 99热网站在线观看| 国产无遮挡羞羞视频在线观看| 51午夜福利影视在线观看| 啦啦啦在线免费观看视频4| 赤兔流量卡办理| 一本—道久久a久久精品蜜桃钙片| 啦啦啦啦在线视频资源| 欧美在线一区亚洲| 午夜日韩欧美国产| 亚洲,欧美精品.| 日本猛色少妇xxxxx猛交久久| 国产乱人偷精品视频| 免费观看a级毛片全部| 一二三四在线观看免费中文在| 精品国产乱码久久久久久男人| 国产av国产精品国产| 一区二区日韩欧美中文字幕| e午夜精品久久久久久久| 女的被弄到高潮叫床怎么办| 亚洲av国产av综合av卡| 亚洲四区av| 大陆偷拍与自拍| 久久99热这里只频精品6学生| 欧美变态另类bdsm刘玥| 国产男人的电影天堂91| 人人妻,人人澡人人爽秒播 | 国产又色又爽无遮挡免| 中国三级夫妇交换| 自线自在国产av| 一区二区三区精品91| 黄色视频在线播放观看不卡| 99精国产麻豆久久婷婷| 日韩成人av中文字幕在线观看| 国产精品国产av在线观看| 一区二区av电影网| 久久人人爽人人片av| 亚洲av日韩在线播放| av国产久精品久网站免费入址| 国产在线视频一区二区| 国产精品三级大全| 少妇精品久久久久久久| 欧美日韩一级在线毛片| 青春草国产在线视频| 老鸭窝网址在线观看| 婷婷色综合www| 波野结衣二区三区在线| 亚洲 欧美一区二区三区| 久久婷婷青草| 精品久久蜜臀av无| 亚洲av在线观看美女高潮| 欧美日韩一级在线毛片| 黄网站色视频无遮挡免费观看| 天天操日日干夜夜撸| 精品久久久精品久久久| 黄色一级大片看看| 观看美女的网站| 欧美激情 高清一区二区三区| 亚洲国产中文字幕在线视频| 乱人伦中国视频| 日韩精品有码人妻一区| 飞空精品影院首页| 中文欧美无线码| 日韩一区二区三区影片| 国产精品 国内视频| 色婷婷久久久亚洲欧美| 免费观看av网站的网址| 丰满少妇做爰视频| 看十八女毛片水多多多| 久久97久久精品| 高清视频免费观看一区二区| 国产人伦9x9x在线观看| 精品人妻在线不人妻| 久久影院123| 啦啦啦 在线观看视频| 国产成人a∨麻豆精品| 欧美激情极品国产一区二区三区| 黄色怎么调成土黄色| 高清视频免费观看一区二区| 久久久久人妻精品一区果冻| 国产不卡av网站在线观看| 免费av中文字幕在线| 午夜免费男女啪啪视频观看| 亚洲欧洲日产国产| 午夜福利乱码中文字幕| 午夜激情久久久久久久| 日韩视频在线欧美| 777米奇影视久久| 国产淫语在线视频| 国产亚洲午夜精品一区二区久久| 精品视频人人做人人爽| 一级片'在线观看视频| 青春草视频在线免费观看| 高清av免费在线| 两个人免费观看高清视频| 日本猛色少妇xxxxx猛交久久| 久久狼人影院| 免费av中文字幕在线| 国产午夜精品一二区理论片| 激情五月婷婷亚洲| 久热这里只有精品99| 亚洲精品国产色婷婷电影| 另类亚洲欧美激情| 亚洲国产精品一区三区| 国产毛片在线视频| 久久青草综合色| 黄色视频不卡| 色播在线永久视频| 波野结衣二区三区在线| 美国免费a级毛片| 肉色欧美久久久久久久蜜桃| 一级毛片我不卡| 精品国产乱码久久久久久男人| 亚洲四区av| 纵有疾风起免费观看全集完整版| 97人妻天天添夜夜摸| 国产黄色免费在线视频| 国产亚洲一区二区精品| 亚洲精品视频女| 国产精品一二三区在线看| 日本午夜av视频| 高清在线视频一区二区三区| 亚洲av福利一区| 国产精品国产av在线观看| 国产精品女同一区二区软件|