• <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)用
    交换朋友夫妻互换小说| 国产xxxxx性猛交| 我的亚洲天堂| 如日韩欧美国产精品一区二区三区| 日本精品一区二区三区蜜桃| 亚洲成人免费电影在线观看| 欧美精品av麻豆av| 美女大奶头黄色视频| 男人添女人高潮全过程视频| 免费一级毛片在线播放高清视频 | 精品亚洲乱码少妇综合久久| 欧美日韩亚洲国产一区二区在线观看 | 成人国语在线视频| 人妻一区二区av| 国产精品久久久久久精品电影小说| 精品国产一区二区久久| 国产在线视频一区二区| 亚洲熟女毛片儿| 亚洲伊人久久精品综合| 在线观看人妻少妇| 一本色道久久久久久精品综合| 99精品欧美一区二区三区四区| 777久久人妻少妇嫩草av网站| 日韩中文字幕视频在线看片| 国产在线视频一区二区| 超碰97精品在线观看| 国产成人av教育| videosex国产| 午夜福利乱码中文字幕| 一区在线观看完整版| 丝袜人妻中文字幕| 五月天丁香电影| 成人18禁高潮啪啪吃奶动态图| 国产一区二区三区在线臀色熟女 | 最近最新中文字幕大全免费视频| 成年av动漫网址| 亚洲精品av麻豆狂野| 国产成人影院久久av| 美女高潮喷水抽搐中文字幕| 久久久久久免费高清国产稀缺| 精品一区二区三区av网在线观看 | 亚洲欧美激情在线| 一级毛片电影观看| 超色免费av| 99精国产麻豆久久婷婷| 久久久欧美国产精品| 老司机靠b影院| 天堂8中文在线网| 视频区图区小说| a级毛片在线看网站| 人人妻人人爽人人添夜夜欢视频| 亚洲三区欧美一区| 国产片内射在线| 成人手机av| 亚洲九九香蕉| 免费黄频网站在线观看国产| 大片免费播放器 马上看| 操美女的视频在线观看| 这个男人来自地球电影免费观看| 国产在线观看jvid| netflix在线观看网站| 可以免费在线观看a视频的电影网站| 91av网站免费观看| 男女午夜视频在线观看| 不卡av一区二区三区| 丰满迷人的少妇在线观看| 亚洲美女黄色视频免费看| 亚洲久久久国产精品| 亚洲精品国产av蜜桃| 99国产精品一区二区三区| 国产男人的电影天堂91| 一区在线观看完整版| 在线精品无人区一区二区三| 精品高清国产在线一区| 一区二区av电影网| av在线app专区| 我要看黄色一级片免费的| 国产成人免费观看mmmm| 麻豆av在线久日| 久久久久精品人妻al黑| 免费黄频网站在线观看国产| av免费在线观看网站| 丰满迷人的少妇在线观看| 少妇的丰满在线观看| 欧美日韩视频精品一区| 久久国产精品大桥未久av| 国产淫语在线视频| 久久久久久人人人人人| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久精品人妻al黑| 亚洲 国产 在线| 窝窝影院91人妻| 欧美另类亚洲清纯唯美| 久久精品亚洲av国产电影网| 亚洲视频免费观看视频| 午夜两性在线视频| 丝袜美足系列| 男女高潮啪啪啪动态图| 热99re8久久精品国产| 免费女性裸体啪啪无遮挡网站| 男女无遮挡免费网站观看| av免费在线观看网站| 国产精品免费大片| 亚洲人成电影观看| 岛国毛片在线播放| 午夜福利视频在线观看免费| 亚洲精品第二区| 国产精品一区二区在线不卡| 麻豆乱淫一区二区| 我的亚洲天堂| 黄色怎么调成土黄色| 18禁黄网站禁片午夜丰满| 成人亚洲精品一区在线观看| 成年动漫av网址| 丝瓜视频免费看黄片| 欧美成人午夜精品| 精品少妇内射三级| 亚洲精品国产色婷婷电影| 午夜精品久久久久久毛片777| 男女之事视频高清在线观看| 美女高潮到喷水免费观看| 1024香蕉在线观看| 两个人免费观看高清视频| 国产真人三级小视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 精品少妇一区二区三区视频日本电影| 精品卡一卡二卡四卡免费| 高潮久久久久久久久久久不卡| 91大片在线观看| 日本一区二区免费在线视频| 久久九九热精品免费| 免费看十八禁软件| 91av网站免费观看| 成人亚洲精品一区在线观看| 日韩熟女老妇一区二区性免费视频| 国产一卡二卡三卡精品| 精品人妻在线不人妻| 成人国产一区最新在线观看| 不卡一级毛片| 亚洲熟女毛片儿| 一边摸一边抽搐一进一出视频| 人妻人人澡人人爽人人| 在线观看免费高清a一片| 国产精品一区二区在线不卡| svipshipincom国产片| 中文精品一卡2卡3卡4更新| 亚洲国产欧美在线一区| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕人妻丝袜制服| 90打野战视频偷拍视频| 久久久久精品人妻al黑| kizo精华| 国产精品国产av在线观看| tube8黄色片| 少妇精品久久久久久久| 亚洲av成人不卡在线观看播放网 | 亚洲性夜色夜夜综合| 国产91精品成人一区二区三区 | 日韩 亚洲 欧美在线| 一本大道久久a久久精品| 午夜91福利影院| 免费在线观看日本一区| 亚洲 欧美一区二区三区| 成人影院久久| 一本—道久久a久久精品蜜桃钙片| 菩萨蛮人人尽说江南好唐韦庄| av福利片在线| 成人手机av| 国产欧美日韩一区二区精品| 亚洲av日韩在线播放| 满18在线观看网站| 亚洲黑人精品在线| 啪啪无遮挡十八禁网站| 国内毛片毛片毛片毛片毛片| 男女床上黄色一级片免费看| 丁香六月天网| 免费在线观看影片大全网站| 美女国产高潮福利片在线看| 欧美 亚洲 国产 日韩一| 啦啦啦在线免费观看视频4| e午夜精品久久久久久久| 国产精品偷伦视频观看了| 男女下面插进去视频免费观看| 免费观看人在逋| 日本精品一区二区三区蜜桃| a级片在线免费高清观看视频| 午夜免费观看性视频| 国产又爽黄色视频| 巨乳人妻的诱惑在线观看| 嫩草影视91久久| 丝袜在线中文字幕| 动漫黄色视频在线观看| a级毛片黄视频| 男男h啪啪无遮挡| 黄色怎么调成土黄色| 国产激情久久老熟女| av网站在线播放免费| 爱豆传媒免费全集在线观看| 老司机午夜十八禁免费视频| 免费不卡黄色视频| 国产成人欧美在线观看 | 国产三级黄色录像| 久久狼人影院| 久久久精品免费免费高清| videos熟女内射| 欧美在线黄色| 久热这里只有精品99| 欧美97在线视频| 搡老岳熟女国产| 日韩欧美一区视频在线观看| 久久中文看片网| 18禁国产床啪视频网站| 亚洲精品一二三| 热99久久久久精品小说推荐| av视频免费观看在线观看| 1024香蕉在线观看| www.av在线官网国产| 日韩欧美免费精品| 老鸭窝网址在线观看| 欧美人与性动交α欧美精品济南到| av片东京热男人的天堂| 欧美另类亚洲清纯唯美| 国产精品 欧美亚洲| 黄片播放在线免费| 国产成人精品久久二区二区免费| 九色亚洲精品在线播放| 亚洲欧美精品自产自拍| 亚洲久久久国产精品| 久久久精品免费免费高清| 午夜激情av网站| a级毛片在线看网站| 国产精品二区激情视频| 国产老妇伦熟女老妇高清| 涩涩av久久男人的天堂| 日韩,欧美,国产一区二区三区| 国产成人精品无人区| 免费高清在线观看日韩| 中文欧美无线码| 免费av中文字幕在线| 久久久久网色| 一个人免费看片子| 成人亚洲精品一区在线观看| 国产国语露脸激情在线看| 亚洲精品av麻豆狂野| 悠悠久久av| 精品乱码久久久久久99久播| 桃红色精品国产亚洲av| 久久影院123| 国产视频一区二区在线看| 他把我摸到了高潮在线观看 | av天堂久久9| 国产欧美日韩一区二区三区在线| 1024视频免费在线观看| netflix在线观看网站| kizo精华| 中文精品一卡2卡3卡4更新| 超碰97精品在线观看| 亚洲精品中文字幕一二三四区 | 亚洲天堂av无毛| 亚洲成人免费av在线播放| 热99久久久久精品小说推荐| 五月开心婷婷网| 久久精品成人免费网站| 三级毛片av免费| 午夜日韩欧美国产| www.自偷自拍.com| 99精国产麻豆久久婷婷| 亚洲免费av在线视频| 搡老熟女国产l中国老女人| 天天躁夜夜躁狠狠躁躁| 久久中文字幕一级| 在线av久久热| 在线十欧美十亚洲十日本专区| 热99久久久久精品小说推荐| 久久九九热精品免费| 国产99久久九九免费精品| 伊人久久大香线蕉亚洲五| 精品国产一区二区三区久久久樱花| 老司机午夜十八禁免费视频| 国产精品一二三区在线看| cao死你这个sao货| 午夜福利影视在线免费观看| 91成年电影在线观看| 水蜜桃什么品种好| 久久女婷五月综合色啪小说| 久久99热这里只频精品6学生| 久热这里只有精品99| 汤姆久久久久久久影院中文字幕| 欧美精品啪啪一区二区三区 | 久久性视频一级片| 国产精品一区二区精品视频观看| 亚洲精品中文字幕一二三四区 | 丝袜在线中文字幕| 大香蕉久久成人网| 日韩 亚洲 欧美在线| 香蕉丝袜av| 亚洲精品久久午夜乱码| 黄片小视频在线播放| 久久久久国产一级毛片高清牌| 亚洲专区中文字幕在线| 久久九九热精品免费| 少妇 在线观看| 国内毛片毛片毛片毛片毛片| 9191精品国产免费久久| 十八禁网站免费在线| 欧美另类一区| 天天影视国产精品| 国产成人欧美在线观看 | 欧美成人午夜精品| 中文字幕另类日韩欧美亚洲嫩草| 丝袜美足系列| 大香蕉久久网| 日本av手机在线免费观看| 午夜福利一区二区在线看| 色婷婷av一区二区三区视频| 国产日韩欧美亚洲二区| 午夜精品久久久久久毛片777| 黄色视频在线播放观看不卡| 男女免费视频国产| 亚洲av片天天在线观看| 大香蕉久久网| 一级毛片精品| 十八禁网站网址无遮挡| 丰满迷人的少妇在线观看| 国产av一区二区精品久久| 十分钟在线观看高清视频www| 黑人巨大精品欧美一区二区蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 精品国产超薄肉色丝袜足j| av超薄肉色丝袜交足视频| 岛国在线观看网站| 天堂中文最新版在线下载| 亚洲国产欧美网| 91精品三级在线观看| 正在播放国产对白刺激| 久久天堂一区二区三区四区| 国产高清videossex| 亚洲国产毛片av蜜桃av| 国产成人欧美在线观看 | 婷婷色av中文字幕| 成年人黄色毛片网站| 老司机在亚洲福利影院| 亚洲欧美精品自产自拍| 一边摸一边抽搐一进一出视频| 大陆偷拍与自拍| 在线观看免费高清a一片| 深夜精品福利| 我的亚洲天堂| 精品人妻熟女毛片av久久网站| 男女高潮啪啪啪动态图| 亚洲国产看品久久| 亚洲精品乱久久久久久| 青春草视频在线免费观看| 国产区一区二久久| 久久久国产欧美日韩av| 一本色道久久久久久精品综合| 成人国产一区最新在线观看| 人人妻人人澡人人看| 婷婷成人精品国产| 国产熟女午夜一区二区三区| 欧美 日韩 精品 国产| tocl精华| 久久久久久久国产电影| 久久久国产成人免费| 亚洲精品在线美女| 亚洲天堂av无毛| 老司机靠b影院| 国产欧美日韩一区二区精品| 女人被躁到高潮嗷嗷叫费观| 亚洲熟女精品中文字幕| 国产高清视频在线播放一区 | 无遮挡黄片免费观看| 色婷婷av一区二区三区视频| 欧美av亚洲av综合av国产av| 亚洲免费av在线视频| 亚洲成人免费av在线播放| 国产成人免费观看mmmm| 欧美av亚洲av综合av国产av| 午夜激情久久久久久久| 搡老岳熟女国产| 各种免费的搞黄视频| 久久久国产欧美日韩av| 欧美一级毛片孕妇| 日韩精品免费视频一区二区三区| 亚洲精品国产色婷婷电影| 巨乳人妻的诱惑在线观看| 亚洲综合色网址| av在线老鸭窝| 久久久国产一区二区| 欧美黑人精品巨大| 搡老乐熟女国产| 久久天堂一区二区三区四区| 亚洲情色 制服丝袜| 人人妻人人爽人人添夜夜欢视频| 日日爽夜夜爽网站| 国产精品自产拍在线观看55亚洲 | 91老司机精品| 久久久水蜜桃国产精品网| 精品国内亚洲2022精品成人 | 自线自在国产av| 精品国产乱码久久久久久小说| 在线亚洲精品国产二区图片欧美| 一进一出抽搐动态| 新久久久久国产一级毛片| 两人在一起打扑克的视频| 亚洲成人免费av在线播放| 大码成人一级视频| 精品国产一区二区三区四区第35| 免费av中文字幕在线| 国产精品国产av在线观看| www.熟女人妻精品国产| 18禁国产床啪视频网站| 精品少妇久久久久久888优播| 捣出白浆h1v1| 99久久综合免费| 黄色视频,在线免费观看| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av高清一级| 国产老妇伦熟女老妇高清| 久久免费观看电影| 国产亚洲av高清不卡| 亚洲精品久久久久久婷婷小说| 日韩中文字幕欧美一区二区| 如日韩欧美国产精品一区二区三区| 999久久久精品免费观看国产| 欧美日韩福利视频一区二区| av线在线观看网站| 十分钟在线观看高清视频www| 日韩视频在线欧美| av一本久久久久| 91麻豆精品激情在线观看国产 | 日韩欧美免费精品| 五月天丁香电影| 久久久久久久久久久久大奶| 精品少妇一区二区三区视频日本电影| 国产av又大| 老司机靠b影院| 99久久国产精品久久久| 新久久久久国产一级毛片| 狂野欧美激情性xxxx| 啦啦啦免费观看视频1| √禁漫天堂资源中文www| 18禁国产床啪视频网站| 午夜福利视频精品| 中文字幕人妻丝袜制服| 一边摸一边抽搐一进一出视频| 男女边摸边吃奶| 亚洲三区欧美一区| 久久香蕉激情| 涩涩av久久男人的天堂| 国产激情久久老熟女| 80岁老熟妇乱子伦牲交| 亚洲欧美日韩高清在线视频 | 亚洲av男天堂| 美女午夜性视频免费| 国产成人啪精品午夜网站| av视频免费观看在线观看| 成人影院久久| 亚洲精品国产精品久久久不卡| 黑人操中国人逼视频| 亚洲欧美一区二区三区黑人| 在线亚洲精品国产二区图片欧美| 久久久久久久大尺度免费视频| 汤姆久久久久久久影院中文字幕| 老熟妇乱子伦视频在线观看 | 亚洲精华国产精华精| 老司机在亚洲福利影院| 欧美亚洲 丝袜 人妻 在线| 亚洲国产欧美日韩在线播放| 日本一区二区免费在线视频| 在线观看人妻少妇| 国产亚洲午夜精品一区二区久久| 他把我摸到了高潮在线观看 | 超碰97精品在线观看| 99精品久久久久人妻精品| 成年av动漫网址| 两人在一起打扑克的视频| 亚洲第一av免费看| 精品一区在线观看国产| 一本一本久久a久久精品综合妖精| av国产精品久久久久影院| 大香蕉久久网| 国产成+人综合+亚洲专区| 日韩欧美免费精品| 老司机在亚洲福利影院| 成人手机av| 午夜视频精品福利| 老司机午夜十八禁免费视频| 欧美在线黄色| 91大片在线观看| 日本黄色日本黄色录像| 精品国产国语对白av| 国产成人精品无人区| 妹子高潮喷水视频| 亚洲国产av影院在线观看| 国产精品久久久av美女十八| 高清在线国产一区| 国产精品九九99| 亚洲av成人一区二区三| 亚洲一码二码三码区别大吗| 亚洲欧美激情在线| 99香蕉大伊视频| 精品人妻一区二区三区麻豆| 日本猛色少妇xxxxx猛交久久| 亚洲成人国产一区在线观看| 黄色毛片三级朝国网站| 亚洲九九香蕉| 欧美激情久久久久久爽电影 | 青青草视频在线视频观看| 亚洲人成77777在线视频| 欧美国产精品va在线观看不卡| 一级,二级,三级黄色视频| 啦啦啦啦在线视频资源| 亚洲伊人久久精品综合| 国产av一区二区精品久久| 久久久久久久精品精品| 亚洲人成77777在线视频| 色老头精品视频在线观看| 欧美日韩成人在线一区二区| 99国产极品粉嫩在线观看| 一级片'在线观看视频| 欧美老熟妇乱子伦牲交| 中文字幕高清在线视频| 国产福利在线免费观看视频| 麻豆乱淫一区二区| 人人妻人人澡人人看| 久久精品熟女亚洲av麻豆精品| 欧美精品av麻豆av| 国产亚洲欧美在线一区二区| 老司机亚洲免费影院| 免费看十八禁软件| 成人av一区二区三区在线看 | 另类亚洲欧美激情| 中文字幕av电影在线播放| 亚洲精品中文字幕一二三四区 | 国产一区二区激情短视频 | 99香蕉大伊视频| 国产又爽黄色视频| 9色porny在线观看| 欧美激情高清一区二区三区| 99热全是精品| 日韩人妻精品一区2区三区| 美女扒开内裤让男人捅视频| 午夜两性在线视频| 一二三四社区在线视频社区8| 一区二区三区四区激情视频| 人人妻人人爽人人添夜夜欢视频| 久久人人爽av亚洲精品天堂| kizo精华| 国产成人免费观看mmmm| 亚洲精品国产av蜜桃| 欧美激情高清一区二区三区| 一本大道久久a久久精品| 欧美日韩亚洲综合一区二区三区_| 99国产精品99久久久久| 欧美日韩中文字幕国产精品一区二区三区 | 热99re8久久精品国产| 免费在线观看影片大全网站| 欧美日韩视频精品一区| 欧美精品亚洲一区二区| 欧美午夜高清在线| 久久久精品94久久精品| 国产日韩欧美在线精品| 欧美少妇被猛烈插入视频| www.av在线官网国产| 免费在线观看视频国产中文字幕亚洲 | 天天躁日日躁夜夜躁夜夜| 日韩大码丰满熟妇| 精品亚洲乱码少妇综合久久| 国产片内射在线| 天天添夜夜摸| 亚洲欧美精品综合一区二区三区| 精品一区二区三区四区五区乱码| 国产不卡av网站在线观看| 高清欧美精品videossex| 免费人妻精品一区二区三区视频| 亚洲成人免费av在线播放| 青青草视频在线视频观看| 国产成人精品久久二区二区91| 久久久久久人人人人人| 大片电影免费在线观看免费| 中文字幕人妻丝袜一区二区| 中文字幕精品免费在线观看视频| 制服诱惑二区| 国产又爽黄色视频| 亚洲中文日韩欧美视频| 色精品久久人妻99蜜桃| 成年av动漫网址| 国产高清videossex| 免费女性裸体啪啪无遮挡网站| 老汉色∧v一级毛片| 欧美日韩亚洲国产一区二区在线观看 | 不卡av一区二区三区| 欧美人与性动交α欧美精品济南到| avwww免费| 咕卡用的链子| 一区二区三区精品91| 亚洲欧美色中文字幕在线| 我的亚洲天堂| av超薄肉色丝袜交足视频| 亚洲少妇的诱惑av| 淫妇啪啪啪对白视频 | 一本大道久久a久久精品| 亚洲天堂av无毛| 午夜福利一区二区在线看| 人妻人人澡人人爽人人| 亚洲熟女精品中文字幕| 国产免费av片在线观看野外av| 人妻人人澡人人爽人人| 波多野结衣一区麻豆| 午夜福利乱码中文字幕| 丰满少妇做爰视频| 后天国语完整版免费观看| 一区福利在线观看| 91精品伊人久久大香线蕉|