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

    Impacts of ice cover on local scour around semi-circular bridge abutment*

    2014-06-01 12:29:59WUPengHIRSHFIELDFayeSUIJueyi
    關(guān)鍵詞:王軍

    WU Peng, HIRSHFIELD Faye, SUI Jueyi

    Environmental Engineering, University of Northern British Columbia, Prince George, BC, Canada, E-mail: pwu@unbc.ca

    WANG Jun (王軍), CHEN Pang-pang (陳胖胖)

    School of Civil Engineering, Hefei University of Technology, Hefei 230009, China

    Impacts of ice cover on local scour around semi-circular bridge abutment*

    WU Peng, HIRSHFIELD Faye, SUI Jueyi

    Environmental Engineering, University of Northern British Columbia, Prince George, BC, Canada, E-mail: pwu@unbc.ca

    WANG Jun (王軍), CHEN Pang-pang (陳胖胖)

    School of Civil Engineering, Hefei University of Technology, Hefei 230009, China

    (Received June 6, 2013, Revised October 28, 2013)

    The presence of ice cover in winter can significantly change the flow field around bridge abutments, which can also cause a different local scour pattern. To investigate the impacts of ice cover, results from a recent flume experiments were presented. Smooth and rough ice covers were created to investigate the impacts of ice cover roughness on the scour geometry around the semi-circular abutment. Three bed materials were used, with50D sof 0.58 mm, 0.50 mm, 0.47 mm respectively. Scour volume and scour area were calculated. It was found that the maximum scour depth was located 75oinclined to the flume wall. Under rough ice cover, the scour area and scour depth were the largest. An empirical equation on the maximum scour depth was also developed.

    ice cover, local scour, semi-circular abutment, maximum scour depth, scour pattern

    Introduction

    The protrusion of a bridge abutment into the main channel creates disturbance and obstruction to the sediment transport state in the alluvial channel. The flow accelerates and separates at the upstream face of the abutment which creates a down-flow vortex. The direct result is local scour around bridge abutment. The vortex system and down-flow, along with the turbulence, are the main cause of local scour. Essentially, the local scour phenomenon is a dynamic feedback process between the turbulent flow and bed sediment[1].

    Bridge scour has been identified as the most common cause of highway bridge failures and it accounts for about 60% of all bridge collapses in the United States[2]. According to Kandasamy and Melville[3], 6 of 10 bridge failures which occurred in New Zealand during the Cyclone Bola were related to abutment scour.

    Investigations of bridge failure due to local scour around bridge abutments have been an important topic for hydraulic engineers for many years. In 2011, the National Cooperative Highway Research Program (NCHRP) presented two reports on the local scour around bridge foundations[4,5]. As was reviewed in the reports, several commonly used equations were compared to estimate scour depth around bridge foundations. However, none of these equations are applicable for the local scour estimation under ice cover.

    In the northern region of Canada, rivers can be covered by ice during the winter. Ice cover is a threat to the safety of a bridge and can cause serious problems around local ecosystems. Ice cover presents a different set of geomorphological conditions when compared to that of open flow[6]. The characteristics of flow under ice cover impact the bed-load sediment transport, traverse and vertical mixing and mean flow velocity[7]. However, to date, there is still very limited research on the local scour around bridge infrastructures under ice cover[8-14]. Additionally, most of the previous studies were conducted in small-scale laboratory flumes (0.5 m-1.6 m wide). None of these studies were conducted in a large-scale flume, which can better simulate the scour phenomenon around abutments. In the present study, one large-scale flume (2 mwide, 40 m long) was used to study the local scour around a semi-circular abutment. To fill this gap, the present study was conducted to investigate the scour pattern and maximum scour depth around a semi-circular abutment under ice cover.

    Fig.1 Dimensions of abutment, ice cover and rough ice cover used in the experiment

    1. Methodology

    1.1Experimental setup

    A large flume was used at the Quesnel River Research Centre, Likely, BC, Canada. The flume had a dimensions of 40 m long, 2 m wide, 1.3m deep. The slope of the flume bottom was 0.2%. A holding tank with a volume of 90 m3was located in the upstream portion to keep a constant discharge in the experimental zone. At the end of the holding tank, water overflowed from a rectangular weir to the flume. Figure 1(a) shows the semi-circular abutment dimensions in the flume.

    Two sand boxes were created in the flume, with a distance of 10.2 m from each other. To make sure the sand box was deep enough for the local scour development, the sand boxes were both dug to a 0.30 m depth while other parts of the flume were covered by water treated plywood. The velocity range in sand box #1 was 0.16 m/s-0.26 m/s, while in sand box #2, the range was 0.14 m/s-0.21 m/s. The semi-circular abutment model was made from plexiglass. Three nonuniform sediments were used with50D sof 0.58 mm, 0.50 mm, 0.47 mm respectively. In the present study, since ice cover was the main focus, two types of ice cover were created, namely smooth cover and rough cover. Both types of ice cover were attached in the experimental zone as a fixed ice cover on top of the water surface (see Fig.1(b)). The smooth ice cover was the original styrofoam panels while the rough ice cover was made by attaching small Styrofoam cubes to the bottom of the smooth cover (see Fig.1(c)). The cubic pieces had a dimensions of 0.025 m×0.025 m× 0.025 m. The spacing distance between adjacent cubic pieces was 0.035 m.

    1.2Experiment procedure

    The following steps were strictly followed in the experimental study.

    (1) Before each experiment, the abutment model was leveled and fixed in the sand box to make sure the abutment was vertical to the flume bottom. On the outside surface of the abutment, different measuring lines have been drawn for comparing the scour profile at different locations. In all, 13 measuring lines (P-Q) were made along the semi-circular abutment (see Fig.2).

    Fig.2 Measuring points along the semi-circular abutment

    (2) At the beginning of each experiment, the flume was slowly filled to avoid initial scouring. After the water depth was reached, the required velocity was applied in the flume.

    (3) In front of each sand box, a SonTek IQ was installed to measure the approaching flow velocity andwater depth during the experiment. A 10 Hz SonTek ADV was used to measure the velocity in front of the abutment. An adjustable tailgate was installed at the end of flume to change the water depth. Table 1 summarizes the experimental conditions for each flume experiment.

    (4) After 24 h, the flume was drained slowly. The scour depth was manually measured along the outside lines of the semi-circular abutment. In all, 27 experiments have been carried out. Some of the data can be found in Table 1.

    Table 1 Summary of experimental running conditions

    2. Results and discussion

    Fig.3 The local scour around the abutment and the measurement of the scour

    2.1Local scour pattern

    At the end of each experiment, the local scour was manually measured (see Fig.3). The distance from the abutment outside surface to the boundary of scour

    hole was measured. The sediment deposition ridges around the abutment can be seen from Fig.3. The contour of the local scour hole was mapped in the localcoordinate system by Surfer 10, Golden Software. Based upon the contour mapping, both the volume of the scour hole and the scour area were calculated (see Table 1).

    Due to the narrowing effect created by the abutment, we noticed stronger flow turbulence in the experimental zone. When the velocity in the channel was increased, sediment in the toe areas of the abutment was eroded most quickly.

    Three different non-uniform sediments were used here. During the scouring process, relatively fine particles moved first and sediment in the scour hole was gradually coarsened. An armor layer formed on the surface of the scour hole which prevented the scour hole from scouring further. After 24 h, the armor layer covered the whole area around the bridge abutment. After each experiment, sediment samples were collected at different locations around the abutment. We noticed that at the location from G to I, a secondary scour hole was also developed around the abutment.

    Fig.4 The scour profiles around the abutment under different cover conditions (D50=0.50 mm)

    The scour hole pattern and geometry around the semi-circular abutment exhibits features similar to those found by Zhang et al.[1]. Due to the existence of a primary vortex and wake vortex downstream of the abutment, as well as their interaction, the geometry of the scour area in the upstream is significantly different from that in the downstream. The primary vortex is responsible for the scour hole development, which is analogous to the well-known horseshoe vortex in front of bridge piers[15]. At the upstream of the abutment, an obvious scour hole formed while a fine sediment deposition ridge can be seen in the downstream.

    Figure 4 shows the contour map plotted under different flow cover conditions with the sediment50=D0.50 mm. It can be noted that the maximum scour depth around the semi-circular abutment is located at the upstream surface of the abutment. Additionally, with a decrease in sediment size, maximum scour depth increases correspondingly. Our experiments confirm the conclusion drawn by Ettema et al.[16]that the reductions in the scour depth for large sediment were due to large particles impeding the erosion process inside of the scour hole and dissipating some of the flow energy in the erosion area. This is particularly correct for the non-uniform sediment erosion around bridge abutment. Moreover, for the same bed sediment under the same flow conditions, ice cover results in a larger maximum scour depth.

    Figure 4 also indicates that the scour patterns in the vicinity of the abutment under ice covers were similar to those in open channels. Under rough ice cover, due to the effect of the ice cover, the scour depth was increased.

    Fig.5(a) Cross-section along the semi-circular abutment (D50=0.50 mm)

    2.2Scour profiles along the abutment

    To date, no research has been undertaken for plotting scour profiles along the abutment under ice cover. Hence, scour profiles along the abutment border (From P to Q, refer to Fig.2) were plotted to show the elevations changing along the semi-circular abutment. Figure 5(a) shows the variation in scour depth withdifferent bed sediments under the same flow conditions. Figure 5(b) is the cross section of local scour under different conditions. The following points are noted from the figures:

    (1) For all the cross sections, it can be found that the maximum scour depth is located close to E, 75ofrom the flume wall. This is believed to be caused by the primary vortex, which originates at the upstream of the abutment[17]. The primary vortex is forced by the velocity to drift towards the side of the semi-circular abutment. From Dey’s research[17]on scour in clear water, they mentioned that the velocity and scour depth reached their maxima at 90oto the flume wall. However, from the experimental data, the maximum scour depth occurred not at 90obut rather at 75o from the flume wall. This may be due to the non-uniform sediment used in the present research. The locations of maximum scour depth with or without ice cover were all around 75ofrom the flume wall. A greater number of experiments in ice covered channels will improve the estimation of maximum scour location around the semi-circular abutment.

    (2) Moreover, the authors noted that there is a sudden increase in the elevation from G to H, which corresponds to the second scour hole noted from the contour diagram. In literature the second scour hole has been given little attention because of its relatively shallow scour depth compared to the primary scour hole. However, it may explain the migration of the primary vortex flow along the abutment to the downstream, which may in turn explain the downstream wake vortex. Although the reason for the sudden increase in elevation is not clear, one can still note from Fig.5 that the upstream surface has a steeper slope (from P to F). The local slope of the scour hole in the downstream (from F to K) is smaller than that in the upstream. Unfortunately, there is no clear trend showing the changes of upstream slope corresponding to the change in bed sediments50D. It is also noted the same results were found by Zhang et al.[1]on the local scour around the spur dikes in open channels.

    (3) The ice cover has a strong impact on the scour depth around the abutment. As is shown in Fig.5(b), under rough ice covered conditions, the maximum scour depth increases significantly compared to that in open channel and smooth ice cover. From our understanding, the turbulence caused by rough ice cover moves the maximum velocity closer to the bed compared to that by smooth ice cover, which can be attributed to the deeper scour in the vicinity of the abutment. However, as was expected, there are some inaccuracies of the cross section plot due to the limitation of the measurement and profile.

    2.3Scour volume and scour area

    So far, most of the present research conducted on the local scour focuses on the maximum scour depth, while little attention has been paid on the scour volume and scour area. Based on the scour volume and scour area given in Table 1, the scour volume vs. scour area is plotted in Fig.6, which shows the variation of scour volume and scour area around the semicircular abutment in open and ice covered channels.

    Fig.5(b) Cross-section along the semi-circular abutment under smooth and rough cover (D50=0.50 mm)

    Fig.6 Variation of scour volume around bridge abutment

    From Fig.6, the following three relations can be developed:

    For open channel

    Since the ratio of scour volume to scour area is the average scour depth. Based on Eqs.(1)-(3), the average scour depths for rough cover, smooth cover and open channel are 0.041655 m, 0.030915 m, 0.023515 m respectively. One should also note that the above equations are practical only under certain conditions, otherwise, the scour volume would be negative. The average scour depth followed a similar trend to the maximum scour depth. With smooth ice cover, the average scour depth increases by 31.5% compared to that in open channels, for rough ice cover, the increase in average scour depth is 34.7% compared to that under smooth ice cover.

    Fig.7(a) Variation of maximum scour depth with the Froude number under different sediment compositions

    2.4Maximum scour depth

    To the author’s knowledge, there is still no experimental measurements on the maximum scour under ice cover. For non-uniform sediments, Melville[18]included sediment non-uniformity in his formula of estimating the maximum scour depth around the bridge abutment under open flow condition. By using sediment size factordKas a parameter, the following equations were developed

    in whichLis the projected abutment length anddKis the particle size factor. When the ratio of50/L D>25, the value ofdKequals 1 from Melville’s previous research, which is not practical for the present study. Since the abutment length remains constant, the nonuniform sediments are valued by including the Froude number as defined by the following equation

    wheregis the gravitational acceleration,oUis the approaching velocity, andHis the approaching flow depth.

    Fig.7(b) The comparison of maximum scour depth in open channel and ice covered condition (D50=0.50 mm)

    The experimental data from Fig.7(a) indicates that under the same flow conditions, fine sediment composition can result in a deeper maximum scour depth. With the same sediment composition, the maximum scour depth increases with the Froude number. An imposed ice cover results in an increased composite resistance, so under ice covered conditions, the maximum scour is greater than that in open channels (see Fig.7(b)). To gain a better understanding of the impact of sediment grain size on the maximum depth, regression analysis was conducted. The maximum scour depth around the semi-circular abutment can be described by the following variables

    In all, 27 experiments have been conducted to investigate the relationship between average scour depth and approaching flow depth, in which 9 experiments were in open channels, 9 experiments were under smooth ice cover and 9 experiments were under rough ice cover. By using the regression analysis, the following equations were derived from all the 27 experiments (Fig.8):

    For open channel

    As was reported by Sui et al.[19], with an increase in velocity and particle size, the maximum scour depth will increase. In the present research, regarding the local scour around the semi-circular abutment, the rough ice cover causes the largest average scour depth compared to those under smooth ice cover and in open channel. Hence, we compared the maximum scour depth under different flow conditions and with different composition of bed sediments. It is interesting to note the geometric characteristics of the local scour depend mainly on the approaching flow velocity, bed sediment grain size as well as the cover condition.

    Fig.8 Dependence of maximum scour depth on related variables

    With Fig.8 and Eqs.(7)-(9), the impact of sediment distribution is studied. To study the impact of ice cover roughness on the local scour development around the semi-circular abutment, the following dimensional variables under covered flow were considered whereinis the ice cover roughness andbnis the channel bed roughness. According to the Hydraulic Design Handbook[20], in an un-vegetated alluvial channel, the total roughnessbnconsists of two parts. One is grain roughness ()n' which results from the size of the particle and the other is skin roughness ()n' because of the existence of the bed forms. The total roughness can be expressed as However, there is no reliable method of estimatingn', so in the present research, the grain roughness was used as the channel bed roughness in the analysis. For mixtures of bed material with significant portions of coarse-grain sizes, the following equation from Hager[21]was used

    The presence of ice cover alters the mean flow distribution and flow turbulence characteristics. For smooth ice cover, because the styrofoam panel has a relatively smooth concrete-like surface, by Ref.[20], the value of 0.013 was adopted. The roughness of the ice cover was changed by attaching small cubes with dimensions of 0.025 m×0.025 m×0.025 m with a distance 0.035 m from each other. By using the results of discharge measurements through the ice and supporting field data related to the observed characteristics of the underside of the ice cover, Carey[22]calculated Manning’s roughness coefficient, which was between 0.0100-0.0281. From his calculation, a constant roughness of 0.0251 was used for the winter period. Li[23]reviewed several methods to calculate the Manning’s coefficient for ice cover, the following equation can be used depending on the size of the small cubes

    in whichskis the average roughness height of the ice underside. So in the present research, the roughness coefficient was calculated as 0.021, which is also in the range of Carey’s calculation.

    By using the regression analysis, the following empirical equation was given

    The correlation between maximum scour depth and above three variables is worth mentioning because as is indicated in Fig.9 the regression relationship is strong. Meanwhile, it is also confirmed that the hypothesis for calculation of the ice cover roughness is correct. One can also note that, under the condition of same flow and bed material, the maximum scour depth under rough ice cover is deeper than that under smooth ice cover. However, under the same flow and cover conditions, since the index for50Dis negative (–4.85), with the decrease in sediment grain size, the maximum scour depth will increase. With the same bed material and covered condition, the approaching velocity has a positive impact on the maximum scour depth.

    Fig.9 Dependence of maximum scour depth on related variables under ice cover

    To apply this empirical equation in the hydraulic engineering, the authors assume that during the winter the ice cover can be treated as smooth. While in early spring, with the ice breaking up and ice jamming processes, the ice cover can be treated as rough, therefore increasing the flow velocities and increasing the local scour around bridge foundations. In this case, the sediment transport increases and the safety of bridge infrastructures will be threatened. One can note from Eq.(14) that, in the same river, with the increase in ice cover roughness, the maximum scour depth increases. During the ice break-up period in spring time, due to the accumulation of small ice chunks under side, local bridge scour should be monitored.

    Compared to the research of armor layer development in Ref.[19], the approaching water depth had a stronger impact on the maximum scour depth compared to approaching velocity in open channels. However, for ice covered flow, the authors found that approaching velocity has a stronger impact compared to that of the approaching water depth. One drawback regarding the proposed empirical equation is that the values of roughness of only two ice covers were tested.

    3. Conclusion

    Experiments have been conducted in a largescale flume to study the impact of ice cover roughness and non-uniform sediment on the local scour around semi-circular abutments. The location of the maximum scour depth along the abutment is 75ofrom the flume wall. We noticed that the downstream slope in the scour hole is smaller compared to that in the upstream. In this research, the Froude number is also used to investigate the impacts of non-uniform sediment composition on local scour. The scour volume and scour area have been calculated and compared to those in open channel, smooth and rough cover conditions. Under ice cover, the average scour depth is always greater compared to that in open channels. The average scour depth under rough ice cover is 35% greater than that under smooth ice cover. By using dimensional analysis, an empirical equation of the maximum scour depth has been developed. The equation indicates that with an increase in sediment grain size, the maximum scour depth decreases correspondingly. In conclusion, ice cover roughness plays an important role for the maximum scour depth development.

    Acknowledgements

    This research is funded partially by the Research Project Awards from UNBC and also the Natural Science and Engineering Research Council of Canada (NSERC). The authors are grateful for the support provided by the staff and students at QRRC.

    [1] ZHANG H., NAKAGAWA H. and KAWAIKE K. et al. Experimental and simulation of turbulent flow in local scour around a spur dike[J]. International Journal of Sediment Research, 2009, 24(3): 33-45.

    [2] DENG L., CAI C. S. Bridge scour: Prediction, modeling, monitoring, and countermeasures-Review[J]. Practice Periodical on Structural Design and Construction, 2009, 15(2): 125-134.

    [3] KANDASAMY J. K., MELVILLE B. W. Maximum local scour depth at bridge piers and abutments[J]. Journal of Hydraulic Research, 1998, 36(2): 183-197.

    [4] NCHRP. Web-only Document 175, Evaluation of bridge-scour research: Pier scour processes and predictions[R]. 2011, NCHRP Project 24-27(01).

    [5] NCHRP. Web-only Document 181, Evaluation of bridge-scour research: Abutment and contraction scour processes and prediction[R]. 2011, NCHRP Project 24-27(02).

    [6] HICKS F. An overview of river ice problems[J]. Cold regions Science and Technology, 2009, 55(2): 175-185.

    [7] ANDRE R., THANG T. Mean and turbulent flow fields in a simulated ice-cover channel with a gravel bed: Some laboratory observations[J]. Earth Surface Processes and Landforms, 2012, 37(9): 951-956.

    [8] ACKERMANN N. L., SHEN H. T. and OLSSON P. Local scour around circular piers under ice covers[C]. Proceedings of the 16th IAHR International Symposium on Ice, International Association of Hydraulic Engineering Research. Dunedin, New Zealand, 2002.

    [9] ETTEMA R., DALY S. Sediment transport under ice. ERDC/CRREL TR-04-20[R]. Cold regions research and Engineering Laboratory, US Army Corps of Engineers, 2004.

    [10] HAINS D. B. An experimental study of ice effects on scour at bridge piers[D]. Doctoral Thesis, Bethlehem, Pennsylvania, USA: Lehigh University, 2004.

    [11] WANG Jun, SUI Jueyi and KARNEY Bryan W. Incipient motion of non-cohesive sediment under ice cover–An experimental study[J]. Journal of Hydrodynamics, 2008, 20(1): 177-124.

    [12] MUNTEANU A., FRENETTE R. Scouring around acylindrical bridge pier under ice covered flow condition-experimental analysis[R]. R V Anderson Associates Limited and Oxand Report, 2010.

    [13] SUI J., FARUQUE M. A. A. and BALANCHANDAR R. Local scour caused by submerged square jets under model ice cover[J]. Journal of Hydraulic Engineering, ASCE, 2009, 135(4): 316-319.

    [14] SUI J., WANG J. and HE Y. et al. Velocity profile and incipient motion of frazil particles under ice cover[J]. International Journal of Sediment Research, 2010, 25(1): 39-51.

    [15] MELVILLE B. W. Local scour at bridge abutments[J]. Journal of Hydraulic Engineering, ASCE, 1992, 118(4): 615-631.

    [16] ETTEMA R., NATAKO T. and MUSTE M. Estimation of scour depth at bridge abutments[R]. NCHRP 24-20, Iowa City, USA, The University of Iowa, 2010,.

    [17] DEY S., BARBHUIYA A. K. Turbulent flow field in a scour hole at a semicircular abutment[J]. Canadian Journal of Civil Engineering, 2005, 32(1): 213-232.

    [18] MELVILLE B. W. Pier and abutment scour: Integrated approach[J]. Journal of Hydraulic Engineering, ASCE, 1997, 123(2): 125-136.

    [19] SUI J., AFZALIMEHR H. and SAMANI A. K. et al. Clear-water scour around semi-elliptical abutments with armored beds[J]. International Journal of Sediment Research, 2010, 25(3): 233-244.

    [20] MAYS L. W. Hydraulic design handbook[M]. New York, USA: MaGraw-Hill, 1999.

    [21] HAGER W. H. Wastewater hydraulics[M]. Berlin, Germany: Springer-Verlag, 1999.

    [22] CAREY K. Observed configuration and computed roughness of the underside of river ice St Croix River Wisconsin[J]. Geological Survey Professional Paper, 1966, 550(Part 2): B192-B198.

    [23] LI S. S. Estimates of the Manning’s coefficient for ice covered rivers[J]. Water Management, Proceedings of the Institution of Civil Engineers, 2012, 165(WM9): 495-505.

    10.1016/S1001-6058(14)60002-0

    * Biography: WU Peng (1984-), Male, Ph. D. Candidate

    SUI Jueyi, E-mail:

    E-mail: jueyi.sui@unbc.ca

    猜你喜歡
    王軍
    HITTING PROBABILITIES AND INTERSECTIONS OF TIME-SPACE ANISOTROPIC RANDOM FIELD
    石榴樹(shù)想法妙
    我要好好來(lái)欣賞
    黃陵祭
    不下戰(zhàn)場(chǎng)的士兵——王軍
    活力(2019年19期)2020-01-06 07:34:36
    蜜蜂和油菜花
    可愛(ài)的小丫丫
    生態(tài)景觀在城市規(guī)劃中的應(yīng)用探索
    Revisit submergence of ice blocks in front of ice cover-an experimental study *
    Impact of bridge pier on the stability of ice jam*
    精品午夜福利在线看| 日本wwww免费看| 国产免费现黄频在线看| 国产精品久久久久久av不卡| 赤兔流量卡办理| 一区二区三区免费毛片| 黄片无遮挡物在线观看| 丰满迷人的少妇在线观看| 99视频精品全部免费 在线| 国产69精品久久久久777片| 亚洲一级一片aⅴ在线观看| 国产不卡av网站在线观看| 一本一本综合久久| 美女国产视频在线观看| 亚洲久久久国产精品| 亚洲国产精品成人久久小说| 亚洲怡红院男人天堂| 亚洲久久久国产精品| 97精品久久久久久久久久精品| av免费在线看不卡| www.av在线官网国产| 超碰97精品在线观看| 2018国产大陆天天弄谢| 夜夜看夜夜爽夜夜摸| 亚洲国产成人一精品久久久| av一本久久久久| 午夜福利在线观看免费完整高清在| 女性被躁到高潮视频| 亚洲怡红院男人天堂| 久久久久久久久久成人| 美女cb高潮喷水在线观看| 免费播放大片免费观看视频在线观看| 久久精品国产亚洲av涩爱| 在线免费观看不下载黄p国产| 一边亲一边摸免费视频| 99视频精品全部免费 在线| av国产久精品久网站免费入址| 亚洲av国产av综合av卡| 乱码一卡2卡4卡精品| 在线观看免费高清a一片| 亚洲精品,欧美精品| 色5月婷婷丁香| 久久99热6这里只有精品| 多毛熟女@视频| 国产免费一级a男人的天堂| 亚洲图色成人| 久久精品久久精品一区二区三区| 18在线观看网站| 午夜影院在线不卡| 狂野欧美激情性xxxx在线观看| 欧美精品一区二区大全| 国产成人精品一,二区| av电影中文网址| 亚洲综合色惰| 91成人精品电影| 久久久久久久精品精品| 妹子高潮喷水视频| 桃花免费在线播放| 国产成人精品一,二区| 亚洲国产最新在线播放| 精品久久久久久电影网| 考比视频在线观看| 99热全是精品| 狂野欧美激情性xxxx在线观看| 纵有疾风起免费观看全集完整版| 国产日韩欧美亚洲二区| 国产在视频线精品| 国产精品国产三级国产av玫瑰| .国产精品久久| videos熟女内射| 丰满少妇做爰视频| 成人国产麻豆网| 日本av免费视频播放| 美女视频免费永久观看网站| 在线 av 中文字幕| 日本午夜av视频| 国产午夜精品一二区理论片| 黑人猛操日本美女一级片| av.在线天堂| 久久久a久久爽久久v久久| 热99国产精品久久久久久7| 中文精品一卡2卡3卡4更新| 美女视频免费永久观看网站| 天美传媒精品一区二区| 九九久久精品国产亚洲av麻豆| 久久国产精品大桥未久av| 在线亚洲精品国产二区图片欧美 | 99热全是精品| 一级二级三级毛片免费看| 亚洲精品日韩av片在线观看| 免费大片黄手机在线观看| 自线自在国产av| 91精品国产国语对白视频| 久久久久久久久大av| 国产熟女欧美一区二区| 涩涩av久久男人的天堂| 中文字幕精品免费在线观看视频 | 国产成人精品福利久久| 国产在线一区二区三区精| 老司机亚洲免费影院| 在线观看国产h片| 日产精品乱码卡一卡2卡三| 国产一区二区三区综合在线观看 | 国产精品三级大全| a级毛片在线看网站| 高清毛片免费看| av国产久精品久网站免费入址| 国产精品人妻久久久影院| 男女免费视频国产| 尾随美女入室| 亚洲美女搞黄在线观看| 王馨瑶露胸无遮挡在线观看| 日韩中文字幕视频在线看片| 男女边摸边吃奶| 久久亚洲国产成人精品v| 国产欧美日韩一区二区三区在线 | 九九爱精品视频在线观看| 国模一区二区三区四区视频| 国产精品人妻久久久影院| 亚洲av中文av极速乱| 高清视频免费观看一区二区| 亚洲图色成人| 国产精品一区www在线观看| 国产精品女同一区二区软件| 人人妻人人澡人人看| 亚洲国产精品一区二区三区在线| 国产午夜精品一二区理论片| 国产精品蜜桃在线观看| 午夜福利影视在线免费观看| 97在线视频观看| 亚洲国产av影院在线观看| 寂寞人妻少妇视频99o| 18+在线观看网站| 青青草视频在线视频观看| 国产极品天堂在线| 久久久久久久久久成人| 精品少妇内射三级| 成人手机av| 亚洲性久久影院| av线在线观看网站| 亚洲国产精品成人久久小说| 搡老乐熟女国产| 亚洲国产av新网站| 亚洲欧美一区二区三区黑人 | 亚洲欧美日韩另类电影网站| 国产高清有码在线观看视频| 久久久亚洲精品成人影院| 久久久久久久久久久丰满| 免费观看性生交大片5| 一区二区av电影网| 国产免费一级a男人的天堂| 日日摸夜夜添夜夜添av毛片| 久久精品国产鲁丝片午夜精品| 性色avwww在线观看| 国产亚洲午夜精品一区二区久久| 国产熟女午夜一区二区三区 | 午夜91福利影院| 嫩草影院入口| 亚洲精品一区蜜桃| 性色av一级| 人妻系列 视频| 少妇猛男粗大的猛烈进出视频| 自拍欧美九色日韩亚洲蝌蚪91| 高清黄色对白视频在线免费看| 高清视频免费观看一区二区| 亚洲av在线观看美女高潮| 曰老女人黄片| 亚洲精品久久午夜乱码| 女人精品久久久久毛片| av网站免费在线观看视频| 精品亚洲成a人片在线观看| 51国产日韩欧美| 91精品国产国语对白视频| 男女边吃奶边做爰视频| 精品久久久久久久久亚洲| 日韩熟女老妇一区二区性免费视频| 男女无遮挡免费网站观看| 久热久热在线精品观看| 在线观看美女被高潮喷水网站| 亚洲av成人精品一二三区| 久久免费观看电影| 国产成人精品福利久久| 久久久精品免费免费高清| 人妻人人澡人人爽人人| 国产亚洲欧美精品永久| 欧美bdsm另类| 久久免费观看电影| 免费高清在线观看日韩| 亚洲精华国产精华液的使用体验| 国产亚洲精品第一综合不卡 | 免费大片黄手机在线观看| 在线观看www视频免费| xxx大片免费视频| 亚洲av国产av综合av卡| 欧美最新免费一区二区三区| 少妇的逼水好多| 久久精品国产a三级三级三级| 亚洲国产色片| 精品国产国语对白av| 精品久久蜜臀av无| 纵有疾风起免费观看全集完整版| 久久久久国产网址| 久久精品熟女亚洲av麻豆精品| 欧美成人精品欧美一级黄| 精品一区二区三区视频在线| 毛片一级片免费看久久久久| 中文字幕最新亚洲高清| 亚洲国产av新网站| 汤姆久久久久久久影院中文字幕| 最近中文字幕2019免费版| 国产成人午夜福利电影在线观看| 97在线人人人人妻| 成人二区视频| 在线 av 中文字幕| av国产久精品久网站免费入址| 国产又色又爽无遮挡免| 韩国高清视频一区二区三区| 亚洲国产欧美日韩在线播放| 女的被弄到高潮叫床怎么办| 九九爱精品视频在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产黄片视频在线免费观看| 热re99久久精品国产66热6| 男女无遮挡免费网站观看| 日韩熟女老妇一区二区性免费视频| 亚洲欧美精品自产自拍| 哪个播放器可以免费观看大片| 狂野欧美激情性xxxx在线观看| 免费黄色在线免费观看| 国产成人午夜福利电影在线观看| 18在线观看网站| av一本久久久久| 免费看av在线观看网站| 亚洲国产日韩一区二区| 超色免费av| 激情五月婷婷亚洲| av女优亚洲男人天堂| 欧美日韩亚洲高清精品| 欧美 日韩 精品 国产| 国产片特级美女逼逼视频| 婷婷色av中文字幕| 国产成人精品一,二区| 成年美女黄网站色视频大全免费 | 高清毛片免费看| 美女福利国产在线| 亚洲av不卡在线观看| 国产精品一区二区在线观看99| 五月伊人婷婷丁香| av国产久精品久网站免费入址| 老女人水多毛片| 99精国产麻豆久久婷婷| 99热国产这里只有精品6| 99久久综合免费| 国产色婷婷99| 天天躁夜夜躁狠狠久久av| 男女高潮啪啪啪动态图| 人妻系列 视频| 黄色一级大片看看| 女人精品久久久久毛片| 一本久久精品| 国产乱人偷精品视频| 寂寞人妻少妇视频99o| 少妇的逼好多水| 热99国产精品久久久久久7| 国产精品久久久久久精品古装| 一级毛片aaaaaa免费看小| 久久狼人影院| 亚洲经典国产精华液单| a级毛片在线看网站| av福利片在线| 极品人妻少妇av视频| 伦理电影免费视频| 王馨瑶露胸无遮挡在线观看| 少妇人妻精品综合一区二区| 香蕉精品网在线| 又黄又爽又刺激的免费视频.| a级片在线免费高清观看视频| 看十八女毛片水多多多| 国产黄色免费在线视频| 99热网站在线观看| 久久97久久精品| 最近2019中文字幕mv第一页| 女的被弄到高潮叫床怎么办| 又黄又爽又刺激的免费视频.| 伊人久久精品亚洲午夜| 男的添女的下面高潮视频| 国产精品一国产av| 亚洲一级一片aⅴ在线观看| 黄色怎么调成土黄色| 亚洲av免费高清在线观看| 国产精品蜜桃在线观看| 国产成人精品一,二区| 尾随美女入室| 国产永久视频网站| 亚洲一级一片aⅴ在线观看| 成人漫画全彩无遮挡| 国产免费一区二区三区四区乱码| 一边摸一边做爽爽视频免费| 久久久a久久爽久久v久久| 欧美日本中文国产一区发布| 制服丝袜香蕉在线| 欧美日韩在线观看h| 一级a做视频免费观看| 精品少妇内射三级| 性色av一级| 免费黄频网站在线观看国产| 国产伦理片在线播放av一区| 高清黄色对白视频在线免费看| 亚洲怡红院男人天堂| 久久午夜综合久久蜜桃| 午夜福利视频精品| 蜜桃国产av成人99| kizo精华| 三级国产精品欧美在线观看| 如日韩欧美国产精品一区二区三区 | 青春草视频在线免费观看| 不卡视频在线观看欧美| 成年美女黄网站色视频大全免费 | 综合色丁香网| 午夜91福利影院| 亚洲成人一二三区av| 亚洲第一av免费看| 岛国毛片在线播放| 一区二区三区精品91| 一边亲一边摸免费视频| 丰满少妇做爰视频| 国产69精品久久久久777片| 国产淫语在线视频| 在线精品无人区一区二区三| 欧美精品人与动牲交sv欧美| 日韩三级伦理在线观看| 黄色怎么调成土黄色| 国产成人aa在线观看| 亚洲国产毛片av蜜桃av| 五月开心婷婷网| 在线观看免费日韩欧美大片 | 久久精品国产亚洲av天美| 亚洲在久久综合| 两个人的视频大全免费| 免费不卡的大黄色大毛片视频在线观看| 啦啦啦视频在线资源免费观看| 欧美3d第一页| 日韩一区二区视频免费看| 精品国产露脸久久av麻豆| 精品一区二区免费观看| 国精品久久久久久国模美| 99热全是精品| 亚洲第一av免费看| 寂寞人妻少妇视频99o| 一本—道久久a久久精品蜜桃钙片| 久久国产亚洲av麻豆专区| 韩国av在线不卡| 免费少妇av软件| 边亲边吃奶的免费视频| 国产男人的电影天堂91| 久久午夜福利片| 欧美成人精品欧美一级黄| 看免费成人av毛片| 成人黄色视频免费在线看| 麻豆成人av视频| 欧美日韩成人在线一区二区| 国产精品一区二区在线不卡| av女优亚洲男人天堂| 国产乱人偷精品视频| 国产高清不卡午夜福利| 日韩不卡一区二区三区视频在线| 国产av国产精品国产| 亚洲精品成人av观看孕妇| 在线观看免费视频网站a站| 久久久久久久久久久久大奶| 精品久久久久久久久av| 国产亚洲精品久久久com| 久久精品国产a三级三级三级| 夫妻性生交免费视频一级片| 少妇人妻久久综合中文| 一二三四中文在线观看免费高清| 妹子高潮喷水视频| 亚洲精品美女久久av网站| 3wmmmm亚洲av在线观看| 免费观看性生交大片5| 丰满饥渴人妻一区二区三| 高清视频免费观看一区二区| 亚洲色图 男人天堂 中文字幕 | 亚洲欧美一区二区三区黑人 | 性色av一级| 日韩在线高清观看一区二区三区| 久久99精品国语久久久| 久久午夜综合久久蜜桃| 少妇人妻精品综合一区二区| h视频一区二区三区| 免费高清在线观看视频在线观看| 国产精品一区二区在线不卡| 韩国av在线不卡| 欧美变态另类bdsm刘玥| 亚洲精品,欧美精品| 黄色怎么调成土黄色| 麻豆精品久久久久久蜜桃| 中国美白少妇内射xxxbb| 制服诱惑二区| 亚洲欧美日韩卡通动漫| 国产精品久久久久久精品古装| 日本黄色日本黄色录像| 日韩中字成人| 嫩草影院入口| 蜜桃久久精品国产亚洲av| 国产精品一区二区在线观看99| av免费在线看不卡| 国产在视频线精品| 十八禁高潮呻吟视频| 国产av精品麻豆| 日本猛色少妇xxxxx猛交久久| 欧美精品亚洲一区二区| 青春草视频在线免费观看| 日韩av不卡免费在线播放| 日韩强制内射视频| 丰满饥渴人妻一区二区三| 91精品一卡2卡3卡4卡| 欧美亚洲 丝袜 人妻 在线| av国产精品久久久久影院| 熟妇人妻不卡中文字幕| 女人久久www免费人成看片| 亚洲欧洲国产日韩| 日本欧美国产在线视频| 精品人妻在线不人妻| 天天躁夜夜躁狠狠久久av| 免费大片18禁| 内地一区二区视频在线| 久久久久久久久久人人人人人人| 亚洲国产av影院在线观看| 老司机亚洲免费影院| 99国产精品免费福利视频| 亚州av有码| 国产精品99久久99久久久不卡 | 男人添女人高潮全过程视频| 母亲3免费完整高清在线观看 | 欧美 亚洲 国产 日韩一| 国产色婷婷99| 日本与韩国留学比较| 91精品伊人久久大香线蕉| 极品少妇高潮喷水抽搐| 你懂的网址亚洲精品在线观看| 国产亚洲午夜精品一区二区久久| 丰满饥渴人妻一区二区三| 不卡视频在线观看欧美| 搡女人真爽免费视频火全软件| 国产亚洲精品久久久com| 赤兔流量卡办理| av播播在线观看一区| 好男人视频免费观看在线| 国产淫语在线视频| 欧美日韩成人在线一区二区| 国产成人a∨麻豆精品| 国产精品一区二区在线不卡| 国产av一区二区精品久久| 色哟哟·www| 免费人成在线观看视频色| 国产一区亚洲一区在线观看| 亚洲人成网站在线播| 只有这里有精品99| 欧美日韩av久久| 精品国产露脸久久av麻豆| 午夜福利视频精品| 欧美日本中文国产一区发布| 亚洲情色 制服丝袜| 少妇被粗大的猛进出69影院 | 欧美日韩在线观看h| 97精品久久久久久久久久精品| 亚洲av欧美aⅴ国产| 精品久久久精品久久久| 日韩熟女老妇一区二区性免费视频| 国产成人免费无遮挡视频| 欧美人与善性xxx| 美女大奶头黄色视频| 免费黄色在线免费观看| 国产黄色免费在线视频| 国产亚洲一区二区精品| 午夜av观看不卡| 中文字幕制服av| 国产一区二区在线观看av| 乱码一卡2卡4卡精品| 国产欧美日韩综合在线一区二区| 热re99久久精品国产66热6| 久久99蜜桃精品久久| 日本黄色片子视频| 亚洲精品日韩av片在线观看| 丝袜在线中文字幕| 成人手机av| 成人二区视频| 99热国产这里只有精品6| 亚洲图色成人| 国产伦精品一区二区三区视频9| 国产一级毛片在线| 国产精品秋霞免费鲁丝片| 丰满饥渴人妻一区二区三| 国产精品无大码| 午夜精品国产一区二区电影| 哪个播放器可以免费观看大片| 嘟嘟电影网在线观看| 欧美日韩在线观看h| 久久久久久久大尺度免费视频| 精品国产一区二区久久| 精品亚洲成a人片在线观看| 亚洲一级一片aⅴ在线观看| tube8黄色片| 国产欧美另类精品又又久久亚洲欧美| 免费少妇av软件| 免费av中文字幕在线| 国产毛片在线视频| 啦啦啦视频在线资源免费观看| 亚洲成色77777| 国产精品久久久久成人av| 交换朋友夫妻互换小说| 麻豆精品久久久久久蜜桃| 22中文网久久字幕| 色婷婷久久久亚洲欧美| 精品人妻偷拍中文字幕| 性色av一级| 亚洲欧美成人综合另类久久久| 黑人巨大精品欧美一区二区蜜桃 | 中国三级夫妇交换| 亚洲精品久久成人aⅴ小说 | 女的被弄到高潮叫床怎么办| 极品人妻少妇av视频| av在线app专区| 国产成人免费观看mmmm| 亚洲成人一二三区av| 亚洲精品久久久久久婷婷小说| 下体分泌物呈黄色| 三级国产精品欧美在线观看| 观看美女的网站| 国产又色又爽无遮挡免| 亚洲一级一片aⅴ在线观看| 日本黄色片子视频| 一级,二级,三级黄色视频| 亚洲成人手机| 国产精品不卡视频一区二区| 色哟哟·www| 内地一区二区视频在线| 日本av免费视频播放| 欧美日韩av久久| 色94色欧美一区二区| 久久97久久精品| 成人午夜精彩视频在线观看| 日本vs欧美在线观看视频| 狠狠精品人妻久久久久久综合| 美女大奶头黄色视频| av不卡在线播放| videosex国产| 一级二级三级毛片免费看| 少妇被粗大的猛进出69影院 | 久久av网站| 免费久久久久久久精品成人欧美视频 | 欧美精品高潮呻吟av久久| 国产精品一国产av| 少妇人妻精品综合一区二区| 熟女电影av网| 精品久久久精品久久久| 一个人看视频在线观看www免费| 日韩av免费高清视频| 免费少妇av软件| 精品国产一区二区久久| 中文字幕免费在线视频6| 最新的欧美精品一区二区| 99久国产av精品国产电影| 成人亚洲欧美一区二区av| 久久综合国产亚洲精品| 少妇被粗大的猛进出69影院 | 欧美激情 高清一区二区三区| 2018国产大陆天天弄谢| 精品人妻一区二区三区麻豆| 国产一区二区在线观看日韩| 婷婷成人精品国产| 久久久亚洲精品成人影院| 亚洲欧美一区二区三区国产| 九九在线视频观看精品| 天堂8中文在线网| 久久亚洲国产成人精品v| 国产成人精品婷婷| 欧美97在线视频| 97超视频在线观看视频| 乱码一卡2卡4卡精品| 91aial.com中文字幕在线观看| a 毛片基地| 日韩av不卡免费在线播放| 欧美 亚洲 国产 日韩一| 国产男女超爽视频在线观看| 王馨瑶露胸无遮挡在线观看| 国产探花极品一区二区| 久久久久久久精品精品| 蜜桃在线观看..| 色哟哟·www| 国产av一区二区精品久久| 国产成人精品婷婷| 在线观看免费视频网站a站| 中国三级夫妇交换| 在线观看一区二区三区激情| 日日撸夜夜添| 国产在线一区二区三区精| 久久久久久伊人网av| 午夜精品国产一区二区电影| 久久久久精品久久久久真实原创| 韩国高清视频一区二区三区| 婷婷色麻豆天堂久久| 秋霞在线观看毛片| 久久热精品热| 99热6这里只有精品| 精品人妻偷拍中文字幕| 亚洲欧美日韩另类电影网站| 天天躁夜夜躁狠狠久久av| 美女国产高潮福利片在线看| 国产成人精品无人区| 亚洲四区av| 交换朋友夫妻互换小说| 五月天丁香电影| 一级,二级,三级黄色视频|