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

    精品久久久久久久毛片微露脸| 国产高清有码在线观看视频 | 亚洲自拍偷在线| 黄片小视频在线播放| 精品国产乱码久久久久久男人| 欧美乱妇无乱码| 国内久久婷婷六月综合欲色啪| 久久精品亚洲精品国产色婷小说| 亚洲成人精品中文字幕电影| 亚洲精品中文字幕一二三四区| 老司机午夜十八禁免费视频| 一级毛片女人18水好多| 精品一区二区三区视频在线观看免费| 亚洲av美国av| www.www免费av| 极品人妻少妇av视频| 亚洲全国av大片| 91成人精品电影| 亚洲 欧美 日韩 在线 免费| 99久久久亚洲精品蜜臀av| 久久国产乱子伦精品免费另类| 国产一卡二卡三卡精品| 天堂√8在线中文| 91av网站免费观看| 香蕉国产在线看| 97超级碰碰碰精品色视频在线观看| 亚洲免费av在线视频| av福利片在线| 久久久久国产一级毛片高清牌| 成熟少妇高潮喷水视频| 伦理电影免费视频| 国产亚洲精品一区二区www| 国产精品亚洲一级av第二区| 久久欧美精品欧美久久欧美| 午夜a级毛片| 在线天堂中文资源库| 亚洲av第一区精品v没综合| 久久香蕉激情| 成人三级做爰电影| 色播在线永久视频| 一二三四在线观看免费中文在| 午夜福利高清视频| 女人高潮潮喷娇喘18禁视频| 日韩欧美一区二区三区在线观看| АⅤ资源中文在线天堂| 叶爱在线成人免费视频播放| 视频在线观看一区二区三区| 欧美中文日本在线观看视频| 一级毛片女人18水好多| www日本在线高清视频| 麻豆久久精品国产亚洲av| 日本五十路高清| 天天躁狠狠躁夜夜躁狠狠躁| 久久狼人影院| 久久精品国产清高在天天线| 亚洲中文av在线| 日韩av在线大香蕉| 99在线视频只有这里精品首页| 精品日产1卡2卡| 女人被狂操c到高潮| 日韩 欧美 亚洲 中文字幕| 国产真人三级小视频在线观看| 美女国产高潮福利片在线看| 精品国产国语对白av| 久久青草综合色| 黄色a级毛片大全视频| 欧美乱色亚洲激情| 亚洲人成网站在线播放欧美日韩| av超薄肉色丝袜交足视频| 又大又爽又粗| 精品无人区乱码1区二区| 国产激情久久老熟女| 中文字幕av电影在线播放| www日本在线高清视频| 国产欧美日韩一区二区三区在线| 男人舔女人下体高潮全视频| 欧美激情高清一区二区三区| 香蕉国产在线看| 国内精品久久久久久久电影| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品一区av在线观看| 天天添夜夜摸| 国产精品久久电影中文字幕| 国产亚洲av高清不卡| 日韩精品青青久久久久久| 国产在线观看jvid| 啦啦啦 在线观看视频| 一级片免费观看大全| bbb黄色大片| 日本在线视频免费播放| 又黄又粗又硬又大视频| 丁香欧美五月| www.www免费av| 天堂影院成人在线观看| 一区二区三区高清视频在线| 日韩欧美免费精品| 国产av一区在线观看免费| 日韩中文字幕欧美一区二区| 女性生殖器流出的白浆| 欧美日韩亚洲国产一区二区在线观看| 日韩中文字幕欧美一区二区| 成在线人永久免费视频| 欧洲精品卡2卡3卡4卡5卡区| 九色亚洲精品在线播放| 亚洲国产毛片av蜜桃av| 母亲3免费完整高清在线观看| 老司机午夜福利在线观看视频| 人妻久久中文字幕网| av在线天堂中文字幕| 一级,二级,三级黄色视频| 久99久视频精品免费| 97超级碰碰碰精品色视频在线观看| 免费观看精品视频网站| 国产精品亚洲一级av第二区| 国产亚洲精品第一综合不卡| 亚洲精品在线观看二区| 欧美黄色片欧美黄色片| 最近最新免费中文字幕在线| 国内精品久久久久久久电影| 午夜精品久久久久久毛片777| 日韩大码丰满熟妇| 99riav亚洲国产免费| 亚洲免费av在线视频| 色播在线永久视频| 亚洲专区中文字幕在线| 长腿黑丝高跟| 桃红色精品国产亚洲av| 99久久精品国产亚洲精品| 啦啦啦观看免费观看视频高清 | 精品国产亚洲在线| 在线观看66精品国产| 亚洲自偷自拍图片 自拍| 国产1区2区3区精品| 午夜激情av网站| 在线视频色国产色| 免费高清视频大片| cao死你这个sao货| 亚洲av成人av| 999久久久精品免费观看国产| 久久久久久久久免费视频了| 亚洲人成电影观看| 18禁国产床啪视频网站| 国产亚洲av高清不卡| 亚洲精品国产色婷婷电影| 久久亚洲真实| av欧美777| 在线十欧美十亚洲十日本专区| 桃色一区二区三区在线观看| 天天添夜夜摸| 最好的美女福利视频网| 女人爽到高潮嗷嗷叫在线视频| 亚洲人成伊人成综合网2020| av有码第一页| 免费在线观看影片大全网站| 一级毛片女人18水好多| 欧美人与性动交α欧美精品济南到| 欧美黑人欧美精品刺激| 亚洲精品中文字幕在线视频| 人人澡人人妻人| 成人手机av| 国产精品二区激情视频| 老汉色∧v一级毛片| 国产成人av激情在线播放| 青草久久国产| 日韩三级视频一区二区三区| 12—13女人毛片做爰片一| 看免费av毛片| 国产av又大| 亚洲精品粉嫩美女一区| 欧美+亚洲+日韩+国产| 色综合婷婷激情| 成人18禁高潮啪啪吃奶动态图| 国产亚洲精品av在线| 麻豆久久精品国产亚洲av| 啦啦啦观看免费观看视频高清 | 国产成+人综合+亚洲专区| 首页视频小说图片口味搜索| 久久中文字幕一级| 成人三级黄色视频| 色老头精品视频在线观看| 国产激情久久老熟女| 一边摸一边抽搐一进一小说| 91麻豆精品激情在线观看国产| 日本黄色视频三级网站网址| 我的亚洲天堂| 一进一出抽搐动态| 免费观看人在逋| av在线播放免费不卡| 久久国产精品影院| 亚洲在线自拍视频| 99香蕉大伊视频| 免费高清在线观看日韩| 村上凉子中文字幕在线| 天天一区二区日本电影三级 | 黄色 视频免费看| 精品卡一卡二卡四卡免费| 性欧美人与动物交配| cao死你这个sao货| 欧美黄色片欧美黄色片| www.自偷自拍.com| 又黄又粗又硬又大视频| 午夜成年电影在线免费观看| 青草久久国产| 免费高清在线观看日韩| 男人舔女人下体高潮全视频| 午夜两性在线视频| 91精品国产国语对白视频| 欧美中文日本在线观看视频| 午夜福利在线观看吧| 国产黄a三级三级三级人| 日韩成人在线观看一区二区三区| 9热在线视频观看99| xxx96com| 免费看av在线观看网站| 成人无遮挡网站| 精品久久国产蜜桃| 亚洲精品日韩av片在线观看| 久久香蕉精品热| 免费av观看视频| 精品午夜福利视频在线观看一区| 日韩欧美精品免费久久| 床上黄色一级片| 少妇人妻一区二区三区视频| 国产亚洲精品久久久久久毛片| av专区在线播放| 免费看日本二区| 国产精华一区二区三区| 狂野欧美激情性xxxx在线观看| 亚洲精品色激情综合| 国内精品宾馆在线| 天堂网av新在线| 在线播放无遮挡| 欧美不卡视频在线免费观看| .国产精品久久| 人妻夜夜爽99麻豆av| 午夜福利欧美成人| 亚洲真实伦在线观看| 成人特级黄色片久久久久久久| 搞女人的毛片| 免费一级毛片在线播放高清视频| 欧美日韩亚洲国产一区二区在线观看| 国产人妻一区二区三区在| 欧美激情在线99| 亚洲无线在线观看| 黄色女人牲交| 亚洲美女黄片视频| 12—13女人毛片做爰片一| 亚洲色图av天堂| 精品人妻偷拍中文字幕| 99在线人妻在线中文字幕| 一个人观看的视频www高清免费观看| 嫁个100分男人电影在线观看| 熟女人妻精品中文字幕| 女的被弄到高潮叫床怎么办 | 人人妻,人人澡人人爽秒播| 舔av片在线| 精品久久国产蜜桃| 久久久久久久久久黄片| 精品国产三级普通话版| av中文乱码字幕在线| 12—13女人毛片做爰片一| 午夜福利视频1000在线观看| 制服丝袜大香蕉在线| 成年女人永久免费观看视频| 色吧在线观看| 精品久久久久久,| 成人国产一区最新在线观看| 一区二区三区激情视频| 看免费成人av毛片| 欧美日本视频| 国产午夜精品久久久久久一区二区三区 | 国产熟女欧美一区二区| 伊人久久精品亚洲午夜| 中文字幕av在线有码专区| 国产高清三级在线| 亚洲国产精品久久男人天堂| 男人狂女人下面高潮的视频| 欧美xxxx性猛交bbbb| 国产精品嫩草影院av在线观看 | 成人av在线播放网站| 两个人视频免费观看高清| 国产精品自产拍在线观看55亚洲| 天堂√8在线中文| 亚洲av.av天堂| 久久亚洲精品不卡| 无遮挡黄片免费观看| 国产精品综合久久久久久久免费| 国产熟女欧美一区二区| 国产精品98久久久久久宅男小说| 伊人久久精品亚洲午夜| 日日摸夜夜添夜夜添av毛片 | 国产高清激情床上av| 18禁黄网站禁片免费观看直播| av中文乱码字幕在线| 男女啪啪激烈高潮av片| 国产精品一区二区性色av| 日韩中文字幕欧美一区二区| 亚洲专区国产一区二区| 亚洲va日本ⅴa欧美va伊人久久| 中文资源天堂在线| 在线播放无遮挡| 色综合婷婷激情| 真人做人爱边吃奶动态| 国产欧美日韩精品一区二区| 成人综合一区亚洲| 真实男女啪啪啪动态图| 老女人水多毛片| 久久久久精品国产欧美久久久| 亚洲不卡免费看| 亚洲国产精品久久男人天堂| 欧美高清性xxxxhd video| 校园春色视频在线观看| 欧美极品一区二区三区四区| 久久精品国产99精品国产亚洲性色| aaaaa片日本免费| 亚洲成a人片在线一区二区| 在线免费观看不下载黄p国产 | 国产精品98久久久久久宅男小说| 成年版毛片免费区| 婷婷亚洲欧美| 免费av不卡在线播放| 在线观看免费视频日本深夜| 在线天堂最新版资源| 亚洲中文字幕日韩| 精品一区二区三区视频在线观看免费| 九九在线视频观看精品| 熟女电影av网| 国产精品久久久久久久电影| 动漫黄色视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 99视频精品全部免费 在线| 熟女电影av网| 国产精品一区二区性色av| 一区二区三区高清视频在线| 久久热精品热| 亚洲乱码一区二区免费版| 亚洲精品成人久久久久久| 男人舔女人下体高潮全视频| 国产伦精品一区二区三区视频9| 日日撸夜夜添| 久久人人爽人人爽人人片va| 亚洲熟妇熟女久久| 日韩强制内射视频| 俄罗斯特黄特色一大片| 国内毛片毛片毛片毛片毛片| 两个人视频免费观看高清| 最近最新免费中文字幕在线| 午夜免费激情av| 国产真实乱freesex| 国产日本99.免费观看| 精品午夜福利视频在线观看一区| 蜜桃久久精品国产亚洲av| 老女人水多毛片| 可以在线观看毛片的网站| 国产亚洲精品综合一区在线观看| 国产精品99久久久久久久久| 精品久久久久久久末码| 欧美日韩亚洲国产一区二区在线观看| 日本黄大片高清| 亚洲熟妇熟女久久| 国产亚洲精品久久久com| 久久人妻av系列| 精品不卡国产一区二区三区| 99久久精品热视频| 69人妻影院| 免费观看的影片在线观看| 无人区码免费观看不卡| 九九爱精品视频在线观看| 欧美最黄视频在线播放免费| 伦精品一区二区三区| 波野结衣二区三区在线| 午夜福利高清视频| 少妇丰满av| 哪里可以看免费的av片| 国内毛片毛片毛片毛片毛片| 97热精品久久久久久| 国产麻豆成人av免费视频| 中文字幕免费在线视频6| 中文字幕久久专区| 天天一区二区日本电影三级| 国产精品av视频在线免费观看| 久久午夜福利片| 久久久久精品国产欧美久久久| 免费无遮挡裸体视频| 久久精品夜夜夜夜夜久久蜜豆| 丝袜美腿在线中文| 亚洲男人的天堂狠狠| 婷婷亚洲欧美| 亚洲国产欧洲综合997久久,| 国产黄色小视频在线观看| 日本欧美国产在线视频| 91精品国产九色| 99热精品在线国产| 成年版毛片免费区| 又粗又爽又猛毛片免费看| 国产伦精品一区二区三区视频9| 国国产精品蜜臀av免费| 欧美成人a在线观看| 国内精品久久久久久久电影| 精品人妻一区二区三区麻豆 | 嫩草影院新地址| 在线免费十八禁| 精品久久久久久成人av| 国产激情偷乱视频一区二区| 成年版毛片免费区| 91午夜精品亚洲一区二区三区 | 午夜激情福利司机影院| 亚洲天堂国产精品一区在线| 国产成人一区二区在线| 永久网站在线| 男女下面进入的视频免费午夜| 女生性感内裤真人,穿戴方法视频| 少妇猛男粗大的猛烈进出视频 | 国产欧美日韩精品亚洲av| 久久久久久久精品吃奶| 人妻少妇偷人精品九色| 日本三级黄在线观看| 在线国产一区二区在线| 成人特级黄色片久久久久久久| 97人妻精品一区二区三区麻豆| 长腿黑丝高跟| 精品久久久久久久末码| 免费看av在线观看网站| 日日夜夜操网爽| 国产熟女欧美一区二区| 国产一区二区在线观看日韩| 成人特级av手机在线观看| a在线观看视频网站| 精品日产1卡2卡| netflix在线观看网站| av.在线天堂| 久久婷婷人人爽人人干人人爱| 欧美性猛交╳xxx乱大交人| 在线免费十八禁| 久久人妻av系列| 亚洲av五月六月丁香网| 婷婷精品国产亚洲av在线| 哪里可以看免费的av片| 一本久久中文字幕| 欧美国产日韩亚洲一区| 日本熟妇午夜| 99国产极品粉嫩在线观看| 国产国拍精品亚洲av在线观看| www.www免费av| 亚洲av电影不卡..在线观看| 国产伦在线观看视频一区| 亚洲 国产 在线| 国产国拍精品亚洲av在线观看| 嫁个100分男人电影在线观看| 午夜福利18| 亚洲国产欧洲综合997久久,| 欧美zozozo另类| 少妇人妻精品综合一区二区 | 一本精品99久久精品77| 老熟妇乱子伦视频在线观看| 自拍偷自拍亚洲精品老妇| 琪琪午夜伦伦电影理论片6080| 真人一进一出gif抽搐免费| 天堂影院成人在线观看| 国产精品1区2区在线观看.| 黄片wwwwww| 一个人免费在线观看电影| 久久久久性生活片| 日韩欧美在线乱码| 男人舔奶头视频| 欧美一区二区国产精品久久精品| 99在线视频只有这里精品首页| 观看免费一级毛片| 日本黄色视频三级网站网址| 亚洲av.av天堂| 久久久久九九精品影院| 日韩欧美在线二视频| 国产乱人伦免费视频| 成人国产麻豆网| 毛片一级片免费看久久久久 | 成人美女网站在线观看视频| 日本黄色视频三级网站网址| 日韩欧美三级三区| 国产精华一区二区三区| 国产精品一区二区免费欧美| 国产午夜福利久久久久久| 国产精品国产三级国产av玫瑰| 亚洲va在线va天堂va国产| 久久国产精品人妻蜜桃| 日本-黄色视频高清免费观看| 国产精品免费一区二区三区在线| av女优亚洲男人天堂| 精品一区二区三区人妻视频| 嫁个100分男人电影在线观看| 精品不卡国产一区二区三区| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品久久久久久毛片| 桃红色精品国产亚洲av| 黄色欧美视频在线观看| 身体一侧抽搐| 成人特级av手机在线观看| xxxwww97欧美| 亚洲欧美日韩东京热| 久久精品夜夜夜夜夜久久蜜豆| 亚洲五月天丁香| 男女那种视频在线观看| 一区二区三区高清视频在线| 久久久久久大精品| 级片在线观看| 国产伦精品一区二区三区视频9| 一本一本综合久久| 国产乱人伦免费视频| 乱人视频在线观看| 国产色爽女视频免费观看| 亚洲成人中文字幕在线播放| 可以在线观看的亚洲视频| 欧美区成人在线视频| 99国产精品一区二区蜜桃av| 精品久久久久久久久久免费视频| 日韩精品青青久久久久久| 91在线精品国自产拍蜜月| 又爽又黄a免费视频| 一进一出抽搐gif免费好疼| av在线天堂中文字幕| 内地一区二区视频在线| 三级国产精品欧美在线观看| 噜噜噜噜噜久久久久久91| 欧美另类亚洲清纯唯美| 国产极品精品免费视频能看的| 日韩强制内射视频| 亚洲久久久久久中文字幕| 国内精品美女久久久久久| 成人特级av手机在线观看| 午夜福利成人在线免费观看| av.在线天堂| 99九九线精品视频在线观看视频| 欧美高清性xxxxhd video| 国产精品一区二区性色av| 亚洲aⅴ乱码一区二区在线播放| 日韩国内少妇激情av| 俄罗斯特黄特色一大片| 国产激情偷乱视频一区二区| 五月玫瑰六月丁香| 黄色丝袜av网址大全| 国产一区二区在线观看日韩| 日本撒尿小便嘘嘘汇集6| 日本黄大片高清| 最后的刺客免费高清国语| 黄色配什么色好看| 91在线精品国自产拍蜜月| 欧美激情久久久久久爽电影| 国产精品精品国产色婷婷| 国产 一区 欧美 日韩| 亚洲精品一区av在线观看| 亚洲va在线va天堂va国产| 日韩在线高清观看一区二区三区 | 少妇的逼好多水| 精品99又大又爽又粗少妇毛片 | 久久6这里有精品| 国产欧美日韩精品一区二区| 在线观看一区二区三区| 日韩 亚洲 欧美在线| av天堂中文字幕网| 乱人视频在线观看| av在线老鸭窝| 日本色播在线视频| 午夜精品在线福利| 亚洲国产精品合色在线| 有码 亚洲区| 久久久久久久久久久丰满 | 久久婷婷人人爽人人干人人爱| 99热6这里只有精品| 中文亚洲av片在线观看爽| 国产精品精品国产色婷婷| 十八禁国产超污无遮挡网站| 一级a爱片免费观看的视频| av女优亚洲男人天堂| 黄色日韩在线| 免费看美女性在线毛片视频| 成人综合一区亚洲| 男女做爰动态图高潮gif福利片| 亚洲人与动物交配视频| 国产熟女欧美一区二区| 最近中文字幕高清免费大全6 | 欧美不卡视频在线免费观看| 久久精品综合一区二区三区| 老师上课跳d突然被开到最大视频| 无遮挡黄片免费观看| 熟女人妻精品中文字幕| 国产熟女欧美一区二区| 69av精品久久久久久| 在线观看舔阴道视频| 国产av不卡久久| 天堂影院成人在线观看| 少妇被粗大猛烈的视频| 国产视频内射| 久久热精品热| 国产精品嫩草影院av在线观看 | 中文字幕av成人在线电影| 春色校园在线视频观看| 91狼人影院| 国产精品三级大全| 国产成人a区在线观看| 热99在线观看视频| 干丝袜人妻中文字幕| 大型黄色视频在线免费观看| 毛片女人毛片| 国产乱人伦免费视频| 久久中文看片网| aaaaa片日本免费| 悠悠久久av| 超碰av人人做人人爽久久| 老女人水多毛片| 国产免费av片在线观看野外av| 美女高潮的动态| 精品人妻1区二区| 一进一出抽搐动态| 一区二区三区高清视频在线| 亚洲四区av| 欧美三级亚洲精品|