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

    啦啦啦免费观看视频1| 一本精品99久久精品77| 久久亚洲真实| 亚洲性夜色夜夜综合| 国产一卡二卡三卡精品| 国产成人系列免费观看| 成人av一区二区三区在线看| 亚洲成人久久性| 人妻久久中文字幕网| 国产亚洲欧美在线一区二区| 波多野结衣av一区二区av| 亚洲专区国产一区二区| 欧美另类亚洲清纯唯美| 亚洲成av人片免费观看| 亚洲精品色激情综合| 亚洲自偷自拍图片 自拍| xxx96com| 日本在线视频免费播放| 日本一区二区免费在线视频| 青草久久国产| 老鸭窝网址在线观看| 最近最新中文字幕大全免费视频| 男女之事视频高清在线观看| 极品教师在线免费播放| 亚洲性夜色夜夜综合| a级毛片a级免费在线| 亚洲国产欧洲综合997久久, | 18禁裸乳无遮挡免费网站照片 | 成人免费观看视频高清| 神马国产精品三级电影在线观看 | 精品一区二区三区av网在线观看| 99在线视频只有这里精品首页| 最近最新中文字幕大全免费视频| 一二三四社区在线视频社区8| www日本黄色视频网| 少妇粗大呻吟视频| 亚洲国产日韩欧美精品在线观看 | 免费观看人在逋| 少妇被粗大的猛进出69影院| 成人国产一区最新在线观看| 久久欧美精品欧美久久欧美| 神马国产精品三级电影在线观看 | 成人国产一区最新在线观看| 此物有八面人人有两片| 亚洲精品一区av在线观看| 欧美日韩瑟瑟在线播放| 国产精品电影一区二区三区| 精品第一国产精品| 色av中文字幕| 亚洲avbb在线观看| 久99久视频精品免费| 在线观看免费日韩欧美大片| 一边摸一边抽搐一进一小说| 在线观看午夜福利视频| 性色av乱码一区二区三区2| 欧美日韩乱码在线| 亚洲av第一区精品v没综合| 18禁美女被吸乳视频| 亚洲七黄色美女视频| 大型av网站在线播放| 国内精品久久久久久久电影| 午夜精品久久久久久毛片777| 1024香蕉在线观看| 麻豆久久精品国产亚洲av| 久久国产精品影院| 色尼玛亚洲综合影院| 国产av又大| 丝袜人妻中文字幕| 男女做爰动态图高潮gif福利片| 一二三四在线观看免费中文在| 精品高清国产在线一区| 午夜久久久久精精品| 高清在线国产一区| 精品久久久久久久久久免费视频| 免费在线观看日本一区| 成人欧美大片| 黑人操中国人逼视频| 又黄又爽又免费观看的视频| 美女国产高潮福利片在线看| 国产一级毛片七仙女欲春2 | 久久人妻福利社区极品人妻图片| 大型av网站在线播放| 国产亚洲av嫩草精品影院| 亚洲五月婷婷丁香| 日本一区二区免费在线视频| 国产视频内射| 中文字幕最新亚洲高清| 在线国产一区二区在线| 免费无遮挡裸体视频| 日本免费一区二区三区高清不卡| 亚洲欧美日韩无卡精品| 亚洲中文av在线| 久久天躁狠狠躁夜夜2o2o| 国产色视频综合| 亚洲色图 男人天堂 中文字幕| 一个人观看的视频www高清免费观看 | 中文字幕久久专区| 亚洲精品国产一区二区精华液| aaaaa片日本免费| 中文字幕人妻熟女乱码| 久久九九热精品免费| 人成视频在线观看免费观看| 叶爱在线成人免费视频播放| 叶爱在线成人免费视频播放| 天天添夜夜摸| 免费高清视频大片| 黄色女人牲交| 国产成年人精品一区二区| 亚洲成人国产一区在线观看| 一a级毛片在线观看| 国产精品久久久久久精品电影 | 精品久久蜜臀av无| 观看免费一级毛片| 日日夜夜操网爽| 韩国av一区二区三区四区| 国产黄片美女视频| 日本免费一区二区三区高清不卡| 久久欧美精品欧美久久欧美| 久久香蕉精品热| 在线十欧美十亚洲十日本专区| 美国免费a级毛片| 午夜激情av网站| 黄色 视频免费看| 亚洲熟女毛片儿| 一区二区三区激情视频| 18禁观看日本| 中文字幕人成人乱码亚洲影| 国产91精品成人一区二区三区| 老汉色av国产亚洲站长工具| 久久久国产精品麻豆| 久久中文字幕人妻熟女| 久久久久久久精品吃奶| 亚洲色图 男人天堂 中文字幕| 国产精品亚洲av一区麻豆| 国产激情久久老熟女| 老汉色∧v一级毛片| 久久精品国产综合久久久| 麻豆成人av在线观看| 91老司机精品| 国产精品永久免费网站| 可以免费在线观看a视频的电影网站| 这个男人来自地球电影免费观看| 一级片免费观看大全| www.www免费av| 真人一进一出gif抽搐免费| 国产亚洲欧美98| 久久伊人香网站| 一级毛片精品| 曰老女人黄片| 日韩免费av在线播放| 亚洲第一av免费看| 亚洲精品美女久久久久99蜜臀| 19禁男女啪啪无遮挡网站| 国产亚洲精品久久久久5区| 国产精品久久电影中文字幕| 精品国产乱子伦一区二区三区| 天天一区二区日本电影三级| 国产一级毛片七仙女欲春2 | 久久中文看片网| 波多野结衣高清作品| 亚洲熟妇熟女久久| 亚洲狠狠婷婷综合久久图片| 亚洲人成77777在线视频| 久久中文字幕一级| 亚洲一码二码三码区别大吗| 一级片免费观看大全| 日韩大码丰满熟妇| 欧美成狂野欧美在线观看| 午夜激情福利司机影院| 老司机午夜十八禁免费视频| 人人妻人人澡人人看| 在线av久久热| 久久国产精品人妻蜜桃| 亚洲午夜理论影院| 国产乱人伦免费视频| 12—13女人毛片做爰片一| 999久久久精品免费观看国产| 亚洲国产精品sss在线观看| 色综合亚洲欧美另类图片| 免费看美女性在线毛片视频| 久久午夜亚洲精品久久| 91麻豆av在线| 成人国语在线视频| 久久久久久久久免费视频了| 久久精品国产99精品国产亚洲性色| 欧美zozozo另类| 午夜福利在线观看吧| 精品一区二区三区四区五区乱码| √禁漫天堂资源中文www| 亚洲一区高清亚洲精品| 制服人妻中文乱码| 久久精品国产99精品国产亚洲性色| 国产精品久久电影中文字幕| 一级毛片高清免费大全| 婷婷精品国产亚洲av| 欧美在线黄色| 亚洲精品粉嫩美女一区| 国产精品久久视频播放| 久久久国产成人精品二区| 中文在线观看免费www的网站 | 久久人妻福利社区极品人妻图片| 日韩视频一区二区在线观看| 免费高清在线观看日韩| 深夜精品福利| 亚洲精华国产精华精| 51午夜福利影视在线观看| 操出白浆在线播放| 成人手机av| 欧美另类亚洲清纯唯美| 桃色一区二区三区在线观看| 亚洲国产中文字幕在线视频| 啦啦啦观看免费观看视频高清| 不卡av一区二区三区| 日本免费一区二区三区高清不卡| 亚洲激情在线av| 51午夜福利影视在线观看| 午夜老司机福利片| 黑人欧美特级aaaaaa片| 国产精品乱码一区二三区的特点| 少妇粗大呻吟视频| 国产主播在线观看一区二区| 90打野战视频偷拍视频| 国产成人精品久久二区二区91| 国产成人av教育| 日日夜夜操网爽| 国产亚洲av高清不卡| 亚洲国产精品999在线| 国产真人三级小视频在线观看| 久久中文字幕一级| 亚洲av第一区精品v没综合| 日韩欧美一区二区三区在线观看| 18美女黄网站色大片免费观看| 国产精品永久免费网站| 国产黄片美女视频| 国产av一区二区精品久久| 少妇裸体淫交视频免费看高清 | 男女做爰动态图高潮gif福利片| 亚洲av美国av| 成人一区二区视频在线观看| 婷婷亚洲欧美| 男人舔女人的私密视频| 国产av一区二区精品久久| 长腿黑丝高跟| 久久久久久久精品吃奶| 欧美日韩一级在线毛片| 中出人妻视频一区二区| 中文字幕最新亚洲高清| 久久久久久久久中文| av免费在线观看网站| 亚洲专区中文字幕在线| avwww免费| videosex国产| 黄片大片在线免费观看| 欧美中文日本在线观看视频| xxx96com| 日韩精品中文字幕看吧| 高潮久久久久久久久久久不卡| 国产精品国产高清国产av| 亚洲精品在线观看二区| 国产精品精品国产色婷婷| 中出人妻视频一区二区| 久热爱精品视频在线9| 国产成年人精品一区二区| 免费在线观看成人毛片| 国产精品一区二区免费欧美| 午夜两性在线视频| 国产精品香港三级国产av潘金莲| 90打野战视频偷拍视频| 国产熟女午夜一区二区三区| 级片在线观看| 757午夜福利合集在线观看| 亚洲国产精品999在线| 国产精品精品国产色婷婷| 国产麻豆成人av免费视频| 看免费av毛片| 国产高清视频在线播放一区| 国产欧美日韩精品亚洲av| 99久久无色码亚洲精品果冻| 女人被狂操c到高潮| 久久久久久久久免费视频了| 国产精品免费视频内射| 美女免费视频网站| 51午夜福利影视在线观看| 国产熟女午夜一区二区三区| 亚洲人成网站高清观看| 欧美乱色亚洲激情| 国产亚洲精品久久久久久毛片| 97人妻精品一区二区三区麻豆 | 哪里可以看免费的av片| 一个人观看的视频www高清免费观看 | 最近最新免费中文字幕在线| 欧美黄色淫秽网站| 黑人操中国人逼视频| 国产精品久久久久久亚洲av鲁大| 久久精品影院6| 91国产中文字幕| 亚洲在线自拍视频| 午夜两性在线视频| 亚洲成人国产一区在线观看| 在线观看免费日韩欧美大片| 国产av一区在线观看免费| 欧美精品亚洲一区二区| 夜夜看夜夜爽夜夜摸| 久久精品国产99精品国产亚洲性色| 国产片内射在线| netflix在线观看网站| 岛国在线观看网站| 999久久久精品免费观看国产| 一级a爱片免费观看的视频| 国产免费av片在线观看野外av| 天天添夜夜摸| 日韩一卡2卡3卡4卡2021年| 精品卡一卡二卡四卡免费| 日韩精品青青久久久久久| 欧美激情极品国产一区二区三区| 国产国语露脸激情在线看| 国产亚洲精品一区二区www| 精品电影一区二区在线| 听说在线观看完整版免费高清| 欧美日韩一级在线毛片| 丁香六月欧美| 成人国产一区最新在线观看| 日韩视频一区二区在线观看| 亚洲成a人片在线一区二区| 无限看片的www在线观看| 97碰自拍视频| cao死你这个sao货| 欧美国产日韩亚洲一区| 午夜视频精品福利| 亚洲熟妇熟女久久| 一本一本综合久久| 高潮久久久久久久久久久不卡| 中文亚洲av片在线观看爽| 亚洲片人在线观看| 99久久无色码亚洲精品果冻| 黄色片一级片一级黄色片| 免费一级毛片在线播放高清视频| 国产99白浆流出| 黄片小视频在线播放| 99riav亚洲国产免费| 精品久久久久久久久久久久久 | 欧美日韩黄片免| 男女床上黄色一级片免费看| 国产亚洲av嫩草精品影院| 最新美女视频免费是黄的| 老司机在亚洲福利影院| 久久久国产欧美日韩av| 国内毛片毛片毛片毛片毛片| 亚洲欧美激情综合另类| 亚洲中文日韩欧美视频| 18禁观看日本| 女性生殖器流出的白浆| 国产麻豆成人av免费视频| 亚洲精品av麻豆狂野| 婷婷精品国产亚洲av在线| 久久人妻福利社区极品人妻图片| 精品久久蜜臀av无| 亚洲午夜精品一区,二区,三区| 国产在线观看jvid| 午夜福利在线观看吧| 人妻丰满熟妇av一区二区三区| 久久香蕉精品热| 黄片小视频在线播放| 精品免费久久久久久久清纯| 国产精品亚洲av一区麻豆| 日韩大尺度精品在线看网址| 男人舔奶头视频| 国内毛片毛片毛片毛片毛片| 国产成人欧美在线观看| 国内精品久久久久精免费| 国产午夜福利久久久久久| 91成年电影在线观看| 99热这里只有精品一区 | 亚洲av电影不卡..在线观看| 99国产精品99久久久久| 欧美日本亚洲视频在线播放| 99久久无色码亚洲精品果冻| xxx96com| av天堂在线播放| 色播亚洲综合网| 精品久久久久久久久久免费视频| 欧美+亚洲+日韩+国产| 日本免费a在线| 日本免费a在线| 俺也久久电影网| 国产成人精品久久二区二区免费| 亚洲熟女毛片儿| 亚洲av电影在线进入| 制服诱惑二区| 亚洲人成网站高清观看| 亚洲av五月六月丁香网| 午夜影院日韩av| 久久午夜亚洲精品久久| 又黄又爽又免费观看的视频| 最好的美女福利视频网| 又大又爽又粗| 日日摸夜夜添夜夜添小说| 国产一级毛片七仙女欲春2 | 性欧美人与动物交配| 国产精品99久久99久久久不卡| 国产精品一区二区三区四区久久 | 欧美一级a爱片免费观看看 | 精品电影一区二区在线| 性色av乱码一区二区三区2| 中文字幕最新亚洲高清| 欧美又色又爽又黄视频| 狠狠狠狠99中文字幕| 高清在线国产一区| 亚洲专区中文字幕在线| 国产成人精品久久二区二区91| 成人欧美大片| 久久精品国产清高在天天线| 午夜精品在线福利| 日日爽夜夜爽网站| 丁香欧美五月| 色播亚洲综合网| 制服人妻中文乱码| 999精品在线视频| 欧美日韩中文字幕国产精品一区二区三区| 美女 人体艺术 gogo| 国产精品亚洲美女久久久| 国内揄拍国产精品人妻在线 | www国产在线视频色| 美女国产高潮福利片在线看| 亚洲成人国产一区在线观看| 可以在线观看的亚洲视频| 两性午夜刺激爽爽歪歪视频在线观看 | 可以在线观看毛片的网站| 久久久久亚洲av毛片大全| 亚洲七黄色美女视频| √禁漫天堂资源中文www| 国产精品亚洲一级av第二区| 欧美激情极品国产一区二区三区| 亚洲成国产人片在线观看| 亚洲一区中文字幕在线| 88av欧美| 日本撒尿小便嘘嘘汇集6| 国产高清有码在线观看视频 | 久久久久久九九精品二区国产 | 久久精品国产亚洲av高清一级| 丰满人妻熟妇乱又伦精品不卡| 好男人电影高清在线观看| 最近最新免费中文字幕在线| 一区二区三区国产精品乱码| 久热这里只有精品99| 好男人在线观看高清免费视频 | 久久久精品欧美日韩精品| 久久午夜综合久久蜜桃| 日韩 欧美 亚洲 中文字幕| 老鸭窝网址在线观看| 伦理电影免费视频| 女人高潮潮喷娇喘18禁视频| 精品乱码久久久久久99久播| 久久草成人影院| 激情在线观看视频在线高清| 国产主播在线观看一区二区| 午夜精品久久久久久毛片777| 真人做人爱边吃奶动态| 午夜免费观看网址| 国产97色在线日韩免费| 男人舔奶头视频| 国产真实乱freesex| 在线观看www视频免费| 淫秽高清视频在线观看| 欧美日韩乱码在线| 日日爽夜夜爽网站| 免费电影在线观看免费观看| 在线播放国产精品三级| 亚洲国产精品sss在线观看| 99久久精品国产亚洲精品| 男人舔奶头视频| bbb黄色大片| 老司机深夜福利视频在线观看| 国产一区二区三区在线臀色熟女| www.熟女人妻精品国产| av在线播放免费不卡| 精品无人区乱码1区二区| 久久婷婷成人综合色麻豆| 久久久国产成人免费| 黄色 视频免费看| 欧美乱妇无乱码| 亚洲欧美精品综合一区二区三区| 久久人妻福利社区极品人妻图片| 中文字幕av电影在线播放| 亚洲男人天堂网一区| xxxwww97欧美| 少妇熟女aⅴ在线视频| 啦啦啦 在线观看视频| 50天的宝宝边吃奶边哭怎么回事| 精品久久久久久,| 亚洲第一欧美日韩一区二区三区| 日韩成人在线观看一区二区三区| 在线av久久热| 亚洲成人精品中文字幕电影| 国产成人欧美在线观看| 亚洲国产精品sss在线观看| 夜夜夜夜夜久久久久| 婷婷精品国产亚洲av在线| 久久久久久久久久黄片| 久久久国产成人精品二区| 亚洲中文字幕一区二区三区有码在线看 | 两个人免费观看高清视频| 人妻久久中文字幕网| 国产精品影院久久| 无人区码免费观看不卡| 不卡一级毛片| 国产高清有码在线观看视频 | 日本撒尿小便嘘嘘汇集6| 欧美日韩瑟瑟在线播放| 制服丝袜大香蕉在线| 国产片内射在线| 久久久久久人人人人人| 丝袜美腿诱惑在线| 亚洲国产精品久久男人天堂| 精品国产亚洲在线| 午夜福利免费观看在线| 久久久久国产一级毛片高清牌| 国语自产精品视频在线第100页| 日韩欧美三级三区| 午夜老司机福利片| 女性被躁到高潮视频| 黄色视频不卡| 中文字幕av电影在线播放| 久久久水蜜桃国产精品网| 精品久久久久久久人妻蜜臀av| 色综合站精品国产| 亚洲专区字幕在线| 狠狠狠狠99中文字幕| 国产精品一区二区精品视频观看| 久久久久国内视频| 一区二区三区激情视频| 国内毛片毛片毛片毛片毛片| 色尼玛亚洲综合影院| 日韩免费av在线播放| а√天堂www在线а√下载| xxx96com| 国产野战对白在线观看| 亚洲国产欧美一区二区综合| 精品久久久久久久久久免费视频| 老熟妇仑乱视频hdxx| 18禁美女被吸乳视频| 国产在线精品亚洲第一网站| 亚洲精品一区av在线观看| 1024视频免费在线观看| 一级毛片高清免费大全| 午夜a级毛片| 中文字幕人妻熟女乱码| 757午夜福利合集在线观看| 亚洲,欧美精品.| 亚洲成a人片在线一区二区| 国产精品免费视频内射| 亚洲无线在线观看| 国产激情久久老熟女| 午夜视频精品福利| 亚洲黑人精品在线| 国产99久久九九免费精品| 午夜福利欧美成人| 可以在线观看的亚洲视频| 亚洲午夜精品一区,二区,三区| 在线天堂中文资源库| 午夜福利成人在线免费观看| 老司机靠b影院| 满18在线观看网站| 成人欧美大片| 日韩欧美一区二区三区在线观看| x7x7x7水蜜桃| 精品国产亚洲在线| 满18在线观看网站| 真人做人爱边吃奶动态| 国产蜜桃级精品一区二区三区| 午夜激情福利司机影院| 777久久人妻少妇嫩草av网站| 久久久久久大精品| 最新美女视频免费是黄的| 欧美av亚洲av综合av国产av| 久久久国产成人免费| 国产精品一区二区精品视频观看| 一个人免费在线观看的高清视频| 国产高清videossex| 久久人人精品亚洲av| 免费看十八禁软件| 日韩国内少妇激情av| 国产精品日韩av在线免费观看| 级片在线观看| 啦啦啦免费观看视频1| 午夜激情福利司机影院| 久久精品91蜜桃| 色播亚洲综合网| 91九色精品人成在线观看| 国产精品一区二区三区四区久久 | 亚洲人成网站高清观看| avwww免费| 国产精品亚洲av一区麻豆| 午夜福利高清视频| 亚洲五月婷婷丁香| 成年版毛片免费区| 免费看美女性在线毛片视频| 可以在线观看毛片的网站| 1024香蕉在线观看| 国产亚洲欧美精品永久| 欧美日本亚洲视频在线播放| 亚洲人成电影免费在线| 精品一区二区三区视频在线观看免费| 男女视频在线观看网站免费 | 午夜久久久久精精品| 国产真人三级小视频在线观看| aaaaa片日本免费| 亚洲专区中文字幕在线| 亚洲色图 男人天堂 中文字幕| 在线观看66精品国产| 免费搜索国产男女视频| 欧美性长视频在线观看| 国产精品av久久久久免费| 男男h啪啪无遮挡| 亚洲精品粉嫩美女一区|