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

    FLAC3D中巖石能量耗散模型的開發(fā)與應用

    2021-09-01 00:52:58宋子楓鄭冬杰神文龍勾攀峰韋四江
    煤炭學報 2021年8期
    關鍵詞:圍巖變形模型

    王 猛,宋子楓,鄭冬杰,神文龍,勾攀峰,韋四江,2

    (1.河南理工大學 能源科學與工程學院,河南 焦作 454000; 2.中國礦業(yè)大學 深部煤炭資源開采教育部重點實驗室,徐州 江蘇 221116; 3.煤炭安全生產河南省協同創(chuàng)新中心,河南 焦作 454000)

    準確描述巖石強度與變形是評價工程巖體安全和穩(wěn)定性的基礎,現有應力強度理論與破壞準則很難有效描述巖石復雜的強度變化與整體破壞行為[1]。巖石破壞是能量驅動下的一種狀態(tài)失穩(wěn)現象[2],受載巖石變形過程伴隨能量演化,峰前以能量積聚為主,峰后主要表現為能量耗散與釋放,巖石變形破壞可看成是不同形式能量相互轉換的結果。因此,從能量角度研究巖石變形和失穩(wěn),更接近其破壞本質[3-5]。

    基于能量守恒理論,國內外學者對巖石加卸載過程能量演化規(guī)律開展了大量研究,取得了豐富的研究成果。具有代表性的有:謝和平等[6-8]闡述了巖石變形破壞過程能量耗散與損傷、能量釋放與整體破壞的關系,定義了基于能量耗散與釋放原理的巖石強度失效準則和巖體整體破壞準則。張志鎮(zhèn)等[9-11]研究了不同巖石類型能量演化的異同點,建立了巖石能量演化隨加載條件的自我抑制模型,揭示了巖石能量耗散的圍壓效應。宮鳳強、李淼等[12-13]分別研究了動靜態(tài)劈裂過程巖石能量構成與耗散特征,揭示了動靜荷載作用下巖石能量耗散機制及各向異性特征。與此同時,許江等[14]、楊磊等[15]、尹光志和蘇國韶[16-17]、李楊楊等[18]分別研究了含孔隙水砂巖、煤巖組合體、不同加載速率以及循環(huán)加卸載等條件下巖石能量積聚與耗散特征,豐富了巖石能量理論的研究成果。

    現有研究常聚焦于室內試驗,多是針對特定加卸載條件、巖樣類型以及其他條件下(含水率、溫度[19]等)巖石能量演化特征的描述,研究對象著重關注巖石屈服、峰值等特征點,對于貫穿巖石變形破壞全過程的能量實時演化及分配規(guī)律認識不足[9]。巖石能量積聚、耗散及其狀態(tài)間的轉化決定了其內部裂隙擴展和變形破壞,雖然從能量角度反演巖石破壞路徑對實驗結果意義不大,但對于評估工程巖體區(qū)域穩(wěn)定、精準定位加固位置至關重要。目前,實驗室常用巖石破壞監(jiān)測手段有聲發(fā)射定位[20],CT掃描[21],SEM電鏡掃描[22],CCD相機[23]等,其中,僅有聲發(fā)射可以實時定位受載巖樣變形破壞路徑,但該方法多用于單軸壓縮試驗,且存在操作復雜、誤差大等缺點,很難推廣應用于工程巖體的監(jiān)測。分析工程巖體穩(wěn)定性時多是基于圍巖應力場、位移場及支護場的耦合作用關系進行綜合評判。

    為增加能量理論模型分析工程問題時的適用性,筆者基于能量守恒和有限差分理論,推導巖石能量耗散有限差分方程式,采用FISH語言將其應用于FLAC3D應變軟化本構模型。以室內試驗結果驗證采用耗散能演化描述受載巖石變形破壞的合理性,進一步討論能量耗散模型在深部高應力巷道中的應用,揭示工程巖體峰后變形破壞全過程耗散能演化特征與分配規(guī)律,為從能量耗散角度評估工程巖體安全和穩(wěn)定性提供理論依據,彌補FLAC3D軟件采用塑性區(qū)定性判斷地下工程圍巖破壞時的不足。

    1 巖石能量耗散的室內實驗

    1.1 巖石單、三軸壓縮試驗

    試驗巖樣取自徐礦集團三河尖煤礦吳莊區(qū)運輸大巷,選取巖性均勻、結構完整的取芯巖柱加工巖樣,采用GCTS-RTX3000巖石力學試驗系統進行單軸和常規(guī)三軸壓縮試驗,測試巖石基本力學參數,圍壓等級設置為0,5,15和25 MPa,測試結果分別如圖1和表1所示。

    表1 常規(guī)單、三軸試驗結果Table 1 Results of uniaxial and triaxial tests

    圖1 巖石應力-應變曲線Fig.1 Stress-strain curves of the rock samples

    1.2 巖石能量耗散特征

    巖樣在試驗機外力作用下產生變形,假設該物理過程與外界沒有熱交換,即封閉系統,外力功所產生的總輸入能量為Wz,根據熱力學第一定律[6]得到

    Wz=We+Wd

    (1)

    式中,We為巖石可釋放彈性應變能;Wd為耗散能。

    由式(1)可知,巖石總輸入能量Wz一部分作為可釋放的彈性能儲存在其內部,另一部分則伴隨巖石變形破壞耗散和釋放,此過程不可逆[8-9]。巖石損傷程度越高,耗散能越大,殘余彈性能越少。

    對于單軸壓縮試驗,試樣總輸入能量Wz等于軸向輸入能量Wa;而常規(guī)三軸壓縮試驗,軸向輸入能量Wa應包含試樣總輸入能量Wz以及環(huán)向膨脹對液壓油做功所釋放能量Wh兩部分,為此,常三軸壓縮條件下巖石總輸入能量[6]可表示為

    (2)

    式中,ε1,ε3分別為軸向應變和環(huán)向應變;σ1,σ3分別為軸向應力和圍壓。

    根據彈性力學相關理論,巖石可釋放彈性能We[6]可由下式計算:

    (3)

    則巖石耗散能Wd可由下式求得

    Wd=Wz-We

    (4)

    采用式(2)~(4)分別計算巖樣總輸入能量Wz、彈性能We和耗散能Wd,如圖2所示。由圖2可知,巖樣總輸入能量隨軸向應變增加逐漸增大,應力達到峰值前主要以彈性能形式儲存在巖樣內部,耗散能變化不大。當巖樣加載至峰后階段,巖樣總輸入能量依然持續(xù)增加,但主要伴隨峰后破壞而耗散,此外,巖樣內部儲存的大部分彈性能也將轉化成耗散能釋放,能量轉化程度與巖樣變形破壞程度正相關。

    圖2 室內試驗巖樣能量演化規(guī)律Fig.2 Evolution law of energy of rock samples in laboratory test

    2 FLAC3D能量耗散模型的開發(fā)

    σ1<σ2<σ3

    (5)

    相應的主應變增量Δεi[25]可分解為

    (6)

    單元體彈性應變增量和應力增量之間的關系[25]可表示為

    (7)

    (8)

    式中,α1,α2為切變模量G和體積模量K控制的巖石材料常數,其中,α1=K+4G/3,α2=K-2G/3。

    基于Mohr-Coulomb的應變軟化模型在峰后軟化過程控制強度參數描述材料的峰后軟化行為,對于剪切破壞,非關聯流動法則[25]改寫為

    (9)

    式中,λs為塑乘因子;gs為剪切勢函數,其表達式為

    (10)

    式中,ψ為剪脹角。

    因此,聯立式(6),(7)和(9),由總應變表示的應力增量Δσi表達式[25]為

    (12)

    (13)

    式中,φ為內摩擦角。

    根據FLAC3D應變軟化模型屈服準則,巖體材料剪切破壞面上應力[25]滿足:

    (14)

    式中,c為黏聚力。

    拉伸破壞面上應力[25]滿足:

    ft=σ3-σt=0

    (15)

    式中,σt為抗拉強度。

    式(14)和(15)分別為FLAC3D應變軟化模型判定材料壓剪破壞和張拉破壞的失穩(wěn)準則。

    圖3 單元能量計算簡圖Fig.3 Schematic diagram for the energy calculation

    (16)

    假設t時刻對應第n循環(huán),t+Δt時刻對應第n+1循環(huán),則第n+1循環(huán)內總能量增量可表示為

    (17)

    第n循環(huán)至第n+1循環(huán)運算的物理意義表示巖體受外力作用時,某一運算時步內單元網格應力-應變行為的差分運算。t+Δt時刻單元總能量等于n+1循環(huán)總能量增量累加,即

    (18)

    不考慮峰后彈性模量衰減,單元可釋放彈性應變能可由下式計算:

    (19)

    則單元耗散能計算式如下:

    (20)

    圖4 單元能量計算流程Fig.4 Energy calculation flow for the model cells

    3 模型校驗與運算

    3.1 巖樣模型建立

    建立標準巖樣數值模型,尺寸φ×H=50 mm×100 mm,共劃分17 280單元,如圖5所示。模型底部固定位移邊界,四周施加環(huán)向邊界應力模擬圍壓,上邊界通過施加軸向位移對巖樣加壓,加載速率2.5×10-5mm/step,計算時忽略巖石自重影響。

    模型采用基于Mohr-Coulomb的應變軟化本構模型,以表1給出的巖樣力學參數對模型進行初始賦值。巖石峰后軟化階段黏聚力和內摩擦角隨應變的變化數據由試驗結果計算得到,見表2。

    表2 巖石峰后軟化參數Table 2 Post-peak parameters of rock samples

    依據圖5建立的數值模型,表1和2所列巖石力學參數進行數值分析,圖6給出了模擬巖樣應力應變曲線與試驗結果的對比,由圖6可知,不同圍壓下巖樣應力應變曲線模擬結果較為匹配試驗結果,單軸壓縮和低圍壓條件下巖樣脆性特征得到較好描述,表明所選取的巖石參數和數值本構模型合理,可用于后續(xù)能量模塊的模擬分析。

    圖5 巖樣數值模型Fig.5 Numerical model of rock samples

    圖6 應力-應變曲線校驗Fig.6 Verifications of the stress-strain curves

    3.2 巖樣能量模型校驗

    考慮到耗散能是由巖石總能量減去彈性應變能計算得到,限于篇幅,校驗能量模型時只校驗總能量和耗散能兩個指標。如圖7(a)所示,不同圍壓下巖樣總能量數值模擬結果與試驗結果基本吻合,表明FISH語言編寫的能量算法可靠性較強。

    但需要指出的是,模擬巖樣峰后耗散能演化曲線與試驗結果吻合度較好,但峰前存在差異性,如圖7(b)所示。峰前階段,巖樣試驗耗散能主要呈緩慢線性增加;加載至峰值80%時,耗散能近似指數增加;一旦到達峰后某破壞點,耗散能演化為斜率較大的直線增長,最終進入殘余階段。數值模擬可近似反演耗散能演化趨勢,耗散能在峰后增長曲線與試驗結果匹配較好,且總量大致相等;但峰前階段存在顯著差異,模擬峰前耗散能近似等于0,主要原因是能量模型峰前假設不發(fā)生塑性變形,無法準確描述孔隙、松軟巖體初始壓密階段和峰前塑性階段,導致圍壓越大,峰前耗散能相差越大。但是,對于地下工程巖體的分析更關注于峰后變形破壞,由于該模型可較好模擬巖石峰后能量耗散特征,可用于巖石峰后變形破壞的模擬。

    圖7 能量模型校驗Fig.7 Verifications of the energy model

    3.3 巖樣峰后破壞路徑反演

    圖8給出了不同圍壓下巖石破壞模式對比,其中,數值模擬給出的是巖石殘余階段耗散能密度演化云圖,對應試驗巖樣殘余階段破壞照片。由圖8可知,巖樣單軸壓縮試驗呈縱向劈裂破壞,模擬巖樣耗散能集中區(qū)域位于巖樣中線偏上位置,并向周邊演化出一“八”字分支集中線,巖樣耗散能演化沿中線呈軸對稱分布,反映了巖樣單軸壓縮時的縱向劈裂特征,與試驗結果吻合度較高。隨著圍壓增加,試驗巖樣逐漸向剪切破壞演化,破斷角逐漸減小;模擬巖樣在5 MPa時耗散能密度集中成“y”型分布,長邊集中程度顯著大于短邊,為巖樣主破斷面;隨圍壓進一步增加,其破壞模式演化成斜切面,且圍壓越大,耗散能越大,破斷角越小,與試驗結果基本吻合,表明采用耗散能密度演化揭示巖石峰后破壞行為是可行的。

    圖8 巖石破壞模式校驗Fig.8 Verifications of the failure modes of rock samples

    工程現場,準確掌握圍巖變形破壞路徑,精準定位加固位置是控制其穩(wěn)定的前提。為此,反演巖石峰后耗散能演化過程可為揭示其主破壞區(qū)演化、評判巖石穩(wěn)定性提供新途徑。圖9給出了圍壓25 MPa時不同峰后目標點巖樣耗散能密度演化云圖。巖樣初始進入峰后階段(B點),經歷裂隙孕育形成初始主裂隙面,并隨著軸向應變增加,主裂隙面耗散能逐漸增加,一旦巖樣進入殘余階段(F點),主裂隙面耗散能密度增幅減小并趨于穩(wěn)定,表明巖樣此時已發(fā)生宏觀破斷,繼續(xù)加載將導致其結構失穩(wěn)。

    圖9 圍壓25 MPa巖樣峰后破壞路徑Fig.9 Post-peak destruction path of rock sample with a confining pressure of 25 MPa

    4 能量耗散模型的應用

    將能量耗散模型應用于分析吳莊區(qū)運輸大巷,嘗試從能量耗散角度揭示巷道變形破壞路徑,為巷道加固提供依據。建立試驗巷道三維數值模型,尺寸X×Y×Z=60 m×40 m×60 m,水平、底邊界限制位移,上邊界施加20 MPa載荷模擬巷道800 m埋深,依據地應力測試結果,取側壓系數0.8。巷道斷面寬×高=5 m×4 m,模擬巷道錨桿支護參數與現場一致。

    模擬采用應變軟化模型,由于巖石實驗參數無法直接應用于工程巖體,模擬前需校驗巖體參數?;趯嶒灚@取的巖石強度參數衰減規(guī)律,采用RocLab軟件進行巖體參數轉換,以獲取的初始巖體參數進行賦參,以實測巷道變形作為已知值,校驗獲取的巖體參數見表3,4。同時,為防止單元畸變中斷運算,模擬時未設置set large命令。由于巷道周邊圍巖受力多處于三向不等壓狀態(tài),計算單元能量時考慮中間主應力,各應力分量增量計算方法見式(12),再由式(17)~(20)分別計算單元總能量、彈性能和耗散能。

    表3 校驗巖體參數(峰前)Table 3 Calibrated parameters of the rock masses

    圖10給出了試驗巷道模擬結果與現場對比圖。試驗巷道開挖后,采用錨網索支護,巷道頂幫變形得到控制,頂幫變形均在200 mm以內;但底臌嚴重,最大底臌量接近1 000 mm,模擬巷道變形與現場實測數據基本吻合。

    表4 校驗巖體參數(峰后)Table 4 Calibrated parameters of the rock masses

    圖10 試驗巷道變形破壞對比Fig.10 Comparisons on the roadway deformations

    模擬巷道塑性區(qū)主要集中于兩幫和底板,基于軟件拉剪屈服破壞準則,僅可以區(qū)分底幫圍巖變形破壞方式,但很難量化底幫圍巖的破壞程度,也就無法準確定位巷道后期加固區(qū)域。而利用開發(fā)的能量耗散模型,從巷道耗散能集中區(qū)分布可實現對巷道圍巖主失穩(wěn)區(qū)域的實時定位,同時,依據圍巖耗散能集中程度可以量化巖體變形破壞程度。如圖10所示,巷道底板破壞后,底角位置分別存在兩個耗散能集中區(qū),均呈條狀分布,一條貫通與底角,誘發(fā)底角剪切破壞;另一條則向底板中部演化,加劇底臌。

    基于3.3節(jié),研究開發(fā)的能量耗散模型除可定位圍巖主破壞位置,還可以通過控制運算時步反演巷道變形破壞路徑。圖11給出了不同時步對應的巷道圍巖耗散能密度演化云圖,運算平衡時對應的耗散能云圖如圖10所示。由圖11可知,巷道開挖后耗散能集中區(qū)位于巷道兩底角處,表明直角相比于拱形更易產生剪切破壞,這與前人研究結論相吻合。隨著時間增加,巷道底角破壞后,耗散能集中區(qū)域將向深部圍巖轉移,并逐漸演化成帶狀分布,一端向巷道兩幫轉移,一端朝向巷道底板方向,誘發(fā)巷道底臌和兩幫收斂。對于破壞后的巷道,如對底幫尤其是底角圍巖不加以控制,劇烈底臌將進一步引發(fā)巷幫松動破壞,甚至危及頂板穩(wěn)定,最終誘發(fā)巷道整體災變失穩(wěn)。

    圖11 巷道破壞路徑演化模擬Fig.11 Simulated results of the failure paths of roadway

    對于試驗巷道,后期采用了巷幫短錨索和底板錨桿(索)加固,巷道大變形得到有效控制。限于篇幅,具體參數不再贅述。綜上,研究開發(fā)的能量耗散模型對于反演巷道變形破壞路徑和定位巷道主破壞位置提供了一條新途徑,可為巷道失穩(wěn)預測預報和加固設計提供了參考,具有較強的實用價值。

    5 結 論

    (1)基于能量平衡和有限差分理論,推導了巖石耗散能有限差分方程式,采用FISH語言將其寫入FLAC3D應變軟化模型,補充了軟件能量計算模塊。通過與室內試驗對比,該模型可有效描述巖石峰后變形特征及破壞路徑。

    (2)利用開發(fā)的能量耗散模型分析了深部巷道能量演化規(guī)律,采用耗散能演化不僅可定位巷道主破壞位置,同時可有效反演圍巖全過程變形破壞路徑,為巷道變形預測及加固設計提供了一條新途徑。

    (3)巖石能量耗散模型在FLAC3D中的實現,擴展了FLAC3D軟件的適用范圍,一定程度上彌補了常規(guī)采用塑性區(qū)定性判斷地下工程圍巖破壞及其穩(wěn)定的不足,具有重要的理論意義和現場實用價值。

    猜你喜歡
    圍巖變形模型
    一半模型
    重要模型『一線三等角』
    談詩的變形
    中華詩詞(2020年1期)2020-09-21 09:24:52
    重尾非線性自回歸模型自加權M-估計的漸近分布
    隧道開挖圍巖穩(wěn)定性分析
    中華建設(2019年12期)2019-12-31 06:47:58
    “我”的變形計
    軟弱破碎圍巖隧道初期支護大變形治理技術
    江西建材(2018年4期)2018-04-10 12:37:22
    例談拼圖與整式變形
    會變形的餅
    3D打印中的模型分割與打包
    看十八女毛片水多多多| 成人亚洲精品一区在线观看 | 日韩av在线免费看完整版不卡| 1000部很黄的大片| 亚洲欧美一区二区三区黑人 | 亚洲精华国产精华液的使用体验| 赤兔流量卡办理| 欧美日本视频| 赤兔流量卡办理| 日本午夜av视频| 亚洲怡红院男人天堂| 国产欧美另类精品又又久久亚洲欧美| 中国国产av一级| 国产免费视频播放在线视频 | 国产成人freesex在线| 欧美区成人在线视频| 久久热精品热| 国产真实伦视频高清在线观看| 亚洲怡红院男人天堂| 一级a做视频免费观看| 久久6这里有精品| 国内精品宾馆在线| 国产视频内射| 熟妇人妻不卡中文字幕| 夜夜爽夜夜爽视频| 伊人久久精品亚洲午夜| 久久久久久久久久人人人人人人| 边亲边吃奶的免费视频| 亚洲国产精品sss在线观看| 日韩精品有码人妻一区| 日韩欧美三级三区| 亚洲最大成人av| 亚洲高清免费不卡视频| 在线 av 中文字幕| 日本黄色片子视频| 国产在线一区二区三区精| 大又大粗又爽又黄少妇毛片口| 欧美变态另类bdsm刘玥| 久久鲁丝午夜福利片| 美女脱内裤让男人舔精品视频| 成人毛片60女人毛片免费| 国产精品国产三级专区第一集| 久久精品夜夜夜夜夜久久蜜豆| 久久人人爽人人爽人人片va| 最近视频中文字幕2019在线8| 亚洲国产精品成人综合色| 美女cb高潮喷水在线观看| 国产麻豆成人av免费视频| 晚上一个人看的免费电影| 亚洲天堂国产精品一区在线| 80岁老熟妇乱子伦牲交| 特级一级黄色大片| 精品国产露脸久久av麻豆 | 色哟哟·www| 校园人妻丝袜中文字幕| videossex国产| 51国产日韩欧美| 久热久热在线精品观看| 亚洲美女搞黄在线观看| 午夜免费观看性视频| 亚洲欧美日韩东京热| 夜夜看夜夜爽夜夜摸| 亚洲av国产av综合av卡| 在线观看一区二区三区| 中文精品一卡2卡3卡4更新| 少妇裸体淫交视频免费看高清| 婷婷色综合www| 国产精品国产三级国产专区5o| 蜜桃久久精品国产亚洲av| 精品一区二区三卡| 一个人看视频在线观看www免费| 精品不卡国产一区二区三区| 男女啪啪激烈高潮av片| 日韩伦理黄色片| 国产欧美另类精品又又久久亚洲欧美| 99九九线精品视频在线观看视频| 国产成人a∨麻豆精品| 男人爽女人下面视频在线观看| 婷婷色av中文字幕| 日本一本二区三区精品| 亚洲国产日韩欧美精品在线观看| 亚洲成色77777| ponron亚洲| 国产一区二区在线观看日韩| 精品午夜福利在线看| 偷拍熟女少妇极品色| 一本久久精品| 少妇熟女欧美另类| 伊人久久国产一区二区| 赤兔流量卡办理| 成人亚洲精品一区在线观看 | 午夜福利视频精品| 精品一区二区三区视频在线| 国产精品久久久久久久久免| 国产成人精品一,二区| 亚洲最大成人中文| 日韩一本色道免费dvd| 色综合亚洲欧美另类图片| 噜噜噜噜噜久久久久久91| 久久久久九九精品影院| 久久久久久久大尺度免费视频| 亚洲国产精品国产精品| 一个人免费在线观看电影| 久久精品久久精品一区二区三区| 免费观看无遮挡的男女| 日韩 亚洲 欧美在线| 欧美3d第一页| 乱码一卡2卡4卡精品| 亚洲av一区综合| 免费观看性生交大片5| 久99久视频精品免费| 欧美极品一区二区三区四区| 国产成人精品一,二区| 久久精品国产亚洲av涩爱| 三级国产精品欧美在线观看| 中文字幕人妻熟人妻熟丝袜美| 天堂网av新在线| 亚洲最大成人中文| 国产成人精品久久久久久| 街头女战士在线观看网站| 亚洲综合色惰| 免费看美女性在线毛片视频| 精品久久久精品久久久| 熟女人妻精品中文字幕| 亚洲精品乱码久久久v下载方式| 久久这里有精品视频免费| 欧美激情久久久久久爽电影| 亚洲av日韩在线播放| 久久精品国产自在天天线| 欧美精品国产亚洲| 日日干狠狠操夜夜爽| 在线 av 中文字幕| 久久久久久久大尺度免费视频| 成人欧美大片| 看免费成人av毛片| 日韩精品青青久久久久久| 亚洲伊人久久精品综合| 精品一区在线观看国产| 久久久久九九精品影院| 嘟嘟电影网在线观看| 99久久精品一区二区三区| 老司机影院毛片| 国产极品天堂在线| 国产欧美另类精品又又久久亚洲欧美| 天美传媒精品一区二区| 简卡轻食公司| av黄色大香蕉| 免费看光身美女| 欧美精品一区二区大全| 精品国产三级普通话版| 激情 狠狠 欧美| 大香蕉97超碰在线| 欧美极品一区二区三区四区| 80岁老熟妇乱子伦牲交| 插逼视频在线观看| eeuss影院久久| 国产亚洲最大av| 国产精品熟女久久久久浪| 极品教师在线视频| 午夜视频国产福利| 天堂中文最新版在线下载 | 亚洲人成网站在线观看播放| 少妇高潮的动态图| 久久精品人妻少妇| 日日撸夜夜添| av卡一久久| 日本黄色片子视频| 成人午夜精彩视频在线观看| 精品国产露脸久久av麻豆 | 少妇猛男粗大的猛烈进出视频 | 蜜桃亚洲精品一区二区三区| 成人午夜精彩视频在线观看| 高清午夜精品一区二区三区| 亚洲欧洲国产日韩| 日本av手机在线免费观看| 亚洲欧美一区二区三区国产| 午夜精品一区二区三区免费看| 成年版毛片免费区| 国产三级在线视频| 丰满乱子伦码专区| 禁无遮挡网站| 熟妇人妻不卡中文字幕| 免费看光身美女| 国产三级在线视频| 久久精品综合一区二区三区| 免费观看的影片在线观看| 中文字幕免费在线视频6| 国产亚洲精品av在线| 男的添女的下面高潮视频| 日韩欧美三级三区| 男女视频在线观看网站免费| 亚洲一区高清亚洲精品| 中文字幕av在线有码专区| 亚洲欧美一区二区三区黑人 | 高清午夜精品一区二区三区| 九九久久精品国产亚洲av麻豆| 一区二区三区四区激情视频| 亚洲成人一二三区av| 看黄色毛片网站| 国产精品综合久久久久久久免费| 亚洲自偷自拍三级| 婷婷色综合www| av黄色大香蕉| 波野结衣二区三区在线| 又大又黄又爽视频免费| 免费观看无遮挡的男女| 如何舔出高潮| 天天一区二区日本电影三级| 免费高清在线观看视频在线观看| 国产男人的电影天堂91| videossex国产| 99久久人妻综合| 亚洲经典国产精华液单| 日韩伦理黄色片| 久久6这里有精品| 欧美日韩亚洲高清精品| 少妇熟女aⅴ在线视频| 亚洲精品成人久久久久久| 我要看日韩黄色一级片| 亚洲人与动物交配视频| 欧美变态另类bdsm刘玥| 久久热精品热| 夫妻午夜视频| 国产中年淑女户外野战色| 午夜免费观看性视频| 纵有疾风起免费观看全集完整版 | 久久久色成人| 亚洲国产精品专区欧美| 国产黄片视频在线免费观看| 男人狂女人下面高潮的视频| 日韩成人伦理影院| 亚洲av中文av极速乱| 搡老乐熟女国产| 人妻系列 视频| 午夜免费男女啪啪视频观看| 男女边摸边吃奶| 2018国产大陆天天弄谢| 精品一区二区三区视频在线| 在线观看人妻少妇| 肉色欧美久久久久久久蜜桃 | 欧美另类一区| 最近的中文字幕免费完整| 亚洲在线自拍视频| 99视频精品全部免费 在线| www.色视频.com| 午夜老司机福利剧场| 国产色爽女视频免费观看| 亚洲第一区二区三区不卡| 亚洲一级一片aⅴ在线观看| 亚洲精品国产av成人精品| 国产精品麻豆人妻色哟哟久久 | 午夜免费激情av| 国产免费又黄又爽又色| 国产成人精品久久久久久| 久久久精品94久久精品| 美女高潮的动态| 日韩国内少妇激情av| 亚洲欧洲国产日韩| 午夜老司机福利剧场| 国产一区有黄有色的免费视频 | 国产男人的电影天堂91| 国产成人aa在线观看| 午夜免费男女啪啪视频观看| 国产黄频视频在线观看| 国产综合精华液| 欧美极品一区二区三区四区| 看十八女毛片水多多多| 日韩大片免费观看网站| 在线 av 中文字幕| 99久久精品一区二区三区| 午夜日本视频在线| 欧美一区二区亚洲| 又爽又黄无遮挡网站| 精品酒店卫生间| 精品久久久久久电影网| 亚洲欧美清纯卡通| 啦啦啦韩国在线观看视频| 国产爱豆传媒在线观看| 国国产精品蜜臀av免费| 啦啦啦啦在线视频资源| 九草在线视频观看| 久久这里只有精品中国| 18禁裸乳无遮挡免费网站照片| 色尼玛亚洲综合影院| 你懂的网址亚洲精品在线观看| 亚洲av免费在线观看| 好男人在线观看高清免费视频| 亚洲天堂国产精品一区在线| 日本-黄色视频高清免费观看| 久久久国产一区二区| 国产视频内射| 看免费成人av毛片| 99久久精品热视频| 插阴视频在线观看视频| 乱码一卡2卡4卡精品| 麻豆久久精品国产亚洲av| 97在线视频观看| 国产一级毛片在线| 99久久九九国产精品国产免费| 成人国产麻豆网| 好男人视频免费观看在线| 国产成人a区在线观看| 精品人妻偷拍中文字幕| 噜噜噜噜噜久久久久久91| 亚洲精品中文字幕在线视频 | 欧美激情国产日韩精品一区| 日韩中字成人| 国语对白做爰xxxⅹ性视频网站| 国产免费又黄又爽又色| 精品人妻偷拍中文字幕| 偷拍熟女少妇极品色| 丝瓜视频免费看黄片| av在线蜜桃| 国产精品爽爽va在线观看网站| 久久草成人影院| 精品国产一区二区三区久久久樱花 | 亚洲精品亚洲一区二区| 成人漫画全彩无遮挡| 国产精品蜜桃在线观看| 一个人看的www免费观看视频| 欧美精品国产亚洲| 麻豆乱淫一区二区| 久久久久国产网址| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美一区二区三区国产| 特级一级黄色大片| 亚洲国产精品专区欧美| 国产淫片久久久久久久久| 久久久欧美国产精品| 哪个播放器可以免费观看大片| 97热精品久久久久久| 免费看不卡的av| 日韩伦理黄色片| 亚洲精品日本国产第一区| 亚洲无线观看免费| 九九久久精品国产亚洲av麻豆| 亚洲欧美一区二区三区国产| 久久久久久久久久成人| 国产精品av视频在线免费观看| 插阴视频在线观看视频| 精品久久久久久久人妻蜜臀av| 婷婷色综合大香蕉| 欧美成人一区二区免费高清观看| 天堂网av新在线| 亚洲自拍偷在线| 午夜福利在线观看吧| a级一级毛片免费在线观看| 中文欧美无线码| 国产精品久久视频播放| 国产淫语在线视频| 少妇被粗大猛烈的视频| videossex国产| 99久久精品一区二区三区| 亚洲精品成人久久久久久| 亚洲最大成人手机在线| 国产成人一区二区在线| 国产成人91sexporn| 日韩 亚洲 欧美在线| 免费看光身美女| 国产男人的电影天堂91| 汤姆久久久久久久影院中文字幕 | 久久精品久久久久久久性| 九九在线视频观看精品| 中国国产av一级| 内地一区二区视频在线| 亚洲电影在线观看av| 一级二级三级毛片免费看| 亚洲国产精品国产精品| 久久久亚洲精品成人影院| 91久久精品国产一区二区成人| 秋霞伦理黄片| 亚洲欧洲国产日韩| 成人av在线播放网站| 国产毛片a区久久久久| 在线观看一区二区三区| 亚洲精品久久午夜乱码| 校园人妻丝袜中文字幕| 在线观看免费高清a一片| 久久久久久国产a免费观看| 国产一区二区亚洲精品在线观看| 精品一区在线观看国产| 少妇的逼好多水| 一区二区三区四区激情视频| 国产乱人偷精品视频| 三级经典国产精品| 综合色av麻豆| 亚洲不卡免费看| 日日啪夜夜撸| 午夜福利视频1000在线观看| 亚洲,欧美,日韩| 97在线视频观看| 国产乱人偷精品视频| 久久久午夜欧美精品| 亚洲精品影视一区二区三区av| 久久久精品免费免费高清| 搡女人真爽免费视频火全软件| 狂野欧美白嫩少妇大欣赏| 亚洲经典国产精华液单| 国产亚洲最大av| 十八禁国产超污无遮挡网站| 亚洲精品,欧美精品| 午夜福利视频1000在线观看| 熟妇人妻久久中文字幕3abv| xxx大片免费视频| 国产精品伦人一区二区| 日本-黄色视频高清免费观看| 日本av手机在线免费观看| 夜夜爽夜夜爽视频| 最新中文字幕久久久久| 秋霞在线观看毛片| 国产熟女欧美一区二区| 欧美性感艳星| 在线观看人妻少妇| 国产 一区精品| 亚洲四区av| 成年版毛片免费区| 一级毛片电影观看| 国产乱人偷精品视频| 美女大奶头视频| 亚洲人与动物交配视频| 国产黄片美女视频| 久久久久久国产a免费观看| 精品99又大又爽又粗少妇毛片| 男女边摸边吃奶| 一区二区三区乱码不卡18| 中国美白少妇内射xxxbb| 亚洲精品乱久久久久久| 久热久热在线精品观看| 日韩欧美 国产精品| 欧美zozozo另类| 国产成人午夜福利电影在线观看| 久久久久久久久大av| 18禁动态无遮挡网站| 亚洲精品aⅴ在线观看| 一区二区三区免费毛片| 色5月婷婷丁香| 亚洲欧洲日产国产| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩在线观看h| av女优亚洲男人天堂| 三级男女做爰猛烈吃奶摸视频| 观看美女的网站| 久久久欧美国产精品| 爱豆传媒免费全集在线观看| 特大巨黑吊av在线直播| 边亲边吃奶的免费视频| 一区二区三区免费毛片| 久久久久久久久久成人| 亚洲欧洲日产国产| 中文在线观看免费www的网站| 26uuu在线亚洲综合色| 久久精品国产亚洲网站| 永久网站在线| 嫩草影院入口| 特级一级黄色大片| 男女边吃奶边做爰视频| 老女人水多毛片| 少妇被粗大猛烈的视频| 国产乱来视频区| 日韩欧美 国产精品| 亚洲精品国产av成人精品| 亚洲色图av天堂| 国产综合精华液| 赤兔流量卡办理| 国产欧美另类精品又又久久亚洲欧美| 日本免费a在线| 精品人妻熟女av久视频| 波野结衣二区三区在线| 视频中文字幕在线观看| 久久久久久久久中文| 2018国产大陆天天弄谢| 国产色爽女视频免费观看| 日韩伦理黄色片| 亚洲一区高清亚洲精品| 中文天堂在线官网| 最后的刺客免费高清国语| 免费大片黄手机在线观看| 深爱激情五月婷婷| 国产探花在线观看一区二区| 国产 一区精品| 人妻少妇偷人精品九色| 国产精品久久久久久久电影| 亚洲精品影视一区二区三区av| 欧美人与善性xxx| 国国产精品蜜臀av免费| 黄色欧美视频在线观看| 国产 亚洲一区二区三区 | 狂野欧美白嫩少妇大欣赏| 免费观看在线日韩| 日韩欧美精品免费久久| 成人亚洲精品一区在线观看 | 极品少妇高潮喷水抽搐| 久久精品久久久久久噜噜老黄| 亚洲av综合色区一区| 1024视频免费在线观看| 在线 av 中文字幕| 精品国产乱码久久久久久男人| 亚洲成人手机| 国产一区二区激情短视频 | 在线亚洲精品国产二区图片欧美| 日本免费在线观看一区| 免费观看性生交大片5| 亚洲在久久综合| 国产视频首页在线观看| 日韩欧美一区视频在线观看| 成人国语在线视频| 国产精品久久久久久久久免| 青草久久国产| 国产黄色视频一区二区在线观看| 亚洲欧美一区二区三区黑人 | av有码第一页| 中文字幕av电影在线播放| 中文字幕人妻丝袜一区二区 | 国产精品三级大全| 日韩电影二区| 国产男人的电影天堂91| 十八禁高潮呻吟视频| 69精品国产乱码久久久| 黄网站色视频无遮挡免费观看| 欧美+日韩+精品| 国产精品人妻久久久影院| 日日啪夜夜爽| 亚洲精品视频女| 高清不卡的av网站| 欧美日本中文国产一区发布| 国产日韩一区二区三区精品不卡| 精品人妻在线不人妻| 2018国产大陆天天弄谢| 日本欧美视频一区| 色吧在线观看| 97在线人人人人妻| 久久久精品国产亚洲av高清涩受| 免费黄网站久久成人精品| 哪个播放器可以免费观看大片| 丝袜在线中文字幕| 两性夫妻黄色片| 一区二区三区激情视频| 久久99蜜桃精品久久| 在线观看www视频免费| 只有这里有精品99| 性色av一级| 在线观看免费日韩欧美大片| 欧美成人午夜精品| 国产黄频视频在线观看| 性色av一级| 亚洲美女黄色视频免费看| 丰满饥渴人妻一区二区三| freevideosex欧美| 侵犯人妻中文字幕一二三四区| 婷婷色麻豆天堂久久| 国产精品免费视频内射| 婷婷色麻豆天堂久久| 亚洲精品成人av观看孕妇| 2022亚洲国产成人精品| 国产亚洲一区二区精品| 国产黄频视频在线观看| 亚洲五月色婷婷综合| 国产精品国产三级专区第一集| 狂野欧美激情性bbbbbb| 丝袜人妻中文字幕| 国产精品一二三区在线看| 日本色播在线视频| 久久久久久人人人人人| 超色免费av| 黄网站色视频无遮挡免费观看| 国产黄色免费在线视频| 建设人人有责人人尽责人人享有的| 免费人妻精品一区二区三区视频| 亚洲成人一二三区av| 大香蕉久久成人网| 丝袜脚勾引网站| 久久免费观看电影| 韩国高清视频一区二区三区| 国产黄色免费在线视频| 黄网站色视频无遮挡免费观看| 免费人妻精品一区二区三区视频| 妹子高潮喷水视频| 看非洲黑人一级黄片| 熟妇人妻不卡中文字幕| av免费在线看不卡| 国产xxxxx性猛交| av一本久久久久| 亚洲男人天堂网一区| xxxhd国产人妻xxx| 性高湖久久久久久久久免费观看| 少妇猛男粗大的猛烈进出视频| 久久久国产精品麻豆| 高清欧美精品videossex| 肉色欧美久久久久久久蜜桃| 午夜91福利影院| 寂寞人妻少妇视频99o| 国产黄色视频一区二区在线观看| 视频在线观看一区二区三区| 蜜桃国产av成人99| 日韩不卡一区二区三区视频在线| 久久精品国产亚洲av天美| 色吧在线观看| 在线观看免费视频网站a站| 黑人欧美特级aaaaaa片| 看十八女毛片水多多多| 我要看黄色一级片免费的| 交换朋友夫妻互换小说| 三级国产精品片| 欧美人与性动交α欧美精品济南到 | 亚洲欧美清纯卡通| 电影成人av| 国产精品久久久久久精品古装| 最新中文字幕久久久久| 精品亚洲乱码少妇综合久久| 国产男人的电影天堂91| 国产不卡av网站在线观看| xxxhd国产人妻xxx| 国产男女超爽视频在线观看| 亚洲四区av| 国产黄色免费在线视频|