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

    A STUDY OF DRAG COEFFICIENT RELATED WITH VEGETATION BASED ON THE FLUME EXPERIMENT*

    2010-07-02 01:37:53HUIErqingHUXing

    HUI Er-qing, HU Xing-e

    China Three Gorges Corporation, Yichang 443002, China

    State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China, E-mail: hui_erqing@ctgpc.com.cn

    JIANG Chun-bo, MA Fang-kai, ZHU Zhen-duo

    State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China

    A STUDY OF DRAG COEFFICIENT RELATED WITH VEGETATION BASED ON THE FLUME EXPERIMENT*

    HUI Er-qing, HU Xing-e

    China Three Gorges Corporation, Yichang 443002, China

    State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China, E-mail: hui_erqing@ctgpc.com.cn

    JIANG Chun-bo, MA Fang-kai, ZHU Zhen-duo

    State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China

    This study is focused on the effects of ecological factors (diameter and flexibility) and vegetation community composition on the drag coefficient related with vegetation. The single leafy shrub and three mixed communities (including shrub-grass, shrub-reed and reed-grass community) were studied. The flow velocity and water level were measured and used to calculate the drag coefficient based on the Bernoulli’s equation, Darcy drag formula and the expression for the drag coefficient related with Darcy drag factor. The trend of the drag coefficient in the vertical direction was analyzed against flow depth, diameter, diameter Reynolds number, flow depth Reynolds number and relative roughness height in different discharges. The results show that beside the dense leafy shrubs community, the vertical trend of the drag coefficient among other cases against flow depth, diameter, diameter Reynolds number, flow depth Reynolds number and relative roughness height can be approximately expressed by power law functions under different flow discharges. Moreover, in a mixed community with two plants with distinctly different ecological factors, the one with the most distinct variations of ecological factors determines the vertical trend of the drag coefficient; the other one only affects the magnitude of the drag coefficient. Furthermore, if the ecological factors of the vegetation in the vertical direction are kept almost not changed, the drag coefficient can be approximately regarded as a constant.

    vegetated channel, drag coefficient, flume experiment

    1. Introduction

    The flow dynamics conditions are important factors of the riverine wetland, which can affect theecological functions of the wetland. The hydraulic conditions in the wetlands are primarily determined by the flow drag related with vegetation. The factors which can affect the flow drag include the flow depth, profile shape, vegetation composition, sinuousness, bed form, bed load and covered ice, etc.[1]. The riverine wetland is composed of anisomerous grasses, shrubs and arbors. The vegetation includes submerged and unsubmerged vegetation, which considerably affects the hydrodynamics, sediment deposition, riverbed evolvement, nutrients and pollutant transportation[2,3]. All these are related with the blocked and altered course of flow, which is in turn related with the leaves and branches of the vegetation.The factors which affect the drag coefficient include the flexibility, diameter, height, distribution, density, age, health, distribution of branches and leaves, composition of plant and flow depth or velocity, etc.[4]. For large arbors, they not only bring about the drag, but also make the block[5]. The drag coefficient is the key in determining the drag related with vegetation, and also the key to understand the vertical distribution of the velocity. Many researches about the drag related with unsubmerged and submerged vegetation were carried out[6-9].

    In most of previous studies, the artificial plants are considered, which can hardly represent the distribution of branches and leaves of natural plants, and their results can not directly be used for the natural plants. The three dimensional flow fields among the plants concern the understanding of the vertical variations of the drag related with vegetation, in particular and the vertical variations of the drag coefficient, in general. The articles studies the vertical variation trend of the drag coefficient, based on the vertical variations of plant ecological factors (diameter and flexibility), leaves and plant types in the community. The different types of vegetation communities are considered and mounted in the flume in a stagger pattern, including the single leafy shrubs, reeds community and mixed shrubs-grass, shrubs-reeds community. The drag coefficient is calculated based on the Bernoulli’s equation, Darcy drag formula and the drag coefficient expression. The vertical trend of the drag coefficient is analyzed with respect to the water depth, diameter, diameter Reynolds number, water depth Reynolds number and relative roughness height in different discharges. The drag coefficients of a single community and the mixed community are compared. This study supplies some guidance in the determination of the drag coefficient in the vertical direction.

    Fig.1 Front elevation of shrub

    2. The experimental arrangement

    The experiments consider various effective factors that affect the drag coefficient such as plant leaves, composition and the vertical ecological factors (diameter and flexibility), shrubs and reeds. Larger diameter means less flexibility. Four communities are designed, of which two are single communities including leafy shrubs and reeds, other two are mixed communities including shrubs-grasses and shrubs-reeds. Buxus sincia Rehd (Fig.1) and the top segment of the natural reeds are taken to be simulated, which are mounted to a glass wall flume of 22.6 m long, 1.0 m wide and 0.8 m deep (Fig.2), with a slope of 0.0067. On average, the shrubs are 0.275 m in height and 0.2 m in width at the top in the biggest diameter in the longitudinal direction, while the roots are 0.07 m in height. The vertical distribution of the shrub’s diameter follows approximately a power law. The average height of reeds is 0.5 m and the diameter is 3.4 mm. Grasses are very flexible with height of 0.07 m. There are 12 pieces of bamboo leaves and pine needles, divided into three layers: 6 pieces in the bottom layer, 4 pieces in the middle layer and 2 pieces in the top layer. All the bamboo leaves and the pine needles amongst the simulated grasses are fixed in the flume. All plants are mounted in a staggered pattern. Two rows of grasses are placed between shrubs to simulate the mixed community (Fig.3). The mount pattern of the mixed shrub-reed community is the same as the shrub-grass community.

    Fig.2 Schematic diagram of laboratory setup

    Fig.3 Staggered patterns of grasses and shrubs or reeds

    The discharges are conducted by a frequency conversion pump, with frequency value fixed to maintain a constant discharge throughout the process. Upon entering the channel, the flow passes through a series of honeycomb grids to straighten the flow to make the flow uniform across the width of the channel. Four sections are set for measuring velocity and water level. Four vertical lines in each section are set, where are several measurement points in accordance with the flow depth (Fig.4). Water levels are measured by the pressure transducer, which are connected to the software running on a PC. The three components of flow velocity are measured with a three-dimensional Acoustic Doppler Velocimeter (ADV).

    Fig.4 Positions of the measurement points in the plane

    Table 1 Hydraulics and vegetation characteristic parameters in the single plant community

    According to the height of shrubs and reeds and preliminary measurements for the uniform flow in the non-vegetated channel, three types of discharges are determined: unsubmerged, just submerged and submerged, with discharge of 67.4 l/s, 108.6 l/s and 228.6 l/s, respectively. The Reynolds number of this experiment based on the flow depth is within the ranges of 2.20×104to 1.61×105. That based on the plant diameter is within the ranges of 1.0×103to 5.7×104. The Fraud number is less than 0.37, which indicates that the present experiment concerns a subcritical flow. The density of shrubs is 15.71 m–2in the dense conditions and 7.85 m–2in the sparse conditions. The density of the single reeds community is 59.88 m–2. The hydraulic and vegetation characteristic parameters in the single plant community are shown in Table 1 and those in the mixed community are shown in Table 2.

    Table 2 Hydraulic and vegetation characteristic parameters in the mixed plant community

    Fig.5 Schematic diagram of vegetated zone

    3. Theoretical analysis

    3.1The f relation between drag coefficient and Darcy friction factor

    It is assumed that there is a vegetated zone with Δxas the length in the longitudinal direction and Δyas the width in the lateral direction, which are denoted asΔxvand Δyvbetween the planes (Fig.5). The plant height ishv, the flow depth ishand the flow velocity in the cell isu. The forces on the vegetated zone are shown in Fig.6.

    Fig.6 Schematic diagram of the forces on the vegetated zone

    In Fig.6,Fv,Fb,Fsare the drag forces related with the plants, bed and water surface, respectively,Gis the gravity due to the plant, On this vegetated zone,Fvcan be written as:

    whereApis the momentum absorption area,ρis the water density,ρvis the vegetation density andρv=(ΔxvΔyv)?1,Cdis the drag coefficient related with vegetation.

    The volume of vegetation should be considered in a dense plant community, and the gravity forceFgdue to vegetation alongside the slope surface can be written as:

    whereVvegis the volume of vegetation.

    In the dense vegetated zone, the drag force related with vegetation is higher than the friction drag related with the bed roughness and water surface, and the latter can then be neglected. The gravity and the drag force related with vegetation approximately reach a balance under uniform flow conditions. Based on the force balance, it can be deduced thatFv=Fg. According to formulas (1) and (2), the following equation can be attained:

    Formula (4) can be transformed to:

    3.2The vertical distribution of the drag coefficient

    Under the gradually varied flow conditions, in every streamline, the energy lossHfcan be calculated by measuring the head loss based on the Bernoulli’s Eq.(6), and then Darcy friction factorfcan be calculated based on the Darcy-Weisbach Eq.(7). Based on these, the drag coefficientCdcan be calculated by formula (5). Both the potential and velocity heads are incorporated in the calculations. Bernoulli’s equation and Darcy friction factor formula are as follows wherevis the velocity of each measurement point,his the corresponding water level related with the velocity,σis the velocity distribution coefficient representing the flow properties.αapproaches 1.0 under a uniform velocity, which is taken as equal to 1.1 in this study. The subscripts 1 and 2 represent upstream and downstream, respectively.Lis the distance between upstream and downstream. The energy loss due to the flume bottom and the glass wall is much less than that due to the plant, as can be seen by comparing the vegetated and non-vegetated conditions, and it can be neglected. Because it is difficult to determine the streamline positions, whileCdof each measurement point is calculated, the flow velocity values of the upstream and downstream sections are taken as those at the same water level. It is assumed that the energy lossHfin each measurement point is the same as that on the water surface in every section. Because the ecological factors of the plant in the flume are similar, the above-mentioned assumptions are reasonable.

    The plant communities used in the experiment include shrubs, reeds, grasses and mixed communities composed of these plants. The relation between the leafy shrub’s diameter and its height is as follows:

    The relation between the leafless shrub plant’s diameter and its height is as follows:

    The diameter of the mixed plant community is calculated as:

    wheredmis the diameter of the mixed plant community.

    The momentum adsorption areaApof shrubs and grasses is approximately in a form of triangle, theApof reeds is a form of rectangle. TheApof the mixed community is the sum of the area of all plants in the community. The calculation methods of these plants are shown in Table 3. The volumes of shrubs, grasses and reeds are calculated by the formulas as listed in Table 4.

    Table 3 The formulas forAp

    Table 4 The formulas for plant volume

    In Tables 3 and 4,ds,dganddrare the diameter of the shrubs, grasses and reeds, respectively,hs,hgandhrare the height of the shrubs, grasses and reeds, respectively.

    Fig.7 TheCdagainstd,h,Reh,Redandk/hin leafy shrubs community. The symbols▲, ■ and◇represent the discharges of 67.4 l/s, 108.6 l/s and 228.6 l/s, respectively, Figs.7(a1)-7(e1) and Figs.7(a2)-7(e2) are for dense and sparse plants, respectively

    4. The experimental results

    The drag coefficientCdrelated with vegetation, the plant diameterd, the flow depthh, the flow depth Reynolds numberReh, the plant diameter Reynolds numberRedand the relative roughnessk/hin every water height are averaged by four profiles in each conditions. The dependence ofCdond, h,Reh,Redandk/hin different plant communities is analyzed as in the following sections. 4.1Leafy shrubs community

    The experimental results of leafy willows as reported in J?rvel?[10]show thatfincreases with the flow depth, indicating thatCdvalues increase with the flow depth. From Figs.7(a1)-7(e1), it can be seen thatCdvalues increase withd, h,Reh,Redandk/hin power law functions below the plant height. Above the plant height,Cdvalues decrease. Comparing with Juha’s results, it is shown that only if among dense plants and bellow the plant height,Cdvalues increase with the flow depth. Kadlec suggested that the Darcy friction factor of the flexible and emergent plant would decrease with the increase of flow depth[11]. Wu et al.[12]shows that the Manning roughness coefficient of the emergent plant decreases with flow depth increasing. Musleh and Cruise[13]shows that the Chezy coefficient increases with flow depth but nonlinearly. Formula (5) shows thatCdhas a positive correlation with the Darcy friction factor and the Manning coefficient and a negative correlation with the Chezy coefficient. From the above results, it can be deduced thatCddecreases with flow depth increasing. Figs.7(a2)-7(e2) show that the relations ofCdagainstd, h,Reh,Redandk/hcan well be described by declining power law curves, with correlation coefficient up to 0.77.

    4.2Shrub-grass mixed community

    Figure 8 shows thatCddecreases withd, h,Reh,Redandk/hin power law curves. The correlation of the curves is up to 0.92. ThoughCdsees a little difference in different discharges, they stick close to the same curve. Comparing theCdvalues in different plant densities, the scattering ofCdvalues in the sparse plant community is larger than that in the dense plant community. This shows that the effect of discharges (flow velocity) onCdin a sparse plant community is bigger than that in a dense plant community. The drag force related with vegetation is positively correlated with the deflectiondegrees of leaves and branches of the plant. The force of the flow on the vegetation increases with discharges increasing. The leaves and branches of the plant incline to follow the direction of the streamline, therefore, the drag force related with vegetation is reduced. In each discharge, the force of the flow imposed on the plant in the sparse plant community is larger than that in the dense one, which suggests that the deflection of the plant is more distinct and the variations ofCdwith discharges are more prominent in the sparse plant community.

    Fig.8 The plots ofCdagainstd,h,Reh,Redandk/hin leafy shrubs community. The symbols ▲, ■ and◇represent the discharges of 67.4 l/s, 108.6 l/s and 228.6 l/s, respectively, Figs.8(a1)-8(e1) and Figs.8(a2)-8(e2) are for dense and sparse plants, respectively

    4.3Discussions

    Comparing the single leafy shrubs community (Fig.7) with the shrubs-grasses mixed community (Fig.8), the drag coefficient related with vegetation of the latter is higher than that of the former. The trend of the curves of the drag coefficient againstd, h,Reh,Redandk/his different in dense conditions and is nearly the same in sparse conditions. The scattering of the drag coefficient in the shrubs-grasses mixed community is more distinct than that in the single shrubs community. In addition, well above the grasses height, the trend of the curves of the drag coefficient is the same. The effect of grasses on the drag coefficient is only limited within the ranges of a little over the height of grasses.

    The trend of curves of the drag coefficient against the plant diameter varies under different conditions. The drag coefficient of the single shrub community increases with the plant diameter, while the drag coefficient of the shrubs-grasses or shrubs-reeds mixed community decreases. It shows that the trend of the curves of the drag coefficient against the plant diameter is affected by the composition of the plants. The drag coefficient can be regarded as a constant when the vertical ecological factors of the plant (diameter and flexibility) keep little changed. The plant density has much greater effect on the vertical trend of the curves of the drag coefficient againstd, h,Reh,Redandk/h.

    5. The calculation model for the drag coefficient

    As the drag coefficient related with vegetation is a non-dimensional parameter, the average stem and leaf spacingsis divided by plant diameterdor flow depthhfor ensuring the non-dimensional nature of the drag coefficient. Based on the analysis of the experimental data and dimensional analysis of the Jonathan’s research[14], the drag coefficient can be expressed as a function ofd/s, h/s,RehandRed, as follows:

    The relation between every factor and the drag coefficient can be expressed as power or exponential functions. Formula (11) may take 16 different forms. The multiple correlation coefficients of these formsare calculated by the multiple regression analysis method, to show the precision of these forms. A comparison shows that the multiple correlation coefficient of the following form is the biggest:

    whereKis a coefficient,α1,α2,α3andα4are empirical parameters.

    Based on formula (12), the prediction model under each condition including the shrub in this experiment is expressed as follows:

    (1) For leafy shrub community For dense case:

    (2) For leafy shrub-grass mixed community For dense case:

    The statistics parameters for formulas (13)-(18) are listed in Table 5.

    Table 5 The statistics parameters

    In Table 5, the multiple correlation coefficient for formulas (13)-(18) is about 0.96 andFtest value is overF0.05=0.21, which shows that each formula has a high precision.

    6. Conclusions

    Based on the flume experimental data, the Bernoulli’s equation, the Darcy-Weisbach formula and the drag coefficient--the Darcy drag factor expression are used to calculate the drag coefficient in different discharges in various combinations of the flexible and stiff vegetation, such as natural shrubs, reeds and grasses. The dependence of the drag coefficient ond, h,Reh,Redandk/his

    (1) The trend of the curves of the drag coefficient against plant diameter, flow depth, diameter Reynolds number, flow depth Reynolds number and the relative roughness can be expressed as power functions, with different correlations.

    (2) While two types of plants coexist, the one with more distinct variations of the ecological factors analyzed. The following conclusions are reached:determines the trend of the curves of the drag coefficient in a certain range of height, while the other only affects the magnitude of the drag coefficient.

    (3) The drag coefficient is nearly a constant when the vertical ecological factors are kept little changed.

    The present study gives some food of thought in understanding the effects of discharges, the vertical variations of ecological factors and plant composition on the drag coefficient related with vegetation in different submerged conditions. The drag force can be calculated as soon as the numerical simulation results of flow field among the vegetation are available.

    [1] NIKORA V., LARNED S. and NIKORA N. et al. Hydraulic resistance due to aquatic vegetation in small streams: Field study[J].Journal of Hydraulic Engineering,2008, 134(9): 1326-1332.

    [2] NEPF H. M., VIVONI E. R. Flow structure in depthlimited, vegetated flow[J].Journal of Geophysical Research,2000, 105(C12): 28547-28557.

    [3] YEN B. C. Open channel flow resistance[J].Journal of Hydraulic Engineering,2002, 128(1): 20-39.

    [4] DARBY S. E. Effect of riparian vegetation on flow resistance and flood potentia[J].Journal of Hydraulic Engineering,1999, 125(5): 443-454.

    [5] WU Y., FALCONER R. A. and STRUVE J. Mathematical modeling of tidal currents in mangrove forests[J].Environmental Modelling and Software,2001, 16(1): 19-29.

    [6] WU Fu-Sheng. Characteristics of flow resistance in open channels with non-submerged rigid vegetation[J].Journal of Hydrodynamics,2008, 20(2): 239-245.

    [7] LIU Cheng, SHEN Yong-Ming. Flow structure and sediment transport with impacts of aquatic vegetation[J].Journal of Hydrodynamics,2008, 20(4): 461-468.

    [8] GU Feng-feng, NI Han-gen and QI Ding-man. Roughness coefficient for Unsubmerged and submerged reed[J].Journal of Hydrodynamics, Ser. B,2007, 19(4): 421-428.

    [9] ZHANG Jian-tao, SU Xiao-hui. Numerical model for flow motion with vegetation[J].Journal of Hydrodynamics,2008, 20(2): 172-178.

    [10] J?RVEL? J. Flow resistance of flexible and stiff vegetation: A flume study with natural plants[J].Journal of Hydrology,2002, 269(1-2): 44-54.

    [11] KADLEC R. H. Overland flow in wetlands: Vegetation resistance[J].Journal of Hydraulic Engineering,1990, 116(5): 691-706.

    [12] WU F. C., SHEN H. W. and CHOU Y. J. Variation of roughness coefficients for unsubmerged and submerged vegetation[J].Journal of Hydraulic Engineering,1999, 125(9): 934-942.

    [13] MUSLEH F. A., CRUISE J. F. Functional relationships of resistance in wide flood plains with rigid unsubmerged vegetation[J].Journal of Hydraulic Engineering,2006, 132(2): 163-171.

    [14] JONATHAN K. L., LISA C. R. and HARRY L. J. et al. Drag coefficients for modeling flow through emergent vegetation in the Florida Everglades[J].Ecological Engineering,2004, 22(4-5): 237-248.

    July 16, 2009, Revised March 23, 2010)

    * Project supported by the China National Funds for Distinguished Young Scientists (Grant No. 50925932), the National Key Technology R and D Program in the 11th Five Year Plan of China (Grant No. 2009BAB29B09), the State Key Laboratory of Hydroscience and Engineering, Tsinghua University (Grant No. 2008-ZY-5) and the National Key Basic Research Project of China (973 Program, Grant No. 2006CB403304).

    Biography:HUI Er-qing (1977-), Male, Ph. D.

    2010,22(3):329-337

    10.1016/S1001-6058(09)60062-7

    99国产综合亚洲精品| 国产成人一区二区三区免费视频网站 | 久久久国产欧美日韩av| 又黄又粗又硬又大视频| 中文字幕最新亚洲高清| av网站在线播放免费| 亚洲av成人精品一二三区| 免费高清在线观看视频在线观看| 中文字幕人妻丝袜一区二区| 国产淫语在线视频| 亚洲伊人色综图| 国产成人欧美| 秋霞在线观看毛片| a级毛片黄视频| 真人做人爱边吃奶动态| 制服人妻中文乱码| 午夜激情久久久久久久| 国产成人91sexporn| 婷婷色综合大香蕉| 国产精品 国内视频| 国产精品九九99| 国产精品香港三级国产av潘金莲 | 国产精品一区二区在线观看99| 精品国产一区二区三区四区第35| 亚洲国产精品成人久久小说| 亚洲成人手机| av天堂久久9| 狂野欧美激情性xxxx| 精品一区二区三区四区五区乱码 | 最黄视频免费看| 天天操日日干夜夜撸| 一级毛片黄色毛片免费观看视频| 久久精品熟女亚洲av麻豆精品| 欧美在线黄色| 热re99久久精品国产66热6| 啦啦啦啦在线视频资源| 欧美av亚洲av综合av国产av| av欧美777| 国产真人三级小视频在线观看| 男女午夜视频在线观看| 久久精品亚洲熟妇少妇任你| 19禁男女啪啪无遮挡网站| 日韩中文字幕欧美一区二区 | 黑丝袜美女国产一区| 九草在线视频观看| 人人妻人人爽人人添夜夜欢视频| 考比视频在线观看| 亚洲成人国产一区在线观看 | 少妇的丰满在线观看| 成人国语在线视频| 国产男人的电影天堂91| 大码成人一级视频| 国产成人啪精品午夜网站| 欧美日韩亚洲国产一区二区在线观看 | 亚洲熟女毛片儿| 中文乱码字字幕精品一区二区三区| 中文字幕高清在线视频| 国产精品 欧美亚洲| 久久久久国产一级毛片高清牌| 久久久国产一区二区| 熟女少妇亚洲综合色aaa.| 精品一区在线观看国产| 久久影院123| 夫妻午夜视频| 久久人人爽av亚洲精品天堂| 久久国产精品人妻蜜桃| 激情五月婷婷亚洲| 丰满人妻熟妇乱又伦精品不卡| 一级毛片电影观看| 欧美大码av| 日日夜夜操网爽| 亚洲精品av麻豆狂野| 少妇猛男粗大的猛烈进出视频| 啦啦啦视频在线资源免费观看| 久久人人爽av亚洲精品天堂| 亚洲综合色网址| svipshipincom国产片| 青青草视频在线视频观看| 黄色片一级片一级黄色片| 国产成人91sexporn| 国产亚洲欧美精品永久| 熟女少妇亚洲综合色aaa.| 婷婷丁香在线五月| 夫妻午夜视频| 一级毛片女人18水好多 | 免费少妇av软件| 精品视频人人做人人爽| 国产又爽黄色视频| 日韩 欧美 亚洲 中文字幕| 国产野战对白在线观看| 久久人人爽人人片av| 国产精品99久久99久久久不卡| 丝袜脚勾引网站| 好男人视频免费观看在线| 一级毛片电影观看| 欧美xxⅹ黑人| 国产精品.久久久| 午夜福利,免费看| 在线精品无人区一区二区三| 女人久久www免费人成看片| 久久久国产精品麻豆| 熟女av电影| 日本91视频免费播放| 国产亚洲av片在线观看秒播厂| 91成人精品电影| 免费黄频网站在线观看国产| 涩涩av久久男人的天堂| 亚洲国产看品久久| 日韩电影二区| 精品国产一区二区三区四区第35| 一区二区日韩欧美中文字幕| 久久亚洲国产成人精品v| 久热这里只有精品99| 男人添女人高潮全过程视频| 欧美黑人欧美精品刺激| 人人妻人人澡人人爽人人夜夜| 精品久久久久久久毛片微露脸 | 亚洲av成人不卡在线观看播放网 | www日本在线高清视频| 一边摸一边抽搐一进一出视频| 久久久久网色| 黄色视频在线播放观看不卡| 五月天丁香电影| 亚洲专区国产一区二区| 久久久久久人人人人人| 男人舔女人的私密视频| 中文字幕制服av| 狠狠精品人妻久久久久久综合| 亚洲精品一区蜜桃| 在线亚洲精品国产二区图片欧美| 亚洲自偷自拍图片 自拍| 性高湖久久久久久久久免费观看| 亚洲精品在线美女| 亚洲 国产 在线| 97在线人人人人妻| 999久久久国产精品视频| 欧美亚洲日本最大视频资源| 亚洲中文字幕日韩| 一边摸一边做爽爽视频免费| 亚洲三区欧美一区| 一二三四在线观看免费中文在| 久久久久国产精品人妻一区二区| 一边摸一边做爽爽视频免费| 久热爱精品视频在线9| a 毛片基地| 欧美精品啪啪一区二区三区 | 国产福利在线免费观看视频| 国产亚洲精品久久久久5区| 国产亚洲精品久久久久5区| 狠狠婷婷综合久久久久久88av| 美女脱内裤让男人舔精品视频| 午夜福利乱码中文字幕| 日本91视频免费播放| 国产精品久久久av美女十八| 日韩制服丝袜自拍偷拍| 好男人电影高清在线观看| 亚洲精品一区蜜桃| 每晚都被弄得嗷嗷叫到高潮| 国产欧美日韩一区二区三区在线| 麻豆乱淫一区二区| videos熟女内射| 国产有黄有色有爽视频| 国产成人av教育| av网站在线播放免费| 国产一区二区 视频在线| 亚洲国产精品一区三区| 青春草亚洲视频在线观看| 欧美国产精品一级二级三级| 国产精品免费大片| 免费av中文字幕在线| 亚洲国产精品一区三区| 国产99久久九九免费精品| 亚洲成人手机| 成年人午夜在线观看视频| 每晚都被弄得嗷嗷叫到高潮| 男人舔女人的私密视频| 欧美日韩亚洲高清精品| 晚上一个人看的免费电影| 久久精品国产a三级三级三级| 午夜影院在线不卡| 欧美老熟妇乱子伦牲交| 深夜精品福利| 午夜福利,免费看| 国产精品.久久久| 久久人妻福利社区极品人妻图片 | 欧美中文综合在线视频| 亚洲av男天堂| 高清av免费在线| 大片电影免费在线观看免费| www.熟女人妻精品国产| 桃花免费在线播放| 一级片'在线观看视频| 丰满迷人的少妇在线观看| 国产av国产精品国产| 精品久久久精品久久久| 老司机影院成人| 成人影院久久| 成人国产av品久久久| 电影成人av| 国产av精品麻豆| 国产一区二区三区综合在线观看| 香蕉丝袜av| 亚洲精品日韩在线中文字幕| 一级片免费观看大全| 一级,二级,三级黄色视频| 久久精品成人免费网站| 汤姆久久久久久久影院中文字幕| 国产熟女午夜一区二区三区| a级毛片在线看网站| 亚洲激情五月婷婷啪啪| 亚洲黑人精品在线| 欧美日韩视频精品一区| 高清av免费在线| 丁香六月欧美| 午夜激情av网站| 国产在线一区二区三区精| 国产精品国产av在线观看| 大陆偷拍与自拍| 如日韩欧美国产精品一区二区三区| 首页视频小说图片口味搜索 | 巨乳人妻的诱惑在线观看| 另类亚洲欧美激情| 久久中文字幕一级| 亚洲人成电影免费在线| av网站在线播放免费| 黄色a级毛片大全视频| 99国产精品一区二区三区| 亚洲国产欧美日韩在线播放| 午夜免费男女啪啪视频观看| 欧美激情高清一区二区三区| 巨乳人妻的诱惑在线观看| 久久精品aⅴ一区二区三区四区| 精品少妇黑人巨大在线播放| 日本a在线网址| 免费在线观看黄色视频的| 欧美另类一区| 亚洲国产欧美网| 男人爽女人下面视频在线观看| 亚洲国产欧美日韩在线播放| 国产精品麻豆人妻色哟哟久久| 欧美日本中文国产一区发布| 又大又爽又粗| av片东京热男人的天堂| 久久热在线av| 亚洲欧美中文字幕日韩二区| 十八禁高潮呻吟视频| 国产免费又黄又爽又色| 熟女av电影| 亚洲av日韩在线播放| 欧美成狂野欧美在线观看| 可以免费在线观看a视频的电影网站| 黄片播放在线免费| 水蜜桃什么品种好| 午夜福利影视在线免费观看| 国产欧美亚洲国产| 国产精品久久久av美女十八| 侵犯人妻中文字幕一二三四区| 免费看不卡的av| 亚洲欧美清纯卡通| 国产精品一区二区免费欧美 | 精品一区二区三区av网在线观看 | 欧美亚洲日本最大视频资源| 国产精品久久久av美女十八| 蜜桃在线观看..| 国产有黄有色有爽视频| 欧美日韩视频精品一区| xxx大片免费视频| 91老司机精品| 我要看黄色一级片免费的| 亚洲欧美中文字幕日韩二区| 黑人巨大精品欧美一区二区蜜桃| 啦啦啦啦在线视频资源| 亚洲专区中文字幕在线| 少妇被粗大的猛进出69影院| 午夜福利在线免费观看网站| 亚洲一区二区三区欧美精品| 中文欧美无线码| 亚洲综合色网址| 在线观看免费高清a一片| 国产精品免费视频内射| videos熟女内射| 老汉色av国产亚洲站长工具| av网站免费在线观看视频| 午夜老司机福利片| 两个人看的免费小视频| 一区二区三区乱码不卡18| 亚洲成人免费电影在线观看 | 日韩av不卡免费在线播放| 一区二区三区精品91| xxx大片免费视频| av一本久久久久| 人人澡人人妻人| 国语对白做爰xxxⅹ性视频网站| 亚洲五月色婷婷综合| www.自偷自拍.com| 亚洲专区中文字幕在线| 久久久精品区二区三区| 涩涩av久久男人的天堂| 欧美黑人精品巨大| 中文字幕制服av| 少妇被粗大的猛进出69影院| 韩国精品一区二区三区| 亚洲少妇的诱惑av| 日韩电影二区| 美女脱内裤让男人舔精品视频| 日韩一本色道免费dvd| 69精品国产乱码久久久| 日韩大片免费观看网站| 欧美久久黑人一区二区| av电影中文网址| 日韩中文字幕欧美一区二区 | 伊人亚洲综合成人网| 性少妇av在线| 校园人妻丝袜中文字幕| 免费在线观看视频国产中文字幕亚洲 | 精品人妻在线不人妻| 亚洲精品一二三| 亚洲欧美精品自产自拍| 日本色播在线视频| 久久久久精品国产欧美久久久 | 色网站视频免费| 久久99精品国语久久久| 免费不卡黄色视频| 青春草亚洲视频在线观看| 午夜福利影视在线免费观看| 亚洲av综合色区一区| 少妇被粗大的猛进出69影院| 亚洲人成77777在线视频| 欧美激情 高清一区二区三区| 男人舔女人的私密视频| 国产精品久久久人人做人人爽| 成年美女黄网站色视频大全免费| 最黄视频免费看| 丰满迷人的少妇在线观看| 美女国产高潮福利片在线看| 十八禁高潮呻吟视频| 岛国毛片在线播放| 只有这里有精品99| 免费久久久久久久精品成人欧美视频| 老司机深夜福利视频在线观看 | a 毛片基地| 国产免费又黄又爽又色| 国产在线视频一区二区| 大型av网站在线播放| 久久精品久久久久久噜噜老黄| 黑人欧美特级aaaaaa片| 九色亚洲精品在线播放| av欧美777| 天堂中文最新版在线下载| 精品久久久久久电影网| 手机成人av网站| 大香蕉久久网| 波多野结衣一区麻豆| 青春草亚洲视频在线观看| 国产精品免费视频内射| 精品人妻1区二区| 捣出白浆h1v1| 欧美成狂野欧美在线观看| 日韩中文字幕视频在线看片| 日日夜夜操网爽| 97人妻天天添夜夜摸| 亚洲少妇的诱惑av| 亚洲第一青青草原| 免费看十八禁软件| 一本—道久久a久久精品蜜桃钙片| 中文字幕av电影在线播放| 黄网站色视频无遮挡免费观看| 精品亚洲成国产av| 99国产精品免费福利视频| 国产黄色免费在线视频| 午夜久久久在线观看| 老司机影院成人| 亚洲少妇的诱惑av| 久久久国产精品麻豆| 亚洲av男天堂| 在线av久久热| 亚洲成人手机| 伊人久久大香线蕉亚洲五| av在线老鸭窝| 日本a在线网址| 国产精品人妻久久久影院| 国产精品熟女久久久久浪| 精品亚洲成国产av| 亚洲中文av在线| 欧美日韩福利视频一区二区| 欧美中文综合在线视频| 国产免费一区二区三区四区乱码| 亚洲专区国产一区二区| 青春草视频在线免费观看| 高清视频免费观看一区二区| 老鸭窝网址在线观看| 免费不卡黄色视频| 91精品伊人久久大香线蕉| 免费在线观看视频国产中文字幕亚洲 | 国产成人精品无人区| 91国产中文字幕| 熟女少妇亚洲综合色aaa.| 热re99久久国产66热| 啦啦啦啦在线视频资源| 中文字幕人妻丝袜一区二区| 99国产精品一区二区蜜桃av | 99精国产麻豆久久婷婷| 这个男人来自地球电影免费观看| 大片电影免费在线观看免费| 国产片特级美女逼逼视频| 伊人久久大香线蕉亚洲五| 久久久欧美国产精品| 老鸭窝网址在线观看| 亚洲精品国产色婷婷电影| 亚洲欧美精品综合一区二区三区| 久久亚洲精品不卡| 免费看十八禁软件| 亚洲国产中文字幕在线视频| 免费在线观看视频国产中文字幕亚洲 | 国产成人91sexporn| 色播在线永久视频| 欧美97在线视频| 纵有疾风起免费观看全集完整版| 久久久久视频综合| 一区二区三区乱码不卡18| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩亚洲综合一区二区三区_| 欧美日韩亚洲国产一区二区在线观看 | 波多野结衣av一区二区av| 欧美日韩亚洲综合一区二区三区_| 免费不卡黄色视频| 欧美精品亚洲一区二区| 美女午夜性视频免费| 肉色欧美久久久久久久蜜桃| 一级毛片 在线播放| 免费一级毛片在线播放高清视频 | 免费在线观看视频国产中文字幕亚洲 | 真人做人爱边吃奶动态| 天堂俺去俺来也www色官网| 国产一卡二卡三卡精品| 亚洲一卡2卡3卡4卡5卡精品中文| xxxhd国产人妻xxx| 国产欧美日韩一区二区三区在线| 人人妻人人澡人人爽人人夜夜| 日本午夜av视频| 午夜福利影视在线免费观看| 欧美黑人精品巨大| 99精品久久久久人妻精品| 好男人视频免费观看在线| tube8黄色片| 这个男人来自地球电影免费观看| 人成视频在线观看免费观看| 男女免费视频国产| 母亲3免费完整高清在线观看| 在线观看国产h片| 国产在线一区二区三区精| 老鸭窝网址在线观看| 国产99久久九九免费精品| 欧美在线黄色| 亚洲人成电影观看| 亚洲国产中文字幕在线视频| 9191精品国产免费久久| 免费看av在线观看网站| 国产av一区二区精品久久| 精品国产乱码久久久久久男人| 亚洲美女黄色视频免费看| 一级片'在线观看视频| 丁香六月欧美| 少妇 在线观看| 日韩欧美一区视频在线观看| 久久精品人人爽人人爽视色| 欧美精品一区二区免费开放| 亚洲国产中文字幕在线视频| 又黄又粗又硬又大视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲,欧美精品.| 女人久久www免费人成看片| 久久久亚洲精品成人影院| 亚洲伊人久久精品综合| 9色porny在线观看| 久久久久网色| 久久人人97超碰香蕉20202| 亚洲第一青青草原| 纵有疾风起免费观看全集完整版| 久久精品人人爽人人爽视色| 久久精品亚洲av国产电影网| 一级毛片电影观看| 国语对白做爰xxxⅹ性视频网站| av在线老鸭窝| 国产精品国产三级专区第一集| 黑人欧美特级aaaaaa片| 久久久久久人人人人人| 99精国产麻豆久久婷婷| 在线av久久热| 啦啦啦啦在线视频资源| 欧美激情高清一区二区三区| 久久午夜综合久久蜜桃| 免费在线观看影片大全网站 | 秋霞在线观看毛片| 国产一区二区三区综合在线观看| 亚洲成国产人片在线观看| 热re99久久精品国产66热6| 99re6热这里在线精品视频| 国产97色在线日韩免费| 亚洲五月色婷婷综合| 久久久国产精品麻豆| 一级片'在线观看视频| 妹子高潮喷水视频| 男女之事视频高清在线观看 | 一级毛片电影观看| 国产精品久久久久久精品古装| 精品久久久久久久毛片微露脸 | 一级毛片黄色毛片免费观看视频| 2018国产大陆天天弄谢| 国产精品免费视频内射| 久久久国产精品麻豆| 一本色道久久久久久精品综合| 天堂俺去俺来也www色官网| 国产精品二区激情视频| 爱豆传媒免费全集在线观看| 黄色视频在线播放观看不卡| 亚洲精品美女久久久久99蜜臀 | 亚洲熟女毛片儿| 亚洲,欧美,日韩| 大型av网站在线播放| 亚洲欧美日韩高清在线视频 | 纵有疾风起免费观看全集完整版| 黑人欧美特级aaaaaa片| 欧美在线黄色| 操出白浆在线播放| 老司机深夜福利视频在线观看 | 国产一区有黄有色的免费视频| 亚洲激情五月婷婷啪啪| 最近手机中文字幕大全| 99国产精品免费福利视频| 在线观看国产h片| 自线自在国产av| 69精品国产乱码久久久| 中文字幕精品免费在线观看视频| h视频一区二区三区| 91国产中文字幕| 女性被躁到高潮视频| 老汉色av国产亚洲站长工具| 国产男女内射视频| 亚洲精品一区蜜桃| 丝袜人妻中文字幕| 99久久99久久久精品蜜桃| av视频免费观看在线观看| 波多野结衣一区麻豆| 99久久99久久久精品蜜桃| 熟女av电影| 18禁裸乳无遮挡动漫免费视频| 啦啦啦啦在线视频资源| 欧美精品人与动牲交sv欧美| netflix在线观看网站| 乱人伦中国视频| 一级片免费观看大全| 99国产精品99久久久久| 亚洲第一青青草原| 国产亚洲欧美精品永久| 精品亚洲乱码少妇综合久久| 国产色视频综合| 在线 av 中文字幕| 90打野战视频偷拍视频| 亚洲精品国产色婷婷电影| 嫩草影视91久久| 一区二区三区精品91| 2021少妇久久久久久久久久久| 91精品三级在线观看| 成人18禁高潮啪啪吃奶动态图| 国产淫语在线视频| 欧美日韩福利视频一区二区| 欧美日韩精品网址| 国产亚洲精品久久久久5区| 精品高清国产在线一区| 午夜福利,免费看| svipshipincom国产片| 免费黄频网站在线观看国产| bbb黄色大片| 中文字幕制服av| 免费人妻精品一区二区三区视频| 亚洲欧美中文字幕日韩二区| 一二三四社区在线视频社区8| 搡老乐熟女国产| 精品一品国产午夜福利视频| 成年女人毛片免费观看观看9 | 成人18禁高潮啪啪吃奶动态图| 亚洲国产看品久久| 色婷婷av一区二区三区视频| 亚洲国产欧美日韩在线播放| 免费高清在线观看视频在线观看| 五月天丁香电影| 夫妻性生交免费视频一级片| 黄色毛片三级朝国网站| 亚洲精品美女久久久久99蜜臀 | 亚洲国产毛片av蜜桃av| 少妇人妻久久综合中文| 久久精品久久久久久久性| 你懂的网址亚洲精品在线观看| 丝袜美足系列| 一本一本久久a久久精品综合妖精| 精品一区二区三区四区五区乱码 | 久久天堂一区二区三区四区| 国产亚洲欧美精品永久| 青春草视频在线免费观看| 国产成人欧美在线观看 | 国产成人欧美| cao死你这个sao货| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产最新在线播放| 国产三级黄色录像| 校园人妻丝袜中文字幕| 黄片小视频在线播放| www.熟女人妻精品国产| 天天影视国产精品| 久久久国产一区二区| 在线观看一区二区三区激情| 菩萨蛮人人尽说江南好唐韦庄| 嫩草影视91久久| 精品国产一区二区久久|