• <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)力分析研究
    最近最新中文字幕大全电影3 | 日韩有码中文字幕| 中国美女看黄片| 性色av乱码一区二区三区2| 狠狠狠狠99中文字幕| 久久精品国产亚洲av高清一级| 满18在线观看网站| 1024视频免费在线观看| 两性夫妻黄色片| 国产熟女午夜一区二区三区| 99久久国产精品久久久| e午夜精品久久久久久久| 成人三级做爰电影| 精品第一国产精品| 99国产精品99久久久久| 国产黄片美女视频| 国产野战对白在线观看| 国产精品精品国产色婷婷| 亚洲成人久久爱视频| 久久久久久久久免费视频了| 免费无遮挡裸体视频| 成人国产一区最新在线观看| 国产精品av久久久久免费| 精品欧美国产一区二区三| 日本在线视频免费播放| 又紧又爽又黄一区二区| 亚洲成人国产一区在线观看| 禁无遮挡网站| 91老司机精品| 日韩中文字幕欧美一区二区| 亚洲精品美女久久久久99蜜臀| www.www免费av| 久久伊人香网站| www日本在线高清视频| 老熟妇仑乱视频hdxx| 中文字幕人妻丝袜一区二区| 久久九九热精品免费| 亚洲欧洲精品一区二区精品久久久| 久久精品国产亚洲av高清一级| 久99久视频精品免费| 久久青草综合色| 国产99白浆流出| 久久久久久久久久黄片| 亚洲成av人片免费观看| 国产激情偷乱视频一区二区| 午夜免费成人在线视频| 亚洲欧美一区二区三区黑人| 国产精品久久电影中文字幕| 久久国产精品影院| 久久99热这里只有精品18| 欧美黑人巨大hd| 亚洲精品国产一区二区精华液| 欧美绝顶高潮抽搐喷水| 欧美在线黄色| 丝袜人妻中文字幕| 天堂影院成人在线观看| 亚洲中文日韩欧美视频| 免费人成视频x8x8入口观看| 亚洲精品美女久久久久99蜜臀| 亚洲,欧美精品.| 日韩 欧美 亚洲 中文字幕| 香蕉av资源在线| 国产精品久久久久久精品电影 | 国产精品久久久久久人妻精品电影| 免费在线观看视频国产中文字幕亚洲| 午夜久久久久精精品| 我的亚洲天堂| 国产在线精品亚洲第一网站| 真人做人爱边吃奶动态| 亚洲国产日韩欧美精品在线观看 | 美女国产高潮福利片在线看| 老司机靠b影院| 精华霜和精华液先用哪个| 黄片小视频在线播放| 制服丝袜大香蕉在线| 成人手机av| 国产区一区二久久| 欧美大码av| www日本在线高清视频| 伊人久久大香线蕉亚洲五| 亚洲精品色激情综合| 久久久久久久精品吃奶| 免费在线观看亚洲国产| 国产在线精品亚洲第一网站| 国产精品国产高清国产av| 日韩精品免费视频一区二区三区| 人成视频在线观看免费观看| 九色国产91popny在线| 动漫黄色视频在线观看| 91成年电影在线观看| 国产日本99.免费观看| 老司机深夜福利视频在线观看| 亚洲国产精品999在线| 日日干狠狠操夜夜爽| 成人三级做爰电影| 欧美成人午夜精品| 天堂√8在线中文| 免费在线观看黄色视频的| 精品高清国产在线一区| 97人妻精品一区二区三区麻豆 | 中文字幕高清在线视频| 欧美在线黄色| 给我免费播放毛片高清在线观看| 久久 成人 亚洲| 韩国精品一区二区三区| 国产精品久久久久久人妻精品电影| 久久久久精品国产欧美久久久| 国产视频一区二区在线看| 深夜精品福利| 老司机午夜十八禁免费视频| 777久久人妻少妇嫩草av网站| 免费电影在线观看免费观看| 日韩有码中文字幕| 婷婷亚洲欧美| 亚洲成人国产一区在线观看| 中文字幕精品亚洲无线码一区 | 国产一区二区在线av高清观看| 热99re8久久精品国产| 黑人操中国人逼视频| 波多野结衣av一区二区av| 亚洲 国产 在线| 精品久久久久久,| 啪啪无遮挡十八禁网站| 免费一级毛片在线播放高清视频| 亚洲国产精品999在线| 99久久无色码亚洲精品果冻| 可以在线观看毛片的网站| 波多野结衣高清无吗| 少妇裸体淫交视频免费看高清 | 精品国产乱子伦一区二区三区| 午夜福利高清视频| 日日干狠狠操夜夜爽| 国产人伦9x9x在线观看| 色综合站精品国产| 三级毛片av免费| 国产成+人综合+亚洲专区| 国产不卡一卡二| 一区二区三区精品91| 真人做人爱边吃奶动态| 男女做爰动态图高潮gif福利片| 亚洲精品久久国产高清桃花| 啦啦啦韩国在线观看视频| 最好的美女福利视频网| 桃红色精品国产亚洲av| 欧美绝顶高潮抽搐喷水| 天堂√8在线中文| 最新在线观看一区二区三区| 国产亚洲精品久久久久久毛片| 亚洲熟妇中文字幕五十中出| 可以免费在线观看a视频的电影网站| av免费在线观看网站| 亚洲av美国av| 成年女人毛片免费观看观看9| 国产精品99久久99久久久不卡| 久久久久久久久久黄片| 亚洲人成伊人成综合网2020| 精品不卡国产一区二区三区| 免费看美女性在线毛片视频| xxx96com| 97人妻精品一区二区三区麻豆 | 精品国产国语对白av| 国产一区在线观看成人免费| 中文字幕另类日韩欧美亚洲嫩草| 色综合站精品国产| 亚洲一码二码三码区别大吗| 男女之事视频高清在线观看| 国产单亲对白刺激| 一本大道久久a久久精品| 午夜影院日韩av| 国产野战对白在线观看| 亚洲国产欧洲综合997久久, | 成人av一区二区三区在线看| 在线永久观看黄色视频| 日韩av在线大香蕉| 法律面前人人平等表现在哪些方面| 欧美国产精品va在线观看不卡| 大型黄色视频在线免费观看| 久久中文看片网| 国产欧美日韩精品亚洲av| 一级片免费观看大全| 狂野欧美激情性xxxx| 中文字幕高清在线视频| 日韩免费av在线播放| 哪里可以看免费的av片| 成人手机av| 两个人看的免费小视频| 亚洲五月色婷婷综合| 男人舔女人下体高潮全视频| 精品国产乱码久久久久久男人| x7x7x7水蜜桃| 丁香六月欧美| 国产高清视频在线播放一区| 岛国视频午夜一区免费看| 又大又爽又粗| 最新美女视频免费是黄的| 伊人久久大香线蕉亚洲五| 法律面前人人平等表现在哪些方面| 亚洲精品国产区一区二| 99久久无色码亚洲精品果冻| 亚洲精品av麻豆狂野| 97人妻精品一区二区三区麻豆 | 少妇熟女aⅴ在线视频| 久久久久久久久久黄片| 亚洲国产中文字幕在线视频| 99久久久亚洲精品蜜臀av| 国产黄色小视频在线观看| 亚洲天堂国产精品一区在线| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人一区二区三区免费视频网站| 久久久久久免费高清国产稀缺| 久久人妻福利社区极品人妻图片| 级片在线观看| 亚洲第一av免费看| 精品熟女少妇八av免费久了| 女同久久另类99精品国产91| а√天堂www在线а√下载| 我的亚洲天堂| 一a级毛片在线观看| 中文字幕人妻熟女乱码| 欧美成人午夜精品| av视频在线观看入口| 亚洲欧美日韩高清在线视频| 美女午夜性视频免费| 久久久国产成人免费| 亚洲一区中文字幕在线| 俄罗斯特黄特色一大片| 黄片播放在线免费| 国产高清激情床上av| 女警被强在线播放| 亚洲国产高清在线一区二区三 | 亚洲中文日韩欧美视频| 亚洲av成人av| 国产av一区在线观看免费| 波多野结衣高清无吗| 99riav亚洲国产免费| 在线观看舔阴道视频| 国产91精品成人一区二区三区| 99久久精品国产亚洲精品| 欧美成人午夜精品| 日本 欧美在线| 看黄色毛片网站| 国产av一区二区精品久久| 很黄的视频免费| 亚洲中文字幕日韩| 国产精品久久久久久人妻精品电影| 一区二区三区国产精品乱码| 亚洲国产欧洲综合997久久, | 亚洲精品美女久久久久99蜜臀| 色av中文字幕| 亚洲国产毛片av蜜桃av| 久久这里只有精品19| 女同久久另类99精品国产91| 久久亚洲真实| 久久久久国产精品人妻aⅴ院| 黄色视频不卡| 久久热在线av| 国产又色又爽无遮挡免费看| 国产精品亚洲av一区麻豆| 麻豆久久精品国产亚洲av| 国产成人精品久久二区二区91| 波多野结衣巨乳人妻| 欧美午夜高清在线| 亚洲avbb在线观看| 午夜影院日韩av| 麻豆成人午夜福利视频| 动漫黄色视频在线观看| www.www免费av| 欧美 亚洲 国产 日韩一| 女生性感内裤真人,穿戴方法视频| 成人18禁在线播放| 久久精品夜夜夜夜夜久久蜜豆 | 国内久久婷婷六月综合欲色啪| 亚洲精品一区av在线观看| 久久久久国内视频| 成年版毛片免费区| 国产高清有码在线观看视频 | 757午夜福利合集在线观看| 国产精品九九99| 最近最新中文字幕大全电影3 | e午夜精品久久久久久久| 丝袜在线中文字幕| tocl精华| 日本黄色视频三级网站网址| 最近最新中文字幕大全免费视频| av电影中文网址| 国产一区二区激情短视频| 亚洲欧洲精品一区二区精品久久久| 18禁黄网站禁片免费观看直播| 久久这里只有精品19| 不卡一级毛片| 久久九九热精品免费| 亚洲欧美精品综合一区二区三区| 99久久无色码亚洲精品果冻| 亚洲一区中文字幕在线| 国产成人精品久久二区二区91| 欧美日本视频| 欧美 亚洲 国产 日韩一| 美女高潮到喷水免费观看| 中文在线观看免费www的网站 | 欧美成人性av电影在线观看| 不卡av一区二区三区| or卡值多少钱| 精品卡一卡二卡四卡免费| 免费高清视频大片| 日韩大尺度精品在线看网址| 亚洲久久久国产精品| svipshipincom国产片| 美女大奶头视频| 国产成人av教育| 国产午夜福利久久久久久| av超薄肉色丝袜交足视频| 一本大道久久a久久精品| 国产精品爽爽va在线观看网站 | 香蕉av资源在线| 99热这里只有精品一区 | 午夜日韩欧美国产| 国产亚洲欧美在线一区二区| 国产精品电影一区二区三区| 欧美激情极品国产一区二区三区| 啦啦啦韩国在线观看视频| 日韩欧美三级三区| 久久久久久久久免费视频了| 国产91精品成人一区二区三区| √禁漫天堂资源中文www| 国产又爽黄色视频| 欧美日韩亚洲国产一区二区在线观看| 在线观看一区二区三区| 午夜久久久在线观看| 亚洲真实伦在线观看| tocl精华| 午夜视频精品福利| 99精品久久久久人妻精品| 亚洲一区二区三区色噜噜| 一本久久中文字幕| 精品国产超薄肉色丝袜足j| www.www免费av| av免费在线观看网站| 69av精品久久久久久| 国产精品亚洲美女久久久| 黑丝袜美女国产一区| 日日摸夜夜添夜夜添小说| 中文字幕人妻丝袜一区二区| 欧美最黄视频在线播放免费| 久久久久久国产a免费观看| 啦啦啦观看免费观看视频高清| 亚洲,欧美精品.| 中文在线观看免费www的网站 | 嫁个100分男人电影在线观看| 精品福利观看| 久久天躁狠狠躁夜夜2o2o| 国产精品一区二区精品视频观看| 中国美女看黄片| 亚洲欧美精品综合一区二区三区| 日韩高清综合在线| 国产成人系列免费观看| 免费在线观看成人毛片| 欧美又色又爽又黄视频| 十分钟在线观看高清视频www| av中文乱码字幕在线| 在线看三级毛片| 欧美+亚洲+日韩+国产| 黄色片一级片一级黄色片| 91字幕亚洲| 国产成人av教育| 97碰自拍视频| aaaaa片日本免费| 免费看a级黄色片| 99久久精品国产亚洲精品| 国产av在哪里看| 非洲黑人性xxxx精品又粗又长| 国产精品综合久久久久久久免费| 18禁黄网站禁片午夜丰满| 一进一出抽搐gif免费好疼| 欧美日韩一级在线毛片| 欧美又色又爽又黄视频| 国产成人欧美在线观看| 日韩精品免费视频一区二区三区| 亚洲欧洲精品一区二区精品久久久| 香蕉av资源在线| 亚洲天堂国产精品一区在线| 十八禁网站免费在线| 亚洲在线自拍视频| av有码第一页| 99久久99久久久精品蜜桃| 嫩草影院精品99| 搡老熟女国产l中国老女人| 久久中文字幕人妻熟女| 又黄又爽又免费观看的视频| 色精品久久人妻99蜜桃| 久久香蕉精品热| 18禁裸乳无遮挡免费网站照片 | 久久久久久久精品吃奶| 麻豆成人av在线观看| 欧美成人一区二区免费高清观看 | 不卡av一区二区三区| 欧美成人性av电影在线观看| 一级a爱视频在线免费观看| 亚洲五月天丁香| 国产精品免费视频内射| 成年人黄色毛片网站| 国产成人欧美| 久久这里只有精品19| 男女午夜视频在线观看| 精品久久蜜臀av无| 亚洲片人在线观看| 精品第一国产精品| 在线天堂中文资源库| 麻豆一二三区av精品| 亚洲人成77777在线视频| 国产成+人综合+亚洲专区| 最近最新中文字幕大全免费视频| 日本成人三级电影网站| 久久久国产成人免费| www.自偷自拍.com| 久久草成人影院| 中文字幕人成人乱码亚洲影| 中文在线观看免费www的网站 | 在线观看日韩欧美| 国产男靠女视频免费网站| 狂野欧美激情性xxxx| 久久亚洲精品不卡| 九色国产91popny在线| 国产欧美日韩精品亚洲av| 99国产精品99久久久久| 一夜夜www| 国产精品98久久久久久宅男小说| 两性夫妻黄色片| 日韩免费av在线播放| xxxwww97欧美| 国产一区在线观看成人免费| 国产又黄又爽又无遮挡在线| 久久久久久国产a免费观看| 香蕉久久夜色| 亚洲五月天丁香| 动漫黄色视频在线观看| 午夜激情福利司机影院| 亚洲国产精品999在线| 欧洲精品卡2卡3卡4卡5卡区| 欧美黑人欧美精品刺激| 99久久99久久久精品蜜桃| 黄频高清免费视频| 国产亚洲av高清不卡| 好男人在线观看高清免费视频 | svipshipincom国产片| 岛国视频午夜一区免费看| 免费一级毛片在线播放高清视频| 黄频高清免费视频| 国产精品亚洲一级av第二区| 国产精品美女特级片免费视频播放器 | 国产区一区二久久| 国产一区在线观看成人免费| 观看免费一级毛片| 久久久久久大精品| 日本熟妇午夜| www.www免费av| 久久久久久人人人人人| 性欧美人与动物交配| av欧美777| 久久午夜综合久久蜜桃| 成人亚洲精品一区在线观看| 国产精品国产高清国产av| 嫁个100分男人电影在线观看| 午夜视频精品福利| 精品久久蜜臀av无| 香蕉av资源在线| 日日爽夜夜爽网站| 制服丝袜大香蕉在线| 黄色成人免费大全| 久9热在线精品视频| 午夜福利在线在线| 国产精品久久电影中文字幕| 老汉色av国产亚洲站长工具| 中文字幕精品免费在线观看视频| 这个男人来自地球电影免费观看| 男人舔女人的私密视频| 国产亚洲欧美精品永久| 天堂动漫精品| 久久久久免费精品人妻一区二区 | 欧美不卡视频在线免费观看 | 精品国产乱子伦一区二区三区| 亚洲av电影在线进入| 激情在线观看视频在线高清| 成人手机av| 99国产极品粉嫩在线观看| 成人三级做爰电影| 亚洲av成人av| 99国产精品一区二区三区| 十八禁人妻一区二区| 久久久久久九九精品二区国产 | 亚洲午夜精品一区,二区,三区| 91国产中文字幕| 日本五十路高清| 久久伊人香网站| 19禁男女啪啪无遮挡网站| 女人爽到高潮嗷嗷叫在线视频| 久久亚洲真实| 日本免费一区二区三区高清不卡| www国产在线视频色| 黄片小视频在线播放| 亚洲午夜理论影院| 免费看十八禁软件| 国产精品影院久久| 黄色丝袜av网址大全| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩瑟瑟在线播放| 成年人黄色毛片网站| 91成人精品电影| 午夜久久久久精精品| 国产黄a三级三级三级人| 日本五十路高清| 成熟少妇高潮喷水视频| АⅤ资源中文在线天堂| 欧美精品啪啪一区二区三区| 在线观看免费日韩欧美大片| 午夜免费成人在线视频| 中文在线观看免费www的网站 | 美女 人体艺术 gogo| www日本黄色视频网| 欧美久久黑人一区二区| 少妇裸体淫交视频免费看高清 | 免费看a级黄色片| www.熟女人妻精品国产| 欧美性猛交黑人性爽| 欧美精品亚洲一区二区| 美女 人体艺术 gogo| 黄色女人牲交| 热re99久久国产66热| 亚洲五月天丁香| 日韩精品青青久久久久久| 日本三级黄在线观看| 成人免费观看视频高清| 亚洲专区国产一区二区| 在线看三级毛片| 国产一卡二卡三卡精品| 欧美乱码精品一区二区三区| 亚洲五月婷婷丁香| 成人三级黄色视频| 精品一区二区三区av网在线观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲中文字幕一区二区三区有码在线看 | 亚洲性夜色夜夜综合| 国产亚洲av高清不卡| 免费在线观看亚洲国产| 久久中文看片网| 亚洲专区字幕在线| 大香蕉久久成人网| 成人三级黄色视频| 啦啦啦韩国在线观看视频| 人成视频在线观看免费观看| 成人国产一区最新在线观看| 天堂影院成人在线观看| 欧美激情高清一区二区三区| 一区二区三区激情视频| 老司机午夜福利在线观看视频| 长腿黑丝高跟| 亚洲av五月六月丁香网| 99热6这里只有精品| 18禁国产床啪视频网站| 久久久久免费精品人妻一区二区 | av免费在线观看网站| 亚洲人成伊人成综合网2020| 精品乱码久久久久久99久播| 欧美绝顶高潮抽搐喷水| 99精品久久久久人妻精品| 久热这里只有精品99| netflix在线观看网站| 色精品久久人妻99蜜桃| 午夜影院日韩av| 亚洲在线自拍视频| 久久久久久久精品吃奶| 亚洲av成人不卡在线观看播放网| 制服人妻中文乱码| 日本a在线网址| 久久久久久人人人人人| 高清在线国产一区| 欧美激情高清一区二区三区| 99久久久亚洲精品蜜臀av| 午夜成年电影在线免费观看| 国产精品自产拍在线观看55亚洲| 精品国产超薄肉色丝袜足j| 精品日产1卡2卡| 久久香蕉国产精品| 黑人巨大精品欧美一区二区mp4| 成年免费大片在线观看| www.999成人在线观看| 欧美激情久久久久久爽电影| 女警被强在线播放| 亚洲成人国产一区在线观看| 在线十欧美十亚洲十日本专区| 麻豆成人午夜福利视频| 国产蜜桃级精品一区二区三区| 91九色精品人成在线观看| 后天国语完整版免费观看| 国产高清激情床上av| 好看av亚洲va欧美ⅴa在| 久9热在线精品视频| 国产精品久久久人人做人人爽| 欧美乱色亚洲激情| 操出白浆在线播放| 欧美激情 高清一区二区三区| 99国产精品一区二区三区| 黄色片一级片一级黄色片| 免费在线观看完整版高清| 激情在线观看视频在线高清| 在线看三级毛片| АⅤ资源中文在线天堂| 久热爱精品视频在线9| 国产97色在线日韩免费| 亚洲av成人av| 侵犯人妻中文字幕一二三四区| 可以免费在线观看a视频的电影网站| av有码第一页| 精品国产亚洲在线| 国产久久久一区二区三区| av在线天堂中文字幕|