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

    Analysis of fretting wear behavior of unloading valve of gasoline direct injection high-pressure pump

    2022-05-10 10:07:50LiangLUYinpengXUMengruLIQilongXUEManyiZHANGLiangliangLIUZhongyuWU

    Liang LU ,Yin-peng XU ,Meng-ru LI? ,Qi-long XUE ,Man-yi ZHANG ,Liang-liang LIU ,Zhong-yu WU

    1School of Mechanical Engineering,Tongji University,Shanghai 201804,China

    2State Key Laboratory of Fluid Power&Mechatronic Systems,Zhejiang University,Hangzhou 310027,China

    3United Automotive Electronic Systems Co.,Ltd.,Shanghai 200120,China

    4Frontiers Science Center for Intelligent Autonomous Systems,Shanghai 201210,China

    Abstract:The high pressures in gasoline direct injection technology lead to structural damage in some hydraulic components,especially annular damage on the contact area of the valve ball and on the valve seat of the spherical unloading valve in the high-pressure pump.In previous study,the authors have analyzed the damage on the unloading valve and demonstrated that it is caused neither by static damage nor fatigue damage and have put forward the hypothesis of fretting wear.This paper is based on the establishment of the statically indeterminate structure of the unloading valve.The micro friction parameters (stress,friction coefficient,etc.) required for the numerical iterative calculation of fretting wear are calculated.In addition,based on the grid adaptive technology and a modified Archard wear model,the fretting wear is calculated quantitatively and is in good agreement with experimental results.Based on that verification,the wear laws of the valve ball and valve seat under the same hardness,different contact angles,and different assembly stresses,are analyzed in detail,and reasoned suggestions for the structural design and assembly design of the ball valve are given.

    Key words:Fretting wear behavior;Unloading valve;Experimental and numerical analyses;High pressure

    1 Introduction

    Compared with the traditional intake port injection gasoline engine,gasoline direct injection (GDI)technology can reduce fuel cost by 15%and hydrocarbon emission by 30%.Therefore,GDI technology has strong advantages in saving energy consumption and reducing exhaust pollution.High-pressure direct injection technology is one of the directions of development for direct injection vehicles (Zhao et al.,1999).The high-pressure fuel pump is the key to the whole highpressure fuel supply system.The spherical unloading valve,which will open and release high-pressure fuel when the fuel pressure in the high-pressure side of the system is too high,works as the safety valve in the high-pressure direct injection fuel supply system for automobiles.The spherical unloading valve plays a key role in building and maintaining the fuel pressure of the high-pressure fuel supply system.Responding to market demand for high-pressure direct injection systems,the fuel pressure of the high-pressure side of the spherical unloading valve is increased.Consequently,the maximum pressure of the low-pressure side is increased,which causes annular damage on the contact area of the valve ball and the valve seat.The damage causes a premature opening of the spherical valve when the fuel pressure of the high-pressure side is below 35 MPa.Thus,the pressure of the fuel supply system cannot be built and maintained effectively.The system cannot meet the requirements for normal work in respect either of security or of service life.

    Facing the problems of noise,fuel injection and structural strength caused by high pressure,scholars in relevant industries have made targeted analysis and research.Borg et al.(2012) proposed a currentshaping technique and a pulse-skipping control strategy to reduce noise by up to 5 dB in the low-frequency and mid-frequency range and an injector suspension method to reduce noise in the high-frequency range of the gasoline direct injection system.Sharma et al.(2019) investigated the combustion,noise,and vibration characteristics of three oxygenated test fuels(M10,M20,and G100).The results showed that gasohol generates lesser vibration and noise,which makes them suitable for use as oxygenated fuels and for partial replacement of conventional petroleum-based fossil fuels such as gasoline.Liu et al.(2021) built a test platform for ethanol port injection plus gasoline direct injection to explore the effects on the combustion and emissions of a spark ignition (SI) engine.Results clearly show that when the total fuel contains more ethanol,the ignition timing corresponding to the maximum torque is smaller than the ignition timing when the total fuel contains more gasoline.Qian et al.(2019)investigated the effects of combustion system parameters on the scavenging process and proposed an optimum structural shape for the combustion chamber.At present,however,there are few targeted studies on the structural damage of the various components in the high-pressure oil supply system.In this paper,the structural damage of the high-pressure pump unloading valve is analyzed.According to preliminary analysis,it is considered that the damage is caused by fretting wear(Lu et al.,2019).

    Friction problems of hydraulic fluid control components have always been an important research topic.Scholars have done much research on hydraulic pumps,hydraulic slide valves,and other components.Yin et al.(2017)developed a numerical model for predicting the removal rate of corrosive material and the profile evolution of wear in hydraulic spool valves and discussed the main influencing factors,including grain size,differential pressure,reel opening,and flow direction.Their model was validated by comparison with laboratory experiments and field cases.The design and optimization of the key friction pairs are always key and difficult problems in research on the axial piston pump.Xu et al.(2012,2013)established a simulation model based on oil film lubrication theory,analyzed the dynamic characteristics and pressure distribution of radial fretting,and discussed the relationship between radial fretting and bearing capacity,lubrication conditions,and wear.For the high wear rate and rapid loss of the surface modification layer of the piston/cylinder pair in axial piston pumps,Zhang et al.(2021) built a wear degradation model to obtain the variation in wear of the piston/cylinder pair with the operating time.The wear contour at the lowest wear rate was derived to guide the design of a pre-machined surface contour.In addition,additive manufacturing not only provides a new way to manufacture the products with complex internal channels,but also indicates a design revolution for existing products with regular internal channels that are currently limited by conventional fabrication methods (Zhang et al.,2020).For the poor fabricating quality caused by residual stress on the large overhang region when using one type of metal additive manufacturing technology—selective laser melting (SLM),Zhou et al.(2020) developed a new friction factor model to calculate pressure loss in a SLM fabricated channels.The friction loss of the laser powder bed fusion (L-PBF) fabricated fluid channels is greatly affected by the fabrication quality but remains unknown.Zhou et al.(2021)developed a comprehensive model to predict friction factors of the L-PBF fabricated fluid channels and proposed design guidelines of fluid channels by considering both fabrication quality and friction factors.

    Most of the early research on fretting wear was at the level of qualitative understanding.However,since the 1990s,quantitative experimental and simulation research on fretting wear has increased.Goryacheva et al.(2001) developed the analysis technology for the purely contact based wear problem and proposed an analysis method to simulate the fretting wear under the condition of partial sliding in 2D sliding.In this model,the Archard equation was applied to evaluate the wear in one cycle and the distance change in the sliding region.Then an iterative program was used to calculate the change of contact characteristics with the increase of wear cycles.To solve this kind of problem Johansson (1994) proposed a finite element method combined with the Archard equation to evaluate the changes of contact geometry and related contact pressure.?qvist (2001) proposed a numerical simulation study of the slight wear between the cylindrical steel roll and the steel plate in the reciprocating contact fit,obtained the simulation results of the wear profile,and established a good correlation between the numerical prediction and the experimentally measured wear profile over approximately one million cycles.Experiments on fretting wear are very costly,time-consuming,and laborious.The idea of finite element simulation provides a new approach of low cost and easily controlled parameter variables and boundary conditions for the study of fretting wear.The finite element method can provide corresponding research methods for different research objects.In recent years,many scholars have carried out wear research in the fields of car body,car wheel,and bionics (Jia et al.,2020;Sun and Xu,2020;Tao et al.,2021).In this paper,the finite element simulation software ABAQUS is used to study the fretting wear of the unloading valve,to seek the optimal parameters for the improvement of the anti-wear structure,and to provide a new solution to the structural problems caused by the high pressure of the automobile oil supply system.

    This paper is organized as follows.First,the statically indeterminate structure of the unloading valve is analyzed to acquire the essential parameters that would be needed in simulation.Because previous research has proved that the damage of the unloading valve is neither caused by static damage nor fatigue damage but by fretting wear (Lu et al.,2019),experiments to obtain the fretting wear characteristics of the unloading valve are conducted.To quantitatively calculate the wear amount of fretting wear by simulation,a modified Archard model is used together with the adaptive grid technology.A set of comparison results between simulation and experiment verify the rationality of the simulation modeling.Finally,the influence of the contact angle and the assembly stress is discussed in detail to study the law of fretting wear of the unloading valve.

    2 Model and experiment

    2.1 Principle of the unloading valve

    This paper concerns the spherical unloading valve in the high-pressure pump of a direct injection vehicle.The working principle of the high-pressure pump is shown in Fig.1.With the rotation of the cam,the piston rod of the oil cylinder moves downward,and the flow control valve controls the switch to the oneway valve to suck oil from the low-pressure end.With the upward movement of the piston rod,the oil cylinder supplies oil to the high-pressure end.The spherical unloading valve is an important part of the highpressure pump of a direct injection vehicle.As a safety valve,the unloading valve plays a role in building and maintaining the pressure.When the pressure at the high-pressure end is too high,it can release the oil,thus ensuring the normal operation of the whole highpressure pump.

    Fig.2 shows the structure and force state of the valve,where

    f

    represents the friction between valve ball and valve seat,

    N

    represents the normal pressure between valve ball and valve seat,and

    F

    represents the equivalent force of pressure

    p

    .The valve is stressed with an alternating pressure

    p

    of 0 to 45 MPa from the left side and a constant pressure of

    p

    from the right side.The valve is also initially stressed with a preload of the spring from the left side.When the pressure difference exceeds the preload of the spring,the spherical valve opens and high-pressure oil from the right-side flows to the left side to unload.Table 1 shows the structural parameters of the valve and oil pressure parameters.The contact area of the valve ball and valve seat is under an equivalent pressure of 5 to 50 MPa.

    Table 1 Parameters of structure and oil pressure

    2.2 Experiment design

    According to the above hydraulic principle and force state of the unloading valve,experiments were designed to study the fretting wear characteristics of the unloading valve.The contact area between the valve ball and the valve seat receives an alternating pressure of 5-50 MPa,which can be equivalent to a loading force of 6.90-78.75 N.To provide a high frequency alternating force,a piezoelectric ceramic was used as the key component of the experiment device.For the loading scheme of piezoelectric ceramics,the following factors need to be considered in the design process:

    (a)The unloading valve needs to be fixed during the experiment.

    (b)It is necessary to accurately locate and fix the relative position of piezoelectric ceramics and the unloading valve,so as to ensure the accuracy of rigid deformation on the force loading of the unloading valve.

    (c) A micro pressure sensor needs to be added between the piezoelectric ceramic and the unloading valve in order to monitor the force state of the valve ball.

    (d) The frequency demand of amplifier and data collector is very high,and needs to reach more than 20 times the loading frequency.

    To simulate the oil immersion condition of the unloading valve,a transverse flow channel is set up in the lower part of the experimental device.Gasolinelike medium T16 is selected as the test medium,which can simulate the properties of gasoline well and has good aging and corrosion resistance.In addition,the test medium has good stability and is not as flammable and explosive as gasoline,so it can ensure the safety of the experiment.

    The final experimental device composition is shown in Fig.3.The piezoelectric ceramics is connected with the flange plate and the upper base through the shell thread.The flange part is also fixed with the upper base through the four hollowed out strips.The connecting rod No.1,the force sensor,and the connecting rod No.2 are connected in turn to form a “sandwich” structure,and the tested ball valve is fixed through the groove on the lower base.The high frequency alternating load generated by the piezoelectric ceramics can be transmitted to the valve ball through the “sandwich” structure.In view of the heatdissipation problem of the driving power supply of the piezoelectric ceramics under the action of high frequency for a long time,it is necessary to design a heat dissipation channel in the power supply part,and the external ventilation and heat dissipation are carried out by an air compressor.The finally designed system is shown in Fig.4.

    Fig.3 Final experiment device composition.1:piezoelectric ceramics;2:flange plate;3:connecting rod No.1;4:pressure sensor;5:connecting rod No.2;6:tested unloading valve;7-9:bases

    Fig.4 Sketch of the experiment system (a);photo of the experiment system(b)

    2.3 Measurement of wear width

    To compare the experimental data with the later simulation results,it is necessary to measure the wear width of the valve ball and valve seat after the fretting test.The measurement tool is a Leica optical microscope,which can directly measure and read the measurement results with an accuracy of 0.1 μm.The valve ball and valve seat need to be measured differently due to their different geometries.As for the valve seat,it can be directly installed on the fixture vertically in the axial direction.The observing direction of the microscope is also along the axial direction of the valve seat.Since the wear band is located on the sealing slope of the valve seat,it is necessary to convert the measurement data obtained from that direction.Let the width of the wear band obtained be

    s

    ′,and the actual wear band width

    s

    can be obtained by

    s

    =

    s

    ′/sin

    θ

    .As for the valve ball,adding it to the fixture cannot 100% ensure that the wear band is facing the observing direction.Therefore,for the measurement of the wear width of the valve ball,it is necessary to take the average of multiple measurement results to eliminate human error.Fig.5 shows the schematic diagram of seat and ball wear width measurements.

    Fig.5 Schematic diagram of seat and ball wear width measurements

    3 Simulation methodology

    3.1 Theoretical analysis of the statically undetermined structure

    According to the structure and force state of the valve,the valve ball and the valve seat are initially in line contact.To facilitate the analysis of the contact force,the structure of the valve is divided into

    m

    equal divisions as shown in Fig.6.The force state of each micro-element is shown in Fig.7,where Δ

    F

    N

    ,and Δ

    f

    respectively represent the equivalent pressure,the normal pressure,and the friction of the microelement.The force state can be described by

    Fig.6 Micro-division of the structure(Lu et al.,2019)

    Fig.7 Force state of each micro-element(Lu et al.,2019)

    The normal pressure d

    N

    and friction d

    f

    cannot be directly solved by Eq.(1)alone since the force state is statically indeterminate.Therefore,more equations are needed as supplementary conditions.The equivalent friction coefficient

    μ

    between two contact surfaces under static friction can be described as

    where

    k

    represents a constant value of the material;Δ

    x

    represents the micro-slip distance under static friction.The relationship between the normal pressure and friction acting on the micro-element can be described as follows:

    Fig.8 shows the micro-slip model of the valve ball under the preload of the spring,where

    r

    represents the distance between the ball center and valve seat after the micro-slip happens;

    a

    represents the halfwidth of contact between the valve ball and valve seat.Two equations can be concluded from the model as follows:

    Fig.8 Micro-slip model of the micro-element

    The half-width of contact can be solved by Hertz contact theory(Johnson,1985).Fig.9 shows the theoretical model given by Hertz contact theory for the distribution of contact stress on the contact area of two cylinders.The two cylindrical surfaces are in close contact due to the normal pressure

    P

    .In this way,the contact area changes from line contact to surface contact.The shape of the contact stress distribution of the contact area is an ellipse.

    σ

    represents the maximum contact pressure.The half-length of the major axis of the ellipse is the half-width of contact

    a

    .The distribution of the contact stress can be described as follows:

    Fig.9 Hertz contact stress model

    where

    p

    represents the maximum contact stress;

    x

    represents the distance between a point on the contact surface and the initial contact point.The half-width of contact and the maximum contact stress can be described as follows:

    where

    P

    represents the line pressure on two cylinders(N/mm),which can be seen as the normal pressure on the micro-element;

    E

    represents the equivalent elastic modulus,which can be described as

    where

    E

    and

    E

    represent the elastic moduli of two materials,and

    ν

    and

    ν

    represent Poisson’s ratios of the two materials,respectively.

    G

    represents the relative curvature,which can be described as

    where

    G

    and

    G

    respectively represent the relative radius of curvature of two cylinders.If one of the cylinders is changed into a flat,it can be regarded as a cylinder with infinite radius of curvature.Eqs.(1)-(10)are taken into consideration and all the known parameters are substituted in the above equations and give these results:the equivalent friction coefficient

    μ

    =0.0887;the micro-slip distance Δ

    x

    =0.355 μm;the half-width of contact

    a

    =17.15 μm;the maximum contact stress

    p

    =1220 MPa.

    3.2 Archard’s wear model

    According to the Archard equation,the wear volume of the material is directly proportional to the slip distance and positive pressure,but inversely proportional to the hardness of the material.The Archard equation for sliding wear is normally expressed as(Archard,1953):

    where

    V

    is the wear volume,

    S

    is the sliding distance,

    H

    is the hardness of the material,and

    K

    is the wear coefficient.

    After its proposal in the 1950s,the Archard equation was immediately applied to the experimental study of fretting wear.At the same time,it was also applied to analytical and finite element simulation calculations.Fouvry (2001) proposed a simplified model based on the Archard equation using the Archard factor and the corresponding partial Archard factor density,in which the local Archard factor density describes the distribution of the Archard factor on the contact surface and is directly proportional to the wear depth.McColl et al.(2004) simulated and analyzed the conventional fretting wear device in a cylindrical plane based on the Archard wear model.They predicted the change of fretting contact surface profile with the number of cycles and studied in detail the influence of fretting wear on fretting fatigue.

    In the finite element simulation,it is necessary to control the node coordinates of the contact surface and a wear model suitable for fretting wear is needed.However,Archard’s wear model was proposed based on macro wear,so it needs to be modified to be suitable for finite element simulation.In the finite element method,the macro object is divided into finite elements,which are composed of nodes,and then the macro object is simulated and analyzed by the finite element method.Each node represents a local location.If the applicability of the two wear models is extended from the macro object to the local position,the local node coordinates can be controlled scientifically,so as to characterize the wear of materials and realize the quantitative simulation of fretting wear.

    In finite element simulation analysis,differentiating both sides of the Archard equation,we can obtain:

    To simulate the change of contact surface profile with wear times,it is necessary to obtain the wear depth d

    h

    at each position,and set d

    h

    (

    x

    ) as the wear depth at node position

    x

    .For a very small contact surface d

    A

    ,it is assumed that when the relative sliding distance of the contact surface is d

    S

    ,the wear depth of this part is d

    h

    .Dividing the left and right sides of the equation by d

    A

    at the same time:

    where d

    V

    /d

    A

    =d

    h

    is the local wear depth,and d

    P

    /d

    A

    is the local contact compressive stress

    p

    (

    x

    )at position

    x.

    Eq.(13)can be expressed as

    where d

    s

    (

    x

    ) is the relative sliding distance at position

    x

    ,and

    k

    is the wear coefficient,

    k

    =

    K

    /

    H

    .By transforming Eq.(14),the expression of wear depth at position

    x

    can be obtained:

    Eq.(15) is the expression of wear depth when a relative sliding distance d

    s

    (

    x

    ) occurs at position

    x

    .In the finite element simulation calculation,the action process of a cyclic load can be divided into

    n

    increments,and in any incremental step

    i

    ,the contact stress

    p

    (

    x

    ) at position

    x

    remains unchanged,and the relative slip is Δ

    s

    (

    x

    ),where Δ

    s

    (

    x

    )=

    s

    (

    x

    )-

    s

    (

    x

    ),and

    s

    (

    x

    )and

    s

    (

    x

    )are the corresponding sliding distances at the

    i

    th and (

    i

    -1)th incremental steps,respectively.The wear depth at position

    x

    in incremental step

    i

    can be expressed as

    During a complete cyclic load,the wear depth at position

    x

    is the sum of the wear depth of all incremental steps,and can be expressed as

    where

    n

    represents the total incremental steps.

    3.3 Numerical calculation method of cyclic loading

    3.3.1 Numerical iterative acceleration method

    If one step of a finite element simulation calculation only represents the action of one cycle load,the finite element analysis will include thousands of load steps,which will be very time-consuming.In fretting wear,the amount of wear produced by one fretting cycle is very small,and cannot cause significant changes in the contour of the contact surface,so it will not affect the contact stress and other parameters calculated by simulation.Therefore,it is not necessary to update the geometry after each cyclic load in the simulation.It can be assumed that after Δ

    N

    times of cyclic load,the contact surface profile changes significantly.At this time,the geometric morphology of the wear surface needs to be updated,so that the contact stress and relative slip distance after the change of the contact surface profile can be obtained,and provide more accurate data for wear calculation.Therefore,it can be represented by a load step of finite element simulation calculation Δ

    Q

    times of cyclic load,so as to accelerate the calculation and greatly reduce the analysis and calculation time without affecting the calculation accuracy.When the number of accelerations is Δ

    Q

    ,the wear depth at position

    x

    can be expressed as

    3.3.2 Mesh node updating method

    The difficulty of fretting wear simulation is that the contact surface profile changes with the wear process,and leads to a change of contact stress in the contact area,resulting in the real-time change of wear.In this paper,the UMESHMOTION subroutine suitable for finite element mesh generation is written in FORTRAN.The function of the subroutine is to obtain the contact stress and relative slip distance,calculate the wear depth according to the modified Archard equation,and update the nodes of the wear surface of the model after wear.Node updating is realized by adjusting the mesh element of the finite element,which itself will not produce stress and strain.Because of this characteristic,it can be used to simulate fretting wear.The movement of mesh nodes is likely to lead to a huge deformation of the element.In particular,the updated surface of the node may be distorted and unable to continue the calculation,so it also needs the function of automatic mesh adjustment to cooperate,that is,the Arbitrary Lagrangian-Eulerian(ALE) technology in ABAQUS,as shown in Fig.10.This function can make the transition between mesh nodes smooth,so that the mesh can always maintain a more appropriate shape and ensure the convergence of calculation.

    Fig.10 ALE adaptive meshing technology

    3.3.3 Numerical model

    According to Eq.(8),at the end of each step,the UMESHMOTION subroutine realizes the calculation of damage depth and the update of the grid.Therefore,the flow chart of fretting wear numerical simulation is shown in Fig.11,and the specific steps are as follows:

    Fig.11 Wear simulation procedure flowchart

    1.Determining the initial parameters of the finite element model.

    2.Establishing the geometric model in ABAQUS software and generating the INP file.This file contains all the information of the finite element model and can be directly called upon by the solver for analysis.

    3.Calling the modified Archard equation for calculation.

    4.Solving the INP file in ABAQUS to extract the contact stress and relative sliding distance at each node in the contact area.

    5.At the end of each incremental step,the UMESHMOTION subroutine is called upon to calculate the wear depth Δ

    h

    and wear volume Δ

    V

    of each node according to the modified Archard wear model.Then,the node is updated according to the wear depth,and the coordinate at node

    i

    is

    h

    =

    h

    h

    .

    6.After the node is updated,a new surface profile and a new INP file are generated.

    7.Counting the number of fretting cycles and judging whether the total number of cycles is reached.If it is not reached,the judgment result is YES,repeating steps 4-7.The fretting wear continues and the wear amount continues to accumulate;if the total number of cycles is judged to have been reached,the judgment result is NO,and the simulation calculation of fretting wear is completed.

    3.4 Fretting wear model of the unloading valve

    3.4.1 Two-dimensional simplification of computational domain

    The aim of this study is the unloading valve in a high-pressure pump,which has circumferential symmetry in structure and the stress conditions on any section are the same.Considering that the original ball structure and more complex valve seat structure are not conducive to grid division and refinement,a 2D spherical/plane contact structure is adopted for simulation.The established simulation structure is shown in Fig.12.The valve seat and valve ball are equivalent to a 2D plane structure.This simplification has the following advantages:(1) Compared with the original structure,it is conducive to refining the mesh of the contact area;(2)Compared with 3D model and slice model,the number of meshes is less,which improves the simulation speed;(3)The boundary conditions suitable for fretting simulation can be added to the model;(4) The 2D grid structure is more conducive to the updating of the grid and can improve the accuracy of calculation.

    Fig.12 Simplification simulation structure of the unloading valve

    3.4.2 Meshing of computational domain

    In the simulation structure of this paper,the valve ball and valve seat adopt CPE4 grid of neutral axis algorithm,in which the valve ball contains 1598 elements and 1670 nodes,and the valve seat contains 1520 elements and 1590 nodes.As shown in Fig.13a,in order to prevent the fluctuation of contact pressure,the grid elements in the contact area between the valve ball and the valve seat correspond one to one.To accurately simulate the parameter changes in this area(such as relative slip amplitude,contact stress,and stress distribution on the sub surface),the grid in this area is finely divided (Fig.13b),with a size of about 10 μm.Considering that the non-contact part does not play a key role,this saves calculation time;the grid is sparsely divided,and the size is about 100 μm.

    Fig.13 The 2D simulation model of the unloading valve(a);mesh refinement of contact area(b)

    Checking the meshing quality,it is found that large or small element angles (minimum angle 45°,maximum angle 135°)appear on a few elements above the valve ball contact area (Fig.14).However,these elements are not the key parts of simulation and the coverage area is also very small,so their impact can be ignored.The average minimum element angle is 75.21°,and the average maximum element angle is 104.78°,then the average aspect ratio is 1.92.There are no elements with too large an aspect ratio and no intersecting elements and zero volume.Generally speaking,the quality of mesh generation is satisfactory and can meet computational requirements.

    Fig.14 Deviation of the grid angle

    3.4.3 Setting of boundary condition

    The calculation model of fretting wear of a valve ball is established in ABAQUS,and the material parameters of the model are shown in Table 2.Fixed constraints at the bottom of the valve seat and point edge coupling constraints on the upper surface of the valve ball are added to ensure that the valve ball is always level during loading.A periodic normal load

    p

    is applied on the upper surface of the valve ball,and a periodic tangential displacement load

    s

    is applied on the centre line,so that the two specimens are in close contact and slide relatively,resulting in fretting wear(Fig.15).

    Table 2 Material properties

    Fig.15 Boundary conditions of the model

    Fretting wear occurs in the contact area between the two surfaces of the valve ball and valve seat.The lower surface of the valve ball is set as the main surface,and the upper surface of the valve seat is set as the slave surface.In this way,ABAQUS can quickly and accurately calculate the contact parameters (relative slip distance,contact stress,and shear friction)of the contact area.The tangential behavior is set to the isotropic Coulomb friction law calculated by a penalty function,and the friction coefficient adopts the previously calculated value 0.0887.The normal behavior is set to the hard contact constrained by the standard penalty function.The contact algorithm adopts the finite slip algorithm,and the tangential constraint algorithm is set as in the Lagrange multiplier method.In this way,the relative sliding distance and contact stress obtained will be more accurate.

    The two extreme values of valve ball periodic force(5 and 35 MPa)are taken as the boundary conditions of positive pressure and micro slip values,and the positive pressure and micro slip values are added to the simulation load conditions in the form of a periodic alternating cycle.The number of accelerations in the acceleration calculation method Δ

    N

    is set as 1000.Limited by the simulation time,all simulation cases are set as 450 steps.The data of positive pressure and micro-slip distance are shown in Table 3.The two loads are in an alternating state to simulate the respective radial and tangential fretting forms of the valve ball,and thus simulate its fretting wear.

    Table 3 Loading conditions at different contact angles

    3.5 Verification of wear model

    The experimental and simulation results of the wear width of the valve ball under different half contact angles with the hardness of the ball and of the seat as HV700 are shown in Fig.16.

    Each of the experiments was conducted for 25 h at 500 Hz frequency,so 25×60×60×500=4.5×10cycles were performed for each experimental case.The number of accelerations in the acceleration calculation method Δ

    N

    was set as 10000 and all simulation cases were set as 4500 steps,so 10000×4500=4.5×10cycles were performed for each simulation case.The simulation and experiment have the same cycle times.Therefore,if the simulation results are in good agreement with the experimental results,the accuracy and reasonableness of the simulation have been demonstrated.

    It can be seen from Fig.16 that both the experiment and simulation show a decreasing trend of wearwidth as the half contact angle increases.Furthermore,errors between the simulation and experimental results under different half contact angles are all less than ±10°.Considering those experimental and measurement errors,it can be assumed that the simulation calculates the fretting wear of the unloading valve well.

    Fig.16 Comparison between experiment results and simulation results of wear width under different half contact angles.Valve ball material:9Cr18Mo;valve seat material:9Cr18MoV;valve ball hardness:HV700;valve seat hardness:HV700

    4 Fretting wear characteristics of the unloading valve

    4.1 Wear width and depth of different half contact angles

    Based on the verified wear model,it is not difficult to find that the wear width of the valve ball decreases as the half contact angle increases.Therefore,it could be speculated that perhaps both the wear width and the wear depth of the valve ball will decrease with the increase in half contact angle.To test this conjecture,simulation cases of five different half contact angles were conducted under the condition of hardness of HV700 for both ball and seat.Because much time is taken if each simulation case is conducted with the same number of cycles as the experimental case the number of accelerations,Δ

    Q

    ,was set as 1000 and the number of steps was set as 4500.The number of cycles of the simulation is one tenth of that of the experiment.The simulation results are shown in Figs.17 and 18.

    It can be seen from Fig.17 that the wear profile of the valve ball is roughly U-shaped.The wear widths on both sides of the highest point of the U-shape are different because the stress states on both sides of the contact point between the valve ball and the valve seat are different,which is consistent with the wear profile of the valve ball in the experimental result.

    Fig.17 Wear depth profile variation with different half contact angles (valve ball material:9Cr18Mo;valve seat material:9Cr18MoV;valve ball hardness:HV700;valve seat hardness:HV700) (a);valve ball wear profile (2000-h experiment)(b)

    Fig.18 shows the variation of the wear width and depth with different half contact angles.As predicted before,both the wear width and depth decrease with the increase in the half contact angle.This result may be explained by the fact that the wear width is directly proportional to the contact radius between the valve ball and the valve seat,and so the larger the contact area,the larger the wear area.According to Hertz contact theory,the contact radius is directly proportional to the positive pressure and inversely proportional to the contact length.In the process of gradually increasing the contact half angle,the positive pressure between the valve ball and the valve seat increases and the contact length also increases.However,the ratio of the positive pressure and the contact length decreases,so the contact radius decreases with the increase of the contact half angle,and the wear width,which is proportional to the contact radius,also shows the same trend.Therefore,it can be concluded that the wear width and depth of the valve ball will decrease with the increase of the half contact angle of valve ball and valve seat.

    Fig.18 Wear width and depth variation with different half contact angles

    4.2 Wear width and depth of different valve ball radii

    It can be seen from Table 4 that the increase of valve ball radius will lead to the decrease of positive pressure and micro-slip distance.This leads us to consider how the change of valve ball radius will affect the positive pressure,micro slip distance,and fretting wear of the valve ball.Therefore,according to the same method as in Section 3.1,the radius of the valve ball is changed and the statically indeterminate structure of the ball valve is solved correspondingly,with the half contact angle of valve ball fixed as 27.5° and the hardness of valve ball and valve seat set as HV700.The results are shown in Table 4.The increase of valve ball radius also leads to the increase of positive pressure and the decrease of micro-slip distance.Hence,it cannot directly be inferred from the trend of positive pressure and micro-slip distance that that is how the change of valve ball radius will influence the fretting wear of the valve ball.So simulations are needed to confirm the results.

    Table 4 Loading conditions at different contact angles

    In the same way as shown in Section 4.1,the number of accelerations Δ

    Q

    was set as 1000 and the number of steps was set as 4500.The simulation results are shown in Fig.19.

    Fig.19 Wear width and depth variation with different valve ball radii

    As shown in Fig.19,both the wear width and depth decrease as the valve ball radius increases.Although the positive pressure and micro-slip distance show different variation trends compared to the half contact angle situation,the wear width and depth show the same trend.This may be explained by the fact that the product of positive pressure and microslip distance keeps increasing,although the latter is decreasing.Therefore,it can be concluded that the wear width and depth of the valve ball will decrease with the increase of the valve ball radius.

    4.3 Parameter optimization selection

    After determining the influence characteristics of radius and half contact angle,it is essential to select the optimal parameters in combination with the actual working conditions of the unloading valve.

    Currently the contact half angle of the spherical unloading valve in use is 27.5°.According to the influence characteristics of the contact half angle summarized above,this is not the optimal angle to reduce wear.The diameter of the oil pipe of the valve seat is 1.25 mm and the radius of the valve ball is 0.794 mm.From the engineering perspective,priority should be given to meeting its oil passage area.Therefore,according to its structural parameters,as shown in Fig.20,the maximum possible half contact angle is

    θ

    ,which can be calculated by

    Fig.20 Boundary half contact angle analysis

    From the analysis of the fit space between the valve seat and the high-pressure pump body,the minimum half contact angle is

    θ

    when the maximum margin

    v

    is about 5 mm.

    θ

    can be calculated by

    So the range of the half contact angle is 1.94° to 38°.As for the valve ball radius,the value range of 0.704 to 1.963 mm can be calculated in a similar way as the half contact angle.According to the discussion in Sections 4.1 and 4.2,the wear width and depth decrease monotonically with the increase of half contact angle and valve ball radius.Therefore,the maximum values of contact half angle of 38° and ball radius of 1.963 mm can be considered as a better choice for optimizing the unloading valve parameters.

    5 Conclusions

    Increasingly high pressures in gasoline direct injection technology lead to structural damage in spherical unloading valves.In previous studies,the authors have analyzed the stress characteristics of the unloading valve for addressing this problem and have pointed out that it is caused by fretting wear.Taking that as the premise,this essay describes the design and realization of experiments on fretting wear.Furthermore,theoretical analysis of the statically undetermined structure was made to acquire the boundary conditions of simulations.Based on the modified Archard wear equation,a quantitative calculation model for the unloading valve was established.The following conclusions can be obtained:

    1.For 4.5×10cycles fretting wear of the unloading valve,the simulation results are in good agreement with the experiment results.Therefore,the accuracy and rationality of the model can be proved.

    2.For 4.5×10cycles fretting wear of the unloading valve under the same hardness of valve ball and valve seat,the simulation results show that the profile of the valve ball is asymmetrically U-shaped.Besides,both the wear width and the wear depth of the valve ball decrease with the increase of the half contact angle.

    3.For 4.5×10cycles fretting wear of the unloading valve under identical hardness of valve ball and valve seat,the simulation results show that both the wear width and the wear depth of the valve ball decrease with the increase of the valve ball radius.

    Considering the actual application,the half contact angle and valve ball radius have a certain value range.After determining that the value range of the half contact angle is 1.94° to 38° and the value range of valve ball radius is 0.704 to 1.963 mm,according to the law that the wear depth and width decrease with the increase of half contact angle and valve ball radius,suggestions for the optimal selection of unloading valve parameters are made,that is,the valve ball radius should be 1.963 mm and the contact half angle should be 38°.

    Acknowledgments

    This work is supported by the National Key R&D Program of China (No.2019YFB2001502),the National Natural Science Foundation of China (No.52075387),the Open Foundation of the State Key Laboratory of Fluid Power &Mechatronic Systems (No.GZKF-201907),the Shanghai Municipal Science and Technology Major Project (No.2021SHZDZX0100),and the Shanghai Municipal Commission of Science and Technology Project (No.19511132101),China.

    Author contributions

    Liang LU and Qi-long XUE designed the research.Man-yi ZHANG,Liang-liang LIU,and Zhong-yu WU provided the experiment equipment and measured data.Yin-peng XU and Liang LU processed the corresponding data.Yin-peng XU wrote the first draft of the manuscript.Meng-ru LI helped to organize the manuscript.Yin-peng XU revised and edited the final version.

    Conflict of interest

    Liang LU,Yin-peng XU,Meng-ru LI,Qi-long XUE,Man-yi ZHANG,Liang-liang LIU,and Zhong-yu WU declare that they have no conflict of interest.

    a级片在线免费高清观看视频| 成人18禁高潮啪啪吃奶动态图| 日日撸夜夜添| 欧美人与性动交α欧美精品济南到 | 亚洲视频免费观看视频| 午夜免费男女啪啪视频观看| 久久婷婷青草| 国产亚洲精品第一综合不卡| 免费黄网站久久成人精品| 在线观看人妻少妇| videosex国产| 亚洲人成电影观看| 欧美另类一区| 国产亚洲欧美精品永久| 国产亚洲欧美精品永久| 欧美中文综合在线视频| 亚洲欧洲日产国产| 91久久精品国产一区二区三区| 久久ye,这里只有精品| 王馨瑶露胸无遮挡在线观看| 自线自在国产av| 亚洲国产日韩一区二区| 两性夫妻黄色片| 国产人伦9x9x在线观看 | 这个男人来自地球电影免费观看 | 下体分泌物呈黄色| 巨乳人妻的诱惑在线观看| 丝袜美腿诱惑在线| 十分钟在线观看高清视频www| 久久国产亚洲av麻豆专区| 国产成人欧美| 高清视频免费观看一区二区| 一区二区三区精品91| 丰满乱子伦码专区| 黄色视频在线播放观看不卡| 精品视频人人做人人爽| 纵有疾风起免费观看全集完整版| 777米奇影视久久| 一本久久精品| 精品亚洲成国产av| 日本免费在线观看一区| 日日撸夜夜添| 国产在视频线精品| tube8黄色片| 精品国产乱码久久久久久男人| 性少妇av在线| 亚洲国产欧美网| 久久精品国产a三级三级三级| 久久久国产精品麻豆| 日韩伦理黄色片| 最近中文字幕高清免费大全6| 激情视频va一区二区三区| 热99久久久久精品小说推荐| 18在线观看网站| 精品一区在线观看国产| 国产日韩欧美视频二区| 国产伦理片在线播放av一区| 大香蕉久久成人网| 在线观看美女被高潮喷水网站| 日韩中文字幕欧美一区二区 | 久久国产精品男人的天堂亚洲| 中文字幕另类日韩欧美亚洲嫩草| 在线观看三级黄色| 中文乱码字字幕精品一区二区三区| 亚洲第一av免费看| 交换朋友夫妻互换小说| 日本猛色少妇xxxxx猛交久久| 国产老妇伦熟女老妇高清| 美女中出高潮动态图| 最近最新中文字幕大全免费视频 | 国产熟女午夜一区二区三区| 久久久精品国产亚洲av高清涩受| 亚洲av福利一区| 在线天堂中文资源库| 母亲3免费完整高清在线观看 | 欧美少妇被猛烈插入视频| 亚洲精品国产色婷婷电影| 国产精品女同一区二区软件| 午夜日韩欧美国产| 在线观看人妻少妇| 成人毛片60女人毛片免费| 又粗又硬又长又爽又黄的视频| 色婷婷久久久亚洲欧美| 日韩欧美一区视频在线观看| 婷婷色麻豆天堂久久| 国产高清国产精品国产三级| 丰满饥渴人妻一区二区三| 伦精品一区二区三区| 波多野结衣一区麻豆| 国精品久久久久久国模美| 三上悠亚av全集在线观看| 国产精品一区二区在线观看99| 精品国产一区二区久久| 欧美人与性动交α欧美软件| 男人添女人高潮全过程视频| 天天躁夜夜躁狠狠久久av| 91久久精品国产一区二区三区| 国产亚洲一区二区精品| 免费大片黄手机在线观看| 亚洲美女视频黄频| 高清欧美精品videossex| 国产精品麻豆人妻色哟哟久久| 国产成人精品久久二区二区91 | 91在线精品国自产拍蜜月| 欧美日韩亚洲高清精品| 最近中文字幕2019免费版| 亚洲成国产人片在线观看| 午夜老司机福利剧场| tube8黄色片| 交换朋友夫妻互换小说| 亚洲国产精品一区三区| 蜜桃国产av成人99| 人人妻人人澡人人爽人人夜夜| 十分钟在线观看高清视频www| 老司机亚洲免费影院| 午夜福利一区二区在线看| 汤姆久久久久久久影院中文字幕| 看免费av毛片| 高清在线视频一区二区三区| 国产福利在线免费观看视频| 成人18禁高潮啪啪吃奶动态图| 欧美成人精品欧美一级黄| 免费高清在线观看日韩| 欧美精品国产亚洲| 建设人人有责人人尽责人人享有的| 国产探花极品一区二区| 久久久久久伊人网av| 日本欧美国产在线视频| 午夜免费男女啪啪视频观看| 免费在线观看黄色视频的| 精品国产国语对白av| 亚洲欧美色中文字幕在线| 电影成人av| 国产综合精华液| 国产无遮挡羞羞视频在线观看| tube8黄色片| 九草在线视频观看| 男女下面插进去视频免费观看| 亚洲av国产av综合av卡| av福利片在线| 日韩三级伦理在线观看| 亚洲精品久久午夜乱码| av电影中文网址| 日韩精品免费视频一区二区三区| 中文欧美无线码| 毛片一级片免费看久久久久| 在线观看一区二区三区激情| 国产av码专区亚洲av| 日本av免费视频播放| 女人精品久久久久毛片| 亚洲精品久久午夜乱码| 人妻人人澡人人爽人人| 妹子高潮喷水视频| 人妻 亚洲 视频| 国产av码专区亚洲av| 国产精品三级大全| 欧美日韩亚洲高清精品| 狠狠婷婷综合久久久久久88av| 国产激情久久老熟女| 天堂俺去俺来也www色官网| 成人免费观看视频高清| www.精华液| 妹子高潮喷水视频| 免费女性裸体啪啪无遮挡网站| 黑人猛操日本美女一级片| 另类精品久久| 波野结衣二区三区在线| 免费久久久久久久精品成人欧美视频| 999精品在线视频| 99国产综合亚洲精品| 国产男女内射视频| 香蕉精品网在线| 丝袜美腿诱惑在线| 丰满迷人的少妇在线观看| 999久久久国产精品视频| 中文字幕人妻丝袜制服| 久久久久久久久免费视频了| 久久久a久久爽久久v久久| 日韩一本色道免费dvd| 久久久久久久久免费视频了| 欧美av亚洲av综合av国产av | 国产人伦9x9x在线观看 | 欧美日韩亚洲国产一区二区在线观看 | 成人国产麻豆网| 国产精品久久久久成人av| 日韩在线高清观看一区二区三区| 桃花免费在线播放| 一区二区三区激情视频| 中文字幕另类日韩欧美亚洲嫩草| 久久人人爽av亚洲精品天堂| 精品亚洲成国产av| 香蕉国产在线看| 一级毛片我不卡| kizo精华| 女人精品久久久久毛片| 另类精品久久| 永久网站在线| 黄色怎么调成土黄色| 国产野战对白在线观看| 一级,二级,三级黄色视频| 久久 成人 亚洲| 亚洲三级黄色毛片| 中文字幕人妻熟女乱码| 国产精品久久久久久精品电影小说| 中文字幕人妻丝袜一区二区 | 亚洲四区av| 国产成人aa在线观看| 捣出白浆h1v1| 久久久久国产一级毛片高清牌| 国产亚洲一区二区精品| 国产成人精品久久久久久| 黄片播放在线免费| 999精品在线视频| 精品少妇黑人巨大在线播放| 婷婷色麻豆天堂久久| 亚洲成国产人片在线观看| 中文字幕精品免费在线观看视频| 高清在线视频一区二区三区| 免费在线观看完整版高清| 这个男人来自地球电影免费观看 | 一二三四在线观看免费中文在| 亚洲国产精品999| 91午夜精品亚洲一区二区三区| 亚洲av男天堂| 亚洲 欧美一区二区三区| 韩国高清视频一区二区三区| 成人亚洲欧美一区二区av| 一级毛片电影观看| 精品人妻一区二区三区麻豆| 国产精品成人在线| 超碰成人久久| 国产成人精品福利久久| av天堂久久9| 男女午夜视频在线观看| 国产精品无大码| 热99国产精品久久久久久7| 在线精品无人区一区二区三| 国产精品欧美亚洲77777| 伊人久久大香线蕉亚洲五| 亚洲精品日韩在线中文字幕| 美国免费a级毛片| 欧美人与善性xxx| 国产精品人妻久久久影院| 老汉色∧v一级毛片| 午夜福利网站1000一区二区三区| 久久午夜福利片| 国产精品不卡视频一区二区| av又黄又爽大尺度在线免费看| 午夜精品国产一区二区电影| 只有这里有精品99| 中文字幕精品免费在线观看视频| 亚洲中文av在线| www.熟女人妻精品国产| tube8黄色片| 久久精品国产自在天天线| 97精品久久久久久久久久精品| 老司机影院毛片| 久久久久久久久免费视频了| 久热这里只有精品99| 国产成人免费观看mmmm| 在线亚洲精品国产二区图片欧美| 视频在线观看一区二区三区| 日本欧美国产在线视频| 久久国产精品大桥未久av| 国产一区二区 视频在线| 亚洲三级黄色毛片| 性少妇av在线| 女的被弄到高潮叫床怎么办| 久久99精品国语久久久| 国产 精品1| 叶爱在线成人免费视频播放| 亚洲,一卡二卡三卡| 久久久精品免费免费高清| 色婷婷av一区二区三区视频| 女的被弄到高潮叫床怎么办| 国产黄频视频在线观看| 国产精品免费视频内射| 亚洲欧美一区二区三区黑人 | 男的添女的下面高潮视频| 一区二区三区乱码不卡18| 91精品三级在线观看| 97在线视频观看| av在线app专区| 亚洲色图 男人天堂 中文字幕| 国产成人aa在线观看| 在线天堂最新版资源| 啦啦啦视频在线资源免费观看| 久久毛片免费看一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 99国产综合亚洲精品| 欧美人与性动交α欧美精品济南到 | 性少妇av在线| 午夜日韩欧美国产| 各种免费的搞黄视频| 不卡av一区二区三区| 青春草视频在线免费观看| 久久久亚洲精品成人影院| 91在线精品国自产拍蜜月| 十八禁高潮呻吟视频| 18在线观看网站| 女性被躁到高潮视频| 免费高清在线观看日韩| 热re99久久国产66热| 国产探花极品一区二区| 搡老乐熟女国产| 秋霞伦理黄片| 日韩三级伦理在线观看| 久久久国产欧美日韩av| 女人精品久久久久毛片| 97在线人人人人妻| 精品一区二区免费观看| 韩国精品一区二区三区| 观看av在线不卡| 国产乱来视频区| 亚洲av欧美aⅴ国产| 国产在视频线精品| 国产熟女欧美一区二区| av国产久精品久网站免费入址| 不卡视频在线观看欧美| 亚洲精品久久午夜乱码| 99国产精品免费福利视频| 制服人妻中文乱码| 欧美激情 高清一区二区三区| 亚洲欧洲日产国产| 一区二区日韩欧美中文字幕| www.精华液| 欧美在线黄色| xxx大片免费视频| 亚洲少妇的诱惑av| 菩萨蛮人人尽说江南好唐韦庄| 纯流量卡能插随身wifi吗| 91国产中文字幕| 日韩一区二区三区影片| av在线观看视频网站免费| 黄色视频在线播放观看不卡| 精品人妻偷拍中文字幕| 久久精品夜色国产| 老鸭窝网址在线观看| 色婷婷久久久亚洲欧美| 永久免费av网站大全| 中文字幕精品免费在线观看视频| 免费观看av网站的网址| 最近手机中文字幕大全| 免费观看在线日韩| 又黄又粗又硬又大视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲成人一二三区av| 免费高清在线观看视频在线观看| 国产精品一二三区在线看| freevideosex欧美| 韩国高清视频一区二区三区| 如日韩欧美国产精品一区二区三区| 天堂俺去俺来也www色官网| 这个男人来自地球电影免费观看 | 欧美在线黄色| 久久久国产一区二区| 午夜福利视频在线观看免费| 可以免费在线观看a视频的电影网站 | 午夜福利视频精品| 亚洲国产精品999| 最近最新中文字幕免费大全7| a 毛片基地| 日本黄色日本黄色录像| 欧美精品国产亚洲| 纯流量卡能插随身wifi吗| 国产精品国产三级专区第一集| 人妻人人澡人人爽人人| 成人午夜精彩视频在线观看| 91国产中文字幕| 久久午夜福利片| 日韩欧美一区视频在线观看| 欧美日韩精品成人综合77777| 曰老女人黄片| 精品酒店卫生间| 日本爱情动作片www.在线观看| 精品人妻熟女毛片av久久网站| 国产不卡av网站在线观看| 国产片内射在线| 国产精品国产av在线观看| 男女国产视频网站| 一区二区日韩欧美中文字幕| 赤兔流量卡办理| 久久 成人 亚洲| 一区二区日韩欧美中文字幕| 国产一区二区在线观看av| 国产精品免费视频内射| 亚洲av日韩在线播放| 久久精品国产综合久久久| www.av在线官网国产| 国产视频首页在线观看| 午夜精品国产一区二区电影| 欧美xxⅹ黑人| 美女xxoo啪啪120秒动态图| 丝瓜视频免费看黄片| 电影成人av| 午夜福利影视在线免费观看| 国产精品不卡视频一区二区| 日韩精品有码人妻一区| 欧美日韩精品网址| 久久久久国产网址| 亚洲精品国产一区二区精华液| 丝瓜视频免费看黄片| 国产精品成人在线| 99久国产av精品国产电影| 欧美人与善性xxx| 日韩制服骚丝袜av| 少妇 在线观看| 你懂的网址亚洲精品在线观看| 精品少妇久久久久久888优播| 国产xxxxx性猛交| 午夜av观看不卡| 婷婷色综合www| 亚洲第一青青草原| 26uuu在线亚洲综合色| 亚洲国产精品成人久久小说| 亚洲成色77777| 9色porny在线观看| 新久久久久国产一级毛片| 免费不卡的大黄色大毛片视频在线观看| 日韩视频在线欧美| 伦理电影免费视频| 久久久久国产网址| 日韩免费高清中文字幕av| 久久精品国产亚洲av高清一级| 欧美精品亚洲一区二区| 制服诱惑二区| 韩国精品一区二区三区| 91麻豆av在线| 国产欧美日韩一区二区三| 人人澡人人妻人| 五月开心婷婷网| 香蕉久久夜色| av超薄肉色丝袜交足视频| 丰满饥渴人妻一区二区三| av福利片在线| 久久国产乱子伦精品免费另类| 校园春色视频在线观看| 欧美黄色片欧美黄色片| 色综合婷婷激情| 免费女性裸体啪啪无遮挡网站| 亚洲色图av天堂| 超色免费av| 亚洲精品在线美女| 亚洲精华国产精华精| 一本综合久久免费| 女性生殖器流出的白浆| 搡老乐熟女国产| 久久精品国产亚洲av香蕉五月| 老司机午夜福利在线观看视频| 日韩精品青青久久久久久| 后天国语完整版免费观看| 女同久久另类99精品国产91| 亚洲欧美激情在线| 黄片大片在线免费观看| 国产欧美日韩一区二区三区在线| 夜夜夜夜夜久久久久| 在线观看一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 日韩视频一区二区在线观看| 欧美日韩瑟瑟在线播放| 久久性视频一级片| 日韩精品青青久久久久久| 亚洲成av片中文字幕在线观看| 人人澡人人妻人| 90打野战视频偷拍视频| 曰老女人黄片| 国产欧美日韩一区二区精品| 国产av在哪里看| 亚洲男人的天堂狠狠| 人人妻人人爽人人添夜夜欢视频| 我的亚洲天堂| 精品一品国产午夜福利视频| 夜夜夜夜夜久久久久| 久久人人97超碰香蕉20202| 欧美日韩视频精品一区| 亚洲欧美激情在线| 久久天躁狠狠躁夜夜2o2o| 不卡av一区二区三区| 黄色视频,在线免费观看| 黄色怎么调成土黄色| 国产91精品成人一区二区三区| 免费av中文字幕在线| 精品久久久久久,| 操美女的视频在线观看| 国产精品影院久久| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲人成网站在线播放欧美日韩| 中国美女看黄片| 纯流量卡能插随身wifi吗| 可以免费在线观看a视频的电影网站| 一进一出好大好爽视频| 久久 成人 亚洲| 高清黄色对白视频在线免费看| 久久精品91蜜桃| av片东京热男人的天堂| 少妇被粗大的猛进出69影院| 成人精品一区二区免费| 国产成人精品久久二区二区91| 999久久久精品免费观看国产| x7x7x7水蜜桃| 日韩欧美一区二区三区在线观看| 国产成人免费无遮挡视频| 满18在线观看网站| 999精品在线视频| 久久久国产成人精品二区 | 国产欧美日韩综合在线一区二区| 美女午夜性视频免费| 最近最新中文字幕大全电影3 | 午夜影院日韩av| 麻豆久久精品国产亚洲av | 亚洲视频免费观看视频| 色播在线永久视频| 18美女黄网站色大片免费观看| 天天躁夜夜躁狠狠躁躁| 亚洲avbb在线观看| 久久香蕉国产精品| 午夜激情av网站| 在线观看免费午夜福利视频| 在线观看免费视频网站a站| 在线播放国产精品三级| 1024香蕉在线观看| 91老司机精品| 亚洲熟妇熟女久久| 91国产中文字幕| 老熟妇仑乱视频hdxx| 无人区码免费观看不卡| 国产深夜福利视频在线观看| 视频区图区小说| 精品人妻1区二区| 欧美日韩一级在线毛片| 亚洲精品一区av在线观看| 久久久久国产一级毛片高清牌| 大型黄色视频在线免费观看| 夜夜夜夜夜久久久久| 国产亚洲精品第一综合不卡| 国产亚洲欧美精品永久| 日韩欧美一区二区三区在线观看| av天堂久久9| 人人妻人人爽人人添夜夜欢视频| 黄色怎么调成土黄色| 久久精品人人爽人人爽视色| 亚洲国产欧美网| 国产精品av久久久久免费| 久久精品亚洲精品国产色婷小说| 欧美成人性av电影在线观看| 黄色 视频免费看| 欧美日韩亚洲国产一区二区在线观看| 久久香蕉国产精品| 国产精品1区2区在线观看.| 精品国产一区二区久久| 国产一卡二卡三卡精品| 欧美乱色亚洲激情| 波多野结衣一区麻豆| 亚洲国产看品久久| 麻豆一二三区av精品| 久久亚洲真实| 精品人妻1区二区| 精品第一国产精品| 51午夜福利影视在线观看| 免费在线观看影片大全网站| 亚洲国产中文字幕在线视频| 国产91精品成人一区二区三区| 亚洲情色 制服丝袜| 天堂√8在线中文| 一区二区三区国产精品乱码| 99国产综合亚洲精品| 国产精品乱码一区二三区的特点 | av有码第一页| 欧美乱色亚洲激情| 99精国产麻豆久久婷婷| 一本综合久久免费| 亚洲国产精品一区二区三区在线| 亚洲国产精品合色在线| cao死你这个sao货| 黑人巨大精品欧美一区二区mp4| 亚洲一区二区三区欧美精品| 国产亚洲精品久久久久久毛片| 欧美成狂野欧美在线观看| 成人av一区二区三区在线看| 欧美中文日本在线观看视频| 啦啦啦 在线观看视频| 九色亚洲精品在线播放| 两人在一起打扑克的视频| 男女床上黄色一级片免费看| 一边摸一边抽搐一进一出视频| 视频区图区小说| a级片在线免费高清观看视频| 欧美 亚洲 国产 日韩一| 欧美日韩一级在线毛片| 久久香蕉国产精品| 看免费av毛片| 极品教师在线免费播放| 一边摸一边抽搐一进一出视频| 免费看十八禁软件| 亚洲国产欧美日韩在线播放| 久久久精品国产亚洲av高清涩受| 亚洲自拍偷在线| 精品久久久久久成人av| 久久中文看片网| 可以在线观看毛片的网站| 叶爱在线成人免费视频播放| 女同久久另类99精品国产91| 精品久久蜜臀av无| 国产一区二区激情短视频| 在线观看66精品国产| 狂野欧美激情性xxxx| 成人免费观看视频高清| 久久精品91蜜桃| 亚洲男人的天堂狠狠| 叶爱在线成人免费视频播放| 成人av一区二区三区在线看| 亚洲精品国产色婷婷电影| 亚洲精品中文字幕一二三四区| 亚洲av成人不卡在线观看播放网| 成人18禁高潮啪啪吃奶动态图|