• <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.

    色婷婷久久久亚洲欧美| 亚洲精品久久成人aⅴ小说| 男女下面进入的视频免费午夜 | 麻豆av在线久日| 欧美日韩一级在线毛片| 九色国产91popny在线| 欧美乱码精品一区二区三区| 亚洲欧美日韩无卡精品| 久久狼人影院| 中文字幕人成人乱码亚洲影| 国产一区二区三区视频了| 色精品久久人妻99蜜桃| 十八禁网站免费在线| 欧美av亚洲av综合av国产av| 此物有八面人人有两片| 国产一区二区三区视频了| 亚洲一区二区三区色噜噜| 国产午夜福利久久久久久| 欧美三级亚洲精品| 国产午夜福利久久久久久| 999精品在线视频| 妹子高潮喷水视频| 久久久久久人人人人人| 成人国产一区最新在线观看| 日本精品一区二区三区蜜桃| 两性夫妻黄色片| www.精华液| 亚洲一区二区三区不卡视频| 久久人人精品亚洲av| 亚洲精品美女久久av网站| 国产精品久久视频播放| 久久久久精品国产欧美久久久| 搞女人的毛片| 不卡一级毛片| 长腿黑丝高跟| 成年人黄色毛片网站| www.精华液| 成人永久免费在线观看视频| 美女大奶头视频| av有码第一页| 非洲黑人性xxxx精品又粗又长| 亚洲人成网站高清观看| 在线观看www视频免费| 精品国产美女av久久久久小说| 大型av网站在线播放| 在线观看日韩欧美| 淫妇啪啪啪对白视频| 一级a爱视频在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产高清有码在线观看视频 | 免费看十八禁软件| 正在播放国产对白刺激| 国产乱人伦免费视频| 久久久久亚洲av毛片大全| 精品人妻1区二区| АⅤ资源中文在线天堂| 久久久久久久久久黄片| 亚洲国产看品久久| 午夜福利成人在线免费观看| 久久久精品国产亚洲av高清涩受| 国产色视频综合| 久久亚洲真实| 国产精品久久视频播放| 国产精华一区二区三区| 香蕉国产在线看| 麻豆成人午夜福利视频| 国产成人av教育| 性欧美人与动物交配| 欧美中文综合在线视频| 午夜成年电影在线免费观看| 亚洲人成网站高清观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲人成网站在线播放欧美日韩| 成人特级黄色片久久久久久久| av天堂在线播放| 色婷婷久久久亚洲欧美| 国产久久久一区二区三区| 免费看日本二区| av在线天堂中文字幕| 国产精品久久久久久人妻精品电影| 久久香蕉激情| 精品久久久久久成人av| 欧美绝顶高潮抽搐喷水| 国产精品综合久久久久久久免费| 国产在线观看jvid| 国产亚洲欧美在线一区二区| 又大又爽又粗| 久久精品成人免费网站| 亚洲人成电影免费在线| 久久草成人影院| 看片在线看免费视频| 亚洲最大成人中文| 嫩草影视91久久| 香蕉av资源在线| 久久久国产精品麻豆| 亚洲五月天丁香| 桃红色精品国产亚洲av| 精品日产1卡2卡| 精品久久久久久久人妻蜜臀av| 亚洲av五月六月丁香网| 村上凉子中文字幕在线| 哪里可以看免费的av片| av在线天堂中文字幕| 美女 人体艺术 gogo| bbb黄色大片| 操出白浆在线播放| 国产成人一区二区三区免费视频网站| 可以在线观看的亚洲视频| 麻豆国产av国片精品| 国产真实乱freesex| 成年女人毛片免费观看观看9| 非洲黑人性xxxx精品又粗又长| 三级毛片av免费| 久久久久国产一级毛片高清牌| 成人18禁高潮啪啪吃奶动态图| 天堂影院成人在线观看| 国产片内射在线| 性欧美人与动物交配| 久久久久精品国产欧美久久久| 两性夫妻黄色片| 男女之事视频高清在线观看| 深夜精品福利| 啦啦啦 在线观看视频| 欧美av亚洲av综合av国产av| 999久久久国产精品视频| 日韩欧美 国产精品| 黄色毛片三级朝国网站| 搞女人的毛片| 看片在线看免费视频| 国产aⅴ精品一区二区三区波| 欧美丝袜亚洲另类 | 久久天堂一区二区三区四区| 欧美精品亚洲一区二区| 国产精品99久久99久久久不卡| 91在线观看av| 国产精品电影一区二区三区| 国产亚洲欧美98| 青草久久国产| 精品熟女少妇八av免费久了| 在线看三级毛片| 久久精品91无色码中文字幕| 欧美激情高清一区二区三区| 黄色视频不卡| 亚洲国产精品999在线| 99国产综合亚洲精品| 真人一进一出gif抽搐免费| 午夜久久久在线观看| 欧美绝顶高潮抽搐喷水| 久久久久久九九精品二区国产 | 女生性感内裤真人,穿戴方法视频| 俺也久久电影网| 在线播放国产精品三级| 久久国产亚洲av麻豆专区| 午夜两性在线视频| 美女大奶头视频| 啦啦啦 在线观看视频| 欧美黑人巨大hd| 99国产精品一区二区蜜桃av| 日本一区二区免费在线视频| 制服人妻中文乱码| 成人国产综合亚洲| 国产爱豆传媒在线观看 | 法律面前人人平等表现在哪些方面| 窝窝影院91人妻| 成人欧美大片| 国产乱人伦免费视频| 熟女电影av网| 久久久久亚洲av毛片大全| 国产精品久久久人人做人人爽| 亚洲熟妇熟女久久| 国产亚洲精品久久久久5区| 欧美乱妇无乱码| 精品欧美一区二区三区在线| 两人在一起打扑克的视频| 国产1区2区3区精品| 此物有八面人人有两片| 久久中文字幕一级| 国产v大片淫在线免费观看| 在线十欧美十亚洲十日本专区| 国产一区在线观看成人免费| 日韩欧美免费精品| av欧美777| 99久久久亚洲精品蜜臀av| 国语自产精品视频在线第100页| 国产精品一区二区免费欧美| 一级黄色大片毛片| 日韩欧美三级三区| 久久精品aⅴ一区二区三区四区| 在线天堂中文资源库| 亚洲成av人片免费观看| 亚洲精品国产一区二区精华液| 麻豆一二三区av精品| 老司机深夜福利视频在线观看| 韩国精品一区二区三区| 免费电影在线观看免费观看| 俺也久久电影网| 欧洲精品卡2卡3卡4卡5卡区| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久久久久亚洲av鲁大| 美女高潮喷水抽搐中文字幕| 日本免费a在线| 国产极品粉嫩免费观看在线| 久久久久久久久免费视频了| av免费在线观看网站| 午夜福利一区二区在线看| 十八禁网站免费在线| 一本久久中文字幕| 熟妇人妻久久中文字幕3abv| 欧美黑人欧美精品刺激| 在线观看日韩欧美| 国产精品免费视频内射| 国产精品自产拍在线观看55亚洲| 久久午夜综合久久蜜桃| 亚洲自拍偷在线| 1024视频免费在线观看| 国产伦一二天堂av在线观看| 日韩三级视频一区二区三区| 人人妻人人澡人人看| 色婷婷久久久亚洲欧美| 国产欧美日韩精品亚洲av| 狂野欧美激情性xxxx| 欧美三级亚洲精品| 99国产精品一区二区蜜桃av| 精品国产国语对白av| 丰满的人妻完整版| www国产在线视频色| 欧美三级亚洲精品| 大型av网站在线播放| 欧美中文综合在线视频| 动漫黄色视频在线观看| 香蕉丝袜av| 亚洲国产中文字幕在线视频| 久久久久久人人人人人| 国产精品久久视频播放| 亚洲自偷自拍图片 自拍| 亚洲成人国产一区在线观看| 久久草成人影院| 中文亚洲av片在线观看爽| 亚洲人成77777在线视频| 午夜成年电影在线免费观看| 精品久久久久久,| 两性夫妻黄色片| 免费在线观看视频国产中文字幕亚洲| 亚洲熟妇熟女久久| 亚洲国产精品久久男人天堂| 欧美av亚洲av综合av国产av| 美女 人体艺术 gogo| 亚洲精品一卡2卡三卡4卡5卡| 深夜精品福利| 老司机午夜福利在线观看视频| 色哟哟哟哟哟哟| 69av精品久久久久久| 看免费av毛片| 香蕉av资源在线| 国产97色在线日韩免费| 日韩免费av在线播放| 亚洲aⅴ乱码一区二区在线播放 | 老司机在亚洲福利影院| 91av网站免费观看| 波多野结衣巨乳人妻| 精品国产美女av久久久久小说| 欧美色欧美亚洲另类二区| 成人特级黄色片久久久久久久| 国产亚洲精品久久久久5区| 日韩大尺度精品在线看网址| 国产免费av片在线观看野外av| 18禁黄网站禁片午夜丰满| 久久国产精品影院| 亚洲国产精品成人综合色| 看黄色毛片网站| 国产精品免费一区二区三区在线| 色av中文字幕| ponron亚洲| 久久精品成人免费网站| 脱女人内裤的视频| 亚洲精品色激情综合| 91麻豆av在线| 99国产极品粉嫩在线观看| 俺也久久电影网| 欧美精品啪啪一区二区三区| 99在线人妻在线中文字幕| 熟女少妇亚洲综合色aaa.| www.自偷自拍.com| 精品熟女少妇八av免费久了| 亚洲一区二区三区色噜噜| 亚洲国产精品成人综合色| 亚洲av片天天在线观看| 国产一区二区在线av高清观看| 国产成人啪精品午夜网站| 亚洲精品在线观看二区| 在线观看免费午夜福利视频| 老汉色∧v一级毛片| 亚洲国产欧美一区二区综合| 两性午夜刺激爽爽歪歪视频在线观看 | 精品高清国产在线一区| 国产精品二区激情视频| 制服人妻中文乱码| 两个人免费观看高清视频| 99久久精品国产亚洲精品| 免费在线观看影片大全网站| 成人手机av| 免费无遮挡裸体视频| 视频区欧美日本亚洲| 欧美大码av| 一a级毛片在线观看| av片东京热男人的天堂| 大型黄色视频在线免费观看| 久9热在线精品视频| 美国免费a级毛片| 国产一区二区在线av高清观看| 色尼玛亚洲综合影院| 色在线成人网| 最新在线观看一区二区三区| 亚洲 国产 在线| 在线观看免费日韩欧美大片| 久热爱精品视频在线9| 色综合站精品国产| 国产激情偷乱视频一区二区| 亚洲久久久国产精品| 精品日产1卡2卡| 亚洲av成人不卡在线观看播放网| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩中文字幕欧美一区二区| 啪啪无遮挡十八禁网站| 激情在线观看视频在线高清| 欧美激情高清一区二区三区| 国产伦在线观看视频一区| 久久久水蜜桃国产精品网| 中文在线观看免费www的网站 | 久久精品成人免费网站| 日日摸夜夜添夜夜添小说| 91麻豆精品激情在线观看国产| 成人国语在线视频| 亚洲自拍偷在线| 精品无人区乱码1区二区| 欧美日韩黄片免| 1024视频免费在线观看| 少妇 在线观看| 亚洲七黄色美女视频| 欧美乱码精品一区二区三区| 嫩草影院精品99| 亚洲av中文字字幕乱码综合 | 午夜日韩欧美国产| 露出奶头的视频| 一级毛片女人18水好多| 哪里可以看免费的av片| 可以在线观看毛片的网站| 国产一卡二卡三卡精品| √禁漫天堂资源中文www| 巨乳人妻的诱惑在线观看| 十分钟在线观看高清视频www| 一进一出抽搐gif免费好疼| 午夜免费激情av| 变态另类成人亚洲欧美熟女| 在线观看午夜福利视频| 日本成人三级电影网站| 1024视频免费在线观看| 一边摸一边做爽爽视频免费| 视频在线观看一区二区三区| 亚洲人成网站在线播放欧美日韩| 中国美女看黄片| 婷婷精品国产亚洲av| 叶爱在线成人免费视频播放| 欧美乱码精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲精品一区二区www| 亚洲狠狠婷婷综合久久图片| 搡老熟女国产l中国老女人| 成人免费观看视频高清| 99久久99久久久精品蜜桃| 欧美性猛交黑人性爽| 国产精品电影一区二区三区| 女性生殖器流出的白浆| 国产成人精品久久二区二区免费| 亚洲av电影在线进入| 亚洲精品av麻豆狂野| 亚洲熟妇中文字幕五十中出| 可以在线观看毛片的网站| 99在线人妻在线中文字幕| 怎么达到女性高潮| 黄色视频不卡| 在线国产一区二区在线| 成在线人永久免费视频| 搡老熟女国产l中国老女人| 国产一区二区三区在线臀色熟女| 搡老熟女国产l中国老女人| 听说在线观看完整版免费高清| 色综合婷婷激情| 一个人免费在线观看的高清视频| 99久久99久久久精品蜜桃| 91成年电影在线观看| 丁香六月欧美| 亚洲av日韩精品久久久久久密| 看片在线看免费视频| 女性生殖器流出的白浆| 欧美乱码精品一区二区三区| 免费女性裸体啪啪无遮挡网站| 欧美绝顶高潮抽搐喷水| 人成视频在线观看免费观看| 在线观看66精品国产| www.www免费av| 黄网站色视频无遮挡免费观看| 一区二区三区精品91| 国产亚洲av高清不卡| 久久久久九九精品影院| 国产精品精品国产色婷婷| 怎么达到女性高潮| 久久久久久久午夜电影| netflix在线观看网站| 久久久精品欧美日韩精品| 女人爽到高潮嗷嗷叫在线视频| 久久中文字幕一级| 欧美 亚洲 国产 日韩一| 狠狠狠狠99中文字幕| 成人国产综合亚洲| 这个男人来自地球电影免费观看| 在线观看www视频免费| 后天国语完整版免费观看| 国产高清视频在线播放一区| 最近在线观看免费完整版| 欧美亚洲日本最大视频资源| 国产成人av教育| av中文乱码字幕在线| 国产精品亚洲av一区麻豆| 亚洲国产毛片av蜜桃av| 中文字幕精品亚洲无线码一区 | 欧美成人一区二区免费高清观看 | 香蕉久久夜色| 欧美成人性av电影在线观看| 久久国产精品影院| 好男人在线观看高清免费视频 | 欧美一级毛片孕妇| 久久精品夜夜夜夜夜久久蜜豆 | 日本一本二区三区精品| 国产一区二区在线av高清观看| 精品高清国产在线一区| 国产黄色小视频在线观看| 成人国产综合亚洲| 久久久久九九精品影院| 久久久久久九九精品二区国产 | 精品卡一卡二卡四卡免费| 免费高清在线观看日韩| 国产免费男女视频| 免费人成视频x8x8入口观看| 亚洲 欧美 日韩 在线 免费| 可以在线观看的亚洲视频| 午夜免费成人在线视频| www国产在线视频色| avwww免费| 高清在线国产一区| 久久精品国产99精品国产亚洲性色| 白带黄色成豆腐渣| а√天堂www在线а√下载| 中文字幕高清在线视频| 日本精品一区二区三区蜜桃| 一本精品99久久精品77| 在线观看一区二区三区| 欧美一级毛片孕妇| 好男人电影高清在线观看| 免费女性裸体啪啪无遮挡网站| 国产成人啪精品午夜网站| 欧美久久黑人一区二区| 最近最新中文字幕大全免费视频| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品合色在线| 美女免费视频网站| a级毛片a级免费在线| 国产成人啪精品午夜网站| 欧美三级亚洲精品| 19禁男女啪啪无遮挡网站| 三级毛片av免费| 一本大道久久a久久精品| 一a级毛片在线观看| 一级a爱片免费观看的视频| 美女午夜性视频免费| 久久天躁狠狠躁夜夜2o2o| 欧美乱码精品一区二区三区| 国产一区二区三区在线臀色熟女| 99久久99久久久精品蜜桃| 亚洲精品美女久久久久99蜜臀| 人人澡人人妻人| 在线观看一区二区三区| 桃色一区二区三区在线观看| 久久久久久九九精品二区国产 | 国产精品av久久久久免费| 精品久久久久久久久久免费视频| www日本黄色视频网| 在线观看免费午夜福利视频| 精品一区二区三区视频在线观看免费| 人成视频在线观看免费观看| 国产免费男女视频| 成在线人永久免费视频| 精品高清国产在线一区| 国产一区二区在线av高清观看| 日本一本二区三区精品| 香蕉av资源在线| 色尼玛亚洲综合影院| 男女下面进入的视频免费午夜 | 不卡av一区二区三区| 久久中文字幕人妻熟女| 真人做人爱边吃奶动态| 美女高潮喷水抽搐中文字幕| 丰满人妻熟妇乱又伦精品不卡| 一区福利在线观看| 熟女少妇亚洲综合色aaa.| 久久国产精品男人的天堂亚洲| 精品欧美一区二区三区在线| 国产成人系列免费观看| 老司机午夜福利在线观看视频| 亚洲美女黄片视频| 性欧美人与动物交配| 久久久国产欧美日韩av| 首页视频小说图片口味搜索| 国产一区在线观看成人免费| √禁漫天堂资源中文www| 国产蜜桃级精品一区二区三区| 首页视频小说图片口味搜索| 韩国精品一区二区三区| 日本一区二区免费在线视频| 老熟妇仑乱视频hdxx| 欧美成狂野欧美在线观看| 在线观看www视频免费| 在线观看免费日韩欧美大片| 亚洲国产高清在线一区二区三 | 午夜福利成人在线免费观看| 午夜a级毛片| 高清在线国产一区| 宅男免费午夜| 色综合欧美亚洲国产小说| 国产精品野战在线观看| 琪琪午夜伦伦电影理论片6080| 搡老熟女国产l中国老女人| 狠狠狠狠99中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 久久精品成人免费网站| 中文亚洲av片在线观看爽| 十八禁人妻一区二区| 在线观看免费视频日本深夜| 给我免费播放毛片高清在线观看| 国产精品亚洲一级av第二区| 欧美在线一区亚洲| 人人妻人人看人人澡| 亚洲精品粉嫩美女一区| 中文字幕久久专区| videosex国产| 欧美性长视频在线观看| 亚洲av第一区精品v没综合| 少妇熟女aⅴ在线视频| 观看免费一级毛片| 岛国视频午夜一区免费看| 怎么达到女性高潮| 日韩视频一区二区在线观看| 在线十欧美十亚洲十日本专区| 99久久无色码亚洲精品果冻| 精品不卡国产一区二区三区| 日本 欧美在线| 一二三四在线观看免费中文在| а√天堂www在线а√下载| 亚洲va日本ⅴa欧美va伊人久久| а√天堂www在线а√下载| 欧美精品啪啪一区二区三区| 99国产精品一区二区蜜桃av| 亚洲色图av天堂| 1024视频免费在线观看| 中文亚洲av片在线观看爽| 青草久久国产| 757午夜福利合集在线观看| 国产精品久久久久久人妻精品电影| 国产精品一区二区精品视频观看| 中国美女看黄片| 亚洲精品在线美女| 午夜精品久久久久久毛片777| 大型黄色视频在线免费观看| 国产日本99.免费观看| 国产av在哪里看| 国产v大片淫在线免费观看| 亚洲成人精品中文字幕电影| 人妻久久中文字幕网| 久久狼人影院| 丰满人妻熟妇乱又伦精品不卡| 精品乱码久久久久久99久播| 人人澡人人妻人| 亚洲成人国产一区在线观看| 人人澡人人妻人| 狠狠狠狠99中文字幕| ponron亚洲| √禁漫天堂资源中文www| or卡值多少钱| 亚洲 欧美 日韩 在线 免费| 中出人妻视频一区二区| 99在线人妻在线中文字幕| 老司机在亚洲福利影院| 99热6这里只有精品| 国产精品久久久人人做人人爽| 中文资源天堂在线| 国产色视频综合| 午夜福利欧美成人| 亚洲av中文字字幕乱码综合 | 757午夜福利合集在线观看| 精品熟女少妇八av免费久了| 国产免费av片在线观看野外av| 久久久久免费精品人妻一区二区 | 99国产精品一区二区三区| 欧美色视频一区免费| 国内揄拍国产精品人妻在线 | 亚洲精品色激情综合| 在线十欧美十亚洲十日本专区| 亚洲成人久久爱视频| 91麻豆精品激情在线观看国产| 欧美绝顶高潮抽搐喷水| 欧美日韩福利视频一区二区| avwww免费| 91九色精品人成在线观看| 国产欧美日韩一区二区精品| 草草在线视频免费看|