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

    承插式混凝土管道爆破振動動力響應尺寸效應研究

    2021-12-16 20:21黃一文蔣楠周傳波羅學東姚穎康夏宇磬
    振動工程學報 2021年5期

    黃一文 蔣楠 周傳波 羅學東 姚穎康 夏宇磬

    摘要: 為研究爆破振動作用下不同內徑規(guī)格承插式混凝土管道振動響應存在的尺寸效應,通過現(xiàn)場鄰近管道爆破試驗及其振動監(jiān)測數(shù)據(jù)驗證,結合動力有限元軟件建立多尺寸管道爆破動力數(shù)值計算模型,分析承插式混凝土管道爆破振動響應特征;基于量綱分析基本理論,考慮管道尺寸效應影響建立管道爆破振動速度預測模型;結合混凝土的動抗拉強度,提出了不同內徑混凝土管道的爆破振動速度控制標準。研究結果表明:管身和管道承插口峰值應力的出現(xiàn)存在時間差,管身峰值應力大于管道承插口,并隨著管道內徑的增加,兩者峰值應力差值逐漸降低;在爆破地震波作用下,混凝土管道峰值振動速度隨著管內徑的增加而逐漸降低。

    關鍵詞: 爆破振動; 管道尺寸; 動力響應; 承插式混凝土管道; 安全判據(jù)

    中圖分類號: U455.6 文獻標志碼: A 文章編號: 1004-4523(2021)05-0969-10

    DOI:10.16385/j.cnki.issn.1004-4523.2021.05.010

    引 言

    城區(qū)復雜環(huán)境下埋地管道工程交錯縱橫、環(huán)境復雜多變。承插式混凝土管道具有強度高、抗?jié)B性好、抗外壓能力強、生產(chǎn)效率高等優(yōu)點,被廣泛應用于城市給排水工程中。同時,近年來隨著城市交通迅速發(fā)展,地鐵等地下工程建設大量興起,這些工程常與既有地下管道并行、交叉、跨越。爆破作為地下工程建設堅硬巖體重要的開挖手段,施工產(chǎn)生的地震波不可避免地對鄰近混凝土給排水管道的安全性和可用性產(chǎn)生嚴重影響。保證爆破施工鄰近承插式混凝土管道的安全運行,研究爆破振動作用下管道響應特征及其爆破振動安全判據(jù)具有重要的理論研究及實際應用價值。

    目前,對于爆破振動荷載對鄰近管道的影響,眾多學者已開展了相關研究工作[1?3]。其中,在管道爆破振動響應影響因素研究方面,劉建民等[4]利用ANSYNS/LS?DYNA軟件分析了炸藥設置方式、炸藥用量、爆心距和管道管徑等因素對埋地管道損傷的不同作用。鐘冬望等[5]在理論分析的基礎上,通過對埋地無縫鋼管進行現(xiàn)場爆破實驗得出了應變峰值與爆心距和藥量計算公式。在管道爆破振動安全判據(jù)研究方面,郝郁清等[6]結合管道最大容許應變計算出在已知裝藥量的條件下爆破施工時的最小安全距離。張紫劍等[7]通過現(xiàn)場監(jiān)測試驗,并對峰值振速進行擬合分析,得出不同頻率下的爆破振動安全允許峰值振速。張震等[8]采用現(xiàn)場監(jiān)測和動力有限元數(shù)值模擬相結合的研究方法,對超淺埋地鐵站通道爆破開挖附近埋地混凝土管道的動力響應進行研究,提出保證管道安全的地表爆破控制振速。高壇等[9]通過數(shù)值模型計算和理論分析得到相應的管道爆破振動速度安全閾值。張黎明等[10]以管道的最大容許軸向應變峰值為管道的破壞判據(jù),得到埋地鋼管在不同主頻下的安全允許峰值振速。綜合分析已有研究成果可知,混凝土管道的爆破響應數(shù)值模擬研究中往往忽略管道承插口和管道尺寸的影響。事實上,城市市政中采用的承插式混凝土排水管尺寸規(guī)格在主線和支線上大有區(qū)別,在管道尺寸效應影響下管道爆破振動響應往往存在較大差異性。

    為此,本文以武漢市城區(qū)典型土層埋置的給排水管道系統(tǒng)為工程依托,通過全尺度預埋承插式混凝土管道現(xiàn)場爆破試驗,結合ANSYS/LS?DYNA有限元動力軟件建立爆破振動作用下埋地混凝土管道的數(shù)值計算模型,分析了承插式混凝土管道在爆破地震波作用下的應力變化規(guī)律?;诹烤V分析,研究了在管道尺寸效應影響下的管道上方地表土體和管道振速的變化規(guī)律,并依據(jù)混凝土構件的動抗拉強度來確定不同管道直徑的承插式混凝土管道的振動速度安全判據(jù)。

    1 鄰近管道現(xiàn)場爆破試驗

    1.1 試驗場地條件及混凝土管埋設

    為研究承插式混凝土管道的爆破振動動力響應,對武漢市現(xiàn)有的排水管道附近的巖土爆破開挖工程進行大量調研,根據(jù)統(tǒng)計數(shù)據(jù)可知,武漢市采用的排水管道多為混凝土管道,管徑主要在700?1600 mm之間,排水管道多埋置在粉質黏土以及雜填土中。本次試驗選取的試驗場地位于武漢市經(jīng)濟開發(fā)區(qū)硃山路與硃山一路交匯處,爆破所在地地質情況主要為強風化石英砂巖。管道覆土及墊土均為粉質黏土,現(xiàn)場管道埋設及炮孔位置如圖1所示。

    現(xiàn)場試驗使用的承插式混凝土管道符合混凝土和鋼筋混凝土管道國家標準要求[11],管道內徑1 m,外徑1.2 m;采用柔性接頭B型承插口管接頭,接頭處設置止水橡膠圈;共使用4節(jié)管道,每節(jié)管道長2.5 m。為模擬城市承插式混凝土管道埋設條件,現(xiàn)場使用機械挖掘出管道溝槽,溝槽深3.6 m,墊土高0.4 m,上覆黏土厚度為2 m。組裝承插式混凝土管道后,按照給水排水管道工程施工及驗收標準規(guī)范要求分層夯實墊土層至地表標高,并確保槽底至管道頂部以上500 mm范圍內,土中不含大于50 mm的磚、石等硬塊[12]。在試驗過程中,承插式混凝土管道不做加壓處理,只考慮承插式混凝土管道在無水狀態(tài)下爆破地震波作用下的動力響應。

    1.2 管道振動速度監(jiān)測和現(xiàn)場爆破方案

    埋地管道在爆破地震波作用下的動力響應,主要表現(xiàn)為管道的振動特征和應力應變特征變化。本試驗將混凝土管道的振動速度作為主要的監(jiān)測和研究對象。

    根據(jù)研究需要,本試驗中采用爆破振動速度監(jiān)測儀TC?4850,承插式混凝土管道在溝槽內安裝完畢后,將振動速度傳感器沿混凝土管道底部軸向布設,并進行固定,確保傳感器在試驗過程中保持穩(wěn)定,共設置6個監(jiān)測點:D1,D2,D3,D4,D5,D6,如圖2所示,并通過信號輸入抗干擾接頭與爆破振動速度監(jiān)測儀TC?4850相連。

    如圖1所示,在距混凝土管10 m處采用機械進行垂直鉆孔,炮孔直徑9 cm,炮孔深度為6 m,現(xiàn)場采用2#巖石乳化炸藥進行爆破試驗,炸藥總質量為8 kg,采用單孔集中裝藥,并用炮泥進行堵塞,待混凝土管道安裝埋設、監(jiān)測裝置安裝完畢后,采用非電導爆管進行單孔單段起爆。

    2 數(shù)值計算建模及其驗證

    2.1 模型介紹

    由于現(xiàn)場試驗條件有限,未進行不同尺寸承插式混凝土管道的現(xiàn)場爆破試驗,采用ANSYS/LS?DYNA有限元動力軟件建立承插式混凝土管道爆破振動數(shù)值計算模型,對混凝土管道的爆破振動動力響應尺寸效應進行研究。

    根據(jù)現(xiàn)場試驗,建立爆源與混凝土管道水平距離為10 m的數(shù)值計算模型。模型尺寸為21.2 m(長)10 m(寬)8 m(高),通過對幾種不同網(wǎng)格劃分方案進行比選,為保證數(shù)值計算的正常進行以及計算的精度,將土體及巖體的網(wǎng)格尺寸控制在20 cm,將混凝土管道網(wǎng)格按照每環(huán)40個進行劃分,將炸藥塊體等分成10個單元。為模擬承插式混凝土管道的實際工作情況,按照實際尺寸建立承插式混凝土管,分段建立4節(jié)管段,然后設置接觸,并在混凝土管道接觸間設置橡膠,4節(jié)管道總長10 m,爆心距為10 m,炮孔深6 m,裝藥長度為1.06 m,炮泥長度為4.94 m,如圖3所示。為模擬實際情況,管道與土層之間采用面面接觸,設置摩擦系數(shù)為0.15;其他參數(shù)為默認值,在模型頂面定義為自由面,其余面則定義為無反射邊界。

    2.2 材料模型和算法介紹

    模型中的材料包括巖層、土層、混凝土管道等,具體參數(shù)如表1所示。巖體介質是非連續(xù)、不均勻的,為簡化處理,本文中的數(shù)值模擬將巖土體假設為均勻介質的彈塑性材料。根據(jù)現(xiàn)場地質資料,該地區(qū)巖層以強風化石英砂巖為主,采用*MAT_PLASTIC_KINEMATIC材料模型[13]。試驗場地溝槽內覆土及墊土為粉質黏土,采用*MAT_ DRUCKER_PRAGER材料模型[13]。在眾多混凝土動態(tài)模型中,*MAT_JOHNSON_HOLMQUIST_CONCRETE模型可用于模擬大應變、高應變率和高壓下的混凝土[13],并被廣泛應用于混凝土動載問題。在本數(shù)值計算模型中采用該模型對混凝土管道進行模擬。在混凝土承插口處設置橡膠圈,橡膠圈的材料采用*MAT_MOONEY?RIVLIN_RUBBER模型。為模擬現(xiàn)場試驗中炮孔的炮泥堵塞,數(shù)值計算模型中炮泥采用*MAT_ PLASTIC_KINEMATIC材料模型[13]。

    現(xiàn)場試驗中爆破采用2#巖石乳化炸藥,單孔連續(xù)裝藥量為8 kg。采用LS?DYNA軟件自帶的高能炸藥材料*MAT_HIGH EXPLOSIVE_BURN來模擬炸藥模型[13],炸藥密度為1150 kg/m?,炸藥爆速4000 m/s,同時參考相關文獻[14],二號巖石乳化炸藥計算參數(shù)取值如表2所示。JWL方程可以對炸藥的爆轟壓力與相對體積和內能的關系進行描述[15],可通過狀態(tài)方程參數(shù)控制不同的炸藥材料,如下式所示

    式中 V為爆轟產(chǎn)物的相對比容;e為爆轟產(chǎn)物的比內能;A,B,R1,R2,ω為特征參數(shù)。

    炸藥采用ALE算法,ALE算法允許炸藥、巖石和空氣自由穿梭于彼此間,通過物質的流動可以分析爆炸發(fā)生的過程。其能夠很好地避免計算中網(wǎng)格變形過大導致的計算中斷問題。因此,ALE算法能夠很好地應用于爆炸分析計算中[16]。土壤、巖石、混凝土、炮泥、橡膠等結構采用拉格朗日算法,該算法多用于固體結構的應力應變分析,這種方法以物質坐標為基礎,其所描述的網(wǎng)格單元將以類似‘雕刻的方式劃分在用于分析的結構上,主要的優(yōu)點是能夠非常精確地描述結構邊界的運動。炸藥與土體相互作用利用流固耦合算法來實現(xiàn),通過添加流固耦合關鍵字進行處理。

    2.3 數(shù)值模擬結果驗證

    為驗證數(shù)值模擬計算結果是否合理,在管內徑為1 m的模型中選取和現(xiàn)場試驗振速監(jiān)測點相同位置的監(jiān)測點D1,D2,D3,D4,D5和D6,對所選取的監(jiān)測點X,Y和Z方向的振動速度進行監(jiān)測,同時和現(xiàn)場監(jiān)測數(shù)據(jù)進行對比,如表3所示;并將同一條件下現(xiàn)場實測和數(shù)值模擬監(jiān)測點D3的Z方向振速波形對比如圖4所示。

    由表3和圖4可知:數(shù)值計算結果和現(xiàn)場監(jiān)測數(shù)據(jù)誤差率的平均值為7.65%。數(shù)值計算結果和現(xiàn)場監(jiān)測的振動速度時程曲線變化基本吻合,可見數(shù)值模擬以及計算材料參數(shù)、算法的選取較為合理。并且由表3可知,三個方向中的Z方向振動速度最大,X方向次之,Y方向最小。在6個監(jiān)測點位中D3和D4的振動速度最大,即管道中部位置。并且振動速度向混凝土管道兩端逐漸降低。

    3 反映混凝土管道尺寸效應的爆破振動預測模型

    爆破地震波在巖土體內傳播過程中,地震波的衰減受到的影響因素較多,文獻[17?20]研究了爆破地震波在不同條件下的衰減規(guī)律。傳播路徑介質(巖土體的性質、巖體節(jié)理面等)、爆源的強度、爆源到監(jiān)測點的距離、地形地貌是影響爆破地震波傳播衰減的主要因素??紤]現(xiàn)場試驗的條件將混凝土管內徑影響下的爆破地震波傳播衰減規(guī)律所涉及的主要變量進行歸納,如表4所示。

    由量綱分析白金漢定理(π定理)(Langhaar, 1951),地表巖土體質點峰值振動速度v可表示為

    式(10)中:k3,k4,α2,β1分別為函數(shù)變換過程中給定的系數(shù);式(11)中:-α2ln r表示爆破振動速度隨距離r的變化;α2為衰減指數(shù),主要反映場地介質條件的影響;而k3+(α2lnQ)/3綜合反映了傳播路徑介質條件與炸藥量對巖土體質點振動的貢獻。

    由上可以得到

    式中 k為場地影響系數(shù);α2表示為爆破振動衰減系數(shù);β2為混凝土管內徑影響下的效應系數(shù)。

    4 承插式混凝土管道爆破振動響應特性分析

    在城市的混凝土管道設施建設中,市政中排水管道在主線中尺寸較大,而支線管道尺寸較小。為分析承插式混凝土管道不同的直徑對管道上方地表土以及管道爆破振動響應的影響,在只改變管內徑的情況下共進行5種不同工況的數(shù)值模擬(其中D為承插式混凝土管的內徑),如表5所示。

    4.1 不同內徑承插式混凝土管道應力分析

    圖5顯示了管內徑為1 m模型中混凝土管在不同時刻的應力分布。在選取的分析時間中,在3250 μs時,和爆源正對方向的管道承插口和臨近承插口的位置處最先出現(xiàn)較高應力;在7199 μs時,爆源正對方向管道承插口最大應力達到2.234 MPa,隨著時間的變化,混凝土管道應力進一步增大,并從管道中部位置沿軸線向兩端擴展;在8122 μs時,最大應力出現(xiàn)在爆源正對方向管道承插口附近管節(jié)的中部位置,最大應力達到5.324 MPa;在15200 μs時,應力區(qū)向管道環(huán)向進一步擴展,管道承插口處應力降低,管身位置出現(xiàn)高應力區(qū)。

    將其他內徑的混凝土管道在15200 μs時迎爆側應力分布圖與內徑為1 m的模型進行對比,如圖6所示。隨著管徑的增大,高應力區(qū)以及最大應力值減少,在管內徑為1.35 m時,最大應力值急劇減少為1.125 MPa。

    結合圖5和6可知,在爆破地震波作用下,承插式混凝土管道在振動響應過程中高應力區(qū)主要集中在迎爆側。在不同內徑承插式混凝土管道的迎爆側,如圖7所示選取6個監(jiān)測點位,讀取每個監(jiān)測點的峰值應力進行對比,如圖8所示。內徑為0.9 m的承插式混凝土管道管身出現(xiàn)的峰值應力最大為13.51 MPa,承插口出現(xiàn)的峰值應力最大為3.62 MPa,管身和承插口應力差值存在較大差值,隨著承插式混凝土管道管徑的增大,監(jiān)測點的峰值應力急劇降低,承插口與管身的峰值應力差值降低,當管徑達到1.35 m時,承插式混凝土管道管身出現(xiàn)的峰值應力最大為0.91 MPa,承插口出現(xiàn)的峰值應力最大為0.78 MPa,承插口與管身峰值應力基本趨于一致。在爆破地震波作用下,承插式混凝土管道應力最大值出現(xiàn)在監(jiān)測點1和點3位置,由此可確定監(jiān)測點1和點3所在截面為承插式混凝土管道在爆破地震波作用下的危險截面。

    4.2 不同內徑承插式混凝土管道振速分析

    為分析5種不同工況下承插式混凝土管道危險截面處的動力響應特征,將圖7中點1和點3所在的危險截面等分為12部分,用以觀察管道截面不同位置的質點峰值振動速度,如圖9所示,承插式混凝土管道的迎爆側即270°?360°位置處的質點峰值振速最大,隨著與爆源距離的增加,質點的峰值振動速度逐漸降低。隨著管道內徑的增加,危險截面上的質點峰值振動速度逐漸降低。管道內徑為0.9 m時,迎爆側質點峰值振動速度最大為22.40 cm/s;當管徑增加到1.35 m時,迎爆側質點峰值振動速度最大值降低至6.65 cm/s,減少了70.3%。由此可知,在研究爆破地震波作用下管道的動力響應問題時,應該考慮管道尺寸這個影響因素。

    4.3 管道上方地表土體在管道內徑影響下的爆破振動規(guī)律

    在不同管道內徑的數(shù)值計算模型中,如圖2所示,在承插式混凝土管道的危險截面處正上方地表土體選取監(jiān)測點位H1和H2,由于爆源和測點連線方向即Z方向的振動速度最大,故只對測點Z方向的振速進行了統(tǒng)計,不同管內徑模型的爆破振動數(shù)據(jù)如表6所示。由表6可知,隨著混凝土管半徑以及測點和爆源間距離的增大,模型中地表監(jiān)測點Z方向的峰值振動速度依次遞減。

    采用考慮混凝土管半徑影響的地表土體爆破振動速度規(guī)律的數(shù)學模型式(13)對表6中所列測試結果進行回歸分析,得到混凝土管內徑影響的地表土體爆破振動規(guī)律的預測模型,結果如下式所示

    式中 v為管道上方土體峰值振動速度,D為管道內徑,相關系數(shù)R2=0.950,可靠度較高,能比較好地反映v和D之間關系。由式(14)可知,當測點與爆源距離一定時,管道上方土體峰值振速振速隨著管內徑的增加而逐漸降低。

    5 基于管身破壞的管道振動速度安全判據(jù)

    5.1 管道和上方地表土體振速對比分析

    由于管道埋置于地下,現(xiàn)場工程條件下不便于進行開挖揭露。 參考朱斌軍等[21]將管道上方地表振動速度作為更加直觀的爆破振動觀測值,對管道與正上方地表X,Y,Z合方向振速進行線性統(tǒng)計。在每個模型中選取混凝土管正上方土體表面H1,H2處和對應管道截面上的Z方向的峰值振動速度進行線性擬合,結果如圖10所示。

    承插式混凝土管道和正上方地表土體Z方向峰值振動速度統(tǒng)計關系為

    式中 v1為承插式混凝土管道危險截面上最大Z方向的峰值振動速度(cm/s); v為管道正上方地表土體Z方向的峰值振動速度(cm/s)。由圖10可知,混凝土管道危險截面峰值振速大于正上方地表土體,這是由于混凝土管道與爆源間的距離較小,隨著與爆源之間距離的增加,爆破振動響應也隨之降低,管道和土體結構的不同也導致在爆破地震波作用下,兩者的振動響應有所區(qū)別。由式(15)可知,混凝土管道危險截面上最大的Z方向峰值振動速度與正上方地表土體呈線性關系,隨著地表土體峰值振速的增加而增加。結合式(14)可知,當測點與爆源距離一定時,管道內部Z方向的峰值振動速度隨著管內徑的增加而逐漸降低。這是由于隨著混凝土管道內徑的增加,管土間的接觸面積也隨之增加,周圍土體對混凝土管的約束作用也隨之增強,從而導致混凝土管道在爆破地震作用下的峰值振動速度降低。

    5.2 管道單元拉應力與峰值振速關系

    混凝土材料的抗壓強度較大,而抗拉強度較小,在爆破振動響應中,混凝土管道主要是產(chǎn)生受拉破壞[22]。在動態(tài)荷載作用下,不考慮應變速率時,混凝土結構的極限動抗拉強度取2.099 MPa[23],參考混凝土結構設計規(guī)范中給出的C35混凝土軸心抗拉強度標準值?tk=2.20 MPa,混凝土軸心抗拉強度設計值?t=1.57 MPa[24]。取a=?tk/?t=1.40為混凝土結構極限動抗拉強度的修正系數(shù),得到混凝土管道在正常使用條件下的極限動抗拉強度?=1.50 MPa。

    工程實際表明,混凝土管道爆破振動破壞主要為爆破地震波在管道介質中傳播引起的管道應力變化而產(chǎn)生的管道材料本身的失效破壞。根據(jù)波動學理論,波在無限介質中傳播時產(chǎn)生的應力可按照下式計算

    式中 σ為正應力(MPa);ρ為傳播介質密度(g/cm3);c為波在介質中的傳播速度,按照波的類型可以分為縱波和橫波(cm/s);v為波傳播過程中引起的介質質點的振動速度(cm/s)。本研究中,爆破地震波在管道中傳播引起管道單元應力的變化,管道傳播介質一定,因此介質密度ρ以及波的傳播速度c可以近似看作定值,所以ρc可以看作常數(shù)k,故上式可以變?yōu)閼εc速度相關的一次函數(shù)。此外,有學者也以該理論為基礎,通過應力與速度得到一次函數(shù)關系來分析管道的安全性[25]。

    在每種工況中,將承插式混凝土管道危險截面上的峰值拉應力和峰值振速進行線性擬合,統(tǒng)計關系如表7所示。根據(jù)表7中管道峰值拉應力σt和峰值振動速度v2的統(tǒng)計關系,由混凝土管道正常使用條件下的極限動抗拉強度?可到不同內徑混凝土管的振動速度安全判據(jù),如表8所示。

    由表8可知,在保證混凝土管道正常使用的條件下,混凝土管道內徑的不同,管道在爆破地震波作用下的控制振速也有所不同,并隨著管內徑的增加,承插式混凝土管道的振動速度最大安全允許值逐漸減小。

    6 結 論

    本文通過承插式混凝土管道在爆破振動作用下的動力響應現(xiàn)場試驗,并結合數(shù)值計算模型進行工況補充,得到了不同尺寸承插式混凝土管道在爆破振動作用下的動力響應結果如下:

    (1)在爆破地震波的作用下,承插式混凝土管道的管身和管道承插口出現(xiàn)不同的應力響應,峰值應力出現(xiàn)的時間存在差別,管身產(chǎn)生的峰值應力大于管道承插口,并隨著管道內徑的增加,管道產(chǎn)生的峰值應力降低,管身和承插口的峰值應力差值也隨之降低。

    (2)基于量綱分析,結合數(shù)值計算模型的數(shù)據(jù),得到在管道尺寸影響下的管道爆破振動速度的變化規(guī)律?;炷凉艿涝诒频卣鸩ǖ淖饔孟?,產(chǎn)生的峰值振動速度隨著管道內徑的增加而逐漸降低。

    (3)隨著管道內徑的增加,以承插式混凝土管道在爆破振動作用下受拉破壞為判別依據(jù),混凝土管道在爆破振動作用下安全使用的最大振動速度逐漸降低。內徑為0.9,1,1.1,1.2,1.35 m的承插式混凝土管道的爆破振動速度控制閾值分別為23.41,22.95,19.77 ,16.66 ,14.17 cm/s。

    參考文獻:

    [1] AMIR Sajjad Abedi, NADER Hataf, ARSALAN Ghahramani. Analytical solution of the dynamic response of buried pipelines under blast wave[J]. International Journal of Rock Mechanics and Mining Sciences, 2016,88:201-306.

    [2] Zhang L, Liang Z, Zhang J. Mechanical response of a buried pipeline to explosion loading[J]. Journal of Failure Analysis and Prevention, 2016,16(4): 576-582.

    [3] 鄭爽英,楊立中.隧道爆破地震下輸氣管道動力響應數(shù)值試驗[J].西南交通大學學報,2017,52(2):264-271.

    Zheng S Y, Yang L Z. Numerical experiments of dynamic response of buried gas pipeline under the action of seismic waves induced by tunnel blasting[J]. Journal of Southwest Jiaotong University, 2017, 52(2):264-271.

    [4] 劉建民,陳文濤.爆炸荷載下埋地管道動力響應分析研究[J].工程爆破, 2008,14(2):20-24.

    Liu J M, Chen W T. Dynamic response study of buried pipeline subjected to blast loads[J]. Engineering Blasting, 2008,14(2): 20-24.

    [5] 鐘冬望,黃 雄,司劍峰,等.爆破荷載作用下埋地鋼管的動態(tài)響應實驗研究[J].爆破, 2018,35(2):19-25.

    Zhong D W, Huang X, Si J F,et al. Experimental study on dynamic response of buried pipeline under blasting loads[J].Blasting, 2018,35(2):19-25.

    [6] 郝郁清,程 康,趙其達拉圖,等.爆破地震波作用下燃氣管道的安全距離研究[J].爆破, 2018,35(2):137-143.

    Hao Y Q, Cheng K, Zhao Q D L T, et al. Study on safety distance of gas pipeline under action of blasting seismic wave[J]. Blasting, 2018,35(2):137-143.

    [7] 張紫劍,趙昌龍,張黎明,等.埋地管道爆破振動安全允許判據(jù)試驗探究[J].爆破, 2016,33(2):12-16.

    Zhang Z J, Zhao C L, Zhang L M, et al. Experimental investigation of blasting vibration safety criterion on buried pipeline[J]. Blasting, 2016,33(2):12-16.

    [8] 張 震,周傳波,路世偉,等.爆破振動作用下鄰近埋地混凝土管道動力響應特性[J].哈爾濱工業(yè)大學學報, 2017,49(9):79-84.

    Zhang Z, Zhou C B, Lu S W , et al. Dynamic response characteristic of adjacent buried concrete pipeline subjected to blasting vibration[J].Journal of Harbin Institute of Technology, 2017,49(9):79-84.

    [9] 高 壇,周傳波,蔣 楠,等.基坑開挖爆破下鄰近管道振動速度安全閾值研究[J].安全與環(huán)境學報, 2017,17(6):2191-2195.

    Gao T, Zhou C B, Jiang N, et al. Study on the vibration velocity thresh-old of the adjacent pipeline under the blasting excavation of the foundation pit [J]. Journal of Safety and Environment, 2017,17(6):2191-2195.

    [10] 張黎明,趙昌龍,池恩安,等.地下管道爆破振動安全允許峰值振速試驗研究[J].爆破, 2019,36(4):146-149+154.

    Zhang L M, Zhao C L, Chi E A, et al. Experimental investigation on allowed blasting vibration velocity of underground pipeline [J]. Blasting, 2019,36(4):146-149+154.

    [11] 沈麗華,王貫明,王憬山,等. GB/T 11836-2009,混凝土和鋼筋混凝土排水管[S]. 北京: 中國標準出版社,2009.

    Shen Lihua, Wang Guanming, Wang Jingshan, et al. GB/T 11836-2009, Concrete and reinforced concrete sewer pipes[S]. Beijing: Standards Press of China, 2009.

    [12] 中華人民共和國住房和城鄉(xiāng)建設部. GB 50268-2008,給水排水管道工程施工及驗收規(guī)范[S] . 北京: 中國建筑工業(yè)出版社,2008.

    MOHURD. GB 50268-2008,Code for construction and acceptance of water and sewerage pipeline works[S]. Beijing: China Architecture & Building Press, 2008.

    [13] HALLQUIST J. LS-DYNA Keyword Users Manual R8.0[M]. California: Livermore Software Technology Corporation,2015.

    [14] Xie L X, Lu W B, Zhang Q B, et al. Damage evolution mechanisms of rock in deep tunnels induced by cut blasting[J]. Tunnelling and Underground Space Technology Incorporating Trenchless Technology Research, 2016, 58, 257?270.

    [15] 趙 錚,陶 鋼,杜長星.爆轟產(chǎn)物JWL狀態(tài)方程應用研究[J].高壓物理學報, 2009,23(4):277-282.

    Zhao Z, Tao G, Du C X. Application research on JWL equation of state of detonation products[J]. Chinese Journal of High Pressure Physics, 2009,23(4):277-282.

    [16] 王崢崢,張楊生.基于ALE算法的隧道開挖爆破振動特性數(shù)值分析[J].大連理工大學學報,2017,57(3):279-284.

    Wang Z Z, Zhang Y S. Numerical analysis for blasting vibration characteristics in tunnel excavation based on ALE algorithm[J]. Journal of Dalian University of Technology, 2017,57(3):279-284.

    [17] 李繼業(yè),蒲朝欽,何興貴.節(jié)理巖體邊坡爆破振動衰減特征分析[J].爆破, 2018,35(2):32-36.

    Li J Y, Pu C Q, He X G. Analysis of blasting vibration attenuation characteristics of jointed rock mass slope [J]. Blasting, 2018,35(2):32-36.

    [18] 唐 海,李海波,蔣鵬燦,等.地形地貌對爆破振動波傳播的影響實驗研究[J].巖石力學與工程學報, 2007,26(9):1817-1823.

    Tang H, Li H B, Jiang P C, et al. Experimental study on effect of topography on propagation of blasting waves[J]. Chinese Journal of Rock Mechanics and Engineering, 2007,26(9):1817-1823.

    [19] 王智德,江俐敏,祝文化,等.順層巖質邊坡爆破荷載作用下的振動傳播規(guī)律研究[J].爆破, 2019,36(1):55-62.

    Wang Z d, Jiang L M, Zhu W H, et al. Vibration Propagation Characteristics of Bedding Rock Slope under Blasting Loads[J]. Blasting, 2019,36(1):55-62.

    [20] 呂 濤,石永強,黃誠,等.非線性回歸法求解爆破振動速度衰減公式參數(shù)[J].巖土力學, 2007,28(9):1871-1878.

    Lu T, Shi Y Q, Huang C, et al. Study on attenuation parameters of blasting vibration by nonlinear regression analysis[J]. Rock and Soil Mechanics, 2007,28(9):1871-1878.

    [21] 朱 斌,蔣 楠,賈永勝,等.下穿穿燃氣管道爆破振動效應現(xiàn)場試驗研究[J].巖石力學與工程學報,2019,38(12):2582-2592.

    Zhu B, Jiang N, Jia Y S, et al. Field experiment on blasting vibration effect of underpass gas pipelines[J].Chinese Journal of Rock Mechanics and Engineering, 2019,38(12):2582-2592.

    [22] 夏宇磬,蔣 楠,姚穎康,等.粉質黏土層預埋承插式混凝土管道對爆破振動的動力響應[J].爆炸與沖擊, 2020,40(4):73-83.

    XIA Y Q, Jiang N, Yao Y k, et al. Experiment study on dynamic response of spigot concrete pipe buried in silty clay layer subjected to blasting seismic wave[J]. Explosion and Shock Waves, 2020,40(4):73-83.

    [23] Jin X, Hou C, Fan X, et al. Quasi-static and dynamic experimental studies on the tensile strength and failure pattern of concrete and mortar discs [J].Scientific Reports, 2017,7(1):15305.

    [24] 中華人民共和國住房和城鄉(xiāng)建設部. GB50010-2010,混凝土結構設計規(guī)范(2015版) [S]. 北京: 中國建筑工業(yè)出版社,2015.

    MOHURD. GB50010-2010, Code for design of concrete structures (2015) [S]. Beijing: China Architecture & Building Press, 2015.

    [25] 朱 斌,蔣 楠,周傳波,等.基坑開挖爆破作用鄰近壓力燃氣管道動力響應特性研究[J].振動與沖擊, 2020,39(11):201-208.

    Zhu B, Jiang N, Zhou C B, et al. Effect of excavation blast vibration on adjacent buried gas pipeline in a foundation pit[J]. Journal of Vibration and Shock, 2020,39(11):201-208.

    作者簡介: 黃一文(1995-),男,碩士研究生。電話:17635076978;E- mail:h2428948778@163.com

    通訊作者: 蔣 楠(1986-),男,博士,副教授。電話:18164055091;E-mail:happyjohn@foxmail.com

    国产黄片视频在线免费观看| 久久久久久久久久成人| 日韩一本色道免费dvd| 高清视频免费观看一区二区| 人人妻人人澡人人爽人人夜夜| 两个人的视频大全免费| 日本免费在线观看一区| 在线观看av片永久免费下载| 国产精品一及| 国产国拍精品亚洲av在线观看| av卡一久久| a级毛色黄片| 在线观看一区二区三区激情| 观看免费一级毛片| 青春草视频在线免费观看| 久久精品国产亚洲av涩爱| 看免费成人av毛片| 亚洲av欧美aⅴ国产| 精品人妻视频免费看| 一本—道久久a久久精品蜜桃钙片| 精品一区在线观看国产| 日本爱情动作片www.在线观看| 亚洲精品日本国产第一区| 亚洲av不卡在线观看| 又爽又黄a免费视频| 亚洲欧美日韩东京热| 三级国产精品欧美在线观看| 99久久精品热视频| 国产一区亚洲一区在线观看| 成人特级av手机在线观看| 日韩免费高清中文字幕av| 干丝袜人妻中文字幕| 日韩电影二区| 女人久久www免费人成看片| 天天躁夜夜躁狠狠久久av| 成人黄色视频免费在线看| 国产精品一及| 校园人妻丝袜中文字幕| 亚洲国产色片| 国产日韩欧美亚洲二区| 尾随美女入室| 一个人看的www免费观看视频| 亚洲一区二区三区欧美精品| 日韩制服骚丝袜av| 精品少妇久久久久久888优播| 91在线精品国自产拍蜜月| 夜夜骑夜夜射夜夜干| 国产69精品久久久久777片| 精品人妻一区二区三区麻豆| av在线老鸭窝| 少妇高潮的动态图| 黄色一级大片看看| 天堂中文最新版在线下载| 九草在线视频观看| 国产一区二区三区综合在线观看 | 黄色配什么色好看| 午夜视频国产福利| 日韩,欧美,国产一区二区三区| 亚洲欧美日韩东京热| 日韩欧美 国产精品| 免费播放大片免费观看视频在线观看| 噜噜噜噜噜久久久久久91| 亚洲精品国产av成人精品| 香蕉精品网在线| 午夜免费男女啪啪视频观看| 欧美日韩综合久久久久久| 国产精品一区二区性色av| 夫妻性生交免费视频一级片| 国产黄色免费在线视频| 99久久中文字幕三级久久日本| 在线观看国产h片| 日韩免费高清中文字幕av| 亚洲精品久久午夜乱码| 久久ye,这里只有精品| 好男人视频免费观看在线| 在线免费观看不下载黄p国产| 伊人久久国产一区二区| 成人高潮视频无遮挡免费网站| 国产永久视频网站| 久久久久久久久大av| 深爱激情五月婷婷| 国产精品久久久久成人av| 极品教师在线视频| 国产精品久久久久久久久免| 女性被躁到高潮视频| 久久韩国三级中文字幕| 久久影院123| 一级二级三级毛片免费看| 狂野欧美白嫩少妇大欣赏| 99热这里只有是精品50| 国产精品不卡视频一区二区| 狠狠精品人妻久久久久久综合| 九九爱精品视频在线观看| 亚洲,欧美,日韩| 你懂的网址亚洲精品在线观看| 亚洲三级黄色毛片| 国产免费一区二区三区四区乱码| 免费高清在线观看视频在线观看| 精品亚洲成国产av| 久久97久久精品| 亚洲精品日韩在线中文字幕| 精品人妻一区二区三区麻豆| 黄色怎么调成土黄色| 国产高潮美女av| 尤物成人国产欧美一区二区三区| 日韩强制内射视频| 日本-黄色视频高清免费观看| 中文字幕制服av| 成人18禁高潮啪啪吃奶动态图 | 亚洲成人av在线免费| 亚洲精品456在线播放app| 男女边吃奶边做爰视频| 夫妻性生交免费视频一级片| 免费看日本二区| 亚洲精品国产成人久久av| 欧美日韩一区二区视频在线观看视频在线| 国产亚洲一区二区精品| 大片电影免费在线观看免费| 久久久久久久久久人人人人人人| 国产v大片淫在线免费观看| 免费av中文字幕在线| 欧美日韩一区二区视频在线观看视频在线| 欧美日韩亚洲高清精品| 丰满迷人的少妇在线观看| 男的添女的下面高潮视频| 一个人免费看片子| av卡一久久| 黄片wwwwww| 国产伦在线观看视频一区| 欧美xxxx性猛交bbbb| 简卡轻食公司| 国产精品人妻久久久影院| 一个人看的www免费观看视频| 嘟嘟电影网在线观看| 国产精品一区二区在线不卡| 交换朋友夫妻互换小说| 久久久a久久爽久久v久久| 欧美成人a在线观看| 在线观看一区二区三区| 亚洲在久久综合| 国产69精品久久久久777片| 国产精品爽爽va在线观看网站| 91久久精品电影网| 亚洲人成网站高清观看| 能在线免费看毛片的网站| 又大又黄又爽视频免费| 免费看光身美女| 久热久热在线精品观看| 国产精品女同一区二区软件| 国产欧美日韩一区二区三区在线 | 色视频www国产| 噜噜噜噜噜久久久久久91| 午夜日本视频在线| 成人午夜精彩视频在线观看| 亚洲欧美中文字幕日韩二区| 自拍欧美九色日韩亚洲蝌蚪91 | 另类亚洲欧美激情| 亚洲欧美一区二区三区国产| 在线观看一区二区三区| 日韩在线高清观看一区二区三区| 亚洲精品第二区| 国产高清国产精品国产三级 | 91精品伊人久久大香线蕉| 最新中文字幕久久久久| 国产欧美亚洲国产| 中文字幕制服av| 久久国产亚洲av麻豆专区| 日韩不卡一区二区三区视频在线| 亚洲欧美日韩东京热| 久久久久久久国产电影| 91精品伊人久久大香线蕉| 99久久精品热视频| 美女脱内裤让男人舔精品视频| av线在线观看网站| 日日啪夜夜爽| 国产成人精品福利久久| 22中文网久久字幕| 男女边摸边吃奶| av在线观看视频网站免费| 国产精品国产三级专区第一集| 免费高清在线观看视频在线观看| 亚洲欧洲日产国产| 99热网站在线观看| 天天躁夜夜躁狠狠久久av| 国产成人一区二区在线| 久久久久久久国产电影| 中文字幕精品免费在线观看视频 | 在线观看av片永久免费下载| 91狼人影院| 国产91av在线免费观看| 婷婷色综合www| 多毛熟女@视频| 亚洲国产色片| 久久婷婷青草| 色5月婷婷丁香| 精品酒店卫生间| 你懂的网址亚洲精品在线观看| 在线看a的网站| 国产乱人视频| 日本av免费视频播放| 纵有疾风起免费观看全集完整版| 久久av网站| 新久久久久国产一级毛片| 狂野欧美白嫩少妇大欣赏| 久久97久久精品| 在线观看三级黄色| 一区二区三区精品91| 欧美最新免费一区二区三区| 青春草视频在线免费观看| a级毛色黄片| 女人十人毛片免费观看3o分钟| 久久久久久久久大av| 久久精品人妻少妇| 日韩,欧美,国产一区二区三区| videos熟女内射| 51国产日韩欧美| 一区二区三区四区激情视频| 亚洲av.av天堂| 日韩av在线免费看完整版不卡| 免费看光身美女| 中文字幕制服av| 我要看日韩黄色一级片| 国产一区二区在线观看日韩| 精品久久久噜噜| 久久久色成人| 久久久成人免费电影| 午夜福利高清视频| 成人高潮视频无遮挡免费网站| 亚洲国产毛片av蜜桃av| 三级国产精品片| 国模一区二区三区四区视频| 97在线视频观看| 亚洲国产精品成人久久小说| 视频中文字幕在线观看| 高清日韩中文字幕在线| 日日啪夜夜撸| 精品国产一区二区三区久久久樱花 | 视频区图区小说| 永久免费av网站大全| 人妻制服诱惑在线中文字幕| 一级爰片在线观看| 热re99久久精品国产66热6| 婷婷色综合www| 人人妻人人添人人爽欧美一区卜 | 激情五月婷婷亚洲| 久久久久精品久久久久真实原创| 亚洲欧美精品自产自拍| 亚洲,欧美,日韩| 久久国产精品男人的天堂亚洲 | 丰满迷人的少妇在线观看| 丰满人妻一区二区三区视频av| 蜜臀久久99精品久久宅男| av天堂中文字幕网| 尾随美女入室| 国产一区有黄有色的免费视频| 国产精品国产三级国产专区5o| 97热精品久久久久久| 成人漫画全彩无遮挡| 亚洲自偷自拍三级| 精品久久久久久久久亚洲| 人妻夜夜爽99麻豆av| 性高湖久久久久久久久免费观看| 人妻 亚洲 视频| 在线观看人妻少妇| 欧美 日韩 精品 国产| 亚洲国产精品成人久久小说| 欧美成人a在线观看| 超碰av人人做人人爽久久| 又粗又硬又长又爽又黄的视频| 亚洲欧美成人综合另类久久久| 久久久午夜欧美精品| 日本免费在线观看一区| 日韩,欧美,国产一区二区三区| 亚洲av中文字字幕乱码综合| 高清在线视频一区二区三区| 婷婷色av中文字幕| 国产黄色视频一区二区在线观看| 国模一区二区三区四区视频| 色视频在线一区二区三区| 国国产精品蜜臀av免费| 亚洲av男天堂| av免费观看日本| 亚洲精品乱码久久久久久按摩| 插逼视频在线观看| h日本视频在线播放| 人人妻人人爽人人添夜夜欢视频 | 国产亚洲最大av| xxx大片免费视频| 久久久久久伊人网av| 欧美bdsm另类| 亚洲精品,欧美精品| 国内少妇人妻偷人精品xxx网站| 男人舔奶头视频| 18禁在线播放成人免费| 国产成人aa在线观看| 涩涩av久久男人的天堂| 亚洲在久久综合| 久久热精品热| 七月丁香在线播放| 欧美日韩视频高清一区二区三区二| 伦理电影免费视频| a级毛片免费高清观看在线播放| 一级a做视频免费观看| 大码成人一级视频| 七月丁香在线播放| 日日撸夜夜添| 最近最新中文字幕免费大全7| 99久久精品热视频| 草草在线视频免费看| 99热国产这里只有精品6| av国产免费在线观看| 日本午夜av视频| 久久精品国产自在天天线| 亚洲一级一片aⅴ在线观看| 99久国产av精品国产电影| 内地一区二区视频在线| 国产精品99久久久久久久久| 交换朋友夫妻互换小说| 日本wwww免费看| 少妇裸体淫交视频免费看高清| 精品一区二区三区视频在线| 一级爰片在线观看| 少妇的逼水好多| 超碰97精品在线观看| 男人爽女人下面视频在线观看| 精品一区二区三区视频在线| 在线免费十八禁| 久久精品夜色国产| 国产精品国产三级国产专区5o| 日本黄色片子视频| 免费看av在线观看网站| 嘟嘟电影网在线观看| 久久久久久久大尺度免费视频| 久久久久久久国产电影| 久久精品国产a三级三级三级| 国产乱人偷精品视频| 中文字幕人妻熟人妻熟丝袜美| 又爽又黄a免费视频| 少妇精品久久久久久久| 国产免费又黄又爽又色| 久热这里只有精品99| 熟女av电影| www.av在线官网国产| 性高湖久久久久久久久免费观看| 免费av不卡在线播放| 夜夜看夜夜爽夜夜摸| 永久免费av网站大全| 国产极品天堂在线| 欧美精品一区二区免费开放| 日韩伦理黄色片| 亚洲国产精品专区欧美| 亚洲精品色激情综合| 免费av中文字幕在线| 超碰av人人做人人爽久久| 一二三四中文在线观看免费高清| 97热精品久久久久久| 男人爽女人下面视频在线观看| 少妇 在线观看| 亚洲精品456在线播放app| 亚洲美女搞黄在线观看| 亚洲av成人精品一区久久| 亚洲丝袜综合中文字幕| 日韩一区二区视频免费看| 国产成人一区二区在线| 女性生殖器流出的白浆| 午夜福利在线观看免费完整高清在| 在线观看av片永久免费下载| 深夜a级毛片| 国产精品一及| av一本久久久久| 啦啦啦视频在线资源免费观看| 亚洲精品国产av成人精品| 国产黄色视频一区二区在线观看| 一级av片app| 国产日韩欧美亚洲二区| a级一级毛片免费在线观看| 国产在线视频一区二区| 国产精品久久久久久久电影| 中文字幕av成人在线电影| 免费观看a级毛片全部| 国产精品国产三级国产av玫瑰| 亚洲国产欧美人成| 久久久久久久精品精品| 夜夜骑夜夜射夜夜干| 在线免费观看不下载黄p国产| 在线亚洲精品国产二区图片欧美 | 欧美最新免费一区二区三区| 日韩欧美 国产精品| 免费观看a级毛片全部| 久久久久视频综合| 舔av片在线| 丰满少妇做爰视频| 一级黄片播放器| 狠狠精品人妻久久久久久综合| 精品一区二区三卡| 国产一区二区三区综合在线观看 | av在线app专区| 亚洲第一av免费看| 18+在线观看网站| 日本黄色日本黄色录像| 三级国产精品欧美在线观看| 亚洲欧美日韩无卡精品| 国产精品国产三级国产av玫瑰| 青春草亚洲视频在线观看| 好男人视频免费观看在线| 欧美精品亚洲一区二区| 亚洲不卡免费看| 国产精品人妻久久久久久| 国产av精品麻豆| 91久久精品电影网| 菩萨蛮人人尽说江南好唐韦庄| 噜噜噜噜噜久久久久久91| 欧美xxxx黑人xx丫x性爽| 建设人人有责人人尽责人人享有的 | 欧美成人a在线观看| 久久精品国产亚洲网站| 乱码一卡2卡4卡精品| 99久久精品一区二区三区| 久久国产精品大桥未久av | 国产白丝娇喘喷水9色精品| 亚洲精品视频女| 日韩大片免费观看网站| 久久精品久久精品一区二区三区| 精品人妻熟女av久视频| 精品人妻一区二区三区麻豆| 熟女av电影| 亚洲第一av免费看| 国产黄片美女视频| 日本与韩国留学比较| 国产男女内射视频| 97热精品久久久久久| 国产精品伦人一区二区| 小蜜桃在线观看免费完整版高清| 99精国产麻豆久久婷婷| 亚洲成人中文字幕在线播放| 人妻少妇偷人精品九色| 成人亚洲精品一区在线观看 | 高清毛片免费看| 精品人妻视频免费看| 性色avwww在线观看| 男女无遮挡免费网站观看| 只有这里有精品99| av卡一久久| 成人一区二区视频在线观看| 国产午夜精品久久久久久一区二区三区| 97热精品久久久久久| 老熟女久久久| 成人亚洲精品一区在线观看 | 精品午夜福利在线看| 美女内射精品一级片tv| 夫妻午夜视频| 久久久久久伊人网av| 一级黄片播放器| 免费看不卡的av| 午夜福利在线在线| 日韩av免费高清视频| 春色校园在线视频观看| 狠狠精品人妻久久久久久综合| 精品少妇黑人巨大在线播放| 汤姆久久久久久久影院中文字幕| 一级毛片 在线播放| 欧美成人一区二区免费高清观看| 国产欧美日韩一区二区三区在线 | 久久国产乱子免费精品| 中国国产av一级| 久久精品国产亚洲网站| 久久毛片免费看一区二区三区| 欧美精品人与动牲交sv欧美| 又粗又硬又长又爽又黄的视频| 毛片女人毛片| 国产成人a区在线观看| 国产有黄有色有爽视频| 日日撸夜夜添| 日韩欧美精品免费久久| 尤物成人国产欧美一区二区三区| 国产精品一区二区在线观看99| 久久精品国产自在天天线| 一区二区三区四区激情视频| 香蕉精品网在线| 亚洲最大成人中文| 五月伊人婷婷丁香| 国产成人免费无遮挡视频| 亚洲人成网站在线播| 国产精品一区二区三区四区免费观看| 观看av在线不卡| 在线观看免费高清a一片| 高清毛片免费看| 色网站视频免费| 一本—道久久a久久精品蜜桃钙片| 午夜福利影视在线免费观看| av在线老鸭窝| 人妻夜夜爽99麻豆av| 国产亚洲午夜精品一区二区久久| 国产精品99久久99久久久不卡 | 亚洲综合精品二区| 亚洲第一av免费看| 午夜福利高清视频| 噜噜噜噜噜久久久久久91| 热re99久久精品国产66热6| 国产一区二区三区综合在线观看 | 大香蕉久久网| 亚洲国产精品专区欧美| 亚洲精品国产av成人精品| 国产日韩欧美亚洲二区| 在线看a的网站| 亚洲精品国产色婷婷电影| 国产69精品久久久久777片| 久久鲁丝午夜福利片| 精品久久久久久电影网| av不卡在线播放| 成年女人在线观看亚洲视频| 亚洲av国产av综合av卡| 精品少妇久久久久久888优播| 欧美高清性xxxxhd video| 最近中文字幕2019免费版| 一级片'在线观看视频| 一级黄片播放器| 内射极品少妇av片p| 黄片无遮挡物在线观看| 高清不卡的av网站| 久久久久久久大尺度免费视频| 99精国产麻豆久久婷婷| 国产男女内射视频| 成人二区视频| 精品国产一区二区三区久久久樱花 | 成人毛片a级毛片在线播放| 国产黄色视频一区二区在线观看| 蜜桃在线观看..| 日韩不卡一区二区三区视频在线| 久久国内精品自在自线图片| 男女无遮挡免费网站观看| 91aial.com中文字幕在线观看| 亚洲成人手机| 久久婷婷青草| 亚洲精品日韩av片在线观看| 最后的刺客免费高清国语| 亚洲av中文字字幕乱码综合| 一区在线观看完整版| 最后的刺客免费高清国语| 欧美变态另类bdsm刘玥| 少妇被粗大猛烈的视频| 久久久色成人| 国产在线视频一区二区| 国产欧美另类精品又又久久亚洲欧美| 熟女人妻精品中文字幕| 人妻制服诱惑在线中文字幕| 亚洲自偷自拍三级| 欧美3d第一页| 欧美三级亚洲精品| 亚洲国产精品一区三区| 丝瓜视频免费看黄片| 国产成人精品久久久久久| 国产精品精品国产色婷婷| 久久精品国产鲁丝片午夜精品| 国产有黄有色有爽视频| 中国国产av一级| videos熟女内射| 国产永久视频网站| 久久青草综合色| 精品一区二区三卡| 精品久久久久久久久亚洲| 国产淫片久久久久久久久| 九草在线视频观看| 男女边吃奶边做爰视频| 欧美高清性xxxxhd video| 国产乱人视频| 97热精品久久久久久| 日日摸夜夜添夜夜添av毛片| 97超视频在线观看视频| 一级av片app| 亚洲av欧美aⅴ国产| 六月丁香七月| 97超碰精品成人国产| 爱豆传媒免费全集在线观看| 成人二区视频| 久久久久久九九精品二区国产| 毛片一级片免费看久久久久| 在线播放无遮挡| 免费不卡的大黄色大毛片视频在线观看| 国产精品久久久久久精品古装| 两个人的视频大全免费| 大片免费播放器 马上看| 91狼人影院| 久久久久视频综合| 亚洲人成网站在线播| 久久99精品国语久久久| 久久国内精品自在自线图片| www.av在线官网国产| 在线观看国产h片| 成人一区二区视频在线观看| 亚洲综合色惰| 晚上一个人看的免费电影| 亚洲精品自拍成人| 亚洲四区av| 国产av精品麻豆| 欧美性感艳星| 国产欧美日韩精品一区二区| 少妇人妻久久综合中文| 欧美 日韩 精品 国产| 亚洲av欧美aⅴ国产| 免费观看的影片在线观看| 精品少妇黑人巨大在线播放| 91精品伊人久久大香线蕉| 欧美区成人在线视频| 下体分泌物呈黄色| 91午夜精品亚洲一区二区三区| 亚洲精品第二区| 国产精品成人在线| 欧美 日韩 精品 国产| 少妇丰满av| 尤物成人国产欧美一区二区三区| av免费在线看不卡| 日本欧美视频一区| 欧美国产精品一级二级三级 | 特大巨黑吊av在线直播|