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

    Visualizing the dynamic progression of backward erosion piping in a Hele-Shaw cell

    2022-11-30 09:51:36GangZHENGJingboTONGTianqiZHANGZiwuWANGXunLIJiqingZHANGChunyuQIHaizuoZHOUYuDIAO

    Gang ZHENG ,Jing-bo TONG ,Tian-qi ZHANG? ,Zi-wu WANG ,Xun LI ,Ji-qing ZHANG ,Chun-yu QI,Hai-zuo ZHOU,Yu DIAO

    1MOE Key Laboratory of Coast Civil Structure Safety,Tianjin University,Tianjin 300072,China

    2School of Science,Tianjin University,Tianjin 300072,China

    3China Railway Design Corporation,Tianjin 300308,China

    4National Engineering Laboratory for Digital Construction and Evaluation Technology of Urban Rail Transit,Tianjin 300308,China

    Abstract: With the utilization of underground space,backward erosion piping (BEP) has been observed in many underground structures(e.g.,shield tunnels)founded on sandy aquifers.However,due to invisibility,the geometry of the eroded pipe and its spatial evolution with time during the piping process was still not clear.In this study,we developed a Hele-Shaw cell to visualize the dynamic progression of BEP.With imaging process technology,we obtained a typical process of BEP(the erosion process can be divided into a piping progression phase and a piping stabilization phase),quantitatively characterized the formation of erosion pipes,and compared the patterns of erosion(e.g.,the erosion area A and the maximum erosion radius Rmax)that spontaneously develop under different fluxes of water.The most interesting finding is that the sand grains in a thicker Hele-Shaw model are easier to dislodge,which is possibly due to the granular system in a thicker model having more degrees of freedom,reducing the stability of the sand grains.

    Key words: Backward erosion piping(BEP);Dynamic progression;Hele-Shaw cell;Visualize;Imaging processing technology

    1 Introduction

    For shield tunnels built in a confined aquifer of silt or silty (fine) sand,once leakage occurs,the soil particles may fluidify,suspend,and then be carried into the tunnels by underground water,leading to the rapid erosion of the soil,changing the load distribu‐tion around the tunnels,and resulting in the deforma‐tion of the segmental linings.In particular,when the leakage point arises at the bottom of the tunnel,ero‐sion cavities can be generated below the tunnel invert,giving rise to both opening and dislocation at the seg‐ment joints.Thereafter,more leakage points arise and aggravate the erosion (Fig.1).Without prompt con‐trol,progressive failure and damage or even largescale tunnel collapse may take place.As shown in Fig.1,the erosion initiates at the outlet and develops backwards toward the fluid source.This is analogous to the backward erosion piping(BEP)observed in em‐bankments and levees founded on sandy aquifers(Fig.2).Therefore,for ease of interpretation,the erosion caused by leakage in the shield tunnels men‐tioned above is also termed BEP.In other words,the term BEP,in this study,refers to a generalized con‐cept and is not limited to a problem of embankments or levees.

    Fig.1 Erosion caused by leakage in a shield tunnel

    The typical process of BEP in embankments is shown in Fig.2.The erosion generally originates downstream and develops in a direction opposite to the water flow underneath the water-retaining struc‐tures.The critical state that occurs when the erosion pipe reaches the upstream is usually a primary focus because,once the erosion pipe forms a connection be‐tween upstream and downstream,the erosion can ac‐celerate,resulting in failure or collapse of the waterretaining structure.Previous experimental studies have yielded fundamental knowledge on key aspects of BEP either by measuring the critical average hydraulic gradient (Yao et al.,2007;van Beek et al.,2011) or by analysing pipe formation (Richards and Reddy,2012;Vandenboer et al.,2019a,2019b;Akrami et al.,2021;Rosenbrand et al.,2022).These studies have greatly broadened our insight into the mechanisms by which BEP can be triggered and sustained beneath an embankment.

    Apparently different from the BEP observed in embankments or levees,the erosion due to leakage in shield tunnels (Neyer,1984;Zheng et al.,2017) de‐velops in an approximately infinite space (in other words,the distance between the upstream and down‐stream is infinitely long).Therefore,no hydraulic connection between the upstream and downstream can be observed (Fig.2),and a critical state for the progression of BEP may not exist.Under these cir‐cumstances,the concept of a“critical hydraulic gradi‐ent” is meaningless.The geometry (shape and area)of erosion pipes,especially their spatial evolution with time,has become a major concern because these features are closely related to the response of the tunnel structures.For instance,if the erosion pipe is shallow,the deformation of the segmental lining may be characterised by opening,or dislocation might occur due to severe uneven settlement,as shown in Fig.1.However,due to its invisibility,the geometry of the erosion pipe underneath the tunnel is still not clear.

    Fig.2 Process of BEP in a river embankment

    Many attempts have been made to capture the geometry (length,depth,and width) of the eroded pipe beneath an embankment (van Beek et al.,2015;Vandenboer et al.,2018,2019b;Pol et al.,2021;Rice et al.,2021).In most of these studies,comprehensive and precise measurement of the pipe geometry was usually conducted at the end of the tests.It still seems a challenge to precisely capture the dynamic erosion morphology during pipe progression.Among the main factors that lead to difficulties in this measurement,the technical barrier in imaging and digital processing stands out.In conventional BEP experiments,the ac‐quisition of the erosion morphology is usually based on a top view of the sample.However,because the sample is thick while the pipe is shallow,the differ‐ence between an eroded and intact area is not obvious when viewed from the top.This makes the determina‐tion of the pipe contour difficult.

    Because the depth of the shallow pipe is only several times the grain diameter (van Beek et al.,2015;Vandenboer,2019),the use of a thin sample be‐tween two transparent plates for observation is feasi‐ble;this idea coincides with the core concept of the Hele-Shaw cell,a classic device used in physics to study laminar flows (Hele-Shaw,1898) and complex fluid-fluid (Al-Housseiny et al.,2012) or fluid-grain interactions (Cheng et al.,2008;Sandnes et al.,2011;Huang et al.,2012)between two parallel plates with a narrow gap.Inspired by the above concept,in this study a Hele-Shaw cell was designed and manufac‐tured for visualizing the progression of BEP.To mimic the erosion resulting from the leakage in a shield tunnel as much as possible,a radial configuration was adopted.The outlet was set at the centre of a circu‐lar Hele-Shaw cell,whereby the speed of the water would quickly decay with increasing radius.In combi‐nation with imaging processing technology,the geom‐etry of the eroded pipe and its spatial evolution during the piping process were quantitatively characterized.The findings provided intuitive insight into the dynamic progression of BEP.

    2 Experiments

    2.1 Experimental apparatus

    The experimental setup had three key compo‐nents:a radial Hele-Shaw cell,a standard volume con‐troller(SVC),and a high-speed camera(Phantom V7.1,AMETEK,USA;resolution: 2048×2048;frequency:60 Hz),as shown in Fig.3.The granular sample was accommodated in the Hele-Shaw cell.The fluid (pure water in this study)was supplied via the SVC.

    Fig.3 Experimental apparatus

    The whole BEP process was recorded by the high-speed camera.To eliminate the influence of am‐bient light on the test results,the experimental appara‐tus was placed in a photographic studio.

    A more detailed view of the Hele-Shaw cell is shown in Fig.4.The Hele-Shaw cell mainly consisted of two circular polymethylmethacrylate (PMMA)plates.The upper plate was 410 mm in diameter and 30 mm thick.The lower plate was 380 mm in diame‐ter and 29 mm thick.The gap between the two plates,b,was precisely controlled using standard spacers around the perimeter (in this study,bwas 0.5 mm in most experiments,but was changed to 0.25,0.70,1.00,or 1.50 mm in some tests to ascertain the effect of sample thickness).The outlet was a 6 mm-diameter hole,which was set at the centre of the upper plate.The two plates were inserted into a stainless-steel ped‐estal.The contact parts between the pedestal and the two plates were sealed using rubber O-rings.The inlet was set at the bottom of the pedestal,and the cavity that formed between the lower plates and the pedestal served as a reservoir,ensuring that the water was uni‐formly driven through the sample from around the pe‐rimeter.An annular aluminium strip embedded in the lower plate served as a ledge to separate the samples and the reservoir.

    Fig.4 Detailed view of the Hele-Shaw cell

    2.2 Materials

    The fluid used in this study was pure water with a density of 1000 kg/m3and a viscosity of 0.8737×10?3Pa·s.The granular material adopted in this study was Toyoura sand with a specific gravity of 2.65.The grading curve of Toyoura sand can be seen in Fig.5a.Because a good uniformity is essential to trigger BEP(van Beek et al.,2014),the Toyoura sand was sieved,and sand grains in a narrow range of size from 0.15 mm to 0.20 mm were selected for the samples(the scanning electron microscopy (SEM) image in Fig.5b).

    Fig.5 Grain size distribution of the Toyoura sand: (a)grading curve;(b)SEM image

    2.3 Test programme

    Before the test began,the sample was first transferred into the Hele-Shaw.More details on the preparation of the sample can be seen in the electronic supplementary materials(ESM).

    In this study,nine tests(excluding repeated tests)in total were conducted.They can be divided into two groups.In the first group,the gapbwas maintained at 0.5 mm,and the instant fluxJwas set to 140,210,280,350,and 420 mm3/s.In the second group,the ratio of the flux of water to the gapJ/bwas kept the same(280 mm2/s),whereas the gapbwas set to 0.25,0.50,0.70,1.00,and 1.50 mm.

    3 Test results

    3.1 A typical process of BEP

    Once the test began,pure water was injected into the Hele-Shaw cell via the SVC.To mimic the ero‐sion of a porous medium under a sudden flow of water as observed in many shield tunnels,a flux of waterJwas instantly imposed at the start and kept constant throughout the test.

    Fig.6 shows the process of BEP observed in an experiment,in which the flux of waterJwas 420 mm3/s and the gapbwas 0.5 mm.In general,the process of erosion can be divided into two phases: piping pro‐gression and piping stabilization.During the piping progression phase (I in Fig.6),the hydraulic pressurePover the sample(the hydraulic pressurePwas mea‐sured using the piezometer as shown in Fig.3) was shown to increase at first and then to decrease to a smaller value (Psta) after reaching a peak (Pmax).The spatial evolution of the erosion pipe was basically completed in phase I.Then,the BEP entered the piping stabilization phase (II in Fig.6).In phase II,the hydraulic pressure was maintained atPsta.Although the erosion area ceased to expand,transport of sand grains within the formed erosion pipe was still observed.

    From Fig.6,phase I lasted for 16 s in total(from 164 to 180 s).Fig.7 shows the morphology of the ero‐sion after each second.It is clear that under a sudden flow of water,the erosion pipe simultaneously pro‐gressed in multiple directions and exhibited bifurca‐tions as it grew in length(radius)at the start.With the progression of BEP,the incipiently formed pipes drew fluid and eroded faster at the expense of the neigh‐bouring regions.This led to competition among the erosion pipes whereby only one or two pipes eventually dominated (the development of pipe B1 and B2 after 171 s).

    Fig.6 Evolution of the erosion morphology corresponding to variations in the hydraulic pressure drop P measured over the sample during the progression of BEP

    Fig.7 Morphology of the erosion at each second in phase I

    Fig.8 shows the variations in both the erosion areaAand the maximum erosion radiusRmax(the max‐imum distance from the centre of the Hele-Shaw cell to the boundary of the erosion area)with time in phases I and II(the erosion areaAand the maximum erosion radiusRmaxwere obtained based on imaging process‐ing technology;more details can be found in the ESM).As can be seen,the erosion areaAlinearly increases with timeT,while the variations in the maxi‐mum erosion radiusRmaxwith timeTfit well to a para‐bolic curve.

    Fig.8 Variations in the erosion area A(a)and the maximum erosion radius Rmax(b)with time

    3.2 Progression of BEP under different fluxes of water

    In section 3.1,a typical BEP process due to a sudden flow of water is presented.To gain further in‐sight into the morphology of erosion under different fluxes of water,four additional tests were performed,in which the fluxes of water were set to 140,210,280,and 350 mm3/s.A comparison of the stable erosion morphology under different fluxes of water can be seen in Fig.9.In an intuitive sense,a greater flux of water may produce more branched patterns.For instance,whenJ=140 and 210 mm3/s,three main branches were observed,while forJ=350 and 420 mm3/s,four and five main branches were observed,respectively.This occurs because when the initial flux of water is set higher and well above the threshold required to trigger a BEP,more than one site of sand grains can be activated to move,and the erosion pipe may simul‐taneously progress in multiple directions,so the mor‐phology of the erosion exhibits bifurcations and a more branched pattern.

    Fig.9 Comparison of the stable erosion morphology under different fluxes of water

    Figs.10a and 10b further compare the erosion areaAand the maximum erosion radiusRmaxunder dif‐ferent fluxes of water.As expected,both the erosion areaAand the maximum erosion radiusRmaxwere larger when a higher flux of waterJwas employed.However,the increases inAandRmaxwithJwere not linear.It is apparent that whenJwas greater than 350 mm3/s,the increase almost stopped.In contrast,the peak hydrau‐lic pressurePmaxand the stable hydraulic pressurePstaincreased almost linearly with the flux of waterJ,as shown in Figs.10c and 10d.

    4 Effect of sample thickness

    In the above-mentioned experiments,the gap be‐tween the upper and lower plates was maintained at 0.5 mm.This meant that 2 or 3 layers of sand grains were arranged along the vertical direction in these experi‐ments.From Figs.6,7,and 9,it is clear that within the formed erosion area,the intensity of light was not uniform.There were some bright and thin channels from which sand grains were thoroughly removed and some dark areas in which one or two layers of sand grains were still left.This indicated that the BEP ob‐served in the Hele-Shaw cell was not simply a planestrain problem,which is the case in experiments to investigate Saffman-Taylor instability (Cheng et al.,2008;Huang et al.,2012) or the capillary fracturing of a fluid-grain system in a Hele-Shaw cell (Sandnes et al.,2011;Holtzman et al.,2012).As a consequence,the thickness of the sample might have an impact on the response of the erosive medium in the Hele-Shaw cell.To address this effect,four additional tests were performed in which gapbwas set to 0.25(monolayer),0.70(3 or 4 layers),1.00(5 or 6 layers),and 1.50 mm(7–10 layers).In particular,to ensure that the sand grains were eroded under similar water drag,the ratio of the flux of water to the gapJ/bwas kept the same(280 mm2/s) in these experiments (thus,the flux of waterJwas set to 70,140,196,280,and 420 mm3/s when the gaps wereb=0.25,0.50,0.70,1.00,and 1.50 mm,respectively).

    Fig.10 Comparison of A(a),Rmax(b),Pmax(c),and Psta(d)under different fluxes of water

    Fig.11 shows a comparison of the stable ero‐sion morphology obtained with different gaps in the Hele-Shaw cell.Apparently,the thickness of the sam‐ple can have a significant influence on the geometry of the erosion pipe.Basically,a thinner sample can show a greater resistance to erosion.For instance,in the case ofb=0.25 mm,BEP was not triggered.In ad‐dition,the erosion pattern looked more uniform when a larger gapbwas used in the Hele-Shaw cell.For in‐stance,in the case ofb=1.00 and 1.50 mm,the shape of the erosion boundary was well rounded,and minimal branching occurred.

    Fig.11 Comparison of the stable erosion morphologies obtained with different gaps in the Hele-Shaw cell

    Figs.12a and 12b further show a comparison of the erosion areaAand maximum erosion radiusRmaxobtained with different gaps in the Hele-Shaw cell.Both the erosion areaAand maximum erosion radiusRmaxwere greater with a larger gapbin the Hele-Shaw cell.In contrast,both the peak hydraulic pressurePmaxand the stable hydraulic pressurePstameasured over the sample were lower when a larger gapbwas em‐ployed,as shown in Figs.12c and 12d.All these re‐sults suggested that the sand grains in a thicker Hele-Shaw model were easier to dislodge.

    Fig.12 Comparison of A(a),Rmax(b),Pmax(c),and Psta(d)obtained in tests with different gaps

    Fig.13 illustrates the influence of the layer num‐ber on the erodibility of sand grains.For clarity,single grains on the erosion boundaries in a 2-layer sample and a 4-layer sample were selected for comparison.Apparently,the sand grain in a 4-layer sample more easily erodes because the granular system has more degrees of freedom than that of a 2-layer sample.In other words,a sand grain in a 4-layer sample ((1) in Fig.13)has more interrelated“neighbours”(the shaded grains in Fig.13)than a grain in a 2-layer sample((2)in Fig.13).Any perturbations or tiny wobbles of these interrelated “neighbours” can lead to the insta‐bility of the sand grain.

    Fig.13 Influence of the layer number on the erodibility of sand grains

    In addition,once a sand grain in the 4-layer sam‐ple is dislodged,it more easily maintains its motion due to the work done by gravity along a steeper slope(the arrow in Fig.13).Moreover,the pipes are more easily blocked by protruding particles when the gap becomes narrower.For these reasons,the sand grains in a thicker sample may exhibit higher mobility.

    5 Conclusions

    This study presented a series of experiments in which the dynamic progression of BEP was visualized by using a classic Hele-Shaw cell.The following con‐clusions may be drawn from the test observations.

    1.The process of BEP observed in a uniform,po‐rous medium driven by a sudden flow of water can be divided into a piping progression phase and a piping stabilization phase.The spatial evolution of the ero‐sion pipe was essentially accomplished in the piping progression phase;nonetheless,transport of sand grains within the formed erosion pipes could still be observed in the piping stabilization phase.

    2.The morphologies of erosion under different fluxes of water indicated that a higher flux of water may produce more branched patterns when the gap was not too wide (in this study,it is observed in the tests withbof 0.5 mm).As expected,both the erosion area and the maximum erosion radius were larger when a higher flux of water was employed;however,their growth with the flux of water was not linear:they increased very little when the flux was relatively large.In contrast,the hydraulic pressure measured over the sample showed an almost linear increase with the flux of water.

    3.The BEP observed in the Hele-Shaw cell was not a plane-strain problem.The thickness of the sam‐ple may have a significant influence on both the initia‐tion and progression of BEP.The test results suggested that the sand grains in a thicker Hele-Shaw model were easier to dislodge.This occurs because the granu‐lar system of a thicker sample has more degrees of freedom,which increases the mobility of the sand grains.

    In practical engineering,the underlying strata below a tunnel invert can be infinitely thick,while the thickness of the sample in this study was very limited.As shown in Figs.11 and 12,it is seen that the effect of the sample thickness was still significant even when gapbwas set to 1.5 mm(approximately 9 times the mean grain diameter).There should exist a critical thickness of sample,above which the erosion geome‐try may not be affected.However,further increasing the thickness of the sample in the presented Hele-Shaw cell may dim the contour of the erosion pipe,making it hard to distinguish the eroded and intact areas.Therefore,it is necessary to develop a new apparatus to determine the“critical thickness”of the sample and the mechanism by which erosion progresses along the vertical direction.Naturally,this is our future work.

    Acknowledgments

    This work is supported by the National Engineering Lab‐oratory for Digital Construction and Evaluation Technology of Urban Rail Transit (No.2021GY01) and the National Natural Science Foundation of China(No.41630641).

    Author contributions

    Gang ZHENG designed the research.Jing-bo TONG and Tian-qi ZHANG processed the corresponding data.Tian-qi ZHANG wrote the first draft of the manuscript.Zi-wu WANG,Xun LI,Ji-qing ZHANG,Chun-yu QI,and Hai-zuo ZHOU helped to organize the manuscript.Yu DIAO revised and edited the final version.

    Conflict of interest

    Gang ZHENG,Jing-bo TONG,Tian-qi ZHANG,Zi-wu WANG,Xun LI,Ji-qing ZHANG,Chun-yu QI,Hai-zuo ZHOU,and Yu DIAO declare that they have no conflict of interest.

    国产福利在线免费观看视频| 动漫黄色视频在线观看| 欧美在线黄色| 久久久水蜜桃国产精品网| 国产精品久久久久久人妻精品电影 | 中文字幕人妻丝袜一区二区| 最近最新中文字幕大全免费视频| 欧美日韩一级在线毛片| 搡老乐熟女国产| 97人妻天天添夜夜摸| 精品人妻一区二区三区麻豆| 亚洲av电影在线进入| 精品熟女少妇八av免费久了| 一级毛片女人18水好多| 妹子高潮喷水视频| 亚洲第一青青草原| 免费看十八禁软件| 免费av中文字幕在线| 色老头精品视频在线观看| 丁香六月天网| 亚洲九九香蕉| 色婷婷av一区二区三区视频| 成人黄色视频免费在线看| 欧美日韩亚洲高清精品| 50天的宝宝边吃奶边哭怎么回事| 1024香蕉在线观看| 亚洲中文av在线| 人成视频在线观看免费观看| tocl精华| 99热网站在线观看| 亚洲国产中文字幕在线视频| 高潮久久久久久久久久久不卡| 精品国内亚洲2022精品成人 | 精品乱码久久久久久99久播| 国产在线免费精品| 高清黄色对白视频在线免费看| 成人手机av| 国产成人免费观看mmmm| 我的亚洲天堂| 午夜91福利影院| 王馨瑶露胸无遮挡在线观看| 亚洲欧美日韩高清在线视频 | 成人三级做爰电影| 久久热在线av| 免费一级毛片在线播放高清视频 | 日韩中文字幕视频在线看片| 国产激情久久老熟女| 99精品欧美一区二区三区四区| 精品高清国产在线一区| 国产免费福利视频在线观看| 精品免费久久久久久久清纯 | 亚洲精品国产av成人精品| 精品国产超薄肉色丝袜足j| 女人高潮潮喷娇喘18禁视频| 久久天堂一区二区三区四区| 女人久久www免费人成看片| 美女脱内裤让男人舔精品视频| 纵有疾风起免费观看全集完整版| 欧美在线一区亚洲| 在线av久久热| 一级毛片电影观看| 又黄又粗又硬又大视频| 久久久久网色| 999久久久精品免费观看国产| 一个人免费看片子| 中文字幕人妻丝袜一区二区| 天堂中文最新版在线下载| 国产黄色免费在线视频| 制服诱惑二区| 两性夫妻黄色片| 大型av网站在线播放| 一区二区av电影网| 国产国语露脸激情在线看| 亚洲人成电影观看| cao死你这个sao货| 在线观看免费日韩欧美大片| 欧美在线一区亚洲| 青草久久国产| 777米奇影视久久| 国产精品亚洲av一区麻豆| 亚洲成国产人片在线观看| 99热网站在线观看| 91老司机精品| 99国产精品99久久久久| 色婷婷久久久亚洲欧美| 老汉色∧v一级毛片| 欧美 日韩 精品 国产| 一边摸一边抽搐一进一出视频| 九色亚洲精品在线播放| 一区二区三区激情视频| 建设人人有责人人尽责人人享有的| 9热在线视频观看99| 成人av一区二区三区在线看 | 国产成人av教育| 成人国语在线视频| 在线观看免费高清a一片| 黄网站色视频无遮挡免费观看| 在线亚洲精品国产二区图片欧美| av在线老鸭窝| 午夜福利一区二区在线看| 久久女婷五月综合色啪小说| 伊人亚洲综合成人网| 久久性视频一级片| 日本猛色少妇xxxxx猛交久久| 亚洲欧洲精品一区二区精品久久久| 人人妻人人添人人爽欧美一区卜| 一本大道久久a久久精品| av天堂在线播放| 日韩视频在线欧美| 人妻一区二区av| 纵有疾风起免费观看全集完整版| 777米奇影视久久| 国产高清国产精品国产三级| 搡老熟女国产l中国老女人| 伊人久久大香线蕉亚洲五| 国产一区二区三区综合在线观看| 狠狠婷婷综合久久久久久88av| 视频在线观看一区二区三区| 亚洲欧美成人综合另类久久久| 日韩中文字幕视频在线看片| 啦啦啦视频在线资源免费观看| 婷婷色av中文字幕| 99久久人妻综合| 欧美午夜高清在线| 黄色视频不卡| 这个男人来自地球电影免费观看| 国产色视频综合| 欧美少妇被猛烈插入视频| 日韩 亚洲 欧美在线| a级毛片黄视频| 国产精品熟女久久久久浪| 麻豆乱淫一区二区| 精品亚洲成a人片在线观看| 深夜精品福利| 国产伦理片在线播放av一区| 少妇裸体淫交视频免费看高清 | 欧美 日韩 精品 国产| 欧美亚洲日本最大视频资源| 无遮挡黄片免费观看| 色94色欧美一区二区| 亚洲伊人久久精品综合| 欧美乱码精品一区二区三区| 在线观看人妻少妇| 久久99一区二区三区| 午夜久久久在线观看| 狂野欧美激情性bbbbbb| h视频一区二区三区| 高清黄色对白视频在线免费看| 欧美精品一区二区免费开放| 我的亚洲天堂| 欧美 亚洲 国产 日韩一| 成人国产一区最新在线观看| 91成年电影在线观看| 中文字幕人妻丝袜制服| 国产欧美亚洲国产| 国产精品 国内视频| 日本一区二区免费在线视频| 亚洲成国产人片在线观看| 国产高清videossex| 超碰成人久久| 国产精品99久久99久久久不卡| 搡老岳熟女国产| 国产在线视频一区二区| 国产精品 国内视频| 精品福利永久在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 啦啦啦视频在线资源免费观看| 国产精品99久久99久久久不卡| 精品视频人人做人人爽| 岛国毛片在线播放| 又大又爽又粗| 成人影院久久| 高清欧美精品videossex| 18禁观看日本| 午夜福利乱码中文字幕| 欧美xxⅹ黑人| 国产区一区二久久| √禁漫天堂资源中文www| 久久青草综合色| 99国产精品一区二区三区| 十八禁高潮呻吟视频| 国产日韩欧美在线精品| 久久久水蜜桃国产精品网| 美女大奶头黄色视频| 桃花免费在线播放| 亚洲色图综合在线观看| 国产精品久久久久久精品电影小说| 国产亚洲一区二区精品| 十八禁高潮呻吟视频| 亚洲欧美色中文字幕在线| 久久国产精品大桥未久av| 一区二区三区精品91| 中文精品一卡2卡3卡4更新| 日日爽夜夜爽网站| 国产精品熟女久久久久浪| 狂野欧美激情性bbbbbb| 少妇精品久久久久久久| 亚洲中文av在线| 男女边摸边吃奶| 一本色道久久久久久精品综合| 午夜久久久在线观看| 亚洲va日本ⅴa欧美va伊人久久 | 欧美日韩黄片免| 考比视频在线观看| 免费女性裸体啪啪无遮挡网站| 精品国产一区二区久久| 午夜91福利影院| 不卡一级毛片| 91精品三级在线观看| 国产极品粉嫩免费观看在线| 欧美大码av| 秋霞在线观看毛片| h视频一区二区三区| 国产精品av久久久久免费| 18禁黄网站禁片午夜丰满| 一区二区三区四区激情视频| 美女视频免费永久观看网站| 80岁老熟妇乱子伦牲交| 又大又爽又粗| 国产人伦9x9x在线观看| 下体分泌物呈黄色| 精品人妻在线不人妻| 成年人黄色毛片网站| 两个人看的免费小视频| 777久久人妻少妇嫩草av网站| 无遮挡黄片免费观看| 美女脱内裤让男人舔精品视频| 国产一区二区在线观看av| 国产免费视频播放在线视频| 老熟妇乱子伦视频在线观看 | 视频区欧美日本亚洲| 欧美精品人与动牲交sv欧美| 午夜激情久久久久久久| 男男h啪啪无遮挡| 超碰97精品在线观看| 欧美午夜高清在线| 嫁个100分男人电影在线观看| 91麻豆av在线| 亚洲欧美日韩高清在线视频 | 欧美性长视频在线观看| av天堂在线播放| av有码第一页| 视频区图区小说| 看免费av毛片| 韩国高清视频一区二区三区| 两个人免费观看高清视频| 啪啪无遮挡十八禁网站| 一级片'在线观看视频| 少妇 在线观看| 国产伦理片在线播放av一区| 日韩大码丰满熟妇| 丰满人妻熟妇乱又伦精品不卡| 侵犯人妻中文字幕一二三四区| 成年动漫av网址| 国产成人影院久久av| 狂野欧美激情性xxxx| 在线观看舔阴道视频| 啦啦啦免费观看视频1| 日韩大码丰满熟妇| 一级毛片电影观看| 黄网站色视频无遮挡免费观看| 黄色视频不卡| 99精品久久久久人妻精品| 丰满少妇做爰视频| 美女午夜性视频免费| 交换朋友夫妻互换小说| 久久久久国产一级毛片高清牌| 老熟妇乱子伦视频在线观看 | 狠狠婷婷综合久久久久久88av| 少妇人妻久久综合中文| 咕卡用的链子| av视频免费观看在线观看| 久久久久久久久久久久大奶| 国产1区2区3区精品| 纯流量卡能插随身wifi吗| 19禁男女啪啪无遮挡网站| 这个男人来自地球电影免费观看| 淫妇啪啪啪对白视频 | 国产成人精品在线电影| 亚洲伊人久久精品综合| a 毛片基地| 午夜免费观看性视频| 国产成人影院久久av| av网站在线播放免费| 国产视频一区二区在线看| 精品久久久久久久毛片微露脸 | 王馨瑶露胸无遮挡在线观看| 飞空精品影院首页| www.999成人在线观看| 久9热在线精品视频| 精品少妇久久久久久888优播| 下体分泌物呈黄色| 久热这里只有精品99| 国产极品粉嫩免费观看在线| 亚洲va日本ⅴa欧美va伊人久久 | 国产亚洲精品久久久久5区| 国产欧美日韩综合在线一区二区| 丁香六月欧美| 国产亚洲av片在线观看秒播厂| 亚洲精品久久久久久婷婷小说| 丁香六月欧美| www.自偷自拍.com| 日本黄色日本黄色录像| 欧美日韩一级在线毛片| 人妻人人澡人人爽人人| 国产日韩欧美亚洲二区| 久久99一区二区三区| 欧美中文综合在线视频| av免费在线观看网站| 亚洲av美国av| 69精品国产乱码久久久| 夜夜骑夜夜射夜夜干| 99精品欧美一区二区三区四区| 一进一出抽搐动态| 超碰成人久久| 99re6热这里在线精品视频| 日韩制服骚丝袜av| 女性被躁到高潮视频| 午夜免费观看性视频| 在线观看免费日韩欧美大片| 日本av手机在线免费观看| 久久久国产欧美日韩av| 亚洲国产欧美一区二区综合| 色视频在线一区二区三区| 成年av动漫网址| 国产福利在线免费观看视频| 黄片播放在线免费| 男女午夜视频在线观看| 欧美日韩成人在线一区二区| 99九九在线精品视频| 免费少妇av软件| 波多野结衣一区麻豆| 午夜福利视频精品| 午夜福利视频在线观看免费| 亚洲人成电影观看| 一区二区三区乱码不卡18| 男女无遮挡免费网站观看| 免费看十八禁软件| 日日夜夜操网爽| 丁香六月天网| 桃红色精品国产亚洲av| 真人做人爱边吃奶动态| 搡老岳熟女国产| 国产成人影院久久av| 91精品国产国语对白视频| 美女中出高潮动态图| 最黄视频免费看| 色老头精品视频在线观看| 国产在线视频一区二区| 一本色道久久久久久精品综合| 丝袜脚勾引网站| 精品人妻熟女毛片av久久网站| 两人在一起打扑克的视频| 亚洲精品一卡2卡三卡4卡5卡 | 桃红色精品国产亚洲av| 色老头精品视频在线观看| 免费高清在线观看日韩| 法律面前人人平等表现在哪些方面 | 少妇 在线观看| 免费高清在线观看日韩| 多毛熟女@视频| 999久久久精品免费观看国产| 19禁男女啪啪无遮挡网站| 中文欧美无线码| 国产免费视频播放在线视频| 新久久久久国产一级毛片| 狠狠婷婷综合久久久久久88av| av天堂久久9| 十八禁高潮呻吟视频| 国产精品秋霞免费鲁丝片| 999久久久国产精品视频| 精品国产乱码久久久久久小说| 国产欧美日韩一区二区三区在线| 亚洲国产av影院在线观看| 人人妻,人人澡人人爽秒播| 极品人妻少妇av视频| 人人妻,人人澡人人爽秒播| 亚洲欧美成人综合另类久久久| 成人国产一区最新在线观看| 成人国产av品久久久| 国产精品香港三级国产av潘金莲| 久久久国产一区二区| 午夜福利乱码中文字幕| 国产精品一区二区在线不卡| 在线精品无人区一区二区三| 美女福利国产在线| 狠狠婷婷综合久久久久久88av| av不卡在线播放| 男女下面插进去视频免费观看| a级毛片黄视频| 母亲3免费完整高清在线观看| 狂野欧美激情性bbbbbb| 两性夫妻黄色片| 少妇粗大呻吟视频| 精品熟女少妇八av免费久了| 亚洲少妇的诱惑av| 午夜激情久久久久久久| 国产不卡av网站在线观看| 久久午夜综合久久蜜桃| 亚洲国产日韩一区二区| 国产成人精品久久二区二区免费| 久久人妻福利社区极品人妻图片| www日本在线高清视频| 最近最新中文字幕大全免费视频| 97精品久久久久久久久久精品| 久久久久国内视频| 12—13女人毛片做爰片一| bbb黄色大片| 女人高潮潮喷娇喘18禁视频| 成年女人毛片免费观看观看9 | 777久久人妻少妇嫩草av网站| 麻豆av在线久日| 天堂俺去俺来也www色官网| 男人操女人黄网站| 黑人巨大精品欧美一区二区蜜桃| 少妇的丰满在线观看| 一个人免费看片子| 黑人猛操日本美女一级片| 黄片小视频在线播放| 人人妻人人添人人爽欧美一区卜| 午夜福利免费观看在线| 国产精品久久久久成人av| 欧美精品一区二区大全| 丰满饥渴人妻一区二区三| 99香蕉大伊视频| 男人舔女人的私密视频| 国产精品亚洲av一区麻豆| 青春草视频在线免费观看| 久久久久精品国产欧美久久久 | 国产亚洲av高清不卡| 国产一区二区三区在线臀色熟女 | 亚洲精品日韩在线中文字幕| 国产一级毛片在线| av片东京热男人的天堂| 成年动漫av网址| 伦理电影免费视频| 岛国毛片在线播放| 老司机深夜福利视频在线观看 | 在线亚洲精品国产二区图片欧美| 成人国产av品久久久| 天天添夜夜摸| 下体分泌物呈黄色| tocl精华| 国产亚洲av片在线观看秒播厂| 国产又色又爽无遮挡免| 国产一区二区 视频在线| www.av在线官网国产| 在线观看舔阴道视频| 亚洲国产精品一区二区三区在线| 精品少妇一区二区三区视频日本电影| 后天国语完整版免费观看| 亚洲全国av大片| 首页视频小说图片口味搜索| 国产又色又爽无遮挡免| 国产免费视频播放在线视频| 亚洲欧美一区二区三区黑人| 精品少妇一区二区三区视频日本电影| 欧美日韩福利视频一区二区| 久久久国产欧美日韩av| 正在播放国产对白刺激| 国产精品.久久久| 中文字幕精品免费在线观看视频| 中亚洲国语对白在线视频| 婷婷色av中文字幕| 久久久久久亚洲精品国产蜜桃av| 日韩欧美一区二区三区在线观看 | 久久精品久久久久久噜噜老黄| 久久久国产欧美日韩av| 中文字幕最新亚洲高清| 90打野战视频偷拍视频| 精品福利观看| 精品熟女少妇八av免费久了| 日本欧美视频一区| 91老司机精品| 国产精品av久久久久免费| 777久久人妻少妇嫩草av网站| 制服诱惑二区| 波多野结衣av一区二区av| 黑人猛操日本美女一级片| 亚洲男人天堂网一区| 热re99久久国产66热| 51午夜福利影视在线观看| 亚洲五月色婷婷综合| 久久久国产成人免费| tube8黄色片| 国产精品一区二区在线不卡| 亚洲人成77777在线视频| 男人爽女人下面视频在线观看| 青草久久国产| 精品亚洲成a人片在线观看| 丰满人妻熟妇乱又伦精品不卡| 在线天堂中文资源库| 又黄又粗又硬又大视频| 国产黄频视频在线观看| 97在线人人人人妻| 亚洲精品国产色婷婷电影| 国产精品.久久久| 男人爽女人下面视频在线观看| 制服诱惑二区| 在线观看一区二区三区激情| 两个人看的免费小视频| 久久人妻福利社区极品人妻图片| 亚洲欧美激情在线| 国产成人免费观看mmmm| 亚洲精品一二三| 午夜福利视频精品| 电影成人av| 12—13女人毛片做爰片一| 精品欧美一区二区三区在线| 99re6热这里在线精品视频| 久久99热这里只频精品6学生| 十八禁网站免费在线| 中文欧美无线码| 亚洲欧洲精品一区二区精品久久久| 国产又色又爽无遮挡免| 久久久久精品国产欧美久久久 | 国产成人精品在线电影| 一级毛片电影观看| 另类亚洲欧美激情| 50天的宝宝边吃奶边哭怎么回事| 青春草视频在线免费观看| 美女扒开内裤让男人捅视频| av在线老鸭窝| 亚洲精品乱久久久久久| 国产精品影院久久| 中国美女看黄片| 国产亚洲av高清不卡| 中文字幕人妻熟女乱码| 亚洲精品国产av蜜桃| 日本一区二区免费在线视频| 国产欧美日韩一区二区三区在线| 久久久精品国产亚洲av高清涩受| 三上悠亚av全集在线观看| 久久精品久久久久久噜噜老黄| 91字幕亚洲| 亚洲人成电影免费在线| 成人av一区二区三区在线看 | 下体分泌物呈黄色| 国产色视频综合| 日韩制服丝袜自拍偷拍| 午夜老司机福利片| 亚洲全国av大片| 亚洲国产av影院在线观看| 涩涩av久久男人的天堂| 老司机午夜福利在线观看视频 | 热re99久久精品国产66热6| 亚洲熟女精品中文字幕| 18禁观看日本| 久久久久国产一级毛片高清牌| 国产又色又爽无遮挡免| 久久热在线av| 人人妻人人澡人人爽人人夜夜| 一本综合久久免费| 亚洲av欧美aⅴ国产| 天堂俺去俺来也www色官网| 国产欧美日韩一区二区精品| 天天躁夜夜躁狠狠躁躁| 夜夜夜夜夜久久久久| 免费黄频网站在线观看国产| 国产男女内射视频| 老汉色av国产亚洲站长工具| 久久这里只有精品19| 欧美精品高潮呻吟av久久| 亚洲avbb在线观看| 久久精品国产亚洲av香蕉五月 | 久久精品国产亚洲av高清一级| 午夜激情av网站| 午夜福利在线观看吧| 国产高清视频在线播放一区 | 老熟女久久久| 一本色道久久久久久精品综合| 精品国产一区二区三区四区第35| 制服人妻中文乱码| 黄色视频不卡| 99久久99久久久精品蜜桃| 精品福利永久在线观看| 蜜桃国产av成人99| 视频区图区小说| 波多野结衣av一区二区av| 国产激情久久老熟女| 亚洲第一青青草原| 啦啦啦在线免费观看视频4| 最近中文字幕2019免费版| 国产精品1区2区在线观看. | 亚洲国产精品成人久久小说| 亚洲国产欧美一区二区综合| 三上悠亚av全集在线观看| 桃花免费在线播放| 一区福利在线观看| 久久影院123| 久久中文看片网| 真人做人爱边吃奶动态| 国产精品.久久久| 9热在线视频观看99| 五月开心婷婷网| av电影中文网址| 中亚洲国语对白在线视频| 精品国产超薄肉色丝袜足j| 亚洲精品中文字幕一二三四区 | 另类精品久久| 日本vs欧美在线观看视频| 精品人妻在线不人妻| 亚洲成人手机| 日本a在线网址| 侵犯人妻中文字幕一二三四区| 色婷婷av一区二区三区视频| 日韩三级视频一区二区三区| 久久国产精品人妻蜜桃| 免费少妇av软件| 50天的宝宝边吃奶边哭怎么回事| 成在线人永久免费视频| 国产精品一二三区在线看| 老司机福利观看| 老汉色av国产亚洲站长工具| 黄色视频在线播放观看不卡|