999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

負斜樁頂部水平受拉響應及p-y 曲線特征試驗研究

2023-01-31 08:12:34周德泉蔡強王創業朱沁周毅李留璽
湖南大學學報(自然科學版) 2023年1期
關鍵詞:承載力水平模型

周德泉 ,蔡強,王創業,朱沁,周毅,李留璽

(長沙理工大學 土木工程學院,湖南 長沙 410114)

傾斜樁應用于輸電線塔、海洋平臺、海上風機等承受水平荷載的構筑物樁基礎中,但其工作機制與計算方法尚不十分清楚.學術界和工程界在豎直樁水平工作特性研究的基礎上開展了傾斜樁水平工作特性及其p-y曲線特征研究,揭示傾斜樁水平受力響應特征,推動傾斜樁的應用.

豎直樁水平工作特性研究方面,崔新壯等[1]基于粉質黏土中剛性豎直樁水平受荷模型試驗,揭示了剛性直樁的失穩機理,認為土體的軟化是樁失穩破壞的根本原因.趙明華等[2]提出了適合斜坡地基直樁前土體的兩種破壞模式,推導出兩種剛性直樁的極限承載力計算公式,認為內摩擦角及黏聚力增加會提高極限承載力.曹兆虎等[3]設計了一套基于透明土料和PIV 技術的水平加載及光學測量系統來觀測樁、土變形特性,分析發現受荷一側土體的流動呈一定孔擴張的趨勢,隨著深度增加,擴張趨勢逐漸減小至0,樁體在埋深78%位置發生旋轉剛性破壞.張浦陽等[4]、王建華等[5]通過試驗研究發現嵌巖樁水平承載力隨嵌巖深度增加而增加,嵌巖深度超過某一定值后不再發生明顯的變化.周天應等[6]基于理論計算,分析了嵌巖樁的水平承載力,認為在水平荷載下,彈性長樁與較長的嵌巖樁的受力變形特性相似.趙學亮等[7]通過室內模型試驗研究吸力式三筒基礎水平受力特性,發現隨著筒間距的增加,水平承載力有一定提升,水平承載力隨長徑比增大而增大.鄒新軍等[8]通過室內模型試驗與ABAQUS 研究了中砂與軟黏土中豎直樁樁頂豎向荷載對水平承載力的影響,發現中砂中豎向荷載會提高樁的水平承載力,而軟黏土中豎向荷載會減小樁的水平承載力.

傾斜樁水平工作特性研究方面,Zhang 等[9]、Meyerhof 等[10]、Manoppo 等[11]、Hazzar 等[12]、朱 海濤[13]、周德泉等[14]研究發現樁身傾斜角對傾斜樁水平承載力有一定影響,相對于直樁,正斜樁的水平承載力有一定提高,而負斜樁的水平承載力有一定降低.朱照清等[15]在東海附近海域開展大直徑鋼管正、負斜樁與直樁的水平受荷現場試驗,發現正斜樁的承載力最高.樊文甫等[16]通過數值模擬研究帶承臺傾斜單樁水平受荷下的水平承載力,發現正斜樁最大,直樁居中,負斜樁最小,且長徑比對斜樁水平位移影響不大.顧明等[17]開展了水平偏心荷載下斜群樁受力響應的離心機模型試驗,分析了水平荷載偏心距對群樁水平承載力的影響以及斜樁群樁與直樁群樁的差異,發現斜樁群樁優于直樁群樁.

p-y曲線研究方面,龔健等[18]在軟土中開展微型單樁及群樁水平受荷試驗,發現單樁中斜樁比直樁有更好的水平承載力,并應用p-y曲線法對單樁與群樁的水平位移進行計算,計算值能較好地反映實際情況.曹衛平等[19]為得到水平受荷樁p-y曲線,開展室內模型試驗,建立試驗砂土應力應變與斜樁的p-y曲線無量綱化后兩者之間的關系(p/pub=(y/y50)/(a1+b1y/y50)),建立斜樁的p-y曲線并求解,與實測值相差不大.楊曉峰等[20]從非線性位移假設、雙曲線應力-應變關系和樁側剪應力計算三個方面對應變楔模型進行修改,并通過實例驗證了修正模型的合理性.Kim 等[21]通過模型試驗探索樁的安裝方法和樁頭的約束方式對樁p-y曲線的影響,對地基極限土反力和初始地基彈性模量進行模擬,提出新的p-y曲線公式(p=y/(1/kini+y/pub)).Ashour 等[22]基于水平荷載下豎直樁的應變楔模型推求斜樁水平荷載下應變楔模型,該模型考慮了土、樁、傾斜角度、傾斜方向對p-y曲線的影響.楊明輝等[23]開展黏性土中平地和不同邊坡條件下基樁水平加載室內模型試驗,基于樁在邊坡的位置與邊坡角度修正土的初始地基模型和極限土抗力,建立了考慮邊坡處土抗力折減效應的水平受荷樁的p-y曲線計算公式.

