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

    Characteristic Analysis and Harmonic Feature Identification of Micro-Vibration on Flywheels

    2021-04-08 11:08:52YINXianbo尹顯波SHENGXiaowei盛曉偉XUYangSHENYan

    YIN Xianbo(尹顯波), SHENG Xiaowei(盛曉偉), XU Yang(徐 洋), SHEN Yan(申 妍)

    College of Mechanical Engineering, Donghua University, Shanghai 201620, China

    Abstract: To avoid the negative effects of disturbances on satellites, the characteristics of micro-vibration on flywheels are studied. Considering rotor imbalance, bearing imperfections and structural elasticity, the extended model of micro-vibration is established. In the feature extraction of micro-vibration, singular value decomposition combined with the improved Akaike Information Criterion (AIC-SVD) is applied to denoise. More robust and self-adaptable than the peak threshold denoising, AIC-SVD can effectively remove the noise components. Subsequently, the effective harmonic coefficients are extracted by the binning algorithm. The results show that the harmonic coefficients have great identification in frequency domain. Except for the fundamental frequency caused by rotor imbalance, the harmonics are also caused by the coupling of imperfections on bearing components.

    Key words: flywheel; micro-vibration; sigular value decomposition comined with the improved Akaike Information Criterion (AIC-SVD); harmonic coefficient; binning algorithm

    Introduction

    In addition to the space environment, satellites are also disturbed by the internal airborne on orbit[1]. The internal disturbance sources of satellites mainly include flywheel system, solar cell array drive mechanism, cryogenic cooler, camera swing mirror mechanism and thruster,etc. Among them, the micro-vibration of flywheel system has the greatest interference to satellites[2]. Micro-vibration refers to the low-intensity mechanical vibration or disturbance of spacecraft in microgravity environment. It usually occurs in the range of 1 Hz to 1 000 Hz. With a small amplitude, the bandwidth and long duration, micro-vibration is difficult to be suppressed or compensated by the attitude adjustment mechanism[3]. Through the cabin platform, it decreases the pointing accuracy and attitude stability of satellites[4].

    To reduce the adverse effects of micro-vibration, the characteristics of micro-vibration on flywheels are studied and analyzed. Through the experiment method, Mastersonetal.[5]derived the steady-state empirical model of micro-vibration. Considering the resonance, Zhaoetal.[6]introduced the amplification factor to optimize the steady-state empirical model. On this basis, they carried out disturbance comparisons of flywheels under rigid and flexible installation conditions. As for the characteristic analysis of micro-vibration, the binning algorithm is widely used to extract feature parameters[5-7]. Unfortunately, there are still problems such as poor adaptability and large deviation in signal processing. Based on the vibration test of flywheels, Sunetal.[8]improved the identification accuracy of the amplitude coefficient by the energy compensation method. Mbaetal.[9]used stochastic resonance technology to extract the weak pulse characteristics of micro-vibration, which achieved the purpose of fault diagnosis.

    To further discuss the micro-vibration characteristics of flywheels, the dynamic extended model was established based on the micro-vibration mechanism. After preprocessing by singular value decomposition combined with the improved Akaike Information Criterion (AIC-SVD), the binning algorithm achieves more accurate harmonic extraction. Furthermore, the corresponding harmonic disturbance source can be easily identified through the feature parameters.

    1 Micro-vibration Mechanism

    1.1 Flywheel structure

    As shown in Fig. 1, the flywheel is composed of a brushless direct current (DC) motor, a support shaft, a rotor, bearings, a housing and a base. The magnetic torque generated by the coil stator drives the rotor. Through the mounting seat, the base is fixed on the satellite platform. According to the law of conservation of angular momentum, the flywheel controls satellite attitude and balances space disturbances[10].

    Fig. 1 Structure of flywheels

    1.2 Micro-vibration source

    The micro-vibration of flywheels originates from active disturbances, structural disturbances and other disturbances. Active disturbances are mainly caused by rotor imbalance, bearing imperfections and electromagnetic torque fluctuation. Structural disturbances are caused by elasticity of rotor or shell. Meanwhile, flywheels are also affected by friction, lubrication or imperfect factors[11].

    Due to the influence of structural design and material inequality, there are some errors between the rotating shaft and the inertia axis. The errors result in static imbalance and dynamic imbalance of the rotor. As shown in Fig. 2, the equivalent massmsandmdare used to simulate the static imbalance and the dynamic imbalance[12].

    Fig. 2 Illustration of rotor imbalance

    Through the support bearing, the disturbances generated by rotor imbalance produce the fundamental harmonic. In addition, bearing components also produce impact pulses because of their own imperfections during operation, as shown in Fig. 3. Substituting the bearing parameters into the frequency formulas[13], the relationship between the pulse frequencies and the rotational frequencies can be obtained. In Table 1,fois the rotation frequencies of the outer ring, andfiis the rotation frequencies of the inner ring.fIRI,fORI,fCRandfBIcorrespond to the pulse frequencies caused by the inner race, the outer race, the cage and the ball, respectively. In a single operating cycle, the imperfect position of the rolling element contacts with the inner and outer rings respectively. Therefore, the pulse frequencyfBIis twice as much as the rotation frequencyfBS.

    Fig. 3 Illustration of imperfections on bearing parts

    Table 1 Impulse frequency of imperfections on bearing parts

    2 Modeling of Micro-vibration

    Considering the imbalance of the rotor and the bearing imperfections, the micro-vibration of flywheels is assumed to be a series of discrete harmonics proportional to the square of the speed. The steady-state empirical model can capture the harmonic features of flywheels. Lacking the consideration of structural elasticity, it performs in bad agreement with high order harmonics. To reflect the elastic characteristics of the bearing support, the flywheel is modeled according to the dynamic theory. The angular contact ball bearing is equivalent to the stiffness-damper connected to the supporting shaft. Considering the rotor imbalance, the rotor-bearing model with five degrees of freedom is constructed, as shown in Fig. 4.

    Fig. 4 Five degrees of freedom model for rotor-bearing

    Lagrange’s equation is used to derive the system motion differential equation

    (1)

    Considering bearing imperfections, the five degrees of freedom model for rotor-bearing is improved by combining with the steady-state empirical model. The model is given as

    (2)

    (3)

    I always waited for Jennifer and Amy after classes, so we could walk together to our next class. Amy and Jennifer chatted by Jennifer s desk as Jennifer packed her books up and I waited by the door. Sometimes when they left, they d walk right past me. No Thanks for waiting. No Sorry we took so long. It was as if they couldn t even see me. Yet I still waited after every single class.

    According to Eqs. (2) and (3), the main modes of the rotor bearing system can be obtained. Ignoring the system damping, the radial translation frequencyfr, the axial translation frequencyfa, and the swing frequencyf1, 2are respectively

    (4)

    (5)

    (6)

    Substituting the parameters of Table 2 into Eqs. (4)-(6), the relationship between the structural modal frequencies and the rotating speed are obtained, as shown in Fig. 5. In Fig. 5,frandfaare 455 Hz and 197 Hz, respectively, which remain constant with the speed change. Due to the gyroscopic effect of the rotor, the swing mode can be divided into forward precessionf1and reverse precessionf2. Among them,f1increases with the rotating speed, butf2is inversely proportional to the rotating speed. Within the operating speed, the structure mode has no interaction with the fundamental frequency.

    Table 2 Structural parameters of the flywheel

    Fig. 5 Campbell diagram of micro-vibration on the flywheel

    3 Harmonic Feature Extraction

    3.1 Extraction method

    In the extended model of micro-vibration, the harmonic coefficients can directly reflect the components of harmonics. Using the binning algorithm to extract the features, the corresponding disturbance sources can be identified. As shown in Fig. 6, the vibration component signalsΩuunder the speedΩis collected. After noise reduction processing, the peak values (fpeak,Apeak) ofxΩuare extracted in the spectrum. In addition to harmonic disturbances,fpeakalso includes some structural modes. To extract the frequencies of harmonic features accurately, the modal frequencies solved by the extended model are judged to be removed. After normalizing harmonic frequencies, the toleranceεand the peak ratioPkare set to obtain effective harmonic coefficientshui.

    Fig. 6 Flow chart of binning algorithm

    In the binning algorithm, reasonable preprocessing can get double results with half the effort. As we all know, the acquired signals inevitably are mixed with interference noise. For flywheel vibration signals, the current noise reduction processing is mainly based on the peak threshold denoising[5]. Affected by uncertain parametersNσ, the method of peak threshold denoising is prone to over-denoising or under-denoising. As a result, the stability and accuracy also decrease in feature extraction. To extract the harmonic coefficients of vibration effectively, the denoising method named AIC-SVD is used to improve it.

    3.2 Denoising by AIC-SVD

    Based on the singular value decomposition (SVD) of Hankel matrix, the signal can be decomposed into a simple linear superposition of a series of component signals[14]. For additive noise, an advantage of this decomposition is that the main components can be retained directly to achieve the purpose of denoising. According to the principle of matrix orthogonal transformation, the Hankel matrixAcan be obtained by singular value decomposition[15].

    The singular valuesσicontain all the characteristic information of the signal. To effectively separate the main characteristic component, the orderkof the singular values is determined by the AIC. AIC is an estimated measure of the fitting goodness of statistical models. The function of AIC contains the terms of parameter fitness and inserted bias correction. When the sum of the two terms is the smallest, it means that the function balances the two factors to get the best estimate of the effective order[16].

    In addition to white noise with uniform power, there are also colored noises produced by friction or electromagnetic interference in micro-vibration[17]. To smooth the interference components in the background of colored noise, the diagonal loading technique[18]is used to correct the singular value

    (7)

    Substituting the modified singular values into the function of AIC, the improved decision functionAimpand the effective orderkare respectively

    2d(2q-d),

    (8)

    (9)

    whereNrepresents the number of signal sampling points,nis the row of the Hankel matrixA,d=1, 2, …,n-1 denotes the number of sources, andqis the number of non-zero singular values.

    The flow chart of AIC-SVD denoising method is shown in Fig. 7. Firstly, the Hankle matrix of the measured signal is constructed and solved by SVD. Then the improved AIC is performed to determine the effective order of the singular values. Finally, Through the inverse operation of SVD and the average method, the frontkorder components are reconstructed to obtain the denoised signal.

    Fig. 7 Flow chart of signal denoising method using AIC-SVD

    4 Analysis of Micro-vibration

    4.1 Test of micro-vibration

    The micro-vibration test of the flywheel is carried out based on six component force measurement method, as shown in Fig. 8(a). Through the transfer tooling, the disturbances of the flywheel are transmitted to the Kistler force platform. Meanwhile, the built-in sensors of the platform acquire signals in real time. Four piezoelectric three-dimensional force sensors are selected to collect voltage signals. There are evenly distributed between the top plate and the bottom plate of the platform, as shown in Fig. 8(b). Through the charge amplifier, the vibration signals collected by the sensors are transmitted to the data acquisition system. Finally, the tested signals are processed and displayed on the computer.

    Fig. 8 Micro-vibration test of the flywheel

    The disturbance force of environment is less than 0.2 N, and disturbance torque is less than 0.01 N·m. The test result meets the requirements of the force measuring system for the environment. Within the working speed range, the acquisition is performed every 200 r/min interval. The spectrum waterfall diagrams of acquired signals are shown in Fig. 9. It can be seen that the axial disturbance momentMzis very small and could be neglected. Vibration mainly occurs in the radial directionFx, the axial directionFzand the radial directionMx.

    Fig. 9 Waterfall diagrams of the flywheel

    4.2 Harmonic identification

    TakingFxas an example to carry out the analysis of noise reduction, the peak threshold values (DT) are calculated as shown in Fig. 10. In Fig. 10(a), the threshold values increase with the speed. However, there are sharp change occurs at 400 r/min and 1 800 r/min. Analyzing Fig. 10(b), it can be seen that the threshold values increase due to the shell mode of 131.5 Hz at 1 800 r/min. This leads to the loss of some important features, such as the 18 Hz harmonic caused by the bearing cage.

    Fig. 10 Peak threshold denoising in Fx

    To eliminate the adverse effects of noise, the vibration signals of the flywheel are processed by AIC-SVD and shown in Fig. 11. AIC-SVD can effectively retain the weak features in the signal. Among them, 435 Hz is the radial translation frequencyfrof rotor-bearing. The relative error between theory and experiment is 4.40%, which is mainly caused by structural damping. In addition, the axial translational frequency of rotor bearing is 195.5 Hz, and the relative error is 0.76%.

    Fig. 11 AIC-SVD denoising in Fx

    Furthermore, the energy of the denoised signal is analyzed. The energy ratio of the main component to the total signal is calculated, as shown in Fig. 12. After AIC-SVD denoising, the main feature energy ratio is 92.47%, and the corresponding root mean square error (RMSE)Rmis 0.022. Through comparison, the RMSE of AIC-SVD is far less than 0.074 and 0.124 from the peak threshold denoising. It shows that AIC-SVD has better adaptability and stability compared with the traditional method.

    Synthesizing the vibration signals under all measured speeds, the structural modes are removed from the system. Based on AIC-SVD denoising and the binning algorithm, the harmonic coefficients of each direction are extracted. The obtained parameters are shown in Table 3. Combined with the vibration mechanism, it can be known that the denoised signals also contain frequencies offCR,fo,fBI,fORIand multi-harmonic coupling. In the directionxandθ, the flywheel is disturbed by 0.6 times harmonic caused by the bearing cage imperfections and fundamental harmonic caused by the rotor imbalance. The harmonic coefficients of 4.39, 4.98 and 7.13 are caused by the imperfections of rolling element, outer ring and inner ring respectively. Except structural modal factors, the micro-vibration in the directionzis mainly caused by the bearing imperfections.

    Fig. 12 Energy ratio of main features in Fx

    Table 3 Main harmonic coefficients of micro-vibration

    5 Conclusions

    (1) The harmonic disturbance, structural mode and coupling resonance are the main factors of micro-vibration on flywheels, accounting for 92.47% of the total energy. Therefore, the extended model of micro-vibration combining harmonic features and structural elasticity can better express the micro-vibration characteristics.

    (2) The AIC-SVD denoising method is applied to the signal preprocessing of micro-vibration, which can effectively reduce the noise. Compared with the traditional peak threshold denoising, it has better stability and adaptability.

    (3) Processed by AIC-SVD, the principal components of the micro-vibration signals have great identification in frequency domain. The analysis shows that the harmonics along radial direction and around radial direction are caused by rotor imbalance and bearing imperfections. And the harmonics along radial direction is mainly caused by bearing imperfections.

    成人亚洲精品一区在线观看| aaaaa片日本免费| 老司机午夜福利在线观看视频| 国产精品久久久久久人妻精品电影| 亚洲精品中文字幕一二三四区| 大陆偷拍与自拍| 久久精品熟女亚洲av麻豆精品| 少妇粗大呻吟视频| 欧美日韩av久久| 精品久久久久久久久久免费视频 | 国产精品一区二区在线观看99| 欧美日韩亚洲综合一区二区三区_| 色94色欧美一区二区| 午夜福利免费观看在线| 久久精品国产亚洲av香蕉五月 | 91麻豆av在线| av网站在线播放免费| 久久亚洲真实| 欧美在线一区亚洲| 热99久久久久精品小说推荐| 丰满的人妻完整版| 午夜福利视频在线观看免费| 欧美av亚洲av综合av国产av| 久久久久久免费高清国产稀缺| 亚洲国产欧美日韩在线播放| 又黄又粗又硬又大视频| 欧美成人午夜精品| 亚洲欧洲精品一区二区精品久久久| 高潮久久久久久久久久久不卡| 热99re8久久精品国产| 青草久久国产| 国产av又大| 欧美黑人精品巨大| av不卡在线播放| 精品福利观看| 亚洲欧美一区二区三区久久| xxxhd国产人妻xxx| 久久久水蜜桃国产精品网| 亚洲综合色网址| 一级a爱视频在线免费观看| av欧美777| 久久九九热精品免费| 99国产精品99久久久久| 国产一区二区三区视频了| 91成年电影在线观看| 91成人精品电影| 亚洲一区中文字幕在线| 老熟妇乱子伦视频在线观看| 麻豆乱淫一区二区| 热re99久久精品国产66热6| 丰满饥渴人妻一区二区三| 国产成人欧美| 亚洲 国产 在线| 999久久久精品免费观看国产| 妹子高潮喷水视频| 中文字幕制服av| 每晚都被弄得嗷嗷叫到高潮| 欧美精品亚洲一区二区| 亚洲午夜精品一区,二区,三区| 在线国产一区二区在线| 成人18禁在线播放| 18禁裸乳无遮挡动漫免费视频| 久久这里只有精品19| 夜夜夜夜夜久久久久| 69精品国产乱码久久久| 满18在线观看网站| 亚洲av成人av| 女性生殖器流出的白浆| 亚洲性夜色夜夜综合| 欧美日韩乱码在线| 女人久久www免费人成看片| 亚洲国产精品sss在线观看 | 亚洲一区二区三区不卡视频| 高清欧美精品videossex| 国产在视频线精品| 岛国在线观看网站| 亚洲欧美激情在线| 午夜免费成人在线视频| 国产又爽黄色视频| av视频免费观看在线观看| 欧美黑人欧美精品刺激| 日本精品一区二区三区蜜桃| 黄色视频不卡| 国内毛片毛片毛片毛片毛片| 免费女性裸体啪啪无遮挡网站| 狠狠婷婷综合久久久久久88av| 满18在线观看网站| 一区二区三区精品91| 18禁黄网站禁片午夜丰满| 十分钟在线观看高清视频www| 欧美久久黑人一区二区| 精品久久久精品久久久| a级毛片黄视频| 欧美成人午夜精品| 久久精品国产a三级三级三级| 国产免费av片在线观看野外av| 国产成人精品在线电影| 国产黄色免费在线视频| 亚洲人成电影观看| 亚洲人成电影观看| 亚洲国产精品一区二区三区在线| 正在播放国产对白刺激| 国产成人啪精品午夜网站| 国产男女内射视频| 91大片在线观看| www.熟女人妻精品国产| 国产成人免费无遮挡视频| 中文字幕精品免费在线观看视频| 欧美亚洲日本最大视频资源| 十八禁高潮呻吟视频| 国产精品国产av在线观看| 女人久久www免费人成看片| 麻豆av在线久日| 女警被强在线播放| 亚洲一区高清亚洲精品| 正在播放国产对白刺激| 无限看片的www在线观看| 午夜视频精品福利| 国产av一区二区精品久久| 精品视频人人做人人爽| 日本精品一区二区三区蜜桃| 免费观看人在逋| 又大又爽又粗| 性少妇av在线| 亚洲精品国产一区二区精华液| 免费一级毛片在线播放高清视频 | 一本大道久久a久久精品| 精品少妇久久久久久888优播| 国产成人av激情在线播放| 国产精品电影一区二区三区 | 日韩欧美免费精品| a在线观看视频网站| 日本vs欧美在线观看视频| 黄色a级毛片大全视频| 老司机深夜福利视频在线观看| 丰满饥渴人妻一区二区三| 搡老岳熟女国产| 老汉色av国产亚洲站长工具| 精品免费久久久久久久清纯 | 在线观看www视频免费| 在线观看一区二区三区激情| 精品一品国产午夜福利视频| 久久九九热精品免费| 精品国产乱子伦一区二区三区| 纯流量卡能插随身wifi吗| 岛国毛片在线播放| 一级片'在线观看视频| 成年版毛片免费区| 日本精品一区二区三区蜜桃| 成人国语在线视频| 黑人巨大精品欧美一区二区mp4| 久久中文看片网| 国产av又大| 久久 成人 亚洲| 人人妻,人人澡人人爽秒播| 欧美日韩亚洲综合一区二区三区_| 少妇的丰满在线观看| 黄色成人免费大全| 国产免费av片在线观看野外av| 欧美乱妇无乱码| 国产99久久九九免费精品| 12—13女人毛片做爰片一| 欧美日韩中文字幕国产精品一区二区三区 | 国精品久久久久久国模美| 亚洲少妇的诱惑av| 1024视频免费在线观看| 久久国产亚洲av麻豆专区| 91成人精品电影| av超薄肉色丝袜交足视频| 91大片在线观看| 18禁观看日本| 精品福利永久在线观看| 少妇猛男粗大的猛烈进出视频| 十八禁高潮呻吟视频| 精品福利永久在线观看| 99国产精品99久久久久| 国产蜜桃级精品一区二区三区 | 好男人电影高清在线观看| 国产亚洲精品久久久久5区| 黑人猛操日本美女一级片| 久久久久久免费高清国产稀缺| 美女高潮到喷水免费观看| 亚洲免费av在线视频| 97人妻天天添夜夜摸| 亚洲成人国产一区在线观看| 亚洲成a人片在线一区二区| 国产精品国产高清国产av | 看黄色毛片网站| 欧美精品高潮呻吟av久久| 色精品久久人妻99蜜桃| netflix在线观看网站| 一区福利在线观看| 国产国语露脸激情在线看| 久久精品亚洲精品国产色婷小说| 午夜福利,免费看| 人人妻人人澡人人看| 高清欧美精品videossex| 久热这里只有精品99| av福利片在线| 精品福利观看| 91成年电影在线观看| 国产激情久久老熟女| 动漫黄色视频在线观看| 一夜夜www| 亚洲精品国产色婷婷电影| 妹子高潮喷水视频| 亚洲av欧美aⅴ国产| 窝窝影院91人妻| tube8黄色片| 18在线观看网站| 美女福利国产在线| 亚洲男人天堂网一区| 国产精品秋霞免费鲁丝片| 女性被躁到高潮视频| 日本黄色日本黄色录像| 多毛熟女@视频| 亚洲专区中文字幕在线| 欧美日韩成人在线一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 深夜精品福利| 人妻久久中文字幕网| 飞空精品影院首页| 国产精品 国内视频| 真人做人爱边吃奶动态| 久久精品国产清高在天天线| 老司机福利观看| 亚洲精品国产区一区二| 交换朋友夫妻互换小说| 精品国内亚洲2022精品成人 | 99热国产这里只有精品6| 在线十欧美十亚洲十日本专区| 人妻一区二区av| 狂野欧美激情性xxxx| 天堂俺去俺来也www色官网| 国产有黄有色有爽视频| 国产激情欧美一区二区| 国产欧美日韩一区二区三| 精品国产一区二区久久| 欧美成人午夜精品| 人成视频在线观看免费观看| 亚洲成人手机| 久久久水蜜桃国产精品网| 欧美成人免费av一区二区三区 | 久久久久久久久久久久大奶| 欧美日韩国产mv在线观看视频| 国产精品国产av在线观看| 免费不卡黄色视频| 婷婷成人精品国产| 久久精品91无色码中文字幕| 欧美乱码精品一区二区三区| 国产三级黄色录像| 国产日韩一区二区三区精品不卡| 午夜福利乱码中文字幕| 美女 人体艺术 gogo| 国产精品99久久99久久久不卡| 国产精品一区二区在线观看99| 国产亚洲欧美98| 热re99久久国产66热| 天堂俺去俺来也www色官网| 国产免费av片在线观看野外av| 丰满迷人的少妇在线观看| 亚洲欧洲精品一区二区精品久久久| www.熟女人妻精品国产| 人人澡人人妻人| 妹子高潮喷水视频| 欧美激情极品国产一区二区三区| 欧美日韩亚洲综合一区二区三区_| 国产欧美日韩综合在线一区二区| 一区二区三区国产精品乱码| 亚洲伊人色综图| 欧美日韩乱码在线| 国产精品欧美亚洲77777| netflix在线观看网站| 精品亚洲成国产av| 天天添夜夜摸| 欧美日韩视频精品一区| 欧美中文综合在线视频| 黄片播放在线免费| 国产精品99久久99久久久不卡| 激情视频va一区二区三区| 99在线人妻在线中文字幕 | 久热爱精品视频在线9| 国产av一区二区精品久久| 久久九九热精品免费| 丁香六月欧美| 乱人伦中国视频| 男女之事视频高清在线观看| 免费久久久久久久精品成人欧美视频| 亚洲五月色婷婷综合| 老汉色av国产亚洲站长工具| 性少妇av在线| 色婷婷久久久亚洲欧美| 精品视频人人做人人爽| 日韩大码丰满熟妇| 正在播放国产对白刺激| 飞空精品影院首页| 人妻一区二区av| 在线观看午夜福利视频| 69精品国产乱码久久久| 成人精品一区二区免费| 搡老岳熟女国产| 国产成人精品久久二区二区91| 一级片'在线观看视频| 国产精品亚洲一级av第二区| ponron亚洲| 女人久久www免费人成看片| 18禁裸乳无遮挡免费网站照片 | 岛国毛片在线播放| 母亲3免费完整高清在线观看| 69精品国产乱码久久久| 怎么达到女性高潮| 黑人巨大精品欧美一区二区蜜桃| 五月开心婷婷网| 青草久久国产| 国产高清videossex| 黄色成人免费大全| 午夜福利乱码中文字幕| 日韩欧美免费精品| 人人妻人人澡人人看| 两个人看的免费小视频| 国产高清视频在线播放一区| 人人妻人人澡人人爽人人夜夜| 91麻豆av在线| 黄网站色视频无遮挡免费观看| 捣出白浆h1v1| 99热网站在线观看| 欧美日韩精品网址| 欧美老熟妇乱子伦牲交| 中文字幕人妻熟女乱码| 久久精品国产综合久久久| 久久狼人影院| 国产精品偷伦视频观看了| 啦啦啦视频在线资源免费观看| 男女午夜视频在线观看| 精品国产一区二区久久| av线在线观看网站| 午夜老司机福利片| 精品电影一区二区在线| 后天国语完整版免费观看| 变态另类成人亚洲欧美熟女 | 午夜福利,免费看| 黄色丝袜av网址大全| 一级a爱片免费观看的视频| 日韩一卡2卡3卡4卡2021年| 亚洲欧美日韩另类电影网站| 一个人免费在线观看的高清视频| 亚洲在线自拍视频| 一二三四社区在线视频社区8| 久久热在线av| 国产免费av片在线观看野外av| 亚洲精品国产一区二区精华液| 一级a爱片免费观看的视频| 在线av久久热| 午夜福利视频在线观看免费| 久久精品国产99精品国产亚洲性色 | 免费不卡黄色视频| 国产野战对白在线观看| 最新在线观看一区二区三区| 日韩成人在线观看一区二区三区| 亚洲一区高清亚洲精品| 天天影视国产精品| 亚洲精品国产一区二区精华液| 国产一区在线观看成人免费| 老司机福利观看| 在线观看一区二区三区激情| 高清av免费在线| 国产在线观看jvid| 亚洲精品国产一区二区精华液| 天天影视国产精品| 国产精品免费大片| 日韩欧美一区视频在线观看| 中文字幕最新亚洲高清| 日本精品一区二区三区蜜桃| 夜夜躁狠狠躁天天躁| 亚洲久久久国产精品| 夜夜躁狠狠躁天天躁| 欧美黑人欧美精品刺激| 亚洲 欧美一区二区三区| 欧美日本中文国产一区发布| 亚洲中文日韩欧美视频| 亚洲av电影在线进入| 看片在线看免费视频| 在线天堂中文资源库| 窝窝影院91人妻| 亚洲中文av在线| 美女福利国产在线| 久久人妻福利社区极品人妻图片| 老司机亚洲免费影院| 18禁美女被吸乳视频| 国产免费男女视频| 满18在线观看网站| 在线观看66精品国产| 亚洲精华国产精华精| 亚洲色图av天堂| 亚洲精华国产精华精| 麻豆成人av在线观看| 嫁个100分男人电影在线观看| 在线观看免费视频日本深夜| 一区二区三区精品91| 国产无遮挡羞羞视频在线观看| 国产精品偷伦视频观看了| 不卡一级毛片| 精品熟女少妇八av免费久了| 亚洲少妇的诱惑av| 精品午夜福利视频在线观看一区| avwww免费| www.熟女人妻精品国产| 国产精品二区激情视频| 亚洲色图综合在线观看| 色94色欧美一区二区| 91九色精品人成在线观看| 中文字幕人妻丝袜制服| 精品欧美一区二区三区在线| 国产在线观看jvid| 免费av中文字幕在线| 久久精品国产清高在天天线| xxxhd国产人妻xxx| 国产一区二区三区在线臀色熟女 | 男女午夜视频在线观看| 久久久国产一区二区| 久久精品国产亚洲av香蕉五月 | 五月开心婷婷网| 王馨瑶露胸无遮挡在线观看| 日韩欧美免费精品| 老司机在亚洲福利影院| 99riav亚洲国产免费| 妹子高潮喷水视频| 午夜老司机福利片| 亚洲第一av免费看| 国产片内射在线| 99精品欧美一区二区三区四区| 人人妻人人添人人爽欧美一区卜| 午夜福利免费观看在线| 免费日韩欧美在线观看| 最近最新免费中文字幕在线| 日本黄色视频三级网站网址 | 在线观看一区二区三区激情| 丰满人妻熟妇乱又伦精品不卡| 人妻一区二区av| 丁香六月欧美| 久久人妻熟女aⅴ| 欧美乱色亚洲激情| 日韩人妻精品一区2区三区| 少妇 在线观看| 大香蕉久久成人网| 韩国精品一区二区三区| 两人在一起打扑克的视频| 亚洲欧美一区二区三区黑人| 国产精品乱码一区二三区的特点 | 99精国产麻豆久久婷婷| 岛国在线观看网站| 日韩大码丰满熟妇| 欧美日韩亚洲高清精品| 日本撒尿小便嘘嘘汇集6| 免费黄频网站在线观看国产| 高清在线国产一区| 老汉色av国产亚洲站长工具| 亚洲一区二区三区不卡视频| 电影成人av| 999久久久精品免费观看国产| 99精品久久久久人妻精品| 日本五十路高清| 亚洲第一欧美日韩一区二区三区| 80岁老熟妇乱子伦牲交| 黄色怎么调成土黄色| 国产一区二区三区在线臀色熟女 | 亚洲精品久久成人aⅴ小说| 久久人妻av系列| 精品人妻熟女毛片av久久网站| 真人做人爱边吃奶动态| 男女之事视频高清在线观看| 国产亚洲av高清不卡| 亚洲熟女精品中文字幕| 99久久人妻综合| 人妻久久中文字幕网| 亚洲第一青青草原| 精品高清国产在线一区| 亚洲av日韩在线播放| 变态另类成人亚洲欧美熟女 | 十八禁人妻一区二区| 在线十欧美十亚洲十日本专区| 在线观看午夜福利视频| 午夜视频精品福利| 成人国语在线视频| 他把我摸到了高潮在线观看| 久久 成人 亚洲| 精品免费久久久久久久清纯 | 精品无人区乱码1区二区| 免费少妇av软件| 欧美大码av| 视频区欧美日本亚洲| 色婷婷久久久亚洲欧美| 成人手机av| 成人国产一区最新在线观看| 亚洲avbb在线观看| 一二三四社区在线视频社区8| 精品久久久精品久久久| 久久天躁狠狠躁夜夜2o2o| 欧美 日韩 精品 国产| 免费不卡黄色视频| 午夜福利一区二区在线看| 我的亚洲天堂| 欧美日韩乱码在线| 日韩中文字幕欧美一区二区| 99久久精品国产亚洲精品| av在线播放免费不卡| 欧美久久黑人一区二区| 99国产综合亚洲精品| 97人妻天天添夜夜摸| 99精品久久久久人妻精品| 国产精品久久久av美女十八| 人人澡人人妻人| 久久中文字幕一级| 国产精品免费视频内射| 别揉我奶头~嗯~啊~动态视频| 国产精品久久久av美女十八| 黄色a级毛片大全视频| 久久中文字幕一级| 亚洲成人免费电影在线观看| 两人在一起打扑克的视频| 日韩欧美一区二区三区在线观看 | 国产精品免费一区二区三区在线 | 欧美日本中文国产一区发布| 欧美精品啪啪一区二区三区| 一区在线观看完整版| 搡老岳熟女国产| 91老司机精品| 欧美日韩亚洲综合一区二区三区_| 精品少妇一区二区三区视频日本电影| 波多野结衣一区麻豆| 丁香欧美五月| 久久人人爽av亚洲精品天堂| 久久香蕉国产精品| 9191精品国产免费久久| 久久人妻熟女aⅴ| 999久久久精品免费观看国产| 成人国语在线视频| 高清视频免费观看一区二区| 一本综合久久免费| 欧美日韩成人在线一区二区| 免费黄频网站在线观看国产| 我的亚洲天堂| 亚洲一码二码三码区别大吗| 国产一区二区三区视频了| 可以免费在线观看a视频的电影网站| 一本一本久久a久久精品综合妖精| 久久精品亚洲精品国产色婷小说| 午夜精品在线福利| 亚洲国产毛片av蜜桃av| 精品第一国产精品| 欧美日韩亚洲综合一区二区三区_| 建设人人有责人人尽责人人享有的| 国产97色在线日韩免费| 欧美黄色片欧美黄色片| 不卡av一区二区三区| 两性夫妻黄色片| 精品久久久久久,| 黄色怎么调成土黄色| 亚洲av熟女| 亚洲全国av大片| 国内久久婷婷六月综合欲色啪| 女人久久www免费人成看片| av免费在线观看网站| 精品一区二区三区四区五区乱码| 宅男免费午夜| av不卡在线播放| 午夜精品久久久久久毛片777| 巨乳人妻的诱惑在线观看| 亚洲 欧美一区二区三区| 夜夜爽天天搞| 人人妻,人人澡人人爽秒播| 久久人妻av系列| 亚洲av第一区精品v没综合| 日本黄色视频三级网站网址 | 亚洲国产欧美日韩在线播放| 国产欧美日韩一区二区精品| 韩国精品一区二区三区| 久久久久久亚洲精品国产蜜桃av| 色在线成人网| 午夜成年电影在线免费观看| 中文字幕人妻丝袜一区二区| 在线观看免费高清a一片| 久久人妻福利社区极品人妻图片| 天天躁日日躁夜夜躁夜夜| 国产精品自产拍在线观看55亚洲 | av有码第一页| 亚洲精品av麻豆狂野| 国产精品av久久久久免费| 一边摸一边抽搐一进一出视频| 一边摸一边抽搐一进一小说 | 黄色片一级片一级黄色片| 欧美精品高潮呻吟av久久| 国产男女内射视频| 国产精品国产av在线观看| 高潮久久久久久久久久久不卡| 久久久水蜜桃国产精品网| 亚洲午夜理论影院| 亚洲欧美一区二区三区黑人| 变态另类成人亚洲欧美熟女 | 男女午夜视频在线观看| 国产亚洲欧美98| 精品一品国产午夜福利视频| 国产日韩一区二区三区精品不卡| 久久香蕉国产精品| 国产亚洲精品一区二区www | 热99国产精品久久久久久7| 视频在线观看一区二区三区| 亚洲国产欧美日韩在线播放| 两人在一起打扑克的视频| 成年动漫av网址| 欧美精品一区二区免费开放| 国产精品综合久久久久久久免费 | 国产深夜福利视频在线观看| 色综合欧美亚洲国产小说| 免费在线观看完整版高清|