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

    Three-dimensional clogging structures of granular spheres near hopper orifice

    2022-01-23 06:36:44JingYang楊敬DianjinfengGong宮殿錦豐XiaoxueWang汪曉雪ZhichaoWang王志超JianqiLi李建奇BingwenHu胡炳文andChengjieXia夏成杰
    Chinese Physics B 2022年1期
    關(guān)鍵詞:宮殿

    Jing Yang(楊敬), Dianjinfeng Gong(宮殿錦豐), Xiaoxue Wang(汪曉雪), Zhichao Wang(王志超),Jianqi Li(李建奇), Bingwen Hu(胡炳文), and Chengjie Xia(夏成杰)

    Shanghai Key Laboratory of Magnetic Resonance,School of Physics and Electronic Science,East China Normal University,Shanghai 200241,China

    Keywords: granular materials,clogging,magnetic resonance imaging

    1. Introduction

    Granular materials are ubiquitous in industrial applications and daily life, which may exhibit behaviors of both complex fluids and marginal solids under different external conditions.[1,2]The commonly-observed clogging phenomenon of granular particles is essentially a spontaneous jamming transition from a flowing state to a mechanicallystable static state,[3-5]which have aroused great interest from both industry and physics community.Granular material clogs when flowing through a narrow outlet in a statistical way.[6]Clogging and related jamming phenomena are ubiquitous not only in granular systems but also in other systems, including colloidal suspensions,[7,8]self-driven active materials,[9-11]and microbes,[12]etc. A flow of granular particles clogs due to the formation of a clogging structure near the orifice with some typical geometry that sustains its overhead loading. The probability of forming such a clogging structure depends on the hopper geometry,[13-15]properties of granular particles including shapes, frictions, and polydispersity,[16,17]as well as other external control parameters affecting the kinematics of granular flow.[18-22]In a two-dimensional (2D) model system of hard disks, an arch-like clogging structure has been observed, and its formation is mapped onto a restricted selfavoiding random walk to estimate the clogging probability using a pure geometry-based model.[23,24]Furthermore, nonpermanent clogging under external perturbations has been attributed to the existence of “defects” in arch structures,[25]and the associated non-ergodic unclogging dynamics has been modeled as a random exploration out of a potential basin in the landscape of arch configurations.[26-29]All these previous works emphasize the importance of the geometry of clogging structures, but whether the established concepts in these probabilistic models are independent of the spatial dimension remains untested. For three-dimensional (3D) granular materials, the lack of knowledge about the statistical geometric properties of the clogging structure renders a quantitative modeling of their clogging transition and discharging process difficult. Specifically, although it is reasonable to presume a dome-like 3D clogging structure by analogy with its 2D counterpart, the geometry of a 3D clogging structure, which is supposed to be a complex network, is much more complicated to be characterized than the one-dimensional(1D)arch in 2D case. Therefore,a systematic characterization of the arrangements of granular particles blocking the orifice based on various imaging experiments is necessary,[30,31]which should form a structural foundation to validate and generalize current models,and even to help understanding other non-granular 3D systems in the light of a general clogging phase diagram.[32]

    In this work, dome-like clogging structures are reconstructed and identified with the help of magnetic resonance imaging (MRI) techniques, and their geometric and topological properties are quantified. It turns out that 3D clogging structures are distributed on a slightly flattened hemispherical surface, and are majorly composed of chain-like linear structures associated with a certain correlation length, which clearly manifest characteristics of mutual-stabilization and load-bearing. These observations allow us to form a simple phenomenological model to express the average size of clogging structures as a function of outlet size and shape,which is crucial in predicting clogging probability and avalanche sizes during silo discharging.[33-36]

    The current paper is organized as follows. In Section 2,we introduce the clogging experiment and image processing method used to reconstruct the packing. In Section 3, we demonstrate various structural characterizations of the clogging structures. In Section 4, we provide a simple model for the clogging structure sizes. Concluding remarks are given in Section 5.

    2. Materials and experiments

    Granular particles used in this work are hollow plastic spheres filled with edible gelatin. The magnetic resonance signals are contributed mostly by protons in water molecules of the gelatin, so that the internal part of every sphere appears bright in magnetic resonance images, which enables a straightforward labelling of all spheres with image processing methods (Fig. 1(a)). The average diameter of these spheres isd= 15.78±0.07 mm. Hoppers with a zero slop angleαslope=0°and eight different orifice diametersDoutare used:fromDout=31.26 mm≈2dtoDout=59.09 mm≈3.75dat an interval of every 0.25d,as well as another three hoppers with a fixedDout=3dand three different slop angles:αslope=30°,45°, and 60°. The hoppers are made of polyvinyl chloride.The friction coefficient between a spheres is measured to be 0.43±0.02 (with a method introduced in Ref. [37]), and the friction coefficient between a sphere and a hopper surfaces is 0.40±0.09. These hoppers can be attached to a cylindrical container with an inner diameter of 19 cm≈12.0dand a height of 18 cm≈11.4d. During each experimental trial,the hopper and container were first attached to an additional large tube from above. Then,more than 1800 spheres were rapidly poured into the tube,and filled the container and hopper with a closed outlet to a height of about 14d. Next, the outlet was left open to discharge the spheres until a spontaneous clogging occurs. Typically, the number of spheres that flew out of the hopper before clogging was several hundred,and about 1200 spheres on average remained in the hopper and container,corresponding to a packing of about 10din height. During the pouring process,the hopper and container were steadily placed on the ground with no external vibration exerted. Data of two kinds of cases were discarded as the flow field and/or packing states before clogging might not be stationary and representative: (i) only few particles just at the orifice drop out while the remaining particles has already formed a clogged structure with no movements at all;(ii)height of the remaining packing above the orifice after clogging is lower than about 5d. The temperature and relative humidity of the lab was controlled to be around 20°C and 50%.

    Fig.1. (a)Slices on two orthogonal planes of a raw magnetic resonance image,with part of the reconstructed packing structures superimposed.(b)A schematic diagram of the hopper geometry and the method to identify clogging structures: the shaded area near the orifice represents a free volume in which a sphere can move continuously without particle overlapping,and the yellow spheres outlining the upper boundary of this region compose the clogging structure. [(c)and(d)]Two reconstructed packing structures with(c)Dout=3.5d,αslope=0° and(d)Dout=3d,αslope=60°,in which the constituent spheres of clogging structures are highlighted with a golden color. [(e)and(f)]The same reconstructed structures as those in panels(c)and(d)viewed upward through orifices. [(g)and(h)]Average clogging size〈Nclog〉as a function of(g)orifice diameter Dout and (h) hopper slope angle αslope. The error-bar on each data point represents a 75% confidence interval, and are basically proportional to the standard variance of Nclog. The solid lines represent results of theoretical model according to Eqs.(3)and(4). A schematic diagram of the method to estimate sizes of clogging structures in hoppers with non-zero slope angles is shown in the inset of panel(h).

    The clogged packing structures were then scanned with a 3T MRI system (MAGNETOM Prisma Fit; Siemens Healthcare,Erlangen,Germany),using a spoiled gradient echo scanning sequence (Turbo-FLASH),[38]and a parallel acquisition technique(GRAPPA),[39]and were finally reconstructed using image processing methods(Figs.1(c)-1(f)). The spatial resolution of the MRI is 1 mm/pixel,and the average uncertainty of sphere centroids is estimated to be around 6%of the diameter.It is observed that spheres near container boundary form layered structures,which are discarded in the following analyses.For each hopper,the aforementioned trials are repeated for at least six times for better statistics: a total of 75 sets of packing structures were analyzed for eleven different outlets. The average packing fraction of packing structures away from the orifice is 0.599±0.005. Contacts between spheres are identified using an error-function fitting method,[40]and the average contact number is calculated to be 5.30±0.05.

    In reconstructed packing structures, a layer of spheres that block the orifice can be clearly observed as shown in Figs. 1(c)-1(f). But compared with the one-dimensional arches in 2D disk systems which can be identified instinctively, a rigorous definition for this 3D clogging structures is prerequisite to analyze their properties. To this end,we come up with the following operational definition (as illustrated in Fig.1(b)). First,a free volume near the hopper orifice is identified which comprises all points that are at least 1daway from the centroids of any spheres. By definition, this free volume can accommodate another sphere with a diameter ofdwithout particle overlapping. Then, all spheres whose centroids are 1daway from the boundary of this free volume is defined to compose the clogging structure. Apparently, if the positions of all particles in the clogging structure are fixed, there’s no way to move any other spheres in the packing continuously out through the orifice without penetrating these blocking spheres.This method rules out some extremely concave clogging structures which are very likely to be mechanically unstable,unless the interparticle friction coefficient is unrealistically large.For example, as shown in the schematic diagram (Fig. 1(b)), assuming that particles 3 squeezes in between particles 1 and 2 by slightly separating them,it would not be identified as a part of the clogging structure. Therefore, this definition is practically compatible with the arches identified in granular disk 2D systems. Clogging structures of 2D monodispersed disks would be identical to the arches near hopper outlet if identified following the procedures introduced in this work, as long as all angles between connections of adjacent contacts are larger than 120°. The probability of having contact angles smaller than 120°is indeed negligible as shown in Ref.[25].

    3. Experimental results

    We first study the dependences of the structure sizeNclog,i.e.,number of particles within each clogging structure,on orifice diameter and hopper angle. According to a simple geometric consideration, the projection area of clogging structures over the orifice should scale asD2out, and therefore, a plausible quadratic relationship between〈Nclog〉and orifice diameterDout,i.e.,〈Nclog〉∝D2outcan be inferred. However,interestingly, a linear relation has been observed in our system(Fig.1(g)),which are actually more similar to that in the 2D systems.[14,24,41]We propose that this counter-intuitive linear relation stems from an intrinsic correlation within clogging structures due to the mechanical stability requirement,and its length scale is comparable with the orifice diameter. For a strongly correlated clogging structure,a quadratic model associated with covering the orifice with randomly-placed spheres would be unsuitable. For hoppers with non-zero slope angles,a conventional decreasing〈Nclog〉with increasingαslopeis observed (Fig. 1(h)).[14]From a mechanical point of view, this stems from the finite friction coefficient between hopper inner surface and particles: the maximum vertical component of a fully activated frictional force depends onαslope, which sets an upper bound of the size of clogging structure that can be supported by the hopper’s inner wall. In Section 4, we introduce a phenomenological model incorporating both the structural correlation and hopper boundary condition to account for the dependences of〈Nclog〉onDoutandαslope. Furthermore,we notice that the standard variance ofNclogdecreases withαslope, while no clear dependence onDoutcan be observed.We propose that this is because a larger hopper slope imposes a stronger constraint on the shape and size of clogging structures against their intrinsic fluctuations with the mechanical boundary condition.

    To provide more information beyond the average sizes of clogging structures,[31]we further quantify the overall shape of the clogging structures using its inertia tensorM,the component of which is

    Furthermore,to quantify the clogging structure shape,we define an aspect ratio as

    Fig. 2. (a) Radius of gyration and (b) aspect ratio of clogging structures as a function of structure size. Each point represents a clogging structure in hoppers with different orifice diameters and/or slop angles. The solid line in panel (a) represent a scaling relationship: Rg =N1/2clog. The two horizontal dashed lines in panel (b) represent the aspect ratios of a homogeneous hemispherical shell(0.79)and a circular disk(0.71).

    To further study the load-bearing properties which is expected to be encoded in the geometry of a clogging structure,we analyze the relative positions of contacting spheres in the packing by calculating: (i) a contact orientation angleθzbetween the gravity direction and a vector connecting the centroids of each pair of contacting spheres,and(ii)a three-point contact angleθcspanned by one sphere and two of its contacting neighbors(see insets of Fig.3). These angles are crucial for the establishment of statistical models predicting the clogging probability and for the force analyses of clogging structure stability,[23,25]and their variance is a key indicator of arch-breaking in two dimensions.[28,29]In Fig. 3, we demonstrate the distributions ofθzandθcfor two groups of particles:(a: the clogging structure group)those in clogging structures,and only contacting neighbors which are also in the clogging structure are included when calculatingθc, and (b: the background group)all particles in the packing apart from those in contact with the container and hopper. Distributions of clogging structures obtained in hopper with different orifices are averaged, as we have verified that the variations of the distributions whenDoutorαslopechanges are smaller than the fluctuations of different clogging structures over a fixed orifice. For particles in the background group,θzexhibits a uniform distribution, and in comparison, particles in clogging structures have a clear preference for larger values ofθz, that is, there are more contacts lying close to the horizontal plane. Similarly, particles in the clogging structure favorθcwith large values close to 180°, corresponding to a quasi-linear chainlike contact geometry,which is a clear feature of load-bearing structures in granular materials.[43]In addition,the peak at 60°of the distributionp(θc) of the clogging structure group is lower than that of the background group, which corresponds to a locally dense structure of an equilateral triangle formed by three contacting spheres. This further implies that contacting spheres in the clogging structure tend to form branch-like networks. Although similar analyses have been conducted for 2D arches,[41]only with the comparison with particles away from the orifice as proposed here,the uniqueness of clogging structures associated with their special mechanical status is revealed. This structural anisotropy also reminds us of the correspondence between strong force chains and arches identified in 2D systems.[44]

    The aforementioned geometric features along with the implied requirements of load-bearing and mechanical stability suggest a possible close relationship between the clogging structure and mutually-stabilized bridge structures in jammed granular packings identified previously.[45]A bridge structure is a configuration of several contacting spheres satisfying the following three-point mutual-stability condition: one sphere is stabilized by a group of three contact neighbors if a vector from its centroid along the gravity direction goes through the triangle formed by the centroids of the three particles(i.e.,the triplet). Note that the triplet with the lowest centroid of mass is chosen if there exist more than one triplet satisfying above criteria.[46]Two spheres are mutually-stabilized if one of them rests on a triplet including the other one, andvice versa(see inset of Fig. 4(b)). We find that, for all pairs of contacts between spheres in the clogging structure, about 7.8% of them form mutually-stabilized configurations, and in comparison,for all contact pairs in the packing,this fraction is reduced to 5.2%. This difference verifies an enhanced cooperativity associated with the clogging structures. In addition, bridge structures have been identified in our packing,and we analyze the bridge size: number of spheres composing each bridge. Interestingly,the average size of bridges near the orifice is quite close to that of all bridges in the packing,and is basically independent ofDoutandαslope. Therefore,the average bridge size is determined by the global packing fraction alone,[45]and represents an intrinsic length scale associated with the cooperative stability.

    Fig. 4. (a) Genus of clogging structures versus structure size. The solid line represents G=0.19. (b)The probability distribution of contact number between constituent spheres of clogging structures,which have also been verified to be basically independent of hopper dimension. (Inset)A schematic diagram of a typical bridge structure:particle 1(2)is supported by particles 2,3,and 4(1,5,and 6),and therefore particles 1 and 2 are mutually stable.

    4. Theoretical model

    Last but not least, we present a simple model to account for the dependence of clogging structure size on hopper dimension based on above observations.As defined in Section 3,each bridge structure identified in the packing is a collection of mutually-stabilized spheres,while correlation between bridge structures should be relatively weak,since no pairs of spheres from two different bridges support each other. As a crude model, we propose that a clogging structure is composed by a few independent structural motifs of several spheres,such as the bridge structures defined above. On the other hand,a necessary condition for a blocked orifice is that there are at least about (2πDout/d)/2 spheres arranged along its periphery to prevent other spheres from flowing through it. Likewise,these spheres should be only very weakly-correlated since there is also no mutual-stability requirement between them. Therefore,we propose that the average number of independent motifs in one clogging structure is proportional to, if not very close to,the number of spheres at the periphery,so that a simple relationship between the average clogging structure size and orifice diameter:

    is obtained.The proportionality factorncorresponds to the average number of spheres associated with one structural motif in clogging structures,and should be quite close to the average number of spheres in one bridge structure. By fitting the experimental data(Fig.1(g)),we calculatento be 2.6±0.1. We also evaluate the average bridge size, which is the number of spheres composing all bridges in each packing divided by the number of bridges,and obtain a self-consistent average value of 2.4±0.1.

    The foundation of above proportionality hypothesis is the structural correlation of spheres in clogging structures.This linear relationship therefore reinforces the conception that only strongly-correlated structures with load-bearing ability are selected during the flow process to eventually halt the flow.[36]Another direct corollary of this hypothesis is that there should be a typical largest orifice, beyond which clogging could hardly happen,since the area covered by spheres in the clogging structures increases much more slowly withDoutthan the orifice area. For a fixedαslope,this critical maximum diameter should be simply limited by this intrinsic correlation length,that is,the average bridge size.

    This expression is consistent with the experimental results as shown in Fig. 1(h). A crude examination of results on the length of 2D arches suggests that above geometric model also works satisfactorily when modified to apply to the 2D case.[14]

    5. Concluding remarks

    Using MRI techniques and a series of structural quantification methods,the 3D clogging structures of granular spheres near hopper orifice have been systematically quantified. Clear geometric features associated with the load-bearing abilities of clogging structures have been revealed, and dependences of the size and shape of clogging structures on hopper dimensions have been studied. Despite that the overall characteristics is qualitatively similar to that of arches in 2D disk systems, the 3D clogging structure is found to be close to a branching network,so that it involves a complex combination of some elementary linear load-bearing motifs. The requirement of load-bearing for clogging structures inevitably introduces a correlation length,which is arguably the typical size of cooperatively-stable bridge structures. This is consistent with the quasi-linear strong force chains observed with 2D photoelastic granular disks,[47]which we suggest is the underlying mechanical reason why only configurations with specific geometric features identified here are selected out to form the clogging domes. Meanwhile,this correlation also implies that the formation of a clogging structure is not a random collection of independent spheres,but is associated with the dynamical cooperativity of the heterogeneous dense granular flow.Hypothetically,it is possible to identify precursory fluctuating structures in the flow prior to clogging with similar geometric features as the clogging structures, which should explain the stochastic nature of clogging.

    Phenomenologically,a direct consequence of the correlation length of load-bearing structures being comparable to the orifice is the linear dependence of the clogging structure size on the orifice diameter,which is unexpected and non-trivial for 3D systems,and should be tested in future with more accuracy.We also propose that based on our observations,more sophisticated stochastic models may be established,such as a randomwalk percolation model on a hemispherical surface with an additional short-range correlation,to provide a solid structural mechanism for the clogging phenomenon of 3D granular

    materials.

    Acknowledgments

    Project supported by the China Postdoctoral Science Foundation (Grant Nos. 2018M641957 and 2019T120319),the National Natural Science Foundation of China (Grant No.11904102),and the Fundamental Research Funds for the Central Universities,China.

    猜你喜歡
    宮殿
    南漢國一號(hào)宮殿發(fā)掘與初步復(fù)原
    廣州文博(2023年0期)2023-12-21 07:22:52
    神秘宮殿
    去地下宮殿做客
    記憶宮殿
    從古希臘的宮殿走來
    優(yōu)雅(2017年12期)2017-12-08 23:32:12
    阿掘的宮殿
    南漢國宮殿十六獅柱礎(chǔ)石
    大眾考古(2015年5期)2015-06-26 08:21:46
    故宮的宮殿少半間
    宮殿,不是只能建在陸地上。
    海洋世界(2014年4期)2014-12-23 12:58:07
    滾來滾去的宮殿
    日韩有码中文字幕| 99国产精品一区二区蜜桃av| 噜噜噜噜噜久久久久久91| 国产精品自产拍在线观看55亚洲| a级毛片a级免费在线| 51午夜福利影视在线观看| 十八禁人妻一区二区| 国产亚洲欧美98| 国产高清videossex| 国产精品乱码一区二三区的特点| 成人av一区二区三区在线看| 色吧在线观看| 岛国在线免费视频观看| 美女 人体艺术 gogo| 天堂网av新在线| 欧美激情久久久久久爽电影| 亚洲 国产 在线| 一区福利在线观看| 欧美日韩乱码在线| 免费观看人在逋| 精品无人区乱码1区二区| 国产在视频线在精品| 亚洲国产日韩欧美精品在线观看 | 国产不卡一卡二| 亚洲国产精品999在线| 一本综合久久免费| 精品人妻一区二区三区麻豆 | 看免费av毛片| svipshipincom国产片| 午夜激情福利司机影院| 一级黄片播放器| 法律面前人人平等表现在哪些方面| 可以在线观看毛片的网站| 在线十欧美十亚洲十日本专区| 国产亚洲精品一区二区www| 国产成+人综合+亚洲专区| 国产毛片a区久久久久| 亚洲精品一卡2卡三卡4卡5卡| 欧美区成人在线视频| 色哟哟哟哟哟哟| 日韩人妻高清精品专区| 国产精品电影一区二区三区| 国产精品日韩av在线免费观看| 国产主播在线观看一区二区| 国产精品久久久久久亚洲av鲁大| 老熟妇仑乱视频hdxx| 国产探花极品一区二区| 99精品在免费线老司机午夜| 欧美乱码精品一区二区三区| 午夜视频国产福利| 精品久久久久久久久久免费视频| 久久国产精品人妻蜜桃| 一个人观看的视频www高清免费观看| 国产精品日韩av在线免费观看| АⅤ资源中文在线天堂| 亚洲国产欧美人成| 有码 亚洲区| 国产av不卡久久| 久久久久久久午夜电影| 欧美一区二区国产精品久久精品| 亚洲精品粉嫩美女一区| 天天添夜夜摸| 一夜夜www| а√天堂www在线а√下载| 好男人在线观看高清免费视频| 国产真实乱freesex| 熟女人妻精品中文字幕| 欧美乱色亚洲激情| 亚洲 国产 在线| 老汉色av国产亚洲站长工具| 亚洲欧美日韩高清专用| 中文字幕av在线有码专区| 午夜两性在线视频| 久久精品国产亚洲av香蕉五月| 偷拍熟女少妇极品色| 亚洲成av人片免费观看| 日韩大尺度精品在线看网址| 亚洲精品在线美女| 深夜精品福利| 免费观看的影片在线观看| 欧美极品一区二区三区四区| 少妇熟女aⅴ在线视频| 97碰自拍视频| 久久久久久久午夜电影| 国产淫片久久久久久久久 | 18美女黄网站色大片免费观看| 色视频www国产| 99精品在免费线老司机午夜| 国产高潮美女av| 久久久久久久久大av| 亚洲七黄色美女视频| aaaaa片日本免费| 99精品在免费线老司机午夜| 一进一出抽搐gif免费好疼| 国产精品一及| 十八禁人妻一区二区| 亚洲第一电影网av| 99视频精品全部免费 在线| 在线观看一区二区三区| а√天堂www在线а√下载| 国产精品久久久久久久久免 | 色综合婷婷激情| 色哟哟哟哟哟哟| 亚洲美女黄片视频| 国内毛片毛片毛片毛片毛片| 99久国产av精品| 国产精品久久电影中文字幕| 又黄又爽又免费观看的视频| 一区福利在线观看| 日本 欧美在线| 亚洲欧美精品综合久久99| 精品久久久久久久人妻蜜臀av| 又黄又爽又免费观看的视频| 国产精品香港三级国产av潘金莲| 深爱激情五月婷婷| 性欧美人与动物交配| 午夜视频国产福利| 免费看a级黄色片| 欧美日韩中文字幕国产精品一区二区三区| bbb黄色大片| 一级毛片女人18水好多| 搞女人的毛片| 在线观看一区二区三区| 中文字幕熟女人妻在线| 欧美色欧美亚洲另类二区| 男女那种视频在线观看| 免费在线观看日本一区| 欧美一区二区国产精品久久精品| 免费大片18禁| 天天躁日日操中文字幕| 免费看美女性在线毛片视频| 岛国在线免费视频观看| 99riav亚洲国产免费| 国产成人欧美在线观看| 婷婷亚洲欧美| 一区二区三区国产精品乱码| 757午夜福利合集在线观看| 久久香蕉精品热| 精品一区二区三区视频在线观看免费| 日韩欧美 国产精品| 久久香蕉精品热| 国产高潮美女av| 亚洲成人精品中文字幕电影| 久久久成人免费电影| 丁香欧美五月| 中文字幕人妻熟人妻熟丝袜美 | 啦啦啦韩国在线观看视频| 真人做人爱边吃奶动态| 亚洲中文日韩欧美视频| 国产精品99久久久久久久久| 老司机在亚洲福利影院| 人妻夜夜爽99麻豆av| 午夜影院日韩av| 久久精品国产亚洲av香蕉五月| 久久香蕉国产精品| 男女之事视频高清在线观看| 亚洲成av人片免费观看| 久久久国产成人精品二区| 丰满人妻熟妇乱又伦精品不卡| 黄色日韩在线| 少妇的逼好多水| 欧美一区二区国产精品久久精品| 日本撒尿小便嘘嘘汇集6| 9191精品国产免费久久| 国产一区在线观看成人免费| 亚洲av熟女| 亚洲av成人精品一区久久| 午夜福利在线观看免费完整高清在 | 亚洲午夜理论影院| 亚洲精品成人久久久久久| 免费看美女性在线毛片视频| 91麻豆av在线| 在线观看免费视频日本深夜| 国产成人福利小说| 日本黄色片子视频| 国产毛片a区久久久久| av欧美777| 丝袜美腿在线中文| 午夜久久久久精精品| 亚洲中文字幕日韩| 久久久久久大精品| 精品午夜福利视频在线观看一区| 国产成+人综合+亚洲专区| 人人妻人人澡欧美一区二区| 69av精品久久久久久| 色尼玛亚洲综合影院| 婷婷精品国产亚洲av| 他把我摸到了高潮在线观看| 免费人成视频x8x8入口观看| 观看免费一级毛片| 成人国产一区最新在线观看| 国内毛片毛片毛片毛片毛片| 中出人妻视频一区二区| 亚洲色图av天堂| 男女做爰动态图高潮gif福利片| 丝袜美腿在线中文| xxxwww97欧美| 亚洲在线观看片| 国产精品99久久久久久久久| 99久久九九国产精品国产免费| 精品国产亚洲在线| 亚洲真实伦在线观看| 午夜免费观看网址| 欧美性猛交╳xxx乱大交人| 亚洲七黄色美女视频| av女优亚洲男人天堂| 亚洲男人的天堂狠狠| 午夜免费男女啪啪视频观看 | 成人18禁在线播放| 一级毛片高清免费大全| 成人精品一区二区免费| 婷婷精品国产亚洲av| 人人妻,人人澡人人爽秒播| 国产av在哪里看| 97超视频在线观看视频| 99精品久久久久人妻精品| 亚洲午夜理论影院| 搡老岳熟女国产| 国产精品久久久久久人妻精品电影| 国产精品嫩草影院av在线观看 | 国产v大片淫在线免费观看| 黑人欧美特级aaaaaa片| 国产探花极品一区二区| 我的老师免费观看完整版| 国产精品av视频在线免费观看| 国产精品自产拍在线观看55亚洲| 一区二区三区免费毛片| 性欧美人与动物交配| 免费在线观看日本一区| 搞女人的毛片| 黄色成人免费大全| 麻豆成人av在线观看| 在线十欧美十亚洲十日本专区| av在线蜜桃| 丁香欧美五月| 欧美av亚洲av综合av国产av| 国产蜜桃级精品一区二区三区| 在线观看一区二区三区| 欧美黑人欧美精品刺激| 91久久精品国产一区二区成人 | 非洲黑人性xxxx精品又粗又长| 午夜久久久久精精品| 国产探花极品一区二区| 国产欧美日韩一区二区三| 国产麻豆成人av免费视频| 欧美+日韩+精品| 久久精品国产自在天天线| 欧美xxxx黑人xx丫x性爽| 最新美女视频免费是黄的| 两个人的视频大全免费| 乱人视频在线观看| 精品人妻一区二区三区麻豆 | 一区二区三区高清视频在线| 亚洲av电影在线进入| 国产中年淑女户外野战色| 熟妇人妻久久中文字幕3abv| 久久中文看片网| 国产精品美女特级片免费视频播放器| 中文字幕av在线有码专区| 亚洲成人中文字幕在线播放| 亚洲欧美日韩高清专用| 91在线观看av| 色噜噜av男人的天堂激情| 国产亚洲av嫩草精品影院| 亚洲中文日韩欧美视频| 久久亚洲精品不卡| 亚洲欧美日韩高清专用| 91九色精品人成在线观看| 欧美在线黄色| 一本一本综合久久| 在线观看舔阴道视频| 国产一区二区在线观看日韩 | 熟女少妇亚洲综合色aaa.| 在线a可以看的网站| 久久精品影院6| 夜夜夜夜夜久久久久| 国产精品1区2区在线观看.| 国产国拍精品亚洲av在线观看 | 中文字幕av在线有码专区| 一级a爱片免费观看的视频| 午夜亚洲福利在线播放| 最近最新中文字幕大全电影3| 成人av一区二区三区在线看| 丝袜美腿在线中文| 在线国产一区二区在线| 宅男免费午夜| 两人在一起打扑克的视频| 亚洲精品一卡2卡三卡4卡5卡| 一二三四社区在线视频社区8| 欧美成狂野欧美在线观看| 成人永久免费在线观看视频| 亚洲激情在线av| 国产成人av教育| 手机成人av网站| 亚洲18禁久久av| 欧美3d第一页| 在线天堂最新版资源| 亚洲狠狠婷婷综合久久图片| 国产 一区 欧美 日韩| 国内少妇人妻偷人精品xxx网站| 欧美一区二区亚洲| 精品久久久久久久久久免费视频| 丁香欧美五月| 黄色片一级片一级黄色片| 少妇熟女aⅴ在线视频| 99精品欧美一区二区三区四区| 国产日本99.免费观看| av天堂中文字幕网| 久久久国产精品麻豆| 国内精品久久久久久久电影| 熟女电影av网| 精品久久久久久久毛片微露脸| 蜜桃久久精品国产亚洲av| 日韩欧美精品v在线| 在线免费观看不下载黄p国产 | www日本在线高清视频| 极品教师在线免费播放| 在线播放无遮挡| АⅤ资源中文在线天堂| 啦啦啦观看免费观看视频高清| 国产蜜桃级精品一区二区三区| 色综合欧美亚洲国产小说| 人人妻,人人澡人人爽秒播| 成人性生交大片免费视频hd| 熟妇人妻久久中文字幕3abv| 亚洲国产欧洲综合997久久,| 色av中文字幕| 女人高潮潮喷娇喘18禁视频| 精品一区二区三区视频在线 | 十八禁人妻一区二区| 香蕉久久夜色| 日韩欧美三级三区| 老司机午夜福利在线观看视频| 麻豆久久精品国产亚洲av| 欧美激情在线99| 99热精品在线国产| 亚洲美女视频黄频| 精品无人区乱码1区二区| 久久香蕉精品热| 一个人免费在线观看电影| 观看免费一级毛片| 少妇丰满av| 18禁黄网站禁片午夜丰满| av福利片在线观看| 99国产综合亚洲精品| 黑人欧美特级aaaaaa片| 国内揄拍国产精品人妻在线| 非洲黑人性xxxx精品又粗又长| 亚洲精品国产精品久久久不卡| 亚洲成人久久爱视频| 老司机深夜福利视频在线观看| 亚洲精品色激情综合| 成年女人毛片免费观看观看9| 精品一区二区三区人妻视频| а√天堂www在线а√下载| 可以在线观看毛片的网站| 草草在线视频免费看| tocl精华| 变态另类成人亚洲欧美熟女| 国产黄色小视频在线观看| 久久久精品大字幕| 美女免费视频网站| 非洲黑人性xxxx精品又粗又长| 国产精品亚洲一级av第二区| 国产一区二区在线av高清观看| xxxwww97欧美| 91久久精品电影网| 国产成人影院久久av| 久久精品亚洲精品国产色婷小说| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av成人精品一区久久| 色噜噜av男人的天堂激情| 在线免费观看不下载黄p国产 | 美女高潮的动态| 亚洲国产高清在线一区二区三| 国产高清有码在线观看视频| 天堂av国产一区二区熟女人妻| 午夜福利18| 久久午夜亚洲精品久久| 日本成人三级电影网站| 精品免费久久久久久久清纯| 国产亚洲精品av在线| 国产一区在线观看成人免费| 国产成人av教育| 欧美三级亚洲精品| 最近最新免费中文字幕在线| www日本在线高清视频| 久久久久精品国产欧美久久久| 亚洲自拍偷在线| 精品熟女少妇八av免费久了| 999久久久精品免费观看国产| 亚洲aⅴ乱码一区二区在线播放| 少妇的逼水好多| 亚洲成人中文字幕在线播放| 国产69精品久久久久777片| 亚洲一区二区三区色噜噜| 亚洲国产精品成人综合色| 亚洲成av人片在线播放无| 十八禁网站免费在线| 久久草成人影院| 国产三级中文精品| 麻豆国产av国片精品| 在线a可以看的网站| 午夜亚洲福利在线播放| 中文资源天堂在线| 19禁男女啪啪无遮挡网站| 中文字幕人妻熟人妻熟丝袜美 | 熟女人妻精品中文字幕| 亚洲电影在线观看av| 亚洲国产精品成人综合色| 国产伦在线观看视频一区| 哪里可以看免费的av片| 国产aⅴ精品一区二区三区波| 久久精品国产自在天天线| 欧美丝袜亚洲另类 | 国产视频内射| 91字幕亚洲| 精品欧美国产一区二区三| 一个人看的www免费观看视频| 美女黄网站色视频| 午夜久久久久精精品| 亚洲片人在线观看| 99久久精品国产亚洲精品| 两个人的视频大全免费| 午夜精品在线福利| 熟女电影av网| 亚洲精品影视一区二区三区av| 搞女人的毛片| 久久久精品欧美日韩精品| 脱女人内裤的视频| 哪里可以看免费的av片| 国产在线精品亚洲第一网站| 午夜免费观看网址| 欧美大码av| 国产高清激情床上av| www国产在线视频色| 午夜免费成人在线视频| 黄色片一级片一级黄色片| 人妻夜夜爽99麻豆av| 18禁裸乳无遮挡免费网站照片| 99久久九九国产精品国产免费| 超碰av人人做人人爽久久 | 日韩中文字幕欧美一区二区| 日韩av在线大香蕉| 99久久九九国产精品国产免费| 国产免费男女视频| av欧美777| 久久国产精品人妻蜜桃| 久久精品综合一区二区三区| 国产伦在线观看视频一区| 免费人成视频x8x8入口观看| 国产午夜福利久久久久久| 别揉我奶头~嗯~啊~动态视频| 又粗又爽又猛毛片免费看| 亚洲国产欧美网| 男女床上黄色一级片免费看| 99精品在免费线老司机午夜| 听说在线观看完整版免费高清| 欧美日韩瑟瑟在线播放| 禁无遮挡网站| 99国产极品粉嫩在线观看| 色综合欧美亚洲国产小说| 久久精品国产清高在天天线| 国产精品野战在线观看| 51午夜福利影视在线观看| 网址你懂的国产日韩在线| 国产亚洲精品久久久com| 嫩草影视91久久| 99久久综合精品五月天人人| 一进一出好大好爽视频| 一级毛片高清免费大全| 女人十人毛片免费观看3o分钟| 在线观看66精品国产| 亚洲欧美日韩高清专用| 欧美日本亚洲视频在线播放| 国产久久久一区二区三区| 日韩欧美三级三区| 亚洲国产日韩欧美精品在线观看 | 少妇高潮的动态图| 真人一进一出gif抽搐免费| 久久久久亚洲av毛片大全| 欧美一级毛片孕妇| 久久久久久久久久黄片| 91九色精品人成在线观看| 日韩av在线大香蕉| 每晚都被弄得嗷嗷叫到高潮| 狂野欧美激情性xxxx| 国产欧美日韩一区二区三| 99久久精品一区二区三区| 国内毛片毛片毛片毛片毛片| 在线观看免费视频日本深夜| 非洲黑人性xxxx精品又粗又长| 欧美bdsm另类| 99热这里只有是精品50| 国产精品女同一区二区软件 | 在线观看舔阴道视频| 亚洲国产精品成人综合色| 精品一区二区三区视频在线观看免费| 少妇人妻精品综合一区二区 | 亚洲av免费在线观看| 少妇的逼水好多| 亚洲欧美日韩卡通动漫| 日本黄大片高清| 亚洲欧美激情综合另类| 久久这里只有精品中国| 中文字幕久久专区| 性色av乱码一区二区三区2| 亚洲精品一卡2卡三卡4卡5卡| tocl精华| 国产成人系列免费观看| xxx96com| av在线天堂中文字幕| 免费av不卡在线播放| 亚洲欧美日韩高清专用| 国产真实伦视频高清在线观看 | 欧美中文日本在线观看视频| 国产91精品成人一区二区三区| 成人一区二区视频在线观看| 男女之事视频高清在线观看| 日本黄色视频三级网站网址| xxxwww97欧美| 18禁国产床啪视频网站| 亚洲国产欧美网| 欧美xxxx黑人xx丫x性爽| 亚洲精品亚洲一区二区| 国产精品久久电影中文字幕| 国产黄a三级三级三级人| 欧美精品啪啪一区二区三区| 久9热在线精品视频| 亚洲成人中文字幕在线播放| 日韩大尺度精品在线看网址| 91在线精品国自产拍蜜月 | 性欧美人与动物交配| 18禁国产床啪视频网站| 国产午夜精品论理片| 黄片大片在线免费观看| 亚洲成人精品中文字幕电影| 亚洲在线自拍视频| 精品免费久久久久久久清纯| 免费在线观看影片大全网站| 亚洲av不卡在线观看| 国产黄a三级三级三级人| 国产精品一区二区三区四区免费观看 | 久久久久免费精品人妻一区二区| 亚洲精品久久国产高清桃花| av中文乱码字幕在线| 中文字幕人妻熟人妻熟丝袜美 | 久久国产精品影院| 老司机午夜十八禁免费视频| 色综合亚洲欧美另类图片| 亚洲午夜理论影院| 国产成人影院久久av| 日本与韩国留学比较| 男女那种视频在线观看| 一夜夜www| 12—13女人毛片做爰片一| 校园春色视频在线观看| 国产黄a三级三级三级人| 白带黄色成豆腐渣| 美女 人体艺术 gogo| 欧美午夜高清在线| 国产黄片美女视频| 狂野欧美白嫩少妇大欣赏| 亚洲片人在线观看| 午夜日韩欧美国产| 男插女下体视频免费在线播放| 日韩欧美精品v在线| 国产色婷婷99| 两个人视频免费观看高清| ponron亚洲| 高清在线国产一区| 超碰av人人做人人爽久久 | 国模一区二区三区四区视频| 欧美又色又爽又黄视频| 中文字幕人妻熟人妻熟丝袜美 | 波野结衣二区三区在线 | 国产午夜福利久久久久久| 欧美一级毛片孕妇| 成年女人永久免费观看视频| 国产精品1区2区在线观看.| 熟女人妻精品中文字幕| 久久这里只有精品中国| 在线观看舔阴道视频| 成年人黄色毛片网站| eeuss影院久久| 少妇的逼水好多| 黄片大片在线免费观看| 亚洲精品成人久久久久久| 99久久精品一区二区三区| 久久久久国内视频| 国产三级中文精品| 成人18禁在线播放| 日韩欧美免费精品| 99久久久亚洲精品蜜臀av| 精品午夜福利视频在线观看一区| 国产成人福利小说| 在线观看66精品国产| 欧美3d第一页| 国产一区二区在线av高清观看| 在线观看66精品国产| 欧美中文日本在线观看视频| 有码 亚洲区| 女同久久另类99精品国产91| 亚洲精品456在线播放app | 久久久久久久久久黄片| 欧美激情久久久久久爽电影| 一级黄片播放器| www.999成人在线观看| 国产精品98久久久久久宅男小说| e午夜精品久久久久久久| 免费高清视频大片| 一级黄色大片毛片| 国产v大片淫在线免费观看| 九九热线精品视视频播放| 国产午夜福利久久久久久|