綜上所述,已有學者對豎直樁[1-8]、傾斜樁[9-17]的水平工作性狀及其p-y曲線[18-23]進行了研究,而樁頂水平受拉狀態下豎直樁、傾斜樁及嵌巖傾斜樁受力響應差異的對比試驗少見開展,制約了傾斜樁工作機制的研究與計算方法的提出.

本文通過室內模型試驗,在模型砂土中設置0°、10°、20°樁及20°樁底嵌巖樁,樁頂水平受拉,拉力方向與傾斜方向相反,測試不同拉力作用下4 根負斜樁的水平位移、應變片和土壓力盒應變,以期獲得負斜樁樁身水平位移、彎矩、土抗力和p-y曲線隨傾斜角和樁底約束條件的變化規律,促進負斜樁的推廣應用,對于輸電線塔、海洋平臺、海上風機等承受水平荷載的構筑物樁基礎的設計和安全評價也具有重要意義.

1 室內模型試驗

模型試驗在1 420 mm×750 mm×1 200 mm(長×寬×高)模型箱(專利號:ZL201520323607.3)中開展.模型箱主骨架由高強度角鋼焊接成框,箱壁用鋼化玻璃或木板封閉.

模型樁(以下簡稱“樁”)的布置如圖1 所示.模型箱用3 個隔板在長度方向隔成4 個相同隔室(長×寬×高為355 mm×750 mm×1 200 mm).4 根樁分別固定在4 個隔室中線.由樁頂確定樁的位置(離模型箱壁200 mm),模型箱內固定角度從左到右分別為0°、10°、20°、20°.特別注意第4 根20°樁是嵌巖樁,嵌巖方法是:將20°樁插入開口木盒內事先攪拌的混凝土中,混凝土通過鋼釘與地板錨固.當混凝土凝結硬化后,樁與混凝土墩凝結成為一體.

圖1 試驗布置Fig.1 Test layout

模型樁由水泥砂漿與鋼筋籠在模具中澆筑、振搗、養護而成.鋼筋籠截面呈正方形,邊長為20 mm,長800 mm,縱向鋼筋采用4 根直徑為6 mm 的HRB235 鋼筋,布置在正方形截面的四角,箍圈采用直徑2.2 mm的14#鐵絲制作.樁身對稱兩側面用光滑木條預留一個800 mm×30 mm×5 mm(長×寬×深)的凹槽.模型樁總共4 根,每根樁長800 mm,截面呈正方形,邊長40 mm,具體參數見表1.樁身對稱兩側面凹槽中對稱布置應變片,并用環氧樹脂充滿凹槽.應變片總共8 對,第1 對應變片離樁頂70 mm,然后每隔100 mm 對稱布置1 對應變片.封槽后,在樁身凹槽面對稱布置5 對土壓力盒,與樁頂的距離分別為110 mm、260 mm、410 mm、560 mm、710 mm,并在每個土壓力盒下10 mm 處通過環氧樹脂固定一根長度為40 mm、直徑φ20 mm 的PVC 軟管(共5 根,為百分表測試水平位移預留孔洞).

表1 模型樁參數Tab.1 Parameters of the model piles

模型砂土用紗網過篩晾干而成.土體密度1.83 g/cm3,土粒相對密度2.67,最大粒徑5 mm,含水率1%,不均勻系數Cu=5.5,曲率系數Cc=2.7,級配良好.級配曲線見圖2.砂土填滿模型箱前,將4 根樁按預定的位置與角度布置在4 個隔室中線.先固定樁頂,再調整樁身傾斜角0°、10°或者20°,最后固定樁底.樁身出土長度皆為40 mm,方便對樁施加水平荷載(水平荷載施加在離砂土面20 mm 高度處).填砂采用“砂雨法”,每填埋20 cm 后平整砂面、再進行下一次填埋,完成6次填埋后,厚度共計1 200 mm.

