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

NC-UHPC界面抗剪性能試驗及基于內聚力模型的仿真分析

2025-08-04 00:00:00黃維王婷張雨婷斯比努爾·下木希奴周知
華僑大學學報(自然科學版) 2025年4期
關鍵詞:鋼纖維抗剪齡期

中圖分類號:U454 文獻標志碼:A 文章編號: 1000-5013(2025)04-0400-1C

NC-UHPC Interface Shear Performance Test and Simulation Based on Cohesive Zone Model

HUANG Wei1'2,WANG Ting1'2, ZHANG Yuting 1,2 ,SIBINUER·Xiamuxinu1'2,ZHOU Zhi3

(1.Hubei Key Laboratory of Theory and Application of Advanced Materials Mechanics, Wuhan University of Technology,Wuhan 430070,China;2. Department of Mechanics and Engineering Sciences,Wuhan University of Technology,Wuhan 43007o,China;3.School of Transportation and Logistics Engineering,Wuhan Universityof Technology,Wuhan 430o7o,China)

Abstract:To investigate the shear mechanism at the interface between ultra-high performance concrete (UHPC) and normal concrete (NC),nine groups of double-sided shear tests were conducted on NC-UHPC specimens,the effects of steel fiber content,concrete age,and boundary constraint conditions on interfacial mechanical properties were studied. Combining with ABAQUS software,a refined finite element model considering interfacial bond-slipcharacteristics was established,and an interfacial bond stress distribution model was constructed through theoretical derivation. The results show that under unrestrained conditions,specimens with low steel fiber content (0%-1% volume fraction) exhibit single-surface failure characteristics,while fixing the botom of theNC blocks efectively activatesa dual-surface failure mode.The steel fiber content in UHPC positively correlates with the interfacial shear strength,while the interfacial stiffness is lessaffected by the steel fiber content.The increase of UHPC age enhances interfacial stifess,while higher NC age has the inhibitory effect.The proposed finiteelement model effectively simulates the damage evolutionatthe interface, but simulates the asymmetric strain field in double-shear specimens with limitations.Based on the principle of energy conservation,the two-parameter model of interfacial load-displacement is established,and its prediction results conforms well with the test data.

Keywords:NC-UHPC interface;pushout test;digital image correlation;finite element;cohesive zone model

由于普通混凝土(NC)與超高性能混凝土(UHPC)界面強度的影響因素眾多且復雜,從上個世紀至今,國內外眾多學者進行了大量的研究,但至今仍存在可進一步深究的部分,如界面處理方式、界面粘結方式的參數影響規律、界面粘結機理等方面。

