• <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ù)分析
    香蕉丝袜av| 日日摸夜夜添夜夜添小说| 国产亚洲欧美精品永久| 一级作爱视频免费观看| 男女之事视频高清在线观看| 99在线人妻在线中文字幕| 男女床上黄色一级片免费看| 亚洲欧美日韩无卡精品| 波多野结衣高清无吗| 757午夜福利合集在线观看| 久久午夜综合久久蜜桃| av片东京热男人的天堂| 国产成人av教育| 亚洲成av人片免费观看| 久久久久久久精品吃奶| 亚洲av电影在线进入| 淫秽高清视频在线观看| 1024香蕉在线观看| 亚洲成人精品中文字幕电影| 韩国精品一区二区三区| 久久婷婷成人综合色麻豆| 亚洲中文字幕日韩| 黑人巨大精品欧美一区二区mp4| 国产精品国产高清国产av| 亚洲专区字幕在线| 99久久99久久久精品蜜桃| 黑人操中国人逼视频| 亚洲在线自拍视频| 久久伊人香网站| www.精华液| 高清在线国产一区| 久久久久久亚洲精品国产蜜桃av| 午夜精品久久久久久毛片777| 久久久久九九精品影院| 久久久久国内视频| 中亚洲国语对白在线视频| 欧美黄色片欧美黄色片| 国产免费av片在线观看野外av| 免费少妇av软件| 久久国产精品男人的天堂亚洲| 欧美乱码精品一区二区三区| 午夜久久久在线观看| 欧美中文日本在线观看视频| 久久这里只有精品19| 99在线视频只有这里精品首页| 欧美乱妇无乱码| 色播亚洲综合网| 狠狠狠狠99中文字幕| 成在线人永久免费视频| 国产一区二区三区在线臀色熟女| 亚洲熟妇熟女久久| 90打野战视频偷拍视频| 91av网站免费观看| 日本a在线网址| 亚洲精华国产精华精| 桃红色精品国产亚洲av| 久久午夜亚洲精品久久| 日韩高清综合在线| 一级片免费观看大全| 亚洲五月天丁香| 中亚洲国语对白在线视频| 久久天躁狠狠躁夜夜2o2o| 18禁裸乳无遮挡免费网站照片 | 亚洲国产高清在线一区二区三 | 两个人看的免费小视频| www国产在线视频色| 人妻久久中文字幕网| 日韩欧美国产在线观看| 欧美中文综合在线视频| 色综合欧美亚洲国产小说| 最新在线观看一区二区三区| 在线av久久热| 成熟少妇高潮喷水视频| 国产成人啪精品午夜网站| 精品久久蜜臀av无| 精品一区二区三区四区五区乱码| 一区在线观看完整版| 久久香蕉激情| 桃红色精品国产亚洲av| 久久香蕉激情| 午夜精品久久久久久毛片777| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久国产成人免费| 国产成人一区二区三区免费视频网站| 国产av又大| 亚洲色图 男人天堂 中文字幕| 真人一进一出gif抽搐免费| 久久久久九九精品影院| 日本a在线网址| 久久精品国产亚洲av高清一级| 搡老妇女老女人老熟妇| 亚洲七黄色美女视频| 国产精品av久久久久免费| 欧美在线一区亚洲| 丝袜人妻中文字幕| 人成视频在线观看免费观看| 国产一区二区三区在线臀色熟女| 久久久久久久午夜电影| 日韩三级视频一区二区三区| 琪琪午夜伦伦电影理论片6080| 一个人免费在线观看的高清视频| 色综合站精品国产| 国产精品一区二区三区四区久久 | 国产亚洲av高清不卡| 亚洲色图av天堂| 国产精品国产高清国产av| 国产亚洲av高清不卡| 97人妻天天添夜夜摸| 免费无遮挡裸体视频| 亚洲电影在线观看av| 国产欧美日韩一区二区三区在线| 黑人巨大精品欧美一区二区蜜桃| 亚洲 国产 在线| ponron亚洲| 999久久久国产精品视频| 少妇粗大呻吟视频| 此物有八面人人有两片| 在线视频色国产色| 免费无遮挡裸体视频| 日本五十路高清| 老汉色∧v一级毛片| 精品国产一区二区三区四区第35| www国产在线视频色| 国产精品久久视频播放| 91国产中文字幕| 一级作爱视频免费观看| 性少妇av在线| 精品国产一区二区久久| 免费在线观看亚洲国产| 国产97色在线日韩免费| 在线av久久热| 国产成人啪精品午夜网站| 国内精品久久久久精免费| 国产97色在线日韩免费| 在线观看免费日韩欧美大片| 精品熟女少妇八av免费久了| 日韩有码中文字幕| 国产欧美日韩综合在线一区二区| 免费在线观看黄色视频的| 无人区码免费观看不卡| 国产精品亚洲av一区麻豆| 久久久久精品国产欧美久久久| 国内精品久久久久精免费| 好看av亚洲va欧美ⅴa在| 中文亚洲av片在线观看爽| 97碰自拍视频| 免费看a级黄色片| 一a级毛片在线观看| 99热只有精品国产| 久久久久久久午夜电影| 69av精品久久久久久| 欧美日韩一级在线毛片| 国产成人欧美在线观看| 国产av一区在线观看免费| 麻豆一二三区av精品| 激情在线观看视频在线高清| 在线观看免费视频日本深夜| 久久香蕉精品热| 成人国产一区最新在线观看| 国产真人三级小视频在线观看| 一进一出好大好爽视频| 亚洲中文av在线| 欧美日韩乱码在线| 桃红色精品国产亚洲av| 国产97色在线日韩免费| 免费观看精品视频网站| 亚洲午夜理论影院| 不卡av一区二区三区| 久久国产精品影院| 91字幕亚洲| 日本 av在线| 国产午夜精品久久久久久| 99国产精品一区二区蜜桃av| 欧美乱妇无乱码| 午夜福利成人在线免费观看| www.精华液| 人妻丰满熟妇av一区二区三区| 不卡av一区二区三区| 久久久久久人人人人人| 久久久国产欧美日韩av| 国产免费av片在线观看野外av| 亚洲精品美女久久久久99蜜臀| 狠狠狠狠99中文字幕| 日韩大码丰满熟妇| 亚洲成a人片在线一区二区| 亚洲精品国产色婷婷电影| 久久青草综合色| 国产乱人伦免费视频| 日韩免费av在线播放| 国产免费av片在线观看野外av| 久久久精品欧美日韩精品| 视频在线观看一区二区三区| 女人精品久久久久毛片| 精品国内亚洲2022精品成人| 女性被躁到高潮视频| 丝袜在线中文字幕| 日韩成人在线观看一区二区三区| 波多野结衣av一区二区av| 国产一区二区三区综合在线观看| av网站免费在线观看视频| 午夜精品久久久久久毛片777| 十分钟在线观看高清视频www| 精品久久蜜臀av无| 一级毛片精品| 天天一区二区日本电影三级 | 男女下面插进去视频免费观看| 亚洲精品久久国产高清桃花| 91国产中文字幕| 国产成人av教育| 一级毛片女人18水好多| 国产精品亚洲一级av第二区| 日本五十路高清| 久久中文看片网| 欧美激情高清一区二区三区| 真人一进一出gif抽搐免费| 欧美中文综合在线视频| 91九色精品人成在线观看| 亚洲成人久久性| 久久婷婷人人爽人人干人人爱 | 可以在线观看毛片的网站| 少妇熟女aⅴ在线视频| 亚洲成人久久性| 日韩欧美一区视频在线观看| 亚洲人成77777在线视频| 亚洲精品久久国产高清桃花| 精品久久蜜臀av无| 9191精品国产免费久久| 亚洲中文字幕一区二区三区有码在线看 | 男人操女人黄网站| 中文字幕色久视频| 高潮久久久久久久久久久不卡| 午夜激情av网站| 日韩成人在线观看一区二区三区| 天堂√8在线中文| 老鸭窝网址在线观看| 精品欧美一区二区三区在线| 真人做人爱边吃奶动态| 一区二区三区激情视频| 99国产极品粉嫩在线观看| 亚洲激情在线av| 搡老妇女老女人老熟妇| 午夜精品久久久久久毛片777| 此物有八面人人有两片| 亚洲精品久久成人aⅴ小说| 日韩视频一区二区在线观看| 别揉我奶头~嗯~啊~动态视频| 人人澡人人妻人| 天天一区二区日本电影三级 | 久久久久九九精品影院| 成熟少妇高潮喷水视频| 在线国产一区二区在线| 精品一区二区三区视频在线观看免费| 精品国产美女av久久久久小说| 男女之事视频高清在线观看| aaaaa片日本免费| 国产精品电影一区二区三区| 精品一区二区三区视频在线观看免费| 日韩欧美一区视频在线观看| 国产精品自产拍在线观看55亚洲| tocl精华| 高清在线国产一区| 岛国在线观看网站| 国产三级在线视频| 极品教师在线免费播放| 日本一区二区免费在线视频| 九色国产91popny在线| 女同久久另类99精品国产91| 久热这里只有精品99| 欧美日韩乱码在线| 欧美色欧美亚洲另类二区 | 少妇被粗大的猛进出69影院| 免费搜索国产男女视频| 成人三级做爰电影| 国产精品98久久久久久宅男小说| 国产不卡一卡二| 亚洲九九香蕉| 99久久国产精品久久久| 国产1区2区3区精品| 国产亚洲av高清不卡| 日韩欧美一区二区三区在线观看| 99国产精品一区二区蜜桃av| 国产区一区二久久| 日日摸夜夜添夜夜添小说| 精品少妇一区二区三区视频日本电影| 欧美在线黄色| 性色av乱码一区二区三区2| 国产成人欧美在线观看| 一区二区三区国产精品乱码| 身体一侧抽搐| 午夜免费观看网址| 日韩欧美免费精品| 波多野结衣一区麻豆| 久久午夜综合久久蜜桃| 国产视频一区二区在线看| av在线播放免费不卡| 满18在线观看网站| 国产欧美日韩综合在线一区二区| 我的亚洲天堂| 中文亚洲av片在线观看爽| 亚洲第一青青草原| 精品欧美国产一区二区三| 成人18禁在线播放| 黑人巨大精品欧美一区二区mp4| 中文字幕人成人乱码亚洲影| 淫秽高清视频在线观看| 丝袜人妻中文字幕| 成人亚洲精品av一区二区| 老司机深夜福利视频在线观看| 人人妻,人人澡人人爽秒播| 欧美日本中文国产一区发布| 精品欧美国产一区二区三| 91大片在线观看| 一区二区三区激情视频| 国产私拍福利视频在线观看| 一级片免费观看大全| 91大片在线观看| 欧美亚洲日本最大视频资源| 亚洲国产看品久久| cao死你这个sao货| 高清在线国产一区| 一边摸一边做爽爽视频免费| 美女高潮喷水抽搐中文字幕| 在线观看免费日韩欧美大片| 黄片播放在线免费| 欧美日韩黄片免| 日韩三级视频一区二区三区| 99在线视频只有这里精品首页| 黑人欧美特级aaaaaa片| 精品一区二区三区视频在线观看免费| 男女床上黄色一级片免费看| 国产亚洲av嫩草精品影院| 夜夜躁狠狠躁天天躁| 亚洲 欧美 日韩 在线 免费| 亚洲av美国av| 亚洲一卡2卡3卡4卡5卡精品中文| 成人欧美大片| www.自偷自拍.com| 久久久精品欧美日韩精品| 女人精品久久久久毛片| 老司机深夜福利视频在线观看| 国产精品香港三级国产av潘金莲| cao死你这个sao货| 午夜成年电影在线免费观看| 国产精品一区二区精品视频观看| 在线观看免费视频网站a站| 51午夜福利影视在线观看| 国产亚洲欧美98| 50天的宝宝边吃奶边哭怎么回事| 久久久久久免费高清国产稀缺| 国产午夜精品久久久久久| 免费人成视频x8x8入口观看| 国产成人欧美在线观看| 老司机在亚洲福利影院| 久久人人精品亚洲av| 黄频高清免费视频| 少妇的丰满在线观看| 色婷婷久久久亚洲欧美| 国产不卡一卡二| 免费看十八禁软件| 国产精品一区二区在线不卡| 精品一区二区三区视频在线观看免费| 两个人免费观看高清视频| 亚洲五月婷婷丁香| 最近最新中文字幕大全电影3 | 欧美丝袜亚洲另类 | 99re在线观看精品视频| 国产高清有码在线观看视频 | 少妇 在线观看| 日韩欧美一区二区三区在线观看| 国产精品 国内视频| 日韩欧美在线二视频| 中文字幕精品免费在线观看视频| 久久狼人影院| 嫁个100分男人电影在线观看| 一级毛片女人18水好多| 久久欧美精品欧美久久欧美| 日本a在线网址| 亚洲国产精品久久男人天堂| 欧美日韩精品网址| 欧美黑人精品巨大| 免费看十八禁软件| 99精品在免费线老司机午夜| 久久婷婷成人综合色麻豆| 亚洲精品美女久久av网站| 一进一出抽搐动态| 欧美乱码精品一区二区三区| 国产不卡一卡二| 久久九九热精品免费| 欧美最黄视频在线播放免费| 亚洲精品一区av在线观看| 亚洲五月天丁香| 日本在线视频免费播放| 午夜两性在线视频| 午夜福利,免费看| 99久久国产精品久久久| 欧美乱码精品一区二区三区| 男女做爰动态图高潮gif福利片 | 淫秽高清视频在线观看| 满18在线观看网站| 一级毛片女人18水好多| 日韩欧美国产一区二区入口| 满18在线观看网站| 国产aⅴ精品一区二区三区波| 欧美+亚洲+日韩+国产| 国产成人精品久久二区二区91| 黄片播放在线免费| 久久欧美精品欧美久久欧美| 亚洲国产欧美一区二区综合| 精品久久久久久久久久免费视频| 欧美+亚洲+日韩+国产| 69av精品久久久久久| 久久影院123| 嫩草影视91久久| 一级片免费观看大全| 国产99久久九九免费精品| 久久久国产欧美日韩av| 一边摸一边抽搐一进一小说| 亚洲精品av麻豆狂野| 国产av在哪里看| 男人的好看免费观看在线视频 | 亚洲性夜色夜夜综合| 成人三级做爰电影| 精品一区二区三区视频在线观看免费| 香蕉丝袜av| 熟女少妇亚洲综合色aaa.| 在线观看www视频免费| 两个人视频免费观看高清| 免费在线观看日本一区| 国产xxxxx性猛交| 国产欧美日韩一区二区精品| 在线观看舔阴道视频| 久久精品国产亚洲av香蕉五月| 精品熟女少妇八av免费久了| 欧美性长视频在线观看| 99热只有精品国产| 性欧美人与动物交配| 亚洲中文字幕日韩| 一夜夜www| 国产1区2区3区精品| 日韩 欧美 亚洲 中文字幕| 黄色a级毛片大全视频| 日韩欧美国产在线观看| 国产男靠女视频免费网站| 亚洲成av人片免费观看| 免费无遮挡裸体视频| 欧美日本中文国产一区发布| 欧美另类亚洲清纯唯美| 女人高潮潮喷娇喘18禁视频| 亚洲成av人片免费观看| 老熟妇乱子伦视频在线观看| 国产精品精品国产色婷婷| 国产精品一区二区免费欧美| 中文字幕人成人乱码亚洲影| 免费人成视频x8x8入口观看| 中文字幕高清在线视频| 欧美黑人欧美精品刺激| 自拍欧美九色日韩亚洲蝌蚪91| 12—13女人毛片做爰片一| 夜夜爽天天搞| 91在线观看av| 国产成年人精品一区二区| 欧洲精品卡2卡3卡4卡5卡区| 波多野结衣av一区二区av| 可以在线观看的亚洲视频| 一级a爱视频在线免费观看| 精品免费久久久久久久清纯| 久久青草综合色| 91麻豆精品激情在线观看国产| 亚洲成国产人片在线观看| 欧美亚洲日本最大视频资源| 日韩av在线大香蕉| 精品久久久精品久久久| 99国产精品免费福利视频| 成年女人毛片免费观看观看9| 午夜视频精品福利| 禁无遮挡网站| 久热爱精品视频在线9| 性色av乱码一区二区三区2| 亚洲色图综合在线观看| 久久中文看片网| 日本在线视频免费播放| 亚洲国产毛片av蜜桃av| 99久久国产精品久久久| 日韩中文字幕欧美一区二区| 亚洲情色 制服丝袜| 最好的美女福利视频网| 后天国语完整版免费观看| 欧美不卡视频在线免费观看 | 看黄色毛片网站| 亚洲激情在线av| 一卡2卡三卡四卡精品乱码亚洲| av天堂在线播放| 一级作爱视频免费观看| 亚洲成人免费电影在线观看| 国产成人精品无人区| 亚洲专区字幕在线| 亚洲中文av在线| 看免费av毛片| 在线免费观看的www视频| 精品久久久久久成人av| 色在线成人网| 可以在线观看的亚洲视频| 国产99白浆流出| 亚洲七黄色美女视频| 精品国产乱子伦一区二区三区| 嫩草影视91久久| 欧美性长视频在线观看| 国产一区二区三区视频了| 一级片免费观看大全| 亚洲欧洲精品一区二区精品久久久| 亚洲第一av免费看| 麻豆av在线久日| 亚洲av五月六月丁香网| 欧美激情高清一区二区三区| 精品久久久久久久人妻蜜臀av | avwww免费| 亚洲第一欧美日韩一区二区三区| 国产精品综合久久久久久久免费 | 久久久久九九精品影院| 又黄又爽又免费观看的视频| 人人妻人人爽人人添夜夜欢视频| 校园春色视频在线观看| 侵犯人妻中文字幕一二三四区| 国产成人精品久久二区二区91| 日韩欧美在线二视频| 精品久久久久久,| 久久中文字幕一级| 亚洲精品av麻豆狂野| 欧美亚洲日本最大视频资源| 桃红色精品国产亚洲av| 亚洲成人免费电影在线观看| 一夜夜www| 超碰成人久久| 一夜夜www| 欧美精品啪啪一区二区三区| 欧美日韩乱码在线| 久久人妻av系列| 亚洲全国av大片| 亚洲av第一区精品v没综合| 久久久久久亚洲精品国产蜜桃av| 亚洲专区国产一区二区| 免费高清视频大片| 91大片在线观看| 婷婷精品国产亚洲av在线| 一进一出抽搐动态| 一本大道久久a久久精品| 成人亚洲精品一区在线观看| 欧美丝袜亚洲另类 | 曰老女人黄片| 美女免费视频网站| 日本vs欧美在线观看视频| 亚洲自偷自拍图片 自拍| 国产一区二区三区综合在线观看| 午夜精品在线福利| 每晚都被弄得嗷嗷叫到高潮| 69精品国产乱码久久久| 精品一区二区三区视频在线观看免费| 久久久久九九精品影院| 国产激情欧美一区二区| 变态另类丝袜制服| 91精品国产国语对白视频| 国产一区二区激情短视频| 免费在线观看日本一区| av天堂在线播放| 免费在线观看亚洲国产| 亚洲熟女毛片儿| 麻豆av在线久日| 精品国内亚洲2022精品成人| 国产成年人精品一区二区| 99精品在免费线老司机午夜| 国产欧美日韩一区二区三区在线| 欧美日韩福利视频一区二区| 伊人久久大香线蕉亚洲五| 一本久久中文字幕| 亚洲三区欧美一区| 国产一区二区激情短视频| 免费在线观看日本一区| 99在线人妻在线中文字幕| 首页视频小说图片口味搜索| 亚洲最大成人中文| av视频免费观看在线观看| 国产激情久久老熟女| 每晚都被弄得嗷嗷叫到高潮| 69精品国产乱码久久久| 九色亚洲精品在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品久久久人人做人人爽| 免费搜索国产男女视频| 一二三四在线观看免费中文在| 露出奶头的视频| 正在播放国产对白刺激| 国产高清激情床上av| 1024视频免费在线观看| 无限看片的www在线观看| 国产一区在线观看成人免费| 精品欧美国产一区二区三| 窝窝影院91人妻| 99久久久亚洲精品蜜臀av| 日韩欧美一区二区三区在线观看| 国产高清有码在线观看视频 | 性少妇av在线| 国产av在哪里看| 欧美人与性动交α欧美精品济南到| 999久久久国产精品视频| 日韩高清综合在线| 12—13女人毛片做爰片一| 久久久精品国产亚洲av高清涩受| 乱人伦中国视频| 久久久久国产一级毛片高清牌| 久热这里只有精品99| 淫妇啪啪啪对白视频| 日本五十路高清| 首页视频小说图片口味搜索| 精品国产亚洲在线| 成人手机av|