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

    基于IHB法的轉(zhuǎn)子系統(tǒng)松動故障特性研究*

    2014-08-22 11:23:08曹青松熊國良
    機(jī)床與液壓 2014年24期
    關(guān)鍵詞:青松華東工程學(xué)院

    曹青松,向 琴,熊國良

    華東交通大學(xué)機(jī)電工程學(xué)院,南昌 330013

    1.Introduction

    Mechanical fault diagnosis technology,which is usually used to exploit the running states of equipment,includes fault detection,fault identification and fault isolation.Through this diagnosis technology,the performance of equipment could be evaluated as normal or not and the initial cause of failure might be forecast.Since fault characteristics are the basis of fault identification,the prerequisite of fault diagnosis is a kind of in-depth research on fault characteristics.Pedestal looseness is a common fault in a rotor-bearing system and it will reduce holistic stiffness,mechanical damping.All these results could cause violent vibration for the whole system.Therefore,the study of pedestal looseness is very significant for real engineering applications in terms of the safe operation of rotating machinery,the extension of service life,and the improvement of work efficiency.

    In the 90’s,Muszynska and Goldman devoted themselves to the research on dynamics of rotating machinery,and made some contribution to this field.They proposed the bilinear model of a rotor system with one loose pedestal.Their numeric results showed that the vibration response included synchronous and subsynchronous fractional components,even leaded to chaos[1].At present,there are two methods to study the looseness fault characteristics,i.e.,experiment-al method and dynamic modeling method.The former one is limited due to the high cost,so the latter one is adopted widely by domestic and oversea scholars.Dynamic models mainly contain the nonlinear differential equations[2]and finite element model[3].Nonlinear differential equations could be usually solved by using Matlab software and numerical integration method,and the finite element model could be solved based on the Ansys software.The integration methods such as Runge-Kutta and Newmark-βmethod are usually adopted to solve the problem of looseness fault[4].Although the available researches have build the foundation for the looseness fault diagnosis,the looseness models are mainly based on the piecewise-linear model that can not describe the real phenomena of looseness due to the unclear fault mechanism.Besides,most of the available results are limited to graphs and seldom achieve the periodic solution of vibration response.This drawback makes it difficult to distinguish the frequency components and the calculation is a time consuming process.

    IHB is a kind of method which combines the incremental method and harmonic balance method to solve the nonlinear vibration equations[5].Since this method has some advantages such as calculation time-saving,easy to converge and appropriate to a strong non-linear system,it is usually applied to solve strong non-linear problems such as piecewiselinear system[6],Duffing oscillator system[7]and mistuning bladed disks system[8].

    Due to the fact that the vibration equation of rotor system with pedestal looseness and unbalanced force is a strong non-linear and IHB method has excellent performance to solve non-linear equation,this method is introduced into the solution of rotor system with looseness fault to acquire the periodic solution for the vibration response and analyze the fault characteristics of pedestal looseness.The results of present study would provide some theoretical basis for the diagnosis of looseness fault system.

    2.Looseness model

    In this paper,a sample single-disk Jeffcott system was considered.As shown in Figure1,one end of the shaft is loosed.Assume the looseness clearance is far less than the vibration displacement of the system,and then the imbalance caused by looseness clearance could be ignored.

    1)Traditional looseness model.When the vibration amplitude of shaft end is less than or equal to looseness clearance α,stiffness of the system,i.e.,k1,is small.When the vibration amplitude of shaft end is greater than the looseness clearanceα,stiffness of system will be changed to k2due to the impact force from the constrained interface. Therefore,looseness stiffness should be a piecewise-linear function as follows:

    2)Continuous looseness model.The looseness model shown in Eq.(1)is relatively idealized.However,the real phenomena of looseness stiffness could not be such step change.It will be continuously and non-linearly changed with the vibration amplitude of shaft end.Since the looseness mechanism was not clear and the complex change of stiffness was hard to describe accurately,a simple continuous model could be proposed.If the vibration amplitude of shaft end is less than or equal to the looseness clearanceα,the stiffness of the system will be linearly changed.However,if the vibration amplitude of shaft end is greater than the looseness clearanceα,the stiffness of system will be a constant value which is equal to the stiffness of shaft.Therefore,the continuous looseness model could be expressed as follows:

    Then,the differential equation of motion for Jeffcott rotor system will be written as follows:

    Where,m,c,k,k1,k0,F(xiàn)(t)and x denote the total mass of system,viscous damping coefficient,holistic stiffness,stiffness of shaft,looseness stiffness,exciting force of interference outside resulted from eccentric mass and vibration displacement,respectively.

    Figure 1.Schematic diagram of a Jeffcott rotor system with looseness fault

    3.IHB method solution

    IHB method was adopted to solve the Eq.(3).If kx=k0x+f(x),f(x)can be written as follows:

    For convenience,we introduce a phase angleτ given by

    and then the independent variable could be changed from t to τ.Therefore,Eq.(3)could be transformed into the following equation.

    Assume that the current state of vibration corresponding to the excitation parametersω0and F0is denoted by x0(τ),a neighboring state could be expressed by adding the corresponding increments as follows:

    For a small incrementΔx,the function f(x)may be written approximately as follows:

    Where,f'( x0) represents the value of the first derivative of f(x)with respect to x at x0which could be expressed by using the step function.Put the Eq.(7)and Eq.(8)into Eq.(6),and consider the parametersΔω andΔF are equal to zero and ignore the small value of higher order,then a linearized incremental equation could be obtained as the following form:

    Where,

    R is the correction term which will equal to zero if the values of ω0、F0and x0are the precise solutions of Eq.(6).

    The approximate periodic solution x0may be expressed as

    Accordingly,the increment is as follows:

    Where,

    and N is the number of harmonic terms.

    The Galerkin procedure is carried out with the Δan’s and Δbn’s as the generalized coordinate

    Eq.(11)and Eq.(12)are substituted into Eq.(16),one could obtain Eq.(17)after the rearrangement

    Where,

    The explicit expressions for Kij,Ri,Piand ΔQi(i=1,2;j=1,2)could be worked out as follows:

    ①Elements of the matrix K:

    Where,

    ②Elements of the matrix R:

    The nonlinear part of the elements of K and R are expressed by

    4.Simulation analysis

    The initial state of system A0andω0are given to iterate in the program with the following equation

    The updated matrices K(i+1),R(i+1)are then evaluated from Eqs.(18) ~ (25)in terms of the(i+1)th amplitude vector A(i+1)for the next iteration.In order to obtain the next point of solution,the iteration process will be proceeded until the criteria(R)≤ε is satisfied,whereεis the permissible error according to the given accuracy.A will be obtained when the iteration process is stopped,which is the coefficient of the harmonic terms.

    Parameters of the Jeffcott rotor system are as follows:total mass of the system m=3 kg,eccentric mass Δm=5 g,eccentric displacement e=30.8 mm,system span L=630 mm,diameter of shaft d=15.4 mm,and modulus of elasticity E=206 GPa.The stiffness of shaft k2is 1.081 ×105N/m when the system is normal.

    Since there exist only one unbalanced force from eccentric rotor,one pedestal looseness was considered.In this paper,the fault characteristics of rotor system with looseness was studied only in the vertical direction since looseness fault has directivity.Assume the rotate speed ω0is 1 400 r/min or 147 rad/s,looseness clearance α is 5 μm and the number of harmonic N is 12.

    The periodic solution of the rotor system without looseness fault could be obtained by the IHB method,the time domain waveform,amplitude spectrogram,phase diagram and Poincare sectional drawing were shown in Figure 2.Figure 2(a)is the time domain waveform,which is smooth and the vibration displacement is approximate 34μm.From Figure 2(b)one could understand that the component of the vibration frequency is one fold and equal to the rotation frequency of rotor.The rotor motion is periodic since the phase diagram is a closed curve and there is only one point on the Poincare sectional drawing.

    Figure 2.The time domain waveform,amplitude spectrogram,phase diagram and Poincare sectional drawing of the system without looseness

    If the looseness stiffness k0is 3.6 ×104N/m when the fault happened,the periodic solution of the rotor system with looseness fault could be obtained by using the IHB method,the time domain waveform,amplitude spectrogram,phase diagram and Poincare sectional drawing were shown in Figure 3.In this figure,the time domain waveform became abnormal and the maximum vibration displacement was increased to 42 μm.Furthermore,the vibration response has complicated frequency components which consist of 1/3X,3/3X,5/3X,7/3X(X indicates the rotation frequency of rotor)and so on.The phase diagram is still a closed curve but more complicated.The rotor motion was judged to be period-three from the three points on the Poincare sectional drawing.

    If the looseness stiffness k0is 0.6 ×104N/m when the fault happened,the same work could be done using the same method.Time domain waveform,amplitude spectrogram,phase diagram and Poincare sectional drawing were shown in Figure 4,respectively.The time domain waveform is messier than the previous one and the vibration displacement got increased apparently.The frequency components of vibration response became more adjacent which consist of 1/3X,2/3X,3/3X,4/3X and so on.It proved the importance of the Poincare sectional drawing for the phase diagram as shown in Figure 4(c)was too disordered to distinguish the motion state of rotor.The rotor motion is period-three,as shown in Figure 4(d).

    Figure 3.The time domain waveform,amplitude spectrogram,phase diagram and the Poincare sectional drawing of the system when the looseness stiffness k0=3.6 ×104 N/m

    Figure 4.The time domain waveform,amplitude spectrogram,phase diagram and the Poincare sectional drawing of the system when the looseness stiffness k0=0.6 ×104 N/m

    The simulation results showed that not only the vibration amplitude of the failing system is greater than that of the normal one,but also the frequency components become more complex.The phase diagram and the Poincare sectional drawing became more disperse with the decrease of looseness stiffness.All these results almost correspond to the result of reference[4].It also could explain that the IHB method is capable to solve this issue like rotor system with looseness fault,and it also proves the validity of the continuous looseness model.

    5.Conclusion

    In this paper,the nonlinear dynamics and fault characteristics of a Jeffcott rotor system with pedestal looseness was studied.IHB method was exploited to solve a dynamics model with a continuous changed stiffness.Some conclusions could be drawn as follows:

    1)The time domain waveform of vibration response for rotor system is smooth when there was no looseness fault,the vibration amplitude is relatively small,and there exists only 1X amplitude spectrogram,the rotor motion is period-one.

    2)When the looseness fault happened,the time domain waveform of vibration response became messy,the vibration amplitude gets increased,and the degree of mess depends on the reduction of the looseness stiffness.The rotor motion became period-three.

    3)The results also confirmed the feasibility that the IHB method could be applied to solve the issue of looseness system.The new continuous model built some foundation for nonlinear continuous looseness model which could describe the looseness mechanism best.And this research work could provide a theoretical basis for the diagnosis of looseness fault.

    [1] Muszynska A,Goldman P.Chaotic responses of unbalanced rotor/bearing/stator systems with looseness or rubs[J].Chaos,Solitons & Fractals,1995,5(9):1683-1704.

    [2] Goldman P,Muszynska A.Dynamic effects in mechanical structures with gaps and impacting:order and chaos[J].Journal of Vibration and Acoustics,1994,116(4):541-547.

    [3] LIU Xiandong,HE Tian,LI Qihan.Dynamic model of rotor system with support loosening and its diagnosis method[J].Journal of Aerospace Power,2005,20(1):54-59.

    [4] MA Hui,SUN Wei,Xuejun Wang,et al.Characteristic analysis of looseness fault in rotor system[J].Journal of Northeastern University,2009,30(3):400-404.

    [5] LIU Yuanfeng,ZHAO Mei,ZHU Houjun.Chaotic Behavior of a rotor under rub and looseness[J].Journal of Vibration and Shock,2002,21(4):36-38,21.

    [6] YANG Yongfeng,REN Xingmin,QIN Wei-yang.Nonlinear Dynamics Behaviors of Rotor Systems with Crack and Pedestal Looseness Faults[J].Mechanical Science and Technology,2005,24(8):985-987.

    [7] CHEN Shuhui.Quantitative analysis method of strongly nonlinear vibration system[M].Beijing:Science Press,2006.

    [8] Lau SL,ZHANGWS.Nonlinear vibrations of piecewise linear systems by incremental harmonic balance method,Journal of Applied Mechanics,1992,59:153-160.

    [9] Leung A Y T,CHUI SK.Nonlinear vibration of coupled duffing oscillators by an improved incremental harmonic balance method[J].Journal of Sound and Vibration,1995,181(4):619-633.

    [10] WANG Ailun,LONG Qing.Application of IHB method to study the response characteristics of mistuning bladed disks with strong non-linear friction[J].Journal of Mechanical Strength,2012,34(2):159-164.

    猜你喜歡
    青松華東工程學(xué)院
    福建工程學(xué)院
    華東銷售在一線
    福建工程學(xué)院
    Experimental Study on Grinding Force of Electrostatic Coated Grinding Wheel
    相華東:走在欣欣向榮的田野上
    華人時刊(2022年21期)2022-02-15 03:42:36
    豐收歌
    福建工程學(xué)院
    林青松
    寶藏(2018年6期)2018-07-10 02:26:44
    福建工程學(xué)院
    多絲量新品種華東×春晨的引進(jìn)推廣
    蠶桑通報(2015年2期)2015-12-15 00:41:56
    国产在线一区二区三区精| 在线观看免费高清a一片| 久久久亚洲精品成人影院| 亚洲精华国产精华液的使用体验| 亚洲国产成人一精品久久久| a级毛片黄视频| 日日爽夜夜爽网站| 亚洲欧美精品自产自拍| av在线老鸭窝| 中文字幕制服av| av女优亚洲男人天堂| 亚洲欧洲日产国产| 五月伊人婷婷丁香| 永久免费av网站大全| 国产精品不卡视频一区二区| 久久精品国产亚洲av天美| 国产精品久久久久久精品电影小说| 一级片'在线观看视频| 免费观看a级毛片全部| 欧美激情高清一区二区三区 | 国产伦理片在线播放av一区| 香蕉精品网在线| 国产老妇伦熟女老妇高清| 精品少妇久久久久久888优播| 亚洲综合色惰| 汤姆久久久久久久影院中文字幕| 国产精品av久久久久免费| 午夜福利视频在线观看免费| 成人国产av品久久久| 中文精品一卡2卡3卡4更新| 在线看a的网站| 一区二区日韩欧美中文字幕| 丝袜喷水一区| 亚洲图色成人| 美女中出高潮动态图| 国产乱来视频区| 超色免费av| 色网站视频免费| 青春草视频在线免费观看| 美国免费a级毛片| 久久久久久久久免费视频了| 可以免费在线观看a视频的电影网站 | 伦理电影免费视频| 国产成人免费观看mmmm| 久久99精品国语久久久| 色网站视频免费| 午夜福利一区二区在线看| 亚洲少妇的诱惑av| 最近的中文字幕免费完整| 国产精品久久久久久av不卡| 男女高潮啪啪啪动态图| 精品一区二区三区四区五区乱码 | 国产精品免费大片| 少妇人妻 视频| 一级片免费观看大全| 成人免费观看视频高清| 热99久久久久精品小说推荐| 国产日韩欧美亚洲二区| 久久久精品区二区三区| 99久久精品国产国产毛片| 综合色丁香网| 国产精品欧美亚洲77777| 性少妇av在线| 亚洲人成电影观看| 久久99一区二区三区| 90打野战视频偷拍视频| 91精品三级在线观看| 久久精品aⅴ一区二区三区四区 | 欧美人与善性xxx| 亚洲精品久久成人aⅴ小说| 国产日韩欧美亚洲二区| 纯流量卡能插随身wifi吗| 777久久人妻少妇嫩草av网站| 亚洲婷婷狠狠爱综合网| 国产成人av激情在线播放| www.精华液| 国产亚洲欧美精品永久| 免费在线观看黄色视频的| 只有这里有精品99| 纵有疾风起免费观看全集完整版| 成人亚洲精品一区在线观看| 宅男免费午夜| 另类亚洲欧美激情| 国产日韩欧美亚洲二区| 国产精品久久久久久久久免| 欧美 日韩 精品 国产| 黄网站色视频无遮挡免费观看| 人妻人人澡人人爽人人| 久久韩国三级中文字幕| 日本爱情动作片www.在线观看| 一级毛片电影观看| 可以免费在线观看a视频的电影网站 | 99国产精品免费福利视频| 成人午夜精彩视频在线观看| 国产亚洲精品第一综合不卡| 亚洲精品国产av蜜桃| 日韩中文字幕欧美一区二区 | 亚洲av国产av综合av卡| 国产精品成人在线| 精品国产乱码久久久久久小说| 日韩av在线免费看完整版不卡| 日本-黄色视频高清免费观看| 亚洲人成77777在线视频| 欧美人与善性xxx| 99re6热这里在线精品视频| 亚洲欧美成人精品一区二区| 国产av一区二区精品久久| 国产淫语在线视频| tube8黄色片| 亚洲国产精品国产精品| 如何舔出高潮| 国产精品久久久久久久久免| 麻豆精品久久久久久蜜桃| 国产黄色视频一区二区在线观看| 一区二区三区乱码不卡18| 亚洲国产最新在线播放| 少妇猛男粗大的猛烈进出视频| 女性被躁到高潮视频| 色播在线永久视频| 日韩成人av中文字幕在线观看| 我的亚洲天堂| 国产精品女同一区二区软件| 国产亚洲最大av| av又黄又爽大尺度在线免费看| 亚洲成国产人片在线观看| 日韩中文字幕视频在线看片| 精品一区二区三区四区五区乱码 | 中文字幕最新亚洲高清| 欧美精品av麻豆av| 欧美国产精品一级二级三级| av免费在线看不卡| 亚洲欧美成人综合另类久久久| 99香蕉大伊视频| 亚洲美女黄色视频免费看| 制服诱惑二区| 美女主播在线视频| 日本色播在线视频| 国产欧美日韩综合在线一区二区| 一区在线观看完整版| 男人操女人黄网站| 一区在线观看完整版| 午夜av观看不卡| 91国产中文字幕| 制服诱惑二区| 26uuu在线亚洲综合色| 国产欧美日韩综合在线一区二区| 国产极品粉嫩免费观看在线| 国产熟女午夜一区二区三区| av免费观看日本| 日韩精品有码人妻一区| 美女国产视频在线观看| av片东京热男人的天堂| 成人国产av品久久久| 亚洲欧美精品综合一区二区三区 | 日日啪夜夜爽| 极品少妇高潮喷水抽搐| 国产亚洲av片在线观看秒播厂| 极品少妇高潮喷水抽搐| 美女午夜性视频免费| 秋霞伦理黄片| 亚洲精品久久成人aⅴ小说| 巨乳人妻的诱惑在线观看| 成年女人在线观看亚洲视频| 男人添女人高潮全过程视频| 97在线人人人人妻| 少妇精品久久久久久久| 日本91视频免费播放| 久久久久久久亚洲中文字幕| 91成人精品电影| 国产xxxxx性猛交| 欧美97在线视频| 日韩电影二区| 日韩电影二区| 国产淫语在线视频| 精品人妻在线不人妻| 18+在线观看网站| 国产精品二区激情视频| 久久精品国产亚洲av涩爱| 成人亚洲欧美一区二区av| 99热网站在线观看| 美女国产视频在线观看| 久久久久人妻精品一区果冻| 久久午夜综合久久蜜桃| 国产亚洲午夜精品一区二区久久| 美女xxoo啪啪120秒动态图| 国产av国产精品国产| 婷婷色综合www| 九色亚洲精品在线播放| 国产成人欧美| 欧美在线黄色| 国产精品不卡视频一区二区| 老司机影院毛片| 亚洲精品日韩在线中文字幕| 99热网站在线观看| 国产片内射在线| 美女xxoo啪啪120秒动态图| 国产综合精华液| 美女福利国产在线| 久久久久网色| 久久久久久久大尺度免费视频| 亚洲三级黄色毛片| 久久久久精品人妻al黑| av.在线天堂| 免费人妻精品一区二区三区视频| 日韩大片免费观看网站| 在线观看国产h片| 大片免费播放器 马上看| 亚洲精品日本国产第一区| 在线观看www视频免费| 久久精品国产鲁丝片午夜精品| 男女无遮挡免费网站观看| 亚洲在久久综合| 大陆偷拍与自拍| 两性夫妻黄色片| 一区二区av电影网| 最新的欧美精品一区二区| 99久久中文字幕三级久久日本| 精品第一国产精品| av不卡在线播放| www.熟女人妻精品国产| 另类精品久久| 不卡av一区二区三区| 美女视频免费永久观看网站| 男人添女人高潮全过程视频| 久久女婷五月综合色啪小说| 亚洲精品视频女| 精品国产一区二区久久| 男女边吃奶边做爰视频| 嫩草影院入口| 国产精品国产三级国产专区5o| 久久精品国产a三级三级三级| av视频免费观看在线观看| 免费女性裸体啪啪无遮挡网站| 这个男人来自地球电影免费观看 | 90打野战视频偷拍视频| 日韩av在线免费看完整版不卡| 日韩大片免费观看网站| av国产精品久久久久影院| 五月天丁香电影| 亚洲伊人色综图| 黑丝袜美女国产一区| 少妇人妻 视频| 国产高清国产精品国产三级| 不卡av一区二区三区| 午夜精品国产一区二区电影| 少妇的逼水好多| 亚洲第一区二区三区不卡| 久久99蜜桃精品久久| 国产成人精品婷婷| av国产精品久久久久影院| 成年人午夜在线观看视频| 涩涩av久久男人的天堂| 亚洲精品美女久久av网站| a级毛片在线看网站| 国产麻豆69| 男男h啪啪无遮挡| 婷婷色av中文字幕| 高清视频免费观看一区二区| 久久99一区二区三区| 亚洲天堂av无毛| 亚洲欧洲日产国产| 我要看黄色一级片免费的| 男女午夜视频在线观看| 国产免费现黄频在线看| 制服诱惑二区| 日本猛色少妇xxxxx猛交久久| videossex国产| 久久久久久久国产电影| 亚洲欧美日韩另类电影网站| 亚洲av免费高清在线观看| 欧美亚洲日本最大视频资源| 日韩欧美一区视频在线观看| 在线天堂最新版资源| 少妇人妻 视频| 校园人妻丝袜中文字幕| 亚洲欧美成人综合另类久久久| 91久久精品国产一区二区三区| 午夜激情av网站| 亚洲第一av免费看| 日韩熟女老妇一区二区性免费视频| 这个男人来自地球电影免费观看 | 热99久久久久精品小说推荐| 人体艺术视频欧美日本| 色婷婷av一区二区三区视频| 国产在视频线精品| 国产爽快片一区二区三区| www.av在线官网国产| 超色免费av| 97在线视频观看| av国产久精品久网站免费入址| 少妇熟女欧美另类| 熟妇人妻不卡中文字幕| 日韩欧美一区视频在线观看| 中文字幕人妻丝袜制服| 国产男人的电影天堂91| 热re99久久国产66热| 十八禁高潮呻吟视频| 99精国产麻豆久久婷婷| 久久久久精品久久久久真实原创| 嫩草影院入口| 18禁观看日本| 国产一区二区三区综合在线观看| 这个男人来自地球电影免费观看 | 日韩不卡一区二区三区视频在线| 亚洲色图综合在线观看| 黄网站色视频无遮挡免费观看| 久久av网站| 亚洲国产精品国产精品| 精品人妻熟女毛片av久久网站| 亚洲精品,欧美精品| 久久影院123| 丝袜脚勾引网站| 精品一区二区三卡| 国产熟女欧美一区二区| 亚洲美女视频黄频| 国产成人a∨麻豆精品| 麻豆精品久久久久久蜜桃| 亚洲男人天堂网一区| 欧美日韩精品成人综合77777| 多毛熟女@视频| 看免费av毛片| 九色亚洲精品在线播放| 国产免费又黄又爽又色| 久久精品国产亚洲av高清一级| 亚洲av.av天堂| 春色校园在线视频观看| 日本91视频免费播放| 大陆偷拍与自拍| 亚洲av.av天堂| 免费在线观看视频国产中文字幕亚洲 | videosex国产| 久久久精品国产亚洲av高清涩受| 成年女人毛片免费观看观看9 | 亚洲欧美一区二区三区黑人 | 亚洲综合精品二区| 久久精品国产a三级三级三级| 亚洲av电影在线进入| 街头女战士在线观看网站| 如日韩欧美国产精品一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 大香蕉久久网| 亚洲av成人精品一二三区| 91精品国产国语对白视频| 久久精品国产鲁丝片午夜精品| 国产在线免费精品| 又大又黄又爽视频免费| 久久国产亚洲av麻豆专区| 久久精品国产亚洲av天美| 十八禁网站网址无遮挡| 国产 一区精品| 两个人看的免费小视频| 国产精品一二三区在线看| 久久狼人影院| 一级片免费观看大全| av一本久久久久| 亚洲欧美精品综合一区二区三区 | 建设人人有责人人尽责人人享有的| 亚洲三级黄色毛片| 亚洲国产精品999| 91国产中文字幕| 日韩成人av中文字幕在线观看| 久热久热在线精品观看| 99国产精品免费福利视频| 亚洲色图 男人天堂 中文字幕| 搡老乐熟女国产| 妹子高潮喷水视频| 欧美最新免费一区二区三区| 国产精品亚洲av一区麻豆 | 亚洲欧美一区二区三区久久| 日日啪夜夜爽| 日韩熟女老妇一区二区性免费视频| 国产精品偷伦视频观看了| 亚洲 欧美一区二区三区| 飞空精品影院首页| 色94色欧美一区二区| 一本色道久久久久久精品综合| 女性被躁到高潮视频| a级片在线免费高清观看视频| 久久99精品国语久久久| 少妇猛男粗大的猛烈进出视频| 涩涩av久久男人的天堂| 男女免费视频国产| 熟女少妇亚洲综合色aaa.| 日韩 亚洲 欧美在线| 亚洲第一青青草原| 桃花免费在线播放| 热re99久久精品国产66热6| 男的添女的下面高潮视频| 午夜日韩欧美国产| 欧美人与善性xxx| 黑人猛操日本美女一级片| 91aial.com中文字幕在线观看| 亚洲婷婷狠狠爱综合网| 国产淫语在线视频| 国产精品一国产av| 色婷婷久久久亚洲欧美| 不卡视频在线观看欧美| 国产精品亚洲av一区麻豆 | 国语对白做爰xxxⅹ性视频网站| 18+在线观看网站| 五月伊人婷婷丁香| 久久国内精品自在自线图片| 国产av一区二区精品久久| 欧美日韩视频精品一区| 黄色毛片三级朝国网站| 免费少妇av软件| 91国产中文字幕| 9191精品国产免费久久| 午夜影院在线不卡| 汤姆久久久久久久影院中文字幕| 啦啦啦中文免费视频观看日本| 欧美 日韩 精品 国产| 久久精品国产亚洲av高清一级| 多毛熟女@视频| 婷婷色av中文字幕| 最近中文字幕2019免费版| 国产精品欧美亚洲77777| 日韩一区二区三区影片| 久久女婷五月综合色啪小说| 一边亲一边摸免费视频| 另类精品久久| 国产男女超爽视频在线观看| 色哟哟·www| 18禁动态无遮挡网站| 你懂的网址亚洲精品在线观看| 成人漫画全彩无遮挡| 亚洲伊人色综图| 1024视频免费在线观看| 久久97久久精品| 亚洲国产看品久久| 观看av在线不卡| 日韩欧美精品免费久久| 人妻一区二区av| 日韩中文字幕视频在线看片| 亚洲三区欧美一区| 夫妻性生交免费视频一级片| √禁漫天堂资源中文www| 999精品在线视频| 成年美女黄网站色视频大全免费| 久久久久人妻精品一区果冻| 91成人精品电影| 精品人妻偷拍中文字幕| 午夜福利视频精品| 国产人伦9x9x在线观看 | 亚洲精品久久成人aⅴ小说| 亚洲人成电影观看| 国语对白做爰xxxⅹ性视频网站| 赤兔流量卡办理| 日韩伦理黄色片| 亚洲国产精品成人久久小说| 久久久久精品人妻al黑| 国产精品免费大片| 亚洲国产精品国产精品| 久久99精品国语久久久| 国产精品熟女久久久久浪| 国产欧美亚洲国产| 午夜久久久在线观看| 伊人久久大香线蕉亚洲五| 日韩av不卡免费在线播放| 欧美日韩国产mv在线观看视频| 久久亚洲国产成人精品v| 爱豆传媒免费全集在线观看| 国产欧美日韩一区二区三区在线| 国产成人精品无人区| 美女福利国产在线| 亚洲精华国产精华液的使用体验| 日本vs欧美在线观看视频| 三级国产精品片| 免费观看在线日韩| 国精品久久久久久国模美| 两个人免费观看高清视频| 亚洲欧美成人综合另类久久久| 国产黄色视频一区二区在线观看| 成人亚洲欧美一区二区av| 国产成人av激情在线播放| 亚洲精品久久成人aⅴ小说| 亚洲伊人久久精品综合| 国产精品一二三区在线看| 亚洲av日韩在线播放| 我的亚洲天堂| 国产精品香港三级国产av潘金莲 | 精品人妻一区二区三区麻豆| 成人手机av| 曰老女人黄片| 国产成人一区二区在线| 国产麻豆69| 热re99久久精品国产66热6| 男女下面插进去视频免费观看| 色播在线永久视频| 免费在线观看完整版高清| 国产精品人妻久久久影院| 热re99久久精品国产66热6| 好男人视频免费观看在线| 青春草视频在线免费观看| 日本wwww免费看| 国产精品久久久久久久久免| 青草久久国产| 麻豆乱淫一区二区| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品久久午夜乱码| 美女午夜性视频免费| 日韩精品有码人妻一区| 天堂俺去俺来也www色官网| 国产女主播在线喷水免费视频网站| 亚洲精品国产av成人精品| 黄色配什么色好看| 日韩成人av中文字幕在线观看| 久久毛片免费看一区二区三区| 91久久精品国产一区二区三区| 国产一区有黄有色的免费视频| 亚洲国产av影院在线观看| 免费黄色在线免费观看| 亚洲av欧美aⅴ国产| 99久国产av精品国产电影| 亚洲av电影在线进入| 国产免费视频播放在线视频| 久久精品国产亚洲av涩爱| 欧美精品av麻豆av| 男女无遮挡免费网站观看| 搡女人真爽免费视频火全软件| 香蕉丝袜av| 欧美精品国产亚洲| 国产亚洲精品第一综合不卡| 亚洲欧洲国产日韩| 亚洲av综合色区一区| 青草久久国产| 麻豆精品久久久久久蜜桃| 国产极品天堂在线| 国产亚洲一区二区精品| 久久免费观看电影| 亚洲精品第二区| 1024视频免费在线观看| 人妻少妇偷人精品九色| 免费女性裸体啪啪无遮挡网站| 亚洲av福利一区| 一个人免费看片子| 亚洲av男天堂| 午夜精品国产一区二区电影| a级毛片在线看网站| 中文乱码字字幕精品一区二区三区| 精品国产露脸久久av麻豆| 亚洲成人手机| 岛国毛片在线播放| 日韩成人av中文字幕在线观看| 亚洲视频免费观看视频| 男女边吃奶边做爰视频| 少妇精品久久久久久久| av免费在线看不卡| 国产乱人偷精品视频| 免费播放大片免费观看视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 久久99蜜桃精品久久| 亚洲欧美日韩另类电影网站| 你懂的网址亚洲精品在线观看| 国产精品久久久久久av不卡| 精品人妻熟女毛片av久久网站| 成人影院久久| 日韩欧美一区视频在线观看| xxx大片免费视频| 丁香六月天网| 国产成人91sexporn| 高清在线视频一区二区三区| 桃花免费在线播放| 欧美97在线视频| 男女边吃奶边做爰视频| 国产不卡av网站在线观看| 成人漫画全彩无遮挡| 日韩,欧美,国产一区二区三区| 人人妻人人添人人爽欧美一区卜| 高清视频免费观看一区二区| 日韩成人av中文字幕在线观看| 精品少妇黑人巨大在线播放| 日韩成人av中文字幕在线观看| 欧美人与善性xxx| 久久综合国产亚洲精品| 亚洲成av片中文字幕在线观看 | 极品人妻少妇av视频| 久久精品国产a三级三级三级| 久久精品久久久久久噜噜老黄| 欧美 日韩 精品 国产| 视频区图区小说| 亚洲国产精品国产精品| 国产成人91sexporn| 精品99又大又爽又粗少妇毛片| 视频区图区小说| 97在线人人人人妻| 一级毛片 在线播放| 在线免费观看不下载黄p国产| 国产精品嫩草影院av在线观看| 日韩一区二区视频免费看| 一本色道久久久久久精品综合| 女的被弄到高潮叫床怎么办| 丝袜美腿诱惑在线| 久久午夜综合久久蜜桃| 少妇的逼水好多| 男女高潮啪啪啪动态图| 国产精品久久久av美女十八| 最黄视频免费看| 伊人久久大香线蕉亚洲五| 极品人妻少妇av视频| videos熟女内射| 99热国产这里只有精品6| 九九爱精品视频在线观看| 七月丁香在线播放| 欧美精品高潮呻吟av久久| 2018国产大陆天天弄谢| 成年女人在线观看亚洲视频| 狠狠婷婷综合久久久久久88av| 2022亚洲国产成人精品| 在线观看www视频免费| 日韩一本色道免费dvd| 亚洲精品成人av观看孕妇| 久久99蜜桃精品久久| 晚上一个人看的免费电影| 精品国产乱码久久久久久小说|