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

    基于固有頻率的錯斷式墜落危巖體穩(wěn)定系數計算模型

    2024-04-11 03:44:28張曉勇謝謨文高世崇李雙全
    工程力學 2024年4期
    關鍵詞:母巖后緣固有頻率

    張曉勇,謝謨文,張 磊,杜 巖,高世崇,李雙全

    (1.北京科技大學土木與資源工程學院,北京 100083;2.北京中關村智連安全科學研究院有限公司,北京 100083)

    危巖體失穩(wěn)崩塌是危巖體自母巖脫離并突然崩落的地質現象,是多山地區(qū)常見的地質災害之一。危巖體失穩(wěn)崩塌往往能夠帶來巨大的生命財產安全損失,而我國平均每年發(fā)生危巖體崩塌災害次數約為2000 余起[1],嚴重制約我國經濟建設發(fā)展。

    針對危巖體崩塌失穩(wěn)的防控,眾多學者從危巖體崩塌機理、防控手段著手進行 研究。BERTRAN[2]基于地質地貌學對阿爾卑斯山脈Claix 崩塌形成機制進行了分析;BRAATHEN 等[3]通過對Norway境內巖質邊坡失穩(wěn)原因進行分析,認為軟質巖層對巖質邊坡穩(wěn)定影響較大;ISHIKAWA 等[4]通過對裂縫寬度、溫度進行監(jiān)測,發(fā)現裂縫擴展受低溫度影響;AREF[5]通過對也門山區(qū)危巖體綜合調查,并對該區(qū)域危巖失穩(wěn)崩塌致災因素進行總結,提出響應的防控措施;陳洪凱等[6]通過斷裂力學理論分析危巖體失穩(wěn)機理,并提出基于強度因子的穩(wěn)定系數計算方法,陳維等[7]通過對凹腔型危巖體研究發(fā)現,危巖體失穩(wěn)受巖腔深度、危巖體厚度和高度、主控裂隙深度等影響。王根龍等[8]認為危巖體失穩(wěn)是由于危巖體后緣裂隙尖端應力集中引起的。危巖體崩塌產生原因多類多樣,但危巖體穩(wěn)定程度直接受其后緣裂隙強度控制[9]。在分析各類危巖體失穩(wěn)機理基礎上,針對各類危巖體監(jiān)控技術應運,常見的監(jiān)測指標包括:位移、應力、應變及環(huán)境因素等[10]。但危巖體由穩(wěn)定到失穩(wěn)時常年累月地損傷積累過程,初期危巖體的形變緩慢而微小,直至臨近失穩(wěn)才會出現較大變形,另外,邊坡危巖體的變形及破壞是一個多因素相互作用的非線性過程[11-12],因此危巖體失穩(wěn)表現出空間位置隨機性、時間不確定性及突發(fā)性的特點,常見的基于位移、應力-應變等監(jiān)測指標在監(jiān)測預警中有一定的時間滯后性及不確定性[13],基于當前的監(jiān)測技術仍無法有效地防控危巖體失穩(wěn)。

    如何引入有效的監(jiān)測指標來定量分析危巖體穩(wěn)定系數,是實現危巖體安全評價的關鍵所在。借鑒結構力學分析方法,巖體本身可作為結構進行分析,危巖體后緣裂隙擴展是巖體結構損傷的表現。當前結構損傷識別主要集中在高危建筑、橋隧及機械等領域[14-17],而基于動力學特征參數的損傷識別及定位是常用方法之一,基于結構動力學的損傷識別方法所采用的主要模態(tài)參數包括:固有頻率、阻尼及振型等。同樣,巖體結構損傷后在時域、頻域上出現明顯的改變[18],眾多學者[19-22]通過實踐及試驗發(fā)現,隨著危巖體損傷程度變化,其振幅、頻率等動力特征有著較為明顯的變化。HAN 等[23]通過設計振動臺對具有不連續(xù)面高陡巖質邊坡的動力響應進行研究;葉陽升等[24]通過室內試驗發(fā)現墜落危巖體的固有頻率隨危巖后緣貫通率的增大而明顯降低;杜巖等[25]通過試驗及理論分析得出結論:隨著危巖體與基巖粘結程度降低,滑移型危巖體的固有頻率降低。VALENTIN 等[26]基于常時微動測量技術對現場危巖體進行測量,并建議采用頻譜特征分析監(jiān)測巖石損傷斷裂情況。DU 等[27]通過模擬凍融的方法模擬了危巖體崩落過程,并提出將基于固有頻率作為危巖體穩(wěn)定程度的表征指標的方法。當前危巖體動力學基本停留在理論及試驗角度以定性的角度得到危巖體損傷程度與固有頻率之間的關系,其主要分析方法是將危巖體作為質點進行考慮,而實際上危巖體與母巖都是可變形體,對危巖振動頻率與其后緣裂隙深度之間精確定量關系鮮有報道[28]。本文基于修正Timoshenko 梁理論,建立錯斷式墜落危巖體穩(wěn)定系數與固有頻率的定量關系,以期為危巖體崩塌防控提供有益借鑒。

    1 錯斷式墜落危巖體穩(wěn)定系數靜力學計算方法

    錯斷式墜落危巖體[29]又稱“懸掛式危巖體”或“壁掛式危巖體”,該類危巖體往往厚度較薄,呈片狀或長條狀,巖體結構豎向節(jié)理發(fā)育。危巖體底面臨空無支撐,危巖體與基巖的粘結面受由自重引起的張拉-剪切復合作用力,當作用力大于主控斷裂面強度時,巖體沿主控斷裂面拉剪崩落。如圖1 所示,錯斷式墜落危巖體形成主要原因為差異風化作用。由于坡體底部受河谷下切或人工開挖擾動,巖體在卸荷回彈作用下在后緣形成陡傾張拉裂縫,在自重及漸進風化作用下裂縫不斷加深,巖橋部位應力集中程度加劇進一步加深,直至巖橋突然破裂,危巖發(fā)生失穩(wěn)崩落。

    圖1 危巖體風化過程示意圖Fig.1 Schematic diagram of weathering process of dangerous rock mass

    錯斷式墜落危巖體受力示意圖如圖2 所示,危巖體模型高度為H,寬度為L, 厚度為B,后緣裂縫深度為h。基于該模型可計算得裂縫處截面所受最大拉應力 σmax和剪應力τ,二者表達式如下:

    圖2 危巖體受力示意圖Fig.2 Force diagram of dangerous rock mass

    式中, ρ /( kg·m-3)為巖體密度。

    如式(1)所示,隨著裂縫深度h增加,裂縫處截面上最大拉應力 σmax和剪應力τ隨之增大。若墜落危巖體以剪切破壞為主,則當前危巖體穩(wěn)定系數為Fs1:

    實際上,裂縫處截面受拉剪復合作用,拉剪復合作用條件下巖石抗拉強度遠小于單軸抗拉強度[30],拉剪復合作用條件下巖石的抗拉強度為:

    式中: σt為拉剪復合作用下巖石抗拉強度;Rt為巖石單軸抗拉強度;c為巖石內聚力;τ為巖石當前所受剪應力。由式(1)與式(3)可計算當前巖石抗拉強度為:

    若危巖體破壞形式為拉剪復合作用下張拉破壞,則當前危巖體穩(wěn)定系數為:

    實際上,危巖體破壞形式受幾何形式、巖體力學參數及后緣裂縫深度等因素,采用單一形式的破壞判據有失客觀,為安全起見,危巖體當前穩(wěn)定系數取上述兩種破壞形式所采用的穩(wěn)定系數最小值:

    式(6)所涉及的巖石材料力學參數及危巖體幾何形式往往是可測的,根據該式可基于巖體材料參數得到危巖體后緣裂隙臨界深度值。然而,隨著危巖體整體力學性質劣化,其后緣裂縫深度h不斷增大,而裂縫深度的動態(tài)變化無法經靜力學進行實時計算,危巖體的穩(wěn)定系數也就無法確定。因此,如何反映穩(wěn)定系數計算必需參數h的動態(tài)變化仍是亟需解決的問題。從動力學角度而言,任何物體的動力特征參數可反映該物體的本質特征,如材料參數、形狀、邊界條件等。固有頻率是動力特征參數之一,能夠反映物體固有的特征參數。隨著危巖體后緣裂隙深度增加,母巖對危巖體的約束程度降低,固有頻率亦隨之變化,為分析固有頻率與危巖體后緣裂隙深度的關系,筆者對危巖體固有頻率計算方法進行分析。

    2 錯斷式墜落危巖體的固有頻率計算方法研究

    2.1 錯斷式墜落危巖體動力學方程建立

    錯斷式墜落危巖體動力學模型如圖2 所示,將危巖體與基巖分離的部分簡化為巖梁,巖梁長度為裂縫深度h, 寬度為L,厚度為B。

    由圖3 可知,隨著裂縫深度增大,巖梁長度h不斷增大,母巖對危巖體的約束長度H-h不斷減小,巖梁的幾何形式及邊界特征變化直接導致其振動特征發(fā)生改變,為獲取巖梁長度h對振動特征的影響,建立基于修正Timoshenko 梁理論的動力學方程[31],修正Timoshenko 梁理論的同時考慮了梁的剪切變形、轉動慣量,且不存在雙頻譜的問題,修正Timoshenko 梁理論方程式如下:

    圖3 危巖體動力學模型Fig.3 Dynamic model of dangerous rock mass

    式中:S為巖梁橫截面積;G為巖石剪切模量;EI為巖梁抗彎剛度;μ為巖梁剪切系數,若截面為矩形,則

    基于分離變量法可獲取方程的通解,獲取振型函數Y(x) 、彎矩引起的轉角函數 Φ(x)及剪切引起的轉角函數B(x):

    式中:

    通解見式(8)~式(10),式中的c1~c4為待求系數,需根據邊界條件進行求解。

    2.2 巖梁約束條件分析

    本文不考慮巖梁在X方向的平移運動,僅考慮梁的橫向彎曲。危巖體與母巖未分離部分受側向約束,橫向彎曲轉角近乎為0°,可忽略。將危巖體未分離部分的運動變形進行分解,如圖4 所示。未分離部分的運動變形包括剛性運動及變形,其中剛性運動又包括剛性平移與剛性旋轉?;谠摬糠值倪\動變形特征可知,該部分受側向母巖的抗搖擺約束及抗平移約束,另外剪切變形也提供抗上部巖梁的搖擺運動作用。

    圖4 運動分解示意圖Fig.4 Motion decomposition diagram

    由于母巖尺度遠大于危巖,因此,將未分離部分與母巖之間關系視為半無限空間地基與基礎的關系,可采用動力基礎半空間理[32]進行分析,如圖5 所示。

    圖5 半無限空間地基與基礎示意圖Fig.5 Semi-infinite space foundation bed and foundation schematic diagram

    半空間提供的抗平移剛度Ky與抗搖擺剛度Kψ計算方法如下:

    式中:G′為半無限空間材料剪切模量;r1為橫向運

    其中系數F1(α1) 計算方法[33]見表1,F2(α2)表達式見式(12):

    表1F1(α1)計算式Table 1 Formulas ofF1(α1)

    1) 抗巖梁搖擺剛度分析

    “基礎”部分可發(fā)生整體縱向剪切變性和局部剪切變形,兩種變形同時存在。首先考慮“基礎”整體的剪切變形,如圖4 所示。

    由材料力學可得“基礎”剪切變形提供的抗搖擺靜剛度為:

    “基礎”局部變形如圖5 所示,由材料力學理論可得“基礎”部分局部變形產生的抗搖擺靜剛度為:

    式中,Kψ1與Kψ2為“彈簧串聯”的關系?!凹羟小弊冃伟殡S著旋轉運動,其參振質量不可忽略,設“基礎”搖擺角度的簡諧運動為 Θ:

    則考慮參振部分轉動慣量的抗“基礎”搖擺動剛度Kψ3為:

    半無限空間的抗搖擺動剛度可由式(11)進行計算,但該方法是基于等效圓面積比擬所得,為更加精確求解“地基”對矩形“基礎”的動約束剛度,本文采用等效圓面積條件下動、靜剛度之比乘以矩形“基礎”靜約束剛度的方法得到矩形“基礎”底部動剛度。等效圓面積條件下“地基”對“基礎”的靜約束剛度為:

    結合式(11)與式(17),可計算得到動、靜剛度比 γ1:

    “地基”對矩形“基礎”的靜搖擺約束剛度為:

    式中:γ1′=-0.0001χ4+0.002χ3-0.0148χ2+0.139χ+0.363,

    結合式(18)與式(19)可得“地基”對矩形“基礎”動約束剛度:

    考慮“基礎”的轉動慣量造成“地基―基礎”對巖梁約束剛度的影響,“地基―基礎”對巖梁搖擺的約束剛度為:

    由于Kψ3與Kψ7是“彈簧串聯”的關系,基于“彈簧串聯”理論最終得到“地基―基礎”對巖梁搖擺的約束剛度:

    2) 抗巖梁橫向位移剛度分析

    根據式(11)可得,不考慮“基礎”參振質量條件下“地基”對與圓形等效面積的“矩形基礎”平移約束剛度。“地基”對圓形等效面積矩形“基礎”的平移靜約束剛度表達式如下:

    結合式(11)與式(23),可計算得到動、靜剛度比 γ2:

    “地基”對矩形“基礎”的平移靜約束剛度為:

    式中: γ2′=-0.0036χ2+0.1231χ+1.9932。

    結合式(24)與式(25)可得“地基”對矩形“基礎”的平移動約束剛度:

    由于“基礎”是質量體,其慣性條件需考慮,類比“地基―基礎”對巖梁搖擺約束剛度的計算方法,可得“地基―基礎”抗巖梁橫向位移剛度為:

    通過上述分析,可根據式(22)與式(27)獲得抗巖梁搖擺約束剛度與抗巖梁橫向位移約束剛度,為巖梁振型的求解提供明確的約束邊界條件。

    2.3 巖梁邊界條件分析及固有頻率求解

    危巖體懸空端受剪力及彎矩皆為0,可將懸空自由端x=h處邊界條件寫作:

    由于巖梁底部未分離部分受基巖限制,其橫向彎曲主要由彎矩引起,因此僅考慮梁未分離部分由彎矩引起的橫向彎曲及橫向平移。值得注意的是,修正Timoshenko 梁理論是基于梁內部相對變形及作用內力推導的,并不考慮整體剛性的位移或旋轉。而“半空間地基”需作為可變形體考慮,巖梁搖擺勢必引起剛性旋轉及剛性平移。基于該性質,可推導x=0時梁的橫向彎曲運動及橫向平移運動的邊界條件,即:

    式(30)中c1~c4為待求解系數,由線性代數理論可知,當c1~c4全部為0 時,式(30)無意義,因此,c1~c4不全為0,故c1~c4的系數組成的4 階矩陣D行列式為零,對式 |D|=0進行求解便可得到巖梁的固有頻率,其中矩陣D見式(31)。由于計算較為復雜,需要借助數值計算。式(31)中未知數僅為含ω項,利用編程軟件采用數值迭代法便可快速求解得到ω值,計算結果形式如圖6 所示,其中|D|第一個零點對應頻率值即為第一階固有頻率。

    圖6 固有頻率計算示意圖Fig.6 Diagram of natural frequency calculation

    3 基于固有頻率的錯斷式墜落危巖體穩(wěn)定系數計算模型

    3.1 固有頻率指標敏感性分析

    危巖體穩(wěn)定系數計算服務于危巖體失穩(wěn)的防控,危巖體失穩(wěn)防控需要明確兩個問題:① 危巖體穩(wěn)定程度是否發(fā)生劣化;② 當前危巖體穩(wěn)定系數是否達到臨界值。若無法實時地確定危巖體穩(wěn)定程度是否發(fā)生劣化,則當前的穩(wěn)定系數的計算不具有任何意義,而反映危巖體穩(wěn)定程度或后緣裂縫深度指標的敏感性是最為關鍵的,因此表征危巖體后緣裂縫深度的關鍵參數不僅要能夠反映當前危巖體穩(wěn)定程度,還要對危巖體后緣裂縫深度變化有足夠的敏感性。

    為了危巖體固有頻率對危巖體后緣裂縫深度變化的敏感性,本文采用數值軟件ABAQUS 數值計算方法進行分析。由于母巖相對危巖尺寸無限大,且母巖的位移相對危巖近乎為0,因此,本次模擬中母巖尺寸為60 cm×60 cm×60 cm,并對母巖模型進行固支。危巖模型尺寸:厚度B為10 cm、寬度L為10 cm、高度H為40 cm。數值模型如圖7 所示。

    考慮危巖體與母巖一般材料相近或一致,本次模擬對母巖和危巖賦值同種材料屬性。數值模擬旨在分析危巖體固有頻率對裂縫深度的敏感性,暫不考慮危巖體的強度力學參數。材料各參數見表2。模擬過程中采用不斷加深后緣裂縫深度的方式模擬裂縫深度變化,每次加深長度為3 cm,裂縫深度h最大為36 cm。本次數值計算重點考察在平面上的彎曲運動,由于模型是三維的,危巖體對應的模態(tài)不一定全部在單一平面做彎曲運動,裂縫深度為3 cm時,考察的平面彎曲運動為第4 階模態(tài);裂縫深度為6 cm~27 cm 時,考察的平面彎曲運動為第2 階模態(tài);裂縫深度為27 cm~36 cm 時,考察的平面彎曲運動為第1 階模態(tài)。危巖體固有頻率理論值與模擬值對比如圖8 所示。

    表2 材料屬性表Table 2 Material property table

    圖8 固有頻率計算Fig.8 Natural frequency calculation

    如圖8 所示,隨著裂縫深度的增加,危巖體在考察平面彎曲運動固有頻率逐步降低,其降低的趨勢為:裂縫深度為3 cm~6 cm 時,固有頻率值降低速度相對較緩;裂隙深度為9 cm~24 cm時,固有頻率值降低速度較快;裂縫深度為27 cm~33 cm 時,固有頻率值降低速度又呈相對較緩的趨勢。裂隙深度由3 cm 逐步變?yōu)?3 cm 時,固有頻率下降幅值約250 Hz,由此可見,固有頻率值對后緣裂縫深度的變化較為敏感,可以有效地表征危巖體穩(wěn)定程度是否發(fā)生劣化。數值模擬與理論值存在些許誤差,原因在于:① 數值模擬過程中計算單元中的偽應變能是不可避免的,偽應變能的存在使得模型變形變大,一定程度上使得固有頻率降低;② 由于理論模型的復雜性,本文采用動力基礎半空間理論時未考慮輻射阻尼,一定程度上使得固有頻率理論計算值比實際要高。

    3.2 基于固有頻率錯斷式墜落危巖體穩(wěn)定系數計算模型建立

    由于振動系統(tǒng)機械能主要集中在一階頻率,本文重點考察一階頻率,由數值模擬驗證了危巖體固有頻率與后緣裂縫深度呈負相關,危巖體固有頻率隨后緣裂縫深度增大而降低。因此,可根據當前危巖體固有頻率反演危巖體后緣裂縫深度。將固有頻率-裂縫深度曲線擬合得到關系式:

    將式(32)代入式(2)、式(5)及式(6)中可得錯斷式墜落危巖體穩(wěn)定系數計算式:

    4 試驗驗證

    為了說明本算法的正確性及適用性,筆者在室內澆筑危巖-母巖模型,考慮到實際條件中危巖體與母巖材料屬性一致,因此試驗中二者采用同種材料進行澆筑,澆筑材料及配比見表3,材料的力學參數見表4,母巖模型尺寸為長度30 cm、寬度30 cm、高度40 cm 的矩形。危巖體模型寬度10 cm、厚度10 cm、高度30 cm,為了模擬更加準確,母巖的頂部進行加載固定,模型示意圖如圖9 所示。

    表3 模型材料配比Table 3 Model material ratio

    表4 模型材料力學參數Table 4 Mechanical parameters of model materials

    圖9 模型示意圖Fig.9 Diagram of model

    如圖10 所示,實驗過程中,沿著危巖體與母巖體的接觸面向下進行切割裂隙。每次切割深度為2 cm,切割完成后靜置5 min,然后每隔3 s 用力錘輕輕在危巖體表面進行敲擊,并利用拾振器采集振動加速度數據,單次采集完成后重復切割步驟,當切割深度為26 cm 時,繼續(xù)切割過程中模型發(fā)生破壞。

    圖10 裂隙切割示意圖Fig.10 Fracture cutting schematic

    裂隙深度增加過程中,固有頻率隨裂隙深度變化數據如表5 所示。固有頻率ω-裂隙深度h曲線如圖11 所示。

    表5 固有頻率與穩(wěn)定系數值Table 5 Natural frequency and stability coefficient

    圖11ω-h關聯曲線Fig.11 The correlation curve ofω-h

    通過對固有頻率ω(理論值)-裂隙深度h進行擬合,可得式(34):

    穩(wěn)定系數Fs與固有頻率ω的關聯曲線如圖12所示,對Fs-ω曲線進行擬合,二者之間近似呈負線性相關的關系。由此可知,錯斷式墜落危巖體的穩(wěn)定系數完全可由固有頻率ω進行表示,可基于危巖體的固有頻率實時計算當前危巖體的穩(wěn)定系數。該研究成果為危巖體失穩(wěn)防控自動化監(jiān)測提供有益借鑒。

    圖12Fs-ω關聯曲線Fig.12 The correlation curve ofFs-ω

    5 結論

    針對錯斷式墜落危巖體特點,將危巖體―母巖概化為基礎―半無限空間地基模型,結合修正Timoshenko 梁理論得到危巖體固有頻率解析,獲取危巖體后緣裂隙深度、固有頻率及穩(wěn)定系數之間的一一對應關系,并通過數值分析及室內試驗驗證了固有頻率指標的敏感性及有效性。結論如下:

    (1) 錯斷式墜落危巖體安全程度受后緣裂縫深度影響,裂縫深度是評判危巖體安全系數的關鍵參數之一,裂縫深度的動態(tài)變化制約了靜力學計算方法的使用,而通過危巖體一階固有頻率的變化可較好地反映裂縫深度的變化。

    (2) 根據動力基礎半空間理論分析危巖體所受的動約束剛度,結合巖梁模型的運動邊界條件,根據修正Timoshenko 梁理論所得的振型通解可計算不同后緣裂縫深度條件下危巖體固有頻率值。

    (3) 通過數值試驗及室內試驗得出結論:錯斷危墜落巖體的一階固有頻率ω隨后緣裂縫深度h變化而變化的趨勢較為明顯,且試驗結果與計算結果吻合度較高,說明危巖體一階固有頻率可作為危巖體安全系數的有效關鍵參數指標。

    (4) 通過分析固有頻率ω與危巖體后緣裂隙深度h的關系,建立了基于ω反 算h的計算方法,形成了基于固有頻率的錯斷式墜落危巖體穩(wěn)定系數計算模型。

    本文僅分析矩形梁在平面內的彎曲振動,未分析其阻尼特征,且未分析更復雜幾何形式的危巖體動力學特征。

    猜你喜歡
    母巖后緣固有頻率
    水耕條件下兩類富鈣母巖發(fā)育土壤的系統(tǒng)分類歸屬及成因探討*
    土壤學報(2022年3期)2022-08-27 08:41:30
    現場測定大型水輪發(fā)電機組軸系的固有頻率
    大電機技術(2021年2期)2021-07-21 07:28:38
    紫色母巖作基質或覆蓋材料對景觀水體氮磷去除效果研究*
    土壤學報(2019年2期)2019-04-25 06:39:58
    機翼后緣連續(xù)變彎度對客機氣動特性影響
    柔性后緣可變形機翼氣動特性分析
    貴州母巖(母質)對土壤類型及分布的影響
    TNF-α和PGP9.5在椎體后緣離斷癥軟骨終板的表達及意義
    總溫總壓測頭模態(tài)振型變化規(guī)律研究
    現代機械(2015年1期)2015-01-15 03:07:54
    沂蒙山區(qū)不同母巖發(fā)育土壤物理性質
    A novel functional electrical stimulation-control system for restoring motor function of post-stroke hemiplegic patients
    久久韩国三级中文字幕| 伊人久久精品亚洲午夜| 日产精品乱码卡一卡2卡三| 日韩欧美三级三区| 亚洲婷婷狠狠爱综合网| 欧洲精品卡2卡3卡4卡5卡区| 国产免费男女视频| 最好的美女福利视频网| 欧美性猛交╳xxx乱大交人| 久99久视频精品免费| 国产高清三级在线| 夜夜看夜夜爽夜夜摸| 久久精品夜夜夜夜夜久久蜜豆| 免费看美女性在线毛片视频| 91精品国产九色| 在线播放无遮挡| 99久久久亚洲精品蜜臀av| 久久中文看片网| 欧美激情在线99| 欧美3d第一页| 久久久久国产网址| 久久亚洲精品不卡| 国模一区二区三区四区视频| 国产精品久久电影中文字幕| 九草在线视频观看| 色噜噜av男人的天堂激情| 亚洲av一区综合| 亚洲成a人片在线一区二区| 国产伦精品一区二区三区四那| 天堂中文最新版在线下载 | 成人av在线播放网站| 亚洲成人久久爱视频| 久久精品久久久久久久性| 99热6这里只有精品| 免费看日本二区| av黄色大香蕉| 亚洲欧美日韩高清在线视频| 欧美变态另类bdsm刘玥| 国产成人一区二区在线| 99视频精品全部免费 在线| 国产白丝娇喘喷水9色精品| 麻豆成人av视频| 精品少妇黑人巨大在线播放 | 永久网站在线| 蜜臀久久99精品久久宅男| av在线蜜桃| 国产亚洲5aaaaa淫片| 只有这里有精品99| 九九久久精品国产亚洲av麻豆| 床上黄色一级片| 又粗又硬又长又爽又黄的视频 | 性色avwww在线观看| 久久久国产成人免费| 久99久视频精品免费| 少妇人妻精品综合一区二区 | 国产乱人偷精品视频| 国产精品麻豆人妻色哟哟久久 | 99久久精品国产国产毛片| 国产麻豆成人av免费视频| 99热这里只有精品一区| 亚洲第一电影网av| 精品久久久久久久人妻蜜臀av| 中文字幕av在线有码专区| 欧美不卡视频在线免费观看| 精品一区二区三区人妻视频| 国内少妇人妻偷人精品xxx网站| 一区二区三区高清视频在线| 精品人妻视频免费看| 国产男人的电影天堂91| 日本免费a在线| 亚洲va在线va天堂va国产| 99热这里只有是精品在线观看| 内地一区二区视频在线| 成人一区二区视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 夜夜看夜夜爽夜夜摸| 国产三级在线视频| 美女国产视频在线观看| 国产一区二区三区在线臀色熟女| 午夜福利在线观看吧| 1024手机看黄色片| 欧美xxxx性猛交bbbb| 国产精品蜜桃在线观看 | 国产黄片美女视频| 精品一区二区三区人妻视频| 一区福利在线观看| 哪里可以看免费的av片| 欧美激情在线99| 可以在线观看的亚洲视频| 中文字幕制服av| 可以在线观看毛片的网站| 波多野结衣高清作品| 亚洲av免费在线观看| 国产视频首页在线观看| 亚洲精品国产av成人精品| 久久精品国产亚洲网站| 又爽又黄a免费视频| 日韩成人伦理影院| 午夜福利在线观看免费完整高清在 | 亚洲精品自拍成人| avwww免费| 欧美性猛交╳xxx乱大交人| а√天堂www在线а√下载| 99久国产av精品国产电影| 免费看日本二区| 午夜视频国产福利| h日本视频在线播放| 又粗又爽又猛毛片免费看| 久久精品国产亚洲av香蕉五月| 国产一区二区亚洲精品在线观看| 免费搜索国产男女视频| 三级毛片av免费| 中文字幕人妻熟人妻熟丝袜美| 欧美高清成人免费视频www| 国产熟女欧美一区二区| 亚洲国产精品sss在线观看| 国产av不卡久久| 久久6这里有精品| 国产成人aa在线观看| 日本在线视频免费播放| 麻豆成人av视频| ponron亚洲| 菩萨蛮人人尽说江南好唐韦庄 | 99久久精品热视频| 少妇裸体淫交视频免费看高清| 一本久久精品| 特大巨黑吊av在线直播| 欧美变态另类bdsm刘玥| 中文亚洲av片在线观看爽| av国产免费在线观看| 国产一区二区在线av高清观看| 亚洲高清免费不卡视频| 久久99热6这里只有精品| 99热精品在线国产| h日本视频在线播放| 黄片无遮挡物在线观看| 国产av不卡久久| 18禁在线无遮挡免费观看视频| 毛片女人毛片| 午夜精品国产一区二区电影 | 国产精品国产高清国产av| 国产三级中文精品| 亚洲av免费在线观看| 校园人妻丝袜中文字幕| 国产亚洲精品久久久com| 人妻久久中文字幕网| 国产午夜精品论理片| 91午夜精品亚洲一区二区三区| 日韩三级伦理在线观看| 亚洲五月天丁香| 嘟嘟电影网在线观看| 国产伦精品一区二区三区四那| 日本黄色视频三级网站网址| 久久精品久久久久久噜噜老黄 | 99久久无色码亚洲精品果冻| 午夜爱爱视频在线播放| 日本三级黄在线观看| 国产伦在线观看视频一区| 亚洲激情五月婷婷啪啪| 久久国产乱子免费精品| 国产精品爽爽va在线观看网站| 日本一二三区视频观看| 99热全是精品| 大又大粗又爽又黄少妇毛片口| 国产精品一区二区三区四区免费观看| av天堂在线播放| eeuss影院久久| 亚洲av成人av| 中国美女看黄片| 秋霞在线观看毛片| 真实男女啪啪啪动态图| 成人永久免费在线观看视频| 精品一区二区三区人妻视频| 婷婷六月久久综合丁香| 精品免费久久久久久久清纯| 99热这里只有是精品在线观看| 国产av麻豆久久久久久久| 欧美区成人在线视频| 精品久久久久久久久亚洲| 我要搜黄色片| 最近2019中文字幕mv第一页| 精品久久久久久久久久免费视频| 少妇高潮的动态图| 日本免费一区二区三区高清不卡| 日韩欧美三级三区| 午夜老司机福利剧场| 欧美极品一区二区三区四区| 久久久久久久久久成人| 久久欧美精品欧美久久欧美| 长腿黑丝高跟| 一本一本综合久久| 国产综合懂色| av天堂中文字幕网| 久久精品国产自在天天线| 嫩草影院新地址| av国产免费在线观看| 国内精品宾馆在线| 日本黄色片子视频| 熟妇人妻久久中文字幕3abv| 国产午夜福利久久久久久| 26uuu在线亚洲综合色| 有码 亚洲区| 最近手机中文字幕大全| 一区二区三区免费毛片| 亚洲国产色片| 精品人妻一区二区三区麻豆| 国产白丝娇喘喷水9色精品| 日日撸夜夜添| 亚洲欧美中文字幕日韩二区| 久久久久久久久中文| 天天躁夜夜躁狠狠久久av| 亚洲自偷自拍三级| 最近手机中文字幕大全| 久久久a久久爽久久v久久| 国产精品1区2区在线观看.| 看非洲黑人一级黄片| 69人妻影院| 99热这里只有是精品在线观看| 日韩成人伦理影院| 亚洲欧美精品专区久久| 听说在线观看完整版免费高清| 麻豆精品久久久久久蜜桃| kizo精华| 免费看a级黄色片| 中国国产av一级| 美女高潮的动态| 日韩成人伦理影院| 村上凉子中文字幕在线| 免费大片18禁| 99精品在免费线老司机午夜| 亚洲欧洲国产日韩| 又爽又黄无遮挡网站| 三级国产精品欧美在线观看| 18禁在线无遮挡免费观看视频| 18+在线观看网站| 丰满的人妻完整版| 亚洲欧美精品综合久久99| 成人毛片a级毛片在线播放| 亚洲av熟女| 国产伦精品一区二区三区四那| 搡女人真爽免费视频火全软件| 亚洲成人久久爱视频| 三级男女做爰猛烈吃奶摸视频| 久久久精品94久久精品| 精品久久久久久久久久久久久| 日韩一本色道免费dvd| 中国美女看黄片| 国产成人一区二区在线| 欧美性猛交黑人性爽| 亚洲国产精品国产精品| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久精品电影小说 | 成人亚洲精品av一区二区| 国产女主播在线喷水免费视频网站 | 国产午夜精品久久久久久一区二区三区| 一本久久中文字幕| 国语自产精品视频在线第100页| 久久久久久国产a免费观看| 级片在线观看| 婷婷六月久久综合丁香| 午夜福利成人在线免费观看| 最好的美女福利视频网| 午夜激情欧美在线| 日韩欧美 国产精品| 夜夜爽天天搞| 黄色一级大片看看| 黄色视频,在线免费观看| 亚洲欧美日韩高清在线视频| 美女黄网站色视频| 尤物成人国产欧美一区二区三区| 少妇被粗大猛烈的视频| 一级av片app| 91午夜精品亚洲一区二区三区| 亚洲精品国产av成人精品| 成人高潮视频无遮挡免费网站| 国产成人精品婷婷| 精品久久久久久久久久免费视频| 久久精品久久久久久噜噜老黄 | 小蜜桃在线观看免费完整版高清| 99热只有精品国产| 国产综合懂色| 91麻豆精品激情在线观看国产| 如何舔出高潮| 亚洲精品日韩av片在线观看| 熟女人妻精品中文字幕| av在线蜜桃| 97热精品久久久久久| 亚洲,欧美,日韩| 久久久久久久午夜电影| 国产精品野战在线观看| 欧美精品国产亚洲| 99热这里只有是精品50| 国产精品av视频在线免费观看| 啦啦啦啦在线视频资源| 自拍偷自拍亚洲精品老妇| 麻豆成人午夜福利视频| 人妻夜夜爽99麻豆av| 一边摸一边抽搐一进一小说| 国产白丝娇喘喷水9色精品| 一级毛片久久久久久久久女| 久久久久免费精品人妻一区二区| 激情 狠狠 欧美| 欧美一级a爱片免费观看看| 菩萨蛮人人尽说江南好唐韦庄 | 久久久精品大字幕| 91麻豆精品激情在线观看国产| 亚洲人成网站在线播| 久久人人精品亚洲av| 少妇人妻精品综合一区二区 | 能在线免费看毛片的网站| 九草在线视频观看| 岛国在线免费视频观看| 国产探花极品一区二区| 国产私拍福利视频在线观看| 国产精品嫩草影院av在线观看| 欧美最黄视频在线播放免费| 亚洲精品自拍成人| 欧美性感艳星| 欧美色视频一区免费| 天堂av国产一区二区熟女人妻| 亚洲欧美精品自产自拍| 丰满的人妻完整版| 最后的刺客免费高清国语| 国产成人精品一,二区 | 91麻豆精品激情在线观看国产| 国产美女午夜福利| 人妻制服诱惑在线中文字幕| 国产精品一及| 热99re8久久精品国产| 欧美激情久久久久久爽电影| 我的女老师完整版在线观看| 国产黄色视频一区二区在线观看 | 黄色一级大片看看| 黑人高潮一二区| 欧美丝袜亚洲另类| 久久鲁丝午夜福利片| 亚洲av免费在线观看| 韩国av在线不卡| 午夜福利高清视频| 久久午夜亚洲精品久久| 寂寞人妻少妇视频99o| 午夜免费激情av| 在线a可以看的网站| 99久久无色码亚洲精品果冻| 国产精品一区二区在线观看99 | 久久这里只有精品中国| 成年免费大片在线观看| 全区人妻精品视频| kizo精华| 一区二区三区高清视频在线| 亚洲av一区综合| 天堂影院成人在线观看| 久久99精品国语久久久| 18禁裸乳无遮挡免费网站照片| av专区在线播放| 2021天堂中文幕一二区在线观| 性插视频无遮挡在线免费观看| 午夜精品一区二区三区免费看| 永久网站在线| 久久精品影院6| 嘟嘟电影网在线观看| 国产精品久久久久久久电影| 国产一区二区三区在线臀色熟女| 内地一区二区视频在线| 天天躁夜夜躁狠狠久久av| 91久久精品国产一区二区成人| 日韩一区二区视频免费看| 精品久久久久久久久亚洲| 男人和女人高潮做爰伦理| 日韩,欧美,国产一区二区三区 | 亚洲自拍偷在线| 能在线免费看毛片的网站| 中文资源天堂在线| 男女做爰动态图高潮gif福利片| 精品日产1卡2卡| 五月伊人婷婷丁香| 欧美另类亚洲清纯唯美| 99久久精品热视频| 国产在线精品亚洲第一网站| 天堂√8在线中文| 大又大粗又爽又黄少妇毛片口| 国产精品国产三级国产av玫瑰| 午夜久久久久精精品| 我的女老师完整版在线观看| 日韩大尺度精品在线看网址| 成人特级av手机在线观看| 在线国产一区二区在线| 国产黄色小视频在线观看| 欧美+日韩+精品| 亚洲av一区综合| 国产成人freesex在线| 久久6这里有精品| 看十八女毛片水多多多| 亚洲精品久久国产高清桃花| 亚洲中文字幕一区二区三区有码在线看| 久久国内精品自在自线图片| 国内少妇人妻偷人精品xxx网站| 亚洲人与动物交配视频| 久久99精品国语久久久| 成熟少妇高潮喷水视频| 99热网站在线观看| 亚洲欧美日韩高清在线视频| 成人特级av手机在线观看| 97超碰精品成人国产| 青青草视频在线视频观看| 伊人久久精品亚洲午夜| 一区二区三区免费毛片| 精品一区二区三区人妻视频| 99热6这里只有精品| 久久人人精品亚洲av| 特大巨黑吊av在线直播| 久久精品国产清高在天天线| 午夜精品一区二区三区免费看| 少妇的逼好多水| 国产成人午夜福利电影在线观看| 最近中文字幕高清免费大全6| 少妇高潮的动态图| 午夜福利视频1000在线观看| 日韩欧美 国产精品| 一级黄片播放器| 国产三级在线视频| 蜜桃久久精品国产亚洲av| 日韩制服骚丝袜av| 中文字幕久久专区| 日韩强制内射视频| 欧美激情久久久久久爽电影| 一本久久中文字幕| 国产乱人视频| 91精品国产九色| 九九久久精品国产亚洲av麻豆| 国产色爽女视频免费观看| 亚洲国产欧洲综合997久久,| 日韩在线高清观看一区二区三区| 亚洲av熟女| 在线观看免费视频日本深夜| 午夜爱爱视频在线播放| 国产真实伦视频高清在线观看| 久久久精品大字幕| 久久久成人免费电影| 亚洲人成网站高清观看| 日韩强制内射视频| 婷婷亚洲欧美| 午夜精品在线福利| 亚洲精品亚洲一区二区| 一级av片app| 九九久久精品国产亚洲av麻豆| 国产精品麻豆人妻色哟哟久久 | 99热网站在线观看| 国产av麻豆久久久久久久| 99九九线精品视频在线观看视频| 国产黄片美女视频| 久久久国产成人免费| 男人舔奶头视频| 男人狂女人下面高潮的视频| 亚洲av.av天堂| 啦啦啦观看免费观看视频高清| 久久久久久大精品| 老熟妇乱子伦视频在线观看| 一本一本综合久久| 久久99蜜桃精品久久| 91精品国产九色| 久久久久久伊人网av| 免费人成在线观看视频色| 99九九线精品视频在线观看视频| 人人妻人人澡人人爽人人夜夜 | 亚洲av成人av| 伦精品一区二区三区| а√天堂www在线а√下载| av在线亚洲专区| 激情 狠狠 欧美| 欧美性感艳星| 又爽又黄a免费视频| 亚洲真实伦在线观看| 国产高清有码在线观看视频| 欧美潮喷喷水| 国产蜜桃级精品一区二区三区| 精品一区二区三区视频在线| 少妇高潮的动态图| 一区二区三区四区激情视频 | 亚洲无线观看免费| 国产成人福利小说| 欧美一级a爱片免费观看看| 色吧在线观看| 国产在视频线在精品| 精品一区二区免费观看| 久久精品国产亚洲av天美| 亚洲成人中文字幕在线播放| 国产91av在线免费观看| 亚洲国产欧洲综合997久久,| 亚洲丝袜综合中文字幕| 国产国拍精品亚洲av在线观看| 国产视频内射| 日韩大尺度精品在线看网址| 亚洲欧美成人综合另类久久久 | 亚洲最大成人手机在线| 看十八女毛片水多多多| 国产一区二区亚洲精品在线观看| 国产精品人妻久久久久久| 在线观看美女被高潮喷水网站| 国产伦在线观看视频一区| 岛国在线免费视频观看| 网址你懂的国产日韩在线| 男女视频在线观看网站免费| 免费av不卡在线播放| 春色校园在线视频观看| 国产女主播在线喷水免费视频网站 | 一个人看视频在线观看www免费| 嫩草影院精品99| 中文字幕久久专区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人午夜高清在线视频| 亚洲欧美成人精品一区二区| 国产精品三级大全| 国产探花在线观看一区二区| 一级毛片电影观看 | 国产色爽女视频免费观看| 国产成年人精品一区二区| 国内精品久久久久精免费| 99在线人妻在线中文字幕| 欧美潮喷喷水| 欧美性猛交╳xxx乱大交人| 国产成人一区二区在线| 91在线精品国自产拍蜜月| 人妻制服诱惑在线中文字幕| 青青草视频在线视频观看| 中国美白少妇内射xxxbb| 黄色视频,在线免费观看| 色5月婷婷丁香| 亚洲av第一区精品v没综合| 国产精品乱码一区二三区的特点| 狠狠狠狠99中文字幕| 少妇熟女aⅴ在线视频| av在线播放精品| 99久久九九国产精品国产免费| 免费搜索国产男女视频| 夫妻性生交免费视频一级片| 国产精品久久视频播放| 在线观看免费视频日本深夜| 美女cb高潮喷水在线观看| 国产精品久久久久久久电影| 99久国产av精品国产电影| 欧美精品一区二区大全| 狂野欧美白嫩少妇大欣赏| 男人舔奶头视频| 欧美丝袜亚洲另类| 欧美日韩在线观看h| 久久精品国产99精品国产亚洲性色| 精品国产三级普通话版| a级一级毛片免费在线观看| 毛片一级片免费看久久久久| 日韩av不卡免费在线播放| 直男gayav资源| 欧美一区二区国产精品久久精品| av福利片在线观看| 亚洲成a人片在线一区二区| 国产精品国产三级国产av玫瑰| 欧美+日韩+精品| 欧美色视频一区免费| 色噜噜av男人的天堂激情| 一级毛片aaaaaa免费看小| 亚洲av熟女| 日韩av不卡免费在线播放| 国产乱人偷精品视频| 99九九线精品视频在线观看视频| 成年女人看的毛片在线观看| 晚上一个人看的免费电影| 国产精品三级大全| 一级黄色大片毛片| 色视频www国产| 国产精品一区二区三区四区久久| 久久精品国产亚洲av涩爱 | 久久99蜜桃精品久久| 亚洲在线自拍视频| 99国产精品一区二区蜜桃av| 男的添女的下面高潮视频| 国产精品永久免费网站| 偷拍熟女少妇极品色| 国产精品一区二区性色av| 99热这里只有精品一区| 国产一区二区三区在线臀色熟女| 你懂的网址亚洲精品在线观看 | 黄色一级大片看看| 少妇人妻精品综合一区二区 | 亚洲人成网站在线播| a级毛片免费高清观看在线播放| 免费无遮挡裸体视频| 欧美xxxx性猛交bbbb| 我的女老师完整版在线观看| 伊人久久精品亚洲午夜| 日韩 亚洲 欧美在线| 高清日韩中文字幕在线| 免费看光身美女| 最近2019中文字幕mv第一页| 久久中文看片网| 看免费成人av毛片| 狂野欧美白嫩少妇大欣赏| 国产精品免费一区二区三区在线| 欧美激情国产日韩精品一区| 人妻少妇偷人精品九色| 亚洲av男天堂| 亚洲最大成人中文| 午夜激情欧美在线| 深夜精品福利| 成年版毛片免费区| 桃色一区二区三区在线观看| 又粗又硬又长又爽又黄的视频 | 狠狠狠狠99中文字幕| 亚洲精品乱码久久久久久按摩| 国产色爽女视频免费观看| 欧美zozozo另类| 国产亚洲5aaaaa淫片| 青春草视频在线免费观看| 亚洲精品影视一区二区三区av| 黄片无遮挡物在线观看| 国产伦精品一区二区三区四那|