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

        ?

        GA+EA鋼橋面鋪裝復(fù)合結(jié)構(gòu)的高溫性能與力學(xué)特性

        2015-10-10 07:56:53胡靖錢振東楊宇明
        關(guān)鍵詞:變形結(jié)構(gòu)模型

        胡靖,錢振東,楊宇明

        ?

        GA+EA鋼橋面鋪裝復(fù)合結(jié)構(gòu)的高溫性能與力學(xué)特性

        胡靖,錢振東,楊宇明

        (東南大學(xué)智能運輸系統(tǒng)研究中心,江蘇南京,210096)

        研究鋼橋面澆注式瀝青混凝土與環(huán)氧瀝青混凝土(GA+EA)復(fù)合鋪裝結(jié)構(gòu)在荷載和溫度耦合下的高溫性能及力學(xué)特性。分析GA層的高溫流變參數(shù),采用多尺度與子模型有限元技術(shù)建立分析最不利溫度下復(fù)合鋪裝層結(jié)構(gòu)的壓應(yīng)力、剪應(yīng)力分布狀態(tài),并預(yù)估連續(xù)變溫條件下復(fù)合結(jié)構(gòu)的車轍深度及蠕變應(yīng)變隨時間變化情況。研究結(jié)果表明:EA層擴散了荷載中心GA層壓應(yīng)力,但其底受剪應(yīng)力較大。雙輪中央處GA層由于荷載疊加作用處于不利狀態(tài),其層底壓應(yīng)力達0.85 MPa。最不利連續(xù)變溫條件下,EA層變形量較小,GA層占鋪裝結(jié)構(gòu)永久變形90%以上,但總體車轍深度僅為0.32 mm。鋪裝結(jié)構(gòu)永久變形主要產(chǎn)生于夏季10:00?16:00高溫時段。GA+EA結(jié)構(gòu)較好地利用了各自材料的優(yōu)點,具有良好的高溫抗永久變形性能。

        復(fù)合結(jié)構(gòu);高溫性能;力學(xué)特征;溫度作用;多尺度子模型

        鋼橋面瀝青混凝土鋪裝是鋪筑在正交異性鋼橋面板上的薄層構(gòu)造,直接承受荷載、環(huán)境及結(jié)構(gòu)變形等綜合作用,受力條件十分不利?,F(xiàn)階段我國鋼橋面鋪裝大多采用雙層環(huán)氧瀝青混凝土以避免車轍等病害,但雙層環(huán)氧瀝青混凝土鋪裝在使用過程中會因抗疲勞性能不足而產(chǎn)生裂縫[1]。為防止環(huán)氧瀝青混凝土鋪裝層開裂,并考慮橋道系對鋪裝層變形追從性的要求,“上層環(huán)氧瀝青混凝土(EA)+下層澆注式瀝青混凝土(GA)”復(fù)合結(jié)構(gòu)從理論上將同樣具備優(yōu)良的使用性能。單純從鋪裝材料性能的角度分析,GA在高溫及荷載作用下容易產(chǎn)生車轍變形,因此有必要重點研究復(fù)合結(jié)構(gòu)在車載及溫度耦合作用下的力學(xué)響應(yīng)及高溫性能。國內(nèi)外已經(jīng)提出了一系列的本構(gòu)模型用以模擬瀝青混凝土的黏彈性特性,并以此為基礎(chǔ)對瀝青路面的高溫性能開展了大量研究。瀝青混凝土的高溫永久變形主要是由黏塑性變形引起,時間硬化蠕變模型能較好的擬合瀝青混凝土材料的變形特性[2]。利用Maxwell模型及有限元法對車轍深度進行預(yù)估,結(jié)果與加速加載試驗結(jié)果也基本吻合[3]。Mostafa等[4]采用Prony級數(shù)將瀝青混凝土的松弛模量轉(zhuǎn)化為蠕變?nèi)崃?,采用廣義Maxwell模型擬合得到材料參數(shù)。Burgers黏彈性模型用于車轍預(yù)估時,與試驗值也具有較好相關(guān)性[5?6]。鑒于EA已有相關(guān)試驗研究[7],本文作者對GA的黏彈性本構(gòu)進行試驗研究。

        1 復(fù)合結(jié)構(gòu)材料參數(shù)

        1.1 澆注式瀝青混合料動態(tài)模量試驗

        GA采用SBS改性瀝青為基質(zhì)瀝青,按特立尼達湖瀝青(TLA)和SBS改性瀝青質(zhì)量比為3:7配制的混合瀝青為瀝青結(jié)合料。級配選用德國0/8級配,油石比為9.2%。采用簡單性能試驗儀(SPT)對GA的動態(tài)模量進行研究。試件為100 mm×150 mm(直徑×高度)的圓柱體試件。試驗溫度為20,40和60 ℃。利用加載頻率為0.1,0.2,0.5,1,2,5,10,20和25 Hz正矢波進行軸向加載,測定相應(yīng)的應(yīng)力和軸向應(yīng)變,計算相應(yīng)頻率下的動態(tài)模量。圍壓為138 kPa,該圍壓能較好地模擬路面的受力狀態(tài),適合于三軸試驗[8]。按照AASHTO TP62試驗規(guī)程進行試驗,采用控制應(yīng)變方式加載,試驗結(jié)果如圖1所示。

        溫度/℃:1—20;2—40;3—60

        1.2 澆注式瀝青混合料流變性能

        作為典型的黏彈性材料,瀝青混合料的流變性能對溫度和時間具有很強的依賴性。采用SPT進行三軸重復(fù)荷載永久變形試驗獲取GA高溫下流變參數(shù)。試件直徑×高度為100 mm×150 mm。試驗溫度分別為20,40和60 ℃,偏應(yīng)力水平為0.7 MPa,試驗結(jié)果如圖2所示。

        溫度/℃:1—20;2—40;3—60

        1.3 有限元模型數(shù)值參數(shù)

        1.3.1 澆注式瀝青混合料流變學(xué)參數(shù)回歸

        瀝青混合料的永久變形通??梢员硎緸闇囟?、應(yīng)力和時間的函數(shù),即

        當(dāng)采用Bailey?Norton規(guī)律分析瀝青混合料的變形時,時間硬化蠕變模型的表達式為

        式中:為應(yīng)力;1,2和3為依賴于溫度的模型 參數(shù)。

        假定不隨時間變化,有

        則有

        表1 瀝青混合料的材料參數(shù)

        1.3.2 澆注式瀝青混合料流變學(xué)參數(shù)有效性驗證

        為驗證回歸得到的流變參數(shù)的有效性,利用ABAQUS建立三軸永久變形仿真模型,根據(jù)GA的蠕變參數(shù)計算圓柱試件的軸向應(yīng)變。單元類型為線性縮減積分單元C3D8R,限制模型底部的豎向位移,而僅在圓柱底部中心點處限制其側(cè)向位移。數(shù)值結(jié)果如圖3所示。

        圖3 數(shù)值結(jié)果與試驗結(jié)果對比

        由圖3可知:時間硬化蠕變模型對三軸重復(fù)荷載永久變形試驗結(jié)果具有較好的擬合度,因而采用該蠕變參數(shù)進行鋼橋面鋪裝的高溫性能研究是可行的。

        2 復(fù)合結(jié)構(gòu)有限元模型

        通過分析復(fù)合結(jié)構(gòu)在高溫條件下的溫度應(yīng)力,并以此為初始應(yīng)力場導(dǎo)入鋼橋面鋪裝力學(xué)分析模型,研究鋪裝層在溫度場與應(yīng)力場耦合作用下的力學(xué)響應(yīng)。

        本文采用多尺度建模技術(shù),數(shù)值模型的基本參數(shù)如表2所示。根據(jù)某大橋鋪裝形式,GA和EA厚度分別為35 mm和25 mm。按《公路橋涵設(shè)計通用規(guī)范》,取公路?Ⅰ級車輛荷載作用于鋼箱梁段,并采用“子模型”技術(shù)截取局部梁段模型,以獲取精確邊界條件。對于局部梁段模型,采用雙輪矩形均布荷載,單輪接地長×寬為20 cm×25 cm,兩輪中心矩為30 cm,接地壓力為0.7 MPa。最不利荷載中心位于2個U型加勁肋之間的中點,縱向位于相鄰兩橫隔板跨中[10],如圖4所示。

        表2 計算模型的基本參數(shù)

        (a) 鋼箱梁段有限元模型;(b) 局部梁段有限元模型

        鋪裝結(jié)構(gòu)溫度場選取夏季高溫季節(jié)的最不利時段,假設(shè)鋪裝層表面溫度為69 ℃,橋面鋼板溫度為58 ℃作為溫度場計算鋪裝結(jié)構(gòu)力學(xué)響應(yīng)。車轍預(yù)估時,則采用最不利季節(jié)(夏季)24 h溫度變化作為溫度場,鋪裝層熱力學(xué)參數(shù)采用相關(guān)文獻結(jié)果[11],如表3所示。

        表3 熱力學(xué)參數(shù)[11]

        車轍預(yù)估采用荷載作用時間累加的方法,將動載簡化為靜載作用[12]。通過施加持續(xù)恒載來模擬荷載的重復(fù)作用,按式(6)將荷載作用次數(shù)轉(zhuǎn)化為荷載累計作用時間。

        3 高溫條件下復(fù)合結(jié)構(gòu)力學(xué)響應(yīng)與車轍預(yù)估

        3.1 復(fù)合結(jié)構(gòu)力學(xué)響應(yīng)

        鋪裝結(jié)構(gòu)由于荷載作用產(chǎn)生蠕變和橫向剪切變形而形成車轍病害。分析最不利溫度(鋪裝表面溫度為69 ℃,橋面鋼板溫度為58 ℃)條件GA+EA復(fù)合結(jié)構(gòu)在荷載作用下,輪載作用中心與雙輪中央處的豎向壓應(yīng)力和橫向剪切應(yīng)力進行計算研究,最不利高溫條件下鋪裝層壓應(yīng)力和橫向剪應(yīng)力隨深度變化如圖6所示。

        (a) 壓應(yīng)力;(b) 橫向剪應(yīng)力

        由圖6(a)可以看出:左輪與右輪作用中心處的豎向壓應(yīng)力隨深度的變化趨勢幾乎相同,但雙輪中央處的壓應(yīng)力變化情況則完全相反。輪載作用處的EA層承受較大豎向壓應(yīng)力,但由于EA層對荷載的擴散作用,GA層所承受的豎向壓應(yīng)力出現(xiàn)了較大幅度的下降,由800 kPa下降到450 kPa,這表明該種復(fù)合結(jié)構(gòu)能較好的發(fā)揮材料各自的性能。但處于雙輪中央處的鋪裝層次,GA層承受了較大應(yīng)力,這主要是由于雙輪作用下應(yīng)力疊加所致。

        由圖6(b)可知:輪載中心及雙輪中央處的橫向剪應(yīng)力存在較大不同,但都隨著鋪裝深度的增加而增大,表明鋪裝層底部是最易出現(xiàn)剪切變形的區(qū)域。雙輪中心所受橫向剪應(yīng)力方向相反且數(shù)值相差較大,左輪最大剪應(yīng)力為80 kPa,而右輪僅為25 kPa。輪載中央處剪應(yīng)力較右輪中心更為不利,這是正交異性鋼橋板復(fù)雜受力情況導(dǎo)致鋪裝結(jié)構(gòu)受力的不規(guī)律性造成的。

        荷載作用下橫斷面上的應(yīng)力狀態(tài)分布情況也是研究復(fù)合結(jié)構(gòu)高溫性能的關(guān)鍵點,不同鋪裝深度輪載中心橫斷面豎向壓應(yīng)力和橫向剪應(yīng)力的平面分布情況如圖7所示,圖中橫軸代表雙輪荷載相對位置(雙輪荷載分別作用于0.1~0.3 m及0.4~0.5 m范圍內(nèi))。

        (a) 壓應(yīng)力;(b) 橫向剪應(yīng)力

        由圖7(a)可知,各鋪裝層豎向壓應(yīng)力發(fā)展規(guī)律相似,均出現(xiàn)“雙峰”形態(tài)。荷載中心處壓應(yīng)力隨著鋪裝深度逐漸變小,EA層的存在極大的分散了作用于GA層的壓應(yīng)力。同時,雙輪中央的壓應(yīng)力隨深度增加的速率更快,在GA層底部已經(jīng)達到850 kPa,而輪載中心處在GA底部的壓應(yīng)力僅為450 kPa,表明雙輪中央處GA層更容易出現(xiàn)高溫病害。圖7(b)表明:輪載作用下,鋪裝層斷面上橫向剪應(yīng)力變化趨勢類似,鋪裝深度越大,則橫向剪應(yīng)力越大,層底最大接近840 kPa。此外,GA層承受較為不利的橫向剪應(yīng)力狀態(tài),會造成較大的橫向剪切變形,但相鄰區(qū)域橫向剪應(yīng)力方向相反,對橫向剪切變形有一定的抵消作用。

        3.2 復(fù)合鋪裝結(jié)構(gòu)永久變形預(yù)估

        連續(xù)變溫條件下的抗變形能力是鋪裝結(jié)構(gòu)高溫性能的重要表現(xiàn),為了分析復(fù)合結(jié)構(gòu)在抗變形能力中的影響程度,計算最不利季節(jié)(夏季)24 h溫度變化下荷載作用50萬次后各結(jié)構(gòu)層的變形情況,軌跡處永久變形情況如圖8所示。

        圖8 輪跡處永久變形情況

        由圖8可知:EA層幾乎沒有產(chǎn)生車轍,鋪裝層的永久變形主要產(chǎn)生在GA層,其隆起變形近乎與總變形相等,而下陷變形則占總變形90%以上。由于GA的高溫性能相對較弱,在荷載的長期重復(fù)作用下,荷載作用處的鋪裝層產(chǎn)生了明顯的車轍變形,且最大蠕變應(yīng)變出現(xiàn)在鋪裝下層底部,位于U型加勁肋側(cè)肋與鋼橋面板相交處;由于EA層具有良好的整體性及應(yīng)力擴散能力,兩輪輪跡間也沒有產(chǎn)生明顯的凸起現(xiàn)象。此外,在荷載作用50萬次后,鋪裝結(jié)構(gòu)在高溫條件下產(chǎn)生的總車轍量為0.32 mm,表明EA+GA復(fù)合結(jié)構(gòu)具有良好的高溫抗永久變形性能。

        由于EA良好的整體性以及應(yīng)力擴散能力,使得鋪裝層在荷載作用下的蠕變變形及輪載中心處的隆起變形均較小,且增長趨勢十分緩慢,GA+EA鋪裝層連續(xù)變溫下蠕變應(yīng)變?時間關(guān)系曲線如圖9所示。

        1—隆起;2—下降;3—總形變

        在荷載作用3 100 s至5 500 s期間(按交通量日分布計算為10:00~16:00),不論是隆起還是下陷變形,鋪裝層的蠕變應(yīng)變在該時段增長速率都較快,分別占總形變的74%和51%,鋪裝層的永久變形有36%也集中產(chǎn)生在此高溫時段內(nèi),隨后進入蠕變穩(wěn)定階段。這是由于仿真初期鋪裝層溫度相對穩(wěn)定,鋪裝層混合料的蠕變參數(shù)取值較為一致。在加載時間為3 100 s至 5 500 s時,鋪裝層受環(huán)境氣溫影響而快速升溫,鋪裝層的蠕變變形呈現(xiàn)出較快的增長趨勢。在5 500 s之后,由于鋪裝層溫度回落,蠕變應(yīng)變增長十分微小。

        4 結(jié)論

        1) 雙輪荷載作用下剪應(yīng)力分布較壓應(yīng)力復(fù)雜。荷載中心處GA層所受壓應(yīng)力較小但橫向剪應(yīng)力較大。此外,由于雙輪荷載的應(yīng)力疊加作用,輪載中央處壓應(yīng)力較大。

        2) EA層擴散了GA層的壓應(yīng)力,但雙輪中央GA層底壓應(yīng)力達0.85 MPa,處于較不利狀況。GA層承受剪應(yīng)力較大,層底更達到0.84 MPa,雖相鄰區(qū)域剪應(yīng)力方向相反存在抵消作用,但仍需保證與鋼橋面板的黏結(jié)力。

        3) 鋪裝層的永久變形90%以上發(fā)生于GA層,而EA層幾乎不產(chǎn)生車轍,但總車轍量僅為0.32 mm表明EA+GA復(fù)合結(jié)構(gòu)具有良好的高溫抗永久變形性能。

        4) 鋪裝結(jié)構(gòu)的36%的永久變形產(chǎn)生于夏季10:00~16:00之間的高溫時段,在此時段應(yīng)對鋪裝層進行必要的降溫處理,以確保復(fù)合結(jié)構(gòu)良好的高溫性能。

        5) EA+GA復(fù)合結(jié)構(gòu)較好地利用了各自材料的優(yōu)點,EA層保證了路面的平整性,而GA層確保了鋪裝系與橋面系的協(xié)調(diào)變形。由于GA與EA高溫下蠕變應(yīng)變相差較大,因此有必要對高溫與重載耦合作用下GA層與EA層的層間黏結(jié)狀態(tài)及協(xié)調(diào)變形機理做進一步的研究,尋求最優(yōu)化鋪裝結(jié)構(gòu)形式。

        [1] 錢振東, 李智, 陳春紅. 鋼橋面環(huán)氧瀝青混凝土鋪裝層Ⅰ型裂縫的斷裂判據(jù)[J]. 中國公路學(xué)報, 2008, 21(5): 33?38. QIAN Zhendong, LI Zhi, CHEN Chunhong. Fracture criterion for mode Ⅰ crack of epoxy asphalt concrete paving course of steel deck bridge pavement[J]. China Journal of Highway and Transport, 2008, 21(5): 33?38.

        [2] Fwa T F, Tan S A, Zhu L Y. Rutting prediction of asphalt pavement layer using?model[J]. Journal of Transportation Engineering, 2004, 130(5): 675?683.

        [3] Blab R, Harvey J T. Viscoelastic rutting model with improved loading assumptions[C]// Proc 9th Int Conf on Asphalt Pavements, International Society for Asphalt Pavements. White Bear Lake, Minn, 2002: 293?303.

        [4] Mostafa A E, Imad A Q,Yoo P J. Viscoelastic Modeling and field validation of flexible pavements[J]. Journal of Engineering Mechanics, 2006, 132(2): 172?178.

        [5] 何兆益, 雷婷, 王國清,等. 基于動力問題的高等級公路瀝青路面車轍預(yù)估方法[J]. 土木工程學(xué)報, 2007, 40(3): 104?109. HE Zhaoyi, LEI Ting, WANG Guoqing, et al. Calculation method for prediction of highway asphalt pavement rutting under dynamic loading[J]. China Civil Engineering Journal, 2007, 43(3): 104?109.

        [6] 楊軍, 叢菱, 朱浩然,等. 鋼橋面瀝青混合料鋪裝車轍有限元分析[J]. 工程力學(xué), 2009, 26(5): 110?115. YANG Jun, CONG Ling, ZHU Haoran, et al. Analysis on rutting potential of asphalt pavement on the steel deck by finite element method[J]. Engineering Mechanics, 2009, 26(5): 110?115.

        [7] 陳磊磊. 環(huán)氧瀝青混合料動態(tài)參數(shù)試驗研究[D]. 南京: 東南大學(xué)交通學(xué)院, 2009: 12?40. CHEN Leilei. Research on dynamic parameters of epoxy asphalt mixture[D]. Nanjing: Southeast University. School of Transportation, 2009: 12?40.

        [8] Brown S F, Cooper K E. The Mechanical properties of bituminous materials for road bases and base courses[J]. Journal of the Association of Asphalt Paving Technologists, 1984, 53(4): 10?15.

        [9] 沈金安. 高速公路瀝青路面早期損壞分析與防治對策[M]. 北京: 人民交通出版社, 2004: 31?46. SHEN Jinan. Analysis and preventive techniques of premature damage of asphalt pavement in expressway[M]. Beijing: China Communication Press, 2004: 31?46.

        [10] 黃衛(wèi). 大跨徑橋梁鋼橋面鋪裝設(shè)計理論與方法[M]. 北京: 中國建筑工業(yè)出版社, 2006: 51?79. HUANG Wei. Theory and method of deck paving design for long-span bridge[M]. Beijing: China Architecture & Building Press, 2006: 51?79.

        [11] 錢振東, 劉龑. 整橋?溫度?重載耦合作用下鋼橋面黏結(jié)層力學(xué)分析[J]. 東南大學(xué)學(xué)報(自然科學(xué)版), 2012, 42(4): 729?733. QIAN Zhendong, LIU Yan. Mechanical analysis of waterproof bonding layer on steel bridge deck under bridge?temperature? load coupling effect[J]. Journal of Southeast University (Natural Science Edition), 2012, 42(4): 729?733.

        [12] 付凱敏. 瀝青路面結(jié)構(gòu)車轍模擬及抗車轍性能研究[D]. 南京: 東南大學(xué)交通學(xué)院, 2008: 25?50. FU Kai-min. Research on the rutting simulation and resistance characters of asphalt pavement structures[D]. Nanjing: Southeast University. School of Transportation, 2008: 25?50.

        High temperature properties of composite structure of gussasphalt+epoxy asphalt at bridge deck pavement

        HU Jing, QIAN Zhendong, YANG Yuming

        (Intelligent Transportation System Research Center, Southeast University, Nanjing 210096, China)

        The high temperature performance and mechanical properties of composite structure constructed by gussasphalt and epoxy asphalt (GA+EA) considering coupling effect of load and temperature were studied, the rheological parameters of GA were established, and the compressive stress and shear stress of paving structure under worst temperature were analyzed by multi-scale and sub-model finite element technology. Finally, the changes of rutting depth and creep strain at the situation of continuous temperature changing were predicted. The results show that the local vertical compress stress at the center of load of GA is diffused by EA, but the shear stress on bottom of GA is large. The vertical compress stress of GA bottom under the center of double load is 0.85 MPa because of the additive effect of load. Under worse temperature changing condition, the deformation of GA accounts for above 90% of that of paving structure; however, the total rutting depth is only 0.32 mm. The permanent deformation of paving structure mainly appears at 10:00?16:00. The composite paving structure of GA+EA can utilize the advantages of GA and EA, and it has excellent performance of anti-deformation.

        composite structure; high temperature performance; mechanical characteristics; temperature action; multi-scale sub-model

        10.11817/j.issn.1672-7207.2015.05.050

        U443.33;U416.2

        A

        1672?7207(2015)05?1946?07

        2014?03?10;

        2014?06?16

        國家自然科學(xué)基金資助項目(51178114);中央高?;究蒲袠I(yè)務(wù)費專項資金資助項目(CXLX12_0117);東南大學(xué)優(yōu)秀博士學(xué)位論文基金資助項目(YBJJ1318) (Project(51178114) supported by the National Natural Science Foundation of China; Project (CXLX12_0117) supported by the Fundamental Research Funds for the Central University; Project(YBJJ1318) supported by the Foundation for the Author of Excellent Doctoral Dissertation of Southeast University)

        錢振東,教授,從事道路結(jié)構(gòu)與材料研究;E-mail: qianzd@seu.edu.cn

        (編輯 趙俊)

        猜你喜歡
        變形結(jié)構(gòu)模型
        一半模型
        《形而上學(xué)》△卷的結(jié)構(gòu)和位置
        重要模型『一線三等角』
        談詩的變形
        中華詩詞(2020年1期)2020-09-21 09:24:52
        重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
        論結(jié)構(gòu)
        中華詩詞(2019年7期)2019-11-25 01:43:04
        “我”的變形計
        例談拼圖與整式變形
        會變形的餅
        論《日出》的結(jié)構(gòu)
        色欲av永久无码精品无码蜜桃| 中文字幕偷拍亚洲九色| 一区二区三区av资源网| 中文乱码字幕精品高清国产| 久久99精品久久久久久9蜜桃 | 丰满少妇被猛男猛烈进入久久| 妺妺窝人体色www聚色窝韩国| 日韩男女av中文字幕| 亚洲综合av一区二区三区蜜桃| 国产精品无码一本二本三本色| 亚洲中文无码久久精品1| 青青草免费在线手机视频| 日本av亚洲中文字幕| 人人妻人人澡人人爽欧美一区九九| 中文字幕人乱码中文字幕| 人妻丝袜无码国产一区| 伊人亚洲综合网色AV另类| 无码啪啪人妻| 日本免费精品一区二区| 先锋中文字幕在线资源| 日韩毛片基地一区二区三区| 精品免费看国产一区二区白浆| 国产精品国产三级久久| 日本免费a级毛一片| 国产白嫩美女在线观看| 国产av天堂亚洲国产av麻豆| 与最丰满美女老师爱爱视频| 亚洲国产精品无码专区| 中文字幕喷水一区二区| 日本一区二区久久精品亚洲中文无| 国产日本精品一二三四区| 狠狠色成人综合网| 欧美色资源| 午夜理论片日本中文在线| 亚洲精品一品区二品区三品区| 日韩中文字幕欧美亚洲第一区| 日韩丝袜人妻中文字幕| 国产成人综合久久久久久| 蜜臀久久99精品久久久久久小说| 欧美亚洲国产丝袜在线| 手机在线看片国产人妻|