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

    Revisit submergence of ice blocks in front of ice cover-an experimental study *

    2018-05-14 01:43:20JunWang王軍YifanWu吳一帆JueyiSui
    關(guān)鍵詞:王軍

    Jun Wang (王軍), Yi-fan Wu (吳一帆), Jueyi Sui

    1. College of Civil and Hydraulic Engineering, Hefei University of Technology, Hefei 230009, China

    2. Environmental Sciences and Engineering Program, University of Northern British Columbia, Prince George,BC, Canada

    Introduction

    The river ice is an important hydrologic factor in the temperate and polar environments. An ice cover alters the hydraulics of an open channel by imposing an extra boundary for the flow, changing the flow depth and the flow velocity as compared with those under the ice free conditions[1,2]. Under specific conditions such as with a continuous vast amount of incoming frazil ice from upstream due to very low air temperature and in a highly turbulent flow, the ice jam may be formed, to result in a significant increase of the water level, accompanied with the ice flooding,and the damage of hydraulic structures. Without doubt,the changes of the water level caused by the river ice jams attract interest of many researchers and engineers in the world. Generally, the development of an ice jam might be due to the following processes[3]: (1) when the incoming ice blocks/frazils approach the front edge of the ice cover, depending on the flow conditions, some ice blocks and ice pans might be entrained by the water and submerged and then accumulated under the ice cover, (2) the ice jam can be formed by the collapse of the ice cover to result in a congestion of broken ice blocks in the channel,normally in a bend channel or in front of bridge piers or point bars during the river breakup.

    Ashton[4]pointed out that the critical flow velocity for the entrainment of the ice block in front of the ice cover ranged from 0.60 m/s to 0.70 m/s.

    To study the entrainment of the ice block in front of an ice cover,Sinotin and Grentinm[5] carried out physical model experiments in laboratory. They derived the following formula for determining the critical flow velocity (VC) for entrainment of the ice block in front of an ice cover

    whereLis the length of the ice block.Ashton[4]pointed out that the critical Froude number

    Using the data obtained from experiments,(FrC) for the entrainment of the ice block in front of an ice cover only depends on the ratio of the thickness of the ice block to the water depth under the ice cover(t/h) provided that the frontal ice cover is not submerged under the water. Ashton[4]claimed that the critical Froude number for the entrainment of the ice block is hardly related to the ratio of the thickness of the ice block to the length of the ice block (t/L). He developed the following formula for determining the critical flow Froude number for entrainment of the ice block in front of an ice cover:

    Based on the field measurements of the submergence of ice blocks in front of an ice cover in the Hequ Reach of the Yellow River and experiments in laboratory, Sui et al.[6-9]studied the impacts of flow Froude number on the formation of ice jam and ice accumulation. They found that the critical Froude number for the entrainment of the ice block in front of an ice cover depends not only on the ratio of the thickness of the ice block to the water depth under the ice cover, but also on the ratio of the thickness of the ice block to the length of the ice block. They derived the following formula for determining the critical flow Froude number for the entrainment of the ice block in front of an ice cover

    Based on the continuity equation, Pariset et al.[3]studied the critical condition for the entrainment of the ice block in front of an ice cover. With the ice Froude number (Fri) as the variable to determine whether or not an ice block will be entrained in front of an ice cover, it is shown that

    Wang[10-12]carried out experiments to study the critical condition for the entrainment of the ice block in front of an ice cover. Also with the ice Froude number as the variable, he determines if the ice block will be entrained in front of an ice cover. Based on the data obtained in laboratory, he found that the critical ice Froude number is a function of the dimension of the ice block as follows

    By measuring the pressure difference and the pressure distribution under a floating ice block in laboratory, Ambtman et al.[13]used the method of the moment equilibrium to study the stability of the ice block in front of an ice cover.

    Using the experimental data from Ambtman et al.[13], Lian et al.[14]carried out numerical simulations.Based on Ashton’s formula (Eq. (3)), they derived the following formula for assessing whether or not the natural broken ice block can be entrained in front of an ice cover

    The coefficientKranges from 1.15 to 1.35.They pointed out thatKtakes the minimum value if the frontal cross section of the ice block assumes a quasi-rectangular shape.

    Obviously, it is very important to study the entrainment of the ice block in front of an ice cover,since the submergence of the ice block may result in the formation of the river ice jam. In the present study,based on our experimental data and the data obtained from previous studies, using the method of the moment equilibrium, the critical condition for the entrainment of the ice block in front of an ice cover is further re-visited.

    1. Experiment setup

    The experimental study is conducted in the flume at Hefei University of Technology. As shown in Fig.1,the flume has a dimension of 26.68 m long and 0.40 m wide. There are 22 observation cross sections (CS)along the flume with an equal distance of 1.2 m. The ice-hoper for adding the model ice into the flume is placed between CS-1 and CS-2. The model ice cube used in this experimental study is made of polypropylene with the mass density of 0.918×103kg/m3, which is nearly the same as the mass density of the natural ice of 0.917×103kg/m3. The model ice blocks used in this experimental study have the following dimensions(length×width×height): 0.02 m×0.02 m×0.01 m, 0.03 m×0.03 m×0.01 m, 0.04 m×0.04 m×0.01 m and 0.06 m×0.06 m×0.01 m, respectively, the same sizes (the same dimensions) as those used by Healy and Hicks[15].The flow depths are 0.10 m, 0.15 mand 0.20 m, respectively. The average approaching flow velocities are: 0.05 m/s, 0.07 m/s, 0.09 m/s, 0.11 m/s,0.12 m/s, 0.13 m/s, 0.14 m/s, 0.15 m/s, 0.16 m/s, 0.17 m/s,0.18 m/s, 0.20 m/s and 0.22 m/s, respectively.

    Fig.1 The setup of the flume for experiments in laboratory

    In the downstream channel section from CS-21 to CS-22, one styrofoam plate is placed on the water surface as the model ice cover. The blockage of the model ice cover in downstream will initiate the formation of the ice cover which consists of ice blocks/cubes. This ice cover will progress upstream,depending on the initial hydraulic conditions, such as the flow velocity and the ice discharge.

    Experiments show that the amount of ice blocks entrained under the ice cover can be classified into the following four entrainment levels, “a few ice blocks”,“some ice blocks”, “many ice blocks” and “none ice blocks”. If the flow velocity is low enough or the ice block is large enough, the ice blocks can not be entrained, as shown in Fig.2(a). In this case, the ice jam can not be formed. If only a few ice blocks can be entrained under the ice cover, it means that either the flow velocity is low or the size of the ice block is large. In this case, the ice cover is formed by the juxtapositioning blockage process of ice blocks. Normally, the thickness of the ice cover is the thickness of an individual ice block. Forced by the flowing current,the entrainment process of ice blocks is normally a sliding movement from the frontal edge of the ice cover to the bottom of the ice cover, as shown in Fig.2(b). With the increase of the flow velocity or the decrease of the ice block size, the amount of ice blocks entrained under the ice cover increases to achieve the entrainment level of “some ice blocks”.With further increase of the flow velocity or decrease of the ice block size, the movement of the entrained ice blocks is changed from “slide under the ice cover”to “semi-rolled” under the ice cover, as shown in Fig.2(c). When the flow velocity is high enough or the ice block size is small enough, the amount of ice blocks entrained under the ice cover will be further increased to reach the entrainment level of “many ice blocks”.One can see from the experiments that the majority of incoming ice blocks from upstream will be entrained and forced by the water to “roll” under the ice cover,as shown in Fig.2(d). Only small amount of ice blocks can be stopped by the ice cover and accumulated in front of the ice cover.

    2. Forces acting on ice blocks in front of ice cover

    2.1 Rolled entrainment of ice blocks

    Fig.2 Diving/submerge of ice blocks in front of ice cover

    If the ratio of the ice block thickness to the ice block length (t/L) is more than 0.1 but less than 0.8,namely, 0.1<t/L<0.8, the ice block will be rotated at the lower corner (in contact with the ice cover)during the entrainment process, as shown in Fig.3.

    To assess whether or not an ice block will be rotated to submerge under the ice cover, the forces acting on an individual floating ice block in front of the ice cover are analyzed. When an ice block moves with the speed ofVtoward the frontal edge of the ice cover, all forces acting on the ice block are as shown in Fig.4.

    Conceptually, the following forces are acting on an ice block at the frontal edge of the ice cover: the drag force, the lifting force, the gravity force, the collision forces between the ice block and the frontal edge of the ice cover, and the hydrodynamic pressure force, as shown in Fig.4. The method of the moment equilibrium is used to study the stability of an individual ice block at the frontal edge of the ice cover.Special attention should be paid on the moments produced by the hydrodynamic pressure force, the drag force and the collision forces between the ice block and the ice cover.

    (1) The drag force acting on an individual ice block caused by the water can be determined by the following equation

    Fig.3 The “rolled” entrainment process of ice block in front of ice cover

    Fig.4 Forces acting on an ice block

    whereCDis the drag coefficient,Vis the moving speed of the ice block,Ais the cross sectional area of the ice block which is perpendicular to the flow direction. According to Pang[16], the drag force acting on an individual ice block can be further expressed as

    wheret,Landbare the thickness, the length and the width of an individual ice block, respectively.

    The moment generated by the drag force (MD)on an individual ice block is

    (2) If the mass of an individual ice block (M)and the duration time of the collision (T) between an individual ice block and the frontal edge of the ice cover are known, the collision force between an individual ice block and the ice cover can be determined by the following equation

    From the theorem of the kinetic energy, the total work (W) done by the collision force during the collision process between an individual ice block and the frontal edge of the ice cover is

    wherex0is the moving distance of an individual ice block during the collision. Since the collision force can be further expressed as

    Thus, the moment resulted from the collision force between an individual ice block and the frontal edge of the ice cover is

    (3) Dow Ambtman et al.[17]studied the pressure distribution under an ice cover. The hydrodynamic pressure force under the ice cover (FL)mainly consists of the Venturi pressure forces (FV) and the effective stress at the frontal edge of the ice cover(FE):

    wherehis the water depth, φ is the reduction factor,Uis the flow velocity beneath an individual ice block,Vis the flow velocity in front of the frontal edge of the ice block. In their study, Ambtman and Hicks[17]definedX50such: when the difference between the initial pressure and the Venturi pressure is 50%, then,X50is the distance from the frontal edge of the ice block to the point. According to Song et al.[18], the reduction factor changes with bothLandX50, (1) whenL<0.5X50, then φ=1.0, (2)when 0.5X50<L<3.0X50, then φ=0.5, (3) whenL>3.0X50, then φ=0.2. In the present study, φ=0.5. Thus, the moment generated from the pressure force from the bottom of the ice block (ML) can be calculated as follows:

    Based on experiments carried out in laboratory,Uzun er and Kennedy[19]studied the anti-overturning moment (R) based on experiments carried out in laboratories. They pointed out that the magnitude of the anti-overturning moment (R) depends on the angle between the individual ice block and the frontal edge of the ice cover during the overturning process.Here are formulas for calculating the anti-overturning moment (R)

    (1) As shown in Fig.3(c), the angle between an individual ice block and the frontal edge of the ice cover (θ1) can be calculated as follows

    If 0≤θ≤θ1, the anti-overturning moment can be calculated as follows

    (2) As shown in Fig.3(d), the angle between an individual ice block and the frontal edge of the ice cover (θ2) can be calculated as follows

    If θ1≤θ≤θ2, the anti-overturning moment can be calculated as follows

    (3) If θ2≤θ≤ π /2, the anti-overturning moment can be calculated as follows

    Coutermarsh and Randy Mc Gilvary[20]found that the maximum anti-overturning moment (Rmax) for an individual ice block to be entrained can be expressed as follows

    Thus, to determine whether or not an individual ice block will be entrained through the overturning process, the overturning moment (M0) generated from the forces acting on the ice block should meet the following requirement, namely

    Therefore, the following formula will be used to determine if an individual ice block is entrained through the overturning (or rolled) process

    where

    To calculate the critical speed (VC), Eq. (29) is further simplified as the follows:

    The existence of the collision force at the frontal edge of the ice cover is responsible for the instability of the ice block, and an overturning moment will be generated. As a consequence, ice blocks become much easier to be entrained and submerged under the ice cover. The above derived equation for determining whether or not an ice block will be entrained under the ice cover is based on the assumption that the ice block width (b) (which is shorter than the ice block length(L)) is in contact with the frontal edge of the ice cover. When the longer side of an individual ice block(namely, the ice block length) is in contact with the frontal edge of the ice cover, the critical speed (VC)for the entrainment of the ice block is:

    Thus, the critical ice Froude Number for the entrainment of the ice block in front of the ice cover can be expressed as follows

    In the present study, due to the low flow velocity,before the entrainment of the ice blocks, they are assumed to be located immediately in front of the ice cover. Thus, the moments generated by the collision force between the ice block and the frontal edge of the ice cover as well as the drag force can be ignored. It means that only the moment generated by the Venturi pressure force will be considered for the entrainment of the ice blocks in front of the ice cover. Thus, Eq.(33) can be simplified as follows

    Equation (35) is different from Ashton’s formula(Eq. (3)) due to the difference in evaluating the antioverturning moment. In Eq. (3), the anti-overturning moment (R) is determined when θ=θ1, as follows

    In the present study, the maximum anti-overturning moment (Rmax) is used.

    Based on the physical model study, Dow Ambtman et al.[17]found that the difference of the hydrodynamic pressure under the ice block consists of not only the effect of the Venturi pressure but also the effect of the frontal edge of the ice cover. Then, Eq.(33) can be simplified as follows

    Overally, in all existing studies, either it is assumed that the ice blocks in front of the ice cover are under a stationary condition or under specific conditions before their entrainment. In the present study, however, the critical flow velocity for the entrainment of an individual ice block in front of the ice cover is derived by considering the collision force,the drag force and the hydrodynamics pressure force beneath the ice block.

    2.2 Vertical movement of ice blocks during entrainment

    Experiments show that, when the ratio of the ice block thickness (t) to the ice block length (L) is in the range of eithert/L<0.1 ort/L>0.8, the ice block is normally forced (by current) to move vertically first. When the top surface of an ice block reaches the same level of the bottom surface of the frontal ice cover, the ice block will be forced to roll and accumulate under the ice cover, as shown in Fig.5.

    Through the analysis of the forces acting on an individual ice block, the criteria for the ice block entrainment under the ice cover can be described as follows

    where

    From Eqs. (37), (38), we obtain

    Fig.5 Vertical movement of ice block during entrainment process

    In the above section, the entrainment process of ice blocks in front of the ice cover is described by considering the drag force caused by the flowing water, the hydrodynamic pressure force and the collision force between the ice block and the frontal edge of the ice cover. The ice block can be submerged either by the rolled (or rotated) entrainment or the vertical movement entrainment. Formulas for determining the critical flow velocities for the entrainment of the ice block are derived.

    3. Experimental results

    Clearly, the criteria for the entrainment process of ice blocks in front of the ice cover depend on the dimension of the ice block and the flow velocity. In this experimental study, four types of ice blocks with different dimensions are used to investigate the entrainment process of ice blocks in front of an ice cover.As shown in Fig.6, the length of the ice block is one important factor affecting the flow velocity for the entrainment of the ice block. One can see, the longer the ice block, the higher the flow velocity is needed for the entrainment of the ice block in front of the ice cover.

    Fig.6 Relationship between the flow velocity for the entrainment of ice block and the ice block length

    Fig.7 The dependence of critical flow velocity for the entrainment of ice block on the flow depth

    By mean of physical models, Sinotin and Grentinm[5]developed one formula for determining the flow velocity for the entrainment of ice blocks in front of an ice cover. As shown in Fig.6, the calculated flow velocities for the entrainment of ice blocks using the formula proposed by Sinotin and Grentinm[5]are much less than the flow velocities measured in our the experiments in the laboratory. The calculated flow velocities for the entrainment of ice blocks using the formula proposed by Wang et al.[10]are much closer to the experiment results than those using the formula proposed by Sinotin and Grentinm[5].

    It is shown that the criteria for the entrainment process of ice blocks also depend on the flow depth,as shown in Fig.7. One can see that the critical flow velocity for the entrainment of ice blocks decreases with the increase oft/h.

    Ashton developed a formula for calculating the critical flow velocity for the entrainment of ice blocks under specific conditions. In his derivation of the formula, onlyt/his considered, without considera tionb oft/L. It is found that when the ice block is neither so thick, nor so thin, Ashton’s formula could be used to determine the critical flow velocity for the entrainment of ice blocks, with reasonable results.However, in reality, the critical flow velocity for the entrainment of ice blocks not only depends ont/h,but also ont/L.

    In the present experimental study, to determine the critical flow velocity for the entrainment process of ice blocks in front of an ice cover, the following four types of ice blocks are used with different dimensions (length×width×height): 0.02 m×0.02 m×0.01 m, 0.03 m×0.03 m×0.01 m, 0.04 m×0.04 m×0.01m and 0.06 m×0.06 m×0.01 m.

    According to Uzuner and Kennedy[18], if 0.1<t/L<0.8, the ice blocks will be mainly rolled during the entrainment process. If the ice blocks have the dimension of eithert/L<0.1 (thin ice block) ort/L>0.8 (thick ice block), during the entrainment process, the ice block will be first forced by the current to move vertically, and then, move to the downstream. In the present study, all four types of ice blocks used have dimensions of 0.1<t/L<0.8.Therefore, the ice blocks will be mainly forced to roll during the entrainment process, and Eq. (31) will be used for the calculation.

    Table 1 gives the comparison between the calculated flow velocities for the entrainment of ice blocks in front of an ice cover using different formulas and those of measured velocities in laboratory. One can see from Table 1 that the proposed formula in this paper gives the best results.

    4. Conclusions

    In the present study, using four types of ice blocks with different dimensions, experiments are carried out in laboratory to study the stabilities of ice blocks in front of an ice cover. The critical condition for the entrainment of ice blocks in front of an ice cover is studied by considering the drag force caused by the flowing water, the collision force, and the hydraulic pressure force. Following results are obtained:

    Table 1 Comparisons between calculated critical velocities for the entrainment of ice block using different formulas and the measured velocities

    (1) Using the method of the moment equilibrium,the critical condition for the entrainment of ice blocks in front of an ice cover is re-visited. The critical ice Froude number is obtained.

    (2) To determine whether or not an ice block will be entrained under the ice cover, botht/handt/Lare considered in deriving the formula for determining the critical velocity for the entrainment of ice blocks.The calculated results by using the proposed formula are compared with the measured velocities in laboratory. The calculated results by using the proposed formula are also compared with those using formulas developed by other researchers. Results indicate that the values obtained by using the proposed formula are in the best agreement with the experimental results.

    [1] Fu H., Guo X. L., Yang K. L. et al. Ice accumulation and thickness distribution before inverted siphon [J].Journal of Hydrodynamics, 2017, 29(1): 61-67.

    [2] Wang J., Li S. Y., Chen P. P. et al. Impacts of bridge piers on water level during ice jammed period in bend channelan experimental study [J].Journal of Hydrodynamics.2018, 30(1): 160-168.

    [3] Pariset E., Hausser R., Gagnon A. Formation of ice covers and ice jams in rivers [J].Journal of the Hydraulics Division,1966, 92(6): 1-24.

    [4] Ashton G. D. Froude criterion for ice-block stability [J].Journal of Glaciology, 1974, 13(68): 307-313.

    [5] Sinotin V. I., Grentinm Z. A. A study of diving of ice floes under an obstacle [C].Proceedings of the IAHR International Symposium on Ice and Its Action on Hydraulic Structure, Reykjavik, Iceland, 1970.

    [6] Sui J., Fang D., Zhou Y. et al. Study of critical submergence of ice blocks in front of ice cover [J].Journal of Hydraulic Engineering, 1993, (10): 46-51(in Chinese).

    [7] Sui J., Karney B., Fang D. Variation in water level under ice-jammed condition–Field investigation and experimental study [J].Nordic Hydrology (Now: Hydrology Research), 2005, 36(1): 65-84.

    [8] Sui J., Karney B., Sun Z. et al. Field investigation of frazil jam evolution: A case study [J].Journal of Hydraulic Engineering, ASCE, 2002,128(8): 781-787.

    [9] Sui J., Wang J., Balachandar R. et al. Accumulation of frazil ice along a river bend [J].Canadian Journal of Civil Engineering, 2008, 35(2): 158-169.

    [10] Wang J. Study of stability of ice block in front of ice cover[J].Yellow River,1997, (1): 9-12(in Chinese).

    [11] Wang J., Shi F., Chen P. et al. Impacts of bridge abutments on the critical condition for initiation of ice cover-an experimental study [J].Journal of Hydrology and Hydro-mechanics, 2015, 63(4): 327-333.

    [12] Wang J., Shi F. Y., Chen P. P. et al. Impacts of bridge pier on the stability of ice jam [J],Journal of Hydrodynamics,2015,27(6): 865-871.

    [13] Ambtman K. D., Steffler P. M., Hicks F. E. Analysis of the stability of floating ice blocks [J].Journal of Hydraulic Engineering, ASCE,2011, 137(4): 412-422.

    [14] Lian J. J., Luo F. X., Zhao X. Numerical simulation of ice submergence in water conveyance channel [J].Journal of Hydraulic Engineering,2013, 44(6): 673-679(in Chinese).

    [15] Healy D., Hicks F. E. Experimental study of ice jam thickening under dynamic flow conditions [J].Journal of Cold Regions Engineering, 2007, 21(3): 72-91.

    [16] Pang Q. Experimental study of forces acting on a solid subject in flowing water [D]. Master Thesis, Nanjing,China: Hohai University, 2005(in Chinese).

    [17] Dow Ambtman K. E., Hicks F. E., Steffler P. M. Experimental investigation of the pressure distribution beneath a floating ice block [J].Journal of Hydraulic Engineering,ASCE, 2010, 137(4): 399-411.

    [18] Song X., Zhao X., Liu B. Numerical simulation of distribution characteristics of water pressure at bottom of submerged ice blocks of channel [J].Journal of WaterResources and Water Engineering, 2014, 25(1): 117-121.

    [19] Uzuner M. S., Kennedy J. F. Stability of floating ice blocks [J].Journal of the Hydraulics Division, 1972,98(12): 2117-2133.

    [20] Coutermarsh B. A., Randy McGilvary W. Static analysis of floating ice block stability [J].Journal of Hydraulic Research, 1993, 31(2): 147-160.

    猜你喜歡
    王軍
    HITTING PROBABILITIES AND INTERSECTIONS OF TIME-SPACE ANISOTROPIC RANDOM FIELD
    石榴樹想法妙
    我要好好來欣賞
    好孩子畫報(2020年5期)2020-06-27 14:08:05
    黃陵祭
    不下戰(zhàn)場的士兵——王軍
    活力(2019年19期)2020-01-06 07:34:36
    蜜蜂和油菜花
    可愛的小丫丫
    生態(tài)景觀在城市規(guī)劃中的應(yīng)用探索
    Impact of bridge pier on the stability of ice jam*
    Simulations of ice jam thickness distribution in the transverse direction*
    国产高清videossex| 国产精品久久久人人做人人爽| 午夜精品久久久久久毛片777| cao死你这个sao货| 精品欧美国产一区二区三| 免费在线观看亚洲国产| 久久草成人影院| 亚洲欧美日韩高清专用| 亚洲aⅴ乱码一区二区在线播放 | 99热这里只有精品一区 | 久久久国产成人免费| 日韩高清综合在线| 精品国产乱码久久久久久男人| 在线观看一区二区三区| 欧美zozozo另类| 国产亚洲欧美98| 一区二区三区国产精品乱码| 99国产极品粉嫩在线观看| 日韩三级视频一区二区三区| 18禁黄网站禁片午夜丰满| 国产成人av激情在线播放| 757午夜福利合集在线观看| 国产黄片美女视频| 露出奶头的视频| 婷婷精品国产亚洲av在线| 久久香蕉精品热| www日本在线高清视频| 中文字幕高清在线视频| www国产在线视频色| 毛片女人毛片| 免费电影在线观看免费观看| 性色av乱码一区二区三区2| 免费在线观看成人毛片| 国产午夜精品论理片| 天天躁夜夜躁狠狠躁躁| 91成年电影在线观看| 麻豆国产av国片精品| 黑人欧美特级aaaaaa片| 两个人的视频大全免费| 国产精品九九99| 好男人在线观看高清免费视频| www.熟女人妻精品国产| 熟女电影av网| 男女午夜视频在线观看| 搡老岳熟女国产| 亚洲精品在线美女| 亚洲第一电影网av| 麻豆国产av国片精品| 国产精品亚洲av一区麻豆| 人人妻人人看人人澡| 午夜福利高清视频| 国产亚洲欧美在线一区二区| 久99久视频精品免费| 最近最新中文字幕大全电影3| 两性午夜刺激爽爽歪歪视频在线观看 | 国产人伦9x9x在线观看| 亚洲va日本ⅴa欧美va伊人久久| av天堂在线播放| 视频区欧美日本亚洲| 黄色片一级片一级黄色片| 窝窝影院91人妻| 老司机靠b影院| 亚洲精品美女久久久久99蜜臀| 一边摸一边抽搐一进一小说| 成人18禁高潮啪啪吃奶动态图| 国产乱人伦免费视频| 草草在线视频免费看| cao死你这个sao货| 精品一区二区三区av网在线观看| 一进一出抽搐gif免费好疼| 亚洲精品美女久久av网站| 真人一进一出gif抽搐免费| 午夜福利成人在线免费观看| 国产激情久久老熟女| 欧美一级毛片孕妇| 99热6这里只有精品| 色av中文字幕| 国产亚洲欧美98| 精品久久蜜臀av无| www日本在线高清视频| 午夜a级毛片| 久久精品91无色码中文字幕| 久久午夜综合久久蜜桃| 老汉色∧v一级毛片| 日本精品一区二区三区蜜桃| 99热6这里只有精品| 欧美日韩国产亚洲二区| 在线视频色国产色| 国产成人精品久久二区二区91| 免费看日本二区| 黄片小视频在线播放| 深夜精品福利| 桃色一区二区三区在线观看| 好看av亚洲va欧美ⅴa在| 国产野战对白在线观看| av超薄肉色丝袜交足视频| 日本在线视频免费播放| 国产午夜福利久久久久久| 亚洲精品色激情综合| 这个男人来自地球电影免费观看| 99精品久久久久人妻精品| 一本一本综合久久| 男人的好看免费观看在线视频 | 亚洲av中文字字幕乱码综合| 看黄色毛片网站| 亚洲va日本ⅴa欧美va伊人久久| 手机成人av网站| 成人av一区二区三区在线看| 午夜福利免费观看在线| 欧美另类亚洲清纯唯美| 亚洲精品中文字幕在线视频| 极品教师在线免费播放| 亚洲第一电影网av| 夜夜看夜夜爽夜夜摸| 欧美日韩黄片免| 精品日产1卡2卡| 日韩欧美在线二视频| 国产成人精品久久二区二区免费| 精品不卡国产一区二区三区| 无遮挡黄片免费观看| 国产精品,欧美在线| 极品教师在线免费播放| 国内精品久久久久久久电影| 中文字幕最新亚洲高清| 国产精品野战在线观看| 看片在线看免费视频| 久久伊人香网站| 97超级碰碰碰精品色视频在线观看| 亚洲国产欧美一区二区综合| 最近在线观看免费完整版| 亚洲欧美激情综合另类| 欧美国产日韩亚洲一区| 最近最新中文字幕大全电影3| 黄色毛片三级朝国网站| 成人欧美大片| 精品不卡国产一区二区三区| 国产野战对白在线观看| 91麻豆精品激情在线观看国产| 精品免费久久久久久久清纯| 亚洲欧洲精品一区二区精品久久久| 夜夜看夜夜爽夜夜摸| 久久精品综合一区二区三区| 99热这里只有精品一区 | 少妇熟女aⅴ在线视频| 欧美黑人精品巨大| 好男人在线观看高清免费视频| 精品高清国产在线一区| 一区福利在线观看| 女警被强在线播放| av免费在线观看网站| www.999成人在线观看| 热99re8久久精品国产| 女人爽到高潮嗷嗷叫在线视频| 美女免费视频网站| 国产熟女午夜一区二区三区| 757午夜福利合集在线观看| 天堂√8在线中文| 久久九九热精品免费| 国产激情欧美一区二区| 国产一区在线观看成人免费| bbb黄色大片| 超碰成人久久| 天堂√8在线中文| 成人国语在线视频| 全区人妻精品视频| 久久精品综合一区二区三区| 国产高清videossex| 日韩免费av在线播放| 亚洲欧美一区二区三区黑人| 黄色丝袜av网址大全| 老汉色∧v一级毛片| 精品国产乱码久久久久久男人| 亚洲第一电影网av| 男人舔女人的私密视频| 69av精品久久久久久| 国内精品久久久久精免费| 波多野结衣高清无吗| 国产成人精品无人区| 国产精品一区二区免费欧美| 18禁黄网站禁片免费观看直播| 欧美成人一区二区免费高清观看 | 人妻夜夜爽99麻豆av| 舔av片在线| 亚洲国产中文字幕在线视频| 国产男靠女视频免费网站| 色播亚洲综合网| 国产精品电影一区二区三区| av在线播放免费不卡| 欧美人与性动交α欧美精品济南到| netflix在线观看网站| 18禁黄网站禁片免费观看直播| 国产精华一区二区三区| 亚洲最大成人中文| or卡值多少钱| 两个人看的免费小视频| 亚洲国产精品999在线| 亚洲五月婷婷丁香| 香蕉久久夜色| 长腿黑丝高跟| svipshipincom国产片| 12—13女人毛片做爰片一| 久久精品国产亚洲av高清一级| 男女那种视频在线观看| 在线观看一区二区三区| 日本一区二区免费在线视频| 三级男女做爰猛烈吃奶摸视频| 日韩三级视频一区二区三区| 伦理电影免费视频| 亚洲成人精品中文字幕电影| 精品熟女少妇八av免费久了| 日本成人三级电影网站| 午夜亚洲福利在线播放| 久久国产精品人妻蜜桃| 日韩精品免费视频一区二区三区| 国产成人一区二区三区免费视频网站| 黄色成人免费大全| 欧美性猛交╳xxx乱大交人| 成人18禁高潮啪啪吃奶动态图| ponron亚洲| 日本a在线网址| 国产aⅴ精品一区二区三区波| 男人舔女人的私密视频| 青草久久国产| 99国产精品一区二区三区| 长腿黑丝高跟| 最近最新中文字幕大全电影3| 国产亚洲精品第一综合不卡| 成年人黄色毛片网站| 亚洲专区字幕在线| 久久久久久国产a免费观看| 亚洲avbb在线观看| 2021天堂中文幕一二区在线观| 午夜福利在线观看吧| 在线国产一区二区在线| 国产精品久久久人人做人人爽| 国产精品美女特级片免费视频播放器 | 精品国产乱子伦一区二区三区| 亚洲精品在线观看二区| 成人av在线播放网站| 亚洲午夜精品一区,二区,三区| 男女床上黄色一级片免费看| 人妻久久中文字幕网| 好男人电影高清在线观看| 久久久久国内视频| 又黄又爽又免费观看的视频| 国内毛片毛片毛片毛片毛片| 中文字幕av在线有码专区| 国产精品久久久久久精品电影| bbb黄色大片| 亚洲欧美精品综合一区二区三区| 俄罗斯特黄特色一大片| 国产精品电影一区二区三区| 国产成年人精品一区二区| 人妻久久中文字幕网| 日韩精品青青久久久久久| 在线观看午夜福利视频| 精品福利观看| 成人午夜高清在线视频| 一本大道久久a久久精品| 国产精品一区二区三区四区久久| √禁漫天堂资源中文www| 日本免费一区二区三区高清不卡| 嫁个100分男人电影在线观看| 国产亚洲精品一区二区www| 美女高潮喷水抽搐中文字幕| 色尼玛亚洲综合影院| 男女午夜视频在线观看| 中文字幕精品亚洲无线码一区| 麻豆成人av在线观看| 国产精品精品国产色婷婷| 啦啦啦观看免费观看视频高清| 午夜福利欧美成人| 最新在线观看一区二区三区| 亚洲专区中文字幕在线| 别揉我奶头~嗯~啊~动态视频| 老司机在亚洲福利影院| 51午夜福利影视在线观看| 日日夜夜操网爽| 18禁国产床啪视频网站| 国产久久久一区二区三区| 精品人妻1区二区| 亚洲成av人片免费观看| 夜夜看夜夜爽夜夜摸| 亚洲黑人精品在线| 色老头精品视频在线观看| 久久久国产精品麻豆| 久久久精品大字幕| 久久亚洲精品不卡| 国产精品久久电影中文字幕| 美女免费视频网站| 老司机午夜福利在线观看视频| 天天躁夜夜躁狠狠躁躁| 一级作爱视频免费观看| 日本成人三级电影网站| 国模一区二区三区四区视频 | 久久人人精品亚洲av| 欧美日韩福利视频一区二区| 国产一级毛片七仙女欲春2| 一个人观看的视频www高清免费观看 | 久久精品国产99精品国产亚洲性色| 妹子高潮喷水视频| 欧美乱妇无乱码| 三级男女做爰猛烈吃奶摸视频| 日本a在线网址| 国产亚洲精品av在线| 亚洲精品国产精品久久久不卡| 精品久久久久久久久久久久久| 一级a爱片免费观看的视频| 一进一出抽搐动态| 久久精品人妻少妇| 亚洲人成电影免费在线| 国产高清视频在线播放一区| 亚洲成av人片免费观看| 亚洲中文字幕一区二区三区有码在线看 | 国产成人精品久久二区二区91| 在线国产一区二区在线| 搞女人的毛片| 精品久久久久久久毛片微露脸| 欧美一级a爱片免费观看看 | 波多野结衣巨乳人妻| 长腿黑丝高跟| 国产精品电影一区二区三区| 国产午夜精品久久久久久| 91老司机精品| 久久久久久久久久黄片| 一本精品99久久精品77| 免费av毛片视频| 亚洲av电影不卡..在线观看| 五月伊人婷婷丁香| 亚洲第一欧美日韩一区二区三区| 久久久国产欧美日韩av| 亚洲成a人片在线一区二区| 成人国产一区最新在线观看| 操出白浆在线播放| 国产爱豆传媒在线观看 | 村上凉子中文字幕在线| av片东京热男人的天堂| 亚洲真实伦在线观看| 美女 人体艺术 gogo| 999久久久精品免费观看国产| av超薄肉色丝袜交足视频| 日本熟妇午夜| 亚洲欧美一区二区三区黑人| 欧美激情久久久久久爽电影| 大型av网站在线播放| 亚洲专区字幕在线| 欧美成人性av电影在线观看| 岛国在线观看网站| cao死你这个sao货| 三级男女做爰猛烈吃奶摸视频| 亚洲黑人精品在线| 在线a可以看的网站| 色播亚洲综合网| 色老头精品视频在线观看| 最近最新中文字幕大全免费视频| 国内精品一区二区在线观看| 男女视频在线观看网站免费 | 久久久久九九精品影院| 国产精品一区二区精品视频观看| 人成视频在线观看免费观看| 手机成人av网站| 欧美在线一区亚洲| 国产精品99久久99久久久不卡| 亚洲aⅴ乱码一区二区在线播放 | 1024手机看黄色片| 老司机靠b影院| 国产欧美日韩一区二区精品| 日韩欧美 国产精品| 麻豆成人午夜福利视频| 午夜视频精品福利| 国内精品久久久久久久电影| 久久国产乱子伦精品免费另类| 很黄的视频免费| 欧美+亚洲+日韩+国产| 三级毛片av免费| 日韩欧美国产一区二区入口| 最新在线观看一区二区三区| 欧美在线黄色| 高清在线国产一区| 国产成年人精品一区二区| 国产成人av激情在线播放| 波多野结衣高清作品| 亚洲av成人精品一区久久| 国产高清视频在线观看网站| 精品国内亚洲2022精品成人| 久久久久久国产a免费观看| 看片在线看免费视频| 亚洲无线在线观看| 欧美日本视频| 国产高清有码在线观看视频 | 视频区欧美日本亚洲| 人人妻人人看人人澡| 最近视频中文字幕2019在线8| 黑人操中国人逼视频| 舔av片在线| 一本大道久久a久久精品| 18禁美女被吸乳视频| 国产精品国产高清国产av| 夜夜爽天天搞| 在线观看免费视频日本深夜| 桃红色精品国产亚洲av| 久久久久久人人人人人| 中文字幕精品亚洲无线码一区| 日韩欧美国产一区二区入口| 免费看日本二区| 色噜噜av男人的天堂激情| 欧美另类亚洲清纯唯美| 国产精品,欧美在线| 夜夜躁狠狠躁天天躁| 国产午夜精品久久久久久| 亚洲中文日韩欧美视频| 欧美日韩亚洲综合一区二区三区_| 成人永久免费在线观看视频| 人成视频在线观看免费观看| 日韩中文字幕欧美一区二区| 麻豆成人午夜福利视频| 人人妻人人澡欧美一区二区| 欧美一级毛片孕妇| 一进一出好大好爽视频| 亚洲 欧美一区二区三区| 久久久久久大精品| 国产精品乱码一区二三区的特点| 亚洲黑人精品在线| 久久久久久免费高清国产稀缺| 亚洲欧美日韩高清专用| 日本五十路高清| 欧美丝袜亚洲另类 | 视频区欧美日本亚洲| 韩国av一区二区三区四区| 欧美日韩亚洲综合一区二区三区_| 欧美日韩国产亚洲二区| av超薄肉色丝袜交足视频| 亚洲国产高清在线一区二区三| 后天国语完整版免费观看| 伊人久久大香线蕉亚洲五| 在线看三级毛片| 美女免费视频网站| 一本精品99久久精品77| 亚洲国产欧美人成| 日韩国内少妇激情av| 国产精品久久久久久人妻精品电影| 18禁裸乳无遮挡免费网站照片| 久久香蕉精品热| 久久香蕉激情| 国产av在哪里看| 精品日产1卡2卡| 国产精品久久久久久精品电影| 久久久久久久久免费视频了| av有码第一页| 一本大道久久a久久精品| 老鸭窝网址在线观看| 欧美三级亚洲精品| 女警被强在线播放| 五月伊人婷婷丁香| 九色国产91popny在线| 国产一区二区激情短视频| 亚洲午夜理论影院| 男女视频在线观看网站免费 | 99热6这里只有精品| 国产精品久久久av美女十八| 亚洲天堂国产精品一区在线| 麻豆国产97在线/欧美 | 欧美性猛交黑人性爽| 免费在线观看视频国产中文字幕亚洲| 成人午夜高清在线视频| 精品国产亚洲在线| 国产成人av教育| 18禁黄网站禁片免费观看直播| 丰满人妻熟妇乱又伦精品不卡| 天堂动漫精品| 国产亚洲精品综合一区在线观看 | 国产视频一区二区在线看| 国产av又大| 国产av在哪里看| 国产精品久久电影中文字幕| 老汉色av国产亚洲站长工具| 一本大道久久a久久精品| 中国美女看黄片| 老熟妇仑乱视频hdxx| 两人在一起打扑克的视频| 午夜福利欧美成人| 怎么达到女性高潮| 叶爱在线成人免费视频播放| 日本成人三级电影网站| 狂野欧美白嫩少妇大欣赏| 精品久久久久久,| 正在播放国产对白刺激| 午夜福利免费观看在线| 午夜激情av网站| 久9热在线精品视频| 亚洲自拍偷在线| 中文字幕最新亚洲高清| 国产精品爽爽va在线观看网站| 又黄又爽又免费观看的视频| 真人一进一出gif抽搐免费| 人人妻人人澡欧美一区二区| 黄色丝袜av网址大全| 国产精品日韩av在线免费观看| 不卡av一区二区三区| 日本一二三区视频观看| 人人妻人人澡欧美一区二区| 午夜久久久久精精品| 90打野战视频偷拍视频| 黄色成人免费大全| 亚洲人成电影免费在线| 男人舔女人下体高潮全视频| 亚洲欧美日韩高清在线视频| 国产精品av久久久久免费| 久久国产精品影院| 久久这里只有精品19| 国产不卡一卡二| 99热这里只有是精品50| 听说在线观看完整版免费高清| 97碰自拍视频| 99热这里只有精品一区 | 午夜亚洲福利在线播放| 精品一区二区三区视频在线观看免费| 三级毛片av免费| 欧美黑人精品巨大| 一区二区三区国产精品乱码| 亚洲五月天丁香| 伊人久久大香线蕉亚洲五| 亚洲国产中文字幕在线视频| 国产亚洲欧美在线一区二区| 亚洲av五月六月丁香网| bbb黄色大片| 久久精品aⅴ一区二区三区四区| 脱女人内裤的视频| 久久久久久久精品吃奶| 久久婷婷成人综合色麻豆| 精品久久蜜臀av无| 国产av不卡久久| 色老头精品视频在线观看| 久久香蕉激情| av欧美777| 99国产精品99久久久久| 99久久无色码亚洲精品果冻| 久久久久久人人人人人| 法律面前人人平等表现在哪些方面| 国产午夜精品久久久久久| 曰老女人黄片| 午夜老司机福利片| 久久天堂一区二区三区四区| 好看av亚洲va欧美ⅴa在| 这个男人来自地球电影免费观看| 午夜免费成人在线视频| 日韩欧美国产在线观看| 久久人妻福利社区极品人妻图片| 久久精品91蜜桃| www国产在线视频色| 亚洲av成人av| 老司机在亚洲福利影院| 亚洲人成伊人成综合网2020| 亚洲成人国产一区在线观看| 欧美日本视频| 给我免费播放毛片高清在线观看| 人人妻人人澡欧美一区二区| 亚洲av成人精品一区久久| 搡老妇女老女人老熟妇| 少妇粗大呻吟视频| 熟女少妇亚洲综合色aaa.| 亚洲熟妇熟女久久| 婷婷精品国产亚洲av| 免费在线观看成人毛片| 欧美黄色片欧美黄色片| 在线十欧美十亚洲十日本专区| 久久久国产成人免费| 久久亚洲精品不卡| 91在线观看av| 免费人成视频x8x8入口观看| 亚洲中文字幕日韩| 久久香蕉激情| 午夜福利18| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文字幕日韩| 嫩草影院精品99| 亚洲成a人片在线一区二区| 国产av一区二区精品久久| 亚洲精品av麻豆狂野| 男女下面进入的视频免费午夜| 成年女人毛片免费观看观看9| 久久天堂一区二区三区四区| 欧美最黄视频在线播放免费| 舔av片在线| 亚洲专区国产一区二区| 午夜精品久久久久久毛片777| 日韩免费av在线播放| 久久久久久亚洲精品国产蜜桃av| www日本黄色视频网| 男女那种视频在线观看| 极品教师在线免费播放| 欧美精品啪啪一区二区三区| 日韩大码丰满熟妇| 婷婷精品国产亚洲av在线| 亚洲第一电影网av| 欧美日韩中文字幕国产精品一区二区三区| 亚洲成av人片在线播放无| 久久久久久国产a免费观看| 亚洲va日本ⅴa欧美va伊人久久| 久久久久九九精品影院| 好男人电影高清在线观看| 很黄的视频免费| av在线播放免费不卡| 精品久久久久久久久久免费视频| 变态另类丝袜制服| 熟妇人妻久久中文字幕3abv| av视频在线观看入口| 动漫黄色视频在线观看| 非洲黑人性xxxx精品又粗又长| 18禁黄网站禁片免费观看直播| 一本精品99久久精品77| 久久久久久国产a免费观看| 少妇熟女aⅴ在线视频| 亚洲五月婷婷丁香| 亚洲 国产 在线|