• <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*
    无遮挡黄片免费观看| 精品久久久久久,| 国产大屁股一区二区在线视频| 一本一本综合久久| 超碰av人人做人人爽久久| 男女下面进入的视频免费午夜| 欧美不卡视频在线免费观看| 日本在线视频免费播放| 午夜福利视频1000在线观看| 女生性感内裤真人,穿戴方法视频| 12—13女人毛片做爰片一| 午夜日韩欧美国产| 欧美另类亚洲清纯唯美| 在现免费观看毛片| 观看美女的网站| 他把我摸到了高潮在线观看| 亚洲 国产 在线| 日日摸夜夜添夜夜添av毛片 | 日日夜夜操网爽| 国产日本99.免费观看| 毛片一级片免费看久久久久 | 亚洲精品影视一区二区三区av| 免费看美女性在线毛片视频| 超碰av人人做人人爽久久| 一边摸一边抽搐一进一小说| 欧美精品国产亚洲| 久久久久久久久大av| 久久这里只有精品中国| 亚洲美女视频黄频| 欧美日韩精品成人综合77777| 成人国产一区最新在线观看| 在线a可以看的网站| 亚洲专区国产一区二区| 两个人的视频大全免费| av天堂中文字幕网| 日韩在线高清观看一区二区三区 | 日本一本二区三区精品| 老司机福利观看| 日本撒尿小便嘘嘘汇集6| 久9热在线精品视频| 99国产精品一区二区蜜桃av| 噜噜噜噜噜久久久久久91| 97超视频在线观看视频| 婷婷色综合大香蕉| 成人国产一区最新在线观看| 国产极品精品免费视频能看的| 欧美精品啪啪一区二区三区| 赤兔流量卡办理| 可以在线观看毛片的网站| 国产精品永久免费网站| 啦啦啦啦在线视频资源| 在线免费观看不下载黄p国产 | 亚洲va日本ⅴa欧美va伊人久久| 欧美最黄视频在线播放免费| 成人二区视频| 久久6这里有精品| 精品免费久久久久久久清纯| 亚洲狠狠婷婷综合久久图片| 日韩中字成人| 国产激情偷乱视频一区二区| 亚洲男人的天堂狠狠| 成人午夜高清在线视频| 日韩av在线大香蕉| 色吧在线观看| 在线a可以看的网站| 日韩大尺度精品在线看网址| 成年女人毛片免费观看观看9| 高清毛片免费观看视频网站| 免费人成视频x8x8入口观看| 亚洲最大成人手机在线| 亚洲人成伊人成综合网2020| 国产 一区 欧美 日韩| 老司机深夜福利视频在线观看| 一区二区三区高清视频在线| 午夜福利欧美成人| 亚洲最大成人av| 亚洲欧美清纯卡通| 又爽又黄无遮挡网站| 日本黄色片子视频| 国产久久久一区二区三区| 日韩一区二区视频免费看| 内射极品少妇av片p| 精品一区二区三区av网在线观看| 哪里可以看免费的av片| 色噜噜av男人的天堂激情| 美女免费视频网站| 少妇丰满av| 久久精品国产亚洲av香蕉五月| av在线天堂中文字幕| 亚洲精品国产成人久久av| 最近中文字幕高清免费大全6 | 88av欧美| 深夜精品福利| 波多野结衣高清作品| 国产高清视频在线观看网站| 国产精品一区二区三区四区久久| 久久精品国产清高在天天线| 国产精品久久久久久精品电影| 精品国产三级普通话版| 亚洲精品日韩av片在线观看| 日本与韩国留学比较| 亚洲精品在线观看二区| 久久热精品热| 嫩草影视91久久| 日韩人妻高清精品专区| 国产av在哪里看| 欧美日韩精品成人综合77777| 久久精品国产亚洲av香蕉五月| 日本精品一区二区三区蜜桃| 国产高清视频在线观看网站| 亚洲熟妇熟女久久| a级毛片免费高清观看在线播放| 国产精品三级大全| a级毛片免费高清观看在线播放| 有码 亚洲区| 亚洲av日韩精品久久久久久密| av在线天堂中文字幕| 最近中文字幕高清免费大全6 | 男人狂女人下面高潮的视频| 亚洲久久久久久中文字幕| 精品午夜福利视频在线观看一区| 久久久久精品国产欧美久久久| 国产aⅴ精品一区二区三区波| 国产在线精品亚洲第一网站| 国产欧美日韩精品一区二区| 看免费成人av毛片| 亚洲成人免费电影在线观看| 免费在线观看影片大全网站| 欧美一区二区精品小视频在线| 白带黄色成豆腐渣| av福利片在线观看| 亚洲电影在线观看av| 在线天堂最新版资源| 三级国产精品欧美在线观看| 国产男靠女视频免费网站| 亚洲精品日韩av片在线观看| 美女大奶头视频| 亚洲精品日韩av片在线观看| 精品国产三级普通话版| 欧美日韩精品成人综合77777| 最新中文字幕久久久久| АⅤ资源中文在线天堂| 欧美性猛交黑人性爽| 有码 亚洲区| 日本免费a在线| 亚洲黑人精品在线| 精品久久久久久久末码| 国产精品久久久久久av不卡| 欧美激情在线99| av视频在线观看入口| 日本免费a在线| 美女xxoo啪啪120秒动态图| 亚洲狠狠婷婷综合久久图片| 亚洲国产欧美人成| 长腿黑丝高跟| 亚洲熟妇中文字幕五十中出| 日本精品一区二区三区蜜桃| 88av欧美| 国产高清视频在线观看网站| 亚洲18禁久久av| 天堂动漫精品| 国产一区二区激情短视频| 综合色av麻豆| 午夜精品在线福利| 国产精华一区二区三区| 欧美激情在线99| 国模一区二区三区四区视频| 亚洲七黄色美女视频| 偷拍熟女少妇极品色| 毛片女人毛片| 国产精品美女特级片免费视频播放器| 干丝袜人妻中文字幕| 欧美日韩亚洲国产一区二区在线观看| 亚洲人成网站高清观看| 人人妻人人看人人澡| 欧美精品啪啪一区二区三区| 色哟哟哟哟哟哟| 亚洲美女黄片视频| 国产不卡一卡二| 国产一区二区在线av高清观看| 国产黄a三级三级三级人| 日本a在线网址| 桃红色精品国产亚洲av| av女优亚洲男人天堂| 99热这里只有是精品在线观看| 亚洲久久久久久中文字幕| 免费观看精品视频网站| 老司机深夜福利视频在线观看| 女同久久另类99精品国产91| 99在线人妻在线中文字幕| 国产成年人精品一区二区| 综合色av麻豆| 国产成人影院久久av| 在线免费观看不下载黄p国产 | 久久久久久九九精品二区国产| 色在线成人网| 久久人人爽人人爽人人片va| 啪啪无遮挡十八禁网站| 国产亚洲欧美98| 高清在线国产一区| 亚洲午夜理论影院| 亚洲专区国产一区二区| 日本五十路高清| 99热网站在线观看| 变态另类丝袜制服| 69人妻影院| 人妻少妇偷人精品九色| 中文字幕免费在线视频6| 啦啦啦韩国在线观看视频| 亚洲性久久影院| 日本成人三级电影网站| 亚洲av中文av极速乱 | 久99久视频精品免费| 精品99又大又爽又粗少妇毛片 | 免费av毛片视频| 女人十人毛片免费观看3o分钟| 热99re8久久精品国产| 老师上课跳d突然被开到最大视频| 尾随美女入室| 三级国产精品欧美在线观看| 国产精品av视频在线免费观看| 一级毛片久久久久久久久女| 亚洲国产色片| 国产久久久一区二区三区| 欧美日本视频| 全区人妻精品视频| 搡老岳熟女国产| 亚洲性夜色夜夜综合| 国产淫片久久久久久久久| 性欧美人与动物交配| 少妇的逼水好多| 99久久久亚洲精品蜜臀av| 亚洲美女视频黄频| 精品久久国产蜜桃| 久久九九热精品免费| 色av中文字幕| 欧美激情国产日韩精品一区| 成人av一区二区三区在线看| 少妇丰满av| 嫩草影院新地址| 美女高潮的动态| 国产欧美日韩精品一区二区| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久精品电影| 国产精品女同一区二区软件 | 日本黄大片高清| 高清日韩中文字幕在线| 一区二区三区四区激情视频 | 亚洲精品日韩av片在线观看| 我要看日韩黄色一级片| 18禁黄网站禁片午夜丰满| 尤物成人国产欧美一区二区三区| 一个人看视频在线观看www免费| 男人舔奶头视频| 欧美精品国产亚洲| 久久久久国内视频| 日本在线视频免费播放| 亚洲国产精品合色在线| 欧美精品啪啪一区二区三区| 久久欧美精品欧美久久欧美| 欧美精品国产亚洲| 日韩高清综合在线| 99久久成人亚洲精品观看| 网址你懂的国产日韩在线| 欧美bdsm另类| 深爱激情五月婷婷| x7x7x7水蜜桃| 真人做人爱边吃奶动态| 永久网站在线| 国产一区二区在线观看日韩| 深爱激情五月婷婷| 国产麻豆成人av免费视频| 国产大屁股一区二区在线视频| 国产精品亚洲一级av第二区| 国内揄拍国产精品人妻在线| 国产精品自产拍在线观看55亚洲| 午夜精品在线福利| 亚洲美女搞黄在线观看 | 欧美日韩国产亚洲二区| 亚洲av中文av极速乱 | 久久久久久久亚洲中文字幕| 在线观看免费视频日本深夜| 丰满的人妻完整版| 99riav亚洲国产免费| 欧美极品一区二区三区四区| 国产在线精品亚洲第一网站| 国产视频一区二区在线看| 久久久久久久久大av| 一个人观看的视频www高清免费观看| 精品日产1卡2卡| 欧美日本视频| 国产精品久久久久久亚洲av鲁大| 精品久久久久久久久av| 久久久久国内视频| 国产精品永久免费网站| 一个人看的www免费观看视频| 亚洲美女视频黄频| 亚洲五月天丁香| 九九爱精品视频在线观看| av福利片在线观看| 丰满的人妻完整版| 俺也久久电影网| 男女那种视频在线观看| 伦精品一区二区三区| 亚洲乱码一区二区免费版| 偷拍熟女少妇极品色| 深爱激情五月婷婷| 欧美xxxx性猛交bbbb| 色5月婷婷丁香| 国产一区二区三区av在线 | 欧美激情久久久久久爽电影| 深夜a级毛片| 啪啪无遮挡十八禁网站| 日韩欧美国产在线观看| 日本免费一区二区三区高清不卡| 国产精品精品国产色婷婷| 国内精品久久久久久久电影| 日韩在线高清观看一区二区三区 | 日日撸夜夜添| 免费搜索国产男女视频| 久久久久久久久久久丰满 | 久久久久久久久大av| 99riav亚洲国产免费| 男女那种视频在线观看| 国产精品一区二区性色av| 国产精品三级大全| 女人被狂操c到高潮| 欧美日本视频| 国产真实伦视频高清在线观看 | 校园春色视频在线观看| 老熟妇乱子伦视频在线观看| 亚洲成人精品中文字幕电影| 听说在线观看完整版免费高清| 日本黄色片子视频| 国产午夜精品论理片| 久久国产乱子免费精品| 狠狠狠狠99中文字幕| 精品日产1卡2卡| .国产精品久久| 亚洲熟妇熟女久久| 国产美女午夜福利| a级一级毛片免费在线观看| 中文字幕av成人在线电影| 九色成人免费人妻av| 可以在线观看的亚洲视频| 真人做人爱边吃奶动态| 精品人妻1区二区| 国产精品女同一区二区软件 | 身体一侧抽搐| 九色成人免费人妻av| 美女 人体艺术 gogo| 中国美白少妇内射xxxbb| 日韩欧美精品v在线| av视频在线观看入口| 亚洲av不卡在线观看| 成年版毛片免费区| 久久这里只有精品中国| 午夜老司机福利剧场| 国产综合懂色| 我的老师免费观看完整版| 国产精品精品国产色婷婷| 精品无人区乱码1区二区| 国产精品人妻久久久久久| 免费看a级黄色片| 两个人视频免费观看高清| 日本三级黄在线观看| 全区人妻精品视频| 99在线人妻在线中文字幕| 99热这里只有是精品在线观看| 国产探花极品一区二区| 日韩欧美一区二区三区在线观看| 真人一进一出gif抽搐免费| 人人妻人人看人人澡| 18禁黄网站禁片免费观看直播| 天美传媒精品一区二区| 国内久久婷婷六月综合欲色啪| 小说图片视频综合网站| 国产国拍精品亚洲av在线观看| 国产日本99.免费观看| 亚洲国产日韩欧美精品在线观看| 麻豆成人av在线观看| 最新中文字幕久久久久| 男插女下体视频免费在线播放| 看黄色毛片网站| 日本与韩国留学比较| 少妇猛男粗大的猛烈进出视频 | 极品教师在线视频| 99久久九九国产精品国产免费| 精品乱码久久久久久99久播| 免费观看的影片在线观看| 啦啦啦韩国在线观看视频| 男人舔奶头视频| 99国产精品一区二区蜜桃av| 一本精品99久久精品77| 久久热精品热| 在线观看美女被高潮喷水网站| 精品久久久噜噜| 亚洲精品在线观看二区| 男插女下体视频免费在线播放| 九色国产91popny在线| 一进一出抽搐gif免费好疼| 一级毛片久久久久久久久女| 国产色爽女视频免费观看| 97碰自拍视频| 超碰av人人做人人爽久久| 99久久成人亚洲精品观看| 午夜福利高清视频| 国产精品女同一区二区软件 | 22中文网久久字幕| 国产精品爽爽va在线观看网站| 蜜桃亚洲精品一区二区三区| 一区二区三区激情视频| 床上黄色一级片| 国产亚洲欧美98| 白带黄色成豆腐渣| 88av欧美| 久久99热6这里只有精品| 亚洲人成网站在线播放欧美日韩| 亚洲国产欧洲综合997久久,| 韩国av在线不卡| av.在线天堂| 亚洲av一区综合| 久久精品国产鲁丝片午夜精品 | 亚洲黑人精品在线| 免费在线观看日本一区| 精品人妻1区二区| 欧美成人a在线观看| 窝窝影院91人妻| 男女下面进入的视频免费午夜| 真人一进一出gif抽搐免费| 国产精品电影一区二区三区| 无遮挡黄片免费观看| 欧美成人一区二区免费高清观看| 97热精品久久久久久| 亚洲天堂国产精品一区在线| 日韩精品青青久久久久久| 一区二区三区四区激情视频 | 国产精品国产三级国产av玫瑰| 很黄的视频免费| 一区二区三区激情视频| 免费观看人在逋| 国产大屁股一区二区在线视频| 日本免费一区二区三区高清不卡| 男女视频在线观看网站免费| 午夜免费男女啪啪视频观看 | 中亚洲国语对白在线视频| 中文字幕精品亚洲无线码一区| 欧美激情国产日韩精品一区| 国产欧美日韩精品亚洲av| 亚洲国产精品久久男人天堂| 嫩草影视91久久| 欧美成人一区二区免费高清观看| 看黄色毛片网站| 他把我摸到了高潮在线观看| 精品无人区乱码1区二区| 国产精品98久久久久久宅男小说| 亚洲自偷自拍三级| 女的被弄到高潮叫床怎么办 | 欧美成人a在线观看| 亚洲中文日韩欧美视频| 亚洲综合色惰| 欧美3d第一页| 国产极品精品免费视频能看的| 最近中文字幕高清免费大全6 | 成人特级黄色片久久久久久久| 我要看日韩黄色一级片| 老女人水多毛片| 婷婷色综合大香蕉| 国内精品宾馆在线| 国产乱人伦免费视频| 最近视频中文字幕2019在线8| 老师上课跳d突然被开到最大视频| 亚洲乱码一区二区免费版| 午夜精品在线福利| 老熟妇乱子伦视频在线观看| 91久久精品国产一区二区三区| 国产在视频线在精品| 亚洲乱码一区二区免费版| 男人和女人高潮做爰伦理| 黄色一级大片看看| a级毛片免费高清观看在线播放| 中文字幕人妻熟人妻熟丝袜美| 无遮挡黄片免费观看| 最近在线观看免费完整版| 成年版毛片免费区| 国产精品自产拍在线观看55亚洲| 国产高清三级在线| 日韩精品有码人妻一区| 色视频www国产| 国内少妇人妻偷人精品xxx网站| 国模一区二区三区四区视频| 国产伦一二天堂av在线观看| 国产欧美日韩精品亚洲av| 日本在线视频免费播放| 村上凉子中文字幕在线| 国产真实乱freesex| 中文亚洲av片在线观看爽| 亚洲av熟女| 欧美日本视频| 男女边吃奶边做爰视频| 国产成人a区在线观看| 日韩精品有码人妻一区| 亚洲欧美精品综合久久99| 精品人妻1区二区| 成人av在线播放网站| 亚洲专区国产一区二区| 狠狠狠狠99中文字幕| 搡老熟女国产l中国老女人| 国产一区二区在线av高清观看| 观看免费一级毛片| 国产精品av视频在线免费观看| 国产精品久久久久久久久免| 日本与韩国留学比较| 一级黄片播放器| 国产精品综合久久久久久久免费| 久久九九热精品免费| 又紧又爽又黄一区二区| 女人十人毛片免费观看3o分钟| 国产精品亚洲美女久久久| 欧美一区二区国产精品久久精品| 亚洲av成人精品一区久久| 亚洲综合色惰| 别揉我奶头 嗯啊视频| 搡老岳熟女国产| 亚洲欧美精品综合久久99| 成人三级黄色视频| 干丝袜人妻中文字幕| 自拍偷自拍亚洲精品老妇| 美女高潮的动态| 精品人妻1区二区| 又爽又黄无遮挡网站| 91av网一区二区| 很黄的视频免费| 国产乱人视频| 可以在线观看的亚洲视频| 日韩欧美 国产精品| 国内精品宾馆在线| 国产精品亚洲一级av第二区| 欧美成人a在线观看| 精品不卡国产一区二区三区| 69av精品久久久久久| 国产毛片a区久久久久| 美女xxoo啪啪120秒动态图| 亚洲中文日韩欧美视频| 长腿黑丝高跟| 久久久久久久亚洲中文字幕| 深夜精品福利| 免费在线观看日本一区| 日本精品一区二区三区蜜桃| 午夜久久久久精精品| 久久婷婷人人爽人人干人人爱| 亚洲男人的天堂狠狠| 国内揄拍国产精品人妻在线| 亚洲av不卡在线观看| 精品久久久久久久久亚洲 | 国产精品综合久久久久久久免费| 日本欧美国产在线视频| 日韩中文字幕欧美一区二区| 日韩欧美精品免费久久| 在线a可以看的网站| 成人无遮挡网站| 国产精品1区2区在线观看.| 欧美zozozo另类| 国产毛片a区久久久久| 免费搜索国产男女视频| 一个人看的www免费观看视频| 免费人成视频x8x8入口观看| 亚洲最大成人中文| 一夜夜www| 久久久久国产精品人妻aⅴ院| 国产高清有码在线观看视频| 欧美日韩瑟瑟在线播放| 99国产精品一区二区蜜桃av| 欧美激情久久久久久爽电影| 男女视频在线观看网站免费| 欧美最新免费一区二区三区| 亚洲五月天丁香| 精品不卡国产一区二区三区| 一a级毛片在线观看| 精品99又大又爽又粗少妇毛片 | 美女高潮的动态| 精品久久久久久,| 国产v大片淫在线免费观看| av国产免费在线观看| 99热这里只有精品一区| 国产高清激情床上av| 中文字幕人妻熟人妻熟丝袜美| 91麻豆精品激情在线观看国产| 级片在线观看| 88av欧美| 欧美最黄视频在线播放免费| 有码 亚洲区| 真人做人爱边吃奶动态| a级毛片a级免费在线| 久久久久久久久久成人| 亚洲av第一区精品v没综合| 色综合色国产| 91久久精品国产一区二区成人| h日本视频在线播放| 在线观看午夜福利视频| 亚洲成人免费电影在线观看| 精品一区二区三区视频在线| 亚洲精品亚洲一区二区| 欧美丝袜亚洲另类 | 三级毛片av免费| 国产真实伦视频高清在线观看 | 精华霜和精华液先用哪个| 国产久久久一区二区三区| 美女高潮喷水抽搐中文字幕| 成人欧美大片| 我要搜黄色片| 在线观看66精品国产| 成人av在线播放网站| 波多野结衣高清无吗| 老司机午夜福利在线观看视频|