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

    基于FLAC3D的巖質(zhì)邊坡爆破動(dòng)力響應(yīng)規(guī)律研究

    2019-10-20 04:50唐旭方正峰鄒飛
    人民長(zhǎng)江 2019年3期
    關(guān)鍵詞:爆破

    唐旭 方正峰 鄒飛

    摘要:邊坡的動(dòng)力穩(wěn)定一直以來(lái)都是巖土工程及地震工程重點(diǎn)研究的問(wèn)題。利用FLAC3D數(shù)值計(jì)算方法,對(duì)爆破作用下巖質(zhì)邊坡振動(dòng)速度、位移、應(yīng)力等進(jìn)行模擬,并在與實(shí)測(cè)爆破振動(dòng)峰值進(jìn)行對(duì)比驗(yàn)證的基礎(chǔ)上,分析邊坡的動(dòng)力響應(yīng)規(guī)律。結(jié)果表明:地震波具有極強(qiáng)的方向性,主要朝施加荷載的正向傳播;在邊坡坡面和坡腳處因應(yīng)力波在邊坡自由面反射容易產(chǎn)生一定區(qū)域的振動(dòng)加強(qiáng)場(chǎng);應(yīng)力的響應(yīng)特征上,無(wú)論是拉應(yīng)力還是壓應(yīng)力其最大值均出現(xiàn)在坡腳附近,在實(shí)際的邊坡爆破開(kāi)挖過(guò)程中,對(duì)邊坡的坡面尤其是坡腳處應(yīng)該加強(qiáng)保護(hù)措施。研究成果對(duì)巖質(zhì)邊坡開(kāi)挖設(shè)計(jì)及施工具有一定的參考價(jià)值。

    關(guān)鍵詞:動(dòng)力響應(yīng);巖質(zhì)邊坡;FLAC3D;爆破

    中圖法分類號(hào):P642

    文獻(xiàn)標(biāo)志碼:A

    DOI:l0.16232/j.cnki.1001-4179.2019.03.035

    1研究背景

    爆破作為一種最有效、最經(jīng)濟(jì)和最便捷的施工方法已經(jīng)普遍應(yīng)用到山地公路路塹巖質(zhì)高邊坡的形成中。巨大的工程規(guī)模伴隨而來(lái)的是大量的爆破施工作業(yè),而邊坡的動(dòng)力穩(wěn)定一直以來(lái)都是巖土工程及地震工程重點(diǎn)研究的問(wèn)題之一。

    隨著科技的發(fā)展,越來(lái)越多的研究人員借助計(jì)算機(jī)進(jìn)行仿真模擬,提出了各類數(shù)值分析方法,并開(kāi)發(fā)了相應(yīng)的數(shù)值分析軟件,如基于動(dòng)力有限元法的Ansys/Ls-dyna模塊[1],離散元法的UDEC[2],快速拉格朗日法的FLAC3D等[3-4]。

    蔣楠用Ansys/Ls-dyna建立了不同坡度的邊坡爆破振動(dòng)數(shù)值計(jì)算模型,獲得了不同坡度邊坡在同一水平位置的各坡面監(jiān)測(cè)點(diǎn)的質(zhì)點(diǎn)爆破振動(dòng)速度變化趨勢(shì),且同一爆源距的質(zhì)點(diǎn)峰值振速隨高程增大表現(xiàn)出高程放大效應(yīng);在邊坡坡度相近時(shí),質(zhì)點(diǎn)振速隨爆源距及高程的同時(shí)增加,高程效應(yīng)不明顯;不同坡度邊坡,在同一水平質(zhì)點(diǎn)振速隨著邊坡坡度增加以衰減為主,即質(zhì)點(diǎn)振速在爆源距增加及高程增加的同時(shí)作用下,而不再表現(xiàn)為常規(guī)的振動(dòng)速度隨爆源距增加衰減的現(xiàn)象,而是以質(zhì)點(diǎn)振速隨高程增加而呈現(xiàn)增長(zhǎng)的現(xiàn)象[5]。唐海,李海波等應(yīng)用UDEC程序模擬了爆破振動(dòng)波在不同地貌中的傳播特征,得出凸形地貌下放大系數(shù)并不隨臺(tái)階型凸形地貌的增高而單調(diào)增加,在孤立凸形地貌下,放大系數(shù)還與凸形地貌的高寬比有關(guān)[6-7]。言志信、祁生文等利用FLAC3D軟件對(duì)不同坡高、坡角的邊坡動(dòng)力響應(yīng)特征做了數(shù)值計(jì)算分析研究,獲知當(dāng)邊坡坡高較低時(shí),在一定的坡高范圍內(nèi)邊坡內(nèi)的質(zhì)點(diǎn)振動(dòng)速度隨著高程的增大表現(xiàn)出高程放大效應(yīng),且速度放大系數(shù)等值線圖呈近等間距分布;當(dāng)邊坡的坡高高出一定的范圍時(shí),邊坡質(zhì)點(diǎn)振動(dòng)速度不再隨高程的增大表現(xiàn)出線性的反放大,而是表現(xiàn)出節(jié)律性的變化,在坡肩附近時(shí),才再一次表現(xiàn)出放大的跡象[8-9]何錚等結(jié)合西南某水電戰(zhàn)順層巖質(zhì)高邊坡工程實(shí)例,研究了該邊坡地震后工程地質(zhì)特性,恢復(fù)邊坡震前的地形地貌,建立地震前未滑塌邊坡的三維模型,而后利用FLAC3對(duì)模型進(jìn)行動(dòng)力-時(shí)程響應(yīng)分析,繪制出位移、速度、加速度三量放大系數(shù)等值線圖[10]。

    巖質(zhì)高邊坡爆破開(kāi)挖過(guò)程中引起的地震效應(yīng)已成為影響邊坡穩(wěn)定性的重要誘發(fā)因素之一,越來(lái)越多的學(xué)者也開(kāi)始關(guān)注巖質(zhì)高邊坡在爆破荷載作用下的動(dòng)力響應(yīng)特征問(wèn)題。爆破荷載對(duì)邊坡的穩(wěn)定性影響主要是動(dòng)荷載作用下巖體的響應(yīng)特征問(wèn)題[11-12],因此掌握并深入分析巖質(zhì)高邊坡對(duì)爆破振動(dòng)效應(yīng)的響應(yīng)特征規(guī)律,探索爆破開(kāi)挖施工對(duì)邊坡穩(wěn)定性的影響,對(duì)于進(jìn)一步控制和預(yù)測(cè)爆破振動(dòng)效應(yīng)對(duì)巖質(zhì)邊坡的危害,同時(shí)減少爆破作業(yè)對(duì)周邊環(huán)境的影響具有重要的指導(dǎo)和實(shí)踐意義[13-15]。

    本文利用FLAC3D建立公路巖質(zhì)邊坡的數(shù)值分析模型,分析巖質(zhì)邊坡爆破振動(dòng)下振速、位移、應(yīng)力相應(yīng)的響應(yīng)特征,得到并分析巖質(zhì)邊坡爆破振動(dòng)下振速、位移、應(yīng)力相應(yīng)的響應(yīng)特征,細(xì)化研究傳播路徑中地質(zhì)條件、地貌特征的影響。

    2FLAC3D動(dòng)力計(jì)算原理

    FLAC3D的基本原理是采用線性顯式拉格朗日差分算法和混合離散技術(shù)求解運(yùn)動(dòng)方程,將指定計(jì)算區(qū)域劃分成若千個(gè)單元,單元之間由節(jié)點(diǎn)連接,在對(duì)某一節(jié)點(diǎn)施加荷載后,此節(jié)點(diǎn)的運(yùn)動(dòng)方程可以表示為時(shí)間不長(zhǎng)的有限差分形式,在一微小時(shí)間步長(zhǎng)內(nèi)作用于節(jié)點(diǎn)的荷載只對(duì)相鄰的一些節(jié)點(diǎn)產(chǎn)生影響,即根據(jù)某一t時(shí)刻的應(yīng)力狀態(tài)和△t時(shí)間步長(zhǎng)應(yīng)變?cè)隽?,進(jìn)而確定t+△t時(shí)刻的應(yīng)力狀態(tài)。根據(jù)高斯(Gauss)定律,由節(jié)點(diǎn)的速率求解單元的應(yīng)變?cè)隽?再利用應(yīng)力-應(yīng)變關(guān)系,即本構(gòu)方程計(jì)算單元應(yīng)力,并對(duì)其進(jìn)行積分,得到作用在節(jié)點(diǎn)上的應(yīng)力矢量;最后結(jié)合平衡方程,由節(jié)點(diǎn)力進(jìn)一步求解出節(jié)點(diǎn)速率及位移。如此反復(fù)循環(huán)迭代,這一過(guò)程將隨時(shí)間步長(zhǎng)正常擴(kuò)展到整個(gè)計(jì)算范圍直至計(jì)算收斂,從而準(zhǔn)確地模擬巖體及其他材料的塑性破壞和流動(dòng)。其計(jì)算流程如圖1所示。

    在速度和位移的動(dòng)力計(jì)算中,F(xiàn)LAC3D以節(jié)點(diǎn)為計(jì)算對(duì)象,分別將質(zhì)量和力都集中在節(jié)點(diǎn)上,再通過(guò)運(yùn)動(dòng)方程在時(shí)域內(nèi)進(jìn)行求解,節(jié)點(diǎn)運(yùn)動(dòng)方程可以表示為

    公式

    式中,F(xiàn)li(t)為t時(shí)刻l節(jié)點(diǎn)在i方向的平衡力分量,可由虛功原理導(dǎo)出;ml為l節(jié)點(diǎn)的集中質(zhì)量。

    將式(1)左邊用中心差分來(lái)近似,則可以得到節(jié)點(diǎn)的速度:

    公式

    此時(shí),可以進(jìn)一步得到節(jié)點(diǎn)位移:

    公式

    3數(shù)值計(jì)算模型的建立

    3.1計(jì)算模型

    圖2為某高速公路K86+152~K86+600段地形地貌示意圖,本文依托該工程對(duì)巖質(zhì)邊坡爆破振動(dòng)響應(yīng)特征進(jìn)行數(shù)值研究。

    本文采用三維計(jì)算模型,X正向?yàn)檫吰掠覀?cè),Y正向?yàn)榫€路大里程方向,Z軸鉛垂向上,模型計(jì)算區(qū)域?yàn)?52mx230mx137m。。輸入波動(dòng)的最短波長(zhǎng)控制著網(wǎng)格劃分的尺寸,研究表明,當(dāng)網(wǎng)格的最大尺寸為△l時(shí),輸入波動(dòng)的最短波長(zhǎng)為λ時(shí),必須滿足△lt<(1/10~1/8)λ。所建立模型及其網(wǎng)格劃分如圖3所示。

    3.2本構(gòu)模型及材料參數(shù)

    本文模型選取莫爾-庫(kù)倫(Mohr-Coulomb)彈塑性模型,該模型用于描述巖土體的力學(xué)行為。

    該里程段巖性較為均一,僅為灰?guī)r。模型中所需的巖石物理力學(xué)參數(shù)根據(jù)室內(nèi)實(shí)驗(yàn)確定,詳見(jiàn)表1。

    3.3吸收邊界

    在FLA3D中可以設(shè)置吸收邊界(也稱靜止邊界、黏性邊界)和自由場(chǎng)邊界來(lái)減少邊界中波的反射。

    本文將X正、負(fù)向,Y正、負(fù)向4個(gè)側(cè)面及Z負(fù)向地面均設(shè)置為黏性吸收邊界,地表均為自由邊界,詳見(jiàn)圖4。

    3.4等效荷載輸入

    炸藥在爆炸瞬間產(chǎn)生高溫高壓的爆轟氣體,作用于炮孔壁,使其產(chǎn)生破碎區(qū)和塑性變形區(qū),并且瞬間衰減為應(yīng)力波沿巖體繼續(xù)傳播,整個(gè)過(guò)程僅僅持續(xù)數(shù)毫秒。其間影響因素甚多且極其復(fù)雜,目前量測(cè)技術(shù)有限,無(wú)法量化并準(zhǔn)確地獲知爆炸過(guò)程中的每一個(gè)具體細(xì)節(jié),只能通過(guò)相關(guān)經(jīng)驗(yàn)將爆破過(guò)程中荷載作用及變化情況等效處理,通常是將數(shù)值模擬中假定荷載以應(yīng)力時(shí)程的形式輸入。本文擬采用三角脈沖荷載,這種荷載的建立需要確定兩個(gè)基本因素[16-19]。

    3.4.1爆破振動(dòng)峰值荷載

    峰值荷載的確定即確定炸藥爆炸過(guò)程中產(chǎn)生的爆轟氣體作用在炮孔壁上的峰值壓應(yīng)力。根據(jù)爆轟波的C-J理論,炸藥作用在炮孔壁上的平均炮轟壓力為

    公式

    式中,Pj為瞬時(shí)炸藥爆轟平均初始?jí)毫?ρe為炸藥密度;D;為炸藥爆轟速度;γ為炸藥的等熵系數(shù),y的取值與裝藥密度有關(guān)。研究表明[20],當(dāng)ρe<1.2g/cm3時(shí),γ取2.1,當(dāng)ρe≥1.2g/cm3時(shí),y取3。

    在耦合裝藥條件下,炸藥作用在炮孔壁上的初始峰值壓力P0為

    P0=Pj

    在不耦合裝條件下,炸藥作用在炮孔壁上的初始峰值壓力為

    公式

    式中,k為裝藥徑向不耦合系數(shù),

    公式

    分別為炮孔直徑和藥卷直徑。

    研究表明,爆轟氣體作用于炮孔壁時(shí),其壓力將明顯增大,故得到不耦合系數(shù)較小的柱狀藥包爆炸對(duì)炮孔壁的沖擊初始峰值壓力為

    公式

    式中,n為爆炸產(chǎn)物膨脹作用炮孔壁時(shí)的壓力增大系數(shù),n=8~10,本文取10;l為軸向不耦合系數(shù),

    公式

    ;le、lb分別為裝藥長(zhǎng)度和炮孔長(zhǎng)度。

    表2為爆破荷載峰值壓力計(jì)算相關(guān)計(jì)算參數(shù)。

    3.4.2爆破荷載作用時(shí)間

    炸藥爆炸作用時(shí)間非常短,通常認(rèn)為炸藥爆破時(shí)沖擊波作用的持續(xù)時(shí)間約1.0x10-6~0.1s,爆轟氣體壓力作用時(shí)間約為1.0x10-3~0.1s。本次數(shù)值計(jì)算中爆破荷載作用時(shí)間假設(shè)為7ms,其中爆破荷載壓力上升時(shí)間為1ms,'下降時(shí)間為6ms,且荷載按現(xiàn)場(chǎng)實(shí)際爆破設(shè)計(jì)情況單排爆破作用于K86+320~K86+350段第一級(jí)邊坡坡面上,圖5為爆破荷載施加示意圖。

    3.5力學(xué)阻尼的選取

    巖體中的內(nèi)摩擦以及可能存在的接觸表面滑動(dòng)是產(chǎn)生阻尼的主要原因。在動(dòng)力計(jì)算中,動(dòng)力反應(yīng)的阻尼需要在數(shù)值模型計(jì)算分析時(shí)模擬重現(xiàn)自然系統(tǒng)中動(dòng)荷載作用下的阻尼大小。

    FLAC3D的動(dòng)力計(jì)算中提供了3種阻尼形式,分別為滯后阻尼、局部阻尼及瑞利阻尼。實(shí)踐證明,瑞利阻尼因其理論與常規(guī)動(dòng)力分析方法類似,雖然采用瑞利阻尼時(shí)的計(jì)算時(shí)間步長(zhǎng)較小,但是計(jì)算得到的加速度響應(yīng)規(guī)律相較于前二者更加符合實(shí)際情況,故也是動(dòng)力計(jì)算中通常使用的阻尼形式。

    綜合考慮模型計(jì)算時(shí)間及其模擬的準(zhǔn)確性,經(jīng)過(guò)反復(fù)調(diào)試,本文選取最小臨界阻尼比為0.02,最小中心頻率為10Hz。

    4數(shù)值計(jì)算結(jié)果驗(yàn)證

    為了驗(yàn)證數(shù)值計(jì)算結(jié)果,依托某在建高速公路爆破施工工程現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù),實(shí)際爆破參數(shù)設(shè)計(jì)如下。

    (1)鉆孔直徑d=90mm。

    (2)臺(tái)階高度H=10m。

    (3)最小抵抗線W光=K光d=20x120mm=2400mm,式中W光為光面爆破最小抵抗線,K光為光面爆破計(jì)算系數(shù),一般取K=15~25,巖石越硬K取值越小。該工程取K=20。

    (4)炮孔超深h=0.5~1.5m,孔深大和巖石堅(jiān)硬完整者取大值,反之取小值。此工程取h=1.0m。

    (5)孔距光=mW=(0.5~0.8)x2.4=(1.2~1.92)m,式中m為裝藥密集系數(shù)。該工程取光=1.2。

    (6)炮孔長(zhǎng)度L=(H+h)/sin。當(dāng)坡度為1:0.75時(shí),L=13.8m;當(dāng)坡度為1:0.5時(shí),L=12.3m。

    在模型上布設(shè)相應(yīng)監(jiān)測(cè)點(diǎn)來(lái)跟蹤監(jiān)測(cè)其速度時(shí)程,位于K86+335斷面的5個(gè)監(jiān)測(cè)點(diǎn)位置,詳見(jiàn)圖6。因選取的數(shù)值計(jì)算荷載施加相關(guān)計(jì)算參數(shù)與現(xiàn)場(chǎng)監(jiān)測(cè)時(shí)的實(shí)際爆破設(shè)計(jì)裝藥參數(shù)相近,且該工況實(shí)際測(cè)點(diǎn)布置與數(shù)值模型中測(cè)線位于同一斷面,將二者峰值振動(dòng)速度進(jìn)行對(duì)比,得到數(shù)值計(jì)算各監(jiān)測(cè)點(diǎn)時(shí)程曲線如圖7所示,表3為數(shù)值計(jì)算各監(jiān)測(cè)點(diǎn)與現(xiàn)場(chǎng)監(jiān)測(cè)點(diǎn)峰值振動(dòng)速度對(duì)比結(jié)果,并將其繪成折線圖如圖8所示。

    由表3及圖8可知,數(shù)值模擬結(jié)果與現(xiàn)場(chǎng)監(jiān)測(cè)各測(cè)點(diǎn)峰值振動(dòng)速度在爆源近1~3號(hào)測(cè)點(diǎn)區(qū)其相對(duì)誤差均在10%以下,到了爆源遠(yuǎn)區(qū)誤差相對(duì)較大,且數(shù)值模擬峰值振速普遍大于現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù),可能原因是簡(jiǎn)化模型為均勻連續(xù)介質(zhì),忽略了實(shí)際巖體內(nèi)部節(jié)理裂隙等不連續(xù)結(jié)構(gòu)面,故而使爆破振動(dòng)地震波的衰減速度相較與實(shí)際工程略有減小??傮w而言,數(shù)值計(jì)算結(jié)果的精度尤其是在爆源近區(qū)的都在可接受范圍內(nèi),說(shuō)明借助FLAC3D研究該邊坡爆破振動(dòng)動(dòng)力響應(yīng)特征是可行的。

    5動(dòng)態(tài)響應(yīng)分析

    5.1速度響應(yīng)分析

    圖9為爆炸后5,10ms時(shí)模型X向速度分布云圖。由圖可知爆破地震波的作用影響范圍,因爆破振動(dòng)地震波具有方向性,特別是對(duì)于這種深孔臺(tái)階開(kāi)挖,隨著時(shí)間的推進(jìn),雖說(shuō)爆破荷載施加的反方向巖體也有一定的響應(yīng),但爆破振動(dòng)地震波主要朝著施加荷載的右邊坡傳播??偟膩?lái)說(shuō)在這種邊坡爆破中,地震波的傳播可以概化成如圖10所示的橢圓等振線。故可以得到的結(jié)論是,右邊坡的爆破開(kāi)挖對(duì)既有左邊坡影響可以忽略不計(jì)。

    圖11為爆炸后5,10ms時(shí)K86+335斷面X向速度分布云圖??汕逦赜^察到,隨著爆破地震波的傳播,在坡腳右側(cè)某一半徑內(nèi)質(zhì)點(diǎn)振動(dòng)速度在某一區(qū)域得到加強(qiáng),其原因是應(yīng)力波在邊坡自由面反射而產(chǎn)生了一定區(qū)域的振動(dòng)加強(qiáng)場(chǎng);又根據(jù)圖11(b)所顯示的各臺(tái)階周邊速度大小分布情況可知,該時(shí)刻速度值在第2級(jí)臺(tái)階坡腳處形成一個(gè)極大中心,順坡而上其速度值是先減小又增大,并在第4級(jí)臺(tái)階坡腳處再次形成一個(gè)極大中心,隨后減小。

    5.2位移響應(yīng)分析

    為進(jìn)一步研究邊坡在爆破振動(dòng)荷載作用下相關(guān)參數(shù)的響應(yīng)特征分布情況,引入無(wú)量綱位移放大系數(shù),其定義為邊坡爆破振動(dòng)位移響應(yīng)大小與第1級(jí)臺(tái)階邊緣位移大小的比值。

    根據(jù)數(shù)值計(jì)算結(jié)果,繪制出邊坡地表及K86+335斷面在10ms時(shí)Z向位移放大系數(shù)等值線圖如圖12所示。由圖12(a)可知,位移響應(yīng)特征與前述地震波橢圓等振線分布情況基本一致。圖12(b)中的位移放大系數(shù)等值線分布情況再次反映了右邊坡爆破開(kāi)挖未對(duì)左邊坡造成的影響。為進(jìn)一步分析邊坡坡體內(nèi)的位移響應(yīng)特征規(guī)律,將圖12(b)右邊坡局部放大顯示如圖13所示。

    由圖13可知,爆破荷載施加的1級(jí)邊坡其位移在一定深度范圍內(nèi)隨深度的增加而出現(xiàn)極大中心,超過(guò)某一深度后位移放大系數(shù)減小;再觀察位移放大系數(shù)等值線隨高程增加所表現(xiàn)出來(lái)的特征時(shí)發(fā)現(xiàn),其分布規(guī)律為位移放大系數(shù)在1級(jí)邊坡表面先減小,順坡面而上在進(jìn)入2級(jí)臺(tái)階時(shí)位移放大系數(shù)逐漸增大,隨后在2級(jí)邊坡中部位置順坡面而上位移放大系數(shù)開(kāi)始減小。此時(shí)位移放大系數(shù)所表現(xiàn)出的分布規(guī)律與速度響應(yīng)大小的分布規(guī)律近乎一致。

    5.3應(yīng)力響應(yīng)分析

    圖14為K86+335斷面炸藥爆炸后不同時(shí)刻最大主應(yīng)力分布云圖。

    由圖可知,各個(gè)時(shí)刻除爆源附近一定區(qū)域內(nèi)應(yīng)力在進(jìn)一步不斷地調(diào)整,邊坡其他區(qū)域等值線變化不大,其原因是所施加的荷載相較于自重荷載比較小,尚不足以使整個(gè)邊坡的應(yīng)力重分布。因爆破地震波主要以壓縮波在坡體內(nèi)傳播,爆破荷載作用下邊坡受影響區(qū)域多出現(xiàn)正常的壓應(yīng)力,但在局部如邊坡坡腳附近、個(gè)別坡面及臺(tái)階邊緣,會(huì)出現(xiàn)一定區(qū)域的拉應(yīng)力,而出現(xiàn)拉應(yīng)力的區(qū)域基本上都是接近坡面,即自由面,其可能原因是爆破振動(dòng)應(yīng)力波在自由面發(fā)生發(fā)射,進(jìn)而生成一定的拉伸波造成相應(yīng)區(qū)域出現(xiàn)拉應(yīng)力。

    由圖14可知,爆后隨著應(yīng)力的逐步調(diào)整,無(wú)論是拉應(yīng)力還是壓應(yīng)力,最大值均出現(xiàn)在坡腳,即靠近爆源附近,且相應(yīng)的應(yīng)力值均小于巖石相應(yīng)的臨界強(qiáng)度值,在實(shí)際的邊坡開(kāi)挖爆破中,也是坡腳處最容易率先發(fā)生局部的失穩(wěn)破壞。又由6~10ms的應(yīng)力分布云圖可知,壓應(yīng)力主要集中區(qū)域在該時(shí)段內(nèi)逐漸向深處移動(dòng),而因應(yīng)力波逐漸傳播至坡面,同時(shí)在坡面發(fā)生反射,生成一定的拉伸波,該時(shí)段坡面逐漸開(kāi)始產(chǎn)生拉應(yīng)力,且逐漸地在坡面附近產(chǎn)生拉應(yīng)力集中區(qū)域;10ms過(guò)后,坡面周邊拉應(yīng)力集中區(qū)域逐漸消散,僅在邊坡坡腳及1,2級(jí)臺(tái)階突出物保持著較小拉應(yīng)力集中區(qū)。同時(shí),在10~1000ms,受爆破動(dòng)荷載影響,較大的坡腳及1,2級(jí)臺(tái)階附近的應(yīng)力狀態(tài)幾乎不發(fā)生變化,應(yīng)力調(diào)整基本趨于穩(wěn)定,其余部分逐漸向初始狀態(tài)下的應(yīng)力分布進(jìn)行調(diào)整。

    6結(jié)論

    (1)邊坡開(kāi)挖爆破中,地震波的傳播可以概化為橢圓等振線,地震波主要朝施加荷載的正向傳播,對(duì)既有左邊坡的影響較小。

    (2)在巖質(zhì)均一的邊坡中,無(wú)論是速度還是位移,其爆破振動(dòng)的響應(yīng)規(guī)律較為一致,且在坡面和坡腳處因應(yīng)力波在邊坡自由面反射容易產(chǎn)生一定區(qū)域的振動(dòng)加強(qiáng)場(chǎng)。

    (3)在巖質(zhì)邊坡爆破作用影響下,無(wú)論是拉應(yīng)力還是壓應(yīng)力其最大值均出現(xiàn)在坡腳附近。

    (4)在實(shí)際的邊坡爆破開(kāi)挖過(guò)程中,對(duì)邊坡的坡面,尤其是坡腳處應(yīng)該加強(qiáng)保護(hù)措施。

    參考文獻(xiàn):

    [1]Aysilio E,Conte E,Dente G.Seismic stability analysis of reinforcedslopes[J].Soil dynamics and earthquake,2000,19(3):159-172.

    [2]Spyros Sklavounos,F(xiàn)otisRigas.Computer simulation of shock wavestransmission in obstructed Terrains[J].Journal of Loss Prevention inthe Process Industries,2004,(17):407-417.

    [3]Siyahi,Bilge G.Pseudo-static stability analysis innormally consolidated soil slopes subjected to earthquake[J].Teknik Dergi/Technicaljournal of Turksih chamber of civil enginers,1998,9(DEC):457-461.

    [4]UDEC(Universal Distinct Element Code)user' s manual Version 3.0[M].Lundon:Itasca Consulting Group,Inc.,1996.

    [5]蔣楠,周傳波,平雯,等.巖質(zhì)邊坡爆破振動(dòng)高速度高程效應(yīng)[J].中南大學(xué)學(xué)報(bào):自然科學(xué)版,2014,45(1):237-243.

    [6]唐海,李海波,蔣鵬燦,等.地形地貌對(duì)爆破振動(dòng)波傳播的影響實(shí)驗(yàn)研究[J].巖石學(xué)與工程學(xué)報(bào),2007(9):1817-1823.

    [7]唐海,李海波.地形地貌對(duì)爆破振動(dòng)波傳播的影響實(shí)驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào),2007(09):1817-1823

    [8]言志信,蔡漢成,王群敏,等.縱波作用下邊坡動(dòng)力響應(yīng)[J].工程地質(zhì)學(xué)報(bào),2010,18(4):487-505.

    [9]祁生文,伍法權(quán).邊坡動(dòng)力響應(yīng)規(guī)律研究[J].中國(guó)科學(xué)E輯技術(shù)科學(xué),2003,33(s):28-40.

    [10]何錚,徐衛(wèi)亞,石崇,等.順層巖質(zhì)高邊坡地震變形破壞機(jī)制三維數(shù)值反演研究[J].巖土力學(xué),2009,30(11):3512-3518.

    [11]陳明,盧文波,李鵬,等.巖質(zhì)邊坡爆破振動(dòng)速度的高程放大效應(yīng)研究[J].巖石力學(xué)與工程學(xué)報(bào),2011(11):2189-2195.

    [12]李鵬,蘇生瑞,王閆超,等.含軟弱層巖質(zhì)邊坡的動(dòng)力響應(yīng)研究[J].巖土力學(xué),2013(S1):365-370,378.

    [13]孟祥甜.露天邊坡動(dòng)載響應(yīng)特性的數(shù)值模擬研究[D].包頭:內(nèi)蒙古科技大學(xué),2015.

    [14]宋光明,陳壽如,史秀志,等.露天礦邊坡爆破振動(dòng)監(jiān)測(cè)與評(píng)價(jià)方法的研究[J].有色金屬,2000(4):24-27.

    [15]王在泉,陸文興.高邊坡爆破開(kāi)挖振動(dòng)傳統(tǒng)規(guī)律及質(zhì)量控制[J].爆破,1994,11(3):1-4.

    [16]夏詳,李俊如,李海波,等.爆破荷載作用下巖體振動(dòng)特征的數(shù)值模擬[J].巖土力學(xué),2005,26(1):50-56.

    [17]Kyhlemeyer R L,Lysmer J.Finite Element Method Accuracy forWave Propagation Problems[J].Soil Mech & Foundations Div,ASCE,1973,99(5):417-421.

    [18]Josef Henrych.The Dynamics of Explosion and Its Use[M].NewYork:Elsevier Scientific Publishing Company,1979.

    [19]S.G.Chen,J.Zhao.UDEC modeling of a field explosion test[J].International Journal of Blasting and Fragmentation,2000(4):149-163

    [20]Josef Henrych.The Dynamics of Explosion and Its Use[M].NewYork:Elsevier Scientific Publishing Company,1979.

    引用本文:唐旭,方正峰,鄒飛.基于FLAC3的巖質(zhì)邊坡爆破動(dòng)力響應(yīng)規(guī)律研究[J].人民長(zhǎng)江,2019,50(3):198-204.

    Study on dynamic response laws of rock slope blasting based on FLAC3D

    TANG Xu',F(xiàn)ANG Zhengfeng2,ZOU Fei,4

    (1.Guizhou Province Quality and Safety Traffic Engineering Monitoring and Inspection Center Co.,Ltd.,Guiyang 550081,China;2.Guizhou Hongxin Chuangda Engineering Detection & Consultation Co,Ltd.,Guiyang 550014,China;3.China GuizhouTransportation Planning Survey & Design Academe,Guiyang 550081,China;4.School of Ciwil Engineering,Central South University,Changsha 410000,China)

    Abstract:The dynamic stability of slope has always been the focus of geotechnical engineering and seismic engineering.Inthis paper,F(xiàn)LACED numerical calculation method is used to simulate the vibration velocity,displacement and stress of rocky slopeunder blasting action.And on the basis of comparing and verifying the peak value of blasting vibration,the dynamic responselaws of rocky slope was analyzed.The results show that the seismic wave has a strong directionality and propagates mainly to thepositive direction of the applied load.On the slope surface and at the foot of the slope,it is easy to produce the vibrationstrengthening field in a certain area due to the reflection of the force wave on the free surface of the slope.The maximum value ofstress,whether tensile stress or compressive stress,appears near the slope foot.In the actual blasting excavation,the slope surface should be protected,especially the slope foot.The research results have certain reference for the construction design of rockyslope excavation.

    Key words:dynamic response;rock slope;FLAC3D;blasting

    猜你喜歡
    爆破
    緊鄰次高壓燃?xì)夤芫€小間距隧道控制爆破施工技術(shù)及安全防護(hù)
    淺談水利工程施工中高邊坡支護(hù)與開(kāi)挖技術(shù)的應(yīng)用
    安哥拉道碴廠露天爆破施工技術(shù)
    在松散巖土層條件下如何進(jìn)行快速掘進(jìn)爆破技術(shù)探究
    地面減震孔在地鐵側(cè)穿加油站爆破施工中的應(yīng)用
    淺孔臺(tái)階靜力爆破技術(shù)在基坑開(kāi)挖中的應(yīng)用
    1024视频免费在线观看| 在线天堂最新版资源| 天堂中文最新版在线下载| 91aial.com中文字幕在线观看| 国产女主播在线喷水免费视频网站| 亚洲精品日韩在线中文字幕| 男女国产视频网站| 99热全是精品| 考比视频在线观看| 天天躁夜夜躁狠狠久久av| 色吧在线观看| 好男人视频免费观看在线| 亚洲五月色婷婷综合| 欧美97在线视频| 亚洲人成77777在线视频| 亚洲美女黄色视频免费看| 久久久久精品国产欧美久久久 | 欧美激情极品国产一区二区三区| 麻豆av在线久日| 亚洲av福利一区| 一级毛片我不卡| 男人舔女人的私密视频| 午夜福利免费观看在线| 亚洲欧美成人综合另类久久久| 日本猛色少妇xxxxx猛交久久| 我的亚洲天堂| 免费看av在线观看网站| 捣出白浆h1v1| 99香蕉大伊视频| 色视频在线一区二区三区| 亚洲五月色婷婷综合| 两个人免费观看高清视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲美女黄色视频免费看| 欧美另类一区| 一边摸一边抽搐一进一出视频| 九九爱精品视频在线观看| 国产又色又爽无遮挡免| 大香蕉久久成人网| av在线播放精品| 侵犯人妻中文字幕一二三四区| 九九爱精品视频在线观看| 九草在线视频观看| 欧美日韩福利视频一区二区| av国产精品久久久久影院| 妹子高潮喷水视频| 久久久久精品人妻al黑| 亚洲国产av新网站| 最近最新中文字幕免费大全7| 18禁裸乳无遮挡动漫免费视频| 美女福利国产在线| 丁香六月天网| 狠狠婷婷综合久久久久久88av| 黄色毛片三级朝国网站| 免费久久久久久久精品成人欧美视频| 王馨瑶露胸无遮挡在线观看| 韩国精品一区二区三区| videosex国产| 成年女人毛片免费观看观看9 | 亚洲,一卡二卡三卡| 高清视频免费观看一区二区| av免费观看日本| 两个人看的免费小视频| 亚洲一区中文字幕在线| 久久毛片免费看一区二区三区| 天天添夜夜摸| 一本色道久久久久久精品综合| 久久热在线av| 超碰97精品在线观看| 赤兔流量卡办理| 电影成人av| 国产亚洲午夜精品一区二区久久| 9热在线视频观看99| 在线观看免费视频网站a站| 交换朋友夫妻互换小说| av天堂久久9| 国产亚洲av高清不卡| 久久亚洲国产成人精品v| 丝袜人妻中文字幕| 久久国产亚洲av麻豆专区| 亚洲精品国产av蜜桃| 国产精品.久久久| 中文字幕人妻熟女乱码| 老司机深夜福利视频在线观看 | 一级毛片我不卡| 无遮挡黄片免费观看| 51午夜福利影视在线观看| 国产精品人妻久久久影院| 99久国产av精品国产电影| 国产一区有黄有色的免费视频| 九草在线视频观看| 性少妇av在线| 欧美日韩av久久| 考比视频在线观看| 少妇被粗大猛烈的视频| 2021少妇久久久久久久久久久| 丁香六月天网| 国产成人精品久久二区二区91 | 午夜精品国产一区二区电影| 人人妻人人澡人人爽人人夜夜| 国产黄频视频在线观看| 一边亲一边摸免费视频| 激情视频va一区二区三区| 免费看不卡的av| a级毛片在线看网站| 久久久久人妻精品一区果冻| 久久久久国产精品人妻一区二区| 18禁国产床啪视频网站| 日韩伦理黄色片| 99久久99久久久精品蜜桃| 在线观看人妻少妇| 国产精品 国内视频| 久久ye,这里只有精品| 一区二区三区激情视频| 国产免费现黄频在线看| 美女扒开内裤让男人捅视频| 国产乱来视频区| 亚洲久久久国产精品| 亚洲天堂av无毛| 美女福利国产在线| 日韩精品免费视频一区二区三区| 51午夜福利影视在线观看| 2018国产大陆天天弄谢| 国产精品久久久久久精品古装| 国产精品久久久人人做人人爽| 男女下面插进去视频免费观看| 在线观看三级黄色| 777久久人妻少妇嫩草av网站| 国产xxxxx性猛交| 亚洲精品aⅴ在线观看| 美女扒开内裤让男人捅视频| 性色av一级| 国产成人精品福利久久| 亚洲国产欧美在线一区| 水蜜桃什么品种好| 亚洲成av片中文字幕在线观看| 婷婷色麻豆天堂久久| 欧美少妇被猛烈插入视频| 精品国产超薄肉色丝袜足j| 人妻一区二区av| 观看av在线不卡| 午夜91福利影院| 日本欧美视频一区| 日韩av不卡免费在线播放| 丰满迷人的少妇在线观看| 丝袜在线中文字幕| 老司机影院毛片| 一本色道久久久久久精品综合| 制服人妻中文乱码| 韩国精品一区二区三区| 交换朋友夫妻互换小说| av国产精品久久久久影院| 纯流量卡能插随身wifi吗| 国产男女超爽视频在线观看| 夫妻午夜视频| 婷婷成人精品国产| 你懂的网址亚洲精品在线观看| 十分钟在线观看高清视频www| av电影中文网址| 日韩精品免费视频一区二区三区| 亚洲人成电影观看| 99九九在线精品视频| 丁香六月天网| 日韩一本色道免费dvd| 久久久精品94久久精品| 国产乱人偷精品视频| 母亲3免费完整高清在线观看| av网站在线播放免费| 国产亚洲av高清不卡| 在线精品无人区一区二区三| 桃花免费在线播放| 性色av一级| 欧美日韩综合久久久久久| 欧美日韩亚洲国产一区二区在线观看 | 精品午夜福利在线看| 侵犯人妻中文字幕一二三四区| 熟女av电影| 亚洲伊人久久精品综合| 最近2019中文字幕mv第一页| 色94色欧美一区二区| 久久久国产一区二区| 亚洲精品在线美女| 国产伦人伦偷精品视频| 欧美久久黑人一区二区| 国产成人精品在线电影| 少妇的丰满在线观看| 国产精品嫩草影院av在线观看| 日韩一区二区视频免费看| 日本av免费视频播放| 亚洲综合精品二区| 精品人妻熟女毛片av久久网站| 最黄视频免费看| 免费在线观看视频国产中文字幕亚洲 | 色播在线永久视频| 高清视频免费观看一区二区| www.精华液| 美女视频免费永久观看网站| 老鸭窝网址在线观看| 人人澡人人妻人| 91成人精品电影| 大香蕉久久成人网| av女优亚洲男人天堂| 男女床上黄色一级片免费看| 国产麻豆69| 最近的中文字幕免费完整| www.av在线官网国产| 又粗又硬又长又爽又黄的视频| 97人妻天天添夜夜摸| 黄色一级大片看看| 精品第一国产精品| 国产精品香港三级国产av潘金莲 | 成人18禁高潮啪啪吃奶动态图| 女的被弄到高潮叫床怎么办| 一级,二级,三级黄色视频| 久久 成人 亚洲| 男女无遮挡免费网站观看| 欧美日韩一级在线毛片| 秋霞在线观看毛片| 老司机影院毛片| 久久久久精品国产欧美久久久 | 精品卡一卡二卡四卡免费| 日韩欧美一区视频在线观看| 日本wwww免费看| 国产精品99久久99久久久不卡 | 婷婷色综合www| 欧美黑人欧美精品刺激| 国产97色在线日韩免费| 18禁动态无遮挡网站| 男女免费视频国产| 天天躁日日躁夜夜躁夜夜| 18禁国产床啪视频网站| 欧美中文综合在线视频| 欧美av亚洲av综合av国产av | 九九爱精品视频在线观看| 一级片免费观看大全| 黄色毛片三级朝国网站| 国产精品二区激情视频| 欧美少妇被猛烈插入视频| 97在线人人人人妻| 婷婷色综合大香蕉| 韩国高清视频一区二区三区| 精品福利永久在线观看| 老司机深夜福利视频在线观看 | 国产av一区二区精品久久| 少妇人妻 视频| 精品少妇一区二区三区视频日本电影 | 99精品久久久久人妻精品| 久久久久网色| 99精国产麻豆久久婷婷| 啦啦啦视频在线资源免费观看| 在线观看一区二区三区激情| 日韩制服骚丝袜av| 如何舔出高潮| 两性夫妻黄色片| 亚洲国产精品成人久久小说| 狠狠婷婷综合久久久久久88av| 亚洲情色 制服丝袜| 国产一区有黄有色的免费视频| 青草久久国产| 一边摸一边做爽爽视频免费| 中文天堂在线官网| 又黄又粗又硬又大视频| 国产免费现黄频在线看| 国精品久久久久久国模美| 人妻 亚洲 视频| 少妇人妻久久综合中文| 亚洲精品国产av蜜桃| 日韩大码丰满熟妇| 丝袜脚勾引网站| 亚洲人成网站在线观看播放| 成人手机av| 另类亚洲欧美激情| 亚洲色图 男人天堂 中文字幕| 欧美中文综合在线视频| 久久久精品区二区三区| 深夜精品福利| 少妇人妻久久综合中文| 九色亚洲精品在线播放| 国产成人午夜福利电影在线观看| 伊人亚洲综合成人网| 啦啦啦在线免费观看视频4| 国产精品久久久人人做人人爽| 久久精品久久久久久噜噜老黄| 老司机在亚洲福利影院| 国产有黄有色有爽视频| 91精品国产国语对白视频| 久久av网站| 午夜久久久在线观看| 无遮挡黄片免费观看| 免费在线观看黄色视频的| 亚洲欧美激情在线| 精品免费久久久久久久清纯 | 女性被躁到高潮视频| 人人妻人人澡人人爽人人夜夜| 久久久久人妻精品一区果冻| 久久天堂一区二区三区四区| 激情五月婷婷亚洲| 国产免费又黄又爽又色| av又黄又爽大尺度在线免费看| 国产成人欧美| 如日韩欧美国产精品一区二区三区| 五月天丁香电影| 亚洲成人av在线免费| 天美传媒精品一区二区| 国产探花极品一区二区| 人人澡人人妻人| 国产又色又爽无遮挡免| 一边亲一边摸免费视频| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲国产日韩| 欧美成人精品欧美一级黄| 久久99精品国语久久久| 九色亚洲精品在线播放| 男女边摸边吃奶| a级毛片黄视频| 欧美日韩亚洲综合一区二区三区_| 精品一区二区三区av网在线观看 | 操美女的视频在线观看| 亚洲伊人色综图| 少妇被粗大猛烈的视频| 国产爽快片一区二区三区| 哪个播放器可以免费观看大片| 欧美精品亚洲一区二区| 日韩av不卡免费在线播放| 亚洲精品日本国产第一区| 精品国产乱码久久久久久男人| 国产av国产精品国产| 欧美老熟妇乱子伦牲交| 国产片内射在线| 午夜免费观看性视频| 少妇的丰满在线观看| videos熟女内射| 国产精品秋霞免费鲁丝片| 三上悠亚av全集在线观看| 亚洲美女黄色视频免费看| 亚洲欧美激情在线| 街头女战士在线观看网站| 99热全是精品| 青草久久国产| 丝袜人妻中文字幕| 欧美 亚洲 国产 日韩一| 亚洲成人av在线免费| 嫩草影院入口| 美女脱内裤让男人舔精品视频| av国产精品久久久久影院| 中文字幕人妻丝袜一区二区 | 女人久久www免费人成看片| 精品一区在线观看国产| 尾随美女入室| 综合色丁香网| 19禁男女啪啪无遮挡网站| xxxhd国产人妻xxx| 岛国毛片在线播放| 大片免费播放器 马上看| 亚洲av日韩精品久久久久久密 | 中文字幕精品免费在线观看视频| 日韩欧美免费精品| 久久精品亚洲熟妇少妇任你| 国产伦人伦偷精品视频| 欧美不卡视频在线免费观看 | 亚洲久久久国产精品| 黑人欧美特级aaaaaa片| 99在线人妻在线中文字幕| 免费在线观看日本一区| 久久久久久国产a免费观看| 国产精品野战在线观看| 嫩草影院精品99| 91九色精品人成在线观看| 国产精品爽爽va在线观看网站 | 欧美精品啪啪一区二区三区| 啪啪无遮挡十八禁网站| 欧美不卡视频在线免费观看 | 国产精品自产拍在线观看55亚洲| 麻豆一二三区av精品| 国产三级黄色录像| 99久久99久久久精品蜜桃| 亚洲性夜色夜夜综合| 精品久久久久久久人妻蜜臀av | 最新美女视频免费是黄的| 久久中文看片网| 欧美久久黑人一区二区| 一二三四社区在线视频社区8| 久久国产精品人妻蜜桃| 啦啦啦 在线观看视频| 国产精品野战在线观看| 极品教师在线免费播放| www日本在线高清视频| 亚洲av第一区精品v没综合| 中文字幕精品免费在线观看视频| 男女之事视频高清在线观看| 一边摸一边抽搐一进一出视频| а√天堂www在线а√下载| 免费在线观看视频国产中文字幕亚洲| 国产成人精品在线电影| 在线观看舔阴道视频| av在线天堂中文字幕| 又黄又粗又硬又大视频| av网站免费在线观看视频| 久热爱精品视频在线9| 午夜免费激情av| 999精品在线视频| 久9热在线精品视频| 男人舔女人的私密视频| 女人被狂操c到高潮| 97人妻精品一区二区三区麻豆 | 欧美久久黑人一区二区| 国产一区二区激情短视频| 国产成人精品久久二区二区免费| 一个人免费在线观看的高清视频| 久久久久久久久中文| 亚洲狠狠婷婷综合久久图片| www日本在线高清视频| 免费搜索国产男女视频| 在线观看www视频免费| 一级作爱视频免费观看| 涩涩av久久男人的天堂| 精品久久久久久成人av| 88av欧美| 99国产精品99久久久久| 精品免费久久久久久久清纯| 侵犯人妻中文字幕一二三四区| 1024视频免费在线观看| 精品少妇一区二区三区视频日本电影| 桃红色精品国产亚洲av| 天天躁夜夜躁狠狠躁躁| 国产单亲对白刺激| 日韩欧美一区二区三区在线观看| 最近最新中文字幕大全电影3 | 成在线人永久免费视频| 欧美激情久久久久久爽电影 | 久久婷婷成人综合色麻豆| 一进一出抽搐动态| 搡老岳熟女国产| 色尼玛亚洲综合影院| 亚洲男人天堂网一区| 久久 成人 亚洲| 亚洲九九香蕉| 老司机在亚洲福利影院| 后天国语完整版免费观看| 成人国产一区最新在线观看| 久久久久国产一级毛片高清牌| 免费少妇av软件| 国产片内射在线| 99国产精品免费福利视频| 日日摸夜夜添夜夜添小说| 99国产综合亚洲精品| 激情在线观看视频在线高清| 欧美乱妇无乱码| 日韩精品青青久久久久久| 老司机午夜福利在线观看视频| 中文字幕人成人乱码亚洲影| 国产精品一区二区在线不卡| 变态另类丝袜制服| 国产伦一二天堂av在线观看| 夜夜爽天天搞| 12—13女人毛片做爰片一| 久久精品亚洲精品国产色婷小说| 免费女性裸体啪啪无遮挡网站| 中文字幕精品免费在线观看视频| 久久久久久久精品吃奶| 不卡一级毛片| 国产97色在线日韩免费| 他把我摸到了高潮在线观看| 久久国产亚洲av麻豆专区| 日韩欧美免费精品| 在线观看免费视频日本深夜| 黄色片一级片一级黄色片| 亚洲成人精品中文字幕电影| av天堂久久9| 亚洲中文av在线| 又黄又粗又硬又大视频| 午夜免费激情av| 91字幕亚洲| 91麻豆av在线| 中文字幕人妻丝袜一区二区| 香蕉久久夜色| 一边摸一边抽搐一进一小说| 久久人人97超碰香蕉20202| 中出人妻视频一区二区| 午夜福利,免费看| 亚洲av成人一区二区三| 亚洲五月天丁香| 少妇粗大呻吟视频| 国产成人欧美在线观看| 黄片小视频在线播放| 在线国产一区二区在线| 日本 欧美在线| 夜夜看夜夜爽夜夜摸| 一进一出抽搐动态| 久久国产乱子伦精品免费另类| av福利片在线| 亚洲全国av大片| 成人三级做爰电影| 在线永久观看黄色视频| 国产片内射在线| 一进一出好大好爽视频| 色综合婷婷激情| cao死你这个sao货| 老熟妇仑乱视频hdxx| 欧美一级毛片孕妇| 亚洲精品美女久久久久99蜜臀| 在线国产一区二区在线| 日本 欧美在线| 午夜精品国产一区二区电影| 成人亚洲精品av一区二区| 99国产精品免费福利视频| 国产精品亚洲美女久久久| 国产成人影院久久av| 亚洲 国产 在线| 亚洲男人的天堂狠狠| 亚洲精品美女久久久久99蜜臀| netflix在线观看网站| 国产精品 国内视频| 欧美一区二区精品小视频在线| 亚洲免费av在线视频| 18禁国产床啪视频网站| 咕卡用的链子| 久久久久国内视频| 日本五十路高清| 级片在线观看| 免费高清在线观看日韩| 国产亚洲av高清不卡| 日本vs欧美在线观看视频| 精品国产美女av久久久久小说| 看免费av毛片| 两个人视频免费观看高清| 老汉色∧v一级毛片| 女人被狂操c到高潮| 亚洲男人天堂网一区| 无限看片的www在线观看| 国产亚洲av嫩草精品影院| 免费看a级黄色片| 日本在线视频免费播放| 女性生殖器流出的白浆| 亚洲欧美精品综合一区二区三区| 中文字幕人妻熟女乱码| 日韩精品免费视频一区二区三区| 日韩大尺度精品在线看网址 | 在线观看免费视频日本深夜| 免费不卡黄色视频| 亚洲av美国av| 欧美性长视频在线观看| 午夜精品在线福利| 欧美在线黄色| 久久国产亚洲av麻豆专区| 一个人观看的视频www高清免费观看 | 日韩欧美三级三区| 一本大道久久a久久精品| 久久性视频一级片| 国产xxxxx性猛交| 亚洲精品av麻豆狂野| 大型av网站在线播放| 99久久久亚洲精品蜜臀av| 亚洲精品国产一区二区精华液| 中文字幕人妻丝袜一区二区| 亚洲国产精品合色在线| bbb黄色大片| 他把我摸到了高潮在线观看| 女人爽到高潮嗷嗷叫在线视频| 午夜福利在线观看吧| 亚洲伊人色综图| 国产aⅴ精品一区二区三区波| 亚洲av电影在线进入| 黑人巨大精品欧美一区二区mp4| 日本一区二区免费在线视频| 国产精品久久久人人做人人爽| 国产成人免费无遮挡视频| 母亲3免费完整高清在线观看| 成人特级黄色片久久久久久久| 精品国产一区二区久久| 1024香蕉在线观看| 超碰成人久久| 免费在线观看黄色视频的| 色综合亚洲欧美另类图片| 一本综合久久免费| 欧美人与性动交α欧美精品济南到| 搡老岳熟女国产| 一级毛片精品| www.www免费av| 精品久久久久久,| 久久人人爽av亚洲精品天堂| 亚洲av第一区精品v没综合| ponron亚洲| 级片在线观看| 一夜夜www| 大型黄色视频在线免费观看| 久久午夜综合久久蜜桃| 亚洲av第一区精品v没综合| 精品电影一区二区在线| 亚洲天堂国产精品一区在线| 18禁美女被吸乳视频| 91字幕亚洲| 成人18禁高潮啪啪吃奶动态图| 黑丝袜美女国产一区| 999久久久精品免费观看国产| 老司机午夜福利在线观看视频| 国产精品二区激情视频| www日本在线高清视频| 三级毛片av免费| 99热只有精品国产| 手机成人av网站| 亚洲欧美激情综合另类| 十八禁网站免费在线| 一二三四在线观看免费中文在| 热99re8久久精品国产| 午夜福利,免费看| 成人国产综合亚洲| 日本vs欧美在线观看视频| 国产欧美日韩精品亚洲av| 久热爱精品视频在线9| 可以免费在线观看a视频的电影网站| 亚洲av成人av| 麻豆成人av在线观看| 精品午夜福利视频在线观看一区| 亚洲av成人av|