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

    雙空孔間距對(duì)爆破槽腔斷面大小的影響*

    2022-12-02 10:12:02李祥龍張志平王建國(guó)王子琛
    爆炸與沖擊 2022年11期
    關(guān)鍵詞:空孔炮孔裝藥

    李祥龍,張志平,王建國(guó),李 強(qiáng),王子琛

    (1. 昆明理工大學(xué)國(guó)土資源工程學(xué)院,云南 昆明 650093;2. 昆明理工大學(xué)云南省中-德藍(lán)色礦山與特殊地下空間開(kāi)發(fā)利用重點(diǎn)實(shí)驗(yàn)室,云南 昆明 650093;3. 中煤科工集團(tuán)西安研究院有限公司,陜西 西安 710077)

    巷道掘進(jìn)爆破中,一般只有一個(gè)自由面,所受夾制作用較大,通常會(huì)留有空孔為巖石破碎提供補(bǔ)償空間[1]。而空孔數(shù)目、間距等參數(shù)是提高掘進(jìn)進(jìn)尺和掏槽取得成功的關(guān)鍵。目前一致認(rèn)可的是補(bǔ)償空間越大,掏槽效果越好,即增加空孔數(shù)量或空孔直徑,然而大的補(bǔ)償空間需要犧牲的是鉆孔成本。那么,如何讓有限的空孔實(shí)現(xiàn)最佳的爆破空間補(bǔ)償效果,是掏槽爆破研究的熱點(diǎn)問(wèn)題。

    就掏槽爆破參數(shù)優(yōu)化問(wèn)題,主要從數(shù)值模擬、現(xiàn)場(chǎng)實(shí)驗(yàn)等方面進(jìn)行了研究。Li 等[2]利用數(shù)值模擬手段驗(yàn)證了中空孔應(yīng)力集中效應(yīng)的存在;Zhang 等[3]提出了掏槽爆破的一種先進(jìn)的掏槽拋片(CCFT)概念;Sui 等[4]研究了無(wú)限巖體在有自由面條件下裂紋擴(kuò)展過(guò)程;單仁亮等[5]模擬研究了不同空孔孔徑的直眼掏槽爆破效果,得出孔徑對(duì)掏槽爆破的效果和振動(dòng)均有重要影響;汪海波等[6]認(rèn)為增大空孔直徑可提高應(yīng)力峰值和延長(zhǎng)作用時(shí)間,空孔周?chē)逯禐闊o(wú)空孔時(shí)的2.2 倍;Li 等[7]利用數(shù)值模擬手段研究得出兩空孔間距越短,止裂作用越強(qiáng),當(dāng)兩孔間距小于一定值時(shí),可完全阻止外向裂紋;鐘波波等[8]利用數(shù)值模擬方法,分析得出空孔有導(dǎo)向作用,且與空孔大小無(wú)明顯相關(guān);材料的非均勻性,對(duì)裂紋的擴(kuò)展方式有顯著影響;張召冉等[9]針對(duì)含空孔直眼掏槽參數(shù)確定沒(méi)有可靠依據(jù)的問(wèn)題,提出了動(dòng)態(tài)碎脹系數(shù)的概念;Gao 等[10]采用ANSYS/LS-DYNA 軟件對(duì)不同側(cè)壓力系數(shù)的復(fù)雜五孔掏槽爆破進(jìn)行二維優(yōu)化模擬研究,結(jié)果表明,設(shè)置延遲爆轟和孔徑以及優(yōu)化炮眼間距可以優(yōu)化破巖效果;柴修偉等[11]利用數(shù)值模擬軟件LS-DYNA 模擬了不同數(shù)量的空孔情況下裂紋擴(kuò)展情況。

    鑒于以上研究成果,本文以在研的礦山巷道掏槽爆破掘進(jìn)效率優(yōu)化項(xiàng)目為依托,選擇雙大直徑空孔(孔徑100 mm)、側(cè)向兩孔起爆的掏槽方案,研究同等起爆條件下大直徑空孔距離對(duì)掏槽槽腔斷面大小的影響,以期確定最優(yōu)的掏槽布孔方式,保證礦山巷道掘進(jìn)進(jìn)尺的同時(shí),減少空孔鉆鑿數(shù)量。

    1 掏槽爆破的空孔效應(yīng)

    1.1 空孔的自由面作用

    起爆后,沖擊波隨著傳播的距離不斷衰減,應(yīng)力也不斷下降,逐漸變?yōu)閴嚎s波,當(dāng)壓縮波傳到空孔孔壁會(huì)反射形成拉伸波,與應(yīng)力波傳到自由面的作用相同,因此空孔充當(dāng)了自由面的作用,有利于巖石破碎。

    根據(jù)動(dòng)量守恒定律可推出[12]:

    式中:σ 為縱波應(yīng)力,MPa;τ 為橫波應(yīng)力,MPa;ρ 為巖石密度,kg/m3;cp為縱波波速,m/s;cs為橫波波速,m/s;vp和vs分別為質(zhì)點(diǎn)在波的傳播方向和垂直與波的傳播方向上的運(yùn)動(dòng)速度,m/s。

    1.2 空孔的應(yīng)力集中作用

    力學(xué)計(jì)算和實(shí)驗(yàn)研究表明:當(dāng)應(yīng)力波傳到空孔孔壁時(shí),將產(chǎn)生反射拉應(yīng)力波并在空孔附近形成拉應(yīng)力區(qū),最大拉應(yīng)力峰值將出現(xiàn)在的兩孔連心線方向上[13-14]。因此,空孔具有應(yīng)力集中作用[6]。在某測(cè)點(diǎn)的徑向壓應(yīng)力峰值為:

    式中:σrmax為某點(diǎn)的爆炸徑向應(yīng)力,MPa;p0為孔壁初始應(yīng)力峰值,MPa;r0為炮孔半徑,m;r為某點(diǎn)到炮孔距離,m;a為應(yīng)力波衰減指數(shù)。

    切向拉應(yīng)力峰值可通過(guò)徑向壓應(yīng)力峰值得到:

    式中:b為系數(shù),在爆炸近區(qū)b≈1,b值與應(yīng)力波傳播距離成反比,并逐漸趨于穩(wěn)定。

    式中:μ為巖石泊松比。

    耦合裝藥時(shí),孔壁初始應(yīng)力峰值為[15-16]:

    式中:ρe、ρm分別為炸藥和巖石密度,kg/m3;D0為爆速,m/s;cp為縱波波速,m/s。

    計(jì)算可得到某一點(diǎn)切向拉應(yīng)力峰值σθmax,若切向拉應(yīng)力峰值σθmax大于巖石的動(dòng)態(tài)抗拉強(qiáng)度σθ,即滿足巖石拉裂破壞條件,巖石破碎。

    1.3 空孔的卸壓作用

    在復(fù)雜地壓下,根據(jù)厚壁圓筒理論,爆炸與原巖應(yīng)力場(chǎng)耦合為:

    空孔使炮孔在徑向的抵抗線要遠(yuǎn)小于軸向的抵抗線,因此,破碎巖體在爆生氣體的推動(dòng)作用下,優(yōu)先向空孔方向移動(dòng),空孔提供的空間可以消掉一定的碎脹壓力,有利于巖石拋擲出腔,即空孔的卸壓作用。

    由上文可知,大直徑空孔直眼掏槽機(jī)理是應(yīng)力集中、自由面、空孔的卸壓作用共同作用的結(jié)果。當(dāng)爆炸氣體產(chǎn)生時(shí),空孔起到自由面作用,為破碎巖體提供一定的碎脹空間,有利于破碎巖體在爆生氣體的推動(dòng)作用下拋出,形成槽腔,即實(shí)現(xiàn)空孔卸壓作用。下文利用數(shù)值模擬手段研究雙空孔間距對(duì)爆破效果的影響。

    2 數(shù)值分析模型的建立

    2.1 模型尺寸及試驗(yàn)方案

    (1) 補(bǔ)償空間理論

    補(bǔ)償空間應(yīng)滿足礦巖碎脹后所需要的體積,不然就會(huì)發(fā)生擠死的現(xiàn)象,從而影響爆破效果。其應(yīng)滿足如下關(guān)系式:

    式中:S1為預(yù)爆巖體面積;S2為空孔面積;S3為裝藥孔面積;K為巖石碎脹系數(shù),K=1.5。

    根據(jù)圖1,可推導(dǎo)出裝藥孔與空孔距離的關(guān)系:

    圖1 裝藥孔與空孔的距離關(guān)系Fig. 1 Distance between charge hole and empty hole

    式中:L為空孔與裝藥孔的距離;d為裝藥孔直徑;D為空孔直徑。計(jì)算得到L<32.7 cm。

    (2) 炮孔偏斜理論

    為防止炮孔間貫穿,需要滿足:

    式中:H為炮孔深度,a為炮孔偏斜度。

    現(xiàn)場(chǎng)使用KJ311 型全液壓掘進(jìn)鉆車(chē),炮孔偏斜角度在1°以內(nèi),經(jīng)計(jì)算得:L>18.9 cm。所以裝藥孔與空孔之間的距離范圍為:18.9 cm<L<32.7 cm。

    運(yùn)用ANSYS/LSDYNA 建立有限元二維模型,模型尺寸為4000 mm×4000 mm,邊緣定義為無(wú)反射邊界,裝藥炮孔直徑為50 mm,大直徑空孔直徑為100 mm,裝藥方式為耦合裝藥。沿模型中心橫向設(shè)置兩個(gè)裝藥炮孔,同時(shí)起爆,孔距為40 cm,兩空孔按模型中心豎向布置,空孔孔距用dv表示,通過(guò)上式計(jì)算裝藥孔與空孔之間孔距18.9 cm<L<32.7 cm。根據(jù)勾股定理原則可計(jì)算出dv最大為51.7 cm,本文將空孔孔距dv設(shè)置為15、25、35、45 和55 cm,定義兩空孔連心線的中點(diǎn)為P,各炮孔位置及模型網(wǎng)格劃分如圖2 所示。

    圖2 炮孔布置及二維模型網(wǎng)格劃分圖Fig. 2 Layout of blast holes and meshing of the two-dimensional model

    2.2 巖石模型

    結(jié)合相關(guān)文獻(xiàn),本巖石模型選用更加適合巖體的爆炸沖擊數(shù)值模擬的HJC 本構(gòu)模型[17-18]。根據(jù)大紅山銅礦試驗(yàn)區(qū)域參數(shù)[19],確定HJC 本構(gòu)模型中的礦巖參數(shù)見(jiàn)表1。

    表1 大紅山大理巖HJC 本構(gòu)模型參數(shù)Table 1 Dahongshan marble HJC constitutive model parameters

    2.3 炸藥模型參數(shù)及狀態(tài)方程

    礦山作業(yè)使用的是1 號(hào)巖石乳化炸藥,選用*MAT_HIGH_EXPLOSIVE_BURN 模型,采用EOS_JWL狀態(tài)方程模擬炸藥狀態(tài)[20-21]

    式中:p1為爆轟產(chǎn)物內(nèi)部壓力,V為相對(duì)體積,A1、B1、R1、R2和ω 為炸藥相關(guān)常數(shù),E0為初始比內(nèi)能。

    本文采用的炸藥模型參數(shù)見(jiàn)表2。

    表2 炸藥的狀態(tài)方程參數(shù)Table 2 HJC constitutive model parameters of explosive

    3 數(shù)值計(jì)算結(jié)果及分析

    3.1 爆破應(yīng)力場(chǎng)分布情況

    為了研究雙空孔對(duì)直眼掏槽爆炸應(yīng)力場(chǎng)的影響,模擬試驗(yàn)設(shè)置兩個(gè)垂直掏槽孔裝藥起爆,炮孔直徑50 mm,空孔直徑100 mm,僅改變雙空孔孔距,得到dv=15 cm、dv=25 cm、dv=35 cm、dv=45 cm 和dv=55 cm五種情況下不同時(shí)刻的等效壓力云圖,見(jiàn)圖3~圖7,P 點(diǎn)處應(yīng)力時(shí)程曲線如圖8 所示。

    圖3 當(dāng)dv=15 cm 時(shí)模型壓力云圖Fig. 3 Pressure contours of the model when dv=15 cm

    圖4 當(dāng)dv=25 cm 時(shí)模型壓力云圖Fig. 4 Pressure contours of the model when dv=25 cm

    圖5 當(dāng)dv=35 cm 時(shí)模型壓力云圖Fig. 5 Pressure contours of the model when dv=35 cm

    圖6 當(dāng)dv=45 cm 時(shí)模型壓力云圖Fig. 6 Pressure contours of the model when dv=45 cm

    圖7 當(dāng)dv=55 cm 時(shí)模型壓力云圖Fig. 7 Pressure contours of the model when dv=55 cm

    圖8 P 點(diǎn)處的應(yīng)力時(shí)程曲線Fig. 8 Stress time history curve at point P

    分析壓力云圖可知:在應(yīng)力波到達(dá)空孔前,五種模型壓力云圖的應(yīng)力峰值基本一致;在應(yīng)力波傳播至空孔附近時(shí),可明顯觀察出兩空孔連線方向出現(xiàn)應(yīng)力集中現(xiàn)象,應(yīng)力波傳播到中心P點(diǎn)處時(shí)兩炮孔產(chǎn)生的應(yīng)力波相互疊加,應(yīng)力波在空孔壁形成反射拉伸波,隨后出現(xiàn)不同的應(yīng)力波疊加擾動(dòng),直至傳播出視野。

    通過(guò)觀察不同dv模型P點(diǎn)處的應(yīng)力時(shí)程曲線可以看出:受到的壓應(yīng)力峰值均為217 MPa,dv=15 cm、dv=25 cm、dv=35 cm、dv=45 cm 、dv=55 cm 時(shí)P點(diǎn)處拉應(yīng)力峰值分別為197、180、165、160 和156 MPa;即dv由15 cm 增大到25 cm,拉應(yīng)力峰值減小8.62%;dv由25 cm 增大到35 cm,拉應(yīng)力峰值減小8.3%;dv由35 cm 增大到45 cm,拉應(yīng)力峰值減小3.03%;dv由45 cm 增大到55 cm,拉應(yīng)力峰值減小0.25%。說(shuō)明dv越小,P點(diǎn)受到的反射拉應(yīng)力越大,這跟空孔自由面的距離有關(guān)。

    3.2 巖石損傷破壞情況

    只考慮雙大直徑空孔間距對(duì)爆破效果的影響時(shí),空孔間距越小,P點(diǎn)處受到的拉應(yīng)力峰值越大,中心巖體越容易破壞。為方便研究槽腔斷面破壞情況,在k 文件中添加*MAT_ADD_EROSION 關(guān)鍵字[11],設(shè)定壓力強(qiáng)度大于本文實(shí)際巖石力學(xué)參數(shù),即單軸抗壓、單軸抗拉、變形模量等時(shí),模擬巖體的損傷與破壞,如圖9~圖13 所示。

    如圖9 所示,dv=15 cm 時(shí),空孔距離太近,補(bǔ)償空間導(dǎo)向作用不明顯;當(dāng)dv=25 cm 時(shí)(圖10),在t=0.05 ms 時(shí),裝藥孔周?chē)鷰r石單元受壓破壞,應(yīng)力波剛剛達(dá)到空孔壁,在t=0.1 ms 時(shí)應(yīng)力波持續(xù)向外傳播,裝藥孔周?chē)鷰r石破壞向外擴(kuò)展,空孔壁處應(yīng)力波發(fā)生反射開(kāi)始受拉破壞,出現(xiàn)空孔導(dǎo)向現(xiàn)象,使空孔與以破壞單元連接形成初步槽腔輪廓,在t=0.15 ms 時(shí),兩裝藥孔與空孔之間完全貫穿,在t=0.2 ms 時(shí),槽腔輪廓面基本形成,直至0.25 ms 時(shí)刻破壞終止。

    圖9 dv=15 cm 時(shí)損傷破壞圖Fig. 9 Damage failure diagram when dv=15 cm

    圖10 dv=25 cm 時(shí)損傷破壞圖Fig. 10 Damage failure diagram when dv=25 cm

    當(dāng)dv=35 cm 時(shí)(圖11),炮孔周?chē)鸭y發(fā)展更加充分,同時(shí)反射應(yīng)力波對(duì)巖石的拉伸破壞減小,爆破后得到的槽腔輪廓面積較dv=25 cm 時(shí)有所增大,中心區(qū)域雖有少量巖體,但基本貫通形成槽腔輪廓面。

    圖11 dv=35 cm 時(shí)損傷破壞圖Fig. 11 Damage failure diagram when dv=35 cm

    當(dāng)dv=45 cm 時(shí)(圖12),由0.25 ms 時(shí)刻損傷破壞圖可以看出兩裝藥孔之間基本貫穿,但槽腔外斷面仍存在大量未失效單元,說(shuō)明掏槽區(qū)域仍有大量巖體存留,若應(yīng)用于現(xiàn)場(chǎng)能會(huì)出現(xiàn)補(bǔ)償空間不夠的情況,不能保證掏槽效果;當(dāng)dv=55 cm 時(shí)(圖13),空孔距離太大,爆生裂紋無(wú)法貫穿,無(wú)法形成有效槽腔。

    圖12 dv=45 cm 時(shí)損傷破壞圖Fig. 12 Damage failure diagram when dv=45 cm

    圖13 dv=55 cm 時(shí)損傷破壞圖Fig. 13 Damage failure diagram when dv=55 cm

    巖石破碎區(qū)域若存在寬度為1 單元的細(xì)長(zhǎng)區(qū)域,即將該區(qū)域判定為裂隙,統(tǒng)計(jì)破碎區(qū)域面積時(shí)不計(jì)入統(tǒng)計(jì)。對(duì)模擬結(jié)果得到的槽腔外斷面進(jìn)行提取,即CAD 等比例繪圖計(jì)算,得到圖14。根據(jù)圖14 可知,dv=15 cm 模型爆破后形成槽腔斷面寬度為63.03 cm,形成的槽腔斷面積為0.1641 m2;dv=25 cm 模型爆破后形成槽腔斷面寬度為68.6 cm,形成的槽腔斷面積為0.2116 m2;dv=35 cm 模型爆破后形成的槽腔斷面為67.4 cm,形成的槽腔斷面面積為0.2436 m2;dv=45 cm 模型爆破后槽腔斷面內(nèi)存在大量未失效巖石單元,面積為0.1740 m2,dv=55 cm 模型爆破后形成槽腔斷面寬度為62.27 cm,高度僅為16.65 cm,形成的槽腔斷面積為0.0951 m2。可見(jiàn)雙空孔方案中dv由15 cm 增加到25 cm 時(shí),形成的槽腔斷面增大了28.94%;dv由25 cm 增加到35 cm 時(shí),槽腔斷面增大15.1%;但當(dāng)dv由35 cm 增加到45 cm 時(shí),槽腔面積為0.1740 m2,形成的槽腔斷面減小了17.8%;dv增大到55 cm 時(shí),槽腔面積為0.0951 m2,形成的槽腔斷面減小了45.3%;即槽腔面積隨空孔間距的增大先增加后減小,dv=35 cm 時(shí)槽腔面積最大,確定為現(xiàn)場(chǎng)驗(yàn)證方案。

    圖14 炸后槽腔外斷面Fig. 14 The external section of the explosion cavity after explosion

    4 爆破現(xiàn)場(chǎng)驗(yàn)證

    4.1 炮孔布置

    試驗(yàn)地點(diǎn)為大紅山銅礦285 中段、370 水平,根據(jù)大紅山銅礦已有資料顯示,礦巖主要為大理巖,巖石力學(xué)參數(shù)見(jiàn)表1。裝藥孔直徑為50 mm,空孔直徑為100 mm,孔深為3.2 m。由于上述模擬得出dv=35 cm 時(shí)形成的槽腔斷面最大,故對(duì)dv=35 cm 方案進(jìn)行現(xiàn)場(chǎng)驗(yàn)證。掏槽區(qū)域炮孔布置方案如圖2(c)所示。炸藥為1 號(hào)巖石乳化炸藥,采用耦合裝藥,孔底起爆;所有炮孔均用炮泥填塞,炮泥長(zhǎng)度為0.6 m。兩個(gè)裝藥孔共需要1 號(hào)巖石乳化炸藥3.2 kg。

    4.2 爆破成腔效果及分析

    通過(guò)數(shù)值模擬計(jì)算分析,選擇空孔間距dv=35 cm,兩裝藥孔距離40 cm 的掏槽方案進(jìn)行現(xiàn)場(chǎng)試驗(yàn),具體試驗(yàn)結(jié)果見(jiàn)表3,由于篇幅原因本文只對(duì)2 號(hào)試驗(yàn)爆破后形成的槽腔進(jìn)行分析,如圖15 所示。由于未進(jìn)行全斷面開(kāi)挖爆破,實(shí)際上形成錐臺(tái)型槽腔,測(cè)量槽腔底部斷面寬度與高度并繪制槽腔斷面輪廓,以此計(jì)算槽腔面積,并與模擬結(jié)果對(duì)比,如圖16?,F(xiàn)場(chǎng)試驗(yàn)槽腔斷面寬度比模擬結(jié)果小4.0%,槽腔斷面高度比模擬結(jié)果小3.4%,槽腔斷面積比模擬結(jié)果小4.98%,多次現(xiàn)場(chǎng)試驗(yàn)與模擬結(jié)果誤差均在5%以內(nèi)。

    表3 現(xiàn)場(chǎng)試驗(yàn)結(jié)果Table 3 Field test results

    圖15 現(xiàn)場(chǎng)掏槽爆破槽腔斷面圖Fig. 15 Cross-sectional view of the cavity of the on-site cutting blasting cavity

    圖16 現(xiàn)場(chǎng)斷面與模擬結(jié)果對(duì)比Fig. 16 Comparison of the site section with the simulation results

    5 結(jié) 論

    以大紅山銅礦地下某巷道為試驗(yàn)場(chǎng)地,并根據(jù)圍巖力學(xué)性質(zhì)和實(shí)際爆破條件進(jìn)行數(shù)值模擬計(jì)算對(duì)比分析,得到如下結(jié)論:

    (1) 兩等間距掏槽孔同時(shí)起爆條件下,掏槽區(qū)域中心點(diǎn)的壓應(yīng)力峰值不變,反射拉應(yīng)力峰值隨著大直徑空孔間距的增加而減小;

    (2) 當(dāng)空孔間距dv為15 cm、25 cm、35 cm、45 cm、55 cm 時(shí),模擬所得槽腔面積分別為0.1641 m2、0.2116 m2、0.2436 m2、0.1740 m2、0.0951 m2,表明存在一個(gè)最佳的空孔間距;

    (3) 模擬預(yù)測(cè)的槽腔斷面寬度、高度、斷面積分別比現(xiàn)場(chǎng)試驗(yàn)結(jié)果大4.0%、3.4%和4.98%,多次現(xiàn)場(chǎng)試驗(yàn)與模擬結(jié)果誤差均在5%以內(nèi),說(shuō)本文預(yù)測(cè)爆破槽腔斷面積的數(shù)值模擬方法具有一定可靠性。

    因此,可進(jìn)一步將此數(shù)值模擬方法從二維拓展至三維,研究不同掏槽方式下的成腔體積并選擇性開(kāi)展現(xiàn)場(chǎng)試驗(yàn)驗(yàn)證,以期找到預(yù)測(cè)地下巷道掏槽爆破成腔效果的數(shù)值方法。

    猜你喜歡
    空孔炮孔裝藥
    不耦合裝藥下炮孔—空孔距離對(duì)預(yù)裂爆破效果的影響探究
    工程爆破(2023年5期)2023-11-02 12:59:44
    基于智能巖性識(shí)別的炮孔裝藥量計(jì)算*
    爆破(2022年4期)2022-12-17 02:45:26
    《混合炸藥及裝藥??沸蜓?/a>
    孔內(nèi)爆炸填塞效應(yīng)的數(shù)值模擬與分析
    空孔作用下定向漲裂破巖試驗(yàn)研究
    等體積空孔直眼掏槽槽腔形成過(guò)程及其分析*
    爆破(2020年4期)2020-12-16 09:09:52
    某發(fā)射裝藥結(jié)構(gòu)改進(jìn)設(shè)計(jì)的新思路
    深孔不耦合裝藥爆破技術(shù)卸壓效果驗(yàn)證
    3Dmine 在雅滿蘇井下礦扇形中深孔爆破炮孔設(shè)計(jì)中的應(yīng)用
    新疆鋼鐵(2015年3期)2015-11-08 01:59:57
    定向斷裂控制爆破的空孔效應(yīng)實(shí)驗(yàn)分析*
    爆炸與沖擊(2015年3期)2015-06-07 11:38:12
    日本猛色少妇xxxxx猛交久久| 麻豆成人av视频| 熟女人妻精品中文字幕| 久久午夜福利片| 日韩成人av中文字幕在线观看| 十分钟在线观看高清视频www | 亚洲av免费高清在线观看| 黄片无遮挡物在线观看| 国产免费又黄又爽又色| 美女高潮的动态| av女优亚洲男人天堂| 天天躁日日操中文字幕| 舔av片在线| 这个男人来自地球电影免费观看 | 欧美 日韩 精品 国产| 涩涩av久久男人的天堂| 欧美zozozo另类| 精品久久久久久久久av| 日韩精品有码人妻一区| 纵有疾风起免费观看全集完整版| 精品酒店卫生间| 麻豆成人av视频| 午夜免费观看性视频| 久久精品熟女亚洲av麻豆精品| 精品亚洲乱码少妇综合久久| 天堂俺去俺来也www色官网| 男人和女人高潮做爰伦理| 久久97久久精品| 亚洲欧美日韩卡通动漫| 亚洲av电影在线观看一区二区三区| 我的女老师完整版在线观看| 欧美亚洲 丝袜 人妻 在线| 精品熟女少妇av免费看| h日本视频在线播放| 欧美激情国产日韩精品一区| 日日啪夜夜爽| 欧美日韩亚洲高清精品| 亚洲四区av| 日日摸夜夜添夜夜爱| 搡老乐熟女国产| 国产男人的电影天堂91| 男人舔奶头视频| 日韩av不卡免费在线播放| 我的老师免费观看完整版| 日本-黄色视频高清免费观看| 伦理电影大哥的女人| 欧美 日韩 精品 国产| 国产人妻一区二区三区在| 中文字幕精品免费在线观看视频 | 精品少妇久久久久久888优播| 国产男人的电影天堂91| 欧美亚洲 丝袜 人妻 在线| 美女脱内裤让男人舔精品视频| 国产探花极品一区二区| 看免费成人av毛片| 18禁动态无遮挡网站| 国产精品熟女久久久久浪| av在线播放精品| 日韩三级伦理在线观看| 美女主播在线视频| 91aial.com中文字幕在线观看| 1000部很黄的大片| 久久久久久久久久人人人人人人| 久久av网站| 国产伦精品一区二区三区四那| 中文字幕久久专区| 一级爰片在线观看| 高清毛片免费看| 久久精品国产鲁丝片午夜精品| 国产v大片淫在线免费观看| 国产精品嫩草影院av在线观看| 丰满乱子伦码专区| 久久久久视频综合| 免费观看a级毛片全部| 日日啪夜夜撸| freevideosex欧美| 国产免费福利视频在线观看| 插逼视频在线观看| 岛国毛片在线播放| 久久久久久久久久成人| 国产色爽女视频免费观看| 国产一级毛片在线| 免费黄色在线免费观看| 一本久久精品| 亚洲精品国产成人久久av| 精品久久久久久久久av| 一级片'在线观看视频| 亚洲av福利一区| 成人国产av品久久久| 乱系列少妇在线播放| 大片免费播放器 马上看| 亚洲成人手机| 午夜老司机福利剧场| 亚洲四区av| 亚洲欧美日韩另类电影网站 | 99久久精品一区二区三区| 成人午夜精彩视频在线观看| 亚洲av中文字字幕乱码综合| 欧美日本视频| 伊人久久精品亚洲午夜| 国产在线视频一区二区| 九九在线视频观看精品| 永久网站在线| 男人舔奶头视频| av播播在线观看一区| 国产黄频视频在线观看| 一区在线观看完整版| 国产精品av视频在线免费观看| 视频区图区小说| 最近最新中文字幕大全电影3| 成人一区二区视频在线观看| 久久久欧美国产精品| 中文资源天堂在线| 欧美人与善性xxx| 黄色一级大片看看| av福利片在线观看| 成人毛片60女人毛片免费| 日韩人妻高清精品专区| 国产极品天堂在线| 又大又黄又爽视频免费| 久久久久久久久大av| 夫妻午夜视频| 久久青草综合色| 国产在视频线精品| 最近2019中文字幕mv第一页| 日本黄色片子视频| 久久久久久人妻| 夜夜爽夜夜爽视频| 欧美+日韩+精品| 免费av不卡在线播放| 精品熟女少妇av免费看| 中文欧美无线码| 免费不卡的大黄色大毛片视频在线观看| xxx大片免费视频| 啦啦啦中文免费视频观看日本| 亚洲国产精品成人久久小说| 极品教师在线视频| 国产成人精品一,二区| 18+在线观看网站| 三级经典国产精品| 国内少妇人妻偷人精品xxx网站| 国产午夜精品久久久久久一区二区三区| 少妇被粗大猛烈的视频| 自拍欧美九色日韩亚洲蝌蚪91 | 国产国拍精品亚洲av在线观看| 一级毛片久久久久久久久女| 欧美zozozo另类| 国产成人aa在线观看| 欧美极品一区二区三区四区| 高清在线视频一区二区三区| 乱码一卡2卡4卡精品| 男人添女人高潮全过程视频| 成人高潮视频无遮挡免费网站| 亚洲第一区二区三区不卡| 午夜福利在线在线| 免费播放大片免费观看视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 免费久久久久久久精品成人欧美视频 | 午夜福利视频精品| 国内精品宾馆在线| 狂野欧美激情性bbbbbb| 男人狂女人下面高潮的视频| 97在线人人人人妻| 少妇的逼水好多| 亚洲av二区三区四区| 国产永久视频网站| av在线观看视频网站免费| 一级毛片 在线播放| 一级毛片电影观看| 一区二区三区免费毛片| 国产精品99久久99久久久不卡 | 91精品一卡2卡3卡4卡| 夫妻性生交免费视频一级片| 国产爱豆传媒在线观看| 一区二区av电影网| av一本久久久久| 天天躁日日操中文字幕| 韩国av在线不卡| av在线老鸭窝| 国产成人午夜福利电影在线观看| 欧美最新免费一区二区三区| 久久99蜜桃精品久久| 国产精品av视频在线免费观看| 久久ye,这里只有精品| 久久6这里有精品| 国内精品宾馆在线| 午夜精品国产一区二区电影| 国产白丝娇喘喷水9色精品| 我的老师免费观看完整版| 国产一区二区在线观看日韩| 一边亲一边摸免费视频| 卡戴珊不雅视频在线播放| 夜夜骑夜夜射夜夜干| 超碰97精品在线观看| 欧美一级a爱片免费观看看| 亚洲性久久影院| 高清午夜精品一区二区三区| 免费高清在线观看视频在线观看| 亚洲欧美成人综合另类久久久| 一级毛片 在线播放| 男女边吃奶边做爰视频| 免费观看无遮挡的男女| 一本一本综合久久| 久久精品国产鲁丝片午夜精品| 国产伦理片在线播放av一区| 在线天堂最新版资源| 国产一区二区三区av在线| 免费久久久久久久精品成人欧美视频 | 不卡视频在线观看欧美| 成年美女黄网站色视频大全免费 | 欧美高清性xxxxhd video| 久久这里有精品视频免费| 亚洲av二区三区四区| 91aial.com中文字幕在线观看| 国产成人免费观看mmmm| 美女xxoo啪啪120秒动态图| 国产在视频线精品| 久久久久久久大尺度免费视频| 80岁老熟妇乱子伦牲交| 国产av精品麻豆| 久久久久久九九精品二区国产| 精品久久久久久久久亚洲| 亚洲国产欧美在线一区| 色视频在线一区二区三区| 99久国产av精品国产电影| 久久久久久久大尺度免费视频| 80岁老熟妇乱子伦牲交| 人妻夜夜爽99麻豆av| 最新中文字幕久久久久| 深爱激情五月婷婷| 亚洲国产精品成人久久小说| 久久精品国产鲁丝片午夜精品| 内射极品少妇av片p| 亚洲精品久久久久久婷婷小说| 我的女老师完整版在线观看| 免费大片18禁| 麻豆成人午夜福利视频| 熟妇人妻不卡中文字幕| 亚洲欧美成人精品一区二区| 欧美区成人在线视频| 男女下面进入的视频免费午夜| 天堂8中文在线网| 色婷婷久久久亚洲欧美| 内地一区二区视频在线| 欧美+日韩+精品| 久久av网站| 久久人妻熟女aⅴ| 午夜福利在线观看免费完整高清在| 两个人的视频大全免费| 丰满人妻一区二区三区视频av| 国产黄片视频在线免费观看| 亚洲av不卡在线观看| 国产日韩欧美在线精品| 国产精品久久久久久久电影| 亚洲激情五月婷婷啪啪| 97超视频在线观看视频| 成人影院久久| 中文字幕精品免费在线观看视频 | 国产精品一二三区在线看| 大又大粗又爽又黄少妇毛片口| 国产黄片美女视频| 免费人妻精品一区二区三区视频| 又粗又硬又长又爽又黄的视频| 嫩草影院入口| 午夜精品国产一区二区电影| 国产伦理片在线播放av一区| 91久久精品电影网| 啦啦啦在线观看免费高清www| 日韩欧美精品免费久久| 国产91av在线免费观看| 国产毛片在线视频| 国产精品久久久久久久电影| 日韩中字成人| 成人特级av手机在线观看| av不卡在线播放| 欧美高清性xxxxhd video| 大香蕉久久网| 自拍偷自拍亚洲精品老妇| 免费久久久久久久精品成人欧美视频 | 免费黄频网站在线观看国产| 麻豆成人av视频| 人体艺术视频欧美日本| 欧美亚洲 丝袜 人妻 在线| 国产精品免费大片| 国模一区二区三区四区视频| 亚洲国产色片| 亚洲人成网站高清观看| 伊人久久精品亚洲午夜| 中文字幕av成人在线电影| 亚洲精品乱久久久久久| av卡一久久| 国产精品蜜桃在线观看| 亚洲在久久综合| 22中文网久久字幕| 各种免费的搞黄视频| 91久久精品国产一区二区成人| 久久午夜福利片| 国产亚洲一区二区精品| 日本欧美视频一区| 欧美xxⅹ黑人| 日本午夜av视频| 在线观看国产h片| 51国产日韩欧美| 老师上课跳d突然被开到最大视频| 久久av网站| 男人舔奶头视频| 各种免费的搞黄视频| 亚洲,一卡二卡三卡| 最近的中文字幕免费完整| tube8黄色片| 夜夜看夜夜爽夜夜摸| 亚洲综合色惰| 91在线精品国自产拍蜜月| 国产精品国产三级国产专区5o| 联通29元200g的流量卡| 日本-黄色视频高清免费观看| 日韩伦理黄色片| 国产高清有码在线观看视频| 韩国av在线不卡| 国产精品国产av在线观看| 国产乱人视频| 九色成人免费人妻av| 在线 av 中文字幕| 亚洲国产精品国产精品| www.av在线官网国产| 在线观看一区二区三区| 亚洲欧美日韩卡通动漫| 美女视频免费永久观看网站| 日韩成人伦理影院| 亚洲精品久久久久久婷婷小说| 久久人人爽人人爽人人片va| 亚洲欧美中文字幕日韩二区| 国产日韩欧美在线精品| 国产高清不卡午夜福利| 亚洲精品乱久久久久久| 国模一区二区三区四区视频| 久久99热6这里只有精品| 成年女人在线观看亚洲视频| 亚洲国产高清在线一区二区三| 亚洲国产欧美在线一区| 下体分泌物呈黄色| 男人添女人高潮全过程视频| 亚洲成人av在线免费| 3wmmmm亚洲av在线观看| 国产黄片视频在线免费观看| 欧美日韩亚洲高清精品| 99国产精品免费福利视频| 国产真实伦视频高清在线观看| 视频中文字幕在线观看| 久久久久性生活片| 久久婷婷青草| 永久免费av网站大全| 亚洲内射少妇av| 亚洲欧美中文字幕日韩二区| 少妇精品久久久久久久| 国产精品一区二区在线观看99| 精品国产三级普通话版| 精品少妇久久久久久888优播| 三级国产精品片| 日韩在线高清观看一区二区三区| 这个男人来自地球电影免费观看 | 国产乱来视频区| 国产视频首页在线观看| 纯流量卡能插随身wifi吗| 国产亚洲av片在线观看秒播厂| 男女无遮挡免费网站观看| 伊人久久精品亚洲午夜| 黄片无遮挡物在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品乱码久久久v下载方式| 欧美区成人在线视频| 自拍欧美九色日韩亚洲蝌蚪91 | 久久韩国三级中文字幕| 中文字幕亚洲精品专区| 色综合色国产| 色吧在线观看| 国产老妇伦熟女老妇高清| 中国美白少妇内射xxxbb| 精品国产乱码久久久久久小说| 在线观看免费高清a一片| 精品人妻一区二区三区麻豆| 看免费成人av毛片| 日韩制服骚丝袜av| 日韩欧美精品免费久久| av在线app专区| 成人毛片60女人毛片免费| 少妇的逼好多水| 亚洲欧美中文字幕日韩二区| 亚洲欧美成人精品一区二区| 国产一区二区三区av在线| 综合色丁香网| 亚洲精品一区蜜桃| 三级国产精品欧美在线观看| 欧美少妇被猛烈插入视频| av卡一久久| 国产淫语在线视频| 婷婷色麻豆天堂久久| 国产黄片视频在线免费观看| 免费黄网站久久成人精品| 肉色欧美久久久久久久蜜桃| 精品久久久久久久久av| 免费看av在线观看网站| 久久久精品免费免费高清| 日韩人妻高清精品专区| 伦理电影免费视频| 免费黄网站久久成人精品| 老司机影院毛片| 国产一区亚洲一区在线观看| 亚洲婷婷狠狠爱综合网| 亚洲第一av免费看| 国产精品偷伦视频观看了| 91精品伊人久久大香线蕉| 日本欧美视频一区| 免费人妻精品一区二区三区视频| 观看美女的网站| 观看免费一级毛片| 免费人妻精品一区二区三区视频| 久久国产亚洲av麻豆专区| 青春草视频在线免费观看| 自拍偷自拍亚洲精品老妇| 久久婷婷青草| 两个人的视频大全免费| 国产黄片视频在线免费观看| 熟妇人妻不卡中文字幕| 亚洲欧美精品自产自拍| 成人国产av品久久久| 国产熟女欧美一区二区| 精品久久久久久久久av| 久久久久久久国产电影| 亚洲精品自拍成人| 能在线免费看毛片的网站| 国产精品国产av在线观看| 欧美bdsm另类| 亚洲精品日本国产第一区| 国产在线视频一区二区| 一本色道久久久久久精品综合| 国产中年淑女户外野战色| 亚洲色图av天堂| 激情五月婷婷亚洲| 国产乱人偷精品视频| 久久久亚洲精品成人影院| 亚洲aⅴ乱码一区二区在线播放| 在线观看免费视频网站a站| 成人黄色视频免费在线看| 五月天丁香电影| 欧美激情极品国产一区二区三区 | 激情五月婷婷亚洲| 超碰av人人做人人爽久久| 精品亚洲乱码少妇综合久久| 人妻夜夜爽99麻豆av| 久久久久久久精品精品| 一个人免费看片子| 成人黄色视频免费在线看| 亚洲国产最新在线播放| kizo精华| 十分钟在线观看高清视频www | 色婷婷av一区二区三区视频| 久久人人爽人人片av| 一个人看的www免费观看视频| 欧美国产精品一级二级三级 | 亚洲国产成人一精品久久久| 寂寞人妻少妇视频99o| 一区二区三区四区激情视频| 日韩电影二区| 日韩伦理黄色片| 少妇人妻精品综合一区二区| 免费av不卡在线播放| 我要看日韩黄色一级片| 亚洲图色成人| 午夜福利在线在线| 日韩欧美精品免费久久| 亚洲,欧美,日韩| 色哟哟·www| 新久久久久国产一级毛片| 午夜日本视频在线| 夜夜骑夜夜射夜夜干| 免费大片黄手机在线观看| 日韩av免费高清视频| 久久久久久久国产电影| 国产乱人视频| 日本一二三区视频观看| 亚洲国产日韩一区二区| 国产精品人妻久久久影院| 久久久久久人妻| 免费高清在线观看视频在线观看| 国语对白做爰xxxⅹ性视频网站| 内射极品少妇av片p| 日本猛色少妇xxxxx猛交久久| 国产在线一区二区三区精| 国产极品天堂在线| 一本一本综合久久| 日韩,欧美,国产一区二区三区| 高清av免费在线| 久久久久久久久久久丰满| 国产黄频视频在线观看| 久久久色成人| 国产精品一区二区在线观看99| 亚洲精品国产成人久久av| 国产91av在线免费观看| 久久精品人妻少妇| 日韩欧美精品免费久久| 亚洲在久久综合| 丰满人妻一区二区三区视频av| 成人二区视频| 亚洲精品视频女| 高清欧美精品videossex| av福利片在线观看| 日本av手机在线免费观看| 精品一区二区三卡| 色婷婷av一区二区三区视频| 日本黄色片子视频| 高清黄色对白视频在线免费看 | 好男人视频免费观看在线| 伦精品一区二区三区| 久久ye,这里只有精品| 青春草视频在线免费观看| 亚洲av男天堂| 纯流量卡能插随身wifi吗| 一区二区三区精品91| 色婷婷久久久亚洲欧美| 成人亚洲精品一区在线观看 | 狂野欧美白嫩少妇大欣赏| 欧美xxⅹ黑人| 欧美xxxx黑人xx丫x性爽| 女性生殖器流出的白浆| 国产精品精品国产色婷婷| 日韩一区二区视频免费看| 日韩中字成人| 日本黄色片子视频| 国产乱人视频| 精品一区二区三区视频在线| 夫妻性生交免费视频一级片| 亚洲,一卡二卡三卡| 国产成人freesex在线| 天美传媒精品一区二区| 欧美高清成人免费视频www| 国产精品三级大全| 日韩欧美一区视频在线观看 | 亚洲av福利一区| 99热这里只有是精品50| 日韩三级伦理在线观看| 精品人妻一区二区三区麻豆| 亚洲怡红院男人天堂| 特大巨黑吊av在线直播| 青春草国产在线视频| 午夜免费观看性视频| 只有这里有精品99| a 毛片基地| 午夜福利高清视频| 午夜老司机福利剧场| 精品国产露脸久久av麻豆| 97热精品久久久久久| 久久久久久久久久人人人人人人| 欧美另类一区| 国产在线男女| a级毛片免费高清观看在线播放| 小蜜桃在线观看免费完整版高清| 国产亚洲午夜精品一区二区久久| 黄色欧美视频在线观看| 777米奇影视久久| 男女无遮挡免费网站观看| 午夜激情久久久久久久| 久久国产精品男人的天堂亚洲 | 亚洲av福利一区| 国产成人精品一,二区| 色婷婷av一区二区三区视频| 亚洲美女搞黄在线观看| 免费播放大片免费观看视频在线观看| 在线亚洲精品国产二区图片欧美 | 免费少妇av软件| 大码成人一级视频| 日本欧美视频一区| 国产国拍精品亚洲av在线观看| 日本欧美视频一区| 国产精品三级大全| 另类亚洲欧美激情| 91久久精品国产一区二区三区| 国产色婷婷99| 国产精品三级大全| 欧美老熟妇乱子伦牲交| 观看美女的网站| www.av在线官网国产| 久久韩国三级中文字幕| 网址你懂的国产日韩在线| 女人十人毛片免费观看3o分钟| av福利片在线观看| 日本一二三区视频观看| 欧美精品一区二区大全| 国产精品久久久久久精品电影小说 | 亚洲av中文av极速乱| 国产精品一区www在线观看| 日韩av不卡免费在线播放| 亚洲精品成人av观看孕妇| 自拍偷自拍亚洲精品老妇| 人妻系列 视频| 国产 一区精品| 中文字幕免费在线视频6| 菩萨蛮人人尽说江南好唐韦庄| 精品亚洲乱码少妇综合久久| 成人国产av品久久久| 噜噜噜噜噜久久久久久91| 亚洲无线观看免费| 18禁在线播放成人免费| 国产在线视频一区二区| 亚洲,一卡二卡三卡| 精品酒店卫生间| 午夜老司机福利剧场| 日本免费在线观看一区| 久久久久网色| 亚洲精品亚洲一区二区| 国产视频内射| 成人黄色视频免费在线看| 欧美 日韩 精品 国产| 男女下面进入的视频免费午夜| 干丝袜人妻中文字幕|