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

    正常固結非飽和黏性土的三剪邊界面模型研究

    2021-07-06 07:01:40胡小榮蔡曉鋒
    工程力學 2021年6期
    關鍵詞:等量非飽和黏性

    胡小榮,蔡曉鋒,陳 昊,2

    (1.南昌大學建筑工程學院,南昌 330031;2.江西中煤建設集團有限公司,南昌 330001)

    基于傳統(tǒng)彈塑性理論建立的非飽和黏性土本構模型[1]需要解決屈服函數(shù)、流動法則和硬化規(guī)律這三個關鍵問題[2],對于像非飽和黏性土這樣土性復雜的力學介質(zhì),難度在于后兩個問題。而采用邊界面模型方法則可解決這兩個問題,其原理方法也較容易,且適用于土體受到單調(diào)、循環(huán)荷載作用,此外還可解決由于加載誘發(fā)導致土的各向異性問題,因此在各種飽和及非飽和土中得到了一定的應用[3?17]。其中在非飽和土方面,文獻[13]在邊界面模型理論的基礎上引入含水率,提出了重塑非飽和黃土邊界面動力本構模型;文獻[14 ? 15]利用改變屈服面大小來體現(xiàn)土的結構性,并基于非飽和結構性黏土提出循環(huán)荷載下的邊界面模型;文獻[16]把BBM模型和塑性硬化準則相結合并引入邊界面理論,提出了適用黏性土特性的彈塑性雙面模型。文獻[17]在邊界面模型框架下,結合結構損傷理論和非飽和土力學,并引入移動映射中心思想,建立了適用非飽和結構性土的本構模型;建立邊界面模型的關鍵所在就是要確定邊界面方程及其隨荷載變化的演化規(guī)律。如今被普遍應用的邊界面確定方法是:先根據(jù)土體在某一應力水平下的屈服面方程來確定初始邊界面,初始邊界面在隨后的加、卸載過程中,將隨荷載變化而不斷演化得到當前邊界面。越來越多的研究者把邊界面理論與p′?q平面上的臨界狀態(tài)橢圓形屈服面方程相結合。破壞應力比M=q/p是橢圓形邊界面在其演化過程中的關鍵參數(shù),一般是通過土的強度準則來確定。目前,研究者們提出了多種針對土體的強度準則:Druker-Prager準則、Mohr-Coulomb準則、Lade準則、SMP準則、雙剪統(tǒng)一屈服準則和臨界狀態(tài)橢圓屈服準則等。這些準則有的不能反映中間主應力效應,有些不能反映土體拉壓差效應,還有些不能反映土體的多樣性。而三剪統(tǒng)一強度準則可體現(xiàn)土體強度的區(qū)間效應、拉壓差效應和中間主應力效應[18],并且通過改變中主應力影響系數(shù)該準則可對其他任何準則作非線性逼近。三剪統(tǒng)一強度準則已被應用于分析單調(diào)加載作用下的超固結土[19? 20]和結構性土[21]的力學特性,還與邊界面模型理論相結合應用于對飽和含泥砂土的單調(diào)、循環(huán)荷載力學特性分析[22 ?23]。

    本文主要內(nèi)容是基于Bishop單應力變量[24]和Fredlund雙應力變量[25]提出的非飽和黏性土三剪強度準則[26],采用等量代換法和坐標平移法得到非飽和黏性土的三剪破壞應力比,并將其引入非飽和黏性土的修正劍橋模型以確定其為橢圓形屈服面方程,以橢圓形屈服面作為初始邊界面,采用徑向映射法則確定后續(xù)邊界面,根據(jù)插值函數(shù)得到加載面上的塑性模量,再根據(jù)彈塑性理論建立單應力變量和雙應力變量下等量代換法和坐標平移法的4個三剪彈塑性邊界面本構模型。為驗證模型的正確性,采用南昌地區(qū)非飽和紅黏土對該模型進行單調(diào)壓縮荷載下的三軸排水試驗驗證及真三軸數(shù)值模擬,采用非飽和黏質(zhì)粉土進行循環(huán)動三軸試驗驗證及非飽和黃土進行循環(huán)真三軸數(shù)值模擬。

    1 非飽和黏性土的三剪強度準則及三剪破壞應力比

    非飽和黏性土的單應力變量和雙應力變量三剪強度準則為[26]:

    1)單應力變量

    由式(1)和式(2)得等量代換法和坐標平移法的三剪破壞應力比為:

    1)單應力變量下采用等量代換法時的三剪破壞應力比為:

    2)單應力變量下采用坐標平移法時的三剪破壞應力比為:

    3)雙應力變量下采用等量代換法時的三剪破壞應力比為:

    4)雙應力變量下采用坐標平移法時的三剪破壞應力比為:

    圖1為非飽和黏性土三剪強度理論計算值和文獻[27]對壓實的黏土砂(SP-SC)的柔性真三軸試驗結果比較圖。

    由圖1可知,在雙應力變量狀態(tài)下,實驗數(shù)據(jù)與理論計算值基本一致,從而驗證了非飽和黏性土強度理論的正確性。盡管在單個應力變量狀態(tài)下的理論計算值有所出入,但是總體趨勢仍然相同。

    圖1 非飽和黏性土三剪強度準則理論值與柔性真三軸試驗結果比較Fig.1 Result comparisonsof flexible truetriaxial triple-shear test for unsaturated clays

    圖2為非飽和黏性土三剪強度理論計算值和文獻[28]對非飽和黃土的剛柔混合型真三軸試驗結果對比圖。

    由圖2可知,實驗數(shù)據(jù)與理論計算值均吻合良好,驗證了非飽和黏性土三剪強度準則的正確性。在相同條件下,單應力變量狀態(tài)和雙應力變量狀態(tài)下的理論計算值基本一致。

    圖2 非飽和黏性土三剪強度準則理論值與剛柔混合型真三軸試驗結果比較Fig.2 Comparisonsof theoretical valuesof thetriple-shear strength criterion and resultsof rigid-flexible truetriaxial test for unsaturated clays

    2 非飽和黏性土的三剪彈塑性邊界面模型

    將第1節(jié)所得三剪破壞應力比與原有非飽和土屈服面方程相結合得到非飽和黏性土的三剪屈服面方程,再結合邊界面理論建立非飽和黏性土三剪邊界面彈塑性本構模型。最后,結合非飽和黏性土三剪強度準則與邊界面彈塑性理論推導出剛度矩陣各元素的具體表達式。

    2.1 非飽和黏性土的三剪屈服面方程式

    根據(jù)文獻[29]提出的LC屈服線,將曲線方程改寫為:

    故非飽和狀態(tài)時土的修正劍橋模型屈服面方程可寫為:

    以及飽和狀態(tài)時土的應變硬化規(guī)律為:

    1)單應力變量下及雙應力變量下采用等量代換法時

    2.2 非飽和黏性土三剪邊界面彈塑性本構模型

    假定屈服面、邊界面和塑性勢面相似,則邊界面方程為:

    1)單應力變量及雙應力變量下采用等量代換法

    采用相關聯(lián)流動法則,即塑性勢面與屈服面重合,進而可知塑性流動方向Rij和加載方向Lij相一致,則有:

    采用如下徑向映射法則:

    1)單應力變量下及雙應力變量下采用等量代換法時

    其中:

    1)單應力變量下采用等量代換法時

    2)雙應力變量下采用等量代換法時

    3)單應力變量下采用坐標平移法時

    由邊界面方程一致性相容條件可得:

    其中:

    1)單應力變量下采用等量代換法時

    2)雙應力變量下采用等量代換法時

    3)單應力變量下及雙應力變量下采用坐標平移法時

    式中:單應力變量下時三剪破壞應力比M(θ,s)=M1(θ,s)、Λ=Λ1;雙應力變量下時三剪破壞應力比M(θ,s)=M2(θ,s)、 Λ=Λ2。

    將塑性體積應變增量公式和式(30)中損傷變量增量公式代入式(15)可得出:

    將各偏導等式代入式(31)可得出邊界面上像應力點對應的塑性模量如下:

    1)單應力變量下采用等量代換法時

    2)雙應力變量下采用等量代換法時

    3)單應力及雙應力變量下均采用坐標平移法時

    式中:單應力變量下時三剪破壞應力比M(θ,s)=M1(θ,s)、Λ=Λ1;雙應力變量下時三剪破壞應力比M(θ,s)=M2(θ,s)、 Λ=Λ2。

    已知像應力點的塑性模量求法,對應實際應力點的塑性模量則由塑性模量插值函數(shù)[30]確定:

    式中:Kp為當前應力點的塑性模量;Kp為像應力點塑性模量;H0為模型參數(shù),通常采用模擬試算方法取值。

    式(41)中各分量具體表達式見附錄。

    3 非飽和黏性土的三剪彈塑性邊界面模型驗證

    3.1 非飽和三軸試驗驗證

    為驗證單調(diào)荷載下該模型的合理性,用文獻[31? 32]中的南昌非飽和重塑紅土三軸壓縮排水試驗與本文所提本構模型計算結果做了對比。在排水條件下將凈圍壓分別為100 kPa、200 kPa和300 kPa和基質(zhì)吸力分別為0 kPa、20 kPa、50 kPa和80 kPa的模擬結果與試驗結果均進行對比,限于篇幅且對比結果類似,只列出凈圍壓為100 kPa及基質(zhì)吸力為50 kPa的對比結果,如圖3所示。

    圖3 凈圍壓為100 kPa及基質(zhì)吸力為50 kPa的對比曲線圖Fig.3 Curve comparisons of the net confining pressure of 100 kPa and the matric suction of 50 kPa

    由圖3(a)可知,在單應力變量下,等量代換法和與實驗數(shù)據(jù)差值最大不超過12%,而坐標平移法與實驗數(shù)據(jù)差值最大不超過14%。在雙應力變量下,等量代換法和與實驗數(shù)據(jù)差值最大不超過2%,而坐標平移法與實驗數(shù)據(jù)差值最大不超過5%。由圖3(b)可知,在單應力變量下,等量代換法和與實驗數(shù)據(jù)差值最大不超過19%,而坐標平移法與實驗數(shù)據(jù)差值最大不超過24%。在雙應力變量下,等量代換法和與實驗數(shù)據(jù)差值最大不超過10%,而坐標平移法與實驗數(shù)據(jù)差值最大不超過18%。實驗數(shù)據(jù)與四種模擬數(shù)據(jù)均符合較好,在軸向應變超過5%左右時,單應力變量下的模擬效果相對較差。與采用單應力變量的模擬數(shù)據(jù)相比,雙應力變量的模擬數(shù)據(jù)更符合實驗結果,這是因為單應力變量法只考慮單一的有效應力變量使得無法全面表征非飽和土的力學特性,而雙應力變量法考慮把吸力和凈應力作為應力狀態(tài)量能較全面地表征非飽和土的力學特性;與采用坐標平移法得到的模擬數(shù)據(jù)相比,采用等量代換法得到的模擬數(shù)據(jù)更符合實驗結果。驗證了所提本構模型在單調(diào)荷載作用下的適用性。

    3.2 非飽和循環(huán)動三軸試驗驗證

    為驗證循環(huán)荷載下該模型的合理性,用文獻[16]中非飽和黏質(zhì)粉土在基質(zhì)吸力分別為50 kPa和100 kPa下的循環(huán)動三軸試驗與本文所提本構模型計算結果做了對比。限于篇幅且對比結果類似,圖4為排水條件下基質(zhì)吸力為50 kPa時的對比結果圖。

    由圖4(a)可知,在加載初期,單應力變量法下等量代換法和實驗數(shù)據(jù)差值最大不超過10%,坐標平移法和實驗數(shù)據(jù)差值最大不超過16%。在循環(huán)荷載過程中,單應力變量法下等量代換法和實驗數(shù)據(jù)差值最大不超過6%,坐標平移法和實驗數(shù)據(jù)差值最大不超過9%。由圖4(b)可知,在加載初期,雙應力變量法下等量代換法和實驗數(shù)據(jù)差值最大不超過20%,坐標平移法和實驗數(shù)據(jù)差值最大不超過21%。在循環(huán)荷載過程中,雙應力變量法下等量代換法和實驗數(shù)據(jù)差值最大不超過3%,坐標平移法和實驗數(shù)據(jù)差值最大不超過5%。所建本構模型均能夠較好地描述非飽和黏質(zhì)粉土在循環(huán)動荷載下的應力-應變變化規(guī)律。雙應力變量下建立的模型模擬結果相對更接近實驗結果,采用等量代換法建立的模型模擬結果相對更接近實驗結果;循環(huán)荷載加載初期時,土體的塑性變形發(fā)展相對較快,但隨著加載時間的增加,土體產(chǎn)生的塑性變形逐漸變小。

    圖4 凈圍壓為100 kPa及基質(zhì)吸力為50 kPa時偏應力與軸向應變關系試驗與計算結果對比圖Fig.4 Comparisons between experimental and calculated resultsof the deviatoric stressvs. the axial strain at the net confining pressureof 100 kPa and the matric suction of 50 kPa

    4 本構模型的真三軸模擬

    4.1 單調(diào)壓縮荷載下的真三軸模擬

    模擬施加的應力路徑:1)首先對六面體單元的三個方向同時施加最小主應力σ3,固結排水;2)固結完成后,一個方向應力σ3保持不變,其余兩個方向應力增至中間主應力σ2,固結排水;3)一個方向應力σ2繼續(xù)保持不變,另一個方向應力由σ2逐漸增加,即為最大主應力σ1,在σ1增大過程中,對土樣進行排水數(shù)值模擬。限于篇幅且結果類似,只列出了雙應力變量下等量代換法本構模型對南昌非飽和紅黏土在不同中間主應力影響系數(shù)下的模擬結果,如圖5所示,其中凈圍壓和基質(zhì)吸力分別為100 kPa和50 kPa。

    由圖5(a)可知,中間主應力影響系數(shù)b=0時與b=0.25時模擬結果差值最大為2%,b=0時與b=0.5時模擬結果差值最大為5%,b=0時與b=0.75時模擬結果差值最大為10%,b=0時與b=1時模擬結果差值最大為12%;由圖5(b)可知,中間主應力影響系數(shù)b=0時與b=0.25時模擬結果差值最大為2.5%,b=0時與b=0.5時模擬結果差值最大為5%,b=0時與b=0.75時模擬結果差值最大為7.5%,b=0時與b=1時模擬結果差值最大為10%。在真三軸單調(diào)加載初期,剪應力與體應變增長較快,并最終趨于穩(wěn)定狀態(tài);同等條件下,隨著中主應力影響系數(shù)b值增大,剪應力與體應變也隨之增大,表明中主應力影響系數(shù)與土體抗剪強度為正相關關系。

    圖5 凈圍壓為100 kPa及基質(zhì)吸力為50 kPa時中間主應力影響曲線圖Fig.5 Influence curves of the intermediate principal stress when the net confining pressure is100 kPa and thematric suction is50 kPa

    圖6為基質(zhì)吸力50 kPa、中主應力影響系數(shù)b=0.5和凈圍壓100 kPa時四種本構模型的真三軸模擬結果。

    由圖6(a)可知,單應力變量坐標平移法與單應力變量等量代換法模擬數(shù)據(jù)差值最大不超過5%,單應力變量坐標平移法與雙應力變量等量代換法模擬數(shù)據(jù)差值最大不超過17%,單應力變量坐標平移法與雙應力變量坐標平移法模擬數(shù)據(jù)差值最大不超過21%;由圖6(b)可知,單應力變量坐標平移法與單應力變量等量代換法模擬數(shù)據(jù)差值最大不超過6%,單應力變量坐標平移法與雙應力變量等量代換法模擬數(shù)據(jù)差值最大不超過21%,單應力變量坐標平移法與雙應力變量坐標平移法模擬數(shù)據(jù)差值最大不超過27%。相同應力變量下,采用等量代換法和坐標平移法的模擬差異很??;土體在雙應力變量下的模擬抗剪強度相對單應力變量下較大。

    圖6 四種本構模型的真三軸模擬對比結果圖Fig.6 Comparisons of the truetriaxial simulated resultsof thefour constitutivemodels

    4.2 循環(huán)壓縮荷載下的真三軸模擬

    模擬施加的應力路徑:1)首先對六面體單元的2個側向和1個豎向同時施加主應力σ3,固結排水;2)固結完成后一個水平側向應力σ3保持不變,其余兩個方向應力增至σ2,固結排水;3)另一個側向應力σ2繼續(xù)保持不變,豎向施加以σ2為初始值的循環(huán)荷載σ1,固結排水。限于篇幅且雙應力變量等量代換法建立的模型模擬結果相對更接近實驗結果,故本文圖7~圖9主要列出雙應力變量下采用等量代換法的模型模擬結果來分析最小主應力、中間主應力以及應力幅值變化時對土體產(chǎn)生的影響,其余模型規(guī)律與之相似。圖7~圖10的計算參數(shù)引自文獻[33? 34]中的非飽和黃土動三軸試驗。圖7是幅值為75 kPa、σ3=25 kPa、σ2分別為40 kPa和50 kPa的模擬結果;圖8是幅值為75 kPa、σ2=50 kPa、σ3分別為15 kPa和25 kPa的模擬結果;圖9是σ2=50 kPa、σ3=25 kPa、幅值分別為75 kPa和100 kPa的模擬結果。

    圖8 幅值為75 kPa、σ2=50 kPa、σ3分別為15 kPa和25 kPa的模擬結果Fig.8 Simulated results with amplitude of 75 kPa,σ2=50 kPa and σ3 of 15 kPa and 25 kPa,respectively

    由圖7~圖9可知,非飽和黃土的塑性變形先期發(fā)展較快,后逐漸減??;若僅分別增加最小主應力或中間主應力,均會使土體的抗剪強度增大;應力大小不變,最小和中間主應力不變,荷載振幅增大,非飽和黃土的應變隨之增大。

    圖7 幅值為75 kPa、σ3=25 kPa、σ2分別為40 kPa和50 kPa的模擬結果Fig.7 Simulated results with amplitude of 75 kPa,σ3=25 kPa and σ2 of 40 kPa and 50 kPa,respectively

    圖9 σ2=50 kPa、σ3=25 kPa、幅值分別為75 kPa和100 kPa的模擬結果Fig.9 Simulated results with σ2 =50 kPa,σ3= 25 kPa and amplitudesof 75 kPa and 100 kPa, respectively

    圖10為荷載幅值、σ3和σ2分別為75 kPa、25 kPa和50 kPa時4種彈塑性本構模型的真三軸循環(huán)荷載模擬結果。

    由圖10(a)可知,雙應力變量等量代換法與雙應力變量坐標平移法滯回圈模擬數(shù)據(jù)差值最大不超過5%,雙應力變量等量代換法與單應力變量坐標平移法滯回圈模擬數(shù)據(jù)差值最大不超過7%,雙應力變量等量代換法與單應力變量等量代換法滯回圈模擬數(shù)據(jù)差值最大不超過12%;由圖10(b)可知,雙應力變量等量代換法與雙應力變量坐標平移法模擬數(shù)據(jù)差值最大不超過3%,雙應力變量等量代換法與單應力變量坐標平移法模擬數(shù)據(jù)差值最大不超過7%,雙應力變量等量代換法與單應力變量等量代換法模擬數(shù)據(jù)差值最大不超過12%。相同應力變量下,采用等量代換法和坐標平移法的模擬差別并不大;在單應力變量下模擬非飽和黃土的塑性變形大于雙應力變量下非飽和黃土的塑性變形。

    圖10 四種模型在真三軸條件下的模擬對比結果圖Fig.10 Simulated comparisons of the four models under the true triaxial conditions

    5 結論

    提出了非飽和黏性土的三剪強度準則,并采用等量代換法與坐標平移法推導了三剪破壞應力比。結合邊界面理論建立了非飽和黏性土的三剪彈塑性邊界面本構模型并對其做了單調(diào)和循環(huán)壓縮荷載試驗驗證。結論如下:

    (1)所提出的三剪破壞應力比有效克服了原修正劍橋模型中破壞應力比為定值的不足,另外還能反映非飽和黏性土黏聚力與全應力狀態(tài)下土體的應力-應變特性的影響。

    (2)所建立的三剪邊界面模型均能較好反映非飽和黏性土中間主應力效應、拉壓差效應和區(qū)間效應,可用于非飽和黏性土在單調(diào)和循環(huán)荷載作用下的彈塑性分析。

    (3)以南昌地區(qū)非飽和重塑紅黏土為試驗研究對象進行非飽和三軸試驗驗證,并采用文獻中的非飽和黏質(zhì)粉土進行非飽和循環(huán)三軸試驗,結果表明本文所提模型在單調(diào)、循環(huán)荷載下均有較好的預測能力。其中雙應力變量等量代換法與實驗結果較為接近,在單調(diào)荷載下與實驗結果差值最大不超過10%,在循環(huán)荷載下與實驗結果差值最大不超過3%。

    (4)對南昌地區(qū)非飽和重塑紅黏土和文獻中非飽和黃土分別進行單調(diào)、循環(huán)荷載下真三軸模擬試驗驗證,結果也表明在真三軸條件下本文所提模型均能較好地模擬單調(diào)、循環(huán)荷載下的力學特性。其中在單調(diào)加載模擬試驗中,中主應力影響系數(shù)或基質(zhì)吸力均與非飽和土體抗剪強度呈正相關關系;且雙應力變量下的模擬抗剪強度相對較大;在循環(huán)加載模擬實驗中,最小主應力或中間主應力均與土體的抗剪強度呈正相關關系;荷載振幅與非飽和土體的應變成正比;土體在單應力變量下模擬產(chǎn)生的塑性變形相對較大。

    1)單應力變量下采用等量代換法

    2)雙應力變量下采用等量代換法

    3)單應力變量下及雙應力變量下采用坐標平移法

    G、L具體表達式如下:

    1)單應力變量下及雙應力變量下采用等量代換法

    2)單應力變量下采用坐標平移法

    3)雙應力變量下采用坐標平移法

    其中:

    猜你喜歡
    等量非飽和黏性
    等量代換
    玩轉等量同減
    曹沖稱象和等量變換
    非飽和原狀黃土結構強度的試驗研究
    工程與建設(2019年1期)2019-09-03 01:12:24
    富硒產(chǎn)業(yè)需要強化“黏性”——安康能否玩轉“硒+”
    當代陜西(2019年14期)2019-08-26 09:41:56
    如何運用播音主持技巧增強受眾黏性
    傳媒評論(2019年4期)2019-07-13 05:49:28
    等量代換
    玩油灰黏性物成網(wǎng)紅
    華人時刊(2017年17期)2017-11-09 03:12:03
    非飽和多孔介質(zhì)應力滲流耦合分析研究
    非飽和土基坑剛性擋墻抗傾覆設計與參數(shù)分析
    免费人妻精品一区二区三区视频| 中国美白少妇内射xxxbb| 人妻系列 视频| 新久久久久国产一级毛片| 亚洲国产日韩一区二区| 欧美bdsm另类| 日日摸夜夜添夜夜爱| 人妻少妇偷人精品九色| 日韩精品免费视频一区二区三区 | 免费高清在线观看视频在线观看| 欧美人与善性xxx| av线在线观看网站| av不卡在线播放| 亚洲va在线va天堂va国产| 国内精品宾馆在线| 日本色播在线视频| .国产精品久久| 精品人妻一区二区三区麻豆| 亚洲精品国产av蜜桃| 在线亚洲精品国产二区图片欧美 | 狠狠精品人妻久久久久久综合| 亚洲精品一区蜜桃| 免费黄频网站在线观看国产| 青春草视频在线免费观看| 少妇被粗大的猛进出69影院 | 搡老乐熟女国产| 天堂8中文在线网| 王馨瑶露胸无遮挡在线观看| 美女内射精品一级片tv| 国产黄色视频一区二区在线观看| 嘟嘟电影网在线观看| 国产精品麻豆人妻色哟哟久久| 精品视频人人做人人爽| 亚洲性久久影院| 久久久久久久亚洲中文字幕| 亚洲精品,欧美精品| 街头女战士在线观看网站| 国产av码专区亚洲av| 色网站视频免费| 国产成人精品久久久久久| 伊人亚洲综合成人网| 亚洲精品亚洲一区二区| 精品亚洲乱码少妇综合久久| 看非洲黑人一级黄片| 一级毛片aaaaaa免费看小| 中文字幕制服av| 极品教师在线视频| 国产亚洲av片在线观看秒播厂| 亚洲人成网站在线观看播放| 黄色一级大片看看| 熟女av电影| 亚洲第一av免费看| 曰老女人黄片| 国产视频内射| 免费黄网站久久成人精品| 国产精品嫩草影院av在线观看| 日本黄色片子视频| 一级黄片播放器| 又爽又黄a免费视频| 免费黄网站久久成人精品| 国产午夜精品一二区理论片| 青春草视频在线免费观看| 久久久久久久久久久免费av| 国内少妇人妻偷人精品xxx网站| 久久久久久伊人网av| 青青草视频在线视频观看| 王馨瑶露胸无遮挡在线观看| 亚洲国产成人一精品久久久| 国产又色又爽无遮挡免| 国产在线视频一区二区| 欧美日韩综合久久久久久| 久久人妻熟女aⅴ| 国产白丝娇喘喷水9色精品| 三级国产精品片| 亚洲精品中文字幕在线视频 | 国产精品女同一区二区软件| av一本久久久久| 丰满乱子伦码专区| 免费av中文字幕在线| 综合色丁香网| 国内精品宾馆在线| 欧美日韩亚洲高清精品| 国产精品麻豆人妻色哟哟久久| 国产精品.久久久| 高清午夜精品一区二区三区| 街头女战士在线观看网站| 国产亚洲91精品色在线| 免费观看性生交大片5| 久久 成人 亚洲| 日本欧美国产在线视频| 国产精品蜜桃在线观看| 成人美女网站在线观看视频| 亚洲国产日韩一区二区| 日韩av在线免费看完整版不卡| 国产熟女午夜一区二区三区 | 黑人高潮一二区| 国内揄拍国产精品人妻在线| 国产永久视频网站| 亚洲天堂av无毛| 精品人妻偷拍中文字幕| 视频中文字幕在线观看| 天天躁夜夜躁狠狠久久av| 国产在线视频一区二区| 最近中文字幕2019免费版| 男女边摸边吃奶| 日韩不卡一区二区三区视频在线| 一级毛片 在线播放| 亚洲情色 制服丝袜| 精品卡一卡二卡四卡免费| 午夜激情福利司机影院| 老熟女久久久| 97超碰精品成人国产| 日韩在线高清观看一区二区三区| av福利片在线观看| 国产一区二区三区av在线| 久久久国产一区二区| 六月丁香七月| 秋霞在线观看毛片| 黄色欧美视频在线观看| 欧美+日韩+精品| 国产av精品麻豆| 极品人妻少妇av视频| 精品卡一卡二卡四卡免费| 韩国高清视频一区二区三区| 国产中年淑女户外野战色| 免费观看性生交大片5| 欧美性感艳星| 欧美日韩精品成人综合77777| 91久久精品国产一区二区三区| 免费大片18禁| 国产黄色视频一区二区在线观看| 在线观看国产h片| 久久99热这里只频精品6学生| 国产乱来视频区| 免费人成在线观看视频色| 青春草视频在线免费观看| 亚洲精品视频女| 午夜日本视频在线| 人妻制服诱惑在线中文字幕| 亚洲中文av在线| 国产精品久久久久久精品古装| av国产精品久久久久影院| 国产精品久久久久久av不卡| 国产精品99久久99久久久不卡 | 国产成人精品无人区| 日本与韩国留学比较| 日本猛色少妇xxxxx猛交久久| 天堂8中文在线网| 国产乱来视频区| 蜜桃在线观看..| 大话2 男鬼变身卡| 国产成人91sexporn| 中文乱码字字幕精品一区二区三区| 久久久国产精品麻豆| 女的被弄到高潮叫床怎么办| 精品一区二区三卡| 精品人妻熟女毛片av久久网站| 亚洲欧美一区二区三区黑人 | 久久99蜜桃精品久久| 各种免费的搞黄视频| 狂野欧美激情性xxxx在线观看| 久久国内精品自在自线图片| 亚洲人成网站在线播| 全区人妻精品视频| 日韩强制内射视频| 91精品国产九色| 国产精品一区二区三区四区免费观看| 黑人高潮一二区| 国产视频内射| 我的老师免费观看完整版| 一级a做视频免费观看| 国产日韩欧美视频二区| 色网站视频免费| 欧美少妇被猛烈插入视频| 国产精品伦人一区二区| 日本-黄色视频高清免费观看| 免费av中文字幕在线| 久久久久久久大尺度免费视频| 大陆偷拍与自拍| 一边亲一边摸免费视频| 久久精品国产a三级三级三级| 欧美三级亚洲精品| 国产av国产精品国产| 成人美女网站在线观看视频| 最新中文字幕久久久久| 日韩 亚洲 欧美在线| 亚洲精华国产精华液的使用体验| 少妇人妻 视频| 在线 av 中文字幕| 少妇高潮的动态图| 又粗又硬又长又爽又黄的视频| 国产精品国产三级国产专区5o| 夫妻性生交免费视频一级片| 色哟哟·www| 丝袜在线中文字幕| 成人毛片60女人毛片免费| 国产高清国产精品国产三级| av在线app专区| 最新的欧美精品一区二区| 人人妻人人添人人爽欧美一区卜| 国产色婷婷99| 日韩不卡一区二区三区视频在线| 国产中年淑女户外野战色| 欧美最新免费一区二区三区| 成人影院久久| av黄色大香蕉| 少妇人妻一区二区三区视频| 少妇 在线观看| 久久久久久久久久人人人人人人| 欧美日韩视频高清一区二区三区二| 十八禁高潮呻吟视频 | 水蜜桃什么品种好| 国产精品国产av在线观看| 成年美女黄网站色视频大全免费 | 久久久久久久亚洲中文字幕| 日本-黄色视频高清免费观看| 亚洲欧洲国产日韩| 99九九线精品视频在线观看视频| 亚洲av免费高清在线观看| 极品少妇高潮喷水抽搐| 日日摸夜夜添夜夜添av毛片| 黑丝袜美女国产一区| 国内揄拍国产精品人妻在线| 少妇猛男粗大的猛烈进出视频| 久久久久久久大尺度免费视频| 少妇高潮的动态图| 亚洲国产精品专区欧美| 久久国内精品自在自线图片| 丁香六月天网| 国产有黄有色有爽视频| 国产高清三级在线| 亚洲婷婷狠狠爱综合网| 国产爽快片一区二区三区| 在线天堂最新版资源| 97在线视频观看| 日本黄色片子视频| 免费大片黄手机在线观看| 一本色道久久久久久精品综合| 97精品久久久久久久久久精品| 国产毛片在线视频| 最新的欧美精品一区二区| 91精品伊人久久大香线蕉| 蜜桃久久精品国产亚洲av| 国产有黄有色有爽视频| 91久久精品国产一区二区三区| 中文字幕亚洲精品专区| 高清午夜精品一区二区三区| 伦理电影免费视频| 人妻系列 视频| 国产欧美亚洲国产| 日韩不卡一区二区三区视频在线| 好男人视频免费观看在线| 久久婷婷青草| 免费播放大片免费观看视频在线观看| 亚洲婷婷狠狠爱综合网| 久久这里有精品视频免费| 黄片无遮挡物在线观看| 中文字幕av电影在线播放| 蜜桃久久精品国产亚洲av| 看十八女毛片水多多多| av又黄又爽大尺度在线免费看| 亚洲欧美中文字幕日韩二区| 日韩一区二区视频免费看| 久久免费观看电影| 观看美女的网站| 婷婷色av中文字幕| 久久久国产精品麻豆| 亚洲精品乱久久久久久| 亚洲欧美日韩另类电影网站| 80岁老熟妇乱子伦牲交| 五月伊人婷婷丁香| 久久精品久久精品一区二区三区| 国产爽快片一区二区三区| 亚洲av日韩在线播放| √禁漫天堂资源中文www| 亚洲精品456在线播放app| 性高湖久久久久久久久免费观看| 汤姆久久久久久久影院中文字幕| 天美传媒精品一区二区| 少妇精品久久久久久久| 国产黄色免费在线视频| 大码成人一级视频| 2021少妇久久久久久久久久久| av在线老鸭窝| 肉色欧美久久久久久久蜜桃| 国产国拍精品亚洲av在线观看| 男女无遮挡免费网站观看| 国产探花极品一区二区| 深夜a级毛片| av国产精品久久久久影院| 精品人妻一区二区三区麻豆| 边亲边吃奶的免费视频| 免费av中文字幕在线| 国产精品久久久久久久久免| 三级经典国产精品| 亚洲综合色惰| 亚洲国产色片| 18+在线观看网站| 国产 一区精品| 大片免费播放器 马上看| av在线app专区| 亚洲av中文av极速乱| 国产黄片视频在线免费观看| 精品亚洲乱码少妇综合久久| 亚洲中文av在线| 最新中文字幕久久久久| 91久久精品国产一区二区成人| 国产爽快片一区二区三区| 亚洲第一av免费看| 肉色欧美久久久久久久蜜桃| 亚洲内射少妇av| 特大巨黑吊av在线直播| 成年美女黄网站色视频大全免费 | 少妇人妻久久综合中文| 女性被躁到高潮视频| 人妻一区二区av| 国产亚洲91精品色在线| 欧美 日韩 精品 国产| 免费高清在线观看视频在线观看| 久久精品国产亚洲网站| 三级国产精品欧美在线观看| 高清欧美精品videossex| 久久久久久久久久人人人人人人| 亚洲婷婷狠狠爱综合网| 欧美日韩综合久久久久久| 久久韩国三级中文字幕| 女人久久www免费人成看片| 久久人人爽av亚洲精品天堂| 午夜免费男女啪啪视频观看| 我要看日韩黄色一级片| 亚洲性久久影院| 中国国产av一级| 另类亚洲欧美激情| 美女国产视频在线观看| 91成人精品电影| 美女国产视频在线观看| 国产高清三级在线| 大码成人一级视频| 久久精品熟女亚洲av麻豆精品| 少妇人妻久久综合中文| 国产亚洲av片在线观看秒播厂| 少妇人妻久久综合中文| 亚洲国产精品国产精品| 又大又黄又爽视频免费| 欧美日韩亚洲高清精品| 人人妻人人澡人人看| 日韩 亚洲 欧美在线| 亚洲av欧美aⅴ国产| 久久国内精品自在自线图片| 亚洲av日韩在线播放| 欧美 亚洲 国产 日韩一| 久久久久久久久久久免费av| 国产免费一区二区三区四区乱码| 国产av码专区亚洲av| 亚洲国产最新在线播放| 精品熟女少妇av免费看| 热99国产精品久久久久久7| 亚洲第一区二区三区不卡| 在线观看免费高清a一片| 99国产精品免费福利视频| 青青草视频在线视频观看| 国产成人一区二区在线| 一区二区三区四区激情视频| h日本视频在线播放| 卡戴珊不雅视频在线播放| 亚洲国产日韩一区二区| 男男h啪啪无遮挡| 成人亚洲精品一区在线观看| 性色avwww在线观看| 国产成人一区二区在线| 亚洲色图综合在线观看| 国产在线男女| 丰满人妻一区二区三区视频av| 国产在线一区二区三区精| 亚洲伊人久久精品综合| 丁香六月天网| 少妇人妻久久综合中文| 欧美3d第一页| 制服丝袜香蕉在线| 最近中文字幕高清免费大全6| 国产免费又黄又爽又色| 国产成人精品久久久久久| 久久久久人妻精品一区果冻| 欧美成人精品欧美一级黄| 亚洲成人av在线免费| 国产高清国产精品国产三级| 伊人久久国产一区二区| 蜜臀久久99精品久久宅男| 国产综合精华液| 精品亚洲成国产av| 男女边摸边吃奶| 美女cb高潮喷水在线观看| 久久av网站| 男女边吃奶边做爰视频| 亚洲精品中文字幕在线视频 | 精品一区在线观看国产| 成年av动漫网址| 在线观看av片永久免费下载| 中文字幕av电影在线播放| 午夜免费观看性视频| 97超视频在线观看视频| 国产精品国产三级国产专区5o| 人体艺术视频欧美日本| 日本vs欧美在线观看视频 | 免费看不卡的av| 国产极品粉嫩免费观看在线 | 国产美女午夜福利| 欧美97在线视频| 美女主播在线视频| 日韩一区二区三区影片| 午夜福利网站1000一区二区三区| 国模一区二区三区四区视频| 尾随美女入室| 亚洲精品国产成人久久av| 精品99又大又爽又粗少妇毛片| 美女大奶头黄色视频| 免费大片18禁| 一级爰片在线观看| 交换朋友夫妻互换小说| 亚洲,一卡二卡三卡| 男人添女人高潮全过程视频| 免费人妻精品一区二区三区视频| 黄色视频在线播放观看不卡| 日本-黄色视频高清免费观看| 精品卡一卡二卡四卡免费| 黄色配什么色好看| 青青草视频在线视频观看| 蜜桃在线观看..| 自拍欧美九色日韩亚洲蝌蚪91 | 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品中文字幕在线视频 | 美女福利国产在线| av黄色大香蕉| 亚洲美女视频黄频| 晚上一个人看的免费电影| 日韩欧美 国产精品| 国产精品嫩草影院av在线观看| 伊人久久精品亚洲午夜| 美女xxoo啪啪120秒动态图| 国产成人午夜福利电影在线观看| 国产极品粉嫩免费观看在线 | 国产精品久久久久久精品电影小说| 三级经典国产精品| 免费观看av网站的网址| videossex国产| 老司机影院毛片| 亚洲情色 制服丝袜| 国产午夜精品一二区理论片| 国产精品伦人一区二区| 久久人人爽人人爽人人片va| 丝瓜视频免费看黄片| 在线看a的网站| 大香蕉97超碰在线| 久久久久网色| 免费黄网站久久成人精品| 国产有黄有色有爽视频| 国产精品久久久久成人av| 欧美区成人在线视频| 国产精品人妻久久久影院| 一级黄片播放器| 国产黄色免费在线视频| 精品人妻熟女av久视频| 日本91视频免费播放| 亚洲综合色惰| 精品亚洲成国产av| 中国三级夫妇交换| 国产91av在线免费观看| 亚洲天堂av无毛| 在线观看www视频免费| 亚洲av日韩在线播放| 亚洲四区av| 天堂中文最新版在线下载| 国产精品久久久久成人av| 久久99一区二区三区| 国产精品一区www在线观看| 男女啪啪激烈高潮av片| 亚洲自偷自拍三级| 久久久a久久爽久久v久久| 天天躁夜夜躁狠狠久久av| 最近手机中文字幕大全| 自线自在国产av| 亚洲av成人精品一二三区| 国产爽快片一区二区三区| 国产亚洲5aaaaa淫片| 一级毛片我不卡| 免费观看在线日韩| 大香蕉97超碰在线| 国产免费一级a男人的天堂| 最近手机中文字幕大全| 国产av码专区亚洲av| 嫩草影院入口| 女人久久www免费人成看片| xxx大片免费视频| 日韩电影二区| 久久久精品94久久精品| 80岁老熟妇乱子伦牲交| 狂野欧美激情性xxxx在线观看| 亚洲欧美成人精品一区二区| 久久综合国产亚洲精品| 久久国产精品男人的天堂亚洲 | 亚洲av成人精品一区久久| 国产真实伦视频高清在线观看| 一级毛片久久久久久久久女| 一级毛片aaaaaa免费看小| av一本久久久久| 国产欧美日韩精品一区二区| 我要看日韩黄色一级片| 蜜桃在线观看..| 综合色丁香网| 精品酒店卫生间| 久久午夜福利片| 午夜福利在线观看免费完整高清在| 欧美日韩国产mv在线观看视频| 成人免费观看视频高清| 9色porny在线观看| 中文字幕久久专区| 2022亚洲国产成人精品| 久久久久国产精品人妻一区二区| 大话2 男鬼变身卡| 婷婷色综合www| 肉色欧美久久久久久久蜜桃| 下体分泌物呈黄色| 蜜臀久久99精品久久宅男| 亚洲成人一二三区av| 久久久午夜欧美精品| videossex国产| 国产熟女欧美一区二区| 美女内射精品一级片tv| 夜夜看夜夜爽夜夜摸| 女人精品久久久久毛片| 99热6这里只有精品| 国产精品偷伦视频观看了| 美女视频免费永久观看网站| 日韩中字成人| 自线自在国产av| 97超视频在线观看视频| 日韩视频在线欧美| 一级a做视频免费观看| 亚洲精品456在线播放app| 五月伊人婷婷丁香| 自线自在国产av| 在线观看免费日韩欧美大片 | 日韩一区二区视频免费看| 噜噜噜噜噜久久久久久91| 99re6热这里在线精品视频| 好男人视频免费观看在线| 永久免费av网站大全| 国产视频内射| 只有这里有精品99| 美女视频免费永久观看网站| 欧美bdsm另类| 国产深夜福利视频在线观看| 国产一区亚洲一区在线观看| 久热这里只有精品99| 天堂8中文在线网| 伦理电影大哥的女人| 国产成人免费无遮挡视频| 午夜福利在线观看免费完整高清在| 欧美xxxx性猛交bbbb| 日日啪夜夜撸| 老熟女久久久| 国产成人精品婷婷| 色5月婷婷丁香| 在线观看免费高清a一片| 亚洲成人av在线免费| 五月天丁香电影| av在线观看视频网站免费| 国产中年淑女户外野战色| 久久久亚洲精品成人影院| 麻豆乱淫一区二区| 日韩精品有码人妻一区| 蜜桃久久精品国产亚洲av| xxx大片免费视频| 在线播放无遮挡| 少妇猛男粗大的猛烈进出视频| 亚洲欧美一区二区三区黑人 | 日本欧美国产在线视频| 亚洲精品亚洲一区二区| 国产精品久久久久久精品古装| 日本wwww免费看| 国产免费又黄又爽又色| 偷拍熟女少妇极品色| 在线精品无人区一区二区三| 最近中文字幕高清免费大全6| 中文欧美无线码| 一级av片app| 亚洲美女视频黄频| 黄色欧美视频在线观看| 国产精品一区二区在线不卡| 国产欧美日韩精品一区二区| 亚洲精品国产色婷婷电影| 欧美激情极品国产一区二区三区 | 亚洲国产av新网站| 精品一区二区三卡| 一级毛片 在线播放| 亚洲精品一二三| 美女内射精品一级片tv| 99视频精品全部免费 在线| 一区二区三区乱码不卡18| 日韩人妻高清精品专区| 中文字幕亚洲精品专区| 国产成人免费观看mmmm| 国产成人精品一,二区| 久久综合国产亚洲精品| 九九在线视频观看精品| 国产精品国产三级国产专区5o| 亚洲综合色惰| 久久亚洲国产成人精品v| 视频中文字幕在线观看| 亚洲欧美成人精品一区二区| 国产精品久久久久久av不卡| 亚洲欧美精品自产自拍| 国产成人aa在线观看| 我的老师免费观看完整版| 色婷婷久久久亚洲欧美|