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

    Frequency-SPeed Control Model Identification of Ultrasonic Motor Using SteP ResPonse

    2015-11-24 02:39:15ShiJingzhuo史敬灼ZhangCaixia張彩霞
    關(guān)鍵詞:彩霞

    Shi Jingzhuo(史敬灼)'Zhang Caixia(張彩霞)

    College of Electrical Engineering'Henan University of Science and Technology'Luoyang 471023'P.R.China

    Frequency-SPeed Control Model Identification of Ultrasonic Motor Using SteP ResPonse

    Shi Jingzhuo(史敬灼)*'Zhang Caixia(張彩霞)

    College of Electrical Engineering'Henan University of Science and Technology'Luoyang 471023'P.R.China

    Control model of ultrasonic motor is the foundation for high control performance.The frequency of driving voltage is commonly used as control variable in the speed control system of ultrasonic motor.Speed control model with the input frequency can significantly improve speed control performance.Step response of rotating speed is tested.Then'the transfer function model is identified through characteristic point method.Considering time-varying characteristics of the model parameters'the variables are fitted with frequency and speed as the independent variables'and the variable model of ultrasonic motor system is obtained'with consideration of the nonlinearity of ultrasonic motor system.The proposed model can be used in the design and analysis of the speed control system in ultrasonic motor.

    ultrasonic motor(USM);speed control;characteristic point method

    0 Introduction

    Ultrasonic motor(USM)is a kind of special motor'and it uses the inverse piezoelectric effect of the piezoelectric material to make the stator generate mechanical vibration within the ultrasonic frequency band.The rotor is driven through the friction between the stator and the rotor[1-4]. USM contains a series of nonlinear processes such as piezoelectric energy conversion'friction energy transfer'and so on'so it is nonlinear'time-varying and strong coupling.Therefore'it is difficult to realize high-precision motion control.

    The mathematical model of the control object is an important foundation of the control system analysis'design and performance evaluation. More reasonable control strategies can also use the model[5-11].Eut the modeling problem of USM has not been solved mainly because of the complexity of its operation mechanism and its short research history.The existing researches mainly concentrate on theoretical modeling and numerical modeling.These methods are all based on piezoelectricity and friction theoretical knowledge'and trying to establish a model completely describing the running process of USM[12-21]. Great progress has been made'and it has become a powerful tool for the analysis and design of USM.Eut these models are too complex and difficult for direct control applications.On the other hand'since our understanding of USM′s nonlinearity is not thorough enough and the expression of it is not comprehensive'the model still has potential.

    Erom the perspective of control application' the modeling of USM can also adopt other methods'e.g.system identification.The input and output signals of USM can reflect the dynamic characteristics of the motor system'and if we can select the appropriate form of the input signal' the input and the output signals can completely contain the non-linear characteristics of USM which we concern[2'3'13].Therefore'we can usethe input and the output data obtained from the test to simulate the model of USM.And the model can be directly applied to control.

    Speed control is the core of motor motion control.In the USM speed control'frequency of the driving voltage is often used as a control variable to regulate speed.The speed control model with an input frequency can significantly improve speed control performance.We adopt the identification modeling method to establish the speed control model of USM which appropriately contains nonlinear characteristics.

    1 Design of Data Test

    We intend to establish the frequency-speed control model of USM through identification method.Therefore'we need to measure the input and the output signals in the experiments as frequency and speed'respectively.The speed of the motor is the output response at a specific frequency input signal.Adopting different forms of input signal'the output response will be different.We must adopt an appropriate form of the frequency signal as the input to make sure the measured data completely reflect the characteristics of USM.The selected input signal must be sufficient to motivate all the dynamic characteristics of USM.It means that the frequency range of the input signal should be able to cover the part that we concern in the dynamic frequency range of USM.Step input signal is a common signal meeting the above requirements.The form of step signal is simple.Eesides it is easy to achieve'analyze and research.Therefore'it is more appropriate as an input signal.In this paper'we use the frequency step signal as input'and measure the step response data of speed under open loop control of the speed to identify the model of motor.

    Shinsei USR60'a two-phase traveling wave ultrasonic motor'is used as the experiment motor.And the homemade H-bridge phase-shift PWM circuit is used as drive control circuit.DC tacho-generator is coaxial rigid connection with the motor'and it is used to measure the motor′s speed.The frequency of the driving voltage is set to a desired value by adjusting the circuit.In order to ensure the accuracy of the identification' we need to obtain the exact value of actual driving frequency.In the experiment'the waveform of the motor′s driving voltage is measured in realtime'and the frequency value is obtained by processing the recorded waveform.In order to ensure the measurement data to accurately reflect the dynamic characteristics of USM'and the credibility of the motor model which is obtained by identifying'we need to confirm the parameters of input signal such as step amplitude'sampling time'length of data records according to the prior knowledge and the control performance before experiment.

    Operating frequency of USM is generally larger than its mechanical resonant frequency'the higher the frequency'the lower the motor speed. If the step amplitude of frequency input is too large'and the speed is too low'the signal-tonoise ratio of the measured signal will decrease. It is not conducive to improve the accuracy of the identification.If the given frequency is too low' and the desired speed corresponding is too high' the variation of the speed of motor will be large and may suddenly stall the open-loop running USM.The operating frequency range of experimental motor is 41.5—44 k Hz.The selected step range of the frequency input is 42.3—43.3 k Hz after open-loop operation.USM has different performance characteristics under different input frequencies due to its complex nonlinearity.Therefore'we need to measure the input and the output data'respectively'under different input step frequencies in the experiment'and to sufficiently reflect the characteristics of motor.Meanwhile' considering the uncertainties and random perturbations which may occur in the testing process' we need to measure sets of data under each frequency to eliminate those obvious deviations.The seted values of the step input frequency in experiment and the corresponding steady-state speed values of the motor are shown in Table 1.

    Table 1 Tested data of stePPing resPonse

    The data recording time needs to be long enough to contain the complete step response process.Eut it will bring unnecessary burden to the calculation of model identification and may affect the accuracy of identification if the time is too long and data quantity is too big.The experiments show that the step response time of the experimental motor is no longer than 100 ms. Therefore'the selected data recording time is 200 ms.

    In the experiment'the data recording is completed by the A/D sampling'and the sampling time directly affects the accuracy of identification. The sampling frequency should be at least twice that of the cutoff frequency of the object according to sampling theorem.If the sampling time is too large'it will cause excessive loss of information and reduce the accuracy of identification. And in the case of the same recording time'decreasing the sampling time will increase the amount of data.At the same time'because of the limitations of response speed and calculation speed of hardware'the sampling time cannot be too small.Eor the experiment motor'the required control response bandwidth is not greater than 500 Hz.And further taking into account that we need to obtain frequency information from the measured voltage waveform'the selected sampling frequency is 10 M Hz in the experiment to ensure the accuracy of measurement.

    In the experiment'we collect the driving voltage of motor and tacho-generator output voltage signal synchronously'the measured step response curve is shown in Eig.1.It only shows part of the time data to make the waveform clearly visible.A four-character indent is recommended for all paragraphs except the first one that immediately follows a section title or a subsection title.

    Eig.1 Tested step response

    2 Model Identification of Motor Based on SteP ResPonse

    We need to preprocess these data before model identification.In Eig.1'the measured data contains noise.The noise becomes more obvious when the amplitude of the tacho-generator′s output signal is small.Since the required bandwidth of the control response is not greater than 500 Hz'we take low pass filtering for the measured speed signal and the cutoff frequency of filter is 1 000 Hz.On the other hand'the start time of data recording is earlier than the action time of step input to ensure the integrity of the measurement data.Therefore'there is a section of zero speed data at the beginning of the recording data' this section is useless for identification.If we reserve those data'they would affect the accuracy of identification.Therefore'they should be deleted.Meanwhile'the DC component existing in the measurement data will also affect the accuracy of identification and can be removed using the average method.After the preprocessing'the obtained speed step response data are shown in Eig.2.These data can be used for model identification.

    Eig.2 Step response after data pretreatment

    As mentioned before'USM is a kind of nonlinear and time-varying object.In order to achieve high-performance control'the control algorithm should have the ability of online correction.It means that the control algorithm is self-adapting. Therefore'the identification motor model is selected as second-order.

    There are many different methods to identify the model'such as area method'characteristic point method'etc.These methods all utilize step response data.The area method can take full advantage of data information of every point to identify the model'and have strong ability to suppress noise.Eut it is not an optimization'and it is not suitable to the specified order model.If the selected model order is equal or close to the order of the actual object'more desirable results can be obtained using the area method.Eut during the actual control system design'we often use a loworder model to simulate the high-order object to simplify the control process.Then the obtained low-order model will inevitably have a deviation and the model accuracy is not high by using the area method.Since the measured step response curve has attenuation oscillation characteristics' we can adopt characteristic point method which is specific to second-order underdamping model to identify the model of motor.

    The transfer function of the frequency-speed control model of USM is set to

    where K=h1/f'h1is the stable value of speed'f the given value of frequency step signal'andτthe delay time.K andτcan be obtained directly from the measured data.ξandω0are all model parameters to be identified'where'is the damping coefficient andthe natural frequency.Eor the convenience of identification' Eq.(1)normalizes to obtain the standard unit transfer function of the second-order underdamping model

    The time domain expression of this normalized object′s step response is

    Erom Eqs.(4'5)and T0'we can obtain

    Using the characteristic point method and the selected motor model of Eq.(2)to identify and calculate the measured frequency and speed data. The model parameters are shown in Table 2.

    Table 2 Model Parameters identified by characteristic Point method

    The frequency-speed control model of ultrasonic motor system can be obtained by putting the model parameters(shown in Table 2)into the model Eq.(2).Ey adopting this model'the step response of ultrasonic motor system can be calculated.Then the accuracy of the model can be verified by comparing the calculated response with the measured step response.Eigs.3'4 show the comparison of the two kinds of step response.In Eigs.3'4'solid line is the curve of measured step response'and dotted line is the step response which is calculated by identification model.It is thus clear that'the model reflects the main dynamic characteristics of actual system except the irregular fluctuations during the steady-state operation of motor.

    Table 2 shows that the motor model parameters are time-varying.It is caused by the nonlinearity of USM.In order to make the model fully reflect the nonlinear characteristics of motor'the time-varying characteristics of parameters need to be represented in the motor model.Since the motor shows different characteristics under different input frequencies'we can use the variation of model parameters along with the change of frequency to represent the time-varying nonlinearity.And in the control process'the given value of frequency is the output of the speed controller.If we ignore the dynamic process of frequency regulating'we can consider the given value of frequency as the actual value'so the frequency value is known.Therefore'the frequency f can be used as a variable to fit the parametersξandω0'and they are expressed withξ(f)andω0(f)'respectively.To make the fitting function be easy to calculate online by using the control chip'e.g. DSP'we choose the quadratic polynomial function to fit according to the change rule of the model parameters'and the obtained fitting function as shown in Eqs.(8'9).

    Eig.3 Tested and simulated step response with n= 22.4 r/min

    Eig.4 Tested and simulated step response with n= 62.8 r/min

    The comparison of the before and the after fitting model step response curve is shown in Eig. 5.In Eig.5'the dashed line is the response which is obtained by calculating the fitting data.It can be seen that both coincide basically and the model fitting effect is good.The model parameter values which are calculated from the fitting function Eqs.(8'9)are shown in Table 3.In the table' the variance is between the step response data which is calculated by using fitting model and the measured data'the relative error is between the identification value of model parameters and the fitting values of model parameters.

    Eig.5 Step response with f as fitting parameter

    Table 3 Model Parameters fitted by frequency fitting function

    In the USM speed control system'the measured speed value is necessary to constitute the closed-loop control.Since the speed data is known'we can also use the speed n as the independent variable to fit the model parametersξand ω0'expressed withξ(n)andω0(n)'respectively. Adopting the output fitting function is likely to get better control effect for some control strategies.According to the identified model parameters'we choose the quadratic polynomial function to fitξandω0.The fitting functions are shown in Eqs.(10'11).

    The comparison of the before and the after fitting model step response curve is shown in Eig.6. The model fitting parameters which adopt the speed n as the independent variable are shown in Table 4.In the table'the variance is between the step response data which is calculated by fitting model and the measured data'the relative error is between the identification value of model parameters and the fitting values of model parameters.

    Eig.6 Step response with n as fitting parameter

    Table 4 Model Parameters fitted by sPeed fitting function

    3 Conclusions

    In the paper'we adopt the characteristic point method to identify the USM model according to the frequency step input data and the speed output data'the conclusions are as follows:

    (1)We need to select an appropriate form of input signal to sufficiently excite characteristics of motor by using the method of system identification to model the motor.The step signal is a kind of input signal which meets the requirement of identification.

    (2)Characteristic point method is a kind of identification method and it aims at the set motor model(second-order underdamped model)'and the identification effect is better.

    (3)Results of identification show that the model parameters are time-varying.Therefore' the model parameters are fitted to the polynomial function with frequency and speed as independent variables'respectively.And we seek timevarying rules to properly reflect the nonlinear characteristics of the motor.It is simple to implement and results good.

    Acknowledgement

    This work was supported by the National Natural Science Eoundation of China(No.U1304501).

    [1] Guo Mingsen'Hu Junhui'Zhu Hua'et al.Threedegree-of-freedom ultrasonic motor using a 5-mm-diameter piezoelectric ceramic tube[J].IEEE Transactions on Ultrasonics'Eerroelectrics'and Erequency Control'2013'60(7):1446-1452.

    [2] Shi Jingzhuo'Zhao Eujie'Shen Xiaoxi'et al.Chaotic operation and chaos control of travelling wave ultrasonic motor[J].Ultrasonics'2013'53(6):1112-1123.

    [3] Senjyu T'Nakamura M'Urasaki M'et al.Mathematical model of ultrasonic motors for speed control[C]∥IEEE'Applied Power Electronics Conference and Exposition'APEC′06.Dallas'TX'United States:IEEE'2006.

    [4] Maas Jürgen.Model-based control for ultrasonic motors[J].IEEE/ASME Transactions on Mechatronics'2000'5(2):165-180.

    [5] Pirrotta Simone'Sinatra Rosario'Meschini Alberto. A novel simulation model for ring type ultrasonic motor[J].Meccanica'2007'42(2):127-139.

    [6] Giraud Erédéric'Semail Eetty.Analysis and phase control of a piezoelectric traveling-wave ultrasonic motor for haptic stick application[J].IEEE Transactions on Industry Applications'2004'40(6):1541-1549.

    [7] Eazrafshan Eazel'Rasti Eehnood'Mojallali Hamed. Euzzy modeling and position control of a traveling wave ultrasonic motor[C]∥ICCAE'2010 The 2nd International Conference on Computer and Automation Engineering.Singapore:IEEE'2010'5:457-461.

    [8] Lin E J.Euzzy adaptive model-following position control for ultrasonic motor[J].IEEE Transactions on Power Electronics'1997'12(2):261-268.

    [9] Erdal Eekiroglu.Microcontroller-based full control of ultrasonic motor with frequency and voltage adjusting[J].Sensors and Actuators A:Physical'2008'141(1):151-159.

    [10]Chen T C'Yu C H'Tsai M C.A new driver based on dual-mode frequency and phase control for traveling-wave type ultrasonic motor[J].Energy Conversion and Management'2008'49(10):2767-2775.

    [11]Lin R J'Hwang W J'Wai R J.Ultrasonic motor servo-drive with online trained neural-network modelfollowing controller[J].IEEE Proceedings:Electric Power Applications'1998'145(2):105-110.

    [12]Pigache E'Giraud E'Lemaire-Semail E.Modelling and identification of a planar standing wave ultrasonic motor[J].EPJ Applied Physics'2006'34(1):55-65.

    [13]Zhang Jiantao'Zhang Tiemin'Liang Li.Nonlinear modeling and generalized predictive control of ultrasonic motor[J].Electric Machines and Control' 2011'15(6):50-56.

    [14]Shi Jingzhuo'You Dongmei.Characteristic model of travelling wave ultrasonic motor[J].Ultrasonics' 2014'54(2):725-730.

    [15]Zhang Xinliang'Tan Yonghong.Modelling of ultrasonic motor with dead-zone based on Hammerstein model structure[J].Journal of Zhejiang University: Science A'2008'9(1):58-64.(in Chinese)

    [16]Sun Eengyan'Qu Jianjun'Wu Peilian.New visco-elastic contact model of traveling wave ultrasonic motor with stator frictional layer[J].Journal of Harbin Institute of Technology(New Series)'2007'14(5): 705-710.(in Chinese)

    [17]Jung P L'Jong S R'Kyung P Y'et al.Characteristic analysis of a traveling wave ultrasonic motor using an ellipsoidal static contact model[J].Smart Materials and Structures'2009'18(11):1-10.

    [18]Shi Jingzhuo'Lv Lin'Zhang Yu.Dynamic Takagi-Sugeno model for the control of ultrasonic motor[J]. Journal of Control Science and Engineering'2011' 2011:1-9.

    [19]Liew Pay Jun'Shimada Keita1'Mizutani Masayoshi' et al.Eabrication of microstructures on RE-SiC by ultrasonic cavitation assisted micro-electrical discharge machining[J].International Journal of Automation Technology'2013'7(6):621-629.

    [20]Qu Jianjun'Sun Eengyan'Zhao Chunsheng.Performance evaluation of traveling wave ultrasonic motor based on a model with visco-elastic friction layer on stator[J].Ultrasonics'2006'45(1/2/3/4):22-31.

    [21]Wang Jingsi'Shimada Keita'Mizutani Masayoshi'et al.Material removal during ultrasonic machining using smoothed particle hydrodynamics[J].International Journal of Automation Technology'2013'7(6):614-620.

    (Executive editor:Xu Chengting)

    TM383 Document code:A Article ID:1005-1120(2015)02-0192-07

    *CorresPonding author:Shi Jingzhuo'Professor'E-mail:sjznew@163.com.

    How to cite this article:Shi Jingzhuo'Zhang Caixia.Erequency-speed control model identification of ultrasonic motor using step response[J].Trans.Nanjing U.Aero.Astro.'2015'32(2):192-198.

    http://dx.doi.org/10.16356/j.1005-1120.2015.02.192

    (Received 5 January 2015;revised 24 January 2015;accepted 2 Eebruary 2015)

    猜你喜歡
    彩霞
    《流失》之《歲月》國畫
    合歡花和彩霞
    沒有愛不能到達(dá)的地方
    做人與處世(2022年4期)2022-05-26 21:50:38
    Time-Consistent Investment and Reinsurance Problems for Mean-Variance Insurers with Default Risk Under Variance Premium Principle
    《云岡石窟》油畫 張彩霞
    安其的田地
    清而不俗 妍而不媚 馮彩霞制蓮子壺
    繞轉(zhuǎn)小舌頭
    Graphene Based Electrochemical Sensor for the Detection of Volatile Organic Compounds
    A Teaching Plan for Lesson 3, Unit 16
    国产精品综合久久久久久久免费| 日本精品一区二区三区蜜桃| 999精品在线视频| 激情在线观看视频在线高清| 亚洲欧美精品综合一区二区三区| 九九在线视频观看精品| 亚洲欧美精品综合一区二区三区| 亚洲精品在线美女| 亚洲午夜精品一区,二区,三区| 又黄又爽又免费观看的视频| 少妇的逼水好多| 国产亚洲精品一区二区www| 麻豆一二三区av精品| 麻豆一二三区av精品| 熟妇人妻久久中文字幕3abv| 香蕉久久夜色| 久久精品国产综合久久久| 亚洲av成人一区二区三| 亚洲av片天天在线观看| 观看免费一级毛片| 淫妇啪啪啪对白视频| 99热这里只有是精品50| 国产综合懂色| 日本三级黄在线观看| 亚洲第一欧美日韩一区二区三区| 久久久久国产一级毛片高清牌| 狂野欧美白嫩少妇大欣赏| 亚洲成人久久爱视频| 麻豆成人午夜福利视频| 香蕉国产在线看| a在线观看视频网站| 给我免费播放毛片高清在线观看| 这个男人来自地球电影免费观看| 欧美成人一区二区免费高清观看 | 99热这里只有是精品50| 亚洲国产精品久久男人天堂| 老司机午夜福利在线观看视频| 精品久久蜜臀av无| 欧美中文综合在线视频| 日本三级黄在线观看| 国产综合懂色| 日日摸夜夜添夜夜添小说| 亚洲国产日韩欧美精品在线观看 | 怎么达到女性高潮| 一级毛片精品| 午夜日韩欧美国产| 亚洲自偷自拍图片 自拍| 好看av亚洲va欧美ⅴa在| 老司机在亚洲福利影院| 一区二区三区激情视频| 国产av在哪里看| 村上凉子中文字幕在线| 国内揄拍国产精品人妻在线| 丰满人妻一区二区三区视频av | 精品电影一区二区在线| 熟女电影av网| 国产极品精品免费视频能看的| 每晚都被弄得嗷嗷叫到高潮| 大型黄色视频在线免费观看| 伊人久久大香线蕉亚洲五| 丁香欧美五月| 久久这里只有精品中国| 美女高潮喷水抽搐中文字幕| 国产午夜精品久久久久久| 久久久久久国产a免费观看| 国产一区二区三区在线臀色熟女| 少妇裸体淫交视频免费看高清| 欧美激情在线99| 日日干狠狠操夜夜爽| 狠狠狠狠99中文字幕| 欧美黑人巨大hd| 国产高清有码在线观看视频| 国产精品久久久久久人妻精品电影| 欧美在线一区亚洲| 精品久久久久久久末码| 久久久久九九精品影院| 色精品久久人妻99蜜桃| 香蕉国产在线看| 操出白浆在线播放| 99在线视频只有这里精品首页| 日韩有码中文字幕| 两人在一起打扑克的视频| 一卡2卡三卡四卡精品乱码亚洲| 伦理电影免费视频| 国产精品一区二区免费欧美| 国产1区2区3区精品| 搡老熟女国产l中国老女人| 欧美日韩瑟瑟在线播放| 一个人看的www免费观看视频| 精品国产超薄肉色丝袜足j| 日本免费一区二区三区高清不卡| 午夜福利18| 久久精品夜夜夜夜夜久久蜜豆| 欧美黑人欧美精品刺激| 亚洲av电影在线进入| 日本 欧美在线| 琪琪午夜伦伦电影理论片6080| 亚洲熟妇熟女久久| 久久久水蜜桃国产精品网| 国产精品自产拍在线观看55亚洲| 母亲3免费完整高清在线观看| 高清毛片免费观看视频网站| 在线观看美女被高潮喷水网站 | 非洲黑人性xxxx精品又粗又长| 国产亚洲欧美98| 99久久综合精品五月天人人| 夜夜躁狠狠躁天天躁| 9191精品国产免费久久| 亚洲国产欧美网| 免费在线观看视频国产中文字幕亚洲| 黄色成人免费大全| 日本黄色视频三级网站网址| 日韩人妻高清精品专区| 无遮挡黄片免费观看| 国产一区在线观看成人免费| 久久精品国产亚洲av香蕉五月| 色综合亚洲欧美另类图片| 热99在线观看视频| 黄片小视频在线播放| 成年女人永久免费观看视频| 黄频高清免费视频| 99riav亚洲国产免费| 国产精品美女特级片免费视频播放器 | 久久精品国产亚洲av香蕉五月| 91麻豆精品激情在线观看国产| 十八禁网站免费在线| 制服丝袜大香蕉在线| av视频在线观看入口| 99riav亚洲国产免费| 国产成人一区二区三区免费视频网站| 欧美一级毛片孕妇| 亚洲av电影不卡..在线观看| 亚洲天堂国产精品一区在线| 国产精品电影一区二区三区| 欧美激情久久久久久爽电影| 亚洲成人免费电影在线观看| 欧美一区二区精品小视频在线| 午夜福利18| 美女黄网站色视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲九九香蕉| 一夜夜www| 这个男人来自地球电影免费观看| 午夜精品久久久久久毛片777| 最新在线观看一区二区三区| 色综合欧美亚洲国产小说| 亚洲av成人精品一区久久| 人妻丰满熟妇av一区二区三区| 欧美日本视频| 免费无遮挡裸体视频| 亚洲熟女毛片儿| 国产野战对白在线观看| 在线永久观看黄色视频| 色综合站精品国产| 老鸭窝网址在线观看| 后天国语完整版免费观看| 日本a在线网址| 欧美成人性av电影在线观看| 久久精品综合一区二区三区| 最近视频中文字幕2019在线8| 两个人视频免费观看高清| 亚洲五月婷婷丁香| 国产视频内射| 国产 一区 欧美 日韩| 亚洲午夜理论影院| 一个人看的www免费观看视频| 国产三级中文精品| 国产精品久久久久久精品电影| 国产精品女同一区二区软件 | 熟女人妻精品中文字幕| 亚洲欧美精品综合久久99| 中文字幕久久专区| 搞女人的毛片| 麻豆久久精品国产亚洲av| 国产免费av片在线观看野外av| 亚洲av免费在线观看| cao死你这个sao货| bbb黄色大片| svipshipincom国产片| 国产伦一二天堂av在线观看| 夜夜看夜夜爽夜夜摸| 成人三级黄色视频| 19禁男女啪啪无遮挡网站| 国产成人精品无人区| 在线免费观看不下载黄p国产 | 啦啦啦观看免费观看视频高清| 一个人看的www免费观看视频| 欧美又色又爽又黄视频| 亚洲av免费在线观看| 国产精品久久视频播放| 亚洲欧美激情综合另类| 人妻丰满熟妇av一区二区三区| 一个人看的www免费观看视频| 色在线成人网| 色综合站精品国产| 久久亚洲真实| 女警被强在线播放| 一本精品99久久精品77| 国产伦精品一区二区三区四那| 黄色成人免费大全| 在线免费观看不下载黄p国产 | 亚洲成a人片在线一区二区| 久久国产精品人妻蜜桃| 成人高潮视频无遮挡免费网站| ponron亚洲| 国产又色又爽无遮挡免费看| 成年女人永久免费观看视频| 一本久久中文字幕| 欧美一级毛片孕妇| 欧美乱码精品一区二区三区| 夜夜看夜夜爽夜夜摸| 男人舔女人下体高潮全视频| 露出奶头的视频| 精品一区二区三区视频在线 | 国内久久婷婷六月综合欲色啪| 少妇的逼水好多| 中文亚洲av片在线观看爽| 天天躁狠狠躁夜夜躁狠狠躁| 十八禁网站免费在线| 精品久久久久久久末码| 午夜激情欧美在线| 美女扒开内裤让男人捅视频| 国产亚洲精品综合一区在线观看| 97碰自拍视频| 久久久久久久午夜电影| 国产亚洲av嫩草精品影院| 久久草成人影院| 国产激情偷乱视频一区二区| 国产极品精品免费视频能看的| 宅男免费午夜| cao死你这个sao货| 精品久久蜜臀av无| 国产亚洲精品一区二区www| 亚洲成人久久性| 丰满人妻一区二区三区视频av | 日韩欧美三级三区| cao死你这个sao货| 美女高潮的动态| 人人妻人人看人人澡| 国产一级毛片七仙女欲春2| 久久久久久久久久黄片| 国模一区二区三区四区视频 | 国产男靠女视频免费网站| 久久精品国产清高在天天线| 高清毛片免费观看视频网站| 国产精品 欧美亚洲| 国产亚洲精品久久久久久毛片| 久久99热这里只有精品18| 亚洲欧美一区二区三区黑人| 欧美日韩一级在线毛片| 国产成人av激情在线播放| 国产极品精品免费视频能看的| 亚洲片人在线观看| 麻豆av在线久日| 日本三级黄在线观看| 岛国在线观看网站| 国产三级在线视频| 亚洲九九香蕉| 久久久国产成人精品二区| netflix在线观看网站| 久久欧美精品欧美久久欧美| 长腿黑丝高跟| 香蕉国产在线看| 国产精品一区二区免费欧美| 午夜福利在线观看免费完整高清在 | 国产成人系列免费观看| 最近视频中文字幕2019在线8| 一进一出抽搐动态| 久久久久久久久中文| 精品人妻1区二区| 亚洲中文av在线| 亚洲国产色片| 免费观看精品视频网站| 亚洲熟女毛片儿| 久久精品影院6| 久久九九热精品免费| 脱女人内裤的视频| 色哟哟哟哟哟哟| 亚洲精品色激情综合| 成年人黄色毛片网站| 男女那种视频在线观看| 国产一区二区三区视频了| 亚洲av成人av| 午夜精品在线福利| 日韩欧美国产在线观看| 国产精品一区二区精品视频观看| 国产精品免费一区二区三区在线| 一本一本综合久久| 脱女人内裤的视频| 日本免费一区二区三区高清不卡| 色综合站精品国产| 成人鲁丝片一二三区免费| 成人高潮视频无遮挡免费网站| av欧美777| 日韩高清综合在线| 国产人伦9x9x在线观看| 黑人巨大精品欧美一区二区mp4| 国产精品香港三级国产av潘金莲| 色精品久久人妻99蜜桃| 日本在线视频免费播放| 变态另类成人亚洲欧美熟女| 丰满的人妻完整版| 99riav亚洲国产免费| 真实男女啪啪啪动态图| 亚洲熟女毛片儿| 久久中文看片网| 最近在线观看免费完整版| 一区二区三区国产精品乱码| 日本与韩国留学比较| 91av网站免费观看| 久久伊人香网站| 日本黄色视频三级网站网址| 国产午夜精品久久久久久| 成年免费大片在线观看| www日本在线高清视频| 久久精品国产亚洲av香蕉五月| 最近最新免费中文字幕在线| 久久人妻av系列| 九九在线视频观看精品| 三级男女做爰猛烈吃奶摸视频| 男人舔女人的私密视频| 99热只有精品国产| 夜夜爽天天搞| 欧美国产日韩亚洲一区| 国产一区二区三区视频了| 色综合站精品国产| 国产精品亚洲av一区麻豆| 又黄又粗又硬又大视频| 国产高清激情床上av| 免费看十八禁软件| 欧美大码av| 色哟哟哟哟哟哟| or卡值多少钱| 国产精品亚洲av一区麻豆| 老司机深夜福利视频在线观看| 亚洲国产精品sss在线观看| 精品久久久久久久末码| 久久伊人香网站| 国产一区二区三区视频了| 露出奶头的视频| 99热这里只有精品一区 | 丁香欧美五月| 色播亚洲综合网| 两个人看的免费小视频| 日韩精品青青久久久久久| 岛国在线免费视频观看| 99久久成人亚洲精品观看| 999精品在线视频| 国产精品自产拍在线观看55亚洲| 欧美另类亚洲清纯唯美| 精品一区二区三区视频在线观看免费| 亚洲18禁久久av| 桃色一区二区三区在线观看| 熟妇人妻久久中文字幕3abv| 午夜免费观看网址| 久99久视频精品免费| 无限看片的www在线观看| 国产精品久久久久久精品电影| 午夜精品久久久久久毛片777| 亚洲av片天天在线观看| avwww免费| 黄片大片在线免费观看| 九色国产91popny在线| 中文亚洲av片在线观看爽| 久久久久久大精品| 欧美极品一区二区三区四区| 全区人妻精品视频| 嫩草影院精品99| 真实男女啪啪啪动态图| 亚洲真实伦在线观看| 一本久久中文字幕| 搡老熟女国产l中国老女人| 88av欧美| 99久久精品国产亚洲精品| 国产亚洲av高清不卡| 俄罗斯特黄特色一大片| 亚洲男人的天堂狠狠| 国产高清视频在线播放一区| 黄色视频,在线免费观看| 欧美大码av| 嫩草影视91久久| 久久久久久久精品吃奶| 人妻丰满熟妇av一区二区三区| 人人妻,人人澡人人爽秒播| 亚洲狠狠婷婷综合久久图片| 日本黄色视频三级网站网址| 日本 欧美在线| 看免费av毛片| 国产不卡一卡二| 亚洲人成电影免费在线| 熟妇人妻久久中文字幕3abv| 99热只有精品国产| 性欧美人与动物交配| 国产午夜福利久久久久久| 久久中文字幕一级| 最新美女视频免费是黄的| 日韩人妻高清精品专区| 久久天躁狠狠躁夜夜2o2o| 免费一级毛片在线播放高清视频| 日本 av在线| 在线观看免费午夜福利视频| 欧美精品啪啪一区二区三区| 天堂动漫精品| 老司机午夜十八禁免费视频| 一区二区三区激情视频| 精品免费久久久久久久清纯| 日日夜夜操网爽| 男人舔女人下体高潮全视频| 午夜激情福利司机影院| 一区二区三区激情视频| 亚洲av第一区精品v没综合| 国产aⅴ精品一区二区三区波| 免费看a级黄色片| 高清毛片免费观看视频网站| 老司机深夜福利视频在线观看| 亚洲国产欧洲综合997久久,| 中出人妻视频一区二区| 国产高清视频在线播放一区| 欧美zozozo另类| 美女被艹到高潮喷水动态| 日本 av在线| 免费在线观看亚洲国产| 九九在线视频观看精品| а√天堂www在线а√下载| 亚洲av成人精品一区久久| 在线视频色国产色| 欧美色视频一区免费| 亚洲午夜精品一区,二区,三区| 成年女人永久免费观看视频| 在线播放国产精品三级| 日本五十路高清| 欧美成狂野欧美在线观看| 国产免费男女视频| 国产精品九九99| 一本精品99久久精品77| 色在线成人网| 99精品久久久久人妻精品| 国产乱人伦免费视频| 亚洲中文av在线| 日韩欧美国产在线观看| 国产三级在线视频| 好男人在线观看高清免费视频| 国语自产精品视频在线第100页| 综合色av麻豆| svipshipincom国产片| 国产av在哪里看| 欧美色视频一区免费| 国产精品1区2区在线观看.| 免费观看人在逋| 久久久久久久久中文| 国产精品一区二区三区四区久久| 国产激情久久老熟女| 久久久久久国产a免费观看| 国产精品 国内视频| 超碰成人久久| 噜噜噜噜噜久久久久久91| 国产精品一区二区免费欧美| 国产成人系列免费观看| 精品一区二区三区四区五区乱码| 日韩中文字幕欧美一区二区| 老鸭窝网址在线观看| 国产精品99久久久久久久久| 国产欧美日韩一区二区三| 国产成人aa在线观看| 搡老熟女国产l中国老女人| 久久人人精品亚洲av| 欧美性猛交╳xxx乱大交人| 亚洲欧美激情综合另类| 国产极品精品免费视频能看的| a级毛片在线看网站| 婷婷亚洲欧美| 一级毛片高清免费大全| 美女免费视频网站| 99久久精品国产亚洲精品| 三级国产精品欧美在线观看 | 欧美日韩福利视频一区二区| 欧美成人一区二区免费高清观看 | 麻豆一二三区av精品| ponron亚洲| 欧美性猛交黑人性爽| 很黄的视频免费| 国产亚洲精品久久久久久毛片| 国产精品精品国产色婷婷| 日韩大尺度精品在线看网址| 久久精品夜夜夜夜夜久久蜜豆| 日本黄色片子视频| 老熟妇仑乱视频hdxx| 高清毛片免费观看视频网站| 亚洲国产精品合色在线| 国产精品av久久久久免费| 99久久精品一区二区三区| 岛国视频午夜一区免费看| 蜜桃久久精品国产亚洲av| www日本黄色视频网| 欧美色欧美亚洲另类二区| 亚洲专区中文字幕在线| 后天国语完整版免费观看| 全区人妻精品视频| 麻豆av在线久日| 国产精品一区二区三区四区免费观看 | 成年免费大片在线观看| 午夜福利高清视频| 后天国语完整版免费观看| 精品国产亚洲在线| 俺也久久电影网| 精品国产亚洲在线| 亚洲五月婷婷丁香| 免费无遮挡裸体视频| 99国产精品一区二区蜜桃av| 老司机深夜福利视频在线观看| 99国产综合亚洲精品| 日韩欧美在线乱码| 亚洲无线在线观看| 看黄色毛片网站| 中文字幕熟女人妻在线| 天堂av国产一区二区熟女人妻| 亚洲精品美女久久久久99蜜臀| 久久久久精品国产欧美久久久| 97人妻精品一区二区三区麻豆| 在线看三级毛片| 日本与韩国留学比较| 午夜激情欧美在线| 欧美三级亚洲精品| 男人和女人高潮做爰伦理| 国产极品精品免费视频能看的| 免费观看的影片在线观看| 熟女电影av网| 桃红色精品国产亚洲av| 99视频精品全部免费 在线 | 又黄又爽又免费观看的视频| 亚洲九九香蕉| 久久久久免费精品人妻一区二区| 俄罗斯特黄特色一大片| 中文字幕高清在线视频| 国产精品日韩av在线免费观看| 99久久精品一区二区三区| 亚洲精品粉嫩美女一区| 日日夜夜操网爽| 亚洲精品美女久久av网站| 国产探花在线观看一区二区| 亚洲熟妇中文字幕五十中出| 高潮久久久久久久久久久不卡| 十八禁网站免费在线| 脱女人内裤的视频| 琪琪午夜伦伦电影理论片6080| 国产av一区在线观看免费| 日韩免费av在线播放| 亚洲国产色片| 国产极品精品免费视频能看的| 18禁美女被吸乳视频| 美女被艹到高潮喷水动态| 国产午夜精品久久久久久| 99久久精品一区二区三区| 淫秽高清视频在线观看| 国产午夜精品论理片| 观看美女的网站| 国产精品乱码一区二三区的特点| 两个人视频免费观看高清| 网址你懂的国产日韩在线| 18美女黄网站色大片免费观看| 精品不卡国产一区二区三区| 最近最新免费中文字幕在线| 全区人妻精品视频| 别揉我奶头~嗯~啊~动态视频| 国产一区二区三区视频了| 国产真实乱freesex| 国产三级中文精品| 美女高潮的动态| 国产不卡一卡二| 丰满人妻熟妇乱又伦精品不卡| 99热只有精品国产| 国产激情欧美一区二区| 亚洲熟妇熟女久久| 日本一本二区三区精品| 久久久水蜜桃国产精品网| 少妇的逼水好多| 久久中文字幕人妻熟女| 国产熟女xx| 亚洲欧美一区二区三区黑人| 老司机在亚洲福利影院| 中国美女看黄片| 亚洲 欧美 日韩 在线 免费| 国产成人aa在线观看| 午夜影院日韩av| 人妻久久中文字幕网| 久久久成人免费电影| 久久人妻av系列| 日韩三级视频一区二区三区| 人人妻人人看人人澡| 偷拍熟女少妇极品色| 757午夜福利合集在线观看| 亚洲五月天丁香| 国产一区二区三区视频了| 国产麻豆成人av免费视频| 国产精品综合久久久久久久免费| 99久久国产精品久久久| 人人妻人人澡欧美一区二区| 99久久综合精品五月天人人| 草草在线视频免费看| 两个人看的免费小视频| 欧美一区二区国产精品久久精品| 在线观看舔阴道视频| 99国产精品一区二区蜜桃av| 日韩中文字幕欧美一区二区| 欧美色视频一区免费| 亚洲精品在线观看二区| 中出人妻视频一区二区| 国产精品av视频在线免费观看| 99热这里只有是精品50| 九九热线精品视视频播放| 免费看光身美女| 日韩高清综合在线| 久久久久久九九精品二区国产| 桃色一区二区三区在线观看| 国产又色又爽无遮挡免费看|