• <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.

    久久99蜜桃精品久久| 人人妻人人看人人澡| 秋霞在线观看毛片| 免费观看的影片在线观看| 欧美xxxx黑人xx丫x性爽| 日韩av在线大香蕉| 超碰av人人做人人爽久久| 国产精品精品国产色婷婷| 中国美白少妇内射xxxbb| 看非洲黑人一级黄片| 免费看a级黄色片| 亚洲中文字幕日韩| 成人高潮视频无遮挡免费网站| 国产色爽女视频免费观看| 日本与韩国留学比较| 中文字幕免费在线视频6| 爱豆传媒免费全集在线观看| 国产亚洲精品久久久久久毛片| 可以在线观看的亚洲视频| 久久久久久伊人网av| 伦精品一区二区三区| 国产精华一区二区三区| 国产极品天堂在线| 午夜福利成人在线免费观看| 亚洲无线在线观看| 综合色av麻豆| 我要看日韩黄色一级片| 一级毛片我不卡| 精品欧美国产一区二区三| 成年免费大片在线观看| 级片在线观看| 最近最新中文字幕大全电影3| 别揉我奶头 嗯啊视频| 亚洲国产精品久久男人天堂| 如何舔出高潮| 夜夜夜夜夜久久久久| 国产成人精品婷婷| 在线观看66精品国产| 尤物成人国产欧美一区二区三区| 国产精品女同一区二区软件| 国产精品无大码| 在线国产一区二区在线| 中国美白少妇内射xxxbb| 日本一本二区三区精品| 丰满人妻一区二区三区视频av| 在线观看一区二区三区| 色噜噜av男人的天堂激情| 国产精品国产三级国产av玫瑰| 一级毛片电影观看 | 亚洲av中文av极速乱| av免费观看日本| 99久久精品热视频| 欧美一区二区国产精品久久精品| 超碰av人人做人人爽久久| 亚洲高清免费不卡视频| 成人av在线播放网站| 亚洲精品456在线播放app| 丝袜喷水一区| 国产探花极品一区二区| 亚洲成人中文字幕在线播放| 国模一区二区三区四区视频| 国内久久婷婷六月综合欲色啪| 黄色日韩在线| www.色视频.com| 亚洲高清免费不卡视频| 黄色欧美视频在线观看| 久久午夜亚洲精品久久| 日本一本二区三区精品| 欧美+亚洲+日韩+国产| 成年女人永久免费观看视频| 成人二区视频| 成人高潮视频无遮挡免费网站| 精品久久久久久久久亚洲| а√天堂www在线а√下载| 又爽又黄无遮挡网站| 熟妇人妻久久中文字幕3abv| 日韩中字成人| 能在线免费看毛片的网站| 中文在线观看免费www的网站| 插逼视频在线观看| 亚洲精品粉嫩美女一区| 亚洲av免费高清在线观看| 国内精品久久久久精免费| 身体一侧抽搐| 亚洲精品粉嫩美女一区| 女人被狂操c到高潮| 国产精品乱码一区二三区的特点| 久久久久久久久大av| 日韩高清综合在线| 可以在线观看毛片的网站| 成人毛片60女人毛片免费| 国产亚洲91精品色在线| 美女高潮的动态| 欧美+日韩+精品| 亚洲欧美日韩无卡精品| 亚洲婷婷狠狠爱综合网| 免费观看精品视频网站| 亚洲欧美日韩无卡精品| 高清毛片免费观看视频网站| 日韩视频在线欧美| 国产精品无大码| 日本av手机在线免费观看| 国产一级毛片在线| 九九在线视频观看精品| 国产一区二区激情短视频| 少妇人妻精品综合一区二区 | 亚洲精品色激情综合| 1024手机看黄色片| 国产成人一区二区在线| 久久久成人免费电影| 最后的刺客免费高清国语| 亚洲av成人精品一区久久| 欧美精品一区二区大全| 日日干狠狠操夜夜爽| 99久久人妻综合| 99视频精品全部免费 在线| 国产精品国产三级国产av玫瑰| 日韩亚洲欧美综合| 日日摸夜夜添夜夜添av毛片| 搡女人真爽免费视频火全软件| 国模一区二区三区四区视频| 赤兔流量卡办理| 亚洲精品久久久久久婷婷小说 | 99久久无色码亚洲精品果冻| 天堂影院成人在线观看| 日韩欧美 国产精品| 一级毛片电影观看 | 天堂影院成人在线观看| 久久精品人妻少妇| 毛片一级片免费看久久久久| 青春草视频在线免费观看| 12—13女人毛片做爰片一| 日韩强制内射视频| 99热这里只有精品一区| 欧美精品一区二区大全| 亚洲欧美日韩卡通动漫| 免费人成在线观看视频色| 国产三级中文精品| 欧美在线一区亚洲| 成人综合一区亚洲| 国产一区二区三区av在线 | 国产高潮美女av| 嫩草影院精品99| 如何舔出高潮| 人妻夜夜爽99麻豆av| 久久久精品94久久精品| 久久精品国产亚洲av天美| 亚洲电影在线观看av| 久久精品综合一区二区三区| 午夜精品国产一区二区电影 | 日本爱情动作片www.在线观看| 精品99又大又爽又粗少妇毛片| 夜夜看夜夜爽夜夜摸| 国产一区二区激情短视频| 久久久国产成人精品二区| 免费无遮挡裸体视频| 国产精品免费一区二区三区在线| 99热全是精品| 小说图片视频综合网站| 国产一区二区三区在线臀色熟女| 国产久久久一区二区三区| 精品一区二区三区人妻视频| 日本免费a在线| 又爽又黄a免费视频| 国产伦在线观看视频一区| 此物有八面人人有两片| 国产成人a∨麻豆精品| 99热这里只有是精品在线观看| 亚洲性久久影院| 久久亚洲精品不卡| 午夜激情欧美在线| 99视频精品全部免费 在线| 国产成人福利小说| 中出人妻视频一区二区| 边亲边吃奶的免费视频| 欧美在线一区亚洲| 亚洲人成网站在线观看播放| 国产精品av视频在线免费观看| 亚洲四区av| 亚洲无线观看免费| 亚洲熟妇中文字幕五十中出| 欧美成人精品欧美一级黄| 夜夜看夜夜爽夜夜摸| or卡值多少钱| 亚洲人成网站高清观看| 少妇被粗大猛烈的视频| 午夜精品在线福利| 国内精品一区二区在线观看| 欧美性猛交黑人性爽| 人妻久久中文字幕网| 国内精品美女久久久久久| 少妇裸体淫交视频免费看高清| 天堂影院成人在线观看| 成人三级黄色视频| 亚洲人与动物交配视频| 亚洲欧美日韩东京热| 精品午夜福利在线看| 热99re8久久精品国产| 欧美不卡视频在线免费观看| 综合色av麻豆| 国产成人a区在线观看| 亚洲乱码一区二区免费版| 99riav亚洲国产免费| 国产成人a区在线观看| 中文欧美无线码| 一级毛片我不卡| 久久久欧美国产精品| 黄片wwwwww| 久久午夜亚洲精品久久| 人人妻人人澡人人爽人人夜夜 | 精品人妻偷拍中文字幕| 最后的刺客免费高清国语| 日韩欧美在线乱码| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品久久久久久婷婷小说 | 可以在线观看的亚洲视频| 黄色配什么色好看| 国产一区二区三区在线臀色熟女| 国产日韩欧美在线精品| 日韩精品有码人妻一区| av卡一久久| 久久精品国产亚洲av天美| 天天一区二区日本电影三级| 国产精品福利在线免费观看| 久久久久九九精品影院| 国产精品乱码一区二三区的特点| 中文精品一卡2卡3卡4更新| 久久久国产成人精品二区| 亚洲欧美精品自产自拍| 18禁在线无遮挡免费观看视频| 麻豆一二三区av精品| 精品国产三级普通话版| 国产黄色视频一区二区在线观看 | 26uuu在线亚洲综合色| 麻豆国产av国片精品| 国产精品乱码一区二三区的特点| 91久久精品国产一区二区三区| 亚洲精品乱码久久久v下载方式| 91午夜精品亚洲一区二区三区| 亚洲高清免费不卡视频| 国产精品乱码一区二三区的特点| 亚洲最大成人av| 直男gayav资源| 男人舔奶头视频| 一本久久精品| 国产白丝娇喘喷水9色精品| 欧美性猛交黑人性爽| 啦啦啦啦在线视频资源| 成人永久免费在线观看视频| 狠狠狠狠99中文字幕| 午夜福利视频1000在线观看| 亚洲丝袜综合中文字幕| 中文欧美无线码| 日韩av在线大香蕉| 亚洲图色成人| 三级经典国产精品| 99热只有精品国产| 一进一出抽搐动态| 麻豆成人av视频| 男女下面进入的视频免费午夜| 91在线精品国自产拍蜜月| 亚洲国产精品国产精品| 国产极品精品免费视频能看的| 国产精品一区二区三区四区久久| 简卡轻食公司| 日韩大尺度精品在线看网址| 欧美+亚洲+日韩+国产| 我要搜黄色片| 内地一区二区视频在线| 国产探花极品一区二区| 草草在线视频免费看| 给我免费播放毛片高清在线观看| 赤兔流量卡办理| 亚洲成人中文字幕在线播放| 夜夜爽天天搞| 尤物成人国产欧美一区二区三区| 亚洲成人av在线免费| 日日撸夜夜添| 精品一区二区三区视频在线| 赤兔流量卡办理| 亚洲国产精品sss在线观看| 免费搜索国产男女视频| 99热这里只有是精品50| 蜜臀久久99精品久久宅男| 午夜亚洲福利在线播放| 蜜桃久久精品国产亚洲av| 国产熟女欧美一区二区| 久久精品影院6| 韩国av在线不卡| 黄色一级大片看看| 国产精品国产三级国产av玫瑰| 成人av在线播放网站| 亚洲无线观看免费| 国产高清激情床上av| 一本一本综合久久| 日日摸夜夜添夜夜添av毛片| 国产伦精品一区二区三区视频9| 国产精品嫩草影院av在线观看| 国产午夜精品论理片| 五月伊人婷婷丁香| 免费看日本二区| 波多野结衣高清无吗| 禁无遮挡网站| av.在线天堂| av卡一久久| 国产精品乱码一区二三区的特点| 天美传媒精品一区二区| 久久久久久久久久成人| 久久久精品欧美日韩精品| 免费人成在线观看视频色| 久久中文看片网| 欧美极品一区二区三区四区| 大香蕉久久网| 国语自产精品视频在线第100页| 国模一区二区三区四区视频| 精品欧美国产一区二区三| 久久这里有精品视频免费| 日日啪夜夜撸| 岛国在线免费视频观看| 国产亚洲av嫩草精品影院| 男女那种视频在线观看| 午夜福利在线观看免费完整高清在 | 国产av不卡久久| 美女高潮的动态| 日本成人三级电影网站| 色哟哟·www| 一级av片app| 亚洲成人av在线免费| 成年女人永久免费观看视频| 中出人妻视频一区二区| 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品亚洲一区二区| 国产成人a∨麻豆精品| 久久这里有精品视频免费| av国产免费在线观看| 亚洲欧美清纯卡通| 久久久久性生活片| 亚洲av第一区精品v没综合| 最近的中文字幕免费完整| 亚洲久久久久久中文字幕| 免费一级毛片在线播放高清视频| 国产一级毛片七仙女欲春2| 日韩制服骚丝袜av| 久99久视频精品免费| 久久久久久久亚洲中文字幕| 国产黄片美女视频| 哪个播放器可以免费观看大片| 国产国拍精品亚洲av在线观看| 一边摸一边抽搐一进一小说| 精品无人区乱码1区二区| 一个人看视频在线观看www免费| 男女下面进入的视频免费午夜| 好男人视频免费观看在线| 亚洲在线自拍视频| 亚洲一区二区三区色噜噜| 久久精品国产亚洲av香蕉五月| 久久九九热精品免费| 99久久成人亚洲精品观看| 亚洲成av人片在线播放无| 深夜a级毛片| 一夜夜www| 欧美成人精品欧美一级黄| 免费无遮挡裸体视频| avwww免费| 免费观看的影片在线观看| 精品国内亚洲2022精品成人| av天堂在线播放| 两个人视频免费观看高清| 亚洲国产精品久久男人天堂| 1000部很黄的大片| 亚洲av男天堂| 99国产精品一区二区蜜桃av| 一本一本综合久久| 色尼玛亚洲综合影院| 亚洲三级黄色毛片| 黄色配什么色好看| 波野结衣二区三区在线| 少妇被粗大猛烈的视频| 麻豆成人午夜福利视频| 精品久久久久久久久久久久久| 中文资源天堂在线| 国产老妇女一区| 国产av在哪里看| 99riav亚洲国产免费| 老师上课跳d突然被开到最大视频| 精品少妇黑人巨大在线播放 | 日日啪夜夜撸| 18禁裸乳无遮挡免费网站照片| 99久久九九国产精品国产免费| 少妇裸体淫交视频免费看高清| 国产欧美日韩精品一区二区| 亚洲美女视频黄频| 欧美3d第一页| 欧美精品一区二区大全| 18禁裸乳无遮挡免费网站照片| 不卡视频在线观看欧美| 亚洲精品粉嫩美女一区| 久久精品国产亚洲av天美| 久久久午夜欧美精品| 中国美女看黄片| 亚洲性久久影院| 免费不卡的大黄色大毛片视频在线观看 | 国内精品宾馆在线| 中文资源天堂在线| 天天躁夜夜躁狠狠久久av| 美女黄网站色视频| 久久草成人影院| 久久热精品热| 国产亚洲av嫩草精品影院| 看十八女毛片水多多多| 一级二级三级毛片免费看| 99热精品在线国产| 91麻豆精品激情在线观看国产| 国产精品蜜桃在线观看 | 一夜夜www| 久久久欧美国产精品| 国产亚洲91精品色在线| 久久久久久久久久久丰满| 深夜a级毛片| 少妇熟女欧美另类| 三级国产精品欧美在线观看| 性色avwww在线观看| 一级毛片aaaaaa免费看小| 3wmmmm亚洲av在线观看| 韩国av在线不卡| 欧美三级亚洲精品| 观看美女的网站| 日韩一区二区视频免费看| 久久精品国产亚洲av香蕉五月| 干丝袜人妻中文字幕| 日韩一本色道免费dvd| 一本久久中文字幕| 少妇熟女aⅴ在线视频| 一级毛片电影观看 | 久久99热6这里只有精品| 边亲边吃奶的免费视频| 18禁裸乳无遮挡免费网站照片| 九草在线视频观看| av在线天堂中文字幕| 成人一区二区视频在线观看| 我的女老师完整版在线观看| 精品99又大又爽又粗少妇毛片| 欧美日本亚洲视频在线播放| 可以在线观看的亚洲视频| 天堂中文最新版在线下载 | 国内精品一区二区在线观看| 欧美日韩乱码在线| 国产精品1区2区在线观看.| 禁无遮挡网站| 国产精品久久久久久亚洲av鲁大| 日韩三级伦理在线观看| 午夜免费男女啪啪视频观看| 一区二区三区四区激情视频 | 久99久视频精品免费| 亚洲真实伦在线观看| 久久草成人影院| 久久99热这里只有精品18| 精品久久久噜噜| 国产精品爽爽va在线观看网站| 一个人看视频在线观看www免费| 欧美成人a在线观看| 日韩亚洲欧美综合| 国产午夜福利久久久久久| 最近2019中文字幕mv第一页| 亚洲精品日韩av片在线观看| 国产综合懂色| 九九热线精品视视频播放| 国产成人freesex在线| 在线观看午夜福利视频| 在线观看av片永久免费下载| 国产精品一区www在线观看| 成人毛片60女人毛片免费| 亚洲国产色片| 91久久精品电影网| 99久久精品国产国产毛片| 国产黄片美女视频| 一级二级三级毛片免费看| 国产午夜精品一二区理论片| 一区福利在线观看| 午夜爱爱视频在线播放| 欧美日本视频| 精品一区二区免费观看| 国产精品嫩草影院av在线观看| 国产伦理片在线播放av一区 | 激情 狠狠 欧美| 精品99又大又爽又粗少妇毛片| av在线蜜桃| 一级毛片久久久久久久久女| 精品久久久久久久久久久久久| 国产亚洲5aaaaa淫片| 午夜老司机福利剧场| av在线亚洲专区| 免费不卡的大黄色大毛片视频在线观看 | 国产高清不卡午夜福利| 成年女人永久免费观看视频| 欧美成人免费av一区二区三区| 国语自产精品视频在线第100页| 成人特级av手机在线观看| 久久精品国产亚洲av天美| 午夜福利在线在线| 日日摸夜夜添夜夜爱| 大又大粗又爽又黄少妇毛片口| 3wmmmm亚洲av在线观看| 亚洲人与动物交配视频| 国产亚洲精品久久久久久毛片| 噜噜噜噜噜久久久久久91| 国产大屁股一区二区在线视频| 精品一区二区三区人妻视频| 女人十人毛片免费观看3o分钟| 99精品在免费线老司机午夜| 亚洲国产精品sss在线观看| 亚洲最大成人手机在线| www.av在线官网国产| 亚洲精品国产成人久久av| 国产高清三级在线| 欧美日本亚洲视频在线播放| 国产人妻一区二区三区在| 在线a可以看的网站| 少妇裸体淫交视频免费看高清| 99精品在免费线老司机午夜| 精品一区二区免费观看| 又爽又黄无遮挡网站| 亚洲精品日韩av片在线观看| 最近2019中文字幕mv第一页| 在线观看免费视频日本深夜| 免费人成视频x8x8入口观看| 亚洲丝袜综合中文字幕| 床上黄色一级片| 赤兔流量卡办理| 欧美日韩乱码在线| 亚洲最大成人中文| 亚洲婷婷狠狠爱综合网| 亚洲中文字幕日韩| 成人性生交大片免费视频hd| 成人毛片a级毛片在线播放| 高清毛片免费观看视频网站| 精品熟女少妇av免费看| 久久九九热精品免费| 国产熟女欧美一区二区| 久久午夜福利片| 免费无遮挡裸体视频| 丝袜美腿在线中文| 99久久九九国产精品国产免费| 亚洲最大成人av| 欧美最黄视频在线播放免费| 国产精品久久久久久精品电影| 一本精品99久久精品77| 亚洲国产欧洲综合997久久,| 成人综合一区亚洲| 国产视频首页在线观看| 亚洲无线观看免费| 久久6这里有精品| 日本成人三级电影网站| 好男人在线观看高清免费视频| 中文字幕av成人在线电影| 国产69精品久久久久777片| 少妇丰满av| av在线老鸭窝| 麻豆国产97在线/欧美| 91精品一卡2卡3卡4卡| 日日干狠狠操夜夜爽| 一级黄色大片毛片| 国产老妇伦熟女老妇高清| 此物有八面人人有两片| 中出人妻视频一区二区| 日韩欧美一区二区三区在线观看| 久久精品国产鲁丝片午夜精品| 美女大奶头视频| 别揉我奶头 嗯啊视频| 久久中文看片网| 久久久精品大字幕| 精品久久国产蜜桃| 欧美一级a爱片免费观看看| 精品久久久久久久久av| 中国美白少妇内射xxxbb| 在线观看美女被高潮喷水网站| 国产精品国产高清国产av| 欧美xxxx黑人xx丫x性爽| 国产午夜福利久久久久久| 国产美女午夜福利| 国产亚洲5aaaaa淫片| 中文字幕制服av| 亚洲av男天堂| 欧美潮喷喷水| 综合色丁香网| 三级男女做爰猛烈吃奶摸视频| 99久久精品国产国产毛片| 12—13女人毛片做爰片一| 欧美色欧美亚洲另类二区| 日韩中字成人| 老女人水多毛片| a级毛片免费高清观看在线播放| 搞女人的毛片| 18禁在线无遮挡免费观看视频| 国产精品一区二区性色av| 最近2019中文字幕mv第一页| 亚洲精品自拍成人| 日本一本二区三区精品| 最近中文字幕高清免费大全6| 身体一侧抽搐| 别揉我奶头 嗯啊视频| 深爱激情五月婷婷| 99久久久亚洲精品蜜臀av| 午夜福利成人在线免费观看| 色哟哟·www| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久久久免| 色综合站精品国产| 久久鲁丝午夜福利片| 禁无遮挡网站| 2021天堂中文幕一二区在线观| 人妻久久中文字幕网| 小说图片视频综合网站| 久久精品综合一区二区三区| 高清在线视频一区二区三区 |