亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        巖橋角度對鎖固段滑坡破壞影響的模型試驗研究

        2022-11-04 02:10:34唐垚
        人民長江 2022年10期
        關鍵詞:裂紋模型

        賀 春 蕾,朱 星,,劉 俊 峰,唐垚,霍 冬 冬,亓 星

        (1.成都理工大學 地質(zhì)災害防治與地質(zhì)環(huán)境保護國家重點試驗室,四川 成都 610059; 2.成都理工大學 信息科學與技術學院,四川 成都 610059; 3.四川輕化工大學 土木工程學院,四川 自貢 643000)

        0 引 言

        鎖固型巖質(zhì)滑坡是大型巖質(zhì)滑坡的主要類型。隨著鎖固段突然發(fā)生脆性斷裂產(chǎn)生震動,滑坡體獲得加速度,在與自身重力的耦合作用下,滑坡體突然失穩(wěn)啟動,具有極強的隱蔽性和突發(fā)性,經(jīng)常造成重大經(jīng)濟損失和人員傷亡。在巖質(zhì)滑坡鎖固段的研究方面,黃潤秋等[1]根據(jù)鎖固段的存在形式將鎖固段滑坡細分為5類:“三段式”、“擋墻式”、阻滑塊體型、線性多級型和階梯狀多級型?!叭问健蹦J绞菐r質(zhì)滑坡中典型的地質(zhì)-力學模式之一?!叭问健蹦J绞侵告i固段滑坡的變形破壞具有三段式的發(fā)育特征:即前緣緩傾結構面蠕滑、后緣結構面拉裂、中部鎖固段被剪斷[2]。在“三段式”巖質(zhì)滑坡中,鎖固段巖體在滑坡重力作用下產(chǎn)生應力集中,并對巖質(zhì)滑坡的變形破壞和整體失穩(wěn)起絕對控制作用。因此,巖質(zhì)滑坡中鎖固段的相關研究一直是大型巖質(zhì)滑坡防災減災研究中的熱點和難點問題。

        針對鎖固段開展物理模型試驗,研究其變形破壞機理和災變規(guī)律對該類巖質(zhì)滑坡成因機理認知、防災能力提升具有一定的科學指導意義。Chen等[3]采用聲發(fā)射技術對具有不同長度巖橋的直剪試驗和大型滑坡模型試驗進行監(jiān)測,發(fā)現(xiàn)聲發(fā)射技術在這兩類試驗中均能準確定位聲發(fā)射源點,反映了巖樣內(nèi)部裂隙的破裂過程,揭示了巖質(zhì)邊坡的破壞機理。黃潤秋等[4]通過物理試驗研究,發(fā)現(xiàn)在加載過程中應力和能量在鎖固段積累到一定程度時瞬間釋放,導致鎖固段突然發(fā)生脆性破壞,揭示了鎖固段控制作用下的大型巖質(zhì)滑坡的發(fā)生機理。陳竑然等[5]構建了鎖固段力學模型,結合物理試驗與數(shù)值模擬,揭示了天然鎖固段之間的力學作用模式。Guo等[6]對具有不同連續(xù)率、正應力、巖橋布置方式和節(jié)理角的鎖固段試件進行了直剪試驗,試驗結果表明,試件的強度受巖橋的控制,不同巖橋的布置方式對試件的抗剪強度影響較大,試件的最終破壞源于巖橋核心區(qū)的貫通發(fā)展。然而,針對具有不同巖橋角度的“三段式”滑坡鎖固段模型的破壞模式以及力學特性差異方面的研究還相對較少。

        鑒于此,本文開展了基于聲發(fā)射技術的不同巖橋角鎖固段模型破壞試驗研究,旨在研究巖橋角度對鎖固段強度及破壞模式的影響。近年來,聲發(fā)射技術作為巖石力學研究領域的一種有效方法常被用于巖石損傷破壞的研究[7-9]。機器學習在利用聲發(fā)射信號預測材料斷裂模式中具有良好的效果[10]。本次研究采用聲發(fā)射監(jiān)測技術對不同巖橋角鎖固段模型在單軸加載條件下的變形破壞過程進行同步觀測,一方面研究鎖固段受荷載作用下變形破壞過程的聲發(fā)射特征變化規(guī)律,研究不同巖橋角度對鎖固段破壞過程的影響;另一方面,大量研究證明,由于聲發(fā)射信號蘊含了豐富的巖石微破裂信息,利用高斯混合模型(GMM)對鎖固段不同破壞模式所對應聲發(fā)射特征進行聚類分析并建立分類判據(jù),實現(xiàn)通過聲發(fā)射特征分析來預測其對應的破壞模式,可為破壞過程的階段識別和監(jiān)測提供一種技術方法。

        1 試驗方案

        1.1 巖石試樣制備

        試驗選用完整砂巖制作了如圖1所示的4類不同巖橋角度的鎖固段模型,每種模型各制備3個。鎖固段模型前緣、后緣均預制裂縫,前緣裂縫中填充低強度材料石蠟用于模擬鎖固段前緣蠕滑破壞,鎖固段模型的基底和兩側均由涂有油脂的鋼板約束。限于篇幅,本文從12個砂巖制成的鎖固段模型中選出了4個進行分析,巖橋角分別為70°,90°,110°,130°。模型的鎖固段模型高度均為250 mm,頂寬133 mm,底寬250 mm,基座高度為20 mm,模型厚度為3 cm,具體物理模型尺寸參數(shù)見表1。鎖固段巖橋長度均為50 mm,不同巖橋的前緣緩傾裂隙長度和傾角保持一致,α和β分別表示張拉裂隙傾角和緩傾裂隙傾角。

        1.2 測試方法

        本次試驗所用加載裝置為WHY-1000微控壓力試驗機,系統(tǒng)可提供最大軸壓1 000 kN。試驗采用靜

        表1 物理模型尺寸參數(shù)Tab.1 Size and parameters of physical model

        態(tài)加載方式,加載過程分兩階段進行:① 預加載階段采用位移控制,加載速度為0.5 mm/min,目標荷載設定為1 kN;② 正式加載階段,加載速度為0.2 mm/min,由軸向應變控制。試驗通過聲發(fā)射系統(tǒng)對模型變形破壞全過程所釋放的聲發(fā)射信號進行同步觀測。4個聲發(fā)射傳感器均勻布置在鎖固段巖橋部位的左右兩側(見圖1)。試驗之前,應對單軸加載系統(tǒng)和聲發(fā)射采集系統(tǒng)進行檢測調(diào)試,確保試驗的順利開展。同時,為了保證單軸加載系統(tǒng)和聲發(fā)射監(jiān)測系統(tǒng)記錄的數(shù)據(jù)同步,兩個系統(tǒng)應同時開始采集數(shù)據(jù)。

        2 試驗結果分析

        2.1 聲發(fā)射特征參數(shù)統(tǒng)計分析

        2.1.1破壞過程聲發(fā)射能量特征

        聲發(fā)射能量是聲發(fā)射源釋放的彈性能量,聲發(fā)射能量的釋放能直觀反映巖石內(nèi)部的破裂變化。因此,本文選用聲發(fā)射能量來比較具有不同巖橋角度的鎖固段模型的內(nèi)部活動變化。通過對比發(fā)現(xiàn),不同巖橋角鎖固段模型中的聲發(fā)射能量隨著加載的進行整體均呈現(xiàn)一個增長的過程。聲發(fā)射能量、累積能量以及荷載隨著時間的變化發(fā)展趨勢如圖2所示?;趯r間-加載曲線和聲發(fā)射能量變化過程的綜合分析,本文將鎖固段模型的加載過程分為3個階段。

        階段Ⅰ,加載初始階段。在這個階段,荷載隨時間的增長不明顯,幾乎沒有太大的增幅,對4個鎖固段模型而言,荷載均在0~2 kN這個范圍內(nèi)浮動。出現(xiàn)這種現(xiàn)象是因為砂巖試樣內(nèi)部存在一些原始微裂紋和孔隙,在前期加載過程中原始微裂紋被壓實。根據(jù)加載-時間曲線,說明試驗中不同鎖固段模型內(nèi)部微裂縫和孔隙被壓實的荷載水平保持在2 kN左右。此階段的聲發(fā)射信號主要來源于試樣內(nèi)部的摩擦作用,并沒有出現(xiàn)真正的巖石破裂。因此,這個階段僅有很小一部分聲發(fā)射能量釋放。

        階段Ⅱ,裂縫穩(wěn)定發(fā)展階段。將鎖固段模型峰值荷載定義為P,此階段的荷載水平維持在2 kN~85%P這個范圍內(nèi)。在這個階段,荷載隨著時間呈現(xiàn)出線性增長的趨勢,模型內(nèi)部逐漸產(chǎn)生新的裂紋,模型內(nèi)部聲發(fā)射活動逐漸變得活躍起來,圖2中可以清楚看到聲發(fā)射能量釋放有了明顯的增長。

        階段Ⅲ,裂縫不穩(wěn)定發(fā)展及破壞階段。鎖固段模型此階段的荷載水平處于85%P~P之間。這個階段,試樣內(nèi)部裂紋進入不穩(wěn)定發(fā)展階段,隨著負載的增加,裂縫在鎖固段上發(fā)展并且貫通,鎖固段結構面開始發(fā)生破壞,此時巖石內(nèi)部聲發(fā)射活動異常劇烈,釋放出大量的聲發(fā)射能量,其能量數(shù)值遠遠大于前兩個階段,從圖3可以看到聲發(fā)射累積能量曲線出現(xiàn)一個幾乎垂直增長的趨勢。

        2.1.2破壞過程聲發(fā)射計數(shù)分布

        聲發(fā)射計數(shù)指的是在試驗中聲發(fā)射信號超過閾值的振鈴脈沖次數(shù),它可以很好地揭示巖石材料的破裂機理。圖3繪制了4種巖橋角度的砂巖試樣聲發(fā)射瞬時計數(shù)、累積計數(shù)和荷載隨時間的變化趨勢。可以發(fā)現(xiàn)它們之間有著一致的變化規(guī)律:在加載初始階段(階段Ⅰ),僅有很少的聲發(fā)射計數(shù);加載中期階段(階段Ⅱ),聲發(fā)射計數(shù)逐漸增加,累積聲發(fā)射計數(shù)隨著時間呈現(xiàn)出線性增長的趨勢;加載后期階段(階段Ⅲ),由于裂紋的不穩(wěn)定發(fā)展,試件逐漸進入破壞階段,此時可以觀察到聲發(fā)射累計計數(shù)曲線出現(xiàn)陡增趨勢。聲發(fā)射計數(shù)的變化過程和聲發(fā)射能量的變化過程基本同步。

        為了進一步對聲發(fā)射計數(shù)數(shù)據(jù)進行分析,將獲得的聲發(fā)射計數(shù)數(shù)據(jù)分階段作圖,如圖4所示。在階段Ⅰ中,70°,90°,110°,130° 4種模型的聲發(fā)射計數(shù)均較少,計數(shù)主要集中在40次以下,出現(xiàn)這個現(xiàn)象的原因是由于加載前期鎖固段模型內(nèi)部主要發(fā)生的是原始裂縫被壓實活動,此時并未產(chǎn)生新裂縫。在加載階段Ⅱ時,4個鎖固段模型的聲發(fā)射計數(shù)均呈現(xiàn)上升趨勢,計數(shù)次數(shù)主要集中在80以下,計數(shù)次數(shù)相較于階段Ⅰ提高了一倍,4種模型最高聲發(fā)射計數(shù)分別達到了140,180,200和190,鎖固段模型內(nèi)部聲發(fā)射活動劇烈,說明此時新裂紋逐漸產(chǎn)生與發(fā)展。隨著加載的進行,模型內(nèi)部產(chǎn)生的裂紋延伸并且逐漸貫通,模型進入到破壞階段(階段Ⅲ),隨著破裂結構面的產(chǎn)生,聲發(fā)射計數(shù)劇烈增加。在Ⅲ階段,模型70°、90°、110°、130°的最高聲發(fā)射計數(shù)均達到20 000次,模型90°的最高聲發(fā)射計數(shù)有18 000多次,聲發(fā)射活動頻率顯著增大,此時模型內(nèi)部裂縫劇烈發(fā)展,直至模型破壞失效。

        2.2 RA-AF分布特征及破壞模式

        本文采用了高斯混合模型(Gaussian Mixture Model,GMM)方法對4類鎖固段模型3個加載階段的RA值和AF值進行分析,進而判斷出鎖固段模型的破裂模式。GMM是一種常用的機器學習聚類算法,該方法使用高斯概率密度函數(shù)量化事物,將一個事物分解為若干個高斯概率密度函數(shù),理論上能擬合出任意類型的分布,常用于解決同一集合下有多重不同分布的問題[11-13]。GMM的詳細信息如下[14-15]。

        如果一個隨機變量y的概率密度函數(shù)如公式(1) 所示,則y服從正態(tài)分布或高斯分布:

        =N(y|μ,σ2);-∞0

        (1)

        式(1)等效于y~N(μ,σ2),表示隨機變量y服從正態(tài)分布,μ是平均值,σ2是方差。GMM是一種參數(shù)化的概率密度函數(shù),表示為高斯分量密度Q的加權和,如公式(2)所示:

        (2)

        式中:y是一個D維連續(xù)向量;ωi是混合權重;Ni(y|μi,∑i)是高斯密度的分量。

        每個單模態(tài)的高斯分量密度的形式是一個D變量高斯函數(shù),表示為

        (3)

        (4)

        (5)

        (6)

        (7)

        (8)

        (9)

        本文通過采集到的聲發(fā)射參數(shù)信息計算RA和AF值,為了更直觀地分析不同鎖固段模型的RA-AF分布差異,基于GMM繪制了關于RA值和AF值的分布密度圖,如圖5所示。其中,密度圖的顏色色階代表了不同的數(shù)據(jù)密度,色階從紅色到紫色表示數(shù)據(jù)密度從稀疏變?yōu)槊芗?,藍色區(qū)域和紫色區(qū)域的數(shù)據(jù)密度更高,而紅色區(qū)域的數(shù)據(jù)密度較為稀疏。其中,RA和AF的值由式(10)和式(11)來定義:

        (10)

        (11)

        GMM依據(jù)該分布特征將破壞模式分為拉伸破裂和剪切破裂,并賦予所有數(shù)據(jù)相對應的分類標簽。為保證模型的泛化能力,本文隨機選取數(shù)據(jù)集中20%的數(shù)據(jù)作為測試集。以往研究結果證明:剪切裂紋通常具有較大的RA值和較小的AF值,而拉伸裂紋則具有較大的AF值和較小的RA值[16-17]。以圖5(a)為例對RA-AF密度圖進行說明,圖中數(shù)據(jù)密度較高的區(qū)域被橢圓虛線標注出來,由豎放橢圓圈出來的數(shù)據(jù)區(qū)域主要代表剪切裂紋,反之則代表拉伸裂紋。根據(jù)GMM聚類結果,計算了鎖固段模型內(nèi)部剪切裂縫占比情況,如表2所列。70°,90°,110°,130°模型在破壞階段的剪切破裂占比分別為32%,37%,56%,22%,該剪切裂紋占比情況與RA-AF密度圖反映出的剪切和拉伸裂紋分布情況相一致。

        有研究表明,初級裂縫一般為拉伸裂縫并且在原有裂縫的尖端或尖端附近引發(fā)[18-19]。本次試驗通過工業(yè)相機記錄了模型表面宏觀裂紋的形成與發(fā)展過程,不同巖橋角度的鎖固段模型均先在前緣緩傾裂縫的后端部因前緣蠕滑受拉產(chǎn)生一條微裂縫,這條微裂縫與前緣緩傾裂縫呈90°并往模型頂部發(fā)展,這條裂縫被稱為機翼裂縫[20-21]。機翼裂縫發(fā)展形成以后,繼續(xù)加載,各個模型的加載形態(tài)開始有所差異:70°模型在機翼裂縫形成后不久,后緣預制裂縫下端部受拉,出現(xiàn)一條與前緣機翼裂縫平行的微裂縫,鎖固段模型被破壞;90°模型隨著加載的進行,又在后緣預制裂縫下端出現(xiàn)一條朝著模型底部發(fā)展的裂縫,能量不斷在鎖固段被積累,最終該裂縫發(fā)展到與前緣緩傾裂縫后端相連接,模型突然被破壞;110°模型和130°模型的裂縫發(fā)展情況比較相似,兩者產(chǎn)生的機翼裂縫不斷向模型上端發(fā)展,幾乎快要與模型上端部相貫通,當加載達到模型極限承載能力時,鎖固段發(fā)生斷裂,兩者均發(fā)生突發(fā)性破壞。

        表2 不同巖橋角度鎖固段模型在不同加載階段時的剪切裂紋占比Tab.2 Proportions of shear cracks in different slope models at different loading stages

        通過對4種鎖固段模型破壞模式的分析,可以發(fā)現(xiàn)鎖固段模型的巖橋部分是抵抗自重和外荷載、維護鎖固段整體穩(wěn)定性的關鍵部位,鎖固段的破壞往往具有突發(fā)性特征。鎖固段模型最終破壞形式如圖6所示。將鎖固段模型上的裂縫發(fā)展過程與計算出的模型剪切破裂占比相結合,可以總結出這4種模型的最終破壞形式:巖橋角為70°和90°時,發(fā)生是以拉伸破裂為主的拉-剪混合破壞模式;巖橋角為110°時,鎖固段模型是以剪切破裂為主的拉-剪混合破壞模式;巖橋角為130°時,鎖固段模型是以受壓產(chǎn)生的劈裂模式為主。

        2.3 力學特性分析

        根據(jù)以往的研究結果,模型的失效通常與相對應的裂紋啟動模式有關[22]。對比4種鎖固段模型的位移-荷載曲線,這4條曲線在單軸加載后期均可觀察到在某一時刻的荷載瞬時下降了一部分,繼而又往上增長,直到加載達到模型的極限承載能力時,鎖固段模型發(fā)生破壞。導致這種現(xiàn)象的原因是,鎖固段模型在加載過程中,首先是在前緣緩傾裂縫后端部產(chǎn)生一條宏觀裂縫,該裂縫的產(chǎn)生導致荷載曲線出現(xiàn)第一次下降,隨著單軸加載的繼續(xù)進行,荷載再一次呈現(xiàn)上升趨勢,直至鎖固段結構面上主要裂縫的貫通,鎖固段的強度完全丟失,試件宣告破壞,這是一個強烈的應力釋放過程。該破壞為脆性破壞,破壞前沒有明顯征兆,具有突發(fā)性強、危險性大的特點。如圖2所示,在4種鎖固段模型加載過程中,有3種模型(分別為70°,90°,130°模型)在最終破壞前均出現(xiàn)應力回彈現(xiàn)象,并且該現(xiàn)象均發(fā)生于峰值荷載的85%以后。出現(xiàn)此現(xiàn)象的原因可能是因為鎖固段某部分開始破壞,但仍具有一定的鎖固能力。因此,可將此現(xiàn)象作為鎖固段模型臨近破壞的前兆信息。

        圖7繪制了4種鎖固段模型的位移-荷載變化曲線圖,圖中70°,90°,110°,130°模型的峰值荷載分別為5.94,16.28,38.54,59.82 kN。對制備的不同巖橋角的鎖固段模型峰值荷載取平均值,從θ=70°~130°,平均峰值荷載分別為6.18,17.46,38.27,49.76 kN??梢钥闯觯i固段模型的平均峰值荷載隨著巖橋角的增大而不斷增大,巖橋角越小,鎖固段模型越早達到峰值點,平均峰值荷載越小。130°模型的平均峰值荷載是70°模型平均峰值載荷的8倍。影響其峰值荷載的原因有兩個:巖橋角度的大小以及后緣張拉裂隙傾角的大小變化。由表1可以看出,隨著巖橋角度的增大,后緣張拉裂隙傾角也隨之線性增大,鎖固段模型承壓巖石材料增加,最終導致鎖固段模型所能承受的峰值荷載增大。圖8非常直觀地反映出鎖固段模型的平均峰值荷載隨著鎖固段巖橋角度增大而不斷增加的趨勢,幾乎呈線性增長趨勢。由此可以得出結論,在前緣緩傾裂隙長度和傾角一定時,鎖固段巖橋角度的大小對于鎖固段模型的強度有很大影響,巖橋角度越大,后緣張拉裂隙傾角越大,鎖固段模型承壓材料越多,故鎖固段模型的抗壓能力越好。

        因此,根據(jù)本次試驗結果,為鎖固段滑坡的防治提供了借鑒:當“三段式”滑坡的巖橋角處于70°~130°范圍時,巖橋角度越大,鎖固段能承受的荷載越大,越有利于鎖固段穩(wěn)定性;當巖橋角處于較小水平時,鎖固段的危險系數(shù)越大,應針對實際情況采取相應安全防護措施。

        3 結 論

        本文采用聲發(fā)射技術對不同巖橋角度的鎖固段模型進行監(jiān)測,分析了鎖固段模型的聲發(fā)射特性及其變化規(guī)律,了解了不同巖橋角度的鎖固段模型的破裂形式,得出了以下幾點結論。

        (1) 基于聲發(fā)射能量和聲發(fā)射計數(shù)變化趨勢,發(fā)現(xiàn)鎖固段模型的破裂過程有明顯的分段特征,聲發(fā)射能量和計數(shù)變化趨勢基本一致。在加載前期階段(階段Ⅰ),聲發(fā)射能量與計數(shù)并未有太大的變化;加載中期階段(階段Ⅱ),聲發(fā)射總能量與總計數(shù)均呈現(xiàn)線性增長的趨勢;加載后期階段(階段Ⅲ),聲發(fā)射總能量與總計數(shù)均發(fā)生陡增。聲發(fā)射參數(shù)的變化形象地描繪了巖石內(nèi)部裂紋的動態(tài)演化過程。

        (2) 對鎖固段模型整個加載過程進行分析,發(fā)現(xiàn)4種鎖固段模型從開始加載到最終破壞均會產(chǎn)生兩條影響巖石強度的主要裂縫:初級裂縫是由于鎖固段模型前緣蠕滑導致前緣緩傾裂縫尖端或尖端附近引發(fā)的拉伸裂縫,該裂縫與前緣緩傾裂縫呈90°傾角并往模型頂部發(fā)展;第二條裂縫是鎖固段結構面上產(chǎn)生的剪切裂縫,這條裂縫的貫通發(fā)展是導致鎖固段模型破壞的關鍵所在。

        (3) 結合機器學習算法,對聲發(fā)射參數(shù)進行非監(jiān)督聚類分析,總結出4種鎖固段模型的最終破壞形式:巖橋角小于90°時,發(fā)生的是以拉伸破裂為主的拉-剪混合破壞模式;巖橋角為110°時,鎖固段模型發(fā)生的是以剪切破裂為主的拉-剪混合破壞模式;巖橋角為130°時,鎖固段模型主要產(chǎn)生因受壓導致的劈裂破壞模式。

        (4) 鎖固段巖橋角度的大小對鎖固段模型的強度有很大影響,不同巖橋角度的鎖固段模型能夠承受不同的峰值荷載。鎖固段模型的承載能力隨巖橋角度的增大而增大,巖橋角度越小,鎖固段模型越早達到峰值點,平均峰值荷載越小;巖橋角度越大,鎖固段模型的抗壓能力越強。

        猜你喜歡
        裂紋模型
        一半模型
        裂紋長度對焊接接頭裂紋擴展驅動力的影響
        一種基于微帶天線的金屬表面裂紋的檢測
        重要模型『一線三等角』
        重尾非線性自回歸模型自加權M-估計的漸近分布
        Epidermal growth factor receptor rs17337023 polymorphism in hypertensive gestational diabetic women: A pilot study
        微裂紋區(qū)對主裂紋擴展的影響
        3D打印中的模型分割與打包
        FLUKA幾何模型到CAD幾何模型轉換方法初步研究
        預裂紋混凝土拉壓疲勞荷載下裂紋擴展速率
        97女厕偷拍一区二区三区| 亚洲午夜精品久久久久久人妖| 亚洲深夜福利| 在线观看免费的黄片小视频| 丝袜美腿在线观看一区| 中文字幕人妻少妇引诱隔壁| 久久综合第一页无码| 国产日产久久福利精品一区| 国产自拍一区在线视频| 亚洲精品欧美精品日韩精品| 人与嘼交av免费| 国产精品久久久久久久y| 国产精品一二三区亚洲| 99久久婷婷国产综合亚洲| 亚洲中文字幕无码久久2020| 男人深夜影院无码观看| 国产在线一区二区三区四区乱码| 亚洲无线一二三四区手机| 99re热这里只有精品最新| 国产精品女丝袜白丝袜| 亚洲中文字幕精品久久a| 国产又粗又猛又黄又爽无遮挡| 精品性高朝久久久久久久| A亚洲VA欧美VA国产综合| 亚洲一区二区三区重口另类| 中文字幕日本人妻久久久免费| 亚洲精品国产成人AV| 国产一级一片内射在线| 亚洲tv精品一区二区三区| 性一交一乱一乱一视频| 在线视频中文字幕乱人伦| 精品国产亚洲av高清日韩专区| 亚洲av综合av一区| chinese国产乱在线观看| 女同国产日韩精品在线| 国产一品二品三品精品在线| 亚洲av无码不卡久久| 亚洲AV无码国产精品久久l| 国产av无毛无遮挡网站| 精品无码国产自产拍在线观看蜜| 精品午夜久久网成年网|