圖2 模型砂集配曲線Fig.2 Grading curve of model sand

填土前,將直徑φ15 mm的PVC硬管穿過模型箱壁的預留小孔,與樁側已經粘貼好的直徑φ20 mm 的PVC 軟管連接,用透明膠密封,另一端伸出模型箱外50 mm.填土后,靜置一段時間.測試前,將百分表的加長探針穿過PVC硬管與塑料軟管的中心并接觸樁側表面,要求不接觸管壁,然后固定百分表.每根樁安裝6個百分表,其中1個百分表安裝在樁頂.

本次試驗所用加載裝置由鋼絞線、定滑輪、掛鉤、砝碼組成.鋼絞線的一端通過繩卡和AB 膠固定在樁頂,另一端通過模型箱壁頂部定滑輪、連接掛鉤.在掛鉤上分級增加標準砝碼,實現樁頂水平分級受拉.特別注意:必須保證鋼絞線與樁身截面中線平行、定滑輪固定不擺動,因為傾斜樁對樁頂拉力非常敏感.最小、最大荷載分別為40 N、440 N,中間每級荷載增加50 N.試驗采用TDS-540 應變儀及其電腦軟件測量應變與土壓力.加載前,應變片和土壓力盒應變調零,記錄百分表初讀數.加載后立即采集1 次應變、記錄百分表讀數.以后每隔10 min 讀取百分表、采集應變.當百分表2 次讀數計算的水平位移差小于0.01 mm時,進行下一級加載.

2 試驗結果與分析

需要說明,外荷載直接通過樁傳遞給巖土時形成主動樁(見圖3(a)(b)),通過巖土間接傳遞給樁時形成被動樁(見圖3(c)(d)).為表達方便,外力促使傾斜樁傾斜角減小時稱為負斜樁(圖3(a)(c)),外力促使傾斜樁傾斜角增大時稱為正斜樁(圖3(b)(d)).本文模型樁的頂部直接承受拉力后傾斜角減小,屬于負斜主動樁(圖3(a)).

圖3 正斜樁與負斜樁示意圖Fig.3 Schematic diagram of positive and negative batter pile

2.1 樁頂水平位移規律

圖4 反映了Z0、Z1、Z2、Z3 樁頂端水平位移變化規律.

圖4 負斜主動樁樁頂水平位移曲線對比Fig.4 Comparison curves of horizontal displacement of negative batter active piles

分析發現:樁頂水平位移加載曲線呈下凹形,曲線光滑,樁頂水平位移y隨水平拉力F(以下簡稱荷載或者水平荷載)的增大逐漸增大,增速隨水平拉力F的增大逐漸增大.相同水平拉力F作用下,樁頂水平位移由大到小的排列順序為Z2、Z1、Z0、Z3,說明樁頂受拉時,負斜樁樁頂水平位移隨樁身傾斜角的增加而逐漸增大,而樁底嵌巖會大大減小樁頂的水平位移.換言之,樁頂受拉時,負斜樁的水平承載力隨樁身傾斜角的增加而逐漸減小,而樁底嵌巖會大大增加負斜樁的水平承載力.此規律與砂土中完成的傾斜樁離心模型試驗結果[9]相似.說明本次水平位移試驗數據可靠.

2.2 樁身整體水平位移規律

圖5 所示為F=290 N、440 N 時4 根負斜主動樁樁身水平位移y隨入土深度z變化規律.本文約定樁身水平位移方向與樁頂拉力方向相同為正,反之為負.經對比分析發現:

圖5 負斜主動樁水平位移曲線對比Fig.5 Horizontal displacement contrast curves of negative batter active piles

1)負斜主動樁樁身水平位移隨荷載的增大而增大,隨著入土深度增加,Z0、Z1、Z2 樁身水平位移先由最大值減小為0,后再逐漸負向增大,即此時3 根樁繞樁身下部某位置發生了旋轉,這一位置為樁“不動點”(水平位移為零,本文稱“位移零點”),Z0、Z1、Z2 樁“位移零點”入土深度分別在590 mm、580 mm、550 mm 附近.Z3 樁水平位移隨著樁入土深度的增加逐漸減小,樁頂最大,樁底最小,繞樁底旋轉彎曲.說明“位移零點”位置隨傾斜角的減小或者樁底約束程度的提高而下移.

