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

    Experimental study on vibration suppression in a rotor system under base excitation using an integral squeeze film damper

    2020-11-27 09:17:06YanWeiHeLidongZhuGangWangShengliDengZhe
    High Technology Letters 2020年4期

    Yan Wei (閆 偉), He Lidong, Zhu Gang, Wang Shengli, Deng Zhe

    (Beijing Key Laboratory of Health Monitoring and Self-Recovery for High-End Mechanical Equipment, Beijing University of Chemical Technology, Beijing 100029, P.R.China)

    Abstract

    Key words: Jeffcott rotor, dynamic characteristics, base excitation, integral squeeze film damper (ISFD), vibration suppression

    0 Introduction

    Turbomachinery, such as turbines, pumps, and compressors, is installed in transportation systems, including warships, submarines, and space vehicles. In their life cycle, they are often subjected to various sudden impact forces, whose specific impact depends on the working conditions and the external environment. Such impacts can be transmitted directly to the rotor-bearing system of turbomachinery through the base or foundation, causing rotor vibration. In severe cases, this may result in direct collision between rotor and bearing, seals and other stators, causing serious damage. Therefore, it is necessary to consider the rotor-support-base system as a whole in order to study the dynamic characteristics of rotor systems under base excitation.

    Although it is common for rotor systems to be subjected to inertial forces related to the fundamental motion, it seems that there is a lack of literature on suppressing the lateral vibration of rotors caused by the base motion. Lee et al.[1,2]proposed a transient response analysis method for rotor systems, and proved that the transient response of the rotor is sensitive to the duration of the impact. Sharma et al.[3,4]conducted numerical analysis to demonstrate that the electromagnetic actuator can significantly reduce the response amplitude of the rotor-shafting system under seismic excitation and improve its stability. Han et al.[5,6]explored the parameter instability of a flexible rotor-bearing system under time-periodic base angular motions. Wang et al.[7]investigated the response of a nonlinear bearing-rotor system under the interaction between unbalanced forces and base excitations. Duchemin et al.[8]conducted theoretical and experimental studies on the dynamic behavior of a flexible rotor system under base excitation (support motion). Song et al.[9]designed an active control system of a rotor-electromagnetic dynamic vibration absorber that was excited through the base. Under the combined action of the vibration absorber and the base excitation, multi-frequency dynamic vibration absorption in the rotor system was achieved. Fawzi and Fred[10]studied the dynamic characteristics of rigid rotors under the combined action of fundamental excitation and mass imbalance, and provided theoretical guidance for the design of rotating systems.

    A new type of damper with novel structure, named integral squeeze film damper (ISFD), appeared in the early 1990s and can provide excellent vibration control performance and effectively improve system stability. Compared to the traditional squeeze film damper, ISFD can effectively suppress rotor vibration, improve the stability of the rotor system, and achieve better vibration reduction performance[11-13]. Santiago and Andrés[14,15]investigated the unbalanced response characteristics of ISFD rotor systems with and without end seals, and revealed that the stiffness and damping of the ISFD had excellent linear characteristics under varying working conditions. Delgado et al.[16]studied the effect of vibration speed on the damping performance of ISFDs. The experiment proved that ISFDs can provide a higher damping coefficient, have better dynamic characteristics than traditional squeeze film dampers, and can effectively suppress the cavitation phenomenon of the oil film. Ertas et al.[17,18]investigated the effect of different end seal clearances on the dynamic coefficients of an ISFD. The hammering method is used to measure the variation trend of the damping coefficients of ISFD, which are found to be linearly correlated with the vibration velocity without cavitation. The experiments demonstrated that ISFDs can solve vibration problems by providing sufficient damping in turbomachinery, such as supercritical carbon dioxide turbochargers. Lu et al.[19-21]performed several experimental studies on the vibration control of gear shafting using ISFD. The experimental results proved that ISFDs can effectively suppress the impact vibration of gear meshing and can efficiently control the vibration in a wider frequency band.

    In this paper, a single-disk rotor test rig with base excitation is developed to study the response of the rotor system under various frequencies and excitation forces, as well as to explore the effect of the ISFDs on the dynamic characteristics of the rotor. The vibration suppression of a rotor system by ISFDs under the action of mass imbalance and base excitation is experimentally investigated, and the lateral vibration law of the shaft relative to the stator is studied.

    1 Integral squeeze film damper

    Squeeze film dampers (SFDs) are commonly used as vibration absorbers in aero-engines, since they provide good vibration reduction performance, and have been widely used in other rotating machinery. Traditional SFDs usually use a centering elastic support, which produces viscous damping by squeezing the oil film vibration modulation in an annular gap. This support provides excellent vibration reduction effect, reduces the rotor bending, reduces the dynamic load exerted on the bearing when the rotor passes through the critical speed, and improves the rotor system stability[22]. However, under severe conditions, such as heavy load and sudden unbalance, where the design scope of work is exceeded, traditional SFDs will exhibit a high degree of nonlinear characteristics, which will lead to nonlinear vibration phenomena, such as lock up, bistable response, non-synchronous precession, and even chaotic motion[23].

    Compared to traditional SFDs, the design method and performance of ISFDs are significantly different. ISFDs are fabricated using wire cut electrical discharge machining (WEDM) technology. The piecewise design of the squeeze film region prevents the circumferential flow of the squeeze oil film and solves the highly nonlinear problem of traditional SFDs.

    A schematic diagram of the ISFD structure is shown in Fig.1. The structure consists of 2 parts: the outer ring and the inner ring. The outer ring and the bearing seat are secured by a transition fit, while the inner ring and the outer ring of the rolling bearing are interference-fitted. An S-type elastomer connects the outer ISFD ring and the inner ISFD ring, and the gap between the outer ring and the inner ring forms the squeeze film region. In case of failure, where the rotor system vibrates, the vibration will be transmitted from the rolling bearing to the ISFD inner ring, and the squeezing action will generate a squeezing effect, such that the vibration of the rotor will be effectively controlled and the stability of the rotor system will be improved. The basic parameters of the ISFD structure are listed in Table 1.

    Fig.1 Schematic diagram of the ISFD structure

    Table 1 ISFD structure parameters

    Studies have shown that the damping force generated by the extruded oil film in an ISFD can be linearized, and the equivalent physical model can be observed in Fig.2[17].KijandCijis the stiffness and damping coefficient of the equivalent physical model, whereiandjrepresentxandydirections. Different stiffness and damping magnitudes can be obtained by altering the ISFD axial length, radial height, radial thickness, oil film clearance, S-type elastomer distribution angle, and other structural parameters.

    Fig.2 ISFD equivalent physical model

    2 Dynamic model

    Firstly, the Jeffcott rotor system is taken as the research object, and rotor dynamic models of rigid and ISFD support under base excitation are established, providing a theoretical basis for experimental research.

    2.1 Dynamic model of rigid support rotor

    A schematic diagram of the Jeffcott rotor system model of the rigid support structure under base excitation is shown in Fig.3. The two ends of the rotor adopt a rigid support structure, which is composed of a bearing seat, a rigid sleeve, and a deep-groove ball bearing. The excitation signal is output by the exciter, providing base excitation for the rotor system. In Fig.3(b),fBrepresents the nonlinear supporting force of the deep-groove ball bearings.

    Fig.3 Dynamic model of a rotor supported on ball bearings with the rigid sleeve subjected to base excitation

    The differential equations of motion of the system are wheremDis the concentrated mass on the disk;cDis the viscous damping of the disk caused by the aerodynamic effect;kSis the stiffness of the shaft, (xD,yD), (xB,yB), and (xO,yO) are the displacements of the disk, the journal of the deep-groove ball bearing, and the base in a fixed coordinate system, respectively;eμis the unbalanced mass eccentricity of the disk; (·) denotes d/dt;ωis the angular velocity of the rotor;mBis the centralized mass of the journal position and the inner ring position of the deep-groove ball bearings;kOis the support stiffness of the base;mOis the centralized mass of the base and the outer ring position of the deep-groove ball bearings;fBxandfByare the supporting forces of the deep-groove ball bearings in thexandydirections, respectively;FxandFyare the exciting forces of the exciter on the base in thexandydirections, respectively.

    (1)

    2.2 Dynamic model of ISFD supported rotor

    A Jeffcott rotor system model of the ISFD supported structure under base excitation is established, and its schematic diagram is shown in Fig.4. Both ends of the rotor adopt an ISFD support structure, which consists of a bearing seat, the ISFD, and deep-groove ball bearings. The excitation force is output by the exciter, providing base excitation for the rotor system. In Fig.4(b),fBrepresents the nonlinear supporting force of the deep-groove ball bearings.

    Fig.4 Dynamic model of a rotor supported on ball bearings with the ISFD subjected to base excitation

    The differential equations of motion of the system are as follows:

    (2)

    wheremIis the centralized mass of the outer ring of the deep-groove ball bearing and the ISFD inner ring,cIis the damping coefficient of the ISFD,kIis the stiffness of the ISFD, and (xI,yI) is the displacement of the ISFD inner ring and the deep-groove ball bearing outer ring in a fixed coordinate system.

    3 Introduction of experimental devices

    3.1 Support structure

    As mentioned above, two different supporting structures are used in the experiments: the rigid support and the ISFD support. Images of the two support types are shown in Fig.5. Fig.5(a) shows the rigid support, which is the traditional support configuration, consisting of a rigid sleeve, a deep-groove ball bearing, a bearing seat, an end cover, and a rubber O-ring. Fig.5(b) shows the ISFD support consisting of an integral squeeze film damper, a deep-groove ball bearing, a bearing seat, an end cover, and a rubber O-ring. In the ISFD support structure, a closed oil chamber is formed by the rubber O-ring, the end cover, and the bearing seat, which provides the ISFD with enough damping fluid to form a squeeze film.

    3.2 Test rig set-up

    Fig.6 shows the experimental platform with the single disk rotor system under base excitation, which includes a motor, a rotor, a disk, the ISFD support, the base, a base fixing bracket, an exciter, and other main components.

    A Jeffcott rotor is used in the test bench. The span of the shaft is 420 mm, the diameter of the shaft is 10 mm, the thickness of the disk is 15 mm, and the diameter of the disk is 78 mm. The two ends of the rotor are supported either by the rigid or the ISFD support structure, which consists of the bearing base, a rigid sleeve, the ISFD, and a deep-groove ball bearing.

    The exciter, which could exert a 0-100 N, 0-100 Hz sinusoidal base excitation force on the rotor system, is attached to the base fixing bracket through bolts. During the experiments, the LC8008 multi-channel fault diagnosis system for mechanical equipment is used to collect vibration displacement signals. The relative displacement of the rotating shaft near the disk and the base is measured by an eddy current sensor. Eddy current sensor 1 captures the relative displacement in the vertical direction (Y-direction) and eddy current sensor 2 captures the relative displacement in the horizontal direction (X-direction). The layout of the experimental set-up can be observed in Fig.6.

    1-Motor; 2-Coupling; 3-Photoelectric sensor; 4-ISFD support; 5-Shaft;6-Disk; 7-Eddy current sensor 1; 8-Eddy current sensor 2; 9-Exciter;10-Base fixing bracket; 11-Base

    4 Results and discussion

    4.1 Effect of different support structures on rotor system response under sinusoidal base excitation

    In order to investigate the effect of the different supports on the dynamic characteristics of the rotor, the rotor speed is set to 3 300 rpm, the excitation force applied by the exciter to the rotor system is 75 N, and the frequency is 65 Hz. In order to accurately capture the response characteristics, the displacement of the rotor relative to the base with the rigid or the ISFD support is recorded. Amplitude-time and amplitude-frequency comparison diagrams of the rotor vibration responses in theXandYdirections at the measuring points under different support structures are obtained and can be seen in Fig.7-Fig.10.

    Fig.7 Amplitude-time comparison diagram of the rotor in theX-direction under different supporting structures

    Fig.8 Amplitude-time comparison diagram of the rotor in theY-direction under different supporting structures

    Fig.9 Amplitude-frequency comparison diagram of the rotor in the X-direction under different supporting structures

    Fig.10 Amplitude-frequency comparison diagram of the rotor in the Y-direction under different supporting structures

    According to the amplitude-time comparison diagrams of Fig.7 and Fig.8, the vibration in the rotor system is effectively reduced when the ISFD support is used. In Fig.9 and Fig.10, it can be seen that the frequency spectrum mainly includes the rotational frequency of the rotor and the base excitation frequency, in which the amplitude corresponding to the rotational frequency is generated by the unbalanced vibration caused by the uneven mass distribution of the rotating disk. In addition, the amplitude corresponding to the base excitation frequency is generated by the excitation force produced by the exciter and transmitted to the rotor through the bearing seat.Through comparative analysis, it can be deduced that the ISFD support demonstrated a good vibration reduction effect on the base excitation frequency and the rotational frequency of the rotor with the reduction of vibration being more than 30%. More specifically, the amplitude of the rotor in theX-direction at a rotational frequency (55 Hz) decreased from 195.1 μm to 136.5 μm, and at an excitation frequency of 65 Hz, it decreased from 36.0 μm to 24.5 μm. The amplitude of the rotor in theY-direction at a rotational frequency of 55 Hz decreased from 181.6 μm to 127.9 μm, and at an excitation frequency of 65 Hz, it decreased from 22.6 μm to 13.9 μm.

    4.2 Effect of ISFD on rotor system response under sinusoidal excitation conditions with different frequency

    In order to investigate the effect of ISFD on rotor system response under sinusoidal excitation with different frequency, the speed of the rotor is set to 3 300 rpm and the exciter applies sinusoidal excitation of different frequency to the rotor system. The magnitude of the excitation force is 75 N and the frequencies are 0 Hz, 5 Hz, 10 Hz, …, 100 Hz. The response characteristics of the rotor system using the two different support structures is investigated. The vibration responses of the rotor in theXandYdirections under base excitations with different frequency are measured. In order to facilitate the observation of the vibration response trend, after the measured data are processed, spectrum waterfall diagrams are plotted (Fig.11-Fig.14).

    Fig.11 Spectral waterfall diagram of the rigid support rotor response in the X-direction under base excitations with different frequencies

    Fig.12 Spectral waterfall diagram of the rigid support rotor response in the Y-direction under base excitations with different frequencies

    Fig.13 Spectral waterfall diagram of the ISFD support rotor response in the X-direction under base excitations with different frequencies

    Fig.14 Spectral waterfall diagram of the ISFD support rotor response in the Y-direction under base excitations with different frequencies

    As it can be seen in Fig.11 and Fig.12, in the rotor system with rigid support, when the excitation frequency is close to 65 Hz or 90 Hz, there are 2 amplitude peaks in the spectrogram: the one corresponding to the rotor rotational frequency (55 Hz) and the other is the corresponding amplitude of the base excitation frequency. When the excitation frequency is not proximal to 65 Hz or 90 Hz, only the amplitude corresponding to the rotor rotational frequency (55 Hz) could be observed in the spectrogram, while the amplitude corresponding to the base excitation frequency is not present.

    From above-mentioned, under a constant excitation force condition, when the excitation frequency approached the frequency of the first critical speed (3 900 rpm) of the rotor system (65 Hz), the corresponding amplitude of the base excitation frequency in the spectrum diagram increased significantly, indicating that the base excitation can excite the natural frequency of the rotor system, which is very dangerous for the rotor system. When the excitation frequency approached 90 Hz, the amplitude corresponding to the base excitation frequency in the spectrum increased again, possibly due to that the frequency of the base excitation is close to the natural frequency of the rotor system base, which resulted in the resonance phenomenon. Due to that the stiffness of the base bracket is large, when the excitation frequency is far from 65 Hz or 90 Hz, the effect of the excitation force on the response of the rotor system is weak, since the amplitude corresponding to the base excitation frequency is not noticeable in the spectrum. As it can be seen in Fig.13 and Fig.14, similar to the rigid support rotor system, when the base excitation frequency of the ISFD-supported rotor system approached 65 Hz or 90 Hz, the corresponding amplitude of the base excitation frequency increased significantly, indicating that the base excitation of these 2 frequencies has the greatest impact on rotor system response.

    By comparing the responses of the rotor system, it is found that the vibration of the rotor system is significantly decreased after the ISFD support structure is adopted, and the corresponding amplitudes of the rotational and base excitation frequencies are reduced. In order to compare the response of the rotor system under the two different supporting structures, the amplitudes corresponding to the base excitation frequency and the rotational frequency with partial excitation frequencies in the spectrum waterfall diagram are identified and the results are presented in Table 2.

    Table 2 Amplitude comparison between the response of the rigid and the ISFD support rotor system under base excitation conditions with different frequencies

    As it can be seen in Table 2, the maximum amplitude reduction in the rotor rotational frequency, after the ISFD support structure is installed, is 30.6%, indicating that the ISFD can effectively reduce the vibration of a rotor system caused by unbalanced forces. The maximum amplitude reduction in the base excitation frequency with partial excitation frequency, after the ISFD support structure is installed, is about 40%, suggesting that the ISFD support structure can effectively suppress the vibration of the rotor system caused by base excitation. At an individual excitation frequency (such as 80 Hz), the vibration of the rotor is increased, while the vibration amplitude is maintained at a smaller value. To sum up, it can be seen that under base excitations of different frequencies, when the rigid support is replaced by the ISFD support, the vibration of the rotor system is significantly improved and the corresponding amplitudes of the rotational and base excitation frequencies are reduced. This indicates that the ISFD can effectively improve the vibration of rotor systems caused by unbalanced forces and different base excitation frequencies to a certain extent.

    4.3 Effect of ISFD on rotor system response under sinusoidal base excitation conditions with different excitation forces

    According to the spectrum waterfall diagrams in Fig.11 and Fig.12, when the excitation frequency is 65 Hz and 90 Hz, the base excitation has a great effect on the rotor system. In the following experiment, a sinusoidal base excitation of 65 Hz is applied to the rotor system in order to investigate the effect of ISFD on the response of the rotor system under sinusoidal base excitation conditions with different excitation forces. The rotor speed is again set to 3 300 rpm, while the exciter applies sinusoidal excitation with different excitation forces to the rotor system. The magnitude of the excitation frequency is 65 Hz and the applied forces are 0 N, 15 N, 30 N, …, 90 N. In order to observe the variation trend of the vibration response, after the measured data are processed, spectrum waterfall diagrams are plotted (Fig.15-Fig.18).

    Fig.15 Spectral waterfall diagram of the rigid support rotor response in the X-direction under base excitations with different forces

    Fig.16 Spectral waterfall diagram of the rigid support rotor response in the Y-direction under base excitations with different forces

    Fig.17 Spectral waterfall diagram of the ISFD support rotor response in the X-direction under base excitations with different forces

    Fig.18 Spectral waterfall diagram of the ISFD support rotor response in the Y-direction under base excitations with different forces

    As it can be seen in Fig.15-Fig.18, in the rigid support rotor system, as the excitation force increased, the amplitude corresponding to the base excitation frequency kept increasing. However, in the ISFD support rotor system, this trend is not observed due to the damping effect of the squeeze oil film.

    By comparing the response between the rigid support and the ISFD support rotor system, it can be seen that the vibration of the rotor system is effectively reduced after the ISFD support structure is adopted,and the corresponding rotation frequency and base excitation frequency amplitudes are reduced. In order to compare the changes in rotor system response under the 2 configurations, the amplitudes corresponding to base excitation and rotational frequencies under different base excitation forces are identified and the results are listed in Table 3.

    Table 3 Amplitude comparison between the response of the rigid and the ISFD support rotor system under base excitation conditions with different forces

    As it can be seen in Table 3, compared to the rotor system with the rigid support structure, the ISFD support structure could effectively suppress the rotor vibration caused by unbalanced forces, and the reduction in vibration amplitude reaches 32%. At the same time, the corresponding amplitude of the rotor excitation frequency is reduced, and the reduction reaches more than 40%. Based on the above analysis, under the action of different base excitation forces, when the rigid support is replaced by the ISFD support, the vibration of the rotor system is significantly improved, and the corresponding amplitudes of the rotational and base excitation frequencies are reduced. This indicates that the ISFD can effectively improve the vibration of rotor systems caused by unbalanced forces and base excitation of different excitation force magnitude to a certain extent.

    5 Conclusions

    The dynamic characteristics of a single disk rotor-rolling bearing-ISFD system under base excitation are investigated. The response of the rotor system under sinusoidal base excitation with different frequencies and excitation forces is experimentally studied. The effect of the ISFD on the dynamic characteristics of the rotor is investigated. Useful conclusions are drawn, which are as follows:

    (1) Comparative experiments demonstrate that when a certain sinusoidal base excitation is applied to a rotor system, the ISFD can significantly inhibit the vibration related to unbalanced forces and sinusoidal base excitation.

    (2) In the case where sinusoidal base excitation with different frequencies is applied to the rotor system, when the base excitation frequency approached the first critical speed of the rotor system (65 Hz) or the natural frequency of the rotor system (90 Hz), strong vibration in the rotor system appeares, indicating that the base excitation of these 2 frequencies has a greater effect on rotor system response. It indicates that the ISFD can effectively improve the vibration of rotor systems caused by unbalanced forces and different base excitation frequencies to a certain extent.

    (3) In the case where sinusoidal base excitation with different exciting forces is applied to the rotor system, when the base excitation force is increased, the vibration response of the rotor system becomes larger and larger. In addition, it is found that the ISFD has a good vibration reduction effect on sinusoidal base excitation with different exciting forces.

    18禁黄网站禁片免费观看直播| 丰满人妻一区二区三区视频av | 国产一级毛片七仙女欲春2| 免费搜索国产男女视频| 午夜福利欧美成人| 国产精品电影一区二区三区| 女警被强在线播放| 一个人免费在线观看的高清视频| 18美女黄网站色大片免费观看| 老汉色av国产亚洲站长工具| 巨乳人妻的诱惑在线观看| 老司机靠b影院| 黄色成人免费大全| 岛国在线观看网站| 天堂av国产一区二区熟女人妻 | 亚洲专区字幕在线| 欧美高清成人免费视频www| 人妻丰满熟妇av一区二区三区| 伦理电影免费视频| 成熟少妇高潮喷水视频| 国产视频内射| 午夜老司机福利片| 一级片免费观看大全| 久久久精品欧美日韩精品| 色综合站精品国产| 欧美黄色片欧美黄色片| 村上凉子中文字幕在线| 91字幕亚洲| 亚洲 欧美 日韩 在线 免费| 亚洲一区中文字幕在线| 高潮久久久久久久久久久不卡| 一本精品99久久精品77| 国产单亲对白刺激| 一级作爱视频免费观看| 欧美zozozo另类| 亚洲七黄色美女视频| 久久人妻福利社区极品人妻图片| 欧美一区二区国产精品久久精品 | 1024视频免费在线观看| 色精品久久人妻99蜜桃| cao死你这个sao货| 亚洲狠狠婷婷综合久久图片| 国产午夜精品久久久久久| 免费高清视频大片| 日韩三级视频一区二区三区| 国产私拍福利视频在线观看| 草草在线视频免费看| 久久精品夜夜夜夜夜久久蜜豆 | 国产黄a三级三级三级人| 欧美午夜高清在线| 黄色女人牲交| 一区福利在线观看| 一区福利在线观看| 久久亚洲精品不卡| 欧美黑人精品巨大| 在线永久观看黄色视频| 亚洲成人国产一区在线观看| 2021天堂中文幕一二区在线观| 好男人电影高清在线观看| 大型黄色视频在线免费观看| 天天躁夜夜躁狠狠躁躁| 波多野结衣高清无吗| 国产99白浆流出| 国产精品久久久久久久电影 | 免费观看人在逋| 亚洲九九香蕉| 精品国产乱码久久久久久男人| 亚洲人成77777在线视频| 欧美一区二区国产精品久久精品 | 亚洲成av人片免费观看| 久久久久九九精品影院| 高清毛片免费观看视频网站| 这个男人来自地球电影免费观看| 国产精品久久视频播放| 中亚洲国语对白在线视频| 麻豆av在线久日| 久久九九热精品免费| 亚洲人与动物交配视频| 久久这里只有精品19| aaaaa片日本免费| 欧美黄色淫秽网站| 成人国语在线视频| 亚洲欧美精品综合久久99| 久久久久性生活片| www国产在线视频色| 国产区一区二久久| 久久久久国产一级毛片高清牌| 女人被狂操c到高潮| 中出人妻视频一区二区| 又爽又黄无遮挡网站| 女人被狂操c到高潮| 精品欧美一区二区三区在线| 性欧美人与动物交配| 欧美中文综合在线视频| 校园春色视频在线观看| 午夜成年电影在线免费观看| 特大巨黑吊av在线直播| svipshipincom国产片| 国产久久久一区二区三区| 欧美又色又爽又黄视频| 天堂av国产一区二区熟女人妻 | 久久精品91蜜桃| 999久久久国产精品视频| 日韩欧美国产一区二区入口| 午夜免费成人在线视频| 欧美日本亚洲视频在线播放| 精华霜和精华液先用哪个| 亚洲成人中文字幕在线播放| 国产免费av片在线观看野外av| 两个人免费观看高清视频| 亚洲精品色激情综合| 岛国视频午夜一区免费看| 丰满人妻一区二区三区视频av | 国产一区二区三区在线臀色熟女| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av日韩精品久久久久久密| 露出奶头的视频| 久久中文字幕一级| 人成视频在线观看免费观看| 一边摸一边抽搐一进一小说| 亚洲精品久久国产高清桃花| 黄片大片在线免费观看| 久9热在线精品视频| 一进一出抽搐动态| 午夜福利成人在线免费观看| 美女 人体艺术 gogo| 国内精品久久久久久久电影| 无人区码免费观看不卡| 很黄的视频免费| 91麻豆av在线| 欧美极品一区二区三区四区| 美女扒开内裤让男人捅视频| 丰满人妻熟妇乱又伦精品不卡| 一级毛片精品| 亚洲av片天天在线观看| 亚洲欧美日韩东京热| 亚洲最大成人中文| 亚洲 欧美 日韩 在线 免费| 久久久久久人人人人人| 欧美又色又爽又黄视频| 亚洲国产精品sss在线观看| 久久精品aⅴ一区二区三区四区| 国产主播在线观看一区二区| 精品乱码久久久久久99久播| 天堂√8在线中文| 亚洲精品美女久久av网站| 黄色视频不卡| 制服丝袜大香蕉在线| 日韩av在线大香蕉| 成人手机av| 免费在线观看亚洲国产| 色综合欧美亚洲国产小说| 欧美久久黑人一区二区| 亚洲午夜理论影院| 成人国产综合亚洲| 亚洲熟女毛片儿| 亚洲一区二区三区色噜噜| 国产av在哪里看| 国内揄拍国产精品人妻在线| 深夜精品福利| 国产亚洲欧美在线一区二区| 国产精品精品国产色婷婷| 亚洲精品久久国产高清桃花| 国产麻豆成人av免费视频| 蜜桃久久精品国产亚洲av| 正在播放国产对白刺激| 精品国产乱子伦一区二区三区| 国产成人精品久久二区二区免费| 国产精品1区2区在线观看.| 久久婷婷成人综合色麻豆| 国产三级中文精品| 亚洲精品一区av在线观看| 欧美黑人精品巨大| 国产精品九九99| 一个人免费在线观看电影 | 五月玫瑰六月丁香| 久久中文字幕人妻熟女| 欧美在线一区亚洲| 嫁个100分男人电影在线观看| 欧美日韩一级在线毛片| 久久精品国产99精品国产亚洲性色| 日本在线视频免费播放| 久久 成人 亚洲| 欧美精品啪啪一区二区三区| 亚洲 欧美 日韩 在线 免费| 2021天堂中文幕一二区在线观| av欧美777| 久久中文看片网| 久久久久久国产a免费观看| 淫秽高清视频在线观看| 老鸭窝网址在线观看| 国产精品美女特级片免费视频播放器 | 欧美日韩瑟瑟在线播放| 在线天堂最新版资源| 蜜臀久久99精品久久宅男| 草草在线视频免费看| 亚洲av成人av| 一区二区三区四区激情视频 | 黄色一级大片看看| 婷婷精品国产亚洲av| 成人亚洲欧美一区二区av| 99国产极品粉嫩在线观看| 免费av毛片视频| 久久久久久大精品| 校园人妻丝袜中文字幕| 日韩高清综合在线| 欧美色欧美亚洲另类二区| av黄色大香蕉| 免费看光身美女| 国内精品美女久久久久久| 精品国产三级普通话版| 中文字幕熟女人妻在线| 日本欧美国产在线视频| 成熟少妇高潮喷水视频| 人人妻人人看人人澡| 国产精品久久久久久久电影| 五月伊人婷婷丁香| 日韩三级伦理在线观看| 搡老妇女老女人老熟妇| 久久国内精品自在自线图片| 亚洲欧美精品综合久久99| 国产av麻豆久久久久久久| 亚洲av二区三区四区| 国产亚洲5aaaaa淫片| АⅤ资源中文在线天堂| 欧美zozozo另类| 又爽又黄无遮挡网站| 啦啦啦韩国在线观看视频| 免费无遮挡裸体视频| 亚洲成人精品中文字幕电影| 亚洲精品色激情综合| 国产精品女同一区二区软件| 亚洲最大成人手机在线| 免费看日本二区| 校园春色视频在线观看| 18+在线观看网站| 日本-黄色视频高清免费观看| 最新中文字幕久久久久| 校园人妻丝袜中文字幕| 午夜激情欧美在线| 国产 一区 欧美 日韩| 成人高潮视频无遮挡免费网站| 男的添女的下面高潮视频| 国产精品一及| 少妇裸体淫交视频免费看高清| 亚洲四区av| 久久精品国产鲁丝片午夜精品| www.av在线官网国产| 亚洲av中文av极速乱| 天堂√8在线中文| 禁无遮挡网站| 高清在线视频一区二区三区 | 成熟少妇高潮喷水视频| 国产伦精品一区二区三区四那| 99久久九九国产精品国产免费| 只有这里有精品99| 丰满乱子伦码专区| 麻豆国产av国片精品| 真实男女啪啪啪动态图| 看十八女毛片水多多多| 在线观看一区二区三区| 一级二级三级毛片免费看| 熟女人妻精品中文字幕| 一本精品99久久精品77| 白带黄色成豆腐渣| 看非洲黑人一级黄片| 69人妻影院| 女人被狂操c到高潮| 亚洲无线观看免费| 1000部很黄的大片| 五月玫瑰六月丁香| 成人无遮挡网站| 直男gayav资源| 欧美人与善性xxx| 欧美又色又爽又黄视频| 久久精品综合一区二区三区| 日本一二三区视频观看| 在线播放国产精品三级| 成人永久免费在线观看视频| 免费看a级黄色片| 九九热线精品视视频播放| 成人性生交大片免费视频hd| 亚洲欧美日韩高清在线视频| .国产精品久久| avwww免费| 久久6这里有精品| 精品无人区乱码1区二区| 国产午夜精品一二区理论片| 一区二区三区高清视频在线| 日韩三级伦理在线观看| 春色校园在线视频观看| 精品人妻一区二区三区麻豆| 黄色一级大片看看| 欧美日韩精品成人综合77777| 日本免费a在线| 青春草国产在线视频 | 悠悠久久av| 五月玫瑰六月丁香| 美女cb高潮喷水在线观看| 亚洲第一电影网av| 九九在线视频观看精品| 国产精品国产三级国产av玫瑰| 久久久久九九精品影院| 亚洲第一电影网av| .国产精品久久| 亚洲欧洲国产日韩| 91在线精品国自产拍蜜月| 亚洲精品456在线播放app| 看黄色毛片网站| 村上凉子中文字幕在线| 一级二级三级毛片免费看| 国产中年淑女户外野战色| 欧美变态另类bdsm刘玥| 精品人妻偷拍中文字幕| 变态另类丝袜制服| 春色校园在线视频观看| 男女那种视频在线观看| 不卡视频在线观看欧美| 国产高清不卡午夜福利| 中文字幕av成人在线电影| 别揉我奶头 嗯啊视频| 国产精品99久久久久久久久| 久久人人精品亚洲av| 91在线精品国自产拍蜜月| 欧美极品一区二区三区四区| 成年版毛片免费区| 国产精品免费一区二区三区在线| 欧美xxxx黑人xx丫x性爽| av在线老鸭窝| 亚洲第一区二区三区不卡| 国产精品一区二区三区四区久久| 性色avwww在线观看| 丝袜喷水一区| 永久网站在线| 中文精品一卡2卡3卡4更新| 成人av在线播放网站| 99久久精品一区二区三区| 男女啪啪激烈高潮av片| 五月玫瑰六月丁香| 国产老妇伦熟女老妇高清| 亚洲人成网站在线观看播放| 人人妻人人澡欧美一区二区| 观看美女的网站| 国产久久久一区二区三区| 亚洲人成网站在线播放欧美日韩| 国产蜜桃级精品一区二区三区| 日韩高清综合在线| 亚洲在线自拍视频| 亚洲精品自拍成人| 偷拍熟女少妇极品色| 亚洲av电影不卡..在线观看| 三级男女做爰猛烈吃奶摸视频| 国产精品,欧美在线| 麻豆国产av国片精品| a级一级毛片免费在线观看| 国产精品无大码| 精品久久久久久久末码| 国产高清三级在线| 亚洲精品久久久久久婷婷小说 | 97超碰精品成人国产| 91aial.com中文字幕在线观看| 超碰av人人做人人爽久久| 欧美另类亚洲清纯唯美| 九草在线视频观看| 日日撸夜夜添| 国内精品宾馆在线| 成人午夜精彩视频在线观看| 91aial.com中文字幕在线观看| 69av精品久久久久久| 蜜桃久久精品国产亚洲av| 午夜福利在线观看免费完整高清在 | 日本爱情动作片www.在线观看| 男人狂女人下面高潮的视频| 两性午夜刺激爽爽歪歪视频在线观看| 色综合亚洲欧美另类图片| 久久国内精品自在自线图片| 中文资源天堂在线| 久久精品国产亚洲av香蕉五月| 最近手机中文字幕大全| 免费观看精品视频网站| ponron亚洲| 亚洲av一区综合| 天堂影院成人在线观看| 久久精品91蜜桃| 欧美一区二区国产精品久久精品| 国产亚洲av嫩草精品影院| 欧美成人精品欧美一级黄| 成人一区二区视频在线观看| 最好的美女福利视频网| 青青草视频在线视频观看| 久久久a久久爽久久v久久| 一进一出抽搐动态| 高清毛片免费看| 精华霜和精华液先用哪个| 日本撒尿小便嘘嘘汇集6| 亚洲电影在线观看av| 性插视频无遮挡在线免费观看| 欧美高清成人免费视频www| 亚洲内射少妇av| 老师上课跳d突然被开到最大视频| 国产av在哪里看| 一本一本综合久久| 人妻制服诱惑在线中文字幕| 青春草亚洲视频在线观看| 欧美潮喷喷水| 国产精品免费一区二区三区在线| 99久国产av精品| 国产精品综合久久久久久久免费| 美女cb高潮喷水在线观看| h日本视频在线播放| 青青草视频在线视频观看| 看黄色毛片网站| 亚洲欧洲日产国产| 国产高清三级在线| 少妇丰满av| 国产精品伦人一区二区| 最好的美女福利视频网| 国产淫片久久久久久久久| 九九在线视频观看精品| 成人午夜高清在线视频| 国产伦理片在线播放av一区 | 男女做爰动态图高潮gif福利片| 国产精品一区www在线观看| 一个人看的www免费观看视频| 69人妻影院| 人人妻人人看人人澡| 97超视频在线观看视频| 亚洲第一电影网av| 哪个播放器可以免费观看大片| 国产精品美女特级片免费视频播放器| 在线天堂最新版资源| 精品免费久久久久久久清纯| 18+在线观看网站| 少妇的逼好多水| 天堂影院成人在线观看| 国产午夜福利久久久久久| 亚洲五月天丁香| 日韩国内少妇激情av| 精华霜和精华液先用哪个| 久久鲁丝午夜福利片| 日本黄大片高清| 村上凉子中文字幕在线| 99热网站在线观看| 欧美性感艳星| 国产老妇伦熟女老妇高清| 日韩一本色道免费dvd| 久久精品国产亚洲av香蕉五月| 成人性生交大片免费视频hd| 精品久久久久久久久久免费视频| 欧美区成人在线视频| 久久精品影院6| 国产老妇伦熟女老妇高清| 麻豆成人av视频| 97超视频在线观看视频| 免费看日本二区| 亚洲成人久久性| 日日干狠狠操夜夜爽| 可以在线观看毛片的网站| 色噜噜av男人的天堂激情| 亚洲国产色片| 婷婷精品国产亚洲av| 日韩强制内射视频| 色尼玛亚洲综合影院| 22中文网久久字幕| 欧美zozozo另类| 久久这里有精品视频免费| 少妇丰满av| 久久99热6这里只有精品| 熟女电影av网| 嫩草影院精品99| 天堂影院成人在线观看| 亚洲欧美精品专区久久| 国模一区二区三区四区视频| 国内精品久久久久精免费| 亚洲欧美日韩无卡精品| 欧美色视频一区免费| 色综合色国产| 国产精品爽爽va在线观看网站| 久久久a久久爽久久v久久| 日韩欧美精品v在线| 亚洲国产色片| 三级男女做爰猛烈吃奶摸视频| 日韩成人av中文字幕在线观看| 国产在线精品亚洲第一网站| 美女大奶头视频| 国产黄色视频一区二区在线观看 | 天堂网av新在线| 久久亚洲精品不卡| 色视频www国产| 日韩一本色道免费dvd| 国产亚洲av嫩草精品影院| 波野结衣二区三区在线| 日本爱情动作片www.在线观看| 成人欧美大片| 久久这里只有精品中国| 成人鲁丝片一二三区免费| 日韩欧美 国产精品| 高清毛片免费观看视频网站| 欧美人与善性xxx| 久久精品91蜜桃| 国产成人a∨麻豆精品| 国产精品99久久久久久久久| 久久久久久国产a免费观看| 免费观看精品视频网站| 亚洲国产欧美在线一区| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美精品v在线| 非洲黑人性xxxx精品又粗又长| 亚洲一区二区三区色噜噜| 一区二区三区免费毛片| 老熟妇乱子伦视频在线观看| 最近视频中文字幕2019在线8| 麻豆一二三区av精品| 国产亚洲91精品色在线| 亚洲人成网站高清观看| 亚洲美女视频黄频| 国产精品久久久久久久久免| 自拍偷自拍亚洲精品老妇| 三级毛片av免费| 一进一出抽搐gif免费好疼| 美女cb高潮喷水在线观看| 婷婷精品国产亚洲av| 亚洲精品影视一区二区三区av| 国产伦一二天堂av在线观看| 中文字幕熟女人妻在线| 国产伦精品一区二区三区四那| 黑人高潮一二区| 观看免费一级毛片| 国产精品一及| 一本精品99久久精品77| avwww免费| 国产日本99.免费观看| 此物有八面人人有两片| 99在线视频只有这里精品首页| 国语自产精品视频在线第100页| 一个人免费在线观看电影| 国产伦在线观看视频一区| 亚洲最大成人av| 午夜免费男女啪啪视频观看| 久久久久久大精品| 天堂av国产一区二区熟女人妻| 国产高清视频在线观看网站| 日本一本二区三区精品| 美女黄网站色视频| 天堂影院成人在线观看| 在线免费观看的www视频| 国产乱人偷精品视频| 精品一区二区免费观看| 国内久久婷婷六月综合欲色啪| 男人舔女人下体高潮全视频| 一级av片app| 亚洲成人中文字幕在线播放| 麻豆精品久久久久久蜜桃| 欧美又色又爽又黄视频| 久久人人爽人人片av| 亚洲精品乱码久久久久久按摩| 国产美女午夜福利| 日韩欧美 国产精品| 熟女人妻精品中文字幕| 日韩一本色道免费dvd| 在线免费十八禁| 久久热精品热| 久久人人爽人人片av| 午夜福利在线在线| 日本黄色片子视频| 亚洲乱码一区二区免费版| 国产亚洲欧美98| 精品人妻视频免费看| 久久99热6这里只有精品| 人人妻人人澡欧美一区二区| 韩国av在线不卡| 自拍偷自拍亚洲精品老妇| 久久久国产成人精品二区| 亚洲国产精品sss在线观看| 免费看av在线观看网站| 欧美日韩在线观看h| 久久草成人影院| 成人国产麻豆网| 国产精品永久免费网站| 高清日韩中文字幕在线| 国产 一区 欧美 日韩| 黑人高潮一二区| 联通29元200g的流量卡| 男女边吃奶边做爰视频| 国产在线男女| 超碰av人人做人人爽久久| 国产探花极品一区二区| 亚洲欧美日韩卡通动漫| 黄色日韩在线| 美女cb高潮喷水在线观看| 1000部很黄的大片| 国产精品久久久久久精品电影小说 | 九九热线精品视视频播放| 亚洲无线观看免费| 亚洲欧美日韩卡通动漫| 欧美激情久久久久久爽电影| 六月丁香七月| 欧美日韩一区二区视频在线观看视频在线 | 国内少妇人妻偷人精品xxx网站| 国产人妻一区二区三区在| 中文字幕av成人在线电影| 欧美色欧美亚洲另类二区| 亚洲成人精品中文字幕电影| 特大巨黑吊av在线直播| 久久久成人免费电影| 中文在线观看免费www的网站| 国产男人的电影天堂91| 欧美最黄视频在线播放免费| 97在线视频观看| 搞女人的毛片| av卡一久久| 午夜福利在线观看免费完整高清在 | 欧美高清性xxxxhd video| 国产精品,欧美在线| 尤物成人国产欧美一区二区三区|