國內外學者通過大量實驗研究新、舊混凝土界面的粘結性能。郭進軍[1采用Z字形粘結抗剪試驗,建立抗剪強度與界面粗糙度及混凝土強度的關系式。葉果[2研究抗剪性能與界面粗糙度及鋼筋摩擦力的關系。張雷順等[3]通過界面抗剪試驗,得出植筋率、混凝土抗壓強度與粘結強度的關系式。刑強[4]發現界面強度與混凝土強度等級、連接界面面積和鋼筋面積密切相關,并提出抗剪強度公式。Tayeh 等[5]通過劈裂和斜剪試驗,揭示NC-UHPC界面粘結的不同破壞階段。Semendary等[6]則通過斜剪、直接拉伸和剪切試驗,探討混凝土強度、表面粗糙度等因素對粘結性能的影響。Al-Madani等7研究基體處理方式和養護條件對界面粘結的具體影響。Bastien-Masse 等[8」研究纖維分布對界面性能的影響,而Hussein等[9發現鋼纖維摻量對粘結劈拉強度影響不大。

同時,有限元軟件在提升新、舊混凝土界面上也起到了重要作用。黨祺[10使用非線性有限元軟件Marc 模擬了UHPC層與NC 層的混合結構,通過Links-Insert、Contact、Touching 等功能增強了兩層結構的連接和鋼筋結合,提高了整體強度和穩定性。王興旺[11]利用ANSYS研究切槽、鑿毛、植筋等技術對混凝土抗剪性能的影響,發現鋼筋的應用比其他方法更能有效提升承載力。田穩苓等[12]通過剪切實驗和數值模型分析新、舊混凝土的抗拉強度,結果顯示植筋率的增加和三角形排列能提高粘結力和抗剪強度。Harris 等[13]通過劈拉、斜剪和拉拔試驗分析混凝土界面粗糙度對粘結性能的影響,并開發了有限元模型評估UHPC 的應力分布。Dias-da-Costa 等[14]利用DIANA 軟件研究混凝土界面的縱向抗剪強度,通過A型直剪試件和平面應力雙線性有限元模擬,結合Mohr-Coulomb 摩擦定律,深人分析各參數對抗剪性能的影響。Farzad等[15]開發了一種高效的數值模擬方法,預測不同混凝土層間的抗剪強度,并通過調整參數來提升粘結性能。

雖然國內外學者對NC-UHPC界面粘結性能的影響因素與機理開展了大量探索,但仍存在以下關鍵問題亟待解決:界面抗剪強度對鋼纖維摻量與混凝土養護齡期的耦合響應規律尚未量化,現有數值模型因忽略界面粗糙度隨機性而難以精確模擬粘結-滑移過程,且缺乏對界面荷載傳遞動態演化機制的理論解析。針對上述問題,本文通過多變量控制試驗量化了鋼纖維摻量(體積分數 φ 為 0%~2% 與混凝土養護齡期 (7~28d) 對抗剪強度的協同效應;基于內聚力模型(CZM)構建考慮界面形貌隨機特征的精細化有限元模型,模擬了界面損傷萌生、擴展至失效的全過程;最后,結合試驗與模擬結果,提出了界面荷載傳遞的三階段發展理論,揭示了界面應力重分布與能量耗散的動態規律。

1試驗方案

1.1 試驗材料及配合比

C4O 混凝土和超高性能混凝土每立方米材料配合比,分別如表1、2所示。表1、2中: m 為質量。

(單位: kg )

表1C40 混凝土每立方米材料配合比Tab.1Material mix proportion per cubic meter of C4O concrete (unit:kg
表2超高性能混凝土每立方米材料配合比 (單位:kg) Tab.2Material mix proportion per cubic meter of UHPC (unit:kg)

1.2 試件設計

考慮到NC與UHPC的粘結長度、UHPC的鋼纖維摻量(體積分數)、NC與UHPC的齡期及混凝土強度等級對界面粘結強度的影響,設計并制作了9組 NC-UHPC 試件,其中,NC-UHPC 試件由2個150mm×150mm×50mm 的NC試件和1個 150mm×150mm×50mm 的UHPC試件組成,試件尺寸,如圖1所示。

Fig.1 Specimen dimensions (unit:mm)

9組試件對應的工況,如表3所示。表3中: tn 為 NC 齡期; tu 為UHPC齡期。試件U1O-N17-0、U10-N17-1和U10-N17-2在上、下端分別粘結了鋼板固定,其加載方式與其他試件有所不同。

表39組試件對應的工況Tab.3Working conditions corresponding to 9 groups of specimens

1.3 測量方案

采取雙面直剪試驗測定試件的粘結強度。為了準確得到NC-UHPC界面的滑移情況,在界面附近的 NC 上端粘結鋼片,布置位移計,獲得 NC 滑移情況。在界面UHPC 側粘結應變片,以測量界面應變情況,加載裝置采用電子萬能試驗機。采用數字圖像相關(DIC)技術設備對試件U10-N17-0、U10-N17-1和U10-N17-2的界面受力過程進行全程測量。

1.4 試驗現象

由于各試件破壞過程相似,以試件U10-N17-0為例,對加載破壞過程進行說明,其荷載 (F)- 位移 (Δ) 曲線,如圖2所示。試件破壞現象,如圖3所示。

由圖2、3可知:試驗加載初期(點 A ),NC-UHPC界面無明顯裂紋,表明此時NC-UHPC構件界面還未出現滑移或滑移量很小,其界面應變最大值出現在NC-UHPC左側界面的下側(圖3(a));隨著荷載增加至點 B ,NC-UHPC界面仍然無明顯裂紋,應變在NC-UHPC兩側界面下側均達到較大值(圖3(b));隨著荷載的進一步加大,在左側出現微裂縫并迅速增大,但荷載達到峰值點 C 時,右側界面也出現微小裂縫(圖3(c));隨著左側界面開裂,右側界面也迅速開裂,最后達到點 D (圖3(d))。

圖2試件U10-N17-0 的荷載-位移曲線Fig.2Load-displacement curve ofspecimen U10-N17-0
圖3試件的破壞現象Fig.3Phenomenon of specimen failure

1.5 荷載-位移曲線

9組試件的荷載-位移曲線,如圖4所示。各試件的特征參數,如表4所示。表4中: KG 為抗剪剛度; Fn 為峰值荷載; τp 為抗剪強度; Sn 為峰值荷載對應的滑移量。

圖4不同試件的各工況荷載-位移曲線Fig.4Load-displacement curves of different specimens under various working conditions

由圖4和表4可知:當鋼纖維摻量從 0% 增至 2% (試件U7-N14-0、U7-N14-1和U7-N14-2)時,界面峰值荷載與抗剪強度分別提升 32.9% 和 33.3% ,峰值滑移量增加 0.37mm ,但界面剛度下降約9.6% ,表明鋼纖維的增韌效應在提升界面承載能力的同時會削弱其剛性;在NC齡期均為28d條件下,UHPC齡期由7d延長至28d(試件U7-N28-0、U14-N28-0 和 U28-N28-0),可使界面剛度提升61.6% ,峰值荷載與抗剪強度分別增大 364.2% 和 363.4% ,且峰值滑移量增加 1.26mm (試件U28-N28-0較試件U7-N28-O),表明UHPC水化進程對界面力學性能具有主導作用;雙側約束試件U10N17-0因鋼板限制側向位移,其峰值荷載較無約束試件U7-N14-0提高了 6.2% ,但鋼纖維摻量為 1% 的試件U10-N17-1 出現剛度與抗剪強度異常,剛度較同組試件U10-N17-0、U10-N17-2 分別提高了 8.8% 和 8.6% ,抗剪強度分別提高了 37.0% 和 6.1% ,推測與約束條件下纖維分布均勻性改善有關。

表4各試件的特征參數Tab.4 Characteristic parameters of each specimen

1.6UHPC-NC界面抗剪承載力影響因素分析

1.6.1鋼纖維摻量不同約束條件下的抗剪強度-鋼纖維摻量曲線,如圖5所示。由圖5可知:無約束條件下,抗剪強度隨著鋼纖維摻量增加呈線性增長趨勢(圖5(a), R2=0.998? );而雙側約束試件中,鋼纖維摻量為 1% 試件的抗剪強度異常高于鋼纖維摻量為 2% 的試件(圖5(b))。這一現象與試件澆筑工藝密切相關,采用先澆筑NC基體、后澆筑UHPC的分層施工方法,可能導致鋼纖維在界面區域的分布受限(尤其高鋼纖維摻量時纖維易在NC表面堆積),致使纖維橋聯效應未能充分發揮,從而削弱了界面抗剪強度的提升效率。

圖5不同約束條件下的試件抗剪強度-鋼纖維摻量曲線Fig.5Shear strength-steel fiber content curve of specimens under different constraint conditions

1.6.2混凝土齡期進一步分析UHPC與NC齡期對抗剪強度的影響,結果如圖6所示。由圖6可知:UHPC齡期從7d增至28d時,抗剪強度增長呈現顯著線性特征;盡管線性擬合相關系數較低(R2=0.880 ,但齡期 14~28d 的抗剪強度增幅達 11.16MPa ,是 7~14d 增幅 (1.16MPa 的9.62倍,表明UHPC后期水化反應加速了界面過渡區(ITZ)的密實化進程,進而強化了界面力學性能。相比之下,NC 齡期固定為28d時,其抗剪強度發展對界面抗剪性能的貢獻趨于穩定。

圖6不同試件的抗剪強度-混凝土齡期曲線Fig.6Shear strength-concrete age curve of different specimens

2 NC-UHPC界面抗剪性能有限元分析

2.1材料及界面本構關系

UHPC和NC采用ABAQUS軟件中的混凝土損傷塑性(CDP)模型,其壓縮拉伸特性采用國家規范GB50010一2010《混凝土結構設計規范》中提出的單軸受拉和單軸受壓本構,混凝土本構模型應力 (σ)- 應變(e)曲線,如圖7所示。圖7中:fc,r 為抗壓強度; εc,r 為與峰值壓應力相對應的應變; εt,r 為與峰值拉應力相對應的應變; ft,r 為抗拉強度。

2.2 單元選取

圖7不同試件的混凝土的應力-應變曲線Fig.7Stress-strain curves of concretewith different specimens

NC-UHPC構件主體部分為NC和UHPC試塊,因此,在建模時所選擇的單元類型將直接影響計算的準確性和速度,NC和UHPC 單元均選用8節點6面體線性縮減積分單元(C3D8R)進行計算。

NC-UHPC界面粘結強度在實際工程中起重要作用。為準確模擬單推荷載下NC與UHPC界面的粘結-滑移特性,需要在NC與UHPC之間建立相應的特殊連接單元來實現粘結-滑移效果。選擇Cohesive單元進行分析,本構模型為雙線性模型,為剪切牽引力與裂紋滑動位移的關系,如圖8所示。圖8中: τ 為剪切牽引力; δ 為裂紋滑動位移; 分別為峰值切應力及其對應的滑動位移。

2.3模型建立與網格劃分

圖8雙線性本構模型

建立幾何模型(圖9(a)),并在試塊底部施加固定約束,在頂部Fig.8Bilinear constitutive model施加位移載荷(圖9(b))。采取結構化網格劃分,為了方便計算,網格不宜太密,綜合多次試驗結果,網格尺寸選擇 10mm (圖9(c))。試驗采用內聚力接觸的方法定義接觸。定義內聚力接觸需要的主要參數包括材料剛度 (K) 和損傷變量,均可通過試驗數據計算獲得。

圖9 有限元模型

2.4有限元計算結果與試驗結果的對比分析

圖10不同試件的有限元計算結果與試驗結果的對比

從有限元計算結果中提取荷載-位移曲線,與試件U7-N14-0、U7-N14-1、U7-N14-2 和試件U10-N17-0、U10-N17-1、U10-N17-2 的試驗結果進行對比,如圖10所示。由圖10 可知:采用內聚力模型建立的 NC-UHPC粘結界面有限元模型能較好地模擬其界面的受力發展過程,其計算得到的承載力誤差在 5% 以內。

試件U10-N17-0、U10-N17-1和U10-N17-2峰值荷載對應應變分布的有限元計算結果與試驗結果的對比,如圖11所示。

圖11不同試件的峰值荷載對應應變分布的有限元計算結果與試驗結果的對比Fig.11 Comparison between finite element calculation results and test results of strain distribution under peak load of different specimens

由圖11可知:采用文中方法得到的峰值荷載有限元計算結果比試驗結果小;有限元計算結果應變分布對稱,而試驗結果應變偏向一側。產生這種現象的原因是建立的有限元模型為對稱模型,而試驗因為界面存在初始差異,造成界面一側達到最大,而另一側處于相對較小的應變狀態。這也說明采用內聚

力模型建立的有限元模型能較好地反映界面整體的受力過程,但無法精確考慮因實際缺陷產生的偏差。

3 NC-UHPC界面抗剪受力分析

3.1NC-UHPC界面抗剪粘結應力分布

為了分析普通混凝土與超高性能混凝土之間的相對滑移,參考殷小激[16的理論分析結果,對新、舊混凝土界面粘結應力進行推導。根據文獻[17]研究表明,UHPC-NC界面粘結應力 τ 在試驗荷載上升階段沿截面位置 d 呈現負指數分布,其表達形式為

τ(d)=Ae-Bd°

粘結界面任意截面位置 d 處的相對滑移 S(d) 的計算式為

由最小勢能原理和平衡方程可以得

式 (1)~(4) 中: A,B 為待定系數; P 為外力荷載; L 為粘結長度; Ec 為 NC 的彈性模量; Ac 為 NC 的橫截面積; Cs 為粘結面寬度的2倍;系數 AC,其中,E。為 UHPC 的彈性模量,A,為UH-PC 的橫截面積。

可通過加載端位移計測得,外力荷載 P 可通過試驗測得,則通過式(3)、(4)可計算得 A 和 B ,代入式(1)可得粘結應力分布曲線。

3.2NC-UHPC界面粘結應力計算結果與試驗結果對比

計算得到試件U10-N17-0、U10-N17-1和U10-N17-2在0.2、0.4、0.6、0.8、1.0倍峰值荷載時的界面粘結應力分布,如圖12所示。圖12中: τm 為界面粘結應力。由圖12可知:理論結果與有限元計算結果基本吻合,界面粘結應力沿粘結長度呈現負指數分布;隨著外荷載的增大,加載端部的應力增加幅值明顯。

將不同試件DIC測量得到的0.6倍、0.8倍和1.0倍峰值荷載時對應的界面粘結應力分布試驗結果與理論結果、有限元計算結果進行對比,如圖13所示。

由圖13可知:DIC測量得到的界面粘結應力分布試驗結果與理論結果、有限元計算結果變化趨勢相同,但在數值上存在一些差異。產生這種現象的原因主要是試件界面不均勻,試件界面粗糙及微裂縫發展使DIC 測量存在一定的誤差。

圖13界面粘結應力分布的試驗結果、理論結果和有限元計算結果的對比 Fig.l3Comparison of test results,theoretical results and finite element calculation results for interface bonding stress distribution

4結論

針對新、舊混凝土的抗剪性能進行試驗和有限元數值模擬,設計9組超高性能混凝土和普通混凝土組合構件的界面剪切試驗,分析UHPC中鋼纖維摻量和混凝土養護齡期對界面抗剪強度的影響。

1)通過9組NC-UHPC試件界面剪切試驗,研究了超高性能混凝土的鋼纖維摻量、混凝土齡期及邊界約束條件對界面抗剪性能的影響。結果表明:在無邊界約束工況下,當UHPC 鋼纖維摻量為0%~1% 時,試件易發生單面破壞;當NC試塊下端固定時,試件更易出現雙剪破壞;UHPC 鋼纖維摻量與界面抗剪強度呈正相關,而界面剛度受鋼纖維摻量影響較小;界面剛度隨著UHPC 齡期的增加而增加,但隨著NC齡期的增加而減小。

2)在ABAQUS有限元軟件中采用內聚力模型,建立基于界面性能的NC-UHPC試塊精細有限元模型,將有限元計算模型和試驗得到的荷載-位移曲線進行對比,結果表明,試驗結果與有限元計算結算結果相似,峰值荷載相近,誤差較小。將峰值荷載對應的應變云圖與DIC測量得到的試驗結果進行對比發現,因雙剪試驗兩界面存在一定差異,其應變發展并不對稱。采用內聚力模型計算得到的界面荷載-位移曲線從宏觀上反映了界面的發展,但是無法精確模擬界面因構造差異產生的受力不均勻現象。

3)以NC-UHPC試件的單壓試驗為基礎,基于能量分析的方法,運用ABAQUS有限元軟件中的內聚力模型的計算方法,得到彈性范圍適用的NC-UHPC界面荷載-位移曲線沿粘結界面豎向長度的分布及函數表達式。與內聚力模型計算結果和試驗結果進行對比,證明了理論公式的適用性。

參考文獻:

[1]郭進軍.高溫后新老混凝土粘結的力學性能研究[D].大連:大連理工大學,2003.

[2]葉果.新老混凝土界面抗剪性能研究[D].重慶:重慶大學,2011.

[3]張雷順,王二花,閆國新.植筋法新老混凝土粘結面剪切性能試驗研究[J].鄭州大學學報(工學版),2006,27(3):34- 37.DOI:10.3969/j.issn. 1671-6833.2006.03.008.

[4]邢強.新舊混凝土界面的連接方法及受力性能研究[D].西安:西安科技大學,2012.

[5]TAYEH B A,BAKAR B H A,JOHARI M A M,et al.Mechanical and permeability properties of the interface between normal concrete substrate and ultra high performance fiber concrete overlay[J].Construction and Building Materials,2012,36(6):538-548.DO1:10.1016/j.conbuildmat. 2012.06.013.

[6]SEMENDARY A A,HAMID W,KHOURY I,et al.Experimental investigation of direct tension bond performance of high-strength concrete and ultrahigh-performanceconcreteconnectionsJ].Journal of Materials in Civil Engineer ing,2019,31(9):04019171.DOI:10.1061/(ASCE)MT.1943-5533.0002800.

[7]AL-MADANI M K,AL-OSTA M A,AHMAD S,et al.Interfacial bond behavior between ultra high performance concrete and normal concrete substrates[J].Construction and Building Materials,202,32o(7):126229.DOI:10. 1016/j.conbuildmat.2021.126229.

[8]BASTIEN-MASSE M,DENARIEE,BRUHWILER E.Effect of fiber orientation on the in-plane tensile response of UHPFRC reinforcement layers[J]. Cement and Concrete Composites,2016,67(2):11-125.DOI:10.1016/j.cemconcomp.2016.01.001.

[9]HUSSEIN L,AMLEHL.Structural behavior of ultra-high performance fiber reinforced concrete-normal strength concrete or high strength concrete composite members[J].Construction and Building Materials,20l5,93(3):105- 1116.DO1:10.1016/j. conbuildmat. 2015.05.030.

[10]黨祺.UHPC 加固損傷板受彎性能的研究[D].長沙:湖南大學,2016.

[11] 王興旺.UHPC 與普通鋼筋混凝土結構界面抗剪性能研究[D].長沙:湖南大學,2016.

[12] 田穩苓,宋昭,肖成志,等.帶植筋新老混凝土粘結面剪切試驗及有限元模擬分析[J].建筑結構,2019,49(23):127- 133.DOI:10.19701/j. jzjg. 2019.23.024.

[13]HARRIS D K,CARBONELL MUNOZ M A,GHEITASI A,et al.The challenges related to interface bond characterization of ultra-high-performance concrete with implications for bridge rehabilitation practices[J].Advances in Civil Engineering Materials,2015,4(2) :75-101. DOI:10. 1520/ACEM20140034.

[14]DIAS-DA-COSTA D,ALFAIATE J,JULIO E.FE modeling of the interfacial behaviour of composite concrete members[J].Construction and Building Materials,201,26(1):233-243.DOI:10.1016/j.conbuildmat.2011.06. 015.

[15]FARZAD M,SHAFIEIFAR M,AZIZINAMINI A. Experimental and numerical study on bond strength between conventional concrete and ultra high-performance concrete (UHPC)[J].Enginering Structures,20l9,186:297- 305.DOI:10.1016/j. engstruct. 2019.02.030.

[16]殷小激.型鋼混凝土粘結-滑移推出試驗的理論分析[J].結構工程師,2010,26(3):48-53.DOI:10.3969/j.issn. 1005-0159.2010.03.008.

[17]胡海昌.彈性力學的變分原理及其應用[M].北京:科學出版社,1981.

猜你喜歡
鋼纖維抗剪齡期
養護齡期對赤泥基泡沫輕質土性能的影響研究
鋼橋面UHPC基層與SMA面層層間施工處理方法研究
栓釘連接雙鋼板-再生混凝土組合剪力墻抗震性能試驗研究
鋼筋混凝土開洞梁抗剪承載力統一模型及其可靠度分析
灌注樁混凝土硬化熱力響應特性現場試驗研究
基于內聚力模型的配筋鋼纖維混凝土結構承載性能研究
基于內聚力模型的配筋鋼纖維混凝土結構承載性能研究
深海沉積物的靜力特性研究
蠶寶寶奉獻的一生
水穩填充大粒徑碎石材料路用性能研究
主站蜘蛛池模板: 国产精品一线天| 亚洲无码91视频| 香蕉网久久| 免费一级毛片| 免费又黄又爽又猛大片午夜| 国产91丝袜在线播放动漫| 婷婷综合缴情亚洲五月伊| 制服丝袜无码每日更新| 国产真实乱人视频| 日本高清成本人视频一区| 永久免费无码日韩视频| 国产v精品成人免费视频71pao| 无码高清专区| 香蕉视频在线观看www| 欧美日韩北条麻妃一区二区| 综1合AV在线播放| 日韩欧美国产另类| 91www在线观看| 精品国产Av电影无码久久久| 五月婷婷精品| 亚洲综合狠狠| 夜精品a一区二区三区| 狠狠做深爱婷婷久久一区| 日本AⅤ精品一区二区三区日| 99视频在线看| 男人天堂亚洲天堂| 亚洲高清中文字幕在线看不卡| 久久狠狠色噜噜狠狠狠狠97视色| 狂欢视频在线观看不卡| 欧美伊人色综合久久天天| 精品国产一二三区| 色综合久久88色综合天天提莫| 国产女人喷水视频| h网站在线播放| 国产在线观看第二页| 亚洲青涩在线| 国产成人一区免费观看| 午夜一区二区三区| 免费国产不卡午夜福在线观看| 国产小视频免费| 欧美中文字幕无线码视频| 久久综合干| 国产成人精品18| 国产91视频免费观看| 亚洲三级成人| a级毛片网| 白浆视频在线观看| 欧美成人A视频| 精品视频福利| 91国内视频在线观看| 久久综合结合久久狠狠狠97色| 久久人搡人人玩人妻精品 | 直接黄91麻豆网站| 中文字幕 日韩 欧美| 国产精品免费p区| 在线a视频免费观看| 福利小视频在线播放| 亚洲天堂精品在线观看| 成年人免费国产视频| 无码 在线 在线| 国产精品妖精视频| 亚洲精品在线观看91| 亚洲福利视频网址| 久久亚洲国产视频| 99热这里只有精品免费国产| 亚洲成a人片77777在线播放| 国产SUV精品一区二区6| 国产嫩草在线观看| 日韩欧美中文字幕一本| 久久人人妻人人爽人人卡片av| 日韩区欧美国产区在线观看| 专干老肥熟女视频网站| 亚洲欧美日韩天堂| 国产精品综合色区在线观看| 亚洲欧洲日韩国产综合在线二区| 中文纯内无码H| yjizz视频最新网站在线| 青青青国产视频| 91视频精品| 亚洲精品在线影院| 国产欧美在线观看精品一区污| 青草视频久久|