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

    Cage slip characteristics of a cylindrical roller bearing with a trilobe-raceway

    2018-03-21 05:29:03SirDENGYujiaLUWnhuZHANGXuSUNZhnwiLU
    CHINESE JOURNAL OF AERONAUTICS 2018年2期

    Sir DENG,Yujia LU,Wnhu ZHANG,Xu SUN,Zhnwi LU

    aSchool of Mechanical Engineering,Henan University of Science and Technology,Luoyang 471003,China

    bCollaborative Innovation Center of Major Machine Manufacturing in Liaoning,Dalian 116024,China

    cPatent Examination Cooperation Center of the Patent Office,SIPO,Henan,Zhengzhou 450002,China

    dSchool of Mechanical Engineering,Northwestern Polytechnical University,Xi’an 710071,China

    eTechnology Center of LUOYANG LYC BEARING CO.,LTD,Luoyang 471003,China

    1.Introduction

    For high-speed cylindrical roller bearings in the application of aircraft engines,cage slip always happens due to high-speed and light-load conditions,and serious cage slip will lead to early failure of such bearings,which has a profound effect on the aircraft safety.One of the effective solutions to this light–load slip problem is to make the outer raceway of a cylindrical roller bearing to be a trilobe-raceway.1There are three man-made preload locations circumferentially equally-spaced on the outer ring raceway of a cylindrical bearing to increase preload onto more rollers,decrease the rotational speed difference at different positions,reduce bearing skidding,and improve its service life.

    The skidding problem of bearings at high speed with insufficient load has always been a hot concern by scholars both at home and abroad.Harris2developed a skid prediction model for a high-speed roller bearing by means of a quasi static analytical method,and investigated the effects of bearing load,rotating speed,and number of rolling elements on bearing skidding.Poplawski3established a quasi dynamics analysis model for a high-speed roller bearing taking into consideration the guide surface between the cage and rings as well as the friction between the rolling element and the cage pocket,and analyzed the cage and roller slip as well as t4he film thickness and forces in the cage.Cavallaro and Nelias presented an analytical model for high-speed cylindrical roller bearings with flexible rings,and investigated the relationship between ring deformation and bearing load distribution as well as roller slipping speed and load distribution.Takafumi et al.5proposed a theoretical model for cage slip of cylindrical roller bearings with the consideration of non-Newtonian fluid properties and temperature rise of lube oil,and analyzed the cage slip ratio under different speeds and loads.Arthanari and Marappan6experimentally analyzed the effects of rotating speed,radial load,and the number of rollers on cage slip of a cylindrical roller bearing.Takabi and Khonsari7built a dynamic model of a high-speed cylindrical roller bearing,and investigated the influence of different traction models on the sliding velocities and cage wear rate between rollers and races.Chinese scholars also did lots of research on bearing skidding.Li and Wu8established a dynamic analytical model for highspeed cylindrical roller bearings,and dynamically simulated the roller skew,axial movement,roller and cage slip characteristics,etc.Hu et al.9,10established a quasi static calculation model for inter-shaft cylindrical roller bearings,and analyzed the relationship between the roller slip ratio and radial load when bearing rings are under different working conditions.Tu and Shao11considered the acceleration phase of roller bearings,and presented an analytical model to study bearings’skidding characteristics under different forces and inner ring rotational accelerations.Chen et al.12developed a bearing dynamic model under a whirling condition taking into consideration the bearing skidding of a high-speed rolling bearing under the whirling condition,and analyzed the effects of various outer loadings,whirling frequencies,and whirling radii on the bearing skidding.Deng et al.13established a dynamic differential equation for a high-speed cylindrical roller bearing,and analyzed the effects of the clearance ratio of the cage,guiding type,etc.on the cage slip ratio and centroid trajectory.All the studies mentioned above have been focused on the performance analysis of cylindrical roller bearings under highspeed and light-load conditions,whereas studies on a cylindrical roller bearing with a trilobe-raceway have been focused on its raceway process technology,1,14,15and little research has been done on the theoretical research of the dynamics of a cylindrical roller bearing with a trilobe-raceway.

    Therefore,this paper analyzes the dynamics of a cylindrical roller bearing with a trilobe-raceway based on the dynamic analysis of the roller bearing,establishes nonlinear dynamics differential equations for a cylindrical roller bearing with a trilobe-raceway,and then adopts the GSTIFF(gear stiff)integer algorithm with a variable step to solve these equations.The focus of this paper is to study the influences of structural parameters of the raceway and parameters of working conditions as well as the relationship between the outer ring installation method and cage slip.The present paper provides some theoretical basis for the structure design of a cylindrical roller bearing with a trilobe-raceway.

    2.Calculation model of the trilobe-raceway

    An outer ring with a trilobe-raceway can be obtained by a predeformation machining method,1and the principle of the method is shown in Fig.1.

    As shown in Fig.1(a),a uniformly-distributed loadFis applied to the 1/3 symmetrical outer surface of the roughcast of the bearing ring along the circumferential direction,which leads to an elastic pre-deformation shown in Fig.1(b).Keep the applied load constant while grinding the raceway to its theoretical base circle dimension as shown in Fig.1(c),and the required trilobe-raceway contour can be obtained as a result of the deformation resilience after releasing the load,as shown in Fig.1(d).The radius of the trilobe-raceway at different azimuth angles after machining can be expressed as:

    whereReois the original radius before machining;Reg(φ)is the radius of the raceway at different azimuth angles φ after the elastic deformation of the outer ring is generated because of the applied force;Regmis the grinding radius,i.e.,the theoretical base circle radius.

    Fig.1 Sketch map of the pre-deformation machining principle.

    2.1.Finite element analytical model of the outer ring with the trilobe-raceway

    Take a certain type of bearing as an example to study the relationship between the trilobe-raceway’s structural parameters and the application of loads.The finite element analytical model of the outer ring with the trilobe-raceway is established by adopting 20-node hexahedron units.The model established in Ref.1 was simplified based on the geometrical relationship analyzed above,and the deformation caused by the grinding force in unloaded areas of the outer ring as well as the elastic deformation of the mold can be neglected.Based on the loads distribution needed for forming a trilobe-raceway,when under ideal processing conditions,gradient loads are applied to the 1/3 symmetrical outer surface of the roughcast of the bearing ring along the circumferential direction,and the included angles between loaded areas are 60°;one of the loaded areas is rotated at a slight degree φfto simulate non-ideal processing conditions,and by this way,the included angles between loaded areas become 60°,60°+ φf,and 60°- φf.The finite element analytical model of the outer ring with the triloberaceway is shown in Fig.2.

    Fig.3 shows the impact of the loads rotating degree φfon the trilobe-raceway’s structural parameters when the original radius before machining is 34.86 mm and the applied loads are 2.0 MPa.

    As shown in Fig.3,with an increase of the loads rotating degree,the high point of the trilobe-raceway counter increases and the low point decreases,i.e.,the roundness of the base circle at low and high points rises,and the D-value between the highest and lowest points is also increased.This is because when a loaded area is rotated,the magnitudes of loads also take suitable changes to achieve the balance of applied loads;in the area where loads are increased,the deformation of the outer ring under applied loads increases,and when the grinding radius remains unchanged,the radius at the low point of the raceway rises while the radius at the high point reduces after molding;likewise,in the area where loads are decreased,the radius of the low point reduces and the radius of the high point rises after molding.

    Fig.2 Finite element analytical model of the outer ring with the trilobe-raceway.

    Fig.3 Wave shape of the raceway imposed with asymmetrical loads.

    2.2.Mathematical model of the outer ring with the triloberaceway

    As shown in Fig.3,with the changes of the azimuth angle,the raceway radius is in accordance with a sinusoidal curve.In order to facilitate the expression of relationship between the trilobe-raceway radius and the azimuth angle,for an ideal trilobe-raceway,the obtained data can be fitted into a sinusoidal function in which the azimuth angle φ is regarded as the independent variable and the radius of the racewayReis regarded as the dependent variable as follows:

    Meanwhile,as for a non-ideal trilobe-raceway,the obtained data can be fitted into six sinusoidal functions in the same form as Eq.(2),where the azimuth angle φ is regarded as the independent variable,the radius of the racewayReis regarded as the dependent variable,andA,B,C,andDare real constants.

    3.Dynamic analytical model of a cylindrical roller bearing with a Trilobe-raceway

    The schematic diagram of a cylindrical roller bearing with a trilobe-raceway is shown in Fig.4.The outer ring is fixed,and the inner ring rotates at a constant speed ω and bears a radial forceFr.The cage is guided by the outer ring,and the bearing face is lubricated by oil jets.Assume that the surfaces of bearing components are absolutely smooth with an ideal geometric shape,and the center of mass is coincident with the center of form.

    In Fig.4,OiandO′iare the inner ring centers of the bearing before and after loading,respectively;P0dis the radial clearance of the theoretical base circle of the outer ring;φjis the azimuth angle of thejth roller;Pdjis the radial clearance between the theoretical base circle and the outer raceway at the azimuth angle φj,andPdj=Re(φj)-Regm;Recan be obtained from Eq.(2);δris the radial displacement of the inner ring.As can be seen from Fig.4,contact deformation between the raceway and part of rollers occurs before applying the radial force,where the radial clearance between the theoretical base circle and the outer raceway is negative and its absolute value is larger thanP0d.

    Fig.4 Schematic diagram of a cylindrical roller bearing with a trilobe-raceway.

    3.1.Force analysis of bearing parts

    3.1.1.Forces between a roller and raceways

    Tilt always happens when a cylindrical roller bearing with a trilobe-raceway operates under a radial forceFr,and the bending momentMand the forces and deformations between a roller and raceways are shown in Fig.5.

    In Fig.5,{O;X,Y,Z} is the inertial coordinate system of the bearing,which is a fixed coordinate system with theXaxis coincident with the rotation axis of the bearing,and theYZplane is parallel to the radial plane across the bearing center.Other coordinate systems are all established according to this inertial coordinate system; {orj;xrj,yrj,zrj} is the coordinate system of thejth roller center,whose original pointorjis coincident with the mass center of the roller,thexrjaxis is along with the roller axis,and theyrjzrjplane is parallel to the radial plane across the mass center of the roller.Every roller has its own local coordinate system which moves and rotates but doesn’t spin with the roller; α is the tilting angle of the inner ring;θjis the tilting angle of thejth roller at the azimuth angle φj.The classic slice method is applied to handle the contact issue between the crowned roller and raceways,where the roller is sliced toNsections along the axial direction,w(w=Ls/N)is the slice width,andLsis the effective length of the roller.

    The elastic deformations between thekth slice of thejth roller and raceways at the azimuth angle φjare expressed as:

    Fig.5 Forces and deformations between a roller and raceways.

    The normal contact forces between thekth slice of thejth roller and raceways are written as16

    whereA1is the elastic deformation coefficient,andA1=1.36η0.9;η is the combined elasticity modulus of two contact bodies.

    The normal contact forces between thejth roller and raceways are expressed as:

    The additional moments caused by normal contact forces between thejth roller and raceways are written as

    According to Eq.(5),the oil drag forces of thejth roller are written as

    where μjis the drag coefficient of oil film(see Ref.17 for detailed equation).

    The additional moments caused by additional moments are depicted as

    3.1.2.Forces between a roller and the cage

    The forces between a roller and the cage pocket are mainly considered as the hydrodynamic pressure between the roller and the beam of the cage pocket as well as elastic impact force due to contact deformation,which are also computed by the‘slice method”.The relationship of displacement and deformation between the roller and the cage pocket is shown in Fig.6. {opj;xpj,ypj,zpj} is the coordinate system of the cage pocket center,whose original pointopjcoincides with the geometric center of the cage pocket,thexpjaxis is parallel to the rotation axis of the cage,and theypjaxis is along the radial direction when thezpjaxis is along the circumferential direction of the pocket center.Every cage pocket has its own local coordinate system which moves and rotates with the cage.In Fig.6,βjis the skewing angle of thejth roller,Δzcjis the displacement of the roller geometric center under the coordinate system of the cage pocket center.

    The contact forceQcjand additional momentMcjbetween thejth roller and the cage bar are expressed as:

    The friction force between thejth roller and the cage bar is expressed as follows:

    where μcjis the friction coefficient between the roller and the cage pocket(see Ref.18 for detailed equation).

    3.1.3.Forces between the cage and the guide surface

    When the cage is guided by the outer ring,the interaction between the guide surface of the ring and the cylindrical surface of the cage is mainly produced by the hydrodynamic effect of lubricating oil.Since the acting surface of the guide rib and the cylindrical surface of the cage are small and mutually sliding,the infinitely short bearing theory is applied in this paper to determine the acting forcesFcy,Fczand momentMcxbetween the cage and the guide surface of the ring(see Ref.19 for detailed equation).

    Fig.6 Forces and deformations between a roller and the cage.

    3.1.4.Surface resistance and end-surface of the cage

    For a high-speed cylindrical roller bearing with a triloberaceway,cage rotation makes the external sides of the cage affected by the resistance from surrounding air/oil mist mixture shearing.See Ref.19 for the calculation method of the retardation torqueMcoof the cage surface and sides.

    3.2.Dynamics differential equations of the bearing

    As seen in Fig.7,dynamics differential equations of thejth roller are expressed as:

    When the bearing is working,the cage is simultaneously applied by the collision force of rollers,the guiding force of the outer ring,and the combined resistance of air/oil mist mixture to both cage ends and its surface,as shown in Fig.8,where ωrjis the spin angular velocity of thejth roller,and ωmis the actual speed of the cage.

    Fig.7 Schematic diagram of roller forces.

    As seen in Fig.8,dynamics differential equations of the cage are expressed as:

    Dynamics differential equations of the inner ring are expressed as:

    Fig.8 Schematic diagram of cage forces.

    Fig.9 Solution procedure of dynamics differential equations.

    3.3.Solution procedure of dynamics differential equations

    The GSTIFF integer algorithm is based on backward difference formulae(BDF)and is a variable-step and multi-step integrator,20which is applied to solve dynamics differential Eq.(16)of a cylindrical roller bearing with a trilobe-raceway in this study.The solution procedure of dynamics differential equations is shown in Fig.9.

    (1)The solution duration,initial step,and convergence error of dynamics differential equations are set,firstly.In this study,the solution duration was set to 0.5 s,the initial step was set to 0.00001 s,and the convergence error was set to 1×10-3.

    (2)According to initial estimated values of a bearing component’s position and motion constraint,initial conditions of dynamics differential equations,namely,relative positions and motion vectors of various components are obtained by quasi-static analysis.

    (3)Given initial conditions of dynamics differential equations,forces and moments applied on bearing components are obtained according to the expressions in Section 3.1.

    (4)The dynamics differential equations are solved by GSTIFF,and then positions and motion vectors of bearing components are obtained.

    (5)Verify whether the error meets the convergence error set in Step(1).If yes,continue next solution after getting the outputs of motion parameters including displacements,velocities,and accelerated speeds of the inner ring,cage,and roller.If no,choose a smaller step value and repeat Step(4)until the solving error meets the convergence error.

    4.Analysis on cage slip characteristics and basic rating life

    In this paper,a certain type of bearing is used to study the impacts of structural parameters and tolerance of the outer raceway,working conditions,and the outer ring installation method on cage slip characteristics and basic rating life.Main parameters of the bearing are shown in Table 1.The inner ring,outer ring,and roller are made of 8Cr4Mo4V,and the cage is made of 40CrNiMo,with a fixed outer ring and a cage guided by the outer ring.

    The cage slip ratioScand the basic rating lifeL10of a cylindrical roller bearing with a trilobe-raceway are defined as follows:

    4.1.Impact of geometric parameters on cage slip characteristics and rating life

    As known from Ref.1,the raceway contour of the outer ring after machining can be determined by the locations of high and low points,as shown in Fig.10.Therefore,this paper takes the low-radiusReminand the D-value δeas the geometric parameters to analyze their impacts on cage slip characteristics.

    4.1.1.Impact of the low-radius on the cage slip ratio and rating life

    Fig.11(a)shows the influence of the low-radius on the cage slip ratio when the rotating speed of the bearing is set to 45000 r/min,the radial force is 800 N,and the D-value is 0.098 mm.As seen in Fig.11(a),with the increase of the low-radius,the cage slip ratio increases and the amplitude grows.Lhis the basic rating life.This is because when the low-radius increases,as shown in Fig.12,the number of loaded rollers and the contact forces between rollers and the inner ring reduce,so the traction forces between rollers and the inner ring decrease,and the drop of the roller pushing force finally leads to a rise of the cage slip ratio.Meanwhile,with a small low-radius,the preloads applied on rollers from the outer raceway are quite large that the bearing rating life declines sharply,as shown in Fig.11(b).Consequently,there must be a reasonable low-radius that can avoid early failure caused by cage slip and also meet the requirement of the bearing rating life.

    Table 1 Main parameters of a bearing.

    4.1.2.Impact of the D-Value on the cage slip ratio and rating life

    Fig.13(a)shows the impact of the D-value on the cage slip ratio and rating life of a bearing when the rotating speed is 45000 r/min,the radial force is 800 N,and the low-radius is 34.8625 mm.As shown in Fig.13(a),the cage slip ratio rises linearly with the increase of the D-value.When the D-value increases,as shown in Fig.14,the loaded zones of the bearing become larger,and there are more rollers pushing the cage forward because of the increased number of loaded rollers,thus the cage slip ratio reduces,and the bearing rating life also decreases,as shown in Fig.13(b).Therefore,a smaller D-value is preferred on condition that the rating life of the bearing is guaranteed.

    4.2.Impact of working conditions on cage slip characteristics and rating life

    4.2.1.Impact of the bearing rotating speed on the cage slip ratio

    Assume that Fig.15 shows the impact of the inner ring rotating speed on the cage slip ratio when the radial force of the bearing is 800 N,the D-value is 0.098 mm,and the lowradius is 34.8600,34.8625,34.8650,and 34.8675 mm,respectively.

    As shown in Fig.15,the cage slip ratio rises with the acceleration of the bearing rotating speed,and the amplitude grows with the increase of the low-radius.The reason lies in the increased rollers centrifugal forces by the acceleration of the rotation speed,and diminished contact forces between rollers and the inner ring shown in Fig.16 resulting in a rise of the cage slip ratio by inadequate traction forces on rollers by the inner ring.When the low-radius is small,the preloads on rollers from the outer raceway are very effective that the centrifugal forces of rollers are relatively tiny compared to the preloads.Therefore,bearings with a smaller low-radius tend to diminish their sensitivity of cage slip to the rotating speed variation during operation.

    Fig.10 Schematic diagram of trilobe-raceway contours.

    Fig.11 Effects of the low-radius on the cage slip ratio and basic rating life.

    Fig.12 Effect of the low-radius on load distribution.

    4.2.2.Impact of the radial force on the cage slip ratio and rating life

    Fig.17 shows the relationship between the radial force and the cage slip ratio when the rotating speed of a bearing is 45000 r/min,the D-value is 0.098 mm,and the low-radius is 34.8600,34.8625,34.8650,and 34.8675 mm,respectively.

    As shown in Fig.17,the cage slip ratio remains unchanged with the increase of the radial force when the low-radius is 34.8600 mm.When the low-radiusis34.8625 mm and 34.8650 mm,with the increase of the radial force,the cage slip ratio firstly remains unchanged,and then decreases after the radial force exceeds a threshold value.The cage slip ratio will decrease with the increase of the radial force when the lowradius is 34.8675 mm.The reason lies in the fact that the loaded zone of the bearing and load distribution are mainly determined by the preloads from the raceway when the lowradius is small,with unchanged raceway geometric parameters and preloads to the roller by the raceway.When the low-radius increases and the radial force exceeds a threshold value,the radial force has a greater impact on bearing load distribution than preloads by the raceway.As Fig.18 shows,with the increases of the radial force,the contact forces between rollers and the inner ring as well as the number of loaded rollers,the cage slip ratio decreases as a result.Fig.17 also shows that the larger the low-radius is,the more evident the influence of the radial force on the cage slip ratio is.In other words,bearings with a smaller low-radius tend to weaken their sensitivity of cage slip when they are working under a variable radial force.

    4.3.Impact of the outer ring installation angle on cage slip characteristics and rating life

    As shown in Fig.1,when the outer ring is fixed and the direction of the radial force goes through the low point of a raceway contour,we define the installation angle at this time as 0°.Installation angle is defined as φewhen the outer ring is rotated clockwise to an angle φe.When φe=60°,i.e.,the installation angle is 60°,the radial force goes through the high point of a raceway contour,and the direction of the radial force changes over one period.

    Fig.13 Effects of the D-value on the cage slip ratio and basic rating life.

    Fig.19(a)shows the impact of the installation angle on the cage slip ratio when the rotating speed of the bearing is 45000 r/min,the radial force is800 N,the D-valueis 0.098 mm,and the low-radius is set to 34.8600,34.8625,34.8650,and 34.8675 mm,respectively.As shown in Fig.19,when the installation angle increases,the cage slip ratio remains unchanged and the bearing rating life rises up a little.This is because of the tiny dragging forces on the cage by rollers,with almost unchanged contact forces between rollers and the inner ring due to fixed raceway geometric parameters and radial force.As the installation angle changes from 0°to 60°,the loaded zone of the bearing as shown in Fig.20 becomes larger,and the dynamic equivalent load decreases with the increased number of loaded rollers,leading to a slight rise of the bearing rating life.Hence,when installing a cylindrical roller bearing with a trilobe-raceway,if possible,make the direction of the radial force go through the high point of the raceway contour to prolong the bearing rating life.

    Fig.14 Effect of the D-value on load distribution.

    Fig.15 Effect of the inner ring rotating speed on the cage slip ratio.

    Fig.16 Effect of the rotating speed on load distribution.

    Fig.17 Effect of the radial force on the cage slip ratio.

    Fig.18 Effect of the radial force on load distribution.

    4.4.Design tolerance of the trilobe-raceway

    4.4.1.Determination of the tolerance of raceway radius parameters

    Fig.21(a)and(b)shows the cage slip ratio and bearing rating life at different D-values when the rotating speed is 45000 r/min,the radial force is 800 N,and the low-radius is 34.8600,34.8625,34.8650,and 34.8675 mm,respectively.As shown in Fig.21,when the low-radius is smaller than 34.8650 mm,the cage slip ratio and bearing rating life rise with the increase of the D-Value,just like the above-mentioned analysis.When the low-radius is 34.8675 mm,with the increase of the D-value,the cage slip ratio and bearing rating life increase at first and then remain unchanged afterwards.This is because when the low-radius is small,preloads provided by the outer raceway is very small;when the D-value exceeds a threshold value,the area of the loaded zone and load distribution are determined by the applied radial force which is set to be a constant value.Therefore,the cage slip ratio remains unchanged.

    Fig.21(a)also shows that compared to the D-value,the low-radius has a greater impact on the cage slip ratio.Since the low-radius can affect not only the number of loaded rollers,but also the contact forces between rollers and the inner ring,once the low-radius is determined,the D-value has limited influence on the contact forces between rollers and the raceway within the loaded area;instead,it affects the cage slip ratio mostly by changing the number of loaded rollers through its impact on the loaded area.Therefore,to design a cylindrical roller bearing with a trilobe-raceway,it is recommended to determine firstly the range of the low-radius based on the demands of the cage slip ratio and bearing rating life,and then the range of the D-value according to the selected low-radius.

    As for the bearing concerned in this paper,the cage slip ratio under 10%is acceptable in actual practice.In order to meet the bearing rating life requirement of 300 h,the lowradius should be within the range of34.8625 mm to 34.8650 mm.When the D-value is less than 0.095 mm,a bearing with a low-radius as 34.8625 mm can no longer meet the rating life requirement.When the D-value is larger than 0.105 mm,the cage slip ratio of a bearing with a low-radius of 34.865 mm will be over 10%.Therefore,the D-value should be within the range of 0.095–0.105 mm.To reduce the influences of the bearing rotating speed and radial force variation during operation,the low-point diameter of the outer raceway contour can be designed as 69.725+00.005mm,and the D-value can be designed as 0.095+00.01mm.

    Fig.19 Effects of the installation angle on the cage slip ratio and basic rating life.

    Fig.20 Effect of the installation angle on load distribution.

    4.4.2.Impact of the base circle roundness on the cage slip ratio

    Fig.21 Effects of the D-value on the cage slip ratio and basic rating life at different low-radius values.

    Fig.22 shows the impact of the roundness of the base circle where the radius of the lowest point of the raceway contour locates on the cage slip ratio and bearing rating life,when the rotating speed is 45000 r/min,the radial force is 800 N,and the low-radius is 34.8625 mm and 34.8650 mm,respectively.As shown in Fig.22,when the roundness increases,the cage slip ratio increases while the bearing rating life reduces.This is because when the roundness of the base circle increases,the radii of three low points of the raceway contour are no longer the same,so are those of three high points.Therefore,the radius of the lowest point decreases when the difference value between the highest and lowest points of the raceway contour increases.Contact forces between rollers and the inner ring increase in areas near the lowest points but decrease in other areas.The increasing amplitude of the total bearing load caused by the former is smaller than the decreasing amplitude of the total bearing load caused by the latter.As a result,the cage slip ratio tends to increase and the bearing rating life decreases gradually.

    As for the bearing concerned in this paper,when the roundness of the base circle at low points is larger than 0.006,a bearing with a low-radius of 69.725 mm and a D-value of 0.095 mm can no longer meet the rating life requirement,so the roundness of the base circle at low points shall be designed as 0.006.

    Fig.22 Effects of the roundness on the cage slip ratio and basic rating life at different low-radius values.

    5.Comparison validation

    Experimental results in Ref.21 are adopted to verify the validity of the dynamic model built in this paper.Main parameters of a test bearing are shown in Table 2.

    Fig.23 shows the comparison of the cage slip ratio between experimental and simulation results when the applied radial force on the tested bearing is 800 N,and the rotating speed is varying from 15000 to 52000 r/min.As shown in Fig.23,the cage slip ratio rises with the acceleration of the bearing rotating speed,so simulation results are generally in line with experimental results,and the error is not more than 9%.When operating at lower speed,simulation results are smaller than experimental results,which may be caused by different lubricating oil used in the dynamic model and the experiment.In addition,when working at a higher speed,simulation results are larger than experimental results,because the dynamic model built in this paper does not include the impact of temperature on bearing geometric parameters.When the rotating speed is high,the radial internal clearance decreases because of the large heat productivity;hence,the cage slip ratio in the experiment is slightly smaller than that in the simulation.

    Table 2 Main parameters of a test bearing.

    Fig.23 Comparison between experimental and simulation results.

    6.Conclusions

    (1)The cage slip ratio and bearing rating life of a cylindrical roller bearing with a trilobe-raceway reduce with the decrease of the low-radius and D-value,and the lowradius has a greater impact than the D-value.Therefore,on the premise of bearing rating life,a small low-radius and D-value are preferred,and the variation range of the low-radius should be determined before the determination of the D-value.

    (2)The cage slip ratio rises with the acceleration of the inner ring rotating speed and drops with the increase of the radial force,and the changing amplitude increases along with the increase of the low-radius.That’s to say,a bearing with a smaller low-radius tends to diminish its sensitivity of cage slip when working under a variable condition.

    (3)When the installation angle of a cylindrical roller bearing with a trilobe-raceway increases within one period,the cage slip ratio remains unchanged while the bearing rating life rises up a little.Hence,the direction of the radial force should go through the high-point radius of the outer raceway contour when installing the bearing to prolong its rating life.

    (4)With the increase of the roundness of the base circle at three low points of the raceway contour,the cage slip ratio rises and the bearing rating life decreases.Thus,the roundness of the base circle should be reduced to the greatest extent under attainable machining precision.

    Acknowledgements

    This research was financially co-supported by the National Natural Science Foundation of China(U1404514),Henan Outstanding Person Foundation in China(144200510020),and the Collaborative Innovation Center of Major Machine Manufacturing in Liaoning,China.

    1.Wang LQ,Chen GC.Pre-deformation machining technology of the noncircular raceway of aero roller bearings.Chin J Mech Eng2005;41(9):223–7.

    2.Harris TA.Analytical method to predict skidding in high speed roller bearings.ASLE Tran1966;9:229–41.

    3.Poplawski JV.Slip and cage forces in a high-speed roller bearing.J Lubr Tech1972;94(2):143–50.

    4.Cavallaro G,Nelias D.Analysis of high-speed inter-shaft cylindrical roller bearing with flexible rings.Tribol T2005;61(11):38–244.

    5.Takafumi Y,Yasuyoshi T,Hisao M,Masaaki S.Analysis of cage slip in cylindrical roller bearings considering non-newtonian behavior and temperature rise of lubricating oil.Tribol T2008;53(6):647–62.

    6.Arthanari S,Marappan R.Experimental analysis of factors influencing the cage slip in cylindrical roller bearing.Int J Adv Manuf Technol2011;53(5):635–44.

    7.Takabi J,Khonsari MM.On the influence of traction coefficient on the cage angular velocity in roller bearings.Tribol T2014;57(5):793–805.

    8.Li JB,Wu LF.Dynamic analysis of high-speed roller bearings.Chin J Aeronaut1992;13(12):625–32.

    9.Hu X,Luo GH,Gao DP.Quasi-static analysis of cylindrical roller inter-shaft bearing.Chin J Aeronaut2006;21(6):1069–74.

    10.Hu X,Luo GH,Gao DP.Performance analysis of aero-engine inter-shaft bearing.Chin J Aeronaut2007;22(3):439–43.

    11.Tu WB,Shao YM.An analytical model to investigate skidding in rolling element bearings during acceleration.J Mech Sci Technol2012;26(8):2451–8.

    12.Chen W,Li JN,Zhang LB.Skidding analysis of high speed rolling bearing considering whirling of bearing.Chin J Mech Eng2013;49(6):38–43.

    13.Deng SE,Gu JF,Cui YC,Sun CY.Analysis on dynamic characteristics of cage in high-speed cylindrical roller bearing.Chin J Aeronaut2014;29(1):492–8.

    14.Fan DP,Yin ZQ.Grinding technology of high speed anti-skidding bearings with noncircular raceway.Bearing1996;12:18–20.

    15.Lopatukhin I,Khaskel L.Deformation of thin-walled rings under the action of n concentrated and distributed forces.Int J Adv Manuf Technol2010;49:627–34.

    16.Harris TA.The effect of misalignment on the fatigue life of cylindrical roller bearings having crowned rolling members.J Tribol1969;91(2):294.

    17.Wang YS,Yang BY,Wang LQ.Investigation into the traction coefficient in elasto-hydrodynamic lubrication.Tribotest2004;11(2):113–24.

    18.Xiuhai Liu,Sier Deng.Dynamic stability analysis of cages in highspeed oil-lubricated angular contact ball bearings.Trans Tianjin Univ2011;17(1):20–7.

    19.Deng SE,Jia QY,Xue JX.Design principle of rolling bearings.2nd ed.Beijing:Standards Press of China;2014.p.113–24.

    20.Gear CW.Simultaneous numerical solution of differential-algebraic equations.IEEE Trans Circ Syst1971;18(1):89–95.

    21.Zuo JT.Study on a certain bearings of aircraft engine[dissertation].Harbin:Harbin Institute of Technology;2008.

    午夜福利,免费看| 日韩 欧美 亚洲 中文字幕| 亚洲精品国产区一区二| 丰满饥渴人妻一区二区三| 91老司机精品| www.精华液| 国产精品永久免费网站| 两个人看的免费小视频| 免费观看精品视频网站| 国产成人免费无遮挡视频| 高潮久久久久久久久久久不卡| 制服人妻中文乱码| av国产精品久久久久影院| 国产激情欧美一区二区| 国产成人免费无遮挡视频| netflix在线观看网站| 欧美国产精品va在线观看不卡| 水蜜桃什么品种好| 国产亚洲av高清不卡| 黑人巨大精品欧美一区二区mp4| 亚洲精品在线美女| 免费看a级黄色片| 亚洲av美国av| 操出白浆在线播放| 色播在线永久视频| 天天影视国产精品| 欧美精品人与动牲交sv欧美| 亚洲avbb在线观看| 久久香蕉精品热| 日韩欧美一区二区三区在线观看 | 91国产中文字幕| 国产精品亚洲av一区麻豆| 天天躁日日躁夜夜躁夜夜| 好男人电影高清在线观看| 亚洲精品美女久久av网站| 欧美精品一区二区免费开放| 国产精品电影一区二区三区 | 国产麻豆69| 国产成人免费无遮挡视频| 欧美日韩亚洲高清精品| 黑丝袜美女国产一区| 日本wwww免费看| 精品高清国产在线一区| 亚洲欧美色中文字幕在线| 香蕉丝袜av| 操美女的视频在线观看| 岛国在线观看网站| 国产免费av片在线观看野外av| 欧美精品人与动牲交sv欧美| 国产精品欧美亚洲77777| 曰老女人黄片| 亚洲第一青青草原| 欧美日韩一级在线毛片| 天天躁日日躁夜夜躁夜夜| 男女午夜视频在线观看| 最近最新中文字幕大全免费视频| 欧美国产精品一级二级三级| 超碰97精品在线观看| 在线免费观看的www视频| 一级片免费观看大全| 欧美国产精品一级二级三级| 老司机靠b影院| www.999成人在线观看| 中文欧美无线码| 淫妇啪啪啪对白视频| 国产无遮挡羞羞视频在线观看| 99精品欧美一区二区三区四区| 成年人黄色毛片网站| 国产成人av激情在线播放| 午夜两性在线视频| 日韩熟女老妇一区二区性免费视频| 欧美日韩福利视频一区二区| 亚洲成人手机| 国产成人系列免费观看| 久久久久久久久久久久大奶| 91老司机精品| 国产亚洲精品第一综合不卡| av网站免费在线观看视频| www.精华液| 一级黄色大片毛片| 国产91精品成人一区二区三区| 中文欧美无线码| 亚洲精品av麻豆狂野| 大型黄色视频在线免费观看| 亚洲欧美一区二区三区黑人| av天堂久久9| 丰满饥渴人妻一区二区三| www.999成人在线观看| 在线观看免费午夜福利视频| 精品高清国产在线一区| 亚洲片人在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 色播在线永久视频| 巨乳人妻的诱惑在线观看| 国产精品久久久久久人妻精品电影| www.精华液| 欧美成人午夜精品| 女人久久www免费人成看片| 欧美激情久久久久久爽电影 | 日本vs欧美在线观看视频| 国产片内射在线| 又黄又粗又硬又大视频| 丰满饥渴人妻一区二区三| 免费黄频网站在线观看国产| 人妻一区二区av| 免费女性裸体啪啪无遮挡网站| 黄片大片在线免费观看| 色婷婷av一区二区三区视频| 91大片在线观看| 成人永久免费在线观看视频| 麻豆乱淫一区二区| 亚洲欧美色中文字幕在线| 亚洲,欧美精品.| 久久亚洲真实| 在线观看日韩欧美| 国内久久婷婷六月综合欲色啪| 日韩免费高清中文字幕av| 婷婷成人精品国产| 亚洲视频免费观看视频| 国产亚洲精品久久久久5区| 成熟少妇高潮喷水视频| 精品国产超薄肉色丝袜足j| 啦啦啦 在线观看视频| 日本wwww免费看| 欧美性长视频在线观看| 久久精品国产清高在天天线| 国产成人av教育| 精品国产国语对白av| 国产一区在线观看成人免费| 后天国语完整版免费观看| 亚洲成人手机| 国产欧美日韩一区二区三| 制服诱惑二区| 搡老熟女国产l中国老女人| 精品视频人人做人人爽| 日韩有码中文字幕| 精品国产国语对白av| 亚洲中文字幕日韩| 久久中文字幕一级| 国产成人免费观看mmmm| 欧美乱色亚洲激情| 精品国产美女av久久久久小说| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美激情久久久久久爽电影 | 精品人妻在线不人妻| 久久精品亚洲av国产电影网| 欧美人与性动交α欧美软件| 天堂动漫精品| 大码成人一级视频| 欧美亚洲日本最大视频资源| 天堂中文最新版在线下载| 亚洲avbb在线观看| 日韩制服丝袜自拍偷拍| 满18在线观看网站| www.熟女人妻精品国产| ponron亚洲| 免费观看人在逋| 法律面前人人平等表现在哪些方面| 久久久精品区二区三区| 黄色视频,在线免费观看| 欧美另类亚洲清纯唯美| 大型黄色视频在线免费观看| 满18在线观看网站| 免费观看精品视频网站| 一区福利在线观看| 欧美日韩福利视频一区二区| 黄色成人免费大全| 极品教师在线免费播放| 黑人欧美特级aaaaaa片| 啦啦啦免费观看视频1| 亚洲人成电影观看| 国产1区2区3区精品| 欧美日韩黄片免| 午夜福利乱码中文字幕| 国产欧美日韩精品亚洲av| 大香蕉久久成人网| 午夜成年电影在线免费观看| 中文字幕高清在线视频| 高清黄色对白视频在线免费看| 久久香蕉国产精品| 黄片大片在线免费观看| 精品欧美一区二区三区在线| 人人妻人人添人人爽欧美一区卜| 一级a爱视频在线免费观看| 成人国产一区最新在线观看| 日本精品一区二区三区蜜桃| 欧美日韩视频精品一区| 国产成人av激情在线播放| www.自偷自拍.com| 91av网站免费观看| av网站在线播放免费| 欧美精品亚洲一区二区| 一本一本久久a久久精品综合妖精| 欧美亚洲 丝袜 人妻 在线| 在线av久久热| 精品国内亚洲2022精品成人 | 下体分泌物呈黄色| 老熟女久久久| 亚洲av电影在线进入| 99热国产这里只有精品6| 国产人伦9x9x在线观看| 一本综合久久免费| av免费在线观看网站| 国产精品久久久久久精品古装| 日韩三级视频一区二区三区| а√天堂www在线а√下载 | 黑人巨大精品欧美一区二区蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人精品无人区| 狂野欧美激情性xxxx| 国产激情欧美一区二区| 妹子高潮喷水视频| 亚洲第一av免费看| 久久精品人人爽人人爽视色| 午夜福利在线观看吧| 高清av免费在线| 无限看片的www在线观看| 国产男靠女视频免费网站| 精品国产一区二区三区久久久樱花| 18禁黄网站禁片午夜丰满| 一区二区三区国产精品乱码| 身体一侧抽搐| 日本精品一区二区三区蜜桃| 如日韩欧美国产精品一区二区三区| 看免费av毛片| 中出人妻视频一区二区| 亚洲一区中文字幕在线| 国产av一区二区精品久久| 亚洲五月天丁香| 伦理电影免费视频| 国产不卡一卡二| 国产av又大| 两性午夜刺激爽爽歪歪视频在线观看 | 人妻丰满熟妇av一区二区三区 | 亚洲aⅴ乱码一区二区在线播放 | av中文乱码字幕在线| 国产黄色免费在线视频| 久久久久国产精品人妻aⅴ院 | 超碰成人久久| 黑丝袜美女国产一区| 黑人欧美特级aaaaaa片| 国产精品国产av在线观看| 在线av久久热| 午夜免费成人在线视频| 18禁黄网站禁片午夜丰满| av免费在线观看网站| 叶爱在线成人免费视频播放| av国产精品久久久久影院| 一区二区日韩欧美中文字幕| 久久久久视频综合| videos熟女内射| 国产熟女午夜一区二区三区| 在线观看www视频免费| 黑丝袜美女国产一区| 国产成人av激情在线播放| 国产精品香港三级国产av潘金莲| 欧美丝袜亚洲另类 | 亚洲熟女毛片儿| 少妇裸体淫交视频免费看高清 | 激情在线观看视频在线高清 | 国产精品免费视频内射| 黑丝袜美女国产一区| 黄色a级毛片大全视频| 午夜福利欧美成人| 亚洲免费av在线视频| 成人av一区二区三区在线看| 丝袜人妻中文字幕| 亚洲专区国产一区二区| 国产精品久久久久久人妻精品电影| 黑人操中国人逼视频| 亚洲男人天堂网一区| 人妻久久中文字幕网| 好男人电影高清在线观看| 黄片大片在线免费观看| 在线观看免费视频日本深夜| 中文字幕人妻丝袜制服| 久久青草综合色| 国产在线观看jvid| 亚洲一区二区三区欧美精品| 亚洲欧洲精品一区二区精品久久久| 欧美日韩精品网址| 捣出白浆h1v1| 女警被强在线播放| 王馨瑶露胸无遮挡在线观看| 久久久国产欧美日韩av| 母亲3免费完整高清在线观看| 成人18禁在线播放| 飞空精品影院首页| 人人妻人人澡人人爽人人夜夜| 亚洲成国产人片在线观看| 又黄又爽又免费观看的视频| 91精品三级在线观看| 操美女的视频在线观看| 9热在线视频观看99| 中文字幕人妻丝袜一区二区| 色婷婷av一区二区三区视频| 一二三四在线观看免费中文在| www.999成人在线观看| 亚洲欧美激情综合另类| 亚洲精品一卡2卡三卡4卡5卡| 99精品欧美一区二区三区四区| 99热网站在线观看| 亚洲欧美日韩另类电影网站| 久久香蕉国产精品| 亚洲综合色网址| 日韩制服丝袜自拍偷拍| 三级毛片av免费| 不卡一级毛片| 成人黄色视频免费在线看| 在线观看免费视频日本深夜| 一边摸一边抽搐一进一小说 | 久久性视频一级片| 亚洲熟女毛片儿| 最近最新免费中文字幕在线| 九色亚洲精品在线播放| 夜夜爽天天搞| 欧美日韩精品网址| 国产成人精品久久二区二区91| 国产在视频线精品| 久久中文字幕一级| 欧美人与性动交α欧美软件| 天堂中文最新版在线下载| 午夜91福利影院| 国产精品.久久久| 淫妇啪啪啪对白视频| 18禁国产床啪视频网站| 国产日韩一区二区三区精品不卡| 免费在线观看完整版高清| 国产高清videossex| 在线观看免费视频网站a站| 成年女人毛片免费观看观看9 | 一进一出抽搐动态| 嫩草影视91久久| 久久香蕉激情| www日本在线高清视频| 老熟女久久久| 国产野战对白在线观看| 亚洲一区二区三区欧美精品| 韩国av一区二区三区四区| 国产精品久久久人人做人人爽| 99久久精品国产亚洲精品| 建设人人有责人人尽责人人享有的| 成人亚洲精品一区在线观看| 国内久久婷婷六月综合欲色啪| 12—13女人毛片做爰片一| 免费女性裸体啪啪无遮挡网站| 国产99久久九九免费精品| 一a级毛片在线观看| 欧美人与性动交α欧美精品济南到| 免费在线观看视频国产中文字幕亚洲| 亚洲精品美女久久久久99蜜臀| 亚洲精品国产一区二区精华液| 老熟女久久久| 免费看十八禁软件| 亚洲色图综合在线观看| 久久国产乱子伦精品免费另类| 午夜精品在线福利| 日本一区二区免费在线视频| 中文字幕精品免费在线观看视频| 国产主播在线观看一区二区| 中亚洲国语对白在线视频| 国产精品一区二区免费欧美| 亚洲av成人不卡在线观看播放网| 精品免费久久久久久久清纯 | 成人18禁高潮啪啪吃奶动态图| 亚洲五月色婷婷综合| 成年人免费黄色播放视频| 国产成人一区二区三区免费视频网站| 看免费av毛片| 国产一区有黄有色的免费视频| 1024视频免费在线观看| 国产在线精品亚洲第一网站| 18禁观看日本| 久久人妻熟女aⅴ| 精品欧美一区二区三区在线| 国产视频一区二区在线看| 天堂√8在线中文| 亚洲欧美激情综合另类| 女性被躁到高潮视频| 国产91精品成人一区二区三区| 久久久久视频综合| 国产精品久久久久久精品古装| 香蕉丝袜av| 人人妻人人添人人爽欧美一区卜| 精品国产一区二区三区四区第35| 别揉我奶头~嗯~啊~动态视频| 高清欧美精品videossex| 99精品欧美一区二区三区四区| 成人三级做爰电影| 精品电影一区二区在线| 久久久久国内视频| 亚洲专区字幕在线| 久久久水蜜桃国产精品网| 国产精品欧美亚洲77777| av中文乱码字幕在线| 久久九九热精品免费| 欧美不卡视频在线免费观看 | 免费黄频网站在线观看国产| 欧美日韩一级在线毛片| 在线观看免费视频日本深夜| 18禁美女被吸乳视频| 男男h啪啪无遮挡| 深夜精品福利| 真人做人爱边吃奶动态| 精品国产一区二区三区久久久樱花| 一本大道久久a久久精品| 在线看a的网站| 欧美精品亚洲一区二区| 香蕉国产在线看| 纯流量卡能插随身wifi吗| 亚洲人成伊人成综合网2020| 亚洲第一av免费看| 伦理电影免费视频| 亚洲熟妇熟女久久| 天天躁日日躁夜夜躁夜夜| 亚洲av第一区精品v没综合| 欧美不卡视频在线免费观看 | 91老司机精品| 国产免费现黄频在线看| 又黄又粗又硬又大视频| 久久天躁狠狠躁夜夜2o2o| 一边摸一边做爽爽视频免费| 欧美日韩一级在线毛片| 久久国产乱子伦精品免费另类| 成人三级做爰电影| 黄片小视频在线播放| 人人妻人人添人人爽欧美一区卜| 国产精品一区二区在线观看99| 亚洲专区字幕在线| 电影成人av| 一区二区三区激情视频| 免费在线观看黄色视频的| 精品亚洲成a人片在线观看| 可以免费在线观看a视频的电影网站| avwww免费| 丰满迷人的少妇在线观看| 19禁男女啪啪无遮挡网站| 在线观看舔阴道视频| 欧美日韩乱码在线| 女同久久另类99精品国产91| 国产欧美日韩一区二区精品| videos熟女内射| 美女扒开内裤让男人捅视频| 热99久久久久精品小说推荐| 久久青草综合色| 黄色视频,在线免费观看| 亚洲国产精品合色在线| 十八禁人妻一区二区| ponron亚洲| 亚洲一区高清亚洲精品| 两个人免费观看高清视频| 亚洲在线自拍视频| 午夜福利在线免费观看网站| 久久久久国产精品人妻aⅴ院 | 久久精品国产a三级三级三级| 在线av久久热| 91麻豆精品激情在线观看国产 | 久久久久久久国产电影| 国产单亲对白刺激| 在线观看免费视频日本深夜| 18禁裸乳无遮挡免费网站照片 | 中亚洲国语对白在线视频| 精品一区二区三区视频在线观看免费 | 大型av网站在线播放| 女同久久另类99精品国产91| 99热国产这里只有精品6| 欧美av亚洲av综合av国产av| av天堂久久9| av有码第一页| 亚洲va日本ⅴa欧美va伊人久久| 最近最新免费中文字幕在线| 纯流量卡能插随身wifi吗| 亚洲人成伊人成综合网2020| av在线播放免费不卡| 久久影院123| 五月开心婷婷网| 老司机亚洲免费影院| 在线观看免费日韩欧美大片| 日本五十路高清| 校园春色视频在线观看| 一二三四社区在线视频社区8| 天天躁日日躁夜夜躁夜夜| 日韩人妻精品一区2区三区| 欧美日韩黄片免| 最新在线观看一区二区三区| 1024香蕉在线观看| 91在线观看av| 久久久久国内视频| 国产精品影院久久| 国内久久婷婷六月综合欲色啪| 久久久水蜜桃国产精品网| 国产精品久久久久成人av| 欧美日韩成人在线一区二区| 日本五十路高清| 成人18禁在线播放| 五月开心婷婷网| 色老头精品视频在线观看| 青草久久国产| 午夜影院日韩av| 中文字幕最新亚洲高清| 国产男女超爽视频在线观看| 又紧又爽又黄一区二区| 黄色片一级片一级黄色片| 一a级毛片在线观看| 午夜久久久在线观看| 国产av一区二区精品久久| 最新美女视频免费是黄的| 欧美黑人欧美精品刺激| 久久人妻福利社区极品人妻图片| 免费在线观看亚洲国产| 午夜视频精品福利| 69精品国产乱码久久久| 国产成人欧美| 国产欧美日韩一区二区精品| 丝袜在线中文字幕| av福利片在线| 嫩草影视91久久| 别揉我奶头~嗯~啊~动态视频| 日本黄色视频三级网站网址 | 欧美精品啪啪一区二区三区| 久久99一区二区三区| 国产成人啪精品午夜网站| 王馨瑶露胸无遮挡在线观看| 天天添夜夜摸| 女人爽到高潮嗷嗷叫在线视频| 精品国产一区二区三区久久久樱花| 9191精品国产免费久久| 午夜影院日韩av| 亚洲第一青青草原| 黄片播放在线免费| 天堂√8在线中文| 美女福利国产在线| 亚洲人成77777在线视频| 精品欧美一区二区三区在线| 国产在线观看jvid| 欧洲精品卡2卡3卡4卡5卡区| 亚洲欧美色中文字幕在线| 中出人妻视频一区二区| 国产欧美日韩一区二区三| 怎么达到女性高潮| 大陆偷拍与自拍| 国产乱人伦免费视频| 久久精品亚洲av国产电影网| 自线自在国产av| e午夜精品久久久久久久| 免费少妇av软件| 伦理电影免费视频| 99国产精品免费福利视频| 国产99久久九九免费精品| 国产午夜精品久久久久久| 怎么达到女性高潮| 一个人免费在线观看的高清视频| 超碰成人久久| 老鸭窝网址在线观看| 精品国产超薄肉色丝袜足j| 日本wwww免费看| 国产又爽黄色视频| 多毛熟女@视频| 欧美激情 高清一区二区三区| 欧美久久黑人一区二区| 视频区欧美日本亚洲| 亚洲专区字幕在线| 一级,二级,三级黄色视频| 国产亚洲一区二区精品| 人人妻人人澡人人看| 国产野战对白在线观看| 日韩欧美在线二视频 | 欧美乱妇无乱码| 国产亚洲一区二区精品| 午夜精品久久久久久毛片777| 久久人人97超碰香蕉20202| 精品一区二区三区四区五区乱码| 亚洲精品在线观看二区| 亚洲第一欧美日韩一区二区三区| av在线播放免费不卡| 久久人人爽av亚洲精品天堂| 亚洲成av片中文字幕在线观看| 国产精品欧美亚洲77777| 淫妇啪啪啪对白视频| 精品人妻在线不人妻| 久久ye,这里只有精品| 一本一本久久a久久精品综合妖精| 啦啦啦在线免费观看视频4| 午夜福利在线免费观看网站| 99国产精品一区二区三区| 最近最新免费中文字幕在线| 国产无遮挡羞羞视频在线观看| 一进一出抽搐动态| 丁香欧美五月| 国产精品一区二区在线不卡| 成年人免费黄色播放视频| 日韩免费高清中文字幕av| 国产精品 国内视频| 国产成人欧美| 久久久久久亚洲精品国产蜜桃av| 在线观看免费视频网站a站| 久久人人爽av亚洲精品天堂| 91成年电影在线观看| 色老头精品视频在线观看| 桃红色精品国产亚洲av| 亚洲专区中文字幕在线| 首页视频小说图片口味搜索| 国产主播在线观看一区二区| 午夜免费鲁丝| 叶爱在线成人免费视频播放| 黄色成人免费大全| 在线观看日韩欧美| 91成人精品电影| 精品第一国产精品| 最近最新中文字幕大全免费视频| 国产深夜福利视频在线观看| 国产成人av激情在线播放| 亚洲精华国产精华精| 亚洲国产毛片av蜜桃av| 少妇猛男粗大的猛烈进出视频|