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

    The lubrication performance of water lubricated bearing with consideration of wall slip and inertial force*

    2017-03-09 09:09:26ZhongliangXie解忠良NaTa塔娜ZhushiRao饒柱石

    Zhong-liang Xie (解忠良), Na Ta (塔娜), Zhu-shi Rao (饒柱石)

    Laboratory of Vibration, Shock and Noise, Shanghai Jiao Tong University, Shanghai 200240, China, E-mail: xiezlsj2011@sjtu.edu.cn

    (Received July 30, 2015, Revised September 30, 2016)

    The lubrication performance of water lubricated bearing with consideration of wall slip and inertial force*

    Zhong-liang Xie (解忠良), Na Ta (塔娜), Zhu-shi Rao (饒柱石)

    Laboratory of Vibration, Shock and Noise, Shanghai Jiao Tong University, Shanghai 200240, China, E-mail: xiezlsj2011@sjtu.edu.cn

    (Received July 30, 2015, Revised September 30, 2016)

    The lubrication mechanism and the performance parameters with consideration of wall slip and inertial force are studied in this paper. Based on the modified Reynolds equation, the finite difference method is used to study the lubrication mechanism and the performance. Effects of the wall slip and the inertial force on the performance parameters are obtained, and found in good agreement with the results of FLUENT. It is shown that the wall slip and the inertial force do not significantly change the distribution of the pressure, the load capacity and the friction force. The inertial force slightly increases the pressure and the load capacity by 1.2% and 4.8%, while the wall slip reduces them by 8.0% and 17.85%. The wall slip and the inertial force increase the friction by about 15.98%, 2.33%, respectively. Compared with the wall slip, the inertial force is smaller, but cannot be neglected.

    Water lubricated bearing, wall slip, inertial force effect, lubrication performance

    Introduction

    The fluid flow boundary condition is one of the most important factors which determines the fluid hydrodynamic performance and characteristics. In the classical fluid mechanics and the lubrication theory the so-called “no slippage boundary conditions” are often used, namely, no wall slip occurs on the solidliquid interface, and the relative speed of motion between the fluid molecules on the solid surface and the solid interface is equal to zero. This hypothesis is verified by experiments in a macroscopic sense, and has been widely used in the theoretical analysis and the experimental researches. However, in recent years, with the development of the micro-nanometer science and technology, it is found that no slip boundary conditions are no longer valid under certain conditions. Boundary slip will occur in many instances[1,2], therefore affects the fluid behaviors and hydrodynamic cha-racteristics. The influence of the wall slip on the lubrication performance of the water lubricated bearing comes to be an important issue. For instance, Spikes et al.[3]analyzed the influence of slippage on the fluid dynamic behaviors with consideration of the wall slip only on the static slider surface. It is found that the bearing load capacity is exactly half of that with no slippage when the limiting shear stress is equal to zero, but the corresponding friction is reduced by several orders of magnitude.

    Consequently, Spikes[4,5]proposed the thought of designing the bearing of a very low coefficient of friction, believing that the existence of the boundary slip would reduce the friction coefficient of the bearing surface, but might have a more complicated influence on the fluid hydrodynamic characteristics. Ma et al.[6,7]et al. studied the influence of the wall slip on the hydrodynamic lubrication performance of a twodimensional journal bearing (finite length journal bearing) based on the limiting shear stress model. Fortier et al.[8,9], Wang et al.[10]and Wu et al.[11]put forward the concept of heterogeneous slip, that is, the wall slip in different positions of the same solid surface of different slippage properties (i.e., with different limiting shear stresses). It goes beyond the framework of the classical lubrication theory[12], which holds that onlythe convergent gap can make the journal bearing generate the hydrodynamic pressure.

    The water lubricated bearing is widely used in ships, water pumps and other mechanical systems for its various advantages, such as no pollution, wide source, safety and fire resistance. Its lubrication performance, reliability and safety help the safe and stable operation of the mechanical system. Therefore, the research and the improvement of the lubrication performance are of direct significance for the promotion and application of the water lubricated bearing.

    Based on the modified Reynolds equation with consideration of both the wall slip and the inertial force, this paper investigates the lubrication mechanism of the water lubricated bearing through numerical methods, focusing on the pressure and thickness distributions as well as the variations of the bearing capacity and the friction force. The effects of the interface slip and the inertial force on the lubrication performance are analyzed, and the results are found to be consistent with the finite element analysis by the software FLUENT. This will help the structure design and optimization of the water lubricated bearing.

    1. Derivation of modified Reynolds equation

    1.1Theoretical basis of wall slip

    It is shown[13,14]that the wall slip only occurs in the surface with smaller limiting shear stress when there exists a relative motion between two surfaces with different limiting shear stresses.

    Fig.1 Schematic diagram of wall slip model

    Figure 1(a) shows the schematic diagram of a wall slip model of the water lubricated bearing. Assuming that when the speed meets the criteriau=uc, Surface 1 just reaches the critical state of the interface slip, which means that the interface limiting shear stress is reached, correspondingly, the critical shear stressτcisτc=υ(uc/h)(υis the viscosity). With the increase of velocity on Surface 2, the wall slip starts to occur on Surface 1, the wall slip velocityusisus=u?uc,his the film thickness,ucthe critical slip velocity.

    If there is no wall slip on any of the surfaces, the wall slip velocity is equal to zero, i.e.us=0. When the wall slip begins to occur, the wall slip velocity is equal to the difference between the velocity of Surface 2 and the critical slip velocity, i.e.,us=u?uc. It was proposed[15,16]that when the interface tension of the lubricant medium becomes larger than that of the friction pair material, then the interface slip is very likely to happen.

    Table 1 shows the interfacial tension of commonly used materials. The interfacial tension of a water lubricated rubber bearing is much smaller, the adhesion fracture of water molecular bonds occurs more easily. Similarly, the wall slip may also exist on the solid-liquid interface of the water lubricated Teflon bearing. Therefore, the influence of the wall slip on the lubrication performance should be fully considered in the lubrication mechanism study for the water lubricated non-metallic or polymer friction pair material bearing.

    Table 1 Interfacial tension of commonly used materials

    Table 2 The ratio of viscous force to inertial force of commonly used lubricants

    1.2Theoretical basis of inertial force effect

    Generally speaking, Reynolds number is very small for conventional oil based lubricants. Table 2 shows the ratio of the viscous force to the inertial force of commonly used lubricants. It is large for theconventional lubricants, so the inertial force effect is almost negligible. Under the same conditions, this ratio is relatively small for low viscosity lubricants, so the inertial force effect cannot be ignored. When the Reynolds number is close to the critical Reynolds number, the impact of the inertial force on the fluid flow will become more and more significant, the modified Reynolds equation with consideration of the inertial force can more properly describe the hydrodynamic characteristics of the water film.

    Fig.2 Schematic diagram of water lubricated bearing configuration

    1.3Derivation of modified Reynolds equation

    The modified Reynolds equation with consideration of both the wall slip and the inertial force is derived in the Cartesian coordinates based on the above theoretical analysis. The schematic diagram of the water lubricated bearing configuration is shown in Fig.2, and a simplified water lubricated bearing configuration in the Cartesian coordinates is shown in Fig.1(b). Surface 1 moves with speeduin the positivexdirection, while Surface 2 is stationary. The wall slip occurs on Surface 1, the critical limiting shear stress isτc, the critical slip velocity isuc. Assuming that the film thicknesshcan vary in thex,ydirections, but is small enough inzdirection so that the lubrication approximation is valid and the transient effect in the Navier-Stokes equations can be ignored. Thermal effects are not considered in the analysis, as the isothermal condition, and the viscosity and density of water regard as constants, without any changes with the temperature and pressure in the film thickness direction.

    Using the above assumption, thex-component of the Navier-Stokes equation becomes

    With the boundary conditions, on Surface1

    and on Surface 2

    Solving Eq.(1), subject to the boundary conditions Eqs.(2) and (3), we obtain the velocity component inxdirection

    Similarly, the velocity component inydirection is

    whereIx,Iyare the momentum components inx,ydirections,ux,uyare the velocity components inx,ydirections. The volumetric flow rate per unit length inxdirection is

    Similarly, the volumetric flow rate per unit length inydirection is

    whereqx,qyare the volumetric flow rates inx,y-directions.

    Conservation of mass requires

    Inserting Eq.(6) and Eq.(7) into Eq.(8), we obtain

    Equation (9) is made dimensionless by defining the following dimensionless variables. The resulting dimensionless modified Reynolds equation becomes

    Linear combinations of these components represent modified Reynolds equations with different contributions as follows:

    The bearing load capacity of the water film can be calculated as follows:

    The friction force can be calculated as follows

    1.4Numerical solution

    These modified Reynolds equations are solved numerically by the finite difference method[12], with the discretization of the partial differential coefficients using the second order central difference. The natural rupture boundary conditions, namely, the Reynolds boundary conditions, are used here.

    2. Results and discussions

    Through analysis of the above theory, a finite length water lubricated bearing of angle of 360ois analyzed. Calculation starts from the thinnest area of the water film, i.e., assuming that the wall slip first occurs in the thinnest area of the film.

    Fig.3 (Color online) The dimensionless pressure distributions with consideration of wall slip and inertial force of water film

    Table 3 Characteristic parameters of circumferential pressure distribution of water film

    2.1Pressure distribution

    The dimensionless three-dimensional pressure distribution with consideration of the wall slip and the inertial force is shown in Fig.3(a). Figure 3(b) shows the dimensionless circumferential pressure distribution on the water film in four cases, whereθrepresents the circumferential coordinate,Lrepresents the axial coordinate. Pressure-I represents the result with consideration of the wall slip alone, pressure-II represents the result with consideration of both the wall slip and the inertial force, pressure-III represents the result without consideration of the wall slip and the inertial force effect, pressure-IV represents the result with consideration of the inertial force alone.

    From Figs.3(a), 3(b) we can see that the pressure distribution shape is not significantly affected after the consideration of the wall slip and the inertial force. The inertial force slightly increases the peak pressure by about 1.2%, while the interface slip significantly decreases the peak pressure by about 8.0%.

    Figure 3(b) shows the effect of the wall slip on the pressure boundary conditions, and also gives the rupture boundary of the water film. From Fig.3(b) we can obtain the maximum pressure and the rupture angle of the water film in four different cases. Characteristic parameters of the dimensionless circumferential pressure distribution of the water film are presented in Table 3. The water film rupture angle is about258owithout consideration of the wall slip and the inertial force, but the rupture angle increases to about 263owith consideration of the inertial force. Obviously, the pressure boundary expands outward and the entire water film lubrication region expands. Therefore, on one hand, the inertial force increases the film pressure peak, but on the other hand, it decreases the lubrication area and optimizes the lubrication state. To some extent, the inertial force effect is helpful for the lubrication performance of the bearing. The influence of the wall slip is just opposite.

    The inertial force only influences the pressure values to a small margin, as compared with the influence of the wall slip, the influence of the inertial force on the water film hydrodynamic pressure is very limited, but it does influence the water film pressure boundary, the lubrication region and state. Furthermore, the calculation shows that, with the increase of the eccentri city ratio, the corresponding hydrodynamic pressure of the water film increases, and the wall slip effect becomes more significant. So is the reduction of the pressure peak, the maximum reduction amplitude may reach as high as 12%. And, the wall slip is first observed in the vicinity of the area of the minimum film thickness.

    2.2Film thickness distributions

    The dimensionless three-dimensional film thickness distribution with consideration of the wall slip and the inertial force is shown in Fig.4(a). Figure 4(b) shows the dimensionless circumferential film thickness distribution of the water film in four cases.

    Fig.4 (Color online) The dimensionless film thickness of water film with consideration of wall slip and inertial force of water film

    From Figs.4(a), 4(b) we can see that the film thickness distribution shape is not significantly affected after the consideration of both the wall slip and the inertial force. The wall slip slightly increases the film thickness values by about 1.8%, The film thickness curve reaches its minimum value at about 180.88owith consideration of the wall slip, but the value is180owithout its consideration. There is a deviation of0.88obetween these two cases. The esse-nce of the wall slip is the adhesive fracture or the fracture of molecular bonds. Therefore, the wall slip does not evidently change the viscosity and the density of the lubricant as well as the shape of the water film thickness. The wall slip firstly occurs in the maximum hydrodynamic pressure region, i.e. the area of the thinnest water film. Because under the same condition, the shear stress in the region of the thinnest film is the first to reach or even exceed the interface limiting shear stress, then the adhesion fracture happens, at the same time, the wall slip occurs.

    The inertial force does not influence the film thickness values and the curve distribution trend. Essentially speaking, the inertial force is a conservative force and it does not work in the shaft rotation process. There is no energy transfer and no variation of the viscosity of lubricant. So the inertial force does not change the shape and the size of the film thickness distribution curve.

    2.3Bearing load capacity

    Figures 5(a), 5(b) show the dimensionless bearing load capacity and the friction force of the water film in four cases, respectively.εrepresents the eccentricity ratio of the bearing.

    Fig.5 (Color online) The dimensionless bearing load capacity and friction force of the water film

    From Fig.5(a) we can see that the bearing loading capacity increases along with the eccentricity, i.e., the larger the eccentricity, the larger the bearing loading capacity. The water film bearing loading capacity distribution is not significantly affected after consideration of both the wall slip and the inertial force. The inertial force slightly increases the bearing capacity, with the maximum increased amplitude of about 4.8%, while the wall slip evidently decreases the bearing capacity, with the reduced amplitude increasing along with the eccentricity, and with the maximum amplitude of about 17.85%. The bearing capacity with consideration of both the wall slip and the inertial force takes a value between that without consideration of the wall slip and the inertial force and that with consideration of the wall slip alone, but that value is much closer to the former. This fully illustrates the fact that compared to the wall slip, the influence of the inertial force on the lubrication performance is smaller, but cannot be neglected. From Fig.5(b) we can see that the curve shape is similar to that of the bearing load capacity. Mathematically speaking, compared to the wall slip, the influence of the inertial force on the friction force is smaller to some extent, but cannot be neglected.

    2.4Finite element analysis of wall slip

    A certain type of water lubricated bearing is analyzed by using the CFD software FLUNET. First we establish the water film model in the CAD software UG, then import it into the preprocessing module Gambit to mesh and fix the boundary conditions, and then import the mesh files into FLUENT for the finite element analysis. Table 4 shows the specific parameters of the water lubricated bearing. The water properties at 20oC as listed in Table 5 are employed.

    Table 4 Basic parameters of the bearing

    Generally speaking, the mesh plays an important role in the calculation. We have repeated the calculation several times with different meshes, to finally determine a proper mesh. We use a structural grid, for the mesh refinement and optimization. In the mesh, 5 divisions are used on the inlet and outlet surfaces along the radial direction, and 100 divisions are used for the edge on the inlet and outlet surfaces along the circumferential direction. The mesh is automatically divided along the axial water film direction. It is usually closely correlated with the length of the bearing. Altogether, there are about 17 952 nodes and 14 520 elements for the model (see in Figs.6(a), 6(b)).

    Figure 7 shows the hydrodynamic water film model and the boundary conditions. The hydrodynamic flow in the wedge gap is considered to be incompressible, in a steady state and isothermal. In the preprocessing modelling, we only build the flow model between the bearing and the journal surface.

    Table 5 Water properties at 20oC

    Fig.6 (Color online) Mesh models of the hydrodynamic water film

    Fig.7 (Color online) Water film model and boundary conditions

    The finite element software ANSYS 15.0 package is used in the simulation with the CFD models. Pro-processed models are incorporated into WORKBENCH. SIMPLEC algorithm is adopted in the calculation. The limiting shear stress of water is 1 806.5 Pa. The corresponding parameters are set in FLUENT and the simulations of two cases are conducted according to whether considering the interface slip or not. The pressure contours are obtained in the post processing module. Figure 8(a) shows the pressure contour of the water film without considering the wall slip. Figure 8(b) shows the pressure contour of the water film with consideration of the wall slip.

    Fig.8 (Color online) The pressure contour of water film without/with consideration of wall slip

    In order to further study the influence of the wall slip on the circumferential pressure distribution clearly, an iso-surface is plotted in the axial middle point of the film. The pressure data are calculated on the intersection of the iso-surface and the water film surface. Note that we ignore the negative pressure zone, and only consider the positive pressure area data in the investigation. Then the data are incorporated into the MATLAB to plot the pressure distribution curve[17,18].

    Figure 9 hows the circumferential pressure distribution of the water film in four cases.Xaxis represents the circumferential coordinate of the bearing,Yaxis represents the water film pressure; P-I represents the result by the finite difference method without considering the wall slip, P-II represents the result by Fluent without considering the wall slip, P-III represents the result by the finite difference method with consideration of the wall slip, and P-IV represents the result by FLUENT with consideration of the wall slip.From Fig.9 we can see that the water film circumferential pressure increases nonlinearly and reaches its maximum value at about 0.063 m, then decreases rapidly to zero. With consideration of the wall slip, the peak pressure decreases visibly, with the maximum amplitude being reduced by about 3.4%. The wall slip first occurs in the vicinity of the region of the maximum hydrodynamic pressure for the shear stress here is the first to reach the interface limiting shear stress. The finite element simulation results are consistent with the theoretical results[19], with a very small deviation between them, which almost can be neglected. Consistency proves the correctness of the theoretical model[19].

    Fig.9 (Color online) The circumferential pressure distribution

    In the subsequent investigations[20,21], related experiments will be conducted systematically in order to obtain experimental results to verify the theory and the lubrication mechanism for the water lubricated plain journal bearing with consideration of wall slip[22].

    3. Conclusions

    This paper establishes a water lubricated bearing fluid dynamics model, derives the modified Reynolds equation with consideration of the wall slip and the inertial force. By numerical analysis methods, the influence of the wall slip and the inertial force on the lubrication performance of the water lubricated bearing is studied. The results are compared with those of the finite element software Fluent. The following conclusions are drawn:

    (1) The wall slip and the inertial force have a certain impact on the hydrodynamic characteristics of the water film. The lubrication state and the lubrication performance can be more accurately described by the modified Reynolds equation with consideration of the wall slip and the inertial force.

    (2) The wall slip slightly increases the film thickness, but does not change its shape. The inertial force has almost no effect on the film thickness distribution, its shape, value or size.

    (3) The wall slip and the inertial force do not change the trend of the water film pressure, the bearing capacity and the friction force with eccentricity, but only change their values.

    (4) Comparison between the finite element analysis by Fluent and the theoretical analysis is done for a certain actual water lubricated bearing with consideration of the wall slip. Simulation results are consistent with theoretical results.

    [1] Kalin M., Polajnar M. The Effect of wetting and surface energy on the friction and slip in oil-lubricated contacts [J].Tribology letter, 2013, 52(2): 185-194.

    [2] Neto C., Evans D. R., Bonaccurso E. et al. Boundary slip in Newtonian liquids: A review of experimental studies [J].Reports on Progress in Physics, 2005, 68(12): 28-59.

    [3] Spikes H. A. The half-wetted bearing. Part 1: Extended Reynolds equation [J].Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, 2003, 217(1): 1-14.

    [4] Spikes H. A. The half-wetted bearing. Part 2: Potential application in low load contacts [J].Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, 2003, 217(1): 15-26.

    [5] Spikes H., Granick S. Equation for slip of simple liquids at smooth solid surfaces [J].Langmuir, 2003, 19(12): 5065-5071.

    [6] Ma G. J., Wu C. W., Zhou P. Wall slip and hydrodynamics of two-dimensional journal bearing [J].Tribology International, 2007, 40(7): 1056-1066.

    [7] Ma G. J., Wu C. W., Zhou P. Influence of wall slip on the hydrodynamic behavior of a two-dimensional slider bearing [J].Acta Mechanica Sinica, 2007, 23(6): 655-661.

    [8] Fortier A. E., Salant R. F. Numerical analysis of a journal bearing with a heterogeneous slip/no-slip surface [J].Journal of Tribology, 2005, 127(4): 820-825.

    [9] Salant R. F., Fortier A. E. Numerical analysis of a slider bearing with a heterogeneous slip/no-slip surface [J].Tribology Transaction, 2004, 47(3): 328-334.

    [10] Wang L. L., Lu C. H., Wang M. The numerical analysis of the radial sleeve bearing with combined surface slip [J].Tribology International, 2012, 47: 100-104.

    [11] Wu C. W., Wu C. D., Ma G. J. Low friction and high load support capacity of slider bearing with a mixed slip surface [J].Journal of Tribology, 2006, 128(4): 904-907.

    [12] Szeri A. Z. Fluid film lubrication: Theory and design [M]. New York, USA: Cambridge University Press, 2005.

    [13] Guo F., Wong P, Geng M. Occurrence of wall slip in elastohydrodynamic lubrication contacts [J].Tribology Letter, 2009, 34(2): 103-111.

    [14] Wu C. W. Performance of hydrodynamic lubrication journal bearing with a slippage surface [J].Industrial Lubrication Tribology, 2008, 60(6): 293-298.

    [15] Aurelian F., Patrick M., Mohamed H. Wall slip effects in (elasto) hydrodynamic journal bearings [J].Tribology Intertational, 2011, 44(7): 868-877.

    [16] Hatzikiriakos S., Dealy J. Wall slip of molten high density polyethylene. I. Sliding plate rheometer studies [J].Journal of Rheology, 1991, 35(4): 497-523.

    [17] Xie Z., Rao Z., Ta N. et al. Investigations on transitions of lubrication states for water lubricated bearing. Part I: Determination of friction coefficient and film thicknessratios [J].Industrial Lubrication and Tribology, 2016, 68(3): 404-415.

    [18] Xie Z., Rao Z., Ta N. et al. Investigations on transitions of lubrication states for water lubricated bearing. Part II: Further insight into the film thickness ratio lambda [J].Industrial Lubrication and Tribology, 2016, 68(3): 416-429.

    [19] Xie Z., Rao Z., Ta N. et al. Theoretical and experimental research on the friction coefficient of water lubricated bearing with consideration of wall slip effects [J].Mechanics and Industry, 2016, 17(1): 106-119.

    [20] Yang J. Sharp interface direct forcing immersed boundary methods: A summary of some algorithms and applications [J].Journal of Hydrodynamics, 2016, 28(5): 713-730.

    [21] Zhang X., Yin Z., Gao G. et al. Determination of stiffness coefficient of hydrodynamic water-lubricated plain journal bearings [J].Tribology International, 2015, 85: 37-47.

    [22] Illner T., Bartel D., Deters L. Determination of the transition speed in journal bearing consideration of bearing deformation [J].Tribology Internaional, 2015, 82: 58-67.

    * Biography:Zhong-liang Xie (1988-), Male, Ph. D. Candidate

    Zhu-shi Rao, E-mail: zsrao@sjtu.edu.cn

    18+在线观看网站| 久久久久久久亚洲中文字幕| 校园人妻丝袜中文字幕| 欧美bdsm另类| 日本免费在线观看一区| 成年女人在线观看亚洲视频| 国产中年淑女户外野战色| 多毛熟女@视频| 久久精品国产自在天天线| 国产精品久久久久久久电影| 久久狼人影院| 亚洲人成网站在线播| 中文乱码字字幕精品一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 久久久久精品性色| 久久久久久久精品精品| 少妇人妻久久综合中文| 99视频精品全部免费 在线| 国产精品麻豆人妻色哟哟久久| 国产 精品1| 美女国产视频在线观看| 免费观看在线日韩| 一级毛片aaaaaa免费看小| av国产久精品久网站免费入址| kizo精华| 啦啦啦啦在线视频资源| 春色校园在线视频观看| 黑人猛操日本美女一级片| 国产一区二区在线观看日韩| 寂寞人妻少妇视频99o| 久久av网站| 免费少妇av软件| 99久久人妻综合| 精品人妻熟女av久视频| 精品卡一卡二卡四卡免费| 欧美国产精品一级二级三级 | 国产亚洲精品久久久com| 日韩人妻高清精品专区| 看十八女毛片水多多多| 国产精品人妻久久久久久| 国产色爽女视频免费观看| 久久精品国产亚洲网站| 国产在视频线精品| 一区二区三区乱码不卡18| 亚洲va在线va天堂va国产| 国产高清有码在线观看视频| 男女国产视频网站| 亚洲一级一片aⅴ在线观看| 九九久久精品国产亚洲av麻豆| 亚洲av欧美aⅴ国产| 亚洲电影在线观看av| 亚洲av二区三区四区| 欧美激情国产日韩精品一区| 少妇精品久久久久久久| 国产欧美日韩一区二区三区在线 | 777米奇影视久久| 亚洲怡红院男人天堂| 欧美日韩视频精品一区| 成人毛片a级毛片在线播放| 日韩熟女老妇一区二区性免费视频| 少妇被粗大的猛进出69影院 | 亚洲av二区三区四区| www.av在线官网国产| 美女视频免费永久观看网站| 在线 av 中文字幕| 王馨瑶露胸无遮挡在线观看| 成人漫画全彩无遮挡| 亚洲国产精品专区欧美| 在线观看免费日韩欧美大片 | 亚洲国产精品999| 日韩电影二区| 亚洲av电影在线观看一区二区三区| 久久久久精品久久久久真实原创| 少妇人妻一区二区三区视频| 少妇被粗大猛烈的视频| 搡女人真爽免费视频火全软件| 亚洲精品色激情综合| 日日撸夜夜添| 亚洲欧美成人综合另类久久久| 如日韩欧美国产精品一区二区三区 | 色94色欧美一区二区| 国产av码专区亚洲av| 黄色配什么色好看| 麻豆成人av视频| 午夜免费观看性视频| 久久国产精品男人的天堂亚洲 | 深夜a级毛片| 欧美亚洲 丝袜 人妻 在线| 三级国产精品片| 美女内射精品一级片tv| 国产高清有码在线观看视频| 久久ye,这里只有精品| 九九久久精品国产亚洲av麻豆| 黄片无遮挡物在线观看| 黄片无遮挡物在线观看| 男女啪啪激烈高潮av片| 男女啪啪激烈高潮av片| a级毛片在线看网站| 看十八女毛片水多多多| 黄色一级大片看看| 国产伦在线观看视频一区| 97超碰精品成人国产| 久久精品熟女亚洲av麻豆精品| 精品视频人人做人人爽| 美女中出高潮动态图| 老熟女久久久| 国产精品麻豆人妻色哟哟久久| 观看免费一级毛片| 成人美女网站在线观看视频| 男女无遮挡免费网站观看| 大又大粗又爽又黄少妇毛片口| 精品久久久噜噜| 久久毛片免费看一区二区三区| 国产午夜精品一二区理论片| 午夜免费男女啪啪视频观看| 蜜桃在线观看..| 久久av网站| 精品酒店卫生间| 大香蕉久久网| 国产av码专区亚洲av| 九九在线视频观看精品| 大陆偷拍与自拍| 22中文网久久字幕| 伊人亚洲综合成人网| 伊人亚洲综合成人网| 久久久久久久久久成人| 亚洲丝袜综合中文字幕| 日韩欧美一区视频在线观看 | 国产午夜精品久久久久久一区二区三区| 婷婷色麻豆天堂久久| 一级,二级,三级黄色视频| 又大又黄又爽视频免费| 日韩中文字幕视频在线看片| 欧美少妇被猛烈插入视频| 欧美性感艳星| 啦啦啦视频在线资源免费观看| 男人爽女人下面视频在线观看| 91精品国产九色| 亚洲国产精品国产精品| 99久久人妻综合| 老女人水多毛片| 涩涩av久久男人的天堂| tube8黄色片| 日韩电影二区| av视频免费观看在线观看| av视频免费观看在线观看| av视频免费观看在线观看| 国产又色又爽无遮挡免| 免费看不卡的av| 国产精品人妻久久久久久| 国产伦理片在线播放av一区| 国产伦理片在线播放av一区| 久久久久久久久久成人| 精品人妻熟女毛片av久久网站| 一个人看视频在线观看www免费| 亚洲性久久影院| 久久精品国产亚洲av涩爱| 高清视频免费观看一区二区| 内射极品少妇av片p| 九草在线视频观看| 人妻制服诱惑在线中文字幕| 看十八女毛片水多多多| 九九在线视频观看精品| 三级国产精品欧美在线观看| 边亲边吃奶的免费视频| 男人舔奶头视频| 国产极品天堂在线| 欧美日韩视频精品一区| 少妇人妻 视频| 亚洲av成人精品一二三区| 我的女老师完整版在线观看| 街头女战士在线观看网站| 精品国产一区二区久久| 久久精品久久久久久噜噜老黄| 最新的欧美精品一区二区| 亚洲美女黄色视频免费看| 国产一区二区在线观看日韩| 99久久精品热视频| 亚洲一级一片aⅴ在线观看| 2018国产大陆天天弄谢| 搡老乐熟女国产| 最近中文字幕高清免费大全6| 亚洲成人手机| 春色校园在线视频观看| 国产精品久久久久久精品电影小说| 男人和女人高潮做爰伦理| 国产精品嫩草影院av在线观看| 中文乱码字字幕精品一区二区三区| 日韩 亚洲 欧美在线| 91在线精品国自产拍蜜月| 美女大奶头黄色视频| 欧美精品人与动牲交sv欧美| 亚洲精品久久久久久婷婷小说| 免费观看a级毛片全部| 中文欧美无线码| 日本wwww免费看| 我的女老师完整版在线观看| 国产亚洲欧美精品永久| 在线观看国产h片| 免费观看性生交大片5| 一级毛片我不卡| 男人爽女人下面视频在线观看| 99热这里只有是精品在线观看| 国产亚洲一区二区精品| 亚洲精品国产av蜜桃| 日本91视频免费播放| 少妇被粗大的猛进出69影院 | 成年人免费黄色播放视频 | 国产欧美日韩综合在线一区二区 | 99久久综合免费| 国产精品熟女久久久久浪| 特大巨黑吊av在线直播| 国产日韩欧美在线精品| 尾随美女入室| 熟女电影av网| 午夜av观看不卡| 国产男人的电影天堂91| 久久久久久久精品精品| 一区二区三区乱码不卡18| 嘟嘟电影网在线观看| 欧美成人午夜免费资源| 男人添女人高潮全过程视频| 色5月婷婷丁香| 亚洲欧美精品自产自拍| 国产一区二区三区av在线| 久久久久久伊人网av| 中国国产av一级| 国产在视频线精品| 久久精品国产鲁丝片午夜精品| 蜜桃在线观看..| 插逼视频在线观看| 男人舔奶头视频| 美女国产视频在线观看| 国产亚洲午夜精品一区二区久久| 免费观看性生交大片5| 麻豆成人午夜福利视频| 久久久久久久国产电影| 在线播放无遮挡| 日本免费在线观看一区| 午夜av观看不卡| 日韩大片免费观看网站| 色视频在线一区二区三区| 国产精品一区www在线观看| 在线亚洲精品国产二区图片欧美 | 欧美丝袜亚洲另类| av播播在线观看一区| 高清黄色对白视频在线免费看 | 久久久久精品性色| 欧美成人午夜免费资源| 18+在线观看网站| 九九在线视频观看精品| 汤姆久久久久久久影院中文字幕| 一本色道久久久久久精品综合| 精品久久久精品久久久| av在线app专区| 国产有黄有色有爽视频| 看非洲黑人一级黄片| 91久久精品国产一区二区三区| 啦啦啦中文免费视频观看日本| 我要看黄色一级片免费的| 亚洲不卡免费看| 欧美激情国产日韩精品一区| 亚洲av日韩在线播放| 精品久久久久久电影网| a 毛片基地| 欧美高清成人免费视频www| 两个人的视频大全免费| 亚洲伊人久久精品综合| 欧美国产精品一级二级三级 | 熟女av电影| 久久99一区二区三区| 精品国产乱码久久久久久小说| 成人特级av手机在线观看| 夜夜骑夜夜射夜夜干| 精品久久久噜噜| 人妻 亚洲 视频| 色网站视频免费| 日韩伦理黄色片| 18禁在线播放成人免费| 久久久久精品性色| 久久久a久久爽久久v久久| 国产av精品麻豆| av免费在线看不卡| 一二三四中文在线观看免费高清| 狂野欧美激情性bbbbbb| 精品亚洲乱码少妇综合久久| 黄色配什么色好看| 日产精品乱码卡一卡2卡三| 亚洲内射少妇av| 欧美日韩国产mv在线观看视频| 成人二区视频| 亚洲成色77777| 我的女老师完整版在线观看| 亚洲欧美日韩卡通动漫| 男人狂女人下面高潮的视频| 亚洲无线观看免费| 亚洲精品,欧美精品| 久久精品久久精品一区二区三区| 内地一区二区视频在线| 免费大片黄手机在线观看| 在线观看免费日韩欧美大片 | 国产av码专区亚洲av| 婷婷色麻豆天堂久久| 日韩三级伦理在线观看| 精品久久久精品久久久| 熟妇人妻不卡中文字幕| 欧美老熟妇乱子伦牲交| 边亲边吃奶的免费视频| 亚洲精品色激情综合| 欧美日韩国产mv在线观看视频| 国产精品蜜桃在线观看| 久久亚洲国产成人精品v| 日韩一本色道免费dvd| 亚洲av福利一区| 欧美日韩国产mv在线观看视频| 欧美日韩亚洲高清精品| 少妇的逼水好多| 欧美xxⅹ黑人| 国产成人91sexporn| 日日爽夜夜爽网站| 天堂俺去俺来也www色官网| 精品一品国产午夜福利视频| 黑丝袜美女国产一区| 99热国产这里只有精品6| 偷拍熟女少妇极品色| 51国产日韩欧美| av在线老鸭窝| 另类亚洲欧美激情| 国产伦在线观看视频一区| 国产免费一级a男人的天堂| 伊人亚洲综合成人网| 亚洲国产欧美在线一区| 精品卡一卡二卡四卡免费| 中文字幕制服av| 亚洲精品一二三| 这个男人来自地球电影免费观看 | 国产爽快片一区二区三区| 伊人久久精品亚洲午夜| 久久99蜜桃精品久久| 青春草亚洲视频在线观看| 亚洲人与动物交配视频| 欧美97在线视频| 午夜福利网站1000一区二区三区| 丰满人妻一区二区三区视频av| 国产在视频线精品| 精品亚洲成国产av| 高清午夜精品一区二区三区| 国产伦精品一区二区三区视频9| 久久国内精品自在自线图片| 十八禁高潮呻吟视频 | 中国国产av一级| 亚洲成人一二三区av| 久久婷婷青草| 伊人亚洲综合成人网| 夫妻午夜视频| 久久人人爽人人片av| 国产又色又爽无遮挡免| 街头女战士在线观看网站| 视频中文字幕在线观看| 亚洲国产欧美日韩在线播放 | xxx大片免费视频| 国产老妇伦熟女老妇高清| 久久精品熟女亚洲av麻豆精品| 亚洲av欧美aⅴ国产| 2018国产大陆天天弄谢| 久久国内精品自在自线图片| 欧美97在线视频| 日韩欧美 国产精品| 伦理电影免费视频| 国产精品国产三级国产专区5o| 高清视频免费观看一区二区| 成人黄色视频免费在线看| 最近中文字幕高清免费大全6| 纯流量卡能插随身wifi吗| 中文字幕人妻丝袜制服| 3wmmmm亚洲av在线观看| 免费观看性生交大片5| 在线亚洲精品国产二区图片欧美 | 一级爰片在线观看| 人体艺术视频欧美日本| 日韩三级伦理在线观看| 欧美日韩av久久| 亚洲性久久影院| 欧美日韩精品成人综合77777| 卡戴珊不雅视频在线播放| 亚洲第一av免费看| 中文天堂在线官网| 久久综合国产亚洲精品| 91精品国产九色| 国产高清不卡午夜福利| av又黄又爽大尺度在线免费看| 日韩一区二区三区影片| 欧美日韩综合久久久久久| 最近最新中文字幕免费大全7| 日韩成人伦理影院| 哪个播放器可以免费观看大片| 午夜福利视频精品| 在线观看一区二区三区激情| freevideosex欧美| 亚洲高清免费不卡视频| 黑人巨大精品欧美一区二区蜜桃 | 免费av中文字幕在线| 亚洲精品国产av蜜桃| a级毛片在线看网站| 久久人人爽av亚洲精品天堂| 午夜激情福利司机影院| 内射极品少妇av片p| 日本vs欧美在线观看视频 | 久久女婷五月综合色啪小说| 国产亚洲91精品色在线| 熟女人妻精品中文字幕| 久久久久久久久久成人| 亚洲欧洲国产日韩| 亚洲av免费高清在线观看| 日韩一区二区视频免费看| 国产精品成人在线| 国产成人a∨麻豆精品| 亚洲av成人精品一区久久| 亚洲激情五月婷婷啪啪| 少妇裸体淫交视频免费看高清| 欧美成人午夜免费资源| 中文字幕人妻丝袜制服| 亚洲精品国产色婷婷电影| 美女视频免费永久观看网站| 久久女婷五月综合色啪小说| xxx大片免费视频| 99热全是精品| 国产成人精品福利久久| 亚州av有码| 日韩 亚洲 欧美在线| 伊人久久精品亚洲午夜| 亚洲av成人精品一二三区| 2021少妇久久久久久久久久久| 亚洲人成网站在线观看播放| 美女主播在线视频| 亚洲欧美清纯卡通| 午夜免费观看性视频| 深夜a级毛片| 嫩草影院入口| 精品一品国产午夜福利视频| 亚洲av男天堂| 男人爽女人下面视频在线观看| 日韩成人伦理影院| 亚洲欧美成人综合另类久久久| 亚洲国产色片| 亚洲在久久综合| 国产伦精品一区二区三区视频9| 亚洲四区av| 人妻 亚洲 视频| 视频区图区小说| 中文乱码字字幕精品一区二区三区| 80岁老熟妇乱子伦牲交| 国内揄拍国产精品人妻在线| 女性被躁到高潮视频| 天美传媒精品一区二区| 嫩草影院新地址| 80岁老熟妇乱子伦牲交| 亚洲国产精品专区欧美| 五月天丁香电影| 亚洲在久久综合| 亚洲自偷自拍三级| 麻豆精品久久久久久蜜桃| 国产精品一区二区在线观看99| 日本猛色少妇xxxxx猛交久久| 国产一级毛片在线| 国产免费福利视频在线观看| 亚洲情色 制服丝袜| 97在线人人人人妻| 最新的欧美精品一区二区| 亚洲婷婷狠狠爱综合网| 国产在线男女| 草草在线视频免费看| 狂野欧美激情性bbbbbb| 国产白丝娇喘喷水9色精品| 国产精品三级大全| 内地一区二区视频在线| 男女边摸边吃奶| av国产久精品久网站免费入址| 18+在线观看网站| 女人精品久久久久毛片| a级毛色黄片| 欧美三级亚洲精品| 校园人妻丝袜中文字幕| 国产色婷婷99| 少妇丰满av| 午夜激情福利司机影院| 色网站视频免费| 天堂8中文在线网| 国产色婷婷99| 日韩精品有码人妻一区| 国产成人aa在线观看| 亚洲成色77777| av天堂久久9| 国产色爽女视频免费观看| 欧美日韩在线观看h| 久久国内精品自在自线图片| 日本黄色日本黄色录像| 国产精品三级大全| 伊人久久国产一区二区| 久久99精品国语久久久| 日韩亚洲欧美综合| 日韩,欧美,国产一区二区三区| 免费大片18禁| 亚洲欧美成人精品一区二区| 九色成人免费人妻av| 在线观看一区二区三区激情| 少妇被粗大猛烈的视频| 秋霞在线观看毛片| 亚洲欧美一区二区三区黑人 | 国产成人免费观看mmmm| 欧美3d第一页| 22中文网久久字幕| 成人无遮挡网站| 麻豆成人av视频| 午夜久久久在线观看| 涩涩av久久男人的天堂| 欧美日韩在线观看h| 黄色视频在线播放观看不卡| 成人亚洲欧美一区二区av| 天美传媒精品一区二区| 日韩一区二区视频免费看| 成人国产av品久久久| www.av在线官网国产| 一级毛片电影观看| 国产一区有黄有色的免费视频| 国内揄拍国产精品人妻在线| 免费在线观看成人毛片| 人妻夜夜爽99麻豆av| 午夜老司机福利剧场| 日韩,欧美,国产一区二区三区| 男女国产视频网站| 精品久久久噜噜| 美女国产视频在线观看| 特大巨黑吊av在线直播| av在线观看视频网站免费| 国产老妇伦熟女老妇高清| 国产男女超爽视频在线观看| 免费观看性生交大片5| 久久久欧美国产精品| 精品久久久久久久久av| 乱码一卡2卡4卡精品| 国内揄拍国产精品人妻在线| 热re99久久精品国产66热6| 亚洲精品视频女| 中文乱码字字幕精品一区二区三区| 亚洲成人av在线免费| 国产黄片视频在线免费观看| 91精品一卡2卡3卡4卡| 精品久久久久久久久av| 欧美精品高潮呻吟av久久| 亚洲国产av新网站| 香蕉精品网在线| 中文欧美无线码| 国产美女午夜福利| 纵有疾风起免费观看全集完整版| 亚洲欧美清纯卡通| 夫妻性生交免费视频一级片| 国精品久久久久久国模美| 欧美+日韩+精品| 搡老乐熟女国产| av视频免费观看在线观看| 大片电影免费在线观看免费| 国产亚洲av片在线观看秒播厂| 亚洲av成人精品一二三区| 久久久久久久大尺度免费视频| 在线观看免费日韩欧美大片 | 午夜福利影视在线免费观看| 人人妻人人澡人人爽人人夜夜| 91久久精品国产一区二区成人| 成年人免费黄色播放视频 | 国产中年淑女户外野战色| 男女啪啪激烈高潮av片| 黄色视频在线播放观看不卡| 边亲边吃奶的免费视频| 亚洲欧美中文字幕日韩二区| 成人毛片60女人毛片免费| 国产老妇伦熟女老妇高清| 成人毛片60女人毛片免费| 日韩亚洲欧美综合| 日韩欧美 国产精品| 国产淫片久久久久久久久| 熟女电影av网| 欧美 亚洲 国产 日韩一| 制服丝袜香蕉在线| 久久久久网色| 男的添女的下面高潮视频| 欧美精品亚洲一区二区| 91午夜精品亚洲一区二区三区| 久久人人爽人人片av| av在线观看视频网站免费| 啦啦啦啦在线视频资源| 欧美日韩视频精品一区| 18禁动态无遮挡网站| 国产一区二区在线观看日韩| 久久久久久人妻| 99热网站在线观看| 中文字幕人妻熟人妻熟丝袜美| 99热网站在线观看| 老女人水多毛片| 精品午夜福利在线看| 成年女人在线观看亚洲视频| 国产精品欧美亚洲77777| 人妻 亚洲 视频| 国产精品一区二区在线不卡| 日日摸夜夜添夜夜添av毛片| 综合色丁香网| 精品酒店卫生间| 成人无遮挡网站| 日本-黄色视频高清免费观看| 热99国产精品久久久久久7| 国产女主播在线喷水免费视频网站| 亚洲欧美中文字幕日韩二区| 91午夜精品亚洲一区二区三区| 日韩一区二区三区影片| 国产精品无大码|