2)本模型負斜主動樁的位移模式為“轉動+彎曲”.4 樁的水平位移曲線都有一定彎曲.在“位移零點”上下,上部彎曲會大于樁下部的彎曲變形.Z3 樁彎曲變形很明顯比Z0、Z1、Z2樁都要大,這是因為Z3樁樁底嵌巖,底部很難發生位移,從而使得Z3 樁的彎曲變形比其他3根非嵌巖樁要大.

3)負斜主動樁水平位移隨傾斜角增大而增大,隨樁底嵌固程度提高而減小,也就是說,傾斜角的增加會降低樁的水平承載能力,而樁底嵌巖會提高樁的水平承載能力.隨著荷載的增大,Z2樁整體樁身水平位移在逐漸遠離其他3 根樁,說明Z2 樁樁身水平位移增大的速率在4 根樁中最大.分析認為:樁頂水平拉力分解為樁身軸向壓力及其法向力,該樁身法向力頂托樁身上部土體,使負斜樁樁身正面方向的土體變得松弛而隆起,土體對樁身的約束能力隨之減弱,且隨著傾斜角增大,負斜樁樁身上部松弛土體范圍也越來越大,所以,負斜樁傾斜角越大,水平位移就越大.

2.3 樁身彎矩規律

負斜樁樁身彎矩根據公式(1)進行計算.

樁頂水平受拉290 N、440 N 時,Z0、Z1、Z2、Z3 樁身彎矩M-入土深度z曲線對比圖(簡稱M-z曲線)如圖6所示.

圖6 負斜主動樁彎矩對比曲線Fig.6 Bending moment comparison curves of negative batter active piles

分析發現:

1)Z0、Z1、Z2、Z3 樁M-z曲線呈“弓”形,隨入土深度的增加,樁身彎矩先增加到峰值,后逐漸減小.Z1 樁的峰值位置(z約為350 mm)比Z0 樁(z約為290 mm)明顯降低,Z1 樁和Z2 樁的峰值位置接近,Z3 樁彎矩峰值點位置(z約為400 mm)比Z2 樁(z約為350 mm)明顯降低.由此推理,負斜主動樁在頂部受拉后將發生彎曲破壞,破壞位置在樁身中部偏上處,且隨傾斜角增大或者樁底約束程度提高而降低.

2)Z0、Z1、Z2 樁彎矩峰值從大到小依次是Z2、Z1、Z0,彎矩峰值隨傾斜角增大而增大.Z3與Z2 彎矩曲線相交于Z2 樁峰值點下一定距離.交點上部,Z2樁彎矩值大于Z3 樁彎矩值;交點下部,Z3 樁彎矩值大于Z2 樁彎矩值.Z3 樁彎矩峰值小于Z2 樁彎矩峰值.由此推理,負斜主動樁頂部水平受拉時,樁底嵌巖降低樁身彎矩峰值與中上部彎矩值、增加整體抗彎能力;增大傾斜角,彎矩峰值隨之增大,大傾斜角負斜主動樁更容易發生彎曲破壞.

2.4 樁身土抗力規律

負斜樁頂部水平受拉時,樁擠壓土,土體對樁產生約束,樁身受到土抗力.在樁兩側面z=75 mm、225 mm、375 mm、525 mm、675 mm處對稱布置土壓力盒,采用TDS-540 應變采集儀測量土壓力盒的微應變,后期采用每個土壓力盒對應的標定方程計算每級荷載的土抗力值.4 根負斜主動樁頂部受拉290 N、440 N時土抗力p-入土深度z規律曲線(簡稱p-z曲線)如圖7所示.p-z曲線光滑、形態穩定,具有明顯的規律,說明數據可靠.

圖7 負斜主動樁土抗力曲線對比Fig.7 Soil resistance comparison curves of negative batter active piles

分析發現:

