• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    多場耦合下大體積混凝土初次蓄水的溫度應(yīng)力問題研究

    2016-06-07 00:53嚴(yán)俊魏迎奇蔡紅璩愛玉
    關(guān)鍵詞:溫度應(yīng)力大體積混凝土數(shù)學(xué)模型

    嚴(yán)俊+魏迎奇+蔡紅+璩愛玉

    摘 要:在大體積混凝土壩初次蓄水時(shí),溫度較低的庫水必然會(huì)對(duì)壩體溫度場產(chǎn)生較明顯的影響,從而影響壩體的變形,甚至產(chǎn)生溫度裂縫.為分析大體積混凝土初次蓄水的溫度應(yīng)力,本文將混凝土類多孔介質(zhì)視為連續(xù)介質(zhì),綜合運(yùn)用水力學(xué)、熱學(xué)和固體力學(xué)等基本理論,根據(jù)動(dòng)量守恒、質(zhì)量守恒和能量守恒方程建立了以位移、孔隙水壓力、孔隙氣壓力、溫度和孔隙率為未知量的多場耦合數(shù)學(xué)模型,在此基礎(chǔ)上編制了有限元計(jì)算程序,并對(duì)大體積碾壓混凝土塊的滲流場、溫度場和應(yīng)力場進(jìn)行了耦合分析,結(jié)果表明,考慮耦合后塊體溫降幅度及溫度大主應(yīng)力均較不考慮耦合條件下大.

    關(guān)鍵詞:多場耦合;數(shù)學(xué)模型;大體積混凝土;溫度應(yīng)力

    中圖分類號(hào):TU43 文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1674-2974(2016)05-0030-09

    Abstract: During the initial impoundment of mass concrete dam, the reservoir water with low temperature would be critical to the temperature field of dam body, which affects the deformation of dam body and even results in temperature cracks. Therefore, in order to investigate the thermal stress distribution of mass concrete during the initial impoundment, the concrete was assumed as a continuous porous media in this paper. In consideration of the basic theories of solid mechanics, hydraulics, and thermodynamics, the multi-field coupling equations of unsaturated porous media that include momentum conservation, mass conservation, and energy conservation were provided as the function of displacements, pore liquid pressure, pore gas pressure, temperature, and porosity. The finite element analysis program was then developed. A mass roller-compacted concrete block was considered for coupling analyses on the seepage field, temperature field, and stress field. The analysis results show that the temperature reduction and principle thermal stress of the concrete block considering the multi-field coupling process are greater than those without the coupling effect.

    Key words: multi-field coupling; mathematical model; mass concrete; thermal stress

    在水利工程中,尤其是對(duì)大體積混凝土而言,較大的溫度拉應(yīng)力會(huì)引起混凝土的溫度裂縫,因此溫度應(yīng)力是大體積混凝土需要額外關(guān)注的.在水庫初次蓄水時(shí),溫度較低的庫水必然會(huì)對(duì)壩體溫度場產(chǎn)生較明顯的影響,從而影響混凝土壩的變形,甚至產(chǎn)生溫度裂縫.如國際上,在Revslstoke[[1],Dworshak[[2]和Russel[[3[]等重力壩的上游表面曾經(jīng)產(chǎn)生過嚴(yán)重的劈頭裂縫,深入壩內(nèi)幾十米,有的甚至將整個(gè)壩段一分為二,產(chǎn)生嚴(yán)重漏水.目前一般都認(rèn)為是在施工過程中壩體的上游側(cè)如果出現(xiàn)了表面裂縫,水庫蓄水之后,經(jīng)過一段時(shí)間,表面裂縫突然大范圍地?cái)U(kuò)展,成為劈頭裂縫,尤其是通倉澆筑的混凝土重力壩更容易出現(xiàn)這種劈頭縫,分析其原因主要是由于這類壩沒有布置縱縫,不進(jìn)行二期冷卻,在水庫初次蓄水時(shí),壩體溫度還很高,與外界低溫的庫水之間形成較大的內(nèi)外溫差,容易使表面擴(kuò)展成為劈頭縫[[4].由此可見,初次蓄水對(duì)混凝土重力壩的溫度場和應(yīng)力場的影響是需要加以研究的,并為采取有效、合理地滲控措施提供科學(xué)依據(jù).

    事實(shí)上,如果將混凝土壩體也視為多孔介質(zhì),則蓄水時(shí)滲流的產(chǎn)生也是必然的,若考慮庫水入滲,壩體非穩(wěn)定溫度場的變化將更加復(fù)雜.考慮滲流因素的存在,對(duì)大體積混凝土壩體溫控防裂及應(yīng)力狀態(tài)研究都有直接意義.國內(nèi)已經(jīng)有學(xué)者用多場耦合的方法來求解大體積混凝土的滲流場溫度場耦合作用[[5-7[]、滲流場溫度應(yīng)力耦合作用[[8[]等.

    本文在總結(jié)國內(nèi)外學(xué)者對(duì)多孔介質(zhì)多場耦合的機(jī)理、數(shù)學(xué)模型建立及求解方法的研究成果基礎(chǔ)上,將混凝土視為多孔介質(zhì),綜合運(yùn)用固體力學(xué)、水力學(xué)、熱學(xué)等基本理論,結(jié)合多孔介質(zhì)的熱本構(gòu)關(guān)系以及孔隙流體的熱運(yùn)動(dòng)規(guī)律,建立了混凝土多相、多場全耦合數(shù)學(xué)模型方程組,并編制了有限元求解程序,對(duì)某一大體積混凝土模型初次蓄水后一段時(shí)間內(nèi)的滲流場、溫度場、應(yīng)力場進(jìn)行耦合分析,以初步闡釋大體積混凝土初次蓄水時(shí)滲流場對(duì)壩體溫度場、應(yīng)力場的影響過程.

    1 混凝土類多孔介質(zhì)概念模型及熱本

    構(gòu)關(guān)系模型

    1.1 混凝土概念模型

    國內(nèi)外有些學(xué)者提出了關(guān)于多孔介質(zhì)的多場耦合模型[[9-13[].事實(shí)上,非飽和狀態(tài)下的混凝土材料也可視為連續(xù)性多孔介質(zhì),其特征單元體主要由固相、液相和氣相三相構(gòu)成,其中固相介質(zhì)主要為經(jīng)硬化后形成具有堆聚結(jié)構(gòu)的復(fù)合物.特征單元體的體積V可以用式(1)表示:

    1.2 混凝土類多孔介質(zhì)熱本構(gòu)關(guān)系模型

    1.2.1 固相介質(zhì)的熱本構(gòu)關(guān)系

    當(dāng)溫度發(fā)生變化時(shí),混凝土類多孔介質(zhì)將由于受熱發(fā)生膨脹、受冷而發(fā)生收縮,從而產(chǎn)生應(yīng)變.如果多孔介質(zhì)熱變形受到邊界的約束,就會(huì)產(chǎn)生應(yīng)力.

    混凝土的熱本構(gòu)關(guān)系表達(dá)方式比較多[[14-18[],為了利用這些成熟的本構(gòu)關(guān)系,同時(shí)考慮溫度變化對(duì)多孔介質(zhì)的影響,可以給出以增量形式表述的熱本構(gòu)關(guān)系如式(3)所示:

    1.2.4 溫度場的熱傳導(dǎo)定律

    在研究溫度場的問題時(shí),廣義的Fourier定律[[21]作為基本定律,它是指在導(dǎo)熱過程中,單位時(shí)間內(nèi)通過給定截面所傳遞的熱量,正比例垂直于該截面方向上的溫度變化率,而熱量傳導(dǎo)的方向與溫度升高的方向相反,用數(shù)學(xué)表達(dá)式表示為:

    2 混凝土類多孔介質(zhì)多場耦合的控制方程

    混凝土類多孔介質(zhì)非飽和狀態(tài)下多場耦合的控制方程將主要基于上述非飽和多孔介質(zhì)的概念模型及其相應(yīng)的熱本構(gòu)關(guān)系[[22]:

    2.1 變形場控制方程

    對(duì)于一個(gè)非飽和多孔介質(zhì)的單元體來說,其準(zhǔn)靜態(tài)下的應(yīng)力平衡方程如式(7):

    3.2 有限元求解程序編制

    本文基于上述有限元離散成果利用FORTRAN語言研發(fā)了多孔介質(zhì)多場耦合求解程序THM-3D,該程序采用模塊式開發(fā),每塊均具有獨(dú)立、明確的功能含義,以充分滿足實(shí)際工程中問題復(fù)雜的要求,主要算法流程圖如圖1所示,限于篇幅,該程序的驗(yàn)證將在其他文章給出.

    4 工程算例

    假設(shè)有一大體積碾壓混凝土塊,如圖2所示,該混凝土塊高94.0 m,上游側(cè)假定有蓄水至77.0 m,下游側(cè)無水.

    該碾壓混凝土塊為澆筑式施工,在4月1日起開始澆筑,混凝土入倉溫度在氣溫的基礎(chǔ)上加3 ℃,分層澆筑至頂部后30 d拆模.壩體在430 d內(nèi)澆筑完成,之后將在10 d內(nèi)分3段快速蓄水至正常蓄水位:在第3 d蓄水至25.0 m;第6 d蓄水至50.0 m;在第10 d蓄至77.0 m.計(jì)算周期為蓄水至2 a,該混凝土塊體由3種級(jí)配的混凝土澆筑而成,從塊體上游至下游依次為變態(tài)混凝土、二級(jí)配混凝土和三級(jí)配混凝土.考慮到大體積混凝土主要關(guān)注的是溫度應(yīng)力的變化,因此應(yīng)力場將主要考慮溫度應(yīng)力.

    4.2 計(jì)算結(jié)果分析

    4.2.1 混凝土澆筑完成

    混凝土塊澆筑完成時(shí)的溫度和應(yīng)力計(jì)算結(jié)果如下圖3(a),(b)所示,其中溫度單位為℃,應(yīng)力單位為MPa.可以看出,在澆筑完成時(shí),混凝土塊體溫度場的高溫區(qū)位于中下部和頂部位置,其中中下部的最高溫度達(dá)到44 ℃,主要是由于本算例中沒有考慮溫控設(shè)施,該處的澆筑溫度較高,且混凝土仍處于升溫階段,水泥水化熱產(chǎn)生的熱量無法散出;同時(shí),塊體中均為拉應(yīng)力,但整體水平不高.

    4.2.2 初次蓄水一段時(shí)間

    混凝土塊澆筑初次蓄水30 d,2 a的孔隙水壓力、溫度和應(yīng)力計(jì)算結(jié)果如圖4(a)~(c)和5(a)~(c)所示.

    從滲流場的計(jì)算結(jié)果來看,隨著蓄水時(shí)間的延長和庫水溫、氣溫的變化,庫水在塊體內(nèi)逐漸向下游側(cè)滲流,蓄水2 a后庫水仍然沒有在塊體的下表面逸出,說明蓄水后混凝土塊體內(nèi)的滲流場在較長的時(shí)間內(nèi)不會(huì)達(dá)到穩(wěn)定狀態(tài),也不會(huì)有水從下游面逸出.

    從溫度場的計(jì)算結(jié)果來看,在水庫蓄水后,隨著蓄水時(shí)間的延長,塊體上游側(cè)的溫度等值線逐漸向內(nèi)部延伸,這體現(xiàn)出了庫水對(duì)塊體溫度的影響過程.同時(shí),塊體溫度場在庫水作用下,內(nèi)部溫度整體水平下降,并且上游側(cè)降溫比下游側(cè)明顯較快;塊體上游庫水以上部位在庫水和外界氣溫共同影響下溫度下降較快,此處的最高溫度由澆筑完成時(shí)的38 ℃降至蓄水2 a后的18 ℃.

    從耦合條件下塊體溫度應(yīng)力分布的變化過程可以看出,隨著蓄水時(shí)間的延長,塊體上游側(cè)表面處的溫度由于受較低的庫水溫控制,與內(nèi)部相鄰部位的混凝土之間存在著溫差,蓄水后不同時(shí)期上表面出現(xiàn)了不同程度的拉應(yīng)力,局部拉應(yīng)力較大,可能對(duì)該部位的溫控防裂不利.

    4.2.3 特征點(diǎn)溫度變化過程

    特征點(diǎn)1~4在耦合與非耦合條件下的溫度計(jì)算結(jié)果對(duì)比如圖6~9所示.通過對(duì)比可以看出,特征點(diǎn)1在非耦合條件下得到的溫度較耦合條件高,最大差值為1.5 ℃;特征點(diǎn)2在非耦合條件下的得到的溫度較耦合條件高,最大差值為1.8 ℃;特征點(diǎn)3在非耦合條件下得到的溫度較耦合條件高,最大差值為3.3 ℃,;特征點(diǎn)4在非耦合條件與耦合條件下得到的溫度計(jì)算結(jié)果吻合得較好,最大差值約0.56 ℃.可見,非耦合的方法對(duì)于上游面附近庫水浸沒區(qū)域以及其他庫水未浸沒區(qū)域的溫度計(jì)算結(jié)果較好,對(duì)于被庫水浸沒的塊體內(nèi)部區(qū)域的溫度計(jì)算結(jié)果則明顯高于耦合計(jì)算得到的結(jié)果.

    5 結(jié)論與建議

    混凝土類多孔介質(zhì)多場耦合作用的數(shù)學(xué)模型,研究的主要是某一物理場方程中因變量或源匯項(xiàng)受其他物理場作用其變化的數(shù)學(xué)描述,也包括本構(gòu)規(guī)律的影響在控制方程中的反映,因此該類數(shù)學(xué)模型較為復(fù)雜,必定包含多個(gè)控制方程:

    1)本文根據(jù)連續(xù)介質(zhì)方法給出了混凝土類多孔介質(zhì)的概念模型,在給出基本假設(shè)的基礎(chǔ)上,依據(jù)熱本構(gòu)關(guān)系模型、流體運(yùn)動(dòng)的廣義Darcy定律以及溫度場熱傳導(dǎo)的廣義Fourier建立了以位移、孔隙水壓力、孔隙氣壓力、溫度、孔隙率為未知量的非飽和多孔介質(zhì)多相多場全耦合研究的數(shù)學(xué)模型,并給出了相應(yīng)的定解條件,獲得了有限元格式的多場耦合求解方程組,并編制了有限元求解程序.

    2)對(duì)一大體積碾壓混凝土塊的滲流場、溫度場和應(yīng)力場進(jìn)行了耦合分析.結(jié)果表明,考慮耦合后塊體入滲區(qū)域內(nèi)特征點(diǎn)的溫降幅度較不考慮耦合條件下大,最大溫降差達(dá)到3.3 ℃;而且考慮耦合效應(yīng)得到的特征點(diǎn)溫度大主應(yīng)力也較非耦合大,最大差值達(dá)到0.32 MPa.

    3)對(duì)于蓄水初期的實(shí)際大體積混凝土工程而言,影響其真實(shí)工作狀態(tài)的因素是極其復(fù)雜的,其中,作用在其上游側(cè)的低溫庫水則是眾多影響因素中最直接、最重要的因素.因此為了更為準(zhǔn)確地計(jì)算分析蓄水初期壩體溫度場和溫度應(yīng)力分布,應(yīng)該考慮滲流場的影響,本文采用耦合分析的方法為以后大體積混凝土壩蓄水后的溫控防裂研究提供了一種新的思路.

    參考文獻(xiàn)

    [1] BRUNNER W J ,WU K H.Cracking of revelstoke concrete gravity dam mass concrete[R].London, UR:15th International Congress on Large Dams,Vol.Ⅱ,1985:201-206.

    [2] HONGTON D L.Measures being taken for prevention of cracks in mass concrete at dworshak and libby dam[R]. Montreal, Canada:10th International Congress on Large Dams,1979:1062-1066.

    [3] ZHU B F. Allowable temperature difference, cooling capacity, inspection and treatment of cracks, and administration of temperature control[M].Beijing:Tsinghua University Press,2014: 439-467.

    [4] 朱伯芳.重力壩的劈頭裂縫[J].水力發(fā)電學(xué)報(bào),1997(59):85-93.

    [5] 柴軍瑞.混凝土壩滲流場與穩(wěn)定溫度場耦合分析的數(shù)學(xué)模型[J].水力發(fā)電學(xué)報(bào),2009,19(1):27-35.

    [6] 陳建余,朱岳明,張建斌.考慮滲流場影響的混凝土壩溫度場分析[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,2003,3(2):119-123.

    [7] 崔皓東,朱岳明. 蓄水初期的壩體非穩(wěn)定滲流場與溫度場耦合的理論模型及數(shù)值模擬[J].水利學(xué)報(bào),2009,15(2):238-243.

    [8] 李守義,陳堯隆,王長江.碾壓混凝土壩滲流對(duì)溫度應(yīng)力的影響[J].西安理工大學(xué)學(xué)報(bào),1996,12(5):41-46.

    [9] HICKS T W, PINE R J, WILLIS-RICHARDS J,et al. A hydro-thermo-mechanical numerical model for HDR geothermal reservoir[J].International Journal of Rock Mechanics and Sciences, 1996,33(5): 499-511.

    [10]CHAPULIOT S, GOURDIN C, PAYEN T,et al. Hydro-thermal-mechanical analysis of thermal fatigue in a mixing tee [J]. Nuclear Engineering and Design, 2005, 235(2):575-596.

    [11]MILLARD A,DURIN M,STIETEL A,et al.Discrete and continuum approaches to simulate the thermo-hydro-mechanical couplings in a large,fractured rock mass[J]. International Journal of Rock Mechanics and Sciences,1995,32(5):409-434.

    [12]NGUYEN T S, SELVADURAI A P S, ARMAND G. Modelling the FEBEX THM experiment using a state surface approach [J]. International Journal of Rock Mechanics & Mining Sciences,2005,42(5): 639-651.

    [13]RUTQVIST J, BRGESSON L, CHIJIMATSU M,et al. Thermohydromechanics of partially saturated geological media: governing equations and formulation of four finite element models [J]. International Journal of Rock Mechanics & Mining Sciences, 2001, 38(6): 105-127.

    [14]劉亞晨.核廢料貯存圍巖介質(zhì)THM耦合過程的數(shù)值模擬[J]. 地質(zhì)災(zāi)害與環(huán)境保護(hù),2006,17(2):78-82.

    [15]POLTRONIERI F, PICCOLROAZ A,BIGONI D,et al.A simple and robust elastoplastic constitutive model for concrete[J].Engineering Structures,2014, 60(5):81-84.

    猜你喜歡
    溫度應(yīng)力大體積混凝土數(shù)學(xué)模型
    活用數(shù)學(xué)模型,理解排列組合
    淺談構(gòu)建數(shù)學(xué)模型,建立千以內(nèi)數(shù)的數(shù)感
    淺談路橋混凝土的施工溫度與裂縫技術(shù)
    工業(yè)建筑工程中大體積混凝土的施工技術(shù)分析
    混凝土的施工溫度和裂痕控制
    淺談建筑工程大體積混凝土施工質(zhì)量控制
    對(duì)一個(gè)數(shù)學(xué)模型的思考
    “費(fèi)馬點(diǎn)”數(shù)學(xué)模型在中考中的應(yīng)用
    淺談房屋建筑工程大體積混凝土施工技術(shù)
    核主泵試驗(yàn)回路中熱水罐應(yīng)力分析研究
    久久午夜综合久久蜜桃| 性少妇av在线| 亚洲精品aⅴ在线观看| 男女下面插进去视频免费观看| 中文字幕亚洲精品专区| 久久韩国三级中文字幕| av国产久精品久网站免费入址| 国产av国产精品国产| 国产综合精华液| 精品一区在线观看国产| 午夜福利,免费看| 亚洲综合精品二区| av一本久久久久| 欧美精品av麻豆av| 国产日韩一区二区三区精品不卡| 午夜免费鲁丝| 久久午夜福利片| 国产精品国产三级专区第一集| 亚洲三级黄色毛片| 一边摸一边做爽爽视频免费| 日韩中字成人| 欧美亚洲 丝袜 人妻 在线| 亚洲av国产av综合av卡| 亚洲激情五月婷婷啪啪| 亚洲人成电影观看| 青春草视频在线免费观看| 美国免费a级毛片| 搡老乐熟女国产| 国产亚洲午夜精品一区二区久久| 午夜久久久在线观看| 宅男免费午夜| 蜜桃国产av成人99| 午夜老司机福利剧场| 国产高清国产精品国产三级| 色婷婷av一区二区三区视频| 国产在视频线精品| 久久久国产一区二区| 91aial.com中文字幕在线观看| 麻豆av在线久日| 午夜av观看不卡| 欧美老熟妇乱子伦牲交| 中文字幕亚洲精品专区| 9热在线视频观看99| 免费观看性生交大片5| 十分钟在线观看高清视频www| 免费高清在线观看视频在线观看| 久久国产精品大桥未久av| 80岁老熟妇乱子伦牲交| 国产欧美日韩一区二区三区在线| av女优亚洲男人天堂| 久久韩国三级中文字幕| 一级毛片电影观看| 日韩制服骚丝袜av| 国产男女超爽视频在线观看| 日本黄色日本黄色录像| 人人妻人人澡人人爽人人夜夜| 亚洲中文av在线| 国产精品.久久久| 日本91视频免费播放| 人妻系列 视频| 亚洲国产色片| 秋霞在线观看毛片| 国产深夜福利视频在线观看| 免费少妇av软件| 久热这里只有精品99| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲内射少妇av| 如日韩欧美国产精品一区二区三区| 久久久久精品久久久久真实原创| 在线观看三级黄色| 9色porny在线观看| 性色avwww在线观看| 亚洲av电影在线观看一区二区三区| 欧美日韩视频精品一区| 亚洲视频免费观看视频| 亚洲成人一二三区av| 人人澡人人妻人| 久久综合国产亚洲精品| 99九九在线精品视频| 一二三四在线观看免费中文在| 国产成人欧美| 亚洲av免费高清在线观看| 精品一区在线观看国产| 制服人妻中文乱码| 欧美国产精品一级二级三级| 国产男人的电影天堂91| 在线观看三级黄色| 精品一区在线观看国产| 美女福利国产在线| 国产亚洲欧美精品永久| 肉色欧美久久久久久久蜜桃| 精品酒店卫生间| 香蕉丝袜av| 日产精品乱码卡一卡2卡三| 99热国产这里只有精品6| www.av在线官网国产| 欧美在线黄色| 晚上一个人看的免费电影| 日韩中文字幕欧美一区二区 | 亚洲av电影在线进入| 亚洲精品美女久久av网站| 亚洲视频免费观看视频| 国产又色又爽无遮挡免| 亚洲av电影在线进入| 老司机影院毛片| 精品一品国产午夜福利视频| 国产亚洲av片在线观看秒播厂| 久久久久久人人人人人| 日韩 亚洲 欧美在线| 最近手机中文字幕大全| 看免费av毛片| 在线观看国产h片| 男女边吃奶边做爰视频| 大码成人一级视频| av片东京热男人的天堂| 免费在线观看黄色视频的| 2022亚洲国产成人精品| 边亲边吃奶的免费视频| 久久人妻熟女aⅴ| 永久网站在线| 久久鲁丝午夜福利片| 在线观看一区二区三区激情| 亚洲精品久久成人aⅴ小说| av有码第一页| 午夜福利乱码中文字幕| 亚洲成人一二三区av| 老司机影院成人| 精品99又大又爽又粗少妇毛片| 午夜精品国产一区二区电影| 国产免费又黄又爽又色| 成人午夜精彩视频在线观看| 亚洲国产看品久久| 伦精品一区二区三区| 黄色配什么色好看| 纵有疾风起免费观看全集完整版| 永久免费av网站大全| 在线天堂最新版资源| 少妇的逼水好多| 亚洲国产精品国产精品| 久久久a久久爽久久v久久| 欧美激情 高清一区二区三区| 黑丝袜美女国产一区| 国产黄频视频在线观看| 日产精品乱码卡一卡2卡三| 成人毛片a级毛片在线播放| 精品国产乱码久久久久久小说| 欧美少妇被猛烈插入视频| 18禁裸乳无遮挡动漫免费视频| 色婷婷av一区二区三区视频| www.精华液| 免费观看无遮挡的男女| 老司机影院成人| 免费日韩欧美在线观看| 久久午夜综合久久蜜桃| 一边亲一边摸免费视频| 又粗又硬又长又爽又黄的视频| 伦理电影免费视频| av一本久久久久| 亚洲精品国产一区二区精华液| 国产精品免费大片| 亚洲精品国产一区二区精华液| 日韩av免费高清视频| 多毛熟女@视频| 天堂8中文在线网| 色播在线永久视频| 久久99热这里只频精品6学生| videossex国产| 韩国av在线不卡| 国产成人欧美| 蜜桃国产av成人99| 国产精品久久久久久久久免| 婷婷成人精品国产| 晚上一个人看的免费电影| 精品一区二区三卡| 日韩精品免费视频一区二区三区| 日韩制服丝袜自拍偷拍| 国产欧美日韩一区二区三区在线| 久久久久久人妻| 久久久久久人妻| 少妇精品久久久久久久| 国产日韩欧美在线精品| 国产免费福利视频在线观看| 欧美xxⅹ黑人| 免费在线观看黄色视频的| 搡老乐熟女国产| 99久久精品国产国产毛片| 色播在线永久视频| 亚洲国产精品成人久久小说| 久久精品国产鲁丝片午夜精品| 人人妻人人爽人人添夜夜欢视频| 亚洲av欧美aⅴ国产| 日韩制服骚丝袜av| 中文精品一卡2卡3卡4更新| 欧美日韩亚洲国产一区二区在线观看 | 精品国产国语对白av| 最近2019中文字幕mv第一页| 在现免费观看毛片| 亚洲精品一区蜜桃| 在线精品无人区一区二区三| 黄频高清免费视频| 交换朋友夫妻互换小说| 欧美日韩精品成人综合77777| 亚洲国产精品一区三区| 免费日韩欧美在线观看| 日日啪夜夜爽| 中文字幕人妻丝袜一区二区 | 久久精品夜色国产| av电影中文网址| 美女视频免费永久观看网站| 午夜福利一区二区在线看| 久久狼人影院| 国产精品二区激情视频| 久久女婷五月综合色啪小说| 在线看a的网站| 两性夫妻黄色片| 美女国产视频在线观看| 男女高潮啪啪啪动态图| xxxhd国产人妻xxx| 日韩免费高清中文字幕av| 久久久久精品性色| 久久人人97超碰香蕉20202| av网站在线播放免费| 日韩三级伦理在线观看| 欧美精品av麻豆av| 啦啦啦视频在线资源免费观看| 中文字幕精品免费在线观看视频| 亚洲国产精品成人久久小说| 久久久精品区二区三区| 丝瓜视频免费看黄片| 久久精品熟女亚洲av麻豆精品| 欧美黄色片欧美黄色片| 久久精品人人爽人人爽视色| 成人毛片60女人毛片免费| 中文欧美无线码| 午夜日本视频在线| 久久99一区二区三区| 狠狠婷婷综合久久久久久88av| 午夜福利在线免费观看网站| 国产女主播在线喷水免费视频网站| 91国产中文字幕| 国产麻豆69| 丰满乱子伦码专区| 少妇猛男粗大的猛烈进出视频| 国产黄色视频一区二区在线观看| 欧美日韩亚洲高清精品| 久久婷婷青草| 永久免费av网站大全| 各种免费的搞黄视频| 王馨瑶露胸无遮挡在线观看| 人妻人人澡人人爽人人| 欧美 亚洲 国产 日韩一| 久久久久久免费高清国产稀缺| 精品国产一区二区三区四区第35| 精品国产露脸久久av麻豆| 免费在线观看黄色视频的| 看免费av毛片| 侵犯人妻中文字幕一二三四区| 91成人精品电影| 国产精品一区二区在线观看99| 五月伊人婷婷丁香| 精品少妇久久久久久888优播| 国产乱来视频区| 青春草视频在线免费观看| 亚洲国产精品999| 精品酒店卫生间| 美女大奶头黄色视频| 美女国产视频在线观看| 国产精品久久久久久av不卡| 亚洲图色成人| 日韩熟女老妇一区二区性免费视频| 久久久久精品久久久久真实原创| 有码 亚洲区| 啦啦啦视频在线资源免费观看| 最近2019中文字幕mv第一页| 在线观看一区二区三区激情| 看免费成人av毛片| 欧美亚洲 丝袜 人妻 在线| 精品午夜福利在线看| 国产免费视频播放在线视频| 久久人人爽人人片av| 18禁国产床啪视频网站| 免费人妻精品一区二区三区视频| 久久热在线av| 三级国产精品片| 美女中出高潮动态图| 老汉色av国产亚洲站长工具| 欧美亚洲 丝袜 人妻 在线| 激情五月婷婷亚洲| 曰老女人黄片| 18+在线观看网站| 免费人妻精品一区二区三区视频| 欧美人与性动交α欧美精品济南到 | 女性生殖器流出的白浆| 青草久久国产| 欧美少妇被猛烈插入视频| 免费在线观看视频国产中文字幕亚洲 | 久久 成人 亚洲| 少妇人妻精品综合一区二区| 一级爰片在线观看| 久久99一区二区三区| 日本vs欧美在线观看视频| 在线观看国产h片| 免费日韩欧美在线观看| 丝瓜视频免费看黄片| 国产成人一区二区在线| 欧美精品一区二区大全| 日韩一卡2卡3卡4卡2021年| 美女xxoo啪啪120秒动态图| 老司机影院成人| 国产乱人偷精品视频| 亚洲精品美女久久av网站| 亚洲精品自拍成人| 大码成人一级视频| 午夜福利视频精品| 国产无遮挡羞羞视频在线观看| 美女中出高潮动态图| 大陆偷拍与自拍| 国产精品偷伦视频观看了| 多毛熟女@视频| 咕卡用的链子| 精品99又大又爽又粗少妇毛片| 亚洲国产精品999| 精品人妻偷拍中文字幕| 国产97色在线日韩免费| 久久韩国三级中文字幕| 亚洲成人一二三区av| 欧美黄色片欧美黄色片| 天天操日日干夜夜撸| 国产成人91sexporn| 午夜老司机福利剧场| 极品少妇高潮喷水抽搐| 国产av码专区亚洲av| 又大又黄又爽视频免费| 亚洲精品久久午夜乱码| 色婷婷久久久亚洲欧美| 波多野结衣av一区二区av| 成人毛片a级毛片在线播放| 一边摸一边做爽爽视频免费| 亚洲,欧美,日韩| 国产精品一国产av| 久久久国产一区二区| 女性被躁到高潮视频| 成人午夜精彩视频在线观看| 一区二区三区精品91| 精品少妇内射三级| 午夜免费观看性视频| 99热国产这里只有精品6| 男男h啪啪无遮挡| 丝袜在线中文字幕| videos熟女内射| 欧美日本中文国产一区发布| 母亲3免费完整高清在线观看 | 国精品久久久久久国模美| 伊人久久国产一区二区| 久久久久久久久久久久大奶| 黄片播放在线免费| 国产欧美日韩一区二区三区在线| 80岁老熟妇乱子伦牲交| 亚洲精品日韩在线中文字幕| 午夜免费鲁丝| 精品99又大又爽又粗少妇毛片| 日韩,欧美,国产一区二区三区| 中文字幕人妻熟女乱码| av国产久精品久网站免费入址| 美女午夜性视频免费| 青草久久国产| 午夜91福利影院| 在线观看免费视频网站a站| 狠狠婷婷综合久久久久久88av| 成人国产麻豆网| 成人国产av品久久久| kizo精华| 久久人妻熟女aⅴ| 久久久久视频综合| 色婷婷久久久亚洲欧美| 亚洲三级黄色毛片| 日韩av在线免费看完整版不卡| 91久久精品国产一区二区三区| 青春草国产在线视频| 国产一区二区三区综合在线观看| 久久精品久久久久久噜噜老黄| 毛片一级片免费看久久久久| www.自偷自拍.com| 亚洲人成77777在线视频| 国产精品亚洲av一区麻豆 | av国产久精品久网站免费入址| 丝瓜视频免费看黄片| 亚洲人成77777在线视频| 美女国产视频在线观看| 亚洲国产最新在线播放| 免费黄频网站在线观看国产| 日本-黄色视频高清免费观看| 看免费av毛片| 久久久久久人妻| 精品人妻熟女毛片av久久网站| 黄片无遮挡物在线观看| 校园人妻丝袜中文字幕| 三上悠亚av全集在线观看| 波野结衣二区三区在线| 成人国产av品久久久| 亚洲在久久综合| 日韩伦理黄色片| 天天躁日日躁夜夜躁夜夜| 成年动漫av网址| 人成视频在线观看免费观看| 国产片内射在线| 三上悠亚av全集在线观看| 满18在线观看网站| 亚洲精品久久午夜乱码| 国产又爽黄色视频| 国产熟女午夜一区二区三区| 天堂8中文在线网| 国产精品蜜桃在线观看| 久久国产精品大桥未久av| 一本—道久久a久久精品蜜桃钙片| 五月天丁香电影| 人人妻人人爽人人添夜夜欢视频| 亚洲一级一片aⅴ在线观看| 亚洲精品成人av观看孕妇| 色网站视频免费| 中文字幕色久视频| 欧美成人午夜精品| 夫妻性生交免费视频一级片| 国产精品欧美亚洲77777| 男人舔女人的私密视频| 91午夜精品亚洲一区二区三区| xxxhd国产人妻xxx| 亚洲三级黄色毛片| 亚洲精品日韩在线中文字幕| 日韩av免费高清视频| 久久久久精品性色| 人妻一区二区av| 久久这里只有精品19| av电影中文网址| 午夜福利,免费看| 精品午夜福利在线看| 国产精品蜜桃在线观看| 亚洲av国产av综合av卡| 午夜福利视频精品| 乱人伦中国视频| 永久网站在线| 国精品久久久久久国模美| 成年女人在线观看亚洲视频| 国产精品国产三级专区第一集| 久久午夜综合久久蜜桃| 97在线人人人人妻| 汤姆久久久久久久影院中文字幕| 亚洲综合精品二区| 国产在线视频一区二区| 国产男女内射视频| 伊人久久大香线蕉亚洲五| 成人毛片a级毛片在线播放| 青春草亚洲视频在线观看| 妹子高潮喷水视频| 免费观看在线日韩| 黄片播放在线免费| 国产精品香港三级国产av潘金莲 | 中文乱码字字幕精品一区二区三区| 国产成人午夜福利电影在线观看| 又黄又粗又硬又大视频| av线在线观看网站| 高清黄色对白视频在线免费看| 久久影院123| videos熟女内射| 超色免费av| 精品久久久精品久久久| 亚洲视频免费观看视频| 精品国产一区二区久久| 国产成人aa在线观看| 黄色一级大片看看| 熟女少妇亚洲综合色aaa.| 美女视频免费永久观看网站| 免费女性裸体啪啪无遮挡网站| 美女xxoo啪啪120秒动态图| 高清在线视频一区二区三区| 亚洲精品日韩在线中文字幕| 亚洲成人手机| 两个人免费观看高清视频| 免费播放大片免费观看视频在线观看| 老司机影院成人| 国产欧美亚洲国产| 日本猛色少妇xxxxx猛交久久| 久久久久国产精品人妻一区二区| 国产又爽黄色视频| 日韩精品有码人妻一区| 99九九在线精品视频| 亚洲少妇的诱惑av| 久久精品久久久久久噜噜老黄| 国产精品久久久久久久久免| 亚洲人成电影观看| 在线免费观看不下载黄p国产| 交换朋友夫妻互换小说| 免费日韩欧美在线观看| 国产av精品麻豆| 国产1区2区3区精品| 亚洲男人天堂网一区| 亚洲精品国产一区二区精华液| 三上悠亚av全集在线观看| a级毛片在线看网站| 性色avwww在线观看| 成年人免费黄色播放视频| 丰满迷人的少妇在线观看| 国产精品国产av在线观看| 男女边吃奶边做爰视频| 美女中出高潮动态图| 99久久中文字幕三级久久日本| 久久久久久人妻| 久久精品国产自在天天线| 性色av一级| 日韩人妻精品一区2区三区| 国产在视频线精品| 五月开心婷婷网| 日韩成人av中文字幕在线观看| 亚洲人成网站在线观看播放| 一边亲一边摸免费视频| 亚洲三级黄色毛片| 欧美精品国产亚洲| 国产xxxxx性猛交| 爱豆传媒免费全集在线观看| 日韩制服丝袜自拍偷拍| 亚洲色图 男人天堂 中文字幕| 人妻少妇偷人精品九色| 亚洲av电影在线进入| 久久综合国产亚洲精品| 国产成人精品一,二区| 亚洲av成人精品一二三区| 女的被弄到高潮叫床怎么办| 欧美日韩视频精品一区| 日韩中文字幕欧美一区二区 | 午夜激情av网站| 一级片免费观看大全| 九色亚洲精品在线播放| 久久免费观看电影| 欧美国产精品va在线观看不卡| 性高湖久久久久久久久免费观看| 丰满迷人的少妇在线观看| 亚洲视频免费观看视频| 黑人猛操日本美女一级片| 国产又爽黄色视频| 久久久久久伊人网av| 国产欧美日韩一区二区三区在线| 国产一区二区三区av在线| 国产精品一二三区在线看| 日产精品乱码卡一卡2卡三| 秋霞在线观看毛片| 18在线观看网站| 少妇 在线观看| 视频在线观看一区二区三区| 中文乱码字字幕精品一区二区三区| 欧美激情高清一区二区三区 | 久久人人爽av亚洲精品天堂| 欧美变态另类bdsm刘玥| 黑人猛操日本美女一级片| 在线观看免费视频网站a站| www.精华液| 午夜福利一区二区在线看| www.熟女人妻精品国产| 午夜福利乱码中文字幕| 日韩中文字幕视频在线看片| 天天躁狠狠躁夜夜躁狠狠躁| 色婷婷av一区二区三区视频| 丝瓜视频免费看黄片| 国产极品粉嫩免费观看在线| 久热久热在线精品观看| 街头女战士在线观看网站| 激情视频va一区二区三区| 久久99精品国语久久久| 麻豆乱淫一区二区| 欧美亚洲 丝袜 人妻 在线| 日韩一区二区三区影片| 亚洲精品久久久久久婷婷小说| 精品少妇一区二区三区视频日本电影 | 啦啦啦视频在线资源免费观看| 国产探花极品一区二区| 精品第一国产精品| 午夜福利,免费看| 一区二区三区四区激情视频| 91精品三级在线观看| 国产亚洲最大av| 考比视频在线观看| 激情五月婷婷亚洲| 赤兔流量卡办理| 有码 亚洲区| 美女中出高潮动态图| 欧美精品国产亚洲| 巨乳人妻的诱惑在线观看| 菩萨蛮人人尽说江南好唐韦庄| 九九爱精品视频在线观看| 最黄视频免费看| 国产欧美日韩一区二区三区在线| 美女主播在线视频| 欧美亚洲日本最大视频资源| 国产 一区精品| 亚洲国产欧美日韩在线播放| 久久精品亚洲av国产电影网| 97精品久久久久久久久久精品| 美女大奶头黄色视频| 在线免费观看不下载黄p国产| 午夜日本视频在线| 国产无遮挡羞羞视频在线观看| 欧美日韩精品成人综合77777| 在线亚洲精品国产二区图片欧美| 菩萨蛮人人尽说江南好唐韦庄| 巨乳人妻的诱惑在线观看| 国产激情久久老熟女| 日韩精品免费视频一区二区三区| 国产精品人妻久久久影院| 欧美在线黄色| 色哟哟·www| 精品人妻在线不人妻| 久久影院123| 欧美激情极品国产一区二区三区| 亚洲成人手机| 永久网站在线| 最近的中文字幕免费完整|