• <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*
    亚洲av片天天在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品自产拍在线观看55亚洲 | 免费不卡黄色视频| 日日摸夜夜添夜夜添小说| 亚洲五月色婷婷综合| 男女无遮挡免费网站观看| av有码第一页| 桃红色精品国产亚洲av| 久久ye,这里只有精品| 在线 av 中文字幕| 日韩成人在线观看一区二区三区| 捣出白浆h1v1| 久9热在线精品视频| 成人18禁在线播放| 国产精品99久久99久久久不卡| av天堂在线播放| 亚洲精品国产精品久久久不卡| 久久性视频一级片| 最新美女视频免费是黄的| 国产精品一区二区精品视频观看| 久热这里只有精品99| cao死你这个sao货| www.精华液| av不卡在线播放| 女人高潮潮喷娇喘18禁视频| 俄罗斯特黄特色一大片| 中文字幕色久视频| 激情在线观看视频在线高清 | 中文字幕另类日韩欧美亚洲嫩草| 久久久久视频综合| 色老头精品视频在线观看| 欧美精品高潮呻吟av久久| 亚洲性夜色夜夜综合| 亚洲国产av新网站| 国产成+人综合+亚洲专区| 乱人伦中国视频| 国产一区二区激情短视频| tocl精华| 久久99热这里只频精品6学生| 成人特级黄色片久久久久久久 | 久久精品国产99精品国产亚洲性色 | 精品少妇久久久久久888优播| 亚洲三区欧美一区| 久久久久精品人妻al黑| 国产极品粉嫩免费观看在线| 亚洲av成人一区二区三| 欧美精品高潮呻吟av久久| 少妇裸体淫交视频免费看高清 | 日韩欧美免费精品| 欧美成人免费av一区二区三区 | 少妇 在线观看| av欧美777| 国产成人精品久久二区二区91| 久久久久视频综合| 亚洲欧美精品综合一区二区三区| av视频免费观看在线观看| 啪啪无遮挡十八禁网站| 久久ye,这里只有精品| 日韩一区二区三区影片| 国产91精品成人一区二区三区 | 中文字幕制服av| 天堂动漫精品| 午夜老司机福利片| 一本久久精品| 久久精品熟女亚洲av麻豆精品| 夜夜爽天天搞| av一本久久久久| 日韩一区二区三区影片| 国产三级黄色录像| 日韩一卡2卡3卡4卡2021年| 亚洲av成人一区二区三| 精品亚洲成a人片在线观看| 一个人免费在线观看的高清视频| 最新在线观看一区二区三区| 一边摸一边抽搐一进一出视频| 亚洲国产看品久久| 十八禁人妻一区二区| av天堂久久9| 中文亚洲av片在线观看爽 | 夫妻午夜视频| 亚洲三区欧美一区| 美女高潮喷水抽搐中文字幕| 精品免费久久久久久久清纯 | 久久午夜亚洲精品久久| √禁漫天堂资源中文www| 欧美日韩av久久| 国产高清激情床上av| 久久99一区二区三区| 国产亚洲欧美精品永久| 变态另类成人亚洲欧美熟女 | 国产99久久九九免费精品| 另类精品久久| 亚洲熟妇熟女久久| av福利片在线| 99国产精品99久久久久| 国产在线观看jvid| 啪啪无遮挡十八禁网站| 黄片播放在线免费| av线在线观看网站| 狠狠狠狠99中文字幕| 精品欧美一区二区三区在线| 老司机亚洲免费影院| 日本欧美视频一区| 欧美精品人与动牲交sv欧美| 国产日韩一区二区三区精品不卡| 老熟妇仑乱视频hdxx| 国产一区二区三区在线臀色熟女 | 欧美日韩视频精品一区| 精品久久久精品久久久| 80岁老熟妇乱子伦牲交| 69精品国产乱码久久久| 亚洲午夜理论影院| 日本欧美视频一区| 日本av免费视频播放| av国产精品久久久久影院| 亚洲专区中文字幕在线| 男女免费视频国产| 9色porny在线观看| 国产在线一区二区三区精| 国产av精品麻豆| 狂野欧美激情性xxxx| 一边摸一边抽搐一进一出视频| 热99久久久久精品小说推荐| 日韩有码中文字幕| 欧美国产精品va在线观看不卡| h视频一区二区三区| 一边摸一边抽搐一进一小说 | 999精品在线视频| 国产免费现黄频在线看| 国产97色在线日韩免费| 老熟妇乱子伦视频在线观看| 99国产精品免费福利视频| 三级毛片av免费| 免费在线观看视频国产中文字幕亚洲| 国产伦理片在线播放av一区| 国产成人一区二区三区免费视频网站| 成年人午夜在线观看视频| 国产一卡二卡三卡精品| 亚洲国产av新网站| 久久国产精品影院| 亚洲中文av在线| 操美女的视频在线观看| 女同久久另类99精品国产91| 精品久久久久久久毛片微露脸| 一本—道久久a久久精品蜜桃钙片| 欧美精品啪啪一区二区三区| 日本黄色日本黄色录像| 久久久久视频综合| 999精品在线视频| 黄色视频不卡| 大香蕉久久成人网| 色综合欧美亚洲国产小说| 真人做人爱边吃奶动态| 中文亚洲av片在线观看爽 | 久9热在线精品视频| 搡老熟女国产l中国老女人| 成年动漫av网址| 国产成人影院久久av| 窝窝影院91人妻| 新久久久久国产一级毛片| 电影成人av| 人人妻人人添人人爽欧美一区卜| 亚洲国产毛片av蜜桃av| 国产欧美日韩一区二区精品| 老司机亚洲免费影院| 久久中文看片网| 18禁黄网站禁片午夜丰满| 欧美性长视频在线观看| 欧美激情 高清一区二区三区| 欧美日韩一级在线毛片| 欧美日韩精品网址| 国产精品.久久久| 亚洲成国产人片在线观看| 黄色成人免费大全| 国产日韩欧美在线精品| 成在线人永久免费视频| 18禁黄网站禁片午夜丰满| av国产精品久久久久影院| 波多野结衣av一区二区av| 男女边摸边吃奶| 窝窝影院91人妻| 久久国产精品男人的天堂亚洲| 亚洲av国产av综合av卡| 国产亚洲精品一区二区www | 色精品久久人妻99蜜桃| 亚洲 欧美一区二区三区| 欧美精品高潮呻吟av久久| 一区二区三区精品91| 色综合婷婷激情| 日韩精品免费视频一区二区三区| 国产国语露脸激情在线看| 亚洲va日本ⅴa欧美va伊人久久| 黄片小视频在线播放| 精品一区二区三区四区五区乱码| 午夜两性在线视频| 亚洲七黄色美女视频| 国产一区二区 视频在线| 涩涩av久久男人的天堂| 亚洲熟女精品中文字幕| 777久久人妻少妇嫩草av网站| 国产精品九九99| 美女扒开内裤让男人捅视频| 欧美日韩黄片免| 极品教师在线免费播放| 人人妻人人添人人爽欧美一区卜| 少妇粗大呻吟视频| 亚洲精品粉嫩美女一区| 国产成人影院久久av| 99久久99久久久精品蜜桃| 午夜福利乱码中文字幕| 国产欧美日韩综合在线一区二区| 亚洲五月色婷婷综合| 免费在线观看完整版高清| 一边摸一边抽搐一进一小说 | 国产精品久久电影中文字幕 | aaaaa片日本免费| av国产精品久久久久影院| 757午夜福利合集在线观看| 欧美黑人精品巨大| 欧美亚洲 丝袜 人妻 在线| 久久久久精品国产欧美久久久| 国产精品熟女久久久久浪| 日韩制服丝袜自拍偷拍| 久久av网站| 麻豆av在线久日| www.精华液| 下体分泌物呈黄色| avwww免费| 99国产精品一区二区蜜桃av | 搡老乐熟女国产| 女同久久另类99精品国产91| 蜜桃在线观看..| 亚洲一区二区三区欧美精品| 久久久久精品国产欧美久久久| 国产精品久久久久久精品古装| 免费观看人在逋| 精品国产乱子伦一区二区三区| 丁香欧美五月| 日韩有码中文字幕| 国产一区二区三区视频了| 满18在线观看网站| 日日摸夜夜添夜夜添小说| 考比视频在线观看| 老司机午夜福利在线观看视频 | 亚洲一区二区三区欧美精品| 精品少妇久久久久久888优播| 18禁观看日本| 一区二区日韩欧美中文字幕| 国产淫语在线视频| 一边摸一边抽搐一进一小说 | xxxhd国产人妻xxx| 99久久人妻综合| 精品免费久久久久久久清纯 | 91精品国产国语对白视频| 日韩成人在线观看一区二区三区| 中亚洲国语对白在线视频| 我的亚洲天堂| 免费黄频网站在线观看国产| 国产福利在线免费观看视频| 一级毛片女人18水好多| 日韩成人在线观看一区二区三区| 国产黄色免费在线视频| 丝袜美足系列| 久久狼人影院| 亚洲五月婷婷丁香| 精品高清国产在线一区| 精品一品国产午夜福利视频| 青草久久国产| 在线观看舔阴道视频| 午夜福利乱码中文字幕| 久久精品成人免费网站| 50天的宝宝边吃奶边哭怎么回事| 国产精品久久久久久精品古装| 叶爱在线成人免费视频播放| svipshipincom国产片| 欧美日韩黄片免| 丝袜美腿诱惑在线| av网站在线播放免费| 国产精品一区二区精品视频观看| 欧美黑人欧美精品刺激| 黄色丝袜av网址大全| 国产在线观看jvid| 国产成人欧美在线观看 | 日本黄色日本黄色录像| 天堂8中文在线网| 在线观看免费视频网站a站| 日韩大片免费观看网站| 乱人伦中国视频| 国产熟女午夜一区二区三区| 久久久精品区二区三区| 最新的欧美精品一区二区| 免费日韩欧美在线观看| 免费女性裸体啪啪无遮挡网站| 国产av国产精品国产| 18禁国产床啪视频网站| 在线播放国产精品三级| 国产免费av片在线观看野外av| 欧美一级毛片孕妇| 无遮挡黄片免费观看| 免费在线观看黄色视频的| 动漫黄色视频在线观看| 中文欧美无线码| 又紧又爽又黄一区二区| 男女边摸边吃奶| 久久久久网色| 亚洲人成77777在线视频| 色视频在线一区二区三区| 80岁老熟妇乱子伦牲交| 成人18禁在线播放| 少妇被粗大的猛进出69影院| 国产亚洲欧美精品永久| 精品视频人人做人人爽| 老司机福利观看| 91成人精品电影| 国产成人av教育| 飞空精品影院首页| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色 视频免费看| 嫁个100分男人电影在线观看| 美女国产高潮福利片在线看| 日本a在线网址| 9191精品国产免费久久| 精品国产一区二区三区久久久樱花| 看免费av毛片| 国产在线免费精品| 中文字幕av电影在线播放| www.999成人在线观看| 亚洲欧洲精品一区二区精品久久久| 久久精品熟女亚洲av麻豆精品| 中文字幕色久视频| 日韩视频在线欧美| 亚洲精品国产一区二区精华液| 在线观看一区二区三区激情| 成年女人毛片免费观看观看9 | 少妇被粗大的猛进出69影院| 国产高清videossex| 精品国产一区二区三区四区第35| 欧美日韩亚洲国产一区二区在线观看 | 国产精品免费视频内射| 久久久久国内视频| 欧美成人午夜精品| 成年人黄色毛片网站| 亚洲欧美激情在线| 久久九九热精品免费| 亚洲欧美激情在线| 香蕉国产在线看| 午夜两性在线视频| 欧美乱妇无乱码| 久久久精品区二区三区| 69av精品久久久久久 | 亚洲精品美女久久久久99蜜臀| 成人三级做爰电影| 久久久久久久大尺度免费视频| 夜夜爽天天搞| 99九九在线精品视频| 99国产精品免费福利视频| 亚洲avbb在线观看| 老熟妇仑乱视频hdxx| 亚洲精华国产精华精| 亚洲欧美一区二区三区久久| 国产一卡二卡三卡精品| 久久久久久亚洲精品国产蜜桃av| 日韩一区二区三区影片| 国产精品av久久久久免费| xxxhd国产人妻xxx| 在线永久观看黄色视频| 日韩大片免费观看网站| 日韩中文字幕视频在线看片| 色婷婷久久久亚洲欧美| 亚洲中文av在线| 国产精品一区二区在线不卡| 自线自在国产av| 91国产中文字幕| 午夜激情av网站| 90打野战视频偷拍视频| 亚洲精品久久午夜乱码| 国产不卡一卡二| 性少妇av在线| 婷婷成人精品国产| 国产精品 国内视频| 一区二区三区国产精品乱码| 久久久精品免费免费高清| 国产精品美女特级片免费视频播放器 | 99国产综合亚洲精品| 高清在线国产一区| www.自偷自拍.com| 久久久久久久精品吃奶| 啦啦啦在线免费观看视频4| 丝袜喷水一区| 免费观看a级毛片全部| 手机成人av网站| 日韩熟女老妇一区二区性免费视频| www.熟女人妻精品国产| 热99久久久久精品小说推荐| 9色porny在线观看| 国产在线免费精品| 一本久久精品| av网站免费在线观看视频| 亚洲色图av天堂| 麻豆av在线久日| 日韩中文字幕欧美一区二区| 精品熟女少妇八av免费久了| 新久久久久国产一级毛片| 中文字幕色久视频| av视频免费观看在线观看| av天堂久久9| 99精品欧美一区二区三区四区| 青青草视频在线视频观看| 精品人妻熟女毛片av久久网站| 啦啦啦在线免费观看视频4| 国产免费现黄频在线看| 超碰成人久久| 亚洲精品自拍成人| 999久久久精品免费观看国产| 亚洲熟女精品中文字幕| 99久久人妻综合| 首页视频小说图片口味搜索| 欧美黑人精品巨大| 成人18禁高潮啪啪吃奶动态图| 十八禁网站免费在线| 久久国产亚洲av麻豆专区| 亚洲av日韩精品久久久久久密| 天天躁日日躁夜夜躁夜夜| 欧美日韩亚洲综合一区二区三区_| 亚洲va日本ⅴa欧美va伊人久久| 啪啪无遮挡十八禁网站| 久久久久久免费高清国产稀缺| 最黄视频免费看| av天堂在线播放| 欧美一级毛片孕妇| 国产一卡二卡三卡精品| 黄色毛片三级朝国网站| 久久中文字幕一级| av免费在线观看网站| 久久精品亚洲av国产电影网| 中文字幕制服av| e午夜精品久久久久久久| 亚洲欧美色中文字幕在线| 法律面前人人平等表现在哪些方面| 久久中文看片网| 日日爽夜夜爽网站| 久久中文字幕一级| 一级a爱视频在线免费观看| 国产色视频综合| 大型黄色视频在线免费观看| 可以免费在线观看a视频的电影网站| 国产主播在线观看一区二区| 黄色怎么调成土黄色| 国产精品.久久久| 九色亚洲精品在线播放| 亚洲av欧美aⅴ国产| 亚洲欧美日韩另类电影网站| 中文字幕人妻熟女乱码| tube8黄色片| 免费av中文字幕在线| 亚洲久久久国产精品| 亚洲国产欧美网| 男女之事视频高清在线观看| 一边摸一边做爽爽视频免费| 视频区图区小说| 国产精品久久久av美女十八| 制服人妻中文乱码| 成人免费观看视频高清| 我要看黄色一级片免费的| 纵有疾风起免费观看全集完整版| 午夜两性在线视频| 国产成人av教育| 搡老岳熟女国产| 成人18禁在线播放| 国产亚洲欧美在线一区二区| 久久精品亚洲熟妇少妇任你| 免费人妻精品一区二区三区视频| 99久久人妻综合| 国产亚洲欧美精品永久| 中文字幕人妻丝袜一区二区| 精品久久久精品久久久| 中国美女看黄片| 在线亚洲精品国产二区图片欧美| 精品久久蜜臀av无| 宅男免费午夜| 王馨瑶露胸无遮挡在线观看| 69av精品久久久久久 | 两人在一起打扑克的视频| 一区福利在线观看| 色视频在线一区二区三区| 久久国产精品大桥未久av| 下体分泌物呈黄色| 波多野结衣一区麻豆| 亚洲美女黄片视频| 久久ye,这里只有精品| 性高湖久久久久久久久免费观看| 在线十欧美十亚洲十日本专区| 超色免费av| 91字幕亚洲| 亚洲第一欧美日韩一区二区三区 | 一级片'在线观看视频| 美女高潮喷水抽搐中文字幕| 女人爽到高潮嗷嗷叫在线视频| 啪啪无遮挡十八禁网站| 在线亚洲精品国产二区图片欧美| 亚洲综合色网址| 国产精品二区激情视频| 美国免费a级毛片| 中亚洲国语对白在线视频| 中文字幕制服av| 婷婷丁香在线五月| 深夜精品福利| 女人爽到高潮嗷嗷叫在线视频| 久久精品91无色码中文字幕| 午夜视频精品福利| 国产亚洲欧美精品永久| 成年人免费黄色播放视频| 久久精品亚洲熟妇少妇任你| tocl精华| 亚洲第一av免费看| 一级毛片女人18水好多| 最近最新中文字幕大全免费视频| 美女视频免费永久观看网站| 青草久久国产| 中文字幕另类日韩欧美亚洲嫩草| 国产在线精品亚洲第一网站| 在线亚洲精品国产二区图片欧美| 久久久久国内视频| 免费在线观看视频国产中文字幕亚洲| 大香蕉久久网| 国产精品 欧美亚洲| 新久久久久国产一级毛片| 桃红色精品国产亚洲av| 曰老女人黄片| 亚洲黑人精品在线| 欧美激情久久久久久爽电影 | 精品国产一区二区三区四区第35| 成人国语在线视频| 亚洲专区中文字幕在线| 国产精品偷伦视频观看了| 亚洲午夜精品一区,二区,三区| 国产精品av久久久久免费| 啦啦啦视频在线资源免费观看| 国产欧美日韩一区二区三| 自拍欧美九色日韩亚洲蝌蚪91| 人人妻,人人澡人人爽秒播| 制服诱惑二区| 久久国产精品影院| 捣出白浆h1v1| 免费不卡黄色视频| 婷婷成人精品国产| 18在线观看网站| 97人妻天天添夜夜摸| 中文字幕人妻丝袜制服| 精品少妇久久久久久888优播| 国产主播在线观看一区二区| 亚洲第一欧美日韩一区二区三区 | 久久久国产精品麻豆| 国产aⅴ精品一区二区三区波| 国产主播在线观看一区二区| 老司机福利观看| 巨乳人妻的诱惑在线观看| 成人国产av品久久久| 午夜福利免费观看在线| 中文字幕最新亚洲高清| av视频免费观看在线观看| 高清av免费在线| 亚洲欧美色中文字幕在线| 老司机午夜十八禁免费视频| 亚洲国产欧美日韩在线播放| 一级片'在线观看视频| 美女高潮到喷水免费观看| 黄色成人免费大全| 亚洲国产av影院在线观看| 宅男免费午夜| av不卡在线播放| 欧美日韩一级在线毛片| 久久久精品免费免费高清| 亚洲国产欧美网| 大香蕉久久成人网| av天堂在线播放| 国产在线精品亚洲第一网站| 久久精品国产亚洲av高清一级| 成人永久免费在线观看视频 | 高清视频免费观看一区二区| 可以免费在线观看a视频的电影网站| 欧美日韩黄片免| 亚洲七黄色美女视频| 夜夜骑夜夜射夜夜干| 国产成人精品在线电影| 午夜福利视频精品| 一本一本久久a久久精品综合妖精| 亚洲成a人片在线一区二区| 麻豆成人av在线观看| 日本a在线网址| xxxhd国产人妻xxx| 一区福利在线观看| 女警被强在线播放| 亚洲全国av大片| 99热国产这里只有精品6| 国产一区二区三区视频了| 国产av又大| 999精品在线视频| 69精品国产乱码久久久| 国产欧美日韩一区二区三区在线| 99香蕉大伊视频| 亚洲国产成人一精品久久久| 在线观看免费午夜福利视频| 美女视频免费永久观看网站| 99精品久久久久人妻精品| 国产深夜福利视频在线观看| 一区二区av电影网| 91成人精品电影| 精品午夜福利视频在线观看一区 | 成人免费观看视频高清| 精品国内亚洲2022精品成人 | 欧美日韩成人在线一区二区| 一个人免费看片子| 国产99久久九九免费精品|