1)負斜主動樁頂部受拉時,p-z曲線呈“弓”形,有峰值.土抗力峰值隨負斜樁傾斜角增大而減小,峰值位置逐漸下降.土抗力從大到小依次是Z0、Z1、Z2樁,Z0、Z1、Z2 樁峰值點出現的位置分別是z=230 mm、250 mm、310 mm 附近.在每級荷載下,Z0、Z1、Z2、Z3 樁土抗力隨著入土深度增加,土抗力先增加到峰值后再逐漸減小,并且Z0、Z1、Z2 樁在樁底附近出現負值土抗力.這是因為樁繞樁身下部發生了旋轉,樁底向負方向移動,樁底正面遠離樁身正面土體,土體松散.Z0、Z1、Z2 樁土抗力存在土抗力“零點”,都在z=600 mm 附近,且與樁身位移“零點”位置相近.土抗力“零點”位置與位移“零點”位置相近,再次說明測試結果可靠.

引入法。借鑒國外成熟的烹飪技術,有效地對中國菜進行改進,打造出中西合璧的菜品,比如蝦仁食料,將其創新出中菜西做或者西菜中做的美味佳肴。

2)負斜樁頂部受拉時,樁底嵌巖會降低樁身上部土抗力及其峰值和位置、增加樁身下部土抗力.Z2樁土抗力峰值位置上部土抗力大于Z3 樁,峰值位置下部土抗力小于Z3樁,這是因為Z2樁的位移“零點”(z=550 mm)比Z3 樁(z=800 mm)位置較高,使Z2 樁位移“零點”附近土抗力變小.Z2 的土抗力峰值比Z3大,峰值點出現的位置比Z3 高.比較發現,Z3 與Z2的p-z曲線與M-z曲線具有相似的特征.

2.5 實測p-y曲線特征

負斜樁頂部受拉時,負斜樁繞樁身下部第4 測點附近轉動、樁底嵌固樁底部位移為零,故取每級荷載下轉動點上部3 個測點(樁身上部3 個土壓力盒)的水平位移y與土抗力p,作出Z0、Z1、Z2、Z3 樁3 個測點處p-y曲線,見圖8.分析發現:

圖8 負斜主動樁p-y曲線對比Fig.8 Comparison of p-y curves of negative batter active piles

1)負斜主動樁p-y曲線特征與位置有關.位置越高,p-y曲線越容易收斂;位置越低,p-y曲線越不容易收斂.說明負斜主動樁存在淺層土與深層土的劃分深度或者極限水平土抗力轉折點的深度.

2)負斜主動樁p-y曲線特征與傾斜角有關.Z0、Z1、Z2、Z3 負斜樁在相同位置的p-y曲線線形相似,但變化幅度不一樣.擬合后更能定量說明.

2.6 p-y曲線擬合特征

p-y曲線反映了樁在水平荷載作用下產生水平位移時樁周土對樁產生的土抗力,它受到樁身傾斜角、樁徑、樁土剛度比、樁頭約束條件等多因素的影響.目前已有多種擬合曲線方程,本文參考已有p-y曲線方程,寫出本文擬合曲線方程,見式(2).

式中:pub為樁周土極限土抗力,kPa;y50為土抗力等于pub/2時對應的樁身水平位移,mm.

采用雙曲線擬合實測p-y曲線,根據Duncan-Chang 模型獲得各測點的極限土抗力pub,從而得到水平位移值y50,其數值見表2.

表2 負斜主動樁pub、y50值隨測點位置的變化Tab.2 Changes of pub and y50 of negative batter active piles with different measuring points

將樁身3 個測點處各級荷載下的土抗力p除以pub,水平位移y除以y50,用各測點的數據采用公式(2)與最小二乘法對3 個測點數據進行擬合,選取最合適的曲線方程,最終得p-y曲線擬合方程,α與β值見表3,擬合曲線見圖9.

表3 負斜主動樁p-y擬合曲線的α與β值變化Tab.3 Variation of α and β in p-y fitting curves of negative batter active piles

圖9 負斜主動樁p-y曲線擬合對比Fig.9 Comparison of p-y curves fitting of negative batter active piles batter active pile

分析發現:

1)采用p/pub=α(y/y50)β擬合p-y曲線,擬合效果較好.

2)參數α在0.45~0.51 之間變化,參數β在0.27~0.62之間變化.Z0、Z1、Z2、Z3樁參數α變化幅度分別為0.03、0.03、0.06、0.04,參數β變化幅度分別為0.21、0.20、0.24、0.14,說明傾斜角與樁端約束程度對α影響不大,傾斜角對β值變化幅度影響較小,α與β值變化幅度體現p-y曲線隨深度變化的敏感度,兩者敏感度隨傾斜角度減小、樁端約束程度的增加而減小.

3 結 論

1)樁頂水平受拉時,本模型負斜主動樁繞樁身某點轉動并發生彎曲,位移模式為“轉動+彎曲”,轉動點位置隨傾斜角的減小或樁底約束程度的提高而下移.樁頂水平位移隨樁身傾斜角的增加或樁底約束程度的降低而增大,也就是說,增加傾斜角或降低樁底約束程度將減小負斜主動樁的水平承載力.

2)樁頂水平受拉時,砂土中負斜主動樁彎矩、土抗力沿深度變化曲線呈“弓”形.增大傾斜角,彎矩峰值逐漸增大,土抗力峰值卻逐漸減?。粯兜浊稁r,彎矩峰值和土抗力峰值均降低,樁身中上部彎矩、土抗力降低,樁身下部彎矩、土抗力增加,增大整體抗彎能力.

3)負斜主動樁p-y曲線可采用p/pub=α(y/y50)β擬合,參數α在0.45~0.51 之間變化,參數β在0.27~0.62 之間變化.傾斜角與樁端約束程度對α影響不大,傾斜角對β值變化幅度影響較小,β值變化幅度隨樁端約束程度提高而降低,可為p-y曲線法求解負斜主動樁內力提供參考.

猜你喜歡
承載力水平模型
一半模型
張水平作品
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
加強上下聯動 提升人大履職水平
人大建設(2019年12期)2019-05-21 02:55:32
3D打印中的模型分割與打包
CFRP-PCP板加固混凝土梁的抗彎承載力研究
耐火鋼圓鋼管混凝土柱耐火極限和承載力
潛艇極限承載力計算與分析
對受壓加勁板極限承載力計算方法的評述
主站蜘蛛池模板: 日韩无码真实干出血视频| 91久久夜色精品| 性欧美精品xxxx| 亚洲中文字幕av无码区| 亚洲免费毛片| 国产欧美在线| 亚洲天堂久久| 国产精品综合久久久| 国产真实自在自线免费精品| 国产一级精品毛片基地| 欧美一级黄色影院| 国产区人妖精品人妖精品视频| 欧美一级在线看| 国产精品性| 99精品免费欧美成人小视频| 国产黄在线观看| 无码专区在线观看| 国产一区二区三区精品欧美日韩| 亚洲男人的天堂在线| 久久精品中文无码资源站| 欧美黑人欧美精品刺激| 97久久人人超碰国产精品| 国产人前露出系列视频| 拍国产真实乱人偷精品| 美女被操91视频| 91免费片| 亚洲国产成人久久77| 污网站免费在线观看| 国产精品视频免费网站| 日本一区高清| 九九九精品成人免费视频7| 少妇被粗大的猛烈进出免费视频| 最新国产高清在线| 国产美女在线观看| 色婷婷综合激情视频免费看| 一本视频精品中文字幕| 尤物成AV人片在线观看| 午夜福利在线观看入口| 欧美福利在线观看| 国产精品网址你懂的| 成人在线综合| 人妻精品久久无码区| 国产又爽又黄无遮挡免费观看| 亚洲成人在线免费| 国产熟睡乱子伦视频网站| 欧美一级视频免费| 狠狠综合久久久久综| 国产精品刺激对白在线| 国产亚洲第一页| 日本一区二区三区精品视频| 麻豆精品在线播放| 在线观看精品自拍视频| 亚洲 日韩 激情 无码 中出| 午夜三级在线| 国产男女免费视频| 午夜a视频| 国产成人av大片在线播放| 波多野结衣久久高清免费| 国产丝袜第一页| 91久久大香线蕉| 国产视频欧美| 欧美日韩专区| 亚洲成人免费看| 国产不卡国语在线| 亚洲免费福利视频| 亚洲国产高清精品线久久| 91视频青青草| 成人av专区精品无码国产| 国产精品无码AV中文| 久久综合伊人 六十路| 美女内射视频WWW网站午夜| 玖玖精品在线| 国产真实乱人视频| 一区二区在线视频免费观看| 欧美日韩免费观看| 在线不卡免费视频| 国产大片喷水在线在线视频| 国产乱子伦视频在线播放 | 97一区二区在线播放| 国产成人无码综合亚洲日韩不卡| 亚洲成A人V欧美综合天堂| 国产小视频在线高清播放|