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

    Parameters Optimization of the Dynamic Absorber to Control the Axial Vibration of Marine Shafting System

    2012-09-22 07:15:50LILiangweiZHAOYaoLITianyunLUPo
    船舶力學 2012年3期

    LI Liang-wei,ZHAO Yao,LI Tian-yun,LU Po

    (School of Naval Architecture&Ocean Engineering,Huazhong University of Science and Technology,Wuhan 430074,China)

    1 Introduction

    Periodical force and torque due to the propeller in non-uniform flow field will cause shaft torsional vibration,axial vibration,lateral vibration and the coupling vibration.Oscillation caused by a small variation of thrust when the blades rotate through the non-uniform wake results in axial resonance at blade frequency[1].Many researchers have explored the axial vibrations of the marine shafting system,which are directly related to the quietness and comfort for passengers.A series of methods has been used for axial shaft vibration reduction,for example,by changing the stiffness of the marine thrust bearing,increasing damping in shafting foundation with the materials that have high damping properties,or using vibration control.One of the effective and feasible ways for vibration control is to install dynamic absorbers to marine shafting system.

    The dynamic absorber attached to the primary system as a subsystem reduces vibration by introducing opposite force and also by dissipating energy through damping.The key to control the main system vibration is to clarify the parameters of dynamic absorber and locations.With the development of the different structural forms and the increasing needs of vibration control,the theory of parameters optimization has been improved.

    The theory of dynamic vibration employed to one degree of freedom system without damping was first systematically investigated by Ding[2].Then,Randall et al[3]determined the optimal linear dynamic absorber for linear damping system based on minimum-maximum amplitudes criteria to find the optimal parameters.Because of popular application of the dynamic absorber in the field of machine engine,marine,automobiles,buildings,and so on,several researchers started to pay more attention to the methods of oversimplification to represent the original system.Ram and Elhay[4]studied the behavior of multi-degree-of-freedom vibration system with a simple dynamic absorber,and they proposed the recoverage of the physical parameters of a subsystem from the frequency response function measurements.Later on,a closed-form vibration transmissibility of a general unidirectional multi-degree-of-freedom system with multiple dynamic absorbers was discussed by Hsueh[5].Based on the graph model,the force and displacement transmissibility were derived by using the topology scheme.Two examples were shown for the implementation of this method.

    Optimization and numerical algorithms were also presented in several literatures.Kitis[6]published a paper on the design of two dynamic absorbers simultaneously applied to an undamped cantilever beam with the re-analysis technique,and the design parameters were optimized by the method of feasible direction.However,in such parameters optimization of the dynamic absorber,only the optimal damping and stiffness were obtained except mass.A new geometrical design method used to simplify the stiffness matrix of dynamic absorbers for the multi-degree-of-freedom system was presented by Jang and Choi[7].Moreover,the performances of different frequency-tuning methods were examined by Sun et al[8].It was shown that the oneone method almost had the same characteristics as the optimization method,and it did not need time-consuming optimization process.

    Dylejko et al[9]analyzed the propeller-shafting system with the four-pole parameter method and used the genetic algorithm to optimize the virtual resonance changer parameters.However,the transmission matrix parameters of the RC are not derived in detail and also the influence on different locations of the RC to decrease the axial vibration was not concerned.Viana et al[10]reported an analytical solution for the parameters optimization of two different dynamic absorbers by ant colony optimization and defined an objective function used to minimize the vibration amplitude of the primary structure.Furthermore,a comprehensive review of dynamic absorber including configurations developments and applications could be found in the papers by Kela et al[11].

    After reviewing the literatures listed above,one should notice that most articles dealt with the theoretical models to optimize the parameters of dynamic absorber,while the practical en-gineering problems are less discussed.For example,little attention has been devoted to use dynamic absorber for specific resonance attenuation and solve the parameter optimization problems with the finite element method.Such problems do have considerable practical significance and therefore will be discussed in the following.

    In this paper,the marine shafting system is modeled as multiple degree-of-freedom system based on the finite element method.The dynamic absorber is used to control the axial vibration of marine shafting system.It is essential to mention that the marine shafting system simply includes propeller,shafting,coupling,journal bearing and thrust bearing.The dynamic response of axial vibration is obtained by the re-analysis technique,which reduces the calculation time.The force transmissibility and the power transmissibility at the thrust bearing are derived to examine the performance of dynamic vibration.With the methods of the genetic algorithm and the multi-objective optimization algorithms,the dynamic absorber parameters are optimized over a certain response frequency range and the problem of specific resonance attenuation is also solved by adding equality constraint in the optimization process.Dynamic vibration is applied to three different locations at marine shafting system to investigate the influence on the axial vibration with different objective functions.The results here validate the feasibility of the optimization algorithm presented,and show that the vibration transmissibility is able to serve as a good criterion for the marine shafting system design as well as vibration control.

    2 Finite element model

    2.1 The axial vibration model of marine shafting system

    The marine shafting system mainly consists of propeller,shafting,coupling,journal bearing,thrust bearing and their foundation.The components of simplified mathematical model are the lumped mass and axial stiffness,which make an impact on the behaviors of the shafting axial vibration.The axial vibration model of marine shafting system is shown in Fig.1.

    Fig.1 The sketch of axial vibration of the marine shafting system

    In order to get the motion equation of the axial vibration based on the finite element method,the marine shafting system is discretized into several elements.The position of the thrust bearing and the coupling should generate nodes according to the simplified criterion of the finite element model.The transmission shaft is generally simulated as beam element in numerical analysis.The mass of the propeller,the couplings and the thrust bearing are modeled as lump masses averagely assigned to their left and right element mass matrices.Similarly,the stiffness of thrust bearing is equivalent as axial stiffness averagely assigned to its left and right stiffness matrix.

    The analytical expression of the beam element axial displacement is solved by using the appropriate displacement interpolation function combined with the boundary conditions.By assembling the matrix of the element,the motion expression of shaft axial vibration is obtained[12].

    2.2 Model of dynamic vibration

    Suppose a single dynamic absorber attached to the marine shaft system along the axial direction at node i and the schematic is given in Fig.2,here miis the mass of marine shafting in the finite element model at node i.The dynamic absorber is composed of a virtual mass m1,a virtual stiffness k1and a virtual damping c1.The displacement of the masses at node i and the displacement of dynamic absorber are shown as ui(t)and u1(t),respectively.According to the D’Alembert principle,the dynamic absorber equation of motion is established as follows:

    Fig.2 Dynamic absorber system

    Substituting the Eq.(1)into the matrix forms:

    After installing the dynamic absorber to the marine shafting system,Fig.3 shows the finite element model of axial vibration of marine shafting system,which has n nodes.Kthis the axial stiffness of thrust bearing and the hull is equivalent as rigid boundary.

    Fig.3 The finite element model of axial vibration of marine shafting system

    The motion equations of axial vibration of marine shafting system could be written as follows based on the performance of typical vibration system.

    where[M],[K ]and[C ]are the mass,stiffness and damping matrices of this system.{F(t)}is the axial excitation force vectors.{u},{u˙}and{u¨}are the axial displacement,velocity and acceleration vectors,respectively.

    3 Parameters optimization of the dynamic absorber

    3.1 Solution of the vibration transmissibility

    As an objective function to evaluate the performance of dynamic absorber,the force transmissibility and the power transmissibility are widely used.The marine structural vibration through the thrust bearing induced by shafting system is the main reason to analyze the axial vibration.This paper will discuss the effect of a dynamic absorber to reduce the axial vibration,using the force transmissibility and the power transmissibility at the thrust bearing as a criterion.

    Suppose that the marine shafting system is excited by a periodically varying force with amplitude vector {F0},{F(t)}in Eq.(3)will be the form:

    Correspondingly,the axial displacement vector of the shafting system is also written as:

    where {u0}is the amplitude of the axial displacement vector.Substitution of Eq.(4)and Eq.(5)into Eq.(3),then takes the form:

    When the dynamic absorber parameters change in the calculation procedures,the mass matrix,stiffness matrix and damping matrix in Eq.(6)will also regenerate.Define the parameter vector x as follows:

    where q is the number of dynamic absorbers attached to the shafting system.

    This paper will use the re-analysis technology to solve the equation of motion in Eq.(3),which could reduce the structure analysis time[6].When the dynamic absorber parameters change,the new mass matrice[M′],stiffness matrice[K′]and damping matrice[C′]of the vibration system in Eq.(6)could be expressed as follows:

    where the Δ [M(x)],Δ[K(x)]and Δ[C(x)]are the variable quantity of each matrix and they are the function of the vector x.

    In order to analyze the dynamic response described as Eq.(6),the initial values of dynamic absorber are required to start the numerical computation.Then,the matrices of the system are updated according to the Eqs.(8-10).Based on the solution technology dealing with the vibration problems of multiple degree-of-freedom system[13],e.g.,the substructure mode synthesis or model reduction,the displacement,velocity and acceleration from Eq.(6)can be finally obtained.

    In the situation of the axial vibration of the marine shafting system,one can assume that the original system is excited a unit harmonically load at the propeller.With the Eqs.(6-10),the axial displacement as well as the axial velocity at each node could be solved.The force transmissibility and the power transmissibility at the thrust bearing are expressed as:

    As shown in the above equations,T(x,ω)and W(x,ω)are the functions of three main variables at thrust bearing,the force FT,the axial displacement un-2and axial velocity vn-2.The interval ω?[ωmin,ωmax]is the response frequency range of the marine shafting axial vibration.

    3.2 Objective function

    Aiming to attenuate the axial vibration of the marine shafting system over the concerned response frequency range,appropriate algorithms should be explored to solve the design variables x of the dynamic absorber.It should be noted that the marine structure vibration through the thrust bearing is the main vibration transmission path.As a result,the response at thrust bearing is to be minimized over a frequency range ωmin≤ω≤ωmaxby attaching the dynamic absorber.Correspondingly,the objective function based on the force transmissibility and the power transmissibility is given by:

    The physical meaning of Eq.(12)could be expressed as:the issue of dynamic absorber parameters optimization is to minimize the maximum value of the force transmissibility and the power transmissibility at the thrust bearing over the frequency range,and finally achieve the purposes of controlling the axial vibration of the marine shafting system.

    3.3 Constraints

    When the dynamic absorber is used to control the marine shafting system,the variables of dynamic absorber should satisfy the design requirements for parameters optimization.Define the xminand xmaxas the lower and upper limits of the design variables x,respectively,and the range of variables can be written as inequality constraints:

    For the non-linear optimization problem in this paper,by combining the objective function in Eq.(12)with the constraints in Eq.(13)and choosing the external point penalty function method,one can express the regenerate objective function as objk′(x,ω,)R :

    where R is the weight factor between the objective functions and the constraints.I is the number of inequality constraints.

    3.4 Optimization algorithm description

    The objective function described in Eq.(14) could be solved with different optimization algorithm techniques.Genetic algorithm is one of the most suitable methodologies,which is a search technique used to find the global approximate solutions to optimization.The significant advantage of the techniques for design values optimization is that there are no strict analytic objective requirements,and it is also suitable for a variety of variables,especially the discrete ones.The details of procedures for solving the Eq.(14)coupled with the genetic algorithm and multi-objective optimization algorithm are described in the following.

    The response frequency band is first discretized.The vibration transmissibility described in Eq.(11)at the ith discrete frequency is written as Yk(x,ωi),where the ωiis the discrete frequency.Using the objective functions in Eq.(14)as the fitness functions,the maximum valuecan be obtained if the parameters of the genetic algorithm are defined,such as individual numbers,generation numbers,and variable precision.All the maximum values at discrete frequencies will constitute a vector,where the s is the number of discrete frequency point.Then,the maximum value of the vibration transmissibility is minimized with the linear weighting method to solve the multi-objective optimization problem.Finally,by transforming Yk(x,ωi)to the dimensionless form at every discrete frequency,the evaluation function is obtained to solve the objective function Eq.(14)based on the linear weighting sum[14].

    where αiis the weighting coefficient which assigns the relative importance of the criterion and.Optimization problem described as Eq.(11)will change the analysis of the objective function(x,ω,R).Eq.(15)will be used to optimize the dynamic absorber parameters over the frequency range with the genetic algorithm and the multi-objective optimization algorithm.

    4 Numerical example

    4.1 Parameters optimization over 0-300 Hz

    To illustrate the numerical method described above,an example of a marine shafting system is presented and a single dynamic absorber is used to control the axial vibration of the original system.The finite element model is established based on the sketch of the marine shaft system shown in Fig.4.

    Fig.4 The finite element model of the marine shafting system

    This sketch contains 19 nodes,17 beam elements,six lump masses and one axial spring.The three natural frequencies of the axial vibration in this example over 0-300 Hz frequency range are 51 Hz,122 Hz and 256 Hz.

    The force transmissibility and the power transmissibility versus frequency of the original marine shafting system are shown in Fig.5.The following force transmissibility and the power transmissibility are calculated with base-10 logarithm.

    Fig.5 The vibration transmissibility curve of the original marine shafting system

    The paper is to study the characteristics of dynamic absorber at different positions in the marine shafting system.Since the transmission shaft mainly consists of the stern shaft,the intermediate shaft and the thrust shaft,the dynamic absorbers are applied to the three locations of the shaft segments shown in Fig.4,which are the location A,location B and location C.

    The optimization parameters of dynamic absorber are solved with objective functions(x,ω,R)given in Eq.(15)based on the optimization algorithm shown above.In addition,the effect of different locations to control the axial vibration is explored.The lower limit and the upper limit of the dynamic absorber within the optimization process are shown in Tab.1.

    Tab.1 The limits values of the dynamic absorber

    Using the force transmissibility at the thrust bearing as the criterion to solve the objective function in Eq.(15),the optimal values of the dynamic absorber are given in Tab.2,as well as the natural frequency f1.Fig.6(1)shows the force transmissibility versus frequency for the cases with dynamic absorber at three different locations.The resonance frequencies of shafting system and the amplitude of resonance peaks are different for the dynamic absorber at different positions.

    Tab.2 Optimal values of dynamic absorber with the force transmissibility

    Fig.6 The vibration transmissibility of different dynamic absorbers locations(over 0-300 Hz)

    A comparison of the resonance frequencies of shafting system is shown in Fig.6(1).The first natural frequency becomes smaller than that of the shafting system without adding dynamic absorber.It is also clear from the figure that all the other resonance peaks,except the first peak,become lower than the corresponding resonance peaks without the dynamic absorber.As a result,after applying a single dynamic absorber to the original system,the second and third resonance peaks are suppressed when the optimal parameters of dynamic absorber are set.However,the first resonance peak attenuation is not obvious.

    Simultaneously,the influence on the dynamic absorber located at three different positions in the first resonance peak attenuation is almost the same.The amplitude of second resonance peak is decreased obviously when the dynamic absorber is fixed on the intermediate shaft,and the amplitude of third resonance peak reduces more when the dynamic absorber is located at the stern shaft.

    With the power transmissibility at the thrust bearing as the criterion to solve the objective function in Eq.(15),the optimal values of the dynamic absorber parameters are also given in Tab.3.A plot of the power transmissibility versus frequency curve with the dynamic absorber at three different locations is shown in Fig.6(2).

    Tab.3 Optimal values of dynamic absorber with the power transmissibility

    Comparing Fig.6(1)and Fig.6(2),the phenomenon of dynamic absorber on axial vibration absorption based on the two objective functions is almost the same.For example,the amplitude of the second and third resonance peaks become smaller when the dynamic absorber is applied to the three different locations.However,it is observed from the Fig.6(2)that there is an anti-resonance peak at the resonance frequency 122 Hz when a single dynamic absorber is attached at the stern shaft and the intermediate shaft.This value is equal to the second resonance frequency of the original system.It is demonstrated that the resonance peaks could be suppressed when the natural frequency of dynamic absorber is equal to that of the primary system.Therefore,parameter optimization of dynamic absorber is also the process of finding the appropriate parameters to solve the system oscillation at certain frequency.

    4.2 Specific resonance frequency attenuation

    The numerical example discussed above cares about the parameters optimization of the dynamic absorber over the entire frequency range.Providing that the problem of reducing the first resonance frequency has been increasingly a major concern,the method to suppress the first resonance frequency to control the axial vibration of original system will be presented as follows.

    Based on the performance of the dynamic absorber to suppress the original system oscillation,an equality constraint will be added to the objective functions in Eq.(15).The method could enable the natural frequency of the dynamic absorber equal to the first resonant frequency.The equality constraint is expressed as follows:

    where f1is the first resonance frequency of the original marine shafting system.

    Substituting Eq.(16)into Eq.(15)to optimize the parameters of dynamic absorber,the new objective functions(x,ω,R,S)is given in Eq.(17)

    where S is the weight factor of the equality constraints.The optimal values are obtained with the objective functions in Eq.(17),as shown in Tab.4 and Tab.5.

    Tab.4 Optimal values for specific resonance attenuation with the force transmissibility

    As can be seen from Tab.4 and Tab.5,for the dynamic absorbers at the three different locations,the natural frequencies of the dynamic absorber based on the two objective functions are close to the first resonance frequency of the original marine shafting system.It is demonstrated that the specific resonance could be attenuated by adding the equality constraint to the objective functions.However,the natural frequencies of the dynamic absorber are not exactly the same regardless of the different objective functions or different locations are chosen.It is the reason why the penalty factor and the number of genetic generation are used in the genetic algorithm process.With the force transmissibility and the power transmissibility to analyze the effect of specific resonance attenuation,the different curves are shown in Fig.7(1)and Fig.7(2).

    Tab.5 Optimal values for specific resonance attenuation with the power transmissibility

    Fig.7 The vibration transmissibility for specific resonance attenuation

    The amplitude of the first resonance peak shown in the two figures is significantly suppressed.With the method of adding the equality constraint to optimize the parameters,the purpose of reducing the specific peaks could be achieveed.A comparison of the three different locations for the first resonance attenuation is presented.It is shown that the efficiency of attaching dynamic absorber to the stern shaft is better than that to the intermediate shaft.The amplitude reduction of the second and third resonance peaks is almost the same based on the two objective functions in Eq.(17),while the dynamic absorber installed to the stern shaft could be the best third resonance peak suppression.

    5 Conclusions and discussion

    Based on the previous results,we could draw the following comments and conclusions:

    1)With the method of the genetic algorithm and multi-objective optimization algorithm,the problem of parameters optimization to control the axial vibration of the marine shafting system could be solved.Over the 0-300 Hz frequency range,the phenomenon of three resonance frequency reduction is similar with the two objective functions.

    2)Comparing to the optimal results with the two objective functions over the 0-300 Hz frequency range,it is obviously proposed in this paper that the second and three resonance frequencies of the marine shafting system could be better attenuated with the optimization algorithm.However,it is shown that less effect is made to control the first resonance peak.

    3)In order to suppress the first resonance peak,the equality constraint is introduced based on the characteristic of the dynamic absorber.Optimization results show that the method presented here could control the first resonance response.Comparison of the attenuation effect of different locations shows that,the three resonance peaks could be controlled better when dynamic absorber is added to the stern shaft.

    [1]Xu Yunxiu,Zhong Xuetian,He Xuanxuan.Axial vibration of the ship shafting system[M].Beijing:China Communications Press,1985.

    [2]Ding Wenjing.Theory of vibration absorption[M].Beijing:Tsinghua University Press,1988.

    [3]Randall S E,Halsted D M,Taulor D L.Optimum vibration absorbers for linear damped systems[J].Journal of Mechanical Design,1981,103:908-913.

    [4]Ram Y M,Elhay S.The theory of a multi-degree-of-freedom dynamic absorber[J].Journal of Sound and Vibration,1996,195(4):607-615.

    [5]Hsueh W J.Vibration transmissibility of a unidirectional multi-degree-of-freedom system with multiple dynamic absorbers[J].Journal of Sound and Vibration,2000,229(4):793-805.

    [6]Kitis L.Vibration reduction over a frequency range[J].Journal of Sound and Vibration,1983,89(4):559-569.

    [7]Jang S J,Choi Y J.Geometrical design method of multi-degree-of-freedom dynamic vibration absorbers[J].Journal of Sound and Vibration,2007,303(2):343-356.

    [8]Sun H L,Zhang P Q,Chen H B,Zhang K,Gong X L.Application of dynamic vibration absorbers in structural vibration control under multi-frequency harmonic excitations[J].Applied Acoustics,2008,69:1361-1367.

    [9]Dylejko P G,Kessissoglou N J,Tso Y,Norwood C J.Optimization of a resonance changer to minimise the vibration transmission in marine vessels[J].Journal of Sound and Vibration,2007,300:101-116.

    [10]Viana F A C,Kotinda G I,Rade D A,Steffen V.Tuning dynamic vibration absorbers by using ant colony optimization[J].Computers and Structures,2008,86(14):1539-1549.

    [11]Kela L,V?h?oja P.Recent studies of adaptive tuned vibration absorbers/neutralizers[J].Applied Mechanics Reviews,2009,62(6):1-9.

    [12]Zhang Zhihua.Numerical calculation of power equipment vibration[M].Harbin:Harbin Engineering University Press,1994.

    [13]Clough R W,Penzien J.Dynamics of structures[M].New York:McGraw-Hill,1993.

    [14]Zeng Guangwu.Optimization design of ship structure[M].Wuhan:Huazhong University of Science and Technology Press,2004.

    亚洲精品乱码久久久v下载方式| 黄色欧美视频在线观看| 人人妻人人澡人人爽人人夜夜 | 日日啪夜夜撸| 夫妻午夜视频| 最近2019中文字幕mv第一页| 青春草视频在线免费观看| 少妇高潮的动态图| 免费高清在线观看视频在线观看| 97在线视频观看| 国语对白做爰xxxⅹ性视频网站| 能在线免费看毛片的网站| 久久久久精品性色| 最近最新中文字幕大全电影3| 久久精品国产鲁丝片午夜精品| 丝瓜视频免费看黄片| 欧美日韩亚洲高清精品| 亚洲精品456在线播放app| 成人一区二区视频在线观看| 高清午夜精品一区二区三区| 91午夜精品亚洲一区二区三区| 亚洲一级一片aⅴ在线观看| 亚洲国产精品成人综合色| 国产视频内射| 成人美女网站在线观看视频| 国产午夜精品一二区理论片| 在线免费观看不下载黄p国产| 久久精品国产自在天天线| 国产精品一区二区性色av| 国产伦一二天堂av在线观看| 天堂av国产一区二区熟女人妻| 伊人久久国产一区二区| 亚洲国产精品国产精品| 国产免费又黄又爽又色| 一级毛片aaaaaa免费看小| 亚洲三级黄色毛片| 日日摸夜夜添夜夜添av毛片| 午夜免费观看性视频| av线在线观看网站| 狠狠精品人妻久久久久久综合| 男人狂女人下面高潮的视频| 中文字幕亚洲精品专区| 熟妇人妻不卡中文字幕| 国产综合懂色| 最近最新中文字幕免费大全7| 国产成年人精品一区二区| 人人妻人人澡欧美一区二区| 人人妻人人澡人人爽人人夜夜 | ponron亚洲| 国产 一区精品| 亚洲成人av在线免费| 亚洲av在线观看美女高潮| 欧美日韩亚洲高清精品| 久久久久精品性色| 成人av在线播放网站| 国精品久久久久久国模美| 久久99热这里只有精品18| 国产午夜福利久久久久久| av.在线天堂| 熟女电影av网| 免费观看在线日韩| 精品国产露脸久久av麻豆 | 婷婷色麻豆天堂久久| 免费高清在线观看视频在线观看| or卡值多少钱| 91久久精品国产一区二区成人| 亚洲精品乱久久久久久| 在线观看美女被高潮喷水网站| 偷拍熟女少妇极品色| 女人十人毛片免费观看3o分钟| 春色校园在线视频观看| 免费黄色在线免费观看| 国产精品久久视频播放| 简卡轻食公司| 老女人水多毛片| 亚洲不卡免费看| 国产 一区 欧美 日韩| 大又大粗又爽又黄少妇毛片口| av女优亚洲男人天堂| 日韩精品青青久久久久久| 国产精品久久久久久久电影| 你懂的网址亚洲精品在线观看| 国产精品福利在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 国产黄a三级三级三级人| 美女xxoo啪啪120秒动态图| 色综合色国产| 熟女电影av网| 97人妻精品一区二区三区麻豆| 亚洲综合色惰| 日本猛色少妇xxxxx猛交久久| 国产成年人精品一区二区| 全区人妻精品视频| 精品国产露脸久久av麻豆 | 波多野结衣巨乳人妻| 又爽又黄a免费视频| 亚洲欧美清纯卡通| 九九爱精品视频在线观看| 国产探花极品一区二区| 日韩在线高清观看一区二区三区| 欧美不卡视频在线免费观看| 亚洲av日韩在线播放| av播播在线观看一区| 最新中文字幕久久久久| 中文精品一卡2卡3卡4更新| 久久久亚洲精品成人影院| 2021天堂中文幕一二区在线观| 神马国产精品三级电影在线观看| 亚洲国产高清在线一区二区三| 免费无遮挡裸体视频| 亚洲精品自拍成人| 日韩伦理黄色片| 搡女人真爽免费视频火全软件| 欧美日韩在线观看h| 国产精品福利在线免费观看| 亚洲精品久久久久久婷婷小说| 哪个播放器可以免费观看大片| 日韩大片免费观看网站| 日韩成人av中文字幕在线观看| 亚洲四区av| 少妇裸体淫交视频免费看高清| 午夜免费观看性视频| 亚洲欧美中文字幕日韩二区| 大话2 男鬼变身卡| 国产精品一区二区性色av| 搞女人的毛片| 国产成人精品福利久久| 国产有黄有色有爽视频| 日韩av不卡免费在线播放| 男人舔女人下体高潮全视频| 男人舔女人下体高潮全视频| 赤兔流量卡办理| 蜜桃久久精品国产亚洲av| 99热这里只有是精品50| 国产在线一区二区三区精| av在线老鸭窝| 尾随美女入室| 精品久久久久久成人av| 伊人久久国产一区二区| 日日摸夜夜添夜夜添av毛片| 精品国产一区二区三区久久久樱花 | or卡值多少钱| 又粗又硬又长又爽又黄的视频| 人妻制服诱惑在线中文字幕| 夜夜看夜夜爽夜夜摸| 街头女战士在线观看网站| 永久网站在线| 黄片无遮挡物在线观看| 搡女人真爽免费视频火全软件| 中文精品一卡2卡3卡4更新| 亚洲婷婷狠狠爱综合网| 99热这里只有是精品50| 少妇人妻一区二区三区视频| 美女cb高潮喷水在线观看| 啦啦啦韩国在线观看视频| 日韩av不卡免费在线播放| 日本免费a在线| 亚洲图色成人| 在线观看av片永久免费下载| 色播亚洲综合网| 国产91av在线免费观看| 男女边摸边吃奶| 欧美日韩视频高清一区二区三区二| 内地一区二区视频在线| 纵有疾风起免费观看全集完整版 | 卡戴珊不雅视频在线播放| 国产中年淑女户外野战色| 中国美白少妇内射xxxbb| 热99在线观看视频| 最近手机中文字幕大全| 一级av片app| 久久这里有精品视频免费| a级一级毛片免费在线观看| 中文资源天堂在线| 久久鲁丝午夜福利片| 亚洲国产高清在线一区二区三| 亚洲欧美精品自产自拍| 亚洲精品日韩av片在线观看| 日本猛色少妇xxxxx猛交久久| 高清av免费在线| 一级二级三级毛片免费看| 国产免费视频播放在线视频 | 蜜桃久久精品国产亚洲av| 99久久精品一区二区三区| 国产白丝娇喘喷水9色精品| 日韩av不卡免费在线播放| 久久久久久久久大av| 久99久视频精品免费| 国产探花极品一区二区| 国产精品三级大全| 亚洲av在线观看美女高潮| 国产成人精品久久久久久| 亚洲av电影在线观看一区二区三区 | 婷婷色av中文字幕| 黄色日韩在线| 久久久久久久午夜电影| 亚洲av不卡在线观看| 国产精品久久久久久久电影| 成人午夜高清在线视频| 成人午夜精彩视频在线观看| 亚洲av男天堂| 色视频www国产| 色网站视频免费| 色综合亚洲欧美另类图片| 精品人妻熟女av久视频| 天天一区二区日本电影三级| 男女边吃奶边做爰视频| 久久久久久久大尺度免费视频| 天天一区二区日本电影三级| 精品午夜福利在线看| 国产乱来视频区| h日本视频在线播放| 有码 亚洲区| 久热久热在线精品观看| 国产白丝娇喘喷水9色精品| 国产男人的电影天堂91| 五月天丁香电影| 国产一区二区三区综合在线观看 | 国产精品人妻久久久久久| 国产精品久久久久久久久免| 久久这里有精品视频免费| 亚洲成人av在线免费| 人人妻人人澡欧美一区二区| 国产亚洲午夜精品一区二区久久 | 久久久久性生活片| 国产午夜精品一二区理论片| 亚洲av中文av极速乱| 亚洲国产精品sss在线观看| 最近最新中文字幕免费大全7| 能在线免费看毛片的网站| 国产 一区 欧美 日韩| 免费观看av网站的网址| 国产精品国产三级国产av玫瑰| 最近手机中文字幕大全| 日本一本二区三区精品| 国产免费福利视频在线观看| 亚洲av成人av| 国产片特级美女逼逼视频| 免费人成在线观看视频色| av网站免费在线观看视频 | 久久精品久久久久久久性| 91狼人影院| 97人妻精品一区二区三区麻豆| 国产高清有码在线观看视频| 亚洲av不卡在线观看| h日本视频在线播放| 日日撸夜夜添| 禁无遮挡网站| 天堂网av新在线| 男的添女的下面高潮视频| 成人特级av手机在线观看| 看免费成人av毛片| 夜夜看夜夜爽夜夜摸| 99久久精品热视频| 91在线精品国自产拍蜜月| 美女高潮的动态| 熟妇人妻久久中文字幕3abv| 黄片wwwwww| 久久韩国三级中文字幕| 亚洲av在线观看美女高潮| 国产在线一区二区三区精| 小蜜桃在线观看免费完整版高清| 99视频精品全部免费 在线| 国内揄拍国产精品人妻在线| av国产免费在线观看| 嫩草影院精品99| 欧美丝袜亚洲另类| 人体艺术视频欧美日本| 嫩草影院新地址| a级一级毛片免费在线观看| 美女国产视频在线观看| 婷婷色av中文字幕| 国产极品天堂在线| 国产精品一区www在线观看| 大香蕉久久网| 99九九线精品视频在线观看视频| av黄色大香蕉| 男人舔女人下体高潮全视频| 嘟嘟电影网在线观看| 亚洲精品aⅴ在线观看| 一个人观看的视频www高清免费观看| 国产精品人妻久久久久久| 亚洲精品色激情综合| 精品久久久久久成人av| 欧美日韩精品成人综合77777| 国产成人aa在线观看| 精品国产露脸久久av麻豆 | 乱码一卡2卡4卡精品| av免费观看日本| 午夜精品在线福利| 久久久久久久久久久丰满| 久久久久久久久久久丰满| 久久久a久久爽久久v久久| 国产黄色免费在线视频| 日韩av在线大香蕉| 少妇人妻精品综合一区二区| 久久久久久久久中文| 亚洲成色77777| 精品熟女少妇av免费看| 国产一区二区在线观看日韩| 亚洲综合色惰| 国产午夜福利久久久久久| 国产亚洲91精品色在线| 国产片特级美女逼逼视频| 免费观看精品视频网站| 国产精品一区二区三区四区免费观看| 精华霜和精华液先用哪个| 国产国拍精品亚洲av在线观看| 亚洲欧美日韩无卡精品| 91久久精品国产一区二区成人| 2021少妇久久久久久久久久久| 久久久a久久爽久久v久久| 亚洲欧美精品专区久久| 好男人视频免费观看在线| 熟妇人妻不卡中文字幕| 亚洲av日韩在线播放| 欧美精品国产亚洲| 非洲黑人性xxxx精品又粗又长| 青春草视频在线免费观看| 国产熟女欧美一区二区| 亚洲国产精品sss在线观看| 九九久久精品国产亚洲av麻豆| www.色视频.com| 最近2019中文字幕mv第一页| 国产在线男女| 国产高潮美女av| 亚洲av电影在线观看一区二区三区 | 蜜桃亚洲精品一区二区三区| 秋霞伦理黄片| 欧美精品一区二区大全| 亚洲欧洲日产国产| 欧美成人精品欧美一级黄| 国产黄片美女视频| 亚洲综合色惰| 国产精品美女特级片免费视频播放器| 亚洲av中文字字幕乱码综合| 2021天堂中文幕一二区在线观| 国产成人精品婷婷| 久久久精品欧美日韩精品| 青春草国产在线视频| 全区人妻精品视频| 69人妻影院| 日韩视频在线欧美| 国产极品天堂在线| 嫩草影院新地址| 亚洲三级黄色毛片| 亚洲欧美日韩无卡精品| 嫩草影院入口| 国产淫片久久久久久久久| 国模一区二区三区四区视频| 亚洲精品乱码久久久v下载方式| 你懂的网址亚洲精品在线观看| 亚洲自偷自拍三级| 韩国高清视频一区二区三区| 99视频精品全部免费 在线| or卡值多少钱| 午夜久久久久精精品| 神马国产精品三级电影在线观看| 少妇猛男粗大的猛烈进出视频 | 一本一本综合久久| 日韩精品有码人妻一区| 有码 亚洲区| 国产精品蜜桃在线观看| 亚洲乱码一区二区免费版| 亚洲av福利一区| 人妻少妇偷人精品九色| 中国美白少妇内射xxxbb| 日本wwww免费看| 亚洲欧洲日产国产| ponron亚洲| 日本午夜av视频| 熟女人妻精品中文字幕| 欧美日韩综合久久久久久| 国产黄片视频在线免费观看| 亚洲在线自拍视频| 黄片wwwwww| 久久精品久久久久久久性| 麻豆久久精品国产亚洲av| 色播亚洲综合网| 伦精品一区二区三区| 久久久久久久久大av| 最近的中文字幕免费完整| 亚洲自拍偷在线| 亚洲va在线va天堂va国产| 国产午夜精品一二区理论片| 国内揄拍国产精品人妻在线| 久久99热6这里只有精品| 我的老师免费观看完整版| 国产黄色小视频在线观看| 成人美女网站在线观看视频| 成人综合一区亚洲| 网址你懂的国产日韩在线| 久久久精品欧美日韩精品| 又爽又黄无遮挡网站| 国产高清不卡午夜福利| 欧美性感艳星| 三级国产精品片| 国产成人精品福利久久| 一级片'在线观看视频| 中文字幕免费在线视频6| 亚洲av日韩在线播放| 成人午夜高清在线视频| 中文资源天堂在线| 亚洲内射少妇av| 日本一本二区三区精品| 亚洲伊人久久精品综合| 国产精品蜜桃在线观看| 91午夜精品亚洲一区二区三区| 免费观看性生交大片5| 午夜免费男女啪啪视频观看| 五月玫瑰六月丁香| 老师上课跳d突然被开到最大视频| 亚洲精品日韩av片在线观看| 99热这里只有是精品在线观看| 日韩精品青青久久久久久| 久久久久精品性色| 成年免费大片在线观看| 欧美区成人在线视频| 综合色丁香网| av专区在线播放| 日本与韩国留学比较| 亚洲欧美中文字幕日韩二区| .国产精品久久| 亚洲国产欧美人成| 精品久久久久久成人av| 纵有疾风起免费观看全集完整版 | 欧美性感艳星| 夜夜爽夜夜爽视频| 亚洲欧美成人综合另类久久久| 一本一本综合久久| 伊人久久精品亚洲午夜| 欧美另类一区| 国产成人a区在线观看| 久久99热这里只有精品18| 久久久久性生活片| 亚洲欧美一区二区三区黑人 | 成人无遮挡网站| 亚洲精品一二三| 一级a做视频免费观看| 国产免费又黄又爽又色| 亚洲国产av新网站| 欧美三级亚洲精品| 日本免费a在线| 最近最新中文字幕大全电影3| 26uuu在线亚洲综合色| 韩国高清视频一区二区三区| 国产精品综合久久久久久久免费| 亚洲国产精品成人综合色| 精品国内亚洲2022精品成人| 亚洲av男天堂| 日本三级黄在线观看| 欧美日韩精品成人综合77777| videossex国产| 久久久久久久久久久免费av| 欧美性猛交╳xxx乱大交人| 高清视频免费观看一区二区 | 别揉我奶头 嗯啊视频| 国产午夜精品久久久久久一区二区三区| 久久久久精品性色| 日本-黄色视频高清免费观看| 国产高清有码在线观看视频| 欧美+日韩+精品| a级毛片免费高清观看在线播放| 午夜免费观看性视频| 久久这里只有精品中国| 国产精品久久久久久久久免| 国产精品福利在线免费观看| 亚洲四区av| 精品国产一区二区三区久久久樱花 | 国产av国产精品国产| 国产 一区 欧美 日韩| 午夜福利在线在线| 丰满少妇做爰视频| 人人妻人人澡人人爽人人夜夜 | 国产乱来视频区| 日韩人妻高清精品专区| 成人av在线播放网站| 最近中文字幕2019免费版| 日韩成人av中文字幕在线观看| 国产午夜福利久久久久久| 男女视频在线观看网站免费| 亚洲经典国产精华液单| 69人妻影院| 日韩成人伦理影院| 男女边吃奶边做爰视频| 亚洲在线自拍视频| 久久久久精品久久久久真实原创| 日本黄色片子视频| 永久免费av网站大全| 久久国内精品自在自线图片| 国产白丝娇喘喷水9色精品| 视频中文字幕在线观看| 丝袜喷水一区| 日韩伦理黄色片| 蜜桃亚洲精品一区二区三区| freevideosex欧美| 最近的中文字幕免费完整| 精品一区二区三区人妻视频| 午夜福利视频1000在线观看| 内射极品少妇av片p| 日韩av免费高清视频| 欧美日韩亚洲高清精品| 成人毛片a级毛片在线播放| 久久精品久久精品一区二区三区| 观看免费一级毛片| 日韩中字成人| 久久精品国产亚洲网站| 欧美日韩一区二区视频在线观看视频在线 | 天堂俺去俺来也www色官网 | 日韩人妻高清精品专区| 成人无遮挡网站| 一个人看视频在线观看www免费| 国产亚洲午夜精品一区二区久久 | 国产精品熟女久久久久浪| 亚洲在久久综合| 一边亲一边摸免费视频| 欧美成人精品欧美一级黄| 亚洲精品自拍成人| 婷婷色综合www| 黄色配什么色好看| 亚洲欧洲国产日韩| 精品一区二区三区视频在线| 国产伦一二天堂av在线观看| av免费在线看不卡| 久久精品国产鲁丝片午夜精品| 日韩伦理黄色片| 晚上一个人看的免费电影| 国产白丝娇喘喷水9色精品| 亚洲乱码一区二区免费版| 日韩强制内射视频| 亚洲久久久久久中文字幕| 在线观看一区二区三区| 99热这里只有精品一区| 搡女人真爽免费视频火全软件| 国产久久久一区二区三区| 男女下面进入的视频免费午夜| 亚洲精品国产av蜜桃| 成人毛片60女人毛片免费| 男女下面进入的视频免费午夜| 建设人人有责人人尽责人人享有的 | 亚洲在线观看片| 亚洲欧美成人综合另类久久久| 日韩一区二区三区影片| 午夜福利网站1000一区二区三区| 两个人视频免费观看高清| 国产精品不卡视频一区二区| 少妇的逼好多水| 观看免费一级毛片| 久久久久久久久大av| 久热久热在线精品观看| freevideosex欧美| 男人爽女人下面视频在线观看| 免费大片黄手机在线观看| 欧美另类一区| 欧美一级a爱片免费观看看| 高清日韩中文字幕在线| 极品教师在线视频| 一区二区三区高清视频在线| 国模一区二区三区四区视频| 成人高潮视频无遮挡免费网站| 亚洲内射少妇av| av线在线观看网站| 亚洲欧美精品专区久久| 久久久色成人| 成人高潮视频无遮挡免费网站| 有码 亚洲区| 中文字幕久久专区| 久久久久久久久大av| 天堂俺去俺来也www色官网 | 亚洲精品乱码久久久久久按摩| 97热精品久久久久久| 尾随美女入室| 69人妻影院| 大陆偷拍与自拍| 一本久久精品| 边亲边吃奶的免费视频| 97精品久久久久久久久久精品| 欧美日韩在线观看h| 蜜桃亚洲精品一区二区三区| 最新中文字幕久久久久| 午夜久久久久精精品| 国产成人精品婷婷| 久久久欧美国产精品| 欧美日韩视频高清一区二区三区二| 麻豆成人av视频| 日本免费a在线| 欧美日韩综合久久久久久| 国产精品久久久久久精品电影小说 | 卡戴珊不雅视频在线播放| 夜夜看夜夜爽夜夜摸| 精品人妻熟女av久视频| 神马国产精品三级电影在线观看| 如何舔出高潮| 我的女老师完整版在线观看| 日韩欧美一区视频在线观看 | 精品欧美国产一区二区三| 国产精品综合久久久久久久免费| 成人鲁丝片一二三区免费| 日韩不卡一区二区三区视频在线| av在线老鸭窝| 一级毛片电影观看| 中国美白少妇内射xxxbb| 丝袜美腿在线中文| 熟女电影av网| 免费在线观看成人毛片| 91av网一区二区| 男女那种视频在线观看| 亚洲熟妇中文字幕五十中出| 成人综合一区亚洲| 最近视频中文字幕2019在线8| 在线播放无遮挡| 亚洲av一区综合| av福利片在线观看| 高清欧美精品videossex| 亚洲精品视频女|