亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        膨脹土地基中單樁受扭非線性分析

        2020-11-14 06:40:56王順葦歐孝奪
        工程力學(xué) 2020年11期
        關(guān)鍵詞:界面模型

        江 杰,王順葦,歐孝奪,王 智,楊 迪

        (1. 廣西大學(xué)土木建筑工程學(xué)院,廣西,南寧 530004;2. 工程防災(zāi)與結(jié)構(gòu)安全教育部重點(diǎn)實驗室,廣西,南寧 530004;3. 廣西防災(zāi)減災(zāi)與工程安全重點(diǎn)實驗室,廣西,南寧 530004;4. 中鐵建設(shè)集團(tuán)有限公司,北京 100040)

        膨脹土在全球40 多個國家均有分布,范圍十分廣泛,在我國主要分布在云南、廣西、四川等10 多個省份[1]。其脹縮特性是造成地基上方建(構(gòu))筑物發(fā)生傾斜、變形、不均勻沉降的的重要因素,具有十分突出的工程問題[2]。為了保證膨脹土地區(qū)建(構(gòu))筑物的安全,實際工程中往往采用樁基礎(chǔ)形式。大型廣告牌、風(fēng)力發(fā)電塔以及彎梁橋在不均勻水平力的作用下,會對其樁基礎(chǔ)產(chǎn)生扭矩[3 ? 7],當(dāng)豎向荷載較小時,扭矩將會成為樁身承載力的控制性因素,在設(shè)計中若忽略扭矩的作用,可能會造成重大的工程損失[8]。膨脹土浸水后會對樁基礎(chǔ)產(chǎn)生軸向拉力作用[1, 9],而研究表明軸向力的存在會使單樁扭矩承載力降低[10 ? 12],因此有必要對膨脹土地基中的受扭單樁進(jìn)行進(jìn)一步的研究,為工程設(shè)計提供參考。

        國內(nèi)外學(xué)者對膨脹土地區(qū)豎向和水平荷載作用下單樁的承載特性進(jìn)行了相關(guān)研究[13 ? 19],對于受扭單樁,Stoll[20]在砂土中對鋼管樁回填混凝土首次對單樁進(jìn)行了扭矩加載試驗,得到了樁頂扭矩和扭轉(zhuǎn)角的關(guān)系曲線,為后續(xù)研究奠定了基礎(chǔ);Poulos[21]基于彈性連續(xù)介質(zhì)理論求解出了樁頂扭矩和扭轉(zhuǎn)角的關(guān)系,并在粘性土中進(jìn)行了驗證;Zhang 等[22]利用力學(xué)的變分原理推導(dǎo)出了樁的控制微分方程,通過迭代法求解出了樁頂扭矩和扭轉(zhuǎn)角的關(guān)系,并在砂土和黏土地基中進(jìn)行了驗證;鄒新軍等[4]結(jié)合平衡原理和剪切位移法求得了砂土中單樁受扭的彈塑性解;趙明華等[23]假設(shè)樁為彈性梁,采用非線性彈簧模擬樁土間的相互作用關(guān)系,建立了粉土地基中受扭單樁的簡化計算模型;Basack 等[24]采用邊界元法求得砂土地基中單樁在扭矩作用下的數(shù)值解;Li 等[25]利用現(xiàn)有的界面剪切試驗建立荷載傳遞模型,采用有限差分法,得到單樁受扭的非線性計算方法。

        綜上所述:目前國內(nèi)外學(xué)者主要針對砂土、黏土和粉土地基中的受扭單樁進(jìn)行了研究,對于膨脹土地區(qū)受扭單樁的研究較少。膨脹土浸水后樁身軸力的存在,使得常規(guī)的單樁受扭計算理論會過高的估計承載力,偏于不安全。鑒于此,本文對膨脹土浸水后樁側(cè)豎向摩阻力的分布模式進(jìn)行了探討,在此基礎(chǔ)上求出膨脹土地基浸水后受扭單樁的非線性解,通過膨脹土地基中單樁模型試驗,對理論進(jìn)行了驗證,并對浸水前后單樁扭轉(zhuǎn)承載特性以及不同膨脹率的膨脹土浸水對基樁的受力和變形的影響進(jìn)行了分析。

        1 計算模型及假定

        本文的分析基于以下假定:

        1)樁為線彈性等截面樁。

        2)定義被浸潤的膨脹土層厚度為“膨脹影響深度”;假設(shè)膨脹影響深度內(nèi)的膨脹土隆起量隨深度線性變化[15],膨脹影響深度以下的隆起量為0,計算公式如下:

        式中:ws(z)為深度z處的膨脹土隆起量;h0為膨脹影響深度;ws0為膨脹的土地表膨脹量。

        3)膨脹土浸水會對樁身產(chǎn)生上拔力,在樁身中性點(diǎn)(樁土相對位移為0 處)以上對樁身作用豎直向上的摩阻力,在中性點(diǎn)下對樁身作用豎直向下的摩阻力。因此可假設(shè)膨脹土浸水時的樁-土界面模型如圖1 所示,浸水過程中樁-土界面豎向正、負(fù)摩阻力以及扭矩作用下的樁身環(huán)向摩阻力(s(z)≥0段)均符合Kraft 等[26]提出的荷載傳遞函數(shù)形式。

        圖1 膨脹土浸水條件下樁-土界面模型Fig. 1 Load-transfer model for pile-soil interaction under expansive soil swelling

        式 中:s(z)為 樁-土 相 對 位 移s(z) =w(z)?ws(z);ws(z)為樁所處位置深度z處的膨脹土浸水后的隆起量;Gs為土體小應(yīng)變時的初始剪切模量;ψ=τ(z)Rf/τf,Rf為應(yīng)力-應(yīng)變曲線擬合常數(shù),可取0.9~1.0;τf為樁側(cè)極限摩阻力;r0為樁截面半徑;rm為有效影響半徑。

        2 理論分析及推導(dǎo)

        2.1 環(huán)向極限摩阻力的確定

        浸水后對樁頂施加扭矩作用,樁側(cè)土體對樁身產(chǎn)生環(huán)向的摩阻力,樁側(cè)豎向及環(huán)向摩阻力共同作用于樁身。如圖2 所示,采用邊界元法將樁縱向離散為n個長度為L/n的單元,第i單元的樁側(cè)豎向和環(huán)向摩阻力分別為τv(i)和τt(i),假設(shè)兩者共同作用下的樁側(cè)合摩阻力τ(i)不大于極限摩阻力τf(i),可得到:

        圖2 浸水后單樁邊界元離散示意圖Fig. 2 Boundary element discretization of pile after immersion

        式 中:τf(i) 為 第i單 元 的 極 限 摩 阻力;τtf(i)為 第i單元環(huán)向極限摩阻力。

        由式(4)可知要想得到環(huán)向的極限摩阻力,須先求出單樁浸水后的豎向摩阻力。

        2.2 浸水作用下單樁豎向摩阻力 τv(i)的求解

        式(14)的求解過程如下:① 假設(shè)沿樁身節(jié)點(diǎn)位移w=0,由樁側(cè)土體自由位移場ws(z)求樁土相對位移s(z)=w(z)?ws(z),進(jìn)而求出樁身節(jié)點(diǎn)切線剛度;② 由式(14)求出沿樁身節(jié)點(diǎn)的位移wk;③ 用新求出的樁身節(jié)點(diǎn)位移wk求出樁-土相對位移和樁側(cè)土體的切線剛度。利用新的割線剛度求出新的樁身節(jié)點(diǎn)位移wk+1;④ 取|wk+1?wk|作為迭代控制誤差,若誤差大于限定值則重復(fù)②~④直至迭代誤差小于限定值。

        由以上數(shù)值解求出樁土相對位移,結(jié)合式(2)即可求出樁身單元豎向摩阻力。

        2.3 浸水作用下單樁受扭非線性計算

        為得到每個樁單元的樁側(cè)環(huán)向摩阻力,采用中心差分法將樁長離散成n個相等的單元,并在樁頂和樁端各增加一個虛擬等分節(jié)點(diǎn),并代入樁端和樁頂邊界條件對式(18)進(jìn)行差分離散,可以得到方程組:

        式(24)的求解過程如下:① 假設(shè)沿樁身節(jié)點(diǎn)的扭轉(zhuǎn)角θ 為任意非零矩陣,由式(17)求得深度z處的相對扭轉(zhuǎn)位移st(z) ,將st(z)代入式(2)求出τt(z) ,由 式(4)求 出τtf(z) ,若τt(z)≥τtf(z),則 取τt(z)=τtf(z) ,由式(9)求出kθ(z) ,進(jìn)而得到Kt′;② 由式(24)求出沿樁身節(jié)點(diǎn)的扭轉(zhuǎn)角θk;③ 用新求出的樁身節(jié)點(diǎn)扭轉(zhuǎn)角θk求出樁身扭轉(zhuǎn)剛度矩陣Kt′k;利用新求出的樁身扭轉(zhuǎn)剛度矩陣Kt′k求出新的樁身節(jié)點(diǎn)扭轉(zhuǎn)角θk+1;④ 取|θk+1?θk|作為迭代控制誤差,若誤差大于限定值則重復(fù)②~④直至迭代誤差小于限定值。

        3 方法驗證與分析

        通過試驗分別對樁側(cè)豎向摩阻力τv(i)和環(huán)向摩阻力τt(i)的求解進(jìn)行驗證,在此基礎(chǔ)上將浸水條件下單樁受扭的求解與本文的試驗結(jié)果與進(jìn)行對比分析。

        3.1 試驗及測量裝置

        采用自主研制的樁基加載裝置,對膨脹土地基中的單樁進(jìn)行浸水前后的室內(nèi)靜載模型試驗,加載裝置如圖3 所示。試驗裝置主要分為三個部分:模型箱、加載裝置和浸水裝置。模型箱尺寸為1 m×1 m×1 m。加載裝置由加載架、滑輪和力臂組成。浸水裝置由豎向浸水管和水桶組成,模型箱底部填筑約250 mm 厚的砂礫層,浸水管插入砂礫層。通過水管將模型箱底部和水桶連接,利用虹吸原理對水桶注水,以保持水平面與土表面高程一致,實現(xiàn)對膨脹土的全浸水膨脹。

        圖3 加載裝置示意圖Fig. 3 Schematic of loading device

        在膨脹土地基表面設(shè)置百分表測量地表隆起量,通過固定于樁頂?shù)牧拷瞧鳒y量樁頂扭轉(zhuǎn)角,在樁頂設(shè)置一個長為0.1 m 的力臂,通過定滑輪對力臂施加水平力,使其轉(zhuǎn)化為扭矩,兩個定滑輪的摩擦系數(shù)μ1和μ2分別為0.943 和0.914,樁頂扭矩可通過下式確定:

        式中:G1和G2分別為兩端重物的重量; θt為測得的樁頂扭轉(zhuǎn)角,為方便計量,本試驗采用水來代替重物。

        通過應(yīng)變片測量樁身切應(yīng)變和軸向應(yīng)變,樁身扭矩和軸力分別由式(25)、式(26)確定:

        式中:樁身剪切模量Gp=Ep/(2(1+ν));樁身泊松比ν取值為0.3; γ為樁身切應(yīng)變; ρ為應(yīng)變片距中性軸的距離;Ep為樁身彈性模量; ε為樁身軸向應(yīng)變。

        3.2 模型樁

        采用鋁合金空心管制作單樁模型,彈性模量69.7 GPa,長800 mm,埋入膨脹土650 mm,直徑25 mm,壁厚2 mm。在模型樁內(nèi)表面每隔100 mm設(shè)置一組測點(diǎn),每組測點(diǎn)設(shè)置一個BF350–3AA和BHF350-3HA 型應(yīng)變片分別測量樁身軸向應(yīng)變和切應(yīng)變,最后一組應(yīng)變片設(shè)置在樁端以上5 mm處,共設(shè)7 組測點(diǎn),通過TST3822EN 靜態(tài)應(yīng)變測試儀采集應(yīng)變。樁端采用尼龍塞進(jìn)行粘貼封底,待模型箱內(nèi)土體填置100 mm 之后將模型樁垂直插入土體,土體填筑過程中不斷測量模型樁與模型箱之間的距離,對樁實時進(jìn)行糾偏,保持其豎直。

        3.3 土體參數(shù)

        試驗用膨脹土取自廣西南寧水牛研究所附近的灰白色膨脹土,其自由膨脹率為65.5%,屬于中等膨脹土,基本物理力學(xué)性質(zhì)如表1 所示。土料經(jīng)過烘干、粉碎,制成初始含水量為20%的土樣。土樣采用人工夯實進(jìn)行分層填筑,每層控制夯實后的厚度為100 mm,填土總高度為700 mm,填筑土層干密度為1.46 g/cm3。

        表1 膨脹土基本物理力學(xué)性質(zhì)Table 1 Physical and mechanical properties of expansive soil

        3.4 試驗方案

        為了研究膨脹土浸水對單樁扭轉(zhuǎn)承載特性的影響,共設(shè)計3 組試驗,第1 組測量浸水后的樁身軸向應(yīng)變;第2、3 組分別測量浸水前后,不同的樁頂扭矩作用下的樁身切應(yīng)變和樁頂扭轉(zhuǎn)角,具體如表2 所示。

        表2 試驗加載方案及測量數(shù)據(jù)Table 2 Test loading scheme and measurement data

        試驗采用分級的方式加載,每級加載增量為預(yù)計極限荷載的1/10,當(dāng)樁頂扭矩-扭轉(zhuǎn)角曲線(T~θ 曲線)出現(xiàn)突變時加載結(jié)束。對于浸水前的模型樁,夯土結(jié)束后即可進(jìn)行加載試驗,對于浸水后的模型樁,浸水完成后即可進(jìn)行加載試驗。當(dāng)百分表讀數(shù)與應(yīng)變儀數(shù)據(jù)趨于穩(wěn)定后48 h內(nèi)再無明顯波動,則浸水試驗完成。本次試驗浸水歷時14 d,地表隆起0.011 m。試驗裝置如圖4所示。

        圖4 模型樁加載試驗裝置Fig. 4 Loading test device of model pile

        3.5 試驗結(jié)果與分析

        1)浸水后樁側(cè)豎向摩阻力分布

        通過第1 組試驗測得的樁身軸向應(yīng)變求出樁身軸力,利用式(5)對樁身軸力數(shù)據(jù)進(jìn)行處理得到樁身豎向摩阻力,圖5 為浸水完成后樁身豎向摩阻力分布圖,圖中豎向摩阻力向上為正,由圖5可知,樁側(cè)豎向摩阻力呈現(xiàn)出先增大再減小然后反向增大的特點(diǎn)。大約在樁深為0.4 m 以上部分的膨脹土對樁產(chǎn)生向上的上拔力,摩阻力為正值;0.4 m以下部分的膨脹土對樁作用向下的拉力,阻礙單樁的抬升,摩阻力為負(fù)值。即未施加樁頂扭矩時,中性點(diǎn)位置約在樁的中下部0.4 m 處,與張大峰[1]的結(jié)論相似,計算結(jié)果和試驗數(shù)據(jù)相比最大誤差為2.8%相差較小,驗證了本文豎向摩阻力的計算方法。

        圖5 樁身豎向摩阻力分布Fig. 5 Vertical frictional resistance distribution of pile shaft

        2)浸水后樁側(cè)環(huán)向摩阻力分布

        對樁頂施加扭矩之后,樁-土界面之間會產(chǎn)生環(huán)向相對位移,從而產(chǎn)生了環(huán)向摩阻力,為驗證本文環(huán)向摩阻力計算的正確性,分別提取第2 組試驗中樁頂扭矩為1.5 N·m、3.0 N·m、4.5 N·m、6.0 N·m 所對應(yīng)的樁身切應(yīng)變數(shù)據(jù),并求出樁身扭矩,利用式(15)對樁身扭矩數(shù)據(jù)進(jìn)行處理得到樁身環(huán)向摩阻力,計算結(jié)果與試驗數(shù)據(jù)對比如圖6所示,4 組數(shù)據(jù)與試驗結(jié)果均比較接近,證明了環(huán)向摩阻力計算的正確性。

        圖6 浸水后樁側(cè)環(huán)向摩阻力分布Fig. 6 Distribution of lateral circular friction after immersion

        由圖6 可知:在0 m~0.15 m 范圍,由于樁土豎向相對位移較大,樁-土界面已經(jīng)進(jìn)入塑性階段,導(dǎo)致樁側(cè)環(huán)向摩阻力較??;在0.15 m~0.25 m深度范圍,樁-土界面處于彈塑性狀態(tài),樁側(cè)環(huán)向摩阻力快速增加;在0.25 m 深度以下,當(dāng)樁頂扭矩為1.5 N·m 和3.0 N·m 時,由于樁-土界面相對位移較小尚處于處于彈性狀態(tài),環(huán)向摩阻力幾乎沒有變化,當(dāng)樁頂扭矩增加到4.5 N·m 和6 N·m時,樁-土界面漸漸進(jìn)入塑性階段,由于土體的抗力作用越來越明顯,環(huán)向摩阻力也在增加,呈現(xiàn)出非線性的變化。

        3)浸水前后荷載位移對比分析

        對于浸水前單樁,只需令豎向摩阻力為0 ,引入Poulos[21]提出的樁端邊界條件進(jìn)行計算,即可得到T-θ 曲線。將本文計算得到的T-θ 曲線與試驗數(shù)據(jù)進(jìn)行對比,如圖7 所示,當(dāng)樁頂扭矩較大時本文計算方法和試驗結(jié)果均較為吻合,但當(dāng)荷載較小時誤差較大,原因是試驗過程中當(dāng)樁頂扭矩較小時,樁頂扭轉(zhuǎn)角并沒有明顯的變化,導(dǎo)致觀測結(jié)果存在一定的誤差。但整體來看本文計算得到的T-θ 曲線和試驗數(shù)據(jù)較為吻合,證明本文提出的膨脹土地基浸水后單樁受扭非線性計算方法的正確性。此外還可看出,浸水前后單樁的極限扭矩分別為19.82 N·m 和6.88 N·m,即浸水后單樁的極限扭矩降低了65%。

        圖7 浸水前后T~θ 曲線對比Fig. 7 Comparison of T~θ curves before and after immersion

        分別提取第2、3 組試驗中樁頂扭矩為2 N·m和4 N·m 時樁身扭矩試驗數(shù)據(jù)與計算數(shù)據(jù)進(jìn)行對比,如圖8 所示:浸水前后的樁身扭矩分布計算結(jié)果與試驗數(shù)據(jù)相比趨勢均一致且誤差較小,浸水前樁身的扭矩近似線性分布,與文獻(xiàn)[4]趨勢一致,再次證明了本文計算方法的正確性。浸水后單樁的樁身扭矩分布呈現(xiàn)出非線性變化的規(guī)律,在深度為0 m~0.25 m 范圍,樁身扭矩變化幅度較小,是由于在施加扭矩之前,此深度范圍內(nèi)膨脹土的膨脹作用導(dǎo)致樁土豎向剪切位移較大,樁-土界面已經(jīng)發(fā)生了滑移(簡稱“滑移段”),從而使環(huán)向摩阻力的作用不明顯,導(dǎo)致樁身扭矩的增量較?。辉?.25 m 深度以下,樁身扭矩呈現(xiàn)出近似線性快速減小的特點(diǎn),究其原因,樁土豎向相對位移較小,樁-土界面尚未發(fā)生滑移(簡稱“未滑移段”),環(huán)向摩阻力有較大的發(fā)揮空間,土體作用的環(huán)向摩阻力使樁身扭矩快速減小。此外還可以看出當(dāng)樁頂扭矩相同時,浸水前后的樁身扭矩分布有著較大的差異,在實際工程中可根據(jù)樁身扭矩分布特點(diǎn)進(jìn)行合理配筋。

        圖8 浸水前后樁身扭矩分布情況Fig. 8 Distribution of pile torque before and after immersion

        4)膨脹率的影響

        為研究不同膨脹率的膨脹土浸水后對單樁T~θ 曲線的影響,設(shè)不同膨脹率的膨脹土對應(yīng)的地表隆起高度為0.004 m~0.016 m,由圖9 可知:隨著膨脹率的增加,單樁的極限扭矩依次降低,且降低的幅度越來越小;當(dāng)樁頂扭矩為定值時,隨著膨脹土膨脹率的增加,樁頂扭轉(zhuǎn)角增加的幅度越來越大,樁-土體系的加載剛度降低的幅度越來越大。因此,地表隆起會造成樁-土體系加載剛度和極限扭矩大幅度降低,且膨脹土的膨脹率越高,單樁抗扭能力和樁-土體系加載剛度降低的幅度越大。

        圖9 不同膨脹率的膨脹土對應(yīng)的T~θ 曲線Fig. 9 T~θ curve corresponding to expansive soil with different expansion rates

        圖10 為樁頂扭矩為4.5 N·m 時不同膨脹率的膨脹土對應(yīng)的樁身扭矩分布情況,假設(shè)不同膨脹率的膨脹土對應(yīng)的地表隆起高度同上,由圖10 可知,膨脹率對樁身扭矩的分布有著重要的影響,隨著膨脹率的增加“滑移段”的長度也在增加。

        圖10 不同膨脹率的膨脹土對應(yīng)的樁身扭矩Fig. 10 Pile shaft torque corresponding to expansive soil withdifferent expansion rates

        4 結(jié)論

        為探討膨脹土地基中受扭單樁的承載特性,考慮膨脹土浸水隆起對樁側(cè)摩阻力及樁端邊界條件的影響,采用有限差分法和邊界元法,對樁身的內(nèi)力進(jìn)行了分析,主要結(jié)論如下:

        (1)提出了膨脹土地基浸水后單樁樁側(cè)豎向摩阻力的計算方法,在此基礎(chǔ)上,結(jié)合邊界元法提出了膨脹土地基浸水后單樁受扭的非線性分析方法,并通過單樁模型試驗進(jìn)行驗證,結(jié)果較為吻合。

        (2)樁側(cè)豎向摩阻力沿樁身先增大后減小,隨后反向增大,中性點(diǎn)位于樁身的中下部;樁側(cè)環(huán)向摩阻力在樁身中上部接近于0,然后沿樁身快速增加,達(dá)到峰值之后,當(dāng)樁-土界面處于彈性狀態(tài)時趨于穩(wěn)定,處于塑性狀態(tài)時沿樁身緩慢減小。

        (3)常規(guī)的計算方法高估了基樁的極限扭矩和樁-土體系加載剛度,對于浸水后的膨脹土地基并不適用,實際工程中應(yīng)考慮地表隆起造成的單樁極限扭矩和樁-土體系加載剛度的降低。

        (4)在樁頂扭矩的作用下,浸水前后單樁的樁身扭矩分布規(guī)律有著很大的差異,浸水前單樁的樁身扭矩沿樁身近似線性減小,浸水后在“滑移段”樁身扭矩變化較小,“未滑移段”轉(zhuǎn)身扭矩近似線性減小??筛鶕?jù)此特點(diǎn)對膨脹土地區(qū)受扭單樁進(jìn)行合理配筋。

        (5)膨脹土的膨脹率會對單樁承載力和樁身扭矩產(chǎn)生較大的影響,隨著膨脹率的增加,單樁的抗扭能力明顯降低,“滑移段”的長度也在增加。

        (6)本文采用有限單元法對膨脹土地基中的受扭單樁進(jìn)行了分析,由于暫無實測資料,采用室內(nèi)模型試驗進(jìn)行了驗證,并對均質(zhì)地層中樁身受力情況進(jìn)行了分析。

        猜你喜歡
        界面模型
        一半模型
        重要模型『一線三等角』
        國企黨委前置研究的“四個界面”
        重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
        基于FANUC PICTURE的虛擬軸坐標(biāo)顯示界面開發(fā)方法研究
        空間界面
        金秋(2017年4期)2017-06-07 08:22:16
        電子顯微打開材料界面世界之門
        人機(jī)交互界面發(fā)展趨勢研究
        3D打印中的模型分割與打包
        FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
        人妻无码第一区二区三区| 日韩中文字幕无码av| 给我播放的视频在线观看| 久久久精品人妻一区二区三区日本 | 精品国产一区二区三区av 性色| 1精品啪国产在线观看免费牛牛| 国产午夜精品av一区二区三| 国产极品少妇一区二区| 狠狠色噜噜狠狠狠狠7777米奇| 四虎在线播放免费永久视频| 美女被插到高潮嗷嗷叫| 男人的天堂一区二av| 午夜福利院电影| 久久精品免费无码区| 亚洲国产精品成人一区二区三区 | 欧美自拍视频在线| 国产白浆精品一区二区三区| 三个黑人插一个女的视频| 亚洲开心婷婷中文字幕| 日本不卡在线视频二区三区| 国产一区二区三区免费在线视频| 熟女免费观看一区二区| 中国午夜伦理片| 国产午夜无码视频免费网站| 日本精品国产1区2区3区| 亚洲成人免费av影院| 特级av毛片免费观看| 精品一区二区三区久久久| 少妇特殊按摩高潮不断| 日韩人妻精品中文字幕专区| 无码国模国产在线观看| 色伊人国产高清在线| 一区二区三区高清视频在线| 久久精品国产只有精品96| 风流少妇又紧又爽又丰满| 亚洲老女人区一区二视频| 五月激情综合婷婷六月久久| 射死你天天日| 超清无码AV丝袜片在线观看| 亚洲成年国产一区二区| av无码人妻中文字幕|