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

    基于β-屈服函數(shù)的炭質(zhì)泥頁巖統(tǒng)一硬化模型

    2023-02-24 02:16:42崔凱胡斌崔阿能馬利遙劉楊祝鑫
    關(guān)鍵詞:炭質(zhì)劍橋屈服

    崔凱,胡斌,崔阿能,馬利遙,劉楊,祝鑫

    (武漢科技大學(xué)資源與環(huán)境工程學(xué)院;冶金礦產(chǎn)資源高效利用與造塊湖北省重點實驗室,武漢 430081)

    在中國中南、西南地區(qū),如云南、四川、貴州、廣西等省存在大范圍的含緩傾軟弱夾層(巖層和夾層傾角小于25°)的二疊系石灰?guī)r地層[1-3],該地層擁有豐富的優(yōu)質(zhì)石灰石礦產(chǎn)資源,是中國大量基礎(chǔ)設(shè)施建設(shè)的重要建材來源。如圖1所示,這些含緩傾軟弱夾層通常由炭質(zhì)泥頁巖組成,具有抗剪強度小且容易在外部因素(人為擾動或環(huán)境)作用下發(fā)生劣化等特點,因此,含緩傾軟弱夾層被視作礦山高邊坡的薄弱帶[4-5]。深入探究炭質(zhì)泥頁巖的力學(xué)特性并建立相應(yīng)的本構(gòu)關(guān)系,對于含緩傾軟弱夾層的礦山高邊坡的穩(wěn)定性分析具有重要意義。

    圖1 礦山高邊坡中含炭質(zhì)泥頁巖的軟弱夾層Fig.1 Weak intercalation of carbonaceous shale in the high slope of the mine

    在復(fù)雜加載條件下,炭質(zhì)泥頁巖類似于超固結(jié)黏土,會發(fā)生彈塑性破壞,其應(yīng)力—應(yīng)變曲線具有應(yīng)變軟化的特點,并且伴隨著剪脹現(xiàn)象的發(fā)生[6-7]。為此,有學(xué)者[8]建議將炭質(zhì)泥頁巖等強度較小的軟巖視作超固結(jié)比很大的黏土,并采用臨界狀態(tài)土力學(xué)的理論框架來描述其常規(guī)力學(xué)特性。宋麗等[9]建立了軟巖的三維統(tǒng)一彈黏塑性本構(gòu)模型,該模型對于日本石川縣能登半島端部的硅藻質(zhì)泥巖在不同圍壓和加載速率等條件下的不排水三軸壓縮試驗結(jié)果具有較好的描述效果。熊勇林等[10]同時考慮了圍壓和溫度對軟巖殘余強度的影響,建立了軟巖的熱彈黏塑性模型。炭質(zhì)泥頁巖另一個顯著的力學(xué)特性是具有流變性,在保持荷載強度不變的情況下,增加加載時間將導(dǎo)致額外的體積變形[11]。范慶忠等[12]提出了軟巖的非線性蠕變模型,齊明山[13]則系統(tǒng)地研究了大變形軟巖流變性態(tài)及其在隧道工程結(jié)構(gòu)中的應(yīng)用。在所有臨界狀態(tài)本構(gòu)關(guān)系中,修正劍橋模型[14]在巖土工程中應(yīng)用最廣泛,然而其只適用于正常固結(jié)黏土力學(xué)特性的描述,在描述炭質(zhì)泥頁巖等超固結(jié)黏土?xí)r具有明顯的局限性。例如,修正劍橋模型計算得到的超固結(jié)黏土的不排水抗剪強度明顯偏大,且應(yīng)力—應(yīng)變曲線為一條不光滑的曲線。另一方面,修正劍橋模型的剪脹方程在描述某些特定類型土的試驗結(jié)果時存在一定的偏差,同時,無法考慮超固結(jié)比和溫度等因素對土的剪脹特性的影響[10,15]。為了克服上述不足,Yamakawa等[16]提出了下加載面的力學(xué)概念,在加載的任意過程中,假定土體當(dāng)前應(yīng)力點始終位于下加載面上,發(fā)生彈塑性體積變形,物理意義明確,能連續(xù)平滑地刻畫超固結(jié)土的彈塑性應(yīng)力—應(yīng)變關(guān)系。Yao等[17-19]定義了統(tǒng)一硬化參量,建立了適用于超固結(jié)重塑土的統(tǒng)一硬化(Unified hardening,UH)模型并取得了成功。統(tǒng)一硬化模型能夠較好地描述超固結(jié)黏土的應(yīng)變軟化和剪脹特性;同時,在加載過程中,當(dāng)土體由超固結(jié)狀態(tài)轉(zhuǎn)變?yōu)檎9探Y(jié)狀態(tài)時,統(tǒng)一硬化模型能夠計算得到光滑連續(xù)的應(yīng)力—應(yīng)變曲線。

    筆者將炭質(zhì)泥頁巖視作超固結(jié)黏土,并且基于β-屈服函數(shù)建立了能夠描述其應(yīng)變軟化和剪脹等復(fù)雜力學(xué)行為的修正統(tǒng)一硬化模型。通過將模型計算結(jié)果與炭質(zhì)泥頁巖和太古石的排水三軸壓縮試驗結(jié)果進(jìn)行對比分析,驗證了模型的合理性。

    1 修正統(tǒng)一硬化模型

    假定軟巖的體積變形以壓為正,其中,有效應(yīng)力p、剪應(yīng)力q和剪應(yīng)力比η分別定義為

    式中:σ1、σ2和σ3為主應(yīng)力。

    1.1 β-屈服函數(shù)

    在三軸壓縮條件下,修正劍橋模型的塑性功增量dWin可以表示為

    式中:dεpv和dεps分別為塑性體積應(yīng)變增量和塑性剪切應(yīng)變增量;M為臨界狀態(tài)應(yīng)力比。

    根據(jù)上述能量耗散方程,可以分別得到修正劍橋模型的剪脹方程和塑性勢。

    式中:d為剪脹因子,可以定義為d=dεpv/dεps;pc為前期固結(jié)應(yīng)力。由于修正劍橋模型采用相關(guān)聯(lián)的流動法則,式(6)同樣可以用于描述土的屈服面。

    在描述超固結(jié)黏土的剪脹特性時,修正劍橋模型的剪脹方程具有一定的局限性,即低應(yīng)力條件下該剪脹方程預(yù)測的剪脹因子d明顯小于土的實測值。此外,由該剪脹方程得到的屈服函數(shù)假定土的物理屈服面具有橢圓形幾何形狀,與部分土的試驗結(jié)果不符。

    為了克服上述不足,首先將修正劍橋模型的剪脹方程改寫為

    式中:β為材料參數(shù)。令β=2,式(7)將直接退化為修正劍橋模型的剪脹方程。

    如圖2所示,與修正劍橋模型的剪脹方程相比,該剪脹方程能夠較好地描述太古石的試驗結(jié)果[7],尤其是在低應(yīng)力水平條件下。

    圖2 炭質(zhì)泥頁巖剪脹曲線Fig.2 Dilatancy curve of carbonaceous mud shale

    通??梢愿鶕?jù)剪脹方程得到塑性勢并確定土的塑性應(yīng)變率,在p-q應(yīng)力空間,塑性勢具有式(8)所示積分表達(dá)式[14]。

    將式(7)代入式(8)并假定相關(guān)聯(lián)流動法則,可得到式(9)所示的β-屈服函數(shù)。

    為了便于探究材料參數(shù)β對屈服面幾何形狀的影響,可以將β-屈服函數(shù)改寫為

    根據(jù)式(10),當(dāng)β=2時,可以得到修正劍橋模型的屈服函數(shù)。

    材料參數(shù)β對屈服面幾何形狀的影響結(jié)果圖3所示。由圖3可以看出,增大參數(shù)β將導(dǎo)致計算得到的屈服面出現(xiàn)明顯收縮,其彈性區(qū)間會隨之減小。

    圖3 參數(shù)β對屈服面幾何形狀的影響Fig.3 The influence of parameter β on the geometry of the yield surface

    β-屈服函數(shù)關(guān)于有效平均主應(yīng)力p和剪應(yīng)力q的偏導(dǎo)數(shù)分別為

    將式(11)、式(12)代入相關(guān)聯(lián)的流動法則中,可以計算得到炭質(zhì)泥頁巖在剪切過程中的塑性應(yīng)變率。

    1.2 相似因子

    炭質(zhì)泥頁巖具有相對較好的連續(xù)性,其變形特征與典型的彈塑性材料類似,即在圍壓很小的情況下表現(xiàn)為脆性破壞,而隨著圍壓的增大,將逐漸轉(zhuǎn)變?yōu)檠有云茐腫7]。

    在加載的初始階段,炭質(zhì)泥頁巖的應(yīng)力—應(yīng)變曲線具有線性增長的特點,此時試樣僅發(fā)生彈性變形。進(jìn)一步加載將導(dǎo)致試樣內(nèi)部的膠結(jié)作用逐漸破壞并形成微裂紋,這些微裂紋的發(fā)展和貫通將導(dǎo)致試樣發(fā)生屈服,其應(yīng)力—應(yīng)變曲線的非線性變化特點也將逐漸增強,表現(xiàn)為應(yīng)變硬化。在達(dá)到峰值強度后,試樣的剪應(yīng)力會逐漸減小,其應(yīng)力—應(yīng)變曲線表現(xiàn)為應(yīng)變軟化。

    可以看出,在復(fù)雜應(yīng)力條件下,炭質(zhì)泥頁巖的應(yīng)力—應(yīng)變曲線變化特點與超固結(jié)黏土十分類似,均表現(xiàn)為應(yīng)變硬化和應(yīng)變軟化。為此,將炭質(zhì)泥頁巖視作超固結(jié)黏土,并且定義其超固結(jié)比OCR為抗剪強度fc與圍壓pic之間的比值,即OCR=fc/pic。

    參考姚仰平等[19]提出的適用于超固結(jié)黏土的統(tǒng)一硬化模型,假定在任意加載過程中炭質(zhì)泥頁巖當(dāng)前的應(yīng)力狀態(tài)點A(p,q)始終位于下加載面上,其參考應(yīng)力狀態(tài)點B(pˉ,qˉ)則位于參考屈服面上,如圖4所示。

    圖4 下加載面和參考屈服面力學(xué)示意圖Fig.4 Schematic diagram of the mechanics of lower loading surface and reference yield surface

    在p-q平面,下加載面和參考屈服面之間具有幾何相似性,兩者之間的相對位置關(guān)系可以采用相似因子R進(jìn)行描述,R的定義為

    根據(jù)式(13),當(dāng)R=1時,下加載面會與參考屈服面重合,此時炭質(zhì)泥頁巖處于正常固結(jié)狀態(tài)。R與OCR具有類似的物理意義,均可以用于描述炭質(zhì)泥頁巖的超固結(jié)性。對于三軸壓縮試驗,剪切伊始,不同圍壓條件下炭質(zhì)泥頁巖的初始相似因子R0為

    采用本文提出的β-屈服函數(shù),參考應(yīng)力狀態(tài)點B()所在的參考屈服面可以表示為

    當(dāng)應(yīng)力狀態(tài)點位于參考屈服面上時,炭質(zhì)泥頁巖將處于正常固結(jié)狀態(tài),因此,類似于修正劍橋模型,可以采用體積硬化法則來描述硬化參量p0-pc的演化規(guī)律。

    式中:cp=(λ-κ)/(1+e0),其中,λ、κ分別為炭質(zhì)泥頁巖的壓縮模量和回彈模量;e0為炭質(zhì)泥頁巖的初始孔隙比。

    將式(16)代入式(15)中,可以得到

    整理式(17),可以得到任意加載時刻有效平均主應(yīng)力的表示式

    根據(jù)相似因子R的定義,可以進(jìn)一步得到R的表達(dá)式

    可以看出,R與應(yīng)力水平、應(yīng)力路徑和累積塑性體積應(yīng)變等狀態(tài)變量有關(guān),式(19)能夠綜合反映上述影響因素對下加載面大小變化規(guī)律的影響。

    1.3 統(tǒng)一硬化參量

    由于當(dāng)前應(yīng)力狀態(tài)點所在的下加載面幾何相似于參考屈服面,因此,其屈服函數(shù)可以表示為

    式中,統(tǒng)一硬化參量H的定義式為

    根據(jù)Yao等[18]的研究,Ω可以表示為

    式中:Mf為潛在峰值強度剪應(yīng)力比。

    在統(tǒng)一硬化模型[17-19]中,Mf被認(rèn)為是關(guān)于相似因子R和臨界狀態(tài)剪應(yīng)力比M的函數(shù),滿足關(guān)系式

    式中,材料參數(shù)k可表示為

    統(tǒng)一硬化參量Ω和剪應(yīng)力比η之間的關(guān)系如圖5所示,可以看出,η=0時,Ω=1;η=M時,Ω=0;η=Mf時,Ω=-∞。

    圖5 剪應(yīng)力比η與統(tǒng)一硬化參數(shù)Ω的關(guān)系曲線Fig.5 Relationship curve between shear stress ratio η and unified hardening parameter Ω

    在模型數(shù)值實現(xiàn)過程中,需要確定炭質(zhì)泥頁巖的體變模量K和剪切模量G,參考修正劍橋模型,可以采用式(25)計算。

    式中:ν為泊松比。

    2 參數(shù)試驗

    在統(tǒng)一硬化模型的基礎(chǔ)上,采用了β-屈服函數(shù)來描述炭質(zhì)泥頁巖的剪脹特性,由此建立的模型能夠合理地描述應(yīng)變硬化和軟化特性,同時考慮了圍壓對炭質(zhì)泥頁巖力學(xué)行為的影響。

    模型包含的材料參數(shù)為M、e0、ν、λ、κ和β。其中,e0、λ、κ的具體取值可以通過炭質(zhì)泥頁巖的等向壓縮曲線得到,臨界狀態(tài)應(yīng)力比M則可以由炭質(zhì)泥頁巖在p—q平面的臨界狀態(tài)線加以確定。不同類型巖土材料的泊松比ν的取值范圍通常為0.2~0.3,可以由特定圍壓下材料的初始應(yīng)力—應(yīng)變曲線得到。

    2.1 超固結(jié)比對模型計算結(jié)果的影響

    為了初步驗證所建立的適用于炭質(zhì)泥頁巖的統(tǒng)一硬化模型的適用性,并且探究初始超固結(jié)比對炭質(zhì)泥頁巖力學(xué)特性的影響,開展了圍壓為200 kPa的排水三軸壓縮數(shù)值試驗,試樣的初始超固結(jié)比分別為1、2、4和8,材料參數(shù)如表1所示。

    表1 數(shù)值試樣的材料參數(shù)Table 1 Material parameters of carbonaceous mud shale numerical samples

    模型計算結(jié)果如圖6所示。由圖6可以看出,隨著初始超固結(jié)比逐漸增大,試樣的應(yīng)力—應(yīng)變曲線表現(xiàn)出明顯的應(yīng)變軟化特性,同時,其體積變形將逐漸由剪縮變?yōu)榧裘洝?/p>

    圖6 初始超固結(jié)比對炭質(zhì)泥頁巖排水三軸壓縮試驗結(jié)果的影響Fig.6 Influence of initial over-consolidation ratio on the results of drainage triaxial compression test of carbonaceous shale

    2.2 材料參數(shù)β對模型計算結(jié)果的影響

    為了更好地描述炭質(zhì)泥頁巖的剪脹特性,對修正劍橋模型進(jìn)行了適當(dāng)?shù)男拚?,并由此得到?屈服函數(shù)。在上述參數(shù)試驗的基礎(chǔ)上,分別令β為1.5、2、2.5和3,計算得到超固結(jié)比為8的試樣的三軸壓縮排水試驗結(jié)果,其中,加載圍壓為200 kPa。

    根據(jù)圖7,增大β的取值會導(dǎo)致試樣的峰值強度逐漸增大,其對應(yīng)的軸向影響則會相應(yīng)減小。此外,隨著β的逐漸增大,試樣的最終體積剪脹量也將隨之增大。

    圖7 參數(shù)β對炭質(zhì)泥頁巖排水三軸壓縮試驗結(jié)果的影響Fig.7 Effect of parameter β on the results of drainage triaxial compression test of carbonaceous mud shale

    3 模型驗證

    3.1 炭質(zhì)泥頁巖

    為了進(jìn)一步驗證模型的適用性,采用該模型對炭質(zhì)泥頁巖[20]的三軸壓縮試驗結(jié)果進(jìn)行計算,材料參數(shù)如表2所示。假定該炭質(zhì)泥頁巖的前期固結(jié)應(yīng)力為50 MPa,當(dāng)圍壓分別為100、200、300、400 kPa時,其初始超固結(jié)比分別為500、250、165、125。

    表2 炭質(zhì)泥頁巖模型參數(shù)Table 2 Model parameters of carbonaceous mud shale

    在不同圍壓下,試樣的應(yīng)力—應(yīng)變曲線試驗結(jié)果和模型計算結(jié)果如圖8(a)所示。對比結(jié)果表明,該模型能夠較好地描述炭質(zhì)泥頁巖的強度和變形特點,尤其是達(dá)到峰值強度后的應(yīng)力—應(yīng)變曲線,此時炭質(zhì)泥頁巖表現(xiàn)出明顯的應(yīng)變軟化特性。不過,在低圍壓條件下,模型計算得到的峰值強度明顯小于試驗結(jié)果,峰值強度對應(yīng)的軸向應(yīng)變也明顯偏大,可以調(diào)整部分材料參數(shù)的取值來進(jìn)一步提高模型的計算精度。

    圖8(b)給出了模型計算得到的炭質(zhì)泥頁巖的剪脹曲線,可以看出,在加載的初始階段,試樣會發(fā)生剪縮變形,而進(jìn)一步加載將導(dǎo)致試樣發(fā)生剪脹破壞。

    圖8 不同圍壓條件下炭質(zhì)泥頁巖排水三軸壓縮試驗結(jié)果與模型計算結(jié)果對比Fig.8 Comparison of the results of triaxial compression test and model calculation of carbonaceous shale under different confining pressures

    3.2 太古石

    基于β-屈服函數(shù)的統(tǒng)一硬化模型同樣能夠較好地描述其他不同類型軟巖的強度和變形特性。對Adachi等[21]的太古石三軸壓縮排水試驗數(shù)據(jù)進(jìn)行計算,材料參數(shù)如表3所示。

    表3 太古石模型參數(shù)Table 3 Model parameters of soft rock

    圖9給出了太古石三軸壓縮排水試驗結(jié)果與計算結(jié)果之間的對比??梢钥闯觯诓煌瑖鷫合?,該模型能較好地描述太古石三軸壓縮排水試驗應(yīng)力—應(yīng)變曲線和體積變形規(guī)律。

    圖9 太古石三軸排水壓縮試驗結(jié)果與計算結(jié)果的對比Fig.9 Comparison between results of drained triaxial tests and calculation on soft rock

    4 結(jié)論

    建立了適用于炭質(zhì)泥頁巖的β-屈服函數(shù)統(tǒng)一硬化模型,該模型能較好地描述炭質(zhì)泥頁巖的應(yīng)變軟化和剪脹特性。得到以下主要結(jié)論:

    1)炭質(zhì)泥頁巖可以被視作超固結(jié)黏土,在任意加載時刻,其當(dāng)前的應(yīng)力狀態(tài)點始終位于下加載面上,并且發(fā)生彈塑性破壞。

    2)該模型能準(zhǔn)確地描述炭質(zhì)泥頁巖的剪脹特性,其采用的β-屈服函數(shù)能考慮屈服面的幾何形狀對模型計算結(jié)果的影響,從而顯著提高模型的計算精度。

    3)模型包含的材料參數(shù)具有明確的物理意義,能夠通過常規(guī)的室內(nèi)試驗進(jìn)行標(biāo)定,便于實際工程應(yīng)用。模型計算結(jié)果與試驗結(jié)果進(jìn)行對比分析表明,模型取得了良好的擬合效果,說明該模型能準(zhǔn)確、合理地描述炭質(zhì)泥頁巖的常規(guī)力學(xué)行為。

    該模型的主要不足是無法考慮在復(fù)雜應(yīng)力條件下炭質(zhì)泥頁巖內(nèi)部膠結(jié)作用破壞對其強度和變形特性的影響。在后續(xù)研究中,需要同時考慮結(jié)構(gòu)性和超固結(jié)性對炭質(zhì)泥頁巖常規(guī)力學(xué)行為的影響。

    猜你喜歡
    炭質(zhì)劍橋屈服
    牙被拔光也不屈服的史良大律師秘書
    紅巖春秋(2022年1期)2022-04-12 00:37:34
    基于敏感性分析的炭質(zhì)板巖引水隧洞支護結(jié)構(gòu)優(yōu)化研究
    高速公路炭質(zhì)頁巖高邊坡穩(wěn)定性分析及施工方案
    基于連續(xù)-非連續(xù)單元方法的炭質(zhì)板巖隧道圍巖穩(wěn)定分析
    Cambridge Audio(劍橋)CXA80解碼/放大器一體機
    The Classic Lines of A Love so Beautiful
    劍橋能源周的啟示
    能源(2018年4期)2018-05-19 01:53:58
    劍橋是最好的起飛平臺
    商周刊(2017年7期)2017-08-22 03:36:20
    勇敢
    干燥和飽水炭質(zhì)板巖流變力學(xué)特性與模型研究
    日韩亚洲欧美综合| 国产 精品1| 赤兔流量卡办理| 人人妻人人澡人人看| 狂野欧美白嫩少妇大欣赏| 欧美xxⅹ黑人| 国产精品秋霞免费鲁丝片| 亚洲欧美日韩另类电影网站| 一区在线观看完整版| 在线观看免费视频网站a站| 国产精品久久久久久av不卡| 观看av在线不卡| 97超碰精品成人国产| 水蜜桃什么品种好| 国产精品人妻久久久影院| 国产熟女欧美一区二区| 涩涩av久久男人的天堂| 日日摸夜夜添夜夜爱| 看非洲黑人一级黄片| 丝袜喷水一区| 亚洲va在线va天堂va国产| 日韩不卡一区二区三区视频在线| 亚洲天堂av无毛| 国产精品国产av在线观看| 丝袜脚勾引网站| 最近最新中文字幕免费大全7| 国产亚洲av片在线观看秒播厂| 在线观看国产h片| 汤姆久久久久久久影院中文字幕| 热re99久久精品国产66热6| 亚洲怡红院男人天堂| 成人亚洲欧美一区二区av| 欧美 亚洲 国产 日韩一| 国产精品三级大全| 欧美人与善性xxx| 高清毛片免费看| 王馨瑶露胸无遮挡在线观看| 精品久久国产蜜桃| 精品国产一区二区久久| 汤姆久久久久久久影院中文字幕| 国产欧美亚洲国产| 国产伦精品一区二区三区四那| 久久精品久久久久久久性| 午夜老司机福利剧场| 国产男女内射视频| 免费播放大片免费观看视频在线观看| 亚洲电影在线观看av| 伦精品一区二区三区| 一级爰片在线观看| 纯流量卡能插随身wifi吗| 国产精品一区www在线观看| 观看美女的网站| 免费不卡的大黄色大毛片视频在线观看| 久久99一区二区三区| av有码第一页| 赤兔流量卡办理| 国产亚洲av片在线观看秒播厂| 亚洲精品成人av观看孕妇| 亚洲va在线va天堂va国产| 国产黄片视频在线免费观看| 99久久中文字幕三级久久日本| 三级国产精品欧美在线观看| 国产日韩欧美亚洲二区| 在线观看www视频免费| 在现免费观看毛片| 啦啦啦视频在线资源免费观看| 国产国拍精品亚洲av在线观看| 免费人妻精品一区二区三区视频| 美女脱内裤让男人舔精品视频| 日韩欧美一区视频在线观看 | 成人美女网站在线观看视频| 建设人人有责人人尽责人人享有的| 精品国产国语对白av| 欧美日韩视频高清一区二区三区二| 一本—道久久a久久精品蜜桃钙片| 最近中文字幕2019免费版| 内射极品少妇av片p| 日韩一区二区视频免费看| 十八禁高潮呻吟视频 | 国产高清不卡午夜福利| 国产伦在线观看视频一区| 亚洲不卡免费看| 国产黄色视频一区二区在线观看| 国产av国产精品国产| 我要看日韩黄色一级片| 国产男女内射视频| 亚洲精品色激情综合| 波野结衣二区三区在线| 2021少妇久久久久久久久久久| 亚洲精品久久午夜乱码| 亚洲国产精品一区三区| 少妇的逼好多水| 久久久久久伊人网av| 如日韩欧美国产精品一区二区三区 | 亚洲精品国产色婷婷电影| 国产男人的电影天堂91| 成人国产av品久久久| 麻豆成人av视频| 久久毛片免费看一区二区三区| 亚洲国产精品国产精品| 青春草国产在线视频| 亚洲四区av| 国模一区二区三区四区视频| 国产高清有码在线观看视频| 久久午夜综合久久蜜桃| 伦理电影免费视频| 男女国产视频网站| 国产精品免费大片| 777米奇影视久久| 久久久久视频综合| 日日撸夜夜添| 久久久久久久久大av| 视频区图区小说| 国产欧美日韩综合在线一区二区 | 久久精品久久久久久噜噜老黄| 久久热精品热| 我的老师免费观看完整版| 中文天堂在线官网| 免费播放大片免费观看视频在线观看| 国产亚洲5aaaaa淫片| 只有这里有精品99| 免费人妻精品一区二区三区视频| 精品熟女少妇av免费看| 美女脱内裤让男人舔精品视频| 看免费成人av毛片| 亚洲性久久影院| 亚洲,欧美,日韩| 高清视频免费观看一区二区| 日韩在线高清观看一区二区三区| 国产亚洲av片在线观看秒播厂| 亚洲美女搞黄在线观看| 麻豆乱淫一区二区| 少妇人妻精品综合一区二区| 中文乱码字字幕精品一区二区三区| 天美传媒精品一区二区| 我要看日韩黄色一级片| 噜噜噜噜噜久久久久久91| 精品一区在线观看国产| 久久久精品94久久精品| 亚洲国产欧美在线一区| 亚洲精品视频女| 欧美老熟妇乱子伦牲交| a级毛色黄片| 9色porny在线观看| 国产极品天堂在线| 80岁老熟妇乱子伦牲交| 亚洲欧美一区二区三区国产| 欧美国产精品一级二级三级 | 日韩中文字幕视频在线看片| 五月开心婷婷网| 亚洲av不卡在线观看| 中文字幕av电影在线播放| 亚洲欧美精品专区久久| 九九爱精品视频在线观看| 王馨瑶露胸无遮挡在线观看| 日日撸夜夜添| 日本黄大片高清| 伊人久久国产一区二区| 亚洲婷婷狠狠爱综合网| 桃花免费在线播放| 五月玫瑰六月丁香| 黑人高潮一二区| 少妇的逼好多水| 丰满乱子伦码专区| 色哟哟·www| av一本久久久久| 久久午夜综合久久蜜桃| 综合色丁香网| av福利片在线观看| 十分钟在线观看高清视频www | 国产一区二区在线观看av| 日韩制服骚丝袜av| 亚洲精品乱久久久久久| 精品一区二区免费观看| 极品人妻少妇av视频| 18禁裸乳无遮挡动漫免费视频| 3wmmmm亚洲av在线观看| 在线观看免费日韩欧美大片 | 欧美日韩精品成人综合77777| 日本免费在线观看一区| 交换朋友夫妻互换小说| 天天操日日干夜夜撸| 日日啪夜夜撸| 亚洲真实伦在线观看| 五月伊人婷婷丁香| 亚洲av成人精品一二三区| 久久精品久久久久久噜噜老黄| 欧美bdsm另类| 99热国产这里只有精品6| 哪个播放器可以免费观看大片| 免费人妻精品一区二区三区视频| 国产黄色免费在线视频| 深夜a级毛片| 欧美bdsm另类| 日韩三级伦理在线观看| 丰满乱子伦码专区| 观看美女的网站| videossex国产| 国内精品宾馆在线| 国产精品99久久99久久久不卡 | 国产女主播在线喷水免费视频网站| 精品熟女少妇av免费看| 日韩一本色道免费dvd| 男女国产视频网站| 日韩不卡一区二区三区视频在线| 亚洲人与动物交配视频| 99九九线精品视频在线观看视频| 亚洲伊人久久精品综合| 26uuu在线亚洲综合色| 看非洲黑人一级黄片| 免费播放大片免费观看视频在线观看| 欧美日韩一区二区视频在线观看视频在线| 亚洲中文av在线| 国产成人精品福利久久| 久久国产精品大桥未久av | 69精品国产乱码久久久| 免费高清在线观看视频在线观看| 寂寞人妻少妇视频99o| 亚洲精品久久久久久婷婷小说| 国产真实伦视频高清在线观看| 午夜av观看不卡| 亚洲综合色惰| 国内揄拍国产精品人妻在线| 日韩中文字幕视频在线看片| 又粗又硬又长又爽又黄的视频| 国产欧美另类精品又又久久亚洲欧美| 99久久精品热视频| 一本色道久久久久久精品综合| 又粗又硬又长又爽又黄的视频| 啦啦啦在线观看免费高清www| 九九久久精品国产亚洲av麻豆| 日韩精品有码人妻一区| 日韩伦理黄色片| 亚洲精品色激情综合| 国产国拍精品亚洲av在线观看| 亚洲精品色激情综合| 日韩中文字幕视频在线看片| 久久99精品国语久久久| 日本wwww免费看| www.色视频.com| 日本vs欧美在线观看视频 | 精品少妇内射三级| 精品国产国语对白av| 精品久久久精品久久久| 少妇 在线观看| 免费看日本二区| 久久精品国产鲁丝片午夜精品| 亚洲人成网站在线播| 老女人水多毛片| 免费人成在线观看视频色| 日日摸夜夜添夜夜添av毛片| 免费av不卡在线播放| 寂寞人妻少妇视频99o| 中文欧美无线码| 大码成人一级视频| 国产综合精华液| 亚洲精品国产av成人精品| 麻豆乱淫一区二区| 伊人久久国产一区二区| 久久精品国产鲁丝片午夜精品| 在线观看免费视频网站a站| 成年人免费黄色播放视频 | 国产 一区精品| 少妇的逼好多水| 王馨瑶露胸无遮挡在线观看| 久久久久网色| 热re99久久精品国产66热6| 观看免费一级毛片| 蜜桃在线观看..| 性高湖久久久久久久久免费观看| 夜夜看夜夜爽夜夜摸| 亚洲精品,欧美精品| 欧美bdsm另类| 国产精品一区二区在线不卡| 在线观看免费高清a一片| 国产淫片久久久久久久久| 国产亚洲av片在线观看秒播厂| 日韩欧美一区视频在线观看 | 国产男人的电影天堂91| 国产黄片视频在线免费观看| 免费观看的影片在线观看| 99热这里只有是精品在线观看| 国产91av在线免费观看| 赤兔流量卡办理| 国产精品免费大片| 久久精品久久精品一区二区三区| 熟女人妻精品中文字幕| 久久精品国产亚洲av天美| 黄色配什么色好看| 久久热精品热| 深夜a级毛片| 欧美精品国产亚洲| 午夜福利影视在线免费观看| 久久精品国产a三级三级三级| 亚洲精品一区蜜桃| 两个人免费观看高清视频 | 精品一区在线观看国产| 欧美日韩国产mv在线观看视频| 热99国产精品久久久久久7| 亚洲av在线观看美女高潮| 久久久久久久亚洲中文字幕| 丝袜喷水一区| 中文在线观看免费www的网站| 国产日韩欧美视频二区| 国产片特级美女逼逼视频| 中文乱码字字幕精品一区二区三区| 18禁在线播放成人免费| 搡老乐熟女国产| 成人影院久久| 男女边吃奶边做爰视频| www.色视频.com| 人人妻人人澡人人看| 国内揄拍国产精品人妻在线| 国产一区二区三区av在线| 免费观看在线日韩| 看免费成人av毛片| 18+在线观看网站| 午夜免费观看性视频| 99热这里只有精品一区| 亚洲人与动物交配视频| 日韩免费高清中文字幕av| 精品少妇久久久久久888优播| 亚洲综合色惰| 久久久精品94久久精品| av国产久精品久网站免费入址| 亚洲欧美日韩另类电影网站| 一级爰片在线观看| .国产精品久久| 日本午夜av视频| 青春草亚洲视频在线观看| av女优亚洲男人天堂| av黄色大香蕉| 国产成人freesex在线| 午夜日本视频在线| 欧美精品国产亚洲| 99久国产av精品国产电影| 日韩免费高清中文字幕av| 在现免费观看毛片| 乱码一卡2卡4卡精品| 精品酒店卫生间| 国产黄频视频在线观看| 伦理电影免费视频| 下体分泌物呈黄色| 国产精品伦人一区二区| 中文资源天堂在线| 日本黄大片高清| 亚洲一级一片aⅴ在线观看| 久久 成人 亚洲| 婷婷色av中文字幕| 午夜视频国产福利| 久久久亚洲精品成人影院| 91久久精品国产一区二区三区| 久久韩国三级中文字幕| 18+在线观看网站| 汤姆久久久久久久影院中文字幕| 18禁裸乳无遮挡动漫免费视频| 美女福利国产在线| 丝袜在线中文字幕| 免费大片黄手机在线观看| 欧美bdsm另类| 最近中文字幕高清免费大全6| 国产日韩欧美视频二区| 国产成人freesex在线| 亚洲欧美精品专区久久| 久久久久人妻精品一区果冻| 国产中年淑女户外野战色| 亚洲精华国产精华液的使用体验| 99视频精品全部免费 在线| 国产国拍精品亚洲av在线观看| 在现免费观看毛片| 国产欧美日韩精品一区二区| av播播在线观看一区| 校园人妻丝袜中文字幕| 久久精品久久久久久久性| 熟女av电影| 国产极品天堂在线| 激情五月婷婷亚洲| 热re99久久精品国产66热6| 岛国毛片在线播放| 菩萨蛮人人尽说江南好唐韦庄| 成人特级av手机在线观看| 亚洲欧美日韩另类电影网站| 国产精品人妻久久久久久| 如日韩欧美国产精品一区二区三区 | 色视频在线一区二区三区| 亚洲av不卡在线观看| 日韩,欧美,国产一区二区三区| 欧美精品一区二区免费开放| 9色porny在线观看| 如何舔出高潮| 伦理电影大哥的女人| 日韩伦理黄色片| 亚洲欧美清纯卡通| 精品一区二区三卡| 亚洲av中文av极速乱| 久久久久久久精品精品| av线在线观看网站| 精品熟女少妇av免费看| 我要看日韩黄色一级片| 亚洲国产精品国产精品| 午夜福利在线观看免费完整高清在| 日韩制服骚丝袜av| 热re99久久精品国产66热6| 我的老师免费观看完整版| 在线观看美女被高潮喷水网站| 亚洲av在线观看美女高潮| 亚洲av国产av综合av卡| 美女福利国产在线| 欧美高清成人免费视频www| 色吧在线观看| 丰满饥渴人妻一区二区三| 边亲边吃奶的免费视频| 国产高清有码在线观看视频| 在线观看免费视频网站a站| 国产极品粉嫩免费观看在线 | 欧美日韩在线观看h| 毛片一级片免费看久久久久| a级一级毛片免费在线观看| 人人妻人人澡人人看| 亚洲欧美中文字幕日韩二区| 色94色欧美一区二区| 亚洲精品国产成人久久av| 女的被弄到高潮叫床怎么办| 午夜福利在线观看免费完整高清在| 国产精品免费大片| 丝瓜视频免费看黄片| 欧美精品一区二区大全| 嫩草影院入口| 成人美女网站在线观看视频| 看非洲黑人一级黄片| 女性被躁到高潮视频| 九草在线视频观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲人与动物交配视频| 人妻制服诱惑在线中文字幕| 日韩三级伦理在线观看| 中文乱码字字幕精品一区二区三区| 伊人久久精品亚洲午夜| 亚洲av.av天堂| 热re99久久精品国产66热6| 亚洲成人av在线免费| 伦理电影免费视频| 高清毛片免费看| 肉色欧美久久久久久久蜜桃| av一本久久久久| 久热久热在线精品观看| 高清黄色对白视频在线免费看 | 观看免费一级毛片| 亚洲怡红院男人天堂| 韩国av在线不卡| 又黄又爽又刺激的免费视频.| videos熟女内射| 久久久亚洲精品成人影院| 黄色毛片三级朝国网站 | 只有这里有精品99| 2021少妇久久久久久久久久久| 国产成人午夜福利电影在线观看| 亚洲熟女精品中文字幕| 亚洲高清免费不卡视频| 赤兔流量卡办理| 久久精品久久久久久噜噜老黄| 午夜老司机福利剧场| av黄色大香蕉| 亚洲一区二区三区欧美精品| 久久 成人 亚洲| 精品一品国产午夜福利视频| 人人妻人人澡人人看| 国产老妇伦熟女老妇高清| 亚洲在久久综合| av免费在线看不卡| 欧美xxxx性猛交bbbb| 少妇高潮的动态图| 99热国产这里只有精品6| 国产免费一区二区三区四区乱码| tube8黄色片| 国产真实伦视频高清在线观看| 一级毛片我不卡| 久久毛片免费看一区二区三区| 欧美日韩在线观看h| 日本欧美国产在线视频| 欧美激情国产日韩精品一区| 日韩制服骚丝袜av| 51国产日韩欧美| 女人久久www免费人成看片| 免费av中文字幕在线| av黄色大香蕉| 国产精品福利在线免费观看| 一级毛片我不卡| 一个人看视频在线观看www免费| 视频区图区小说| 最近的中文字幕免费完整| 久久国产精品大桥未久av | 在线观看人妻少妇| 如日韩欧美国产精品一区二区三区 | 啦啦啦啦在线视频资源| 欧美精品一区二区大全| 亚洲高清免费不卡视频| 久久亚洲国产成人精品v| 国产精品国产三级国产av玫瑰| 成人二区视频| 亚洲精品视频女| 老司机影院成人| 欧美区成人在线视频| 精品少妇久久久久久888优播| 男人狂女人下面高潮的视频| 免费高清在线观看视频在线观看| 99久国产av精品国产电影| 色网站视频免费| 日日摸夜夜添夜夜爱| a级一级毛片免费在线观看| 五月玫瑰六月丁香| 国产亚洲av片在线观看秒播厂| 欧美亚洲 丝袜 人妻 在线| 精品国产国语对白av| 精品人妻熟女毛片av久久网站| 又粗又硬又长又爽又黄的视频| 天堂俺去俺来也www色官网| 亚洲欧美一区二区三区黑人 | 久久久久精品久久久久真实原创| 99九九在线精品视频 | 丰满人妻一区二区三区视频av| 亚洲婷婷狠狠爱综合网| 在线播放无遮挡| 中文精品一卡2卡3卡4更新| 最近2019中文字幕mv第一页| 国产一区二区三区综合在线观看 | 美女视频免费永久观看网站| 天天操日日干夜夜撸| 少妇猛男粗大的猛烈进出视频| 美女中出高潮动态图| 亚洲欧美成人综合另类久久久| 男人和女人高潮做爰伦理| 精品人妻熟女毛片av久久网站| 国产亚洲91精品色在线| 在线亚洲精品国产二区图片欧美 | 国产精品.久久久| 少妇的逼好多水| 丁香六月天网| 日韩不卡一区二区三区视频在线| 97精品久久久久久久久久精品| 免费大片黄手机在线观看| 亚洲内射少妇av| 乱码一卡2卡4卡精品| 久久久久精品久久久久真实原创| 国产精品99久久99久久久不卡 | 国产黄片视频在线免费观看| 亚洲性久久影院| 国产日韩一区二区三区精品不卡 | 亚洲欧美中文字幕日韩二区| 免费观看av网站的网址| av国产久精品久网站免费入址| 不卡视频在线观看欧美| 成人漫画全彩无遮挡| 80岁老熟妇乱子伦牲交| 久久精品国产自在天天线| 欧美精品亚洲一区二区| 在线观看av片永久免费下载| 大香蕉97超碰在线| 久久人妻熟女aⅴ| 成人特级av手机在线观看| 狠狠精品人妻久久久久久综合| 22中文网久久字幕| 91久久精品国产一区二区三区| 久久99热6这里只有精品| 亚洲av欧美aⅴ国产| 丝袜脚勾引网站| 欧美区成人在线视频| 精品卡一卡二卡四卡免费| 三级经典国产精品| 美女福利国产在线| 女人久久www免费人成看片| 黄色毛片三级朝国网站 | 97超视频在线观看视频| 亚洲精品一区蜜桃| 亚洲欧洲国产日韩| 三上悠亚av全集在线观看 | 一本一本综合久久| 免费观看av网站的网址| 一本大道久久a久久精品| 欧美激情极品国产一区二区三区 | 中文在线观看免费www的网站| 99久久精品热视频| 免费黄网站久久成人精品| 日韩一本色道免费dvd| 最近的中文字幕免费完整| 亚洲精品日韩在线中文字幕| 一级黄片播放器| 日本与韩国留学比较| 欧美精品亚洲一区二区| 一级,二级,三级黄色视频| 99热6这里只有精品| 久久精品久久久久久噜噜老黄| 免费在线观看成人毛片| av在线观看视频网站免费| 麻豆精品久久久久久蜜桃| 国产永久视频网站| 亚洲精品第二区| 亚洲精品自拍成人| 青春草国产在线视频| 亚洲av欧美aⅴ国产| 最近的中文字幕免费完整| 2022亚洲国产成人精品| 国产亚洲一区二区精品| 中文欧美无线码| 亚洲国产欧美日韩在线播放 | 性色avwww在线观看| 人体艺术视频欧美日本| 少妇猛男粗大的猛烈进出视频| 欧美日本中文国产一区发布| 国产成人免费无遮挡视频| 国产成人a∨麻豆精品| 国产黄片美女视频| 成年女人在线观看亚洲视频| 亚洲精品国产色婷婷电影| 日韩,欧美,国产一区二区三区|