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

    Erosion degradation characteristics of a linear electro-hydrostatic actuator under a high-frequency turbulent flow field

    2018-05-17 10:05:24YunLIShopingWANGMiletTOMOVICChoZHANG
    CHINESE JOURNAL OF AERONAUTICS 2018年5期

    Yun LI,Shoping WANG,Milet M.TOMOVIC,Cho ZHANG,*

    aSchool of Automation Science and Electrical Engineering,Beihang University,Beijing 100083,China

    bCollege of Engineering and Technology,Old Dominion University,Norfolk,VA 23529,USA

    1.Introduction

    Integrated Electro-Hydrostatic Actuators(EHAs)have seen increased application in More-Electric Aircraft(MEA),due to their numerous advantages including high reliability,long lifetime,and high efficiency.1–4Traditional power-integrated Rotary Electro-Hydrostatic Actuators(REHAs)are facing many problems,such as severe heating,big inertia,low frequency response,and difficulty in redundancy con figuration.In order to solve these deficiencies,many researchers have proposed direct-drive EHAs.5–9Li et al.have presented a novel LEHA for the pump control system based on a collaborative rectification structure for linear pumps where the fundamental subsystem consists of two Direct Drive Pump Cells(DDPCs).10Their prototype has verified feasibility of the proposed actuator and flexibility of the dual-control strategy.Although the system has a number of advantages,its reliability characteristics,such as degradation and life prediction,need to be addressed and further improved.Applying the novel operating principle,the LEHA has improved these performance issues,but it has also introduced a problem due to the fact that it cannot integrate filters.The novel LEHA is also an integrated closed system which can increase the power-to-weight ratio.A linear resonance motor drives piston cylinder suction and discharges oil.The spool of a rectification slide valve is rigidly connected with the driving cylinder rod,resulting in the active rectification mechanism to be in a high-frequency condition.Due to non- filter design,contaminated particles will continuously cycle in an LEHA during its operation,which will precipitate further erosion and produce more contamination particles.Consequently,compared to a conventional hydraulic system with filters,the rectification mechanism in an LEHA can fail more rapidly,which can lead to a dominant progressive degradation process.

    The rectification mechanism is commonly present in hydraulic systems,and when coupled with oil contamination,it can result in system failure.Some research results have indicated that contamination particles in oil wash out and wear the edges of valve components.11The resulting leakage increase caused by wear accounts for approximately 60%of cases of component failure.According to the research by Zhang et al.,the wear resulting from particulate contaminants causes an increase in the internal leakage,the output current hysteresis and null leakage,the input current threshold,and the pressure gain,as well as a decrease in the gain linearity of the electro hydraulic servo valve.12The Physics-of-Failure(PoF)models of particle erosion wear introduced by Fang et al.show that the erosion wear has significant impacts on the electro hydraulic servo valve’s service life and reliability.13Furthermore,the erosion of an LEHA is even more serious due to the particular type of motion and active rectification under relatively severe conditions.This paper is focused on the performance degradation of an LEHA induced by erosion wear.

    Performance degradation analysis is extensively applied in numerous engineering fields to evaluate safety of machine parts and equipment.14The associated theories can be divided into three categories:(a)failure physics,(b)probability statistics,and(c)artificial intelligence.The first approach studies the structural integrity of an object with respect to operating conditions and mechanical and physical properties of materials used to make the object.Probability statistics models are mainly suitable for analysis of degradation based on stress fatigue,which requires significant amount of test data.Alternately,the artificial intelligence approach is dependent on field performance degradation data.In the field of durability analysis and life prediction,numerous physical experiments are performed in order to evaluate the life and reliability of hydraulic components.Therefore,during the design stage,it is essential to analyze the overall LEHA system based on the physical nature of erosion.

    A significant number of studies have been performed to analyze the effect of erosion on system performance degradation.Fitch and Hong investigated the effects of contaminated oil on erosion in pumps and contaminant lock in servo valves,and proposed a new method to predict service life.15The occurrence of contaminant lock is accidental whereas the erosion caused by contaminant particulates is a continuous process which takes place as long as the system is operational.Vaughan et al.examined the effects of the particle size and concentration,differential pressure across the metering land,spool opening,spool surface, flow direction,as well as fluid characteristics on erosion wear.16Yang et al.adopted a gamma process to describe the internal structure degradation under erosion for electrohydraulic servo valves.17Zhang et al.presented a degradation assessment and life prediction method for electrohydraulic servo valves based on the CFD method and hydraulicsimulation.12In addition,other researchers have predicted structural wear by CFD techniques.18–20The mechanism of erosion for a ductile metal material is a micro-cutting process,which was put forward by Finnie,who also presented an analytical erosion model to calculate erosion rates.21Tilly proposed that erosion of ductile materials could occur in two stages,where the first stage is micro-cutting whereas the second stage is surface fragmentation,and found that resulting estimates gave a good correlation with experimentaldata.22Recent theoretical and experimental studies explored the effects of particle properties,impacting speed and angle,and material properties on the severity of erosion.Among those studies,the Edwards model has been widely accepted to be applicable to erosion for gassolid,liquid-solid,or gas-liquid-solid flow,where particle properties,impacting speed and angle,and material properties are taken into account.23Therefore,the Edwards model is utilized in this study because of its extensive applicability and high prediction accuracy.In addition,erosion due to the rectification mechanism in an LEHA under the influence of contaminant particles in hydraulic oil falls within the model framework.

    This paper proposes a new method for analyzing performance degradation under dynamic erosion wear.The Edwards model is utilized to obtain erosion rates due to its broad applicability and high prediction accuracy.Furthermore,erosion rates of the rectification mechanism for different port openings and different degradation stages are estimated.In addition,in order to obtain the degradation curve under dynamic erosion wear,the concentration of contaminant particles is updated at different stages of degradation.Finally,the proposed degradation model is applied to simulate the wear degradation process in an LEHA under different flow conditions,and results are compared with results from traditional wear studies.

    2.Problem description

    2.1.Working principle of the LEHA

    Fig.1 Hydraulic circuit diagram of the novel LEHA design.

    Fig.2 Functional diagram and control loop of an REHA.

    The schematic representation of the proposed LEHA architecture is shown in Fig.1,whereas Fig.2 represents the hydraulic circuit diagram and control loop of a traditional REHA.The REHA utilizes a motor and a pump to convert electrical energy to hydraulic energy,where the motor has reversible rotation in order to control the flow direction of hydraulic oil.Compared with the REHA,the proposed LEHA has the following distinct characteristics:(1)a linear resonance motor corresponding to the rotating electric machine,(2)a single direct drive piston cylinder corresponding to the traditional rotary piston pump,(3)a rectification model with passive reliance on a valve plate and active collaborative rectification dependence on a three-position four-way slide valve,and(4)direction control achieved through the phase difference between two DDPCs.The linear motor commonly operates in the resonant working mode in order to generate the maximum energy output.A Direct Drive Pump Module(DDPM)consists of DDPC1 and DDPC2,which are characterized with collaborative rectification.A DDPC consists of a cylinder and a valve,where a piston rod and a valve spool are integrated,and is driven directly by the linear resonating motor.X and Y are the interface ports of the DDPCs:one is the inlet and the other is the outlet.The deflection angle of aircraft’s control surface is θ.As a result,the novel pump is more compact and has better controllability.The theoretical research and experimental results conducted on the prototype valve have demonstrated that the DDPC-dis charged kinematic flow volume depends on the phase difference between the two cells.The output flow in one cycle from a DDPM port is represented by

    where V is the output flow in one cycle from a DDPM port,Acis the piston’s effective area,Spis the resonance amplitude of a DDPC,and φ is the phase difference between the two cells.

    The kinematic volume efficiency is maximum for φ=±π/2,and zero for φ =0 or φ = ±π.

    2.2.Dominant mode of gradual failure

    At an earlier stage,a Failure Mode and Effect Analysis(FMEA)was performed for the LEHA.24The results show that resonant spring fatigue damage and rectification mechanism erosion are the two major progressive failure modes affecting the performance of the LEHA.The LEHA failure is defined as a condition when a performance indicator cannot be maintained within a specified range.The effect of rectification mechanism erosion is shown in Fig.3,where P,T,A,and B are the hydraulic oil ports of the three-position four-way valve.The figure indicates possible fluid flow patterns in the rectification slide valve when port P is connected with X and port T is connected with Y.It can be seen that the rectification valve is characterized by reciprocating flushing which is different from single-direction flushing which occurs in a traditional valve.The rectification valve at the middle position is equivalent to a piston pump at the upper or lower dead point.Consequently,wear of the internal structure can cause an increase in internal leakage and result in power loss.

    The major sources of initial contamination include particulates resulting from the manufacturing process,dust inclusions from storage,and contaminant particles mixed in oil.These contaminants can impact the system in three different ways:(a)three-body abrasion,(b)particle impingement erosion,and(c)motion impediment.15Motion impediment is rare in practice and caused by relatively large particles.However,the erosion problem is common,and can become serious as the process has a tendency of having positive feedback.

    Fig.3 Fluid flow in the rectification slide valve when P is connected with X and T is connected with Y.

    The major impact of erosion is reflected on the performance of the LEHA.Namely,there are two ways in which the power/weight ratio of the LEHA can be affected:one is a change in the resonant frequency and the other one is a change in the amplitude.The resonant amplitude is subject to the size of the DDPM,which makes modification in the frequency become the most effective method.The initial clearance between the valve spool and sleeve in new valves is generally in the range of 3–5 μm,and the corner radius of edges is generally below 0.5 μm.Throughout their lifecycles,valves will undergo erosion which will cause an increase in internal leakage,which,on the other hand,will ultimately degrade the system performance.

    3.Mathematical model of the degradation process

    The erosion-induced performance degradation process can be analyzed in three consecutive steps:(1)identification of the potential erosion location and determination of the fluid dynamic boundary,(2)determination of the structural degradation after erosion,and(3)prediction of the erosion rate.

    3.1.Stress-strain equation of the dynamic boundary

    The erosion wear is analyzed in both valve and piston regions as indicated in Fig.4.Pv1,Tv1,Av1,and Bv1are the hydraulic oil ports of the valve in DDPC1.Ap2and Bp2are the cylinder chambers in DDPC2.The displacement of the mover is x1in DDPC1,and the velocity is v1.Similarly,the displacement of the mover is x2in DDPC2,and the velocity is v2.The positive direction is the cylinder pointing to the valve.When the spool in DDPC1 or DDPC2 is at the zero displacement,i.e.,x1=0 or x2=0,there is no connection between P,T and A,B.The piston determines the flow velocity and flow rate through the valve inlet,and in turn,the valve transfers the load pressure to the piston.

    The continuity equations of hydraulic chambers Ap2and Bp2can be written as

    where PdA2is the hydraulic pressure in the left chamber,PdB2is the hydraulic pressure in the right chamber,βeis the effective bulk modulus(including hydraulic oil,connecting pipes,and the cylinder’s mechanical flexibility),VdA2is the volume of the left chamber(including valves and connecting pipes),VdB2is the volume of the right chamber(including valves and connecting pipes),Cipis the cylinder’s internal leakage coefficient,Cepis the cylinder’s external leakage coefficient,qdA2outis the out flow from the cylinder’s left chamber,qdB2inis the in flow to the cylinder’s right chamber,VdA02is the initial volume of the left chamber,and VdB02is the initial volume of the right chamber.

    Fig.4 Schematic diagram of the reciprocating system.

    The fluid flow rate Q from the driven hydraulic cylinder chamber to the valve entrance can be described by the following orifice flow equation:

    where Cdis the flow coefficient of the orifice,Aeis the effective area of the orifice,Pfis the hydraulic pressure in front of the orifice,Pbis the hydraulic pressure behind the orifice,ρ is the fluid density,vpais the average velocity of fluid at the valve inlet,and Apis the cross-sectional area of the valve inlet.

    3.2.Continuous flow model

    The continuous phase modeling uses Reynolds averaged Navier-Stokes equations.Because of enhanced accuracy for rapidly strained and swirling flows,the RNG k-ε model is selected for turbulence,where the continuity equation is described as

    where v is the velocity of fluid.

    3.3.Structural degradation model

    According to a significant number of published results,erosion has the greatest impacts on the radial clearance between the spool and the sleeve,and on the throttle edge radii.An ideal geometry calls for the throttle edges to be perfectly square with no rounding.11–13,25,26The erosion process will cause the throttle edges to wear and develop radii,resulting in an increase of the radial clearance between the spool and the sleeve,as shown in Fig.5.The radial clearance between the spool and the sleeve has a small initial value because of the fluid’s lubrication.

    According to the Pythagorean theorem,the gap between the spool and the sleeve can be determined as

    where δ1is the radial clearance between the spool and the sleeve,δ2is the flow clearance between the spool and the sleeve,r1is the radius of the spool’s throttle edge,r2is the radius of the sleeve’s throttle edge.

    It is assumed that r=r1=r2,since the difference between r1and r2is relatively small.

    Fig.5 Throttle edge after wear.

    3.4.Wear volume model

    The geometry of the cylindrical spool changes during erosion wear as indicated in Fig.6.Assuming the same wear conditions between contact surfaces,then the degrees of wear of the spool and the sleeve can be assumed to be equal.In addition,it is assumed that the wear rate of each orifice is the same.Since the edge radius,r,is small compared to the spool radius,Rs,the radius of the spool,Rs,is considered to remain constant throughout the process.

    The longitudinal cross section of the wear-free spool is assumed to be rectangular with square edges.As wear progresses,the spool edges start to develop a quarter-round geometry of radius r.The total wear volume of a valve’s orifice can be obtained by the following equation:

    where Vois the total wear volume of a valve’s orifice,Vsis the wear volume of a single throttle edge,r is the equivalent radius of the throttling edge,c is the perimeter of the throttle edge,and Rsis the radius of the spool’s boss.

    The total wear volume of the valve can be obtained by the time integration of the erosion wear rate as follows:

    where ρpis the contaminated particle’s density,seis the eroded area,and ρe(t)is the erosion rate function.

    3.5.Contaminant particle concentration model

    During a lifetime of LEHA operation,metal particles generated by wear end up in hydraulic oil.Since the LEHA is a closed system,it is reasonable to assume that all particles produced by wear are accumulated in the working hydraulic oil.The contaminant particulate concentration can be expressed by the following mass conservation equation:

    where η0is the initial particle mass flow rate,η is the particle mass flow rate with the corresponding wear volume Vo,and Vhis the hydraulic oil volume in the LEHA.

    Fig.6 Geometries of worn throttle edges.

    3.6.Discrete phase model

    The dynamics equation of a particle in a Lagrangian reference frame can be expressed as

    where u is the fluid velocity,upis the particle velocity,and upis the scalar form of the particle velocity.FD(u-up)represents the force load on unit particle mass,g is the acceleration of gravity,and Fxis the additional force including Saffman’s lift force,the virtual mass force,the pressure gradient force,the thermophoretic force,the forces in moving reference frames,and the Brownian Force.

    where dpis the particle’s diameter,Re is the relative Reynolds number,μ is the dynamic viscosity of the fluid,and CDis the drag coefficient,which is defined by Haider and Levenspiel27as

    where Φ is the particle shape factor,s is the surface area of a sphere with the same volume as the particle,and S is the actual surface area of the particle,and b1,b2,b3,and b4are coefficients related to Φ.

    It can be assumed that the particle volume fraction or concentration is very low and that the collision between particles can be neglected.

    3.7.Particle-wall collision-rebound model

    When particles impinge the valve wall at a certain angle and speed,they rebound at a certain angle and speed,as shownin Fig.7.The energy loss during collision can be determined from the impulse-momentum principle where normal and tangential restitution coefficients are employed,which are introduced as

    Fig.7 Schematic diagram of particle-to-wall collision-rebound.

    where εNis the normal restitution coefficient,εTis the tangential restitution coefficient,vp1,nand vp1,tare the particle’s normal and tangential velocities with respect to the wall prior to collision,and vp2,nand vp2,tare the particle’s normal and tangential velocities after collision.

    It has been demonstrated by Grant and Forder et al.that restitution coefficients are dependent on the particle impingement angle,αp.28,29Furthermore,Forder et al.provided the expressions of restitution coefficients as follows:

    3.8.Semi-empirical material removal model

    The Edwards model is commonly used to determine the erosion rate for certain particle concentration and wall boundaries,and can be expressed as follows:

    where Npis the number of particles that impact the area represented by Aface,mpis a particle’s mass,C(dp)is the function of a particle’s diameter,vpis the particle velocity relative to the wall,b(vp)is the function of the relative particle velocity,Afaceis the area of the cell face on the wall,and f(αp)is the function of the impact angle,which is determined as

    4.Simulation of erosion wear

    A well-established three-step erosion prediction approach using the Eulerian-Lagrangian approach is used in this work to simulate the process.The model,which has been adopted by many authors,consists of the following steps:(1) flow field prediction using viscous or inviscid models,(2)particle trajectory calculation by modeling discrete particle dynamics,and(3)determination of erosion wear using a material removal model.The development of the Eulerian-Lagrangian model of fluid particle flows is well documented and will not be discussed in detail.However,certain aspects of the erosion model applied in this study are discussed below.

    The model of fluid flow is based on the Reynolds averaged Navier-Stokes equations,while the RNG k-ε model,due to its enhanced accuracy for rapidly strained and swirling flows,is used to describe turbulence.The particle trajectory and its velocity are determined by solving a set of ordinary differential equations formed by the Lagrangian method.A simulation is performed with ANSYS Workbench 15.0.

    The simulation is based on the implementation of the flow chart shown in Fig.8.Firstly,the spool position and the concentration of contaminant particles are determined.A dimensional fluid domain model is established next,followed by meshing,parameter initialization,and calculation of erosion wear.The erosion wear is recalculated when the spool position or contaminated particles’concentration changes.Finally,an erosion degradation curve is obtained by integrating the erosion rate with respect to time.

    The main structural parameters of the LEHA’s rectification mechanism are given in Table 1.The initial flow rate of hydraulic oil with contaminant particles can be obtained by 7th level.12It is assumed that the particles are pure iron,the reciprocating frequency is 100 Hz,and the opening size is±3 mm.The dynamic fluid boundary conditions can be obtained from Eqs.(7)and(8).In order to achieve high computational effectiveness,an optimal algorithm is used to balance the grid quality and calculation time.The resulting computational domain mesh is shown in Fig.9.The standard case simulation parameters are given in Table 2.

    Fig.8 Algorithm for erosion prediction.

    Table 1 Main structural parameters of the LEHA.

    Fig.9 Computational domain and mesh(when the valve opening is 50%).

    In addition,the conventional SIMPLE algorithm is applied to couple pressure and velocity via a second-order upwind scheme for the momentum.The convergence criterion is set so that the residual is smaller than 10-3.

    The simulation results representing the erosion rate contours of the spool valve are shown in Fig.10.The results indicate that the edge of the throttle orifice is experiencing the most serious erosion.The results are consistent with the published results focused on erosion of EHSVs,throttle valves,and needle valves.11–13,17,25,30,31The particle path lines and associated velocity magnitudes are shown in Fig.11.The areas with high-speed pro files are evident in the orifice edge region,where the particle velocity can be in excess of 100 m/s.Having a large number of particles with significant kinetic energy,directly impinging the orifice,will result in a severe erosion of the surface.

    5.Results and discussion

    5.1.Effect of the opening gap on the erosion rate

    The average erosion rate of a sharp edge as a function of the valve opening(dop),during the early stage of degradation,is shown in Fig.12.Correspondingly,Fig.13 depicts the erosion rate distribution in the throttling edge for different valve opening sizes.The results are for valve port opening between 0 mm and 3 mm.It can be observed that the maximum average erosion rate occurs when the EHSV opening is 0.5 mm,and that the average erosion rate decreases rapidly for opening sizes between 0.5 mm and 1.0 mm.

    Table 2 Standard case simulation parameters.

    Fig.10 Erosion rate contours of the valve.

    Fig.11 Particles’path lines of the velocity magnitude.

    Fig.12 Erosion wear rate under different valve opening sizes.

    5.2.Effect of the inlet velocity on the erosion rate

    In this section,the effect of the inlet velocity(vin)on flow erosion is studied.Simulations are performed for different inlet velocities while keeping other parameters the same as those in the standard case.The erosion rate distribution in the throttling edge as a function of the inlet velocity is shown in Fig.14.It is apparent that an increase in the inlet velocity increases the erosion severity of both the sleeve and the spool.It is also apparent that the erosion is maximum at the maximumvelocity locations,which occurs at the gap between the edges of the valve spool and sleeve.The simulation results indicate that the erosion zone on the spool wall increases simultaneously with an increase in the severity of erosion as the inlet velocity increases from 45 m/s to 65 m/s,while the erosion in the sleeve is consistently spread over the entire throttle side.

    The maximum erosion rate on the spool occurs at an inlet velocity of 65 m/s and is about 8.9 times higher than that at an inlet velocity of 45 m/s,while the maximum erosion rate along the sleeve throttle side for an inlet velocity of 65 m/s is 6.1 times higher than that for an inlet velocity of 45 m/s.

    5.3.Effect of the inlet pipe diameter on the erosion rate

    The erosion rates for different inlet pipe diameters but for constant pressure and opening are shown in Fig.15.It can be observed that the maximum erosion wear rate initially increases and then decreases with an increase of the inlet diameter.The magnitude of the maximum erosion rate,for both the spool and the sleeve,is largest at din=4 mm.The erosion rate of the valve spool for an inlet pipe diameter of din=2.5 mm is approximately 8%,whereas for din=4 mm,it is 7.5%.

    The results indicate that increasing the diameter of the inlet pipe results in an increase of the flow through the gap between the spool and the sleeve,which ultimately leads to an increase in the erosion severity.

    5.4.Effect of the load pressure on the erosion rate

    The results indicating the erosion rate distribution on the valve spool and along the sleeve for different levels of load pressure(Pl)are shown in Fig.16.When the load pressure increases from 0 MPa to 5 MPa,the erosion wear rate decreases slowly,whereas the erosion rate decreases rapidly with the load pressure increasing from 5 MPa to 8 MPa.The maximum erosion rates on the spool and along the sleeve throttling edges for a load pressure of 0 MPa are respectively 3.21 times and 3.46 times higher than those for a load pressure of 5 MPa.The erosion rate decreases rapidly when the load pressure increases to 8 MPa.

    A lower load pressure causes a greater pressure drop on the valve path,resulting in a greater velocity of particles impinging on the valve wall,which increases erosion wear.

    5.5.Effect of the particle concentration on the erosion rate

    Fig.13 Erosion rate distribution in the throttling edge under different valve opening sizes.

    The results indicate that the particle concentration (cp),expressed as the particle mass flow rate,plays a significant role in affecting the erosion,as shown in Fig.17.In addition,increasing the particle mass flow rate results in an increase of the erosion-dominated zone as well as the severity of erosion.As the mass flow rate increases from 1.78×10-7kg/s to 3.56×10-7kg/s,the erosion zone extends to the entire throttling edge.The peak erosion rate of the spool at a particle mass flow rate of 3.56×10-7kg/s is approximately 2.64 times of that at a flow rate of 1.78×10-7kg/s.Whereas,the maximum erosion rate along the sleeve at a particle mass flow rate of 3.56×10-7kg/s is about 2.21 times of that at a flow rate of 1.78×10-7kg/s.

    The higher the particle mass flow rate is,the more particles impinge on the spool and the sleeve per unit time.The result is that the same region is impacted by a larger number of particles.

    5.6.Effect of the particle diameter on the erosion rate

    The results indicate that the particle diameter clearly affects the erosion rates of both the spool and the sleeve,as illustrated in Fig.18.As the particle diameter increases from 1 μm to 3 μm,the erosion rate increases rapidly;however,as the particle diameter increases from 3 μm to 20 μm,the erosion rate decreases slowly.The maximum erosion rates of the spool and along the sleeve at a particle diameter of 3 μm are respectively 2.41 and 2.35 times higher than those for a particle diameter of 1 μm.Furthermore,the maximum erosion rates of the spool and along the sleeve for a particle diameter of 3 μm are respectively 16.23 and 15.64 times of those for a particle diameter of 15 μm.The results indicate that the particle diameter has a significant effect on erosion.

    Fig.14 Erosion rate distribution in the throttling edge under different inlet velocities.

    Fig.15 Erosion rate distribution in the throttling edge under different inlet pipe diameters.

    The underlying assumption is that the mass flow rate of particles is constant.Consequently,the larger the diameter of a single particle,the greater the mass of a single particle,and the smaller the number of particles.When the number of particles in the unit volume decreases,the probability of the particles impacting the spool or the sleeve is reduced with the number of particles passing through the gap between the valve spool and the sleeve.For the same velocity,a particle of a larger mass has greater kinetic energy than that of a particle of a smaller mass.As a result,assuming the same incidence angle,a particle of a larger mass causes more severe erosion than a particle of a smaller mass.As the particle diameter increases from 1 μm to 3 μm,the particle inertial force plays a dominant role in the erosion process,while for particle sizes between 3 μm and 15 μm,the number of particles plays a primary role.

    Fig.16 Erosion rate distribution in the throttling edge under different load pressures.

    Fig.17 Erosion rate distribution in the throttling edge under different particle concentrations.

    6.Conclusions

    The paper introduces the working principle of a novel LEHA design,which is an integrated closed system without filters.A study was performed to analyze the effects of contaminant particles,caused by the system erosion wear,on the internal leakage of the valve and ultimately on its performance.According to failure analysis,erosion is the dominant progressive failure mode of the LEHA.Since the system does not include filters,the mass of contaminant particles accumulates over the life of the valve,which can lead to progressive failure of the system.Simulation results provide insight into the effects of valve opening on the erosion rate and the relationship between the erosion rate and structural degradation.A comparison between the simulation results and published wear erosion results for EHSVs indicates that the proposed analysis approach is effective.The following conclusions can be made based on the simulation results:(1)the erosion wear dominates the service life of the LEHA(not considering accidental failure),(2)the throttle edge erodes fastest,caused by dissipation of the kinetic energy of a large number of high-speed particles at that location,and(3)the erosion of the LEHA is an accelerated degradation process.The proposed approach can be used to evaluate the service life of the LEHA and to further guide the design of the LEHA’s rectification valve structure.An experimental setup is under construction,which will provide verification of the proposed analysis model,and enable further studies of erosion phenomena in the LEHA and their effects on performance degradation of the LEHA.

    Fig.18 Erosion rate distribution in the throttling edge under different particle diameters.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China (Nos. 51620105010, 51675019,51575019),the National Basic Research Program of China(No.2014CB046402),and the Program 111 of China.The authors would also like to thank the reviewers and the editor for their comments and suggestions that led to significant improvement in the quality of the paper.

    Appendix A.Supplementary material

    Supplementary data associated with this article can be found,in the online version,at https://doi.org/10.1016/j.cja.2017.12.002.

    References

    1.Van den Bossche D.The A380 flight control electrohydrostatic actuators,achievements and lessons learnt.25th International congress of the aeronautical sciences;2006 Sep 3–8;Hamburg,Germany.New York:Curran Associates Inc.;2006.p.3383–90.

    2.MaréJ-C,Fu J.Review on signal-by-wire and power-by-wire actuation for more electric aircraft.Chin J Aeronaut 2017;30(3):857–70.

    3.Naayagi RT.A review of more electric aircraft technology.In International conference on energy efficient technologies for sustainability;2013 April 10–12;Nagercoil,India.Piscataway,NJ:IEEE Computer Society;2013.p.750–3.

    4.Wang SP,Cui XY,Shi J,Tomovic MM,Jiao ZX.Modeling of reliability and performance assessment of a dissimilar redundancy actuation system with failure monitoring.Chin J Aeronaut 2016;29(3):799–813.

    5.Anderson EH,Lindler JE,Regel brugge ME.Smart material actuator with long stroke and high power output.43rd AIAA/ASME/ASCE/AHS/ASC structures,structuraldynamics,and materials conference;2002 April 22–25;Denver,US.Reston,VA:AIAA Inc.;2002.

    6.Lindler JE,Anderson EH,Regelbrugge ME.Design and testing of piezoelectric-hydraulic actuators.Smart structures and materials 2003:industrial and commercial applications of smart structures technologies;2003 March 4–6;San Diego,US.Bellingham,WA:SPIE;2003.p.96–107.

    7.J?nker P,Claeyssen F,Grohmann B,Christmann M,Lorkowski T,LeLetty R,et al.New actuators for aircraft and space applications.11th International conference on new actuators;2008 June 9–11;Bremen,Germany.Bremen:ACTUATOR Interest Group;2008 p.325–30.

    8.Zhang R,Jiao ZX,Yan L,Wu S,Liang HS,Zheng JL.Design and numerical simulation of a continuable and bidirectional piezohydraulic servo pump.2011 International conference on fluid power and mechatronics;2011 August 17–20;Beijing China.Piscataway,NJ:IEEE Computer Society;2011 p.498–504.

    9.Kang RJ,Jiao ZX,Wang SP,Chen LS.Design and simulation of electro-hydrostatic actuator with a built-in power regulator.Chin J Aeronaut 2009;22(6):700–6.

    10.Li Y,Jiao ZX,Yan L,Dong WH.Conceptual design and composition principles analysis of a novel collaborative rectification structure pump.J Dyn Syst Meas Control Trans ASME 2014;136(5):054507–054507-8.

    11.Nystad BH,Gola G,Hulsund JE,Roverso D.Technical condition assessment and remaining useful life estimation of choke valves subject to erosion.Annual conference of the prognostics and health management society;2010 October 13–16;Portland,OR,US.College Park,MD:Prognostics and Health Management Society;2010.p.11–3.

    12.Zhang K,Yao JY,Jiang TM.Degradation assessment and life prediction of electro-hydraulic servo valve under erosion wear.Eng Fail Anal 2014;36(1):284–300.

    13.Fang X,Yao JY,Yin XZ,Chen X,Zhang CH.Physics-offailure models of erosion wear in electrohydraulic servovalve,and erosion wear life prediction method.Mechatronics 2013;23(8):1202–14.

    14.Zhang XL,Chen XF,Li B,He ZJ.Review of life prediction for mechanical majorequipments.J Mech Eng 2011;47(11):100–16.

    15.Fitch EC,Hong IT.Hydraulic system design for service assurance.Stillwater,OK:BarDyne Inc.;2004.p.125–78.

    16.Vaughan ND,Pomeroy PE,Tilley DG.The contribution of erosive wear to the performance degradation of sliding spool servovalves.Proc Inst Mech Eng Part J:J Eng Tribol 1998;212(6):437–51.

    17.Yang YJ,Peng WW,Meng DB,Zhu SP,Huang HZ.Reliability analysis of direct drive electrohydraulic servo valves based on a wear degradation process and individual differences.Proc Inst Mech Eng Part O:J Risk Reliab 2014;228(6):621–30.

    18.Frawley P,Corish J,Niven A,Geron M.Combination of CFD and DOE to analyse solid particle erosion in elbows.Int J Comput Fluid Dyn 2009;23(5):411–26.

    19.Tang P,Yang J,Zheng JY,Wong I,He SZ,Ye JJ,et al.Failure analysis and prediction of pipes due to the interaction between multiphase flow and structure.Eng Fail Anal 2009;16(5):1749–56.

    20.El-Behery SM,Hamed MH,Ibrahim KA,El-Kadi MA.CFD evaluation of solid particles erosion in curved ducts.J Fluids Eng Trans ASME 2010;132(7):071303.

    21.Finnie I.Some observations on the erosion of ductile metals.Wear 1972;19(1):81–90.

    22.Tilly GP.A two stage mechanism of ductile erosion.Wear 1973;23(1):87–96.

    23.Edwards JK.Development,validation,and application of a threedimensional,CFD-based erosion prediction procedure[dissertation].Tulsa(OK):The University of Tulsa;2001.

    24.Li Y,Wang SP,Shi J,Tomovic MM.Dynamic contaminated particles concentration-based degradation model of linear electrohydrostatic actuator.2017 Annual reliability and maintainability symposium;2017 January 23–26;Orlando,FL,US.Piscataway,NJ:IEEE Inc.;2017.

    25.Pan XD,Wang GL,Lu ZS,Liu ZH.Simulation research on effect of diametral clearance of spool valve to valve orifice discharge characteristic.Key Eng Mater 2009;392–394:184–8.

    26.Zhang H,Xiong SP,Liang YW,Xiong XY.Analyses of erosion wear characteristic and structure research on hydraulic valve.J China Coal Soc 2008;33(2):214–7.

    27.Haider A,Levenspiel O.Drag coefficient and terminal velocity of spherical and nonspherical particles.Powder Technol 1989;58(1):63–70.

    28.Grant G,Tabakoff W.Erosion prediction in turbomachinery resulting from environmental solid particles.J Aircr 1975;12(5):471–8.

    29.Forder A,Thew M,Harrison D.A numerical investigation of solid particle erosion experienced within oil field control valves.Wear 1998;216(2):184–93.

    30.Wallace MS,Dempster WM,Scanlon T,Peters J,McCulloch S.Prediction of impact erosion in valve geometries.Wear 2004;256(9–10):927–36.

    31.Zhu HJ,Pan Q,Zhang WL,Feng G,Li X.CFD simulations of flow erosion and flow-induced deformation of needle valve:effects of operation,structure and fluid parameters.Nucl Eng Des 2014;273(1):396–411.

    久久精品熟女亚洲av麻豆精品| 99国产综合亚洲精品| 大片免费播放器 马上看| 十分钟在线观看高清视频www| 亚洲av.av天堂| 欧美人与性动交α欧美精品济南到 | 欧美另类一区| 国产国拍精品亚洲av在线观看| 国产国拍精品亚洲av在线观看| 一级a做视频免费观看| 青春草亚洲视频在线观看| 青春草亚洲视频在线观看| 美女内射精品一级片tv| 日日摸夜夜添夜夜爱| 中文字幕av电影在线播放| 免费观看无遮挡的男女| 欧美 亚洲 国产 日韩一| 天堂俺去俺来也www色官网| 精品亚洲乱码少妇综合久久| 国产日韩欧美亚洲二区| 秋霞伦理黄片| 97超碰精品成人国产| 免费看光身美女| 人人妻人人澡人人爽人人夜夜| 精品亚洲成a人片在线观看| 大香蕉97超碰在线| 高清毛片免费看| 亚洲第一av免费看| 日韩成人伦理影院| 肉色欧美久久久久久久蜜桃| 欧美亚洲 丝袜 人妻 在线| 国产伦精品一区二区三区视频9| 亚洲av福利一区| 国精品久久久久久国模美| 成人18禁高潮啪啪吃奶动态图 | 久久99热6这里只有精品| 免费黄网站久久成人精品| 久热久热在线精品观看| 国产av一区二区精品久久| 亚洲性久久影院| 伊人亚洲综合成人网| 全区人妻精品视频| 交换朋友夫妻互换小说| 最黄视频免费看| 啦啦啦在线观看免费高清www| 18禁在线播放成人免费| 激情五月婷婷亚洲| 中文精品一卡2卡3卡4更新| 亚洲无线观看免费| 五月开心婷婷网| 亚洲国产精品一区二区三区在线| 天堂中文最新版在线下载| 大话2 男鬼变身卡| 日产精品乱码卡一卡2卡三| 国产成人免费观看mmmm| 国产高清国产精品国产三级| 春色校园在线视频观看| 午夜91福利影院| 久久ye,这里只有精品| 黑人高潮一二区| 中国国产av一级| 王馨瑶露胸无遮挡在线观看| 精品人妻熟女av久视频| 少妇熟女欧美另类| 亚洲av不卡在线观看| 赤兔流量卡办理| 少妇高潮的动态图| 日韩av在线免费看完整版不卡| 亚洲精品aⅴ在线观看| 九九久久精品国产亚洲av麻豆| 少妇的逼好多水| 在线观看三级黄色| 18禁在线播放成人免费| 一区二区三区免费毛片| 高清欧美精品videossex| 丝袜脚勾引网站| 国产成人aa在线观看| 日本黄色片子视频| 日本91视频免费播放| 欧美+日韩+精品| 国产免费一级a男人的天堂| 69精品国产乱码久久久| 亚洲国产av影院在线观看| 丰满迷人的少妇在线观看| 人人妻人人添人人爽欧美一区卜| 18禁在线无遮挡免费观看视频| 一区二区三区乱码不卡18| 99热6这里只有精品| 蜜桃久久精品国产亚洲av| 亚洲,一卡二卡三卡| 欧美 亚洲 国产 日韩一| 国产成人aa在线观看| 99热这里只有精品一区| 涩涩av久久男人的天堂| 麻豆成人av视频| 亚洲精品亚洲一区二区| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久av不卡| 曰老女人黄片| 亚洲国产精品国产精品| 国产伦理片在线播放av一区| 99国产综合亚洲精品| 这个男人来自地球电影免费观看 | 亚洲精品久久久久久婷婷小说| 欧美日韩在线观看h| 激情五月婷婷亚洲| 熟女人妻精品中文字幕| 日韩电影二区| 国产成人精品无人区| 欧美人与性动交α欧美精品济南到 | 免费av中文字幕在线| 日日摸夜夜添夜夜爱| 久久精品国产亚洲网站| 欧美97在线视频| 亚洲成色77777| 免费播放大片免费观看视频在线观看| 国产黄片视频在线免费观看| 久久精品久久久久久久性| 美女中出高潮动态图| 午夜免费男女啪啪视频观看| 午夜福利在线观看免费完整高清在| 91国产中文字幕| 国产一区亚洲一区在线观看| 日本爱情动作片www.在线观看| 亚洲欧美中文字幕日韩二区| 亚洲国产欧美日韩在线播放| 亚洲图色成人| 亚洲性久久影院| 中文天堂在线官网| 哪个播放器可以免费观看大片| 夫妻午夜视频| 亚洲欧美日韩卡通动漫| 91精品国产国语对白视频| 久久午夜福利片| 免费观看的影片在线观看| 国产亚洲精品久久久com| 一区二区三区四区激情视频| 国产精品久久久久久精品电影小说| 亚洲国产欧美日韩在线播放| 在线播放无遮挡| 成年人免费黄色播放视频| 在线看a的网站| av视频免费观看在线观看| 国产探花极品一区二区| 十八禁高潮呻吟视频| 国产不卡av网站在线观看| 日韩,欧美,国产一区二区三区| 亚洲成色77777| 欧美亚洲 丝袜 人妻 在线| 桃花免费在线播放| 亚洲综合色惰| 免费观看无遮挡的男女| 亚洲人与动物交配视频| 欧美成人午夜免费资源| 亚洲无线观看免费| 日韩视频在线欧美| 少妇精品久久久久久久| 涩涩av久久男人的天堂| 亚洲美女搞黄在线观看| 又黄又爽又刺激的免费视频.| 国产精品人妻久久久久久| 看十八女毛片水多多多| 黄色一级大片看看| 天天操日日干夜夜撸| 国产女主播在线喷水免费视频网站| 久久热精品热| 91aial.com中文字幕在线观看| 麻豆成人av视频| 永久免费av网站大全| 国产一区二区在线观看日韩| 国产国语露脸激情在线看| 日韩av免费高清视频| 性色avwww在线观看| 嘟嘟电影网在线观看| 少妇被粗大的猛进出69影院 | 老司机影院成人| 99久国产av精品国产电影| 热re99久久精品国产66热6| 中文字幕久久专区| 日本黄色日本黄色录像| 一级,二级,三级黄色视频| 国产成人免费观看mmmm| 日本av手机在线免费观看| 狂野欧美激情性bbbbbb| 黄色欧美视频在线观看| 九九在线视频观看精品| 国产精品一国产av| 国产精品三级大全| av女优亚洲男人天堂| 久久精品国产鲁丝片午夜精品| 久久精品人人爽人人爽视色| 26uuu在线亚洲综合色| 亚洲精品久久午夜乱码| 另类精品久久| 欧美激情国产日韩精品一区| 七月丁香在线播放| 欧美日韩亚洲高清精品| 午夜日本视频在线| 中文字幕精品免费在线观看视频 | 亚洲一区二区三区欧美精品| 啦啦啦视频在线资源免费观看| 午夜影院在线不卡| 国产国语露脸激情在线看| 亚洲精品日本国产第一区| 99re6热这里在线精品视频| 免费不卡的大黄色大毛片视频在线观看| 色婷婷av一区二区三区视频| 少妇猛男粗大的猛烈进出视频| 久久国产精品男人的天堂亚洲 | 日韩成人av中文字幕在线观看| videos熟女内射| 寂寞人妻少妇视频99o| xxx大片免费视频| 国产精品久久久久久av不卡| 狂野欧美白嫩少妇大欣赏| 在线亚洲精品国产二区图片欧美 | 国产一区亚洲一区在线观看| av国产久精品久网站免费入址| 黑人巨大精品欧美一区二区蜜桃 | 大话2 男鬼变身卡| 国产精品一国产av| 黄片播放在线免费| 午夜激情av网站| 欧美一级a爱片免费观看看| 欧美三级亚洲精品| 亚洲国产欧美在线一区| 热99国产精品久久久久久7| 老女人水多毛片| 人妻夜夜爽99麻豆av| 成人漫画全彩无遮挡| 18禁观看日本| 国产高清不卡午夜福利| 亚洲精品日韩在线中文字幕| 黑人高潮一二区| av不卡在线播放| 男人添女人高潮全过程视频| 亚洲精品aⅴ在线观看| 国产欧美日韩综合在线一区二区| 亚洲精品久久午夜乱码| 蜜桃国产av成人99| 亚洲精品美女久久av网站| av福利片在线| 91精品伊人久久大香线蕉| 亚洲精品国产av蜜桃| 国产永久视频网站| 亚洲精品国产av成人精品| 伦精品一区二区三区| 一级片'在线观看视频| 亚洲国产成人一精品久久久| 亚洲精品aⅴ在线观看| 亚洲高清免费不卡视频| 人妻制服诱惑在线中文字幕| 国国产精品蜜臀av免费| 久久久国产欧美日韩av| 国产欧美日韩综合在线一区二区| 少妇熟女欧美另类| 男女啪啪激烈高潮av片| av女优亚洲男人天堂| 亚洲成色77777| 免费大片18禁| 国产精品人妻久久久久久| 国产精品人妻久久久影院| 成人二区视频| 男女边摸边吃奶| 国产亚洲精品久久久com| 岛国毛片在线播放| 91在线精品国自产拍蜜月| 亚洲人成网站在线观看播放| 精品午夜福利在线看| 亚洲美女搞黄在线观看| 哪个播放器可以免费观看大片| 插阴视频在线观看视频| 色94色欧美一区二区| 各种免费的搞黄视频| 欧美 日韩 精品 国产| 亚洲av中文av极速乱| 99九九在线精品视频| 美女中出高潮动态图| 视频区图区小说| 日本黄色片子视频| 亚洲精品,欧美精品| 简卡轻食公司| 天堂中文最新版在线下载| 亚洲色图综合在线观看| 不卡视频在线观看欧美| 18禁在线播放成人免费| 国产伦理片在线播放av一区| 亚洲内射少妇av| 极品少妇高潮喷水抽搐| 午夜av观看不卡| 亚洲精品第二区| 久久国产精品大桥未久av| 一区二区三区四区激情视频| 麻豆成人av视频| 免费av中文字幕在线| 午夜影院在线不卡| 大码成人一级视频| 久久久a久久爽久久v久久| 草草在线视频免费看| 丝袜在线中文字幕| 亚洲美女黄色视频免费看| a 毛片基地| 亚洲成人av在线免费| av国产精品久久久久影院| 成年女人在线观看亚洲视频| 国产精品人妻久久久影院| av黄色大香蕉| 一区二区av电影网| 免费黄网站久久成人精品| 男人操女人黄网站| 婷婷色av中文字幕| 亚洲天堂av无毛| 人妻系列 视频| 天天躁夜夜躁狠狠久久av| 欧美 亚洲 国产 日韩一| 久久国产精品大桥未久av| 免费观看无遮挡的男女| 女性生殖器流出的白浆| 美女福利国产在线| 日韩熟女老妇一区二区性免费视频| 日韩成人av中文字幕在线观看| 国产精品人妻久久久久久| 色哟哟·www| 免费黄频网站在线观看国产| 日本wwww免费看| 欧美日韩亚洲高清精品| 久久精品久久久久久噜噜老黄| 在线观看www视频免费| 最新的欧美精品一区二区| 曰老女人黄片| 日韩av不卡免费在线播放| 色婷婷久久久亚洲欧美| 人成视频在线观看免费观看| 99久久精品国产国产毛片| 欧美最新免费一区二区三区| 少妇猛男粗大的猛烈进出视频| 亚洲国产日韩一区二区| 91久久精品国产一区二区成人| 国产成人freesex在线| 亚洲国产精品一区三区| 国产精品99久久久久久久久| 99九九线精品视频在线观看视频| 日日摸夜夜添夜夜添av毛片| 男女国产视频网站| 国产精品不卡视频一区二区| 久久影院123| 国产成人午夜福利电影在线观看| 日韩强制内射视频| 在线播放无遮挡| 国产亚洲最大av| 午夜免费观看性视频| 又粗又硬又长又爽又黄的视频| 久久国产精品男人的天堂亚洲 | 国产精品久久久久成人av| 亚洲成人av在线免费| 王馨瑶露胸无遮挡在线观看| 成人亚洲欧美一区二区av| 王馨瑶露胸无遮挡在线观看| 成人亚洲欧美一区二区av| 国产免费一区二区三区四区乱码| 99精国产麻豆久久婷婷| 一本久久精品| 两个人的视频大全免费| 最近的中文字幕免费完整| 伊人亚洲综合成人网| 国产精品一区www在线观看| 午夜福利网站1000一区二区三区| 国产精品久久久久久精品电影小说| 人妻夜夜爽99麻豆av| 午夜福利在线观看免费完整高清在| 丝袜美足系列| 夫妻性生交免费视频一级片| videossex国产| 黑丝袜美女国产一区| 18禁在线无遮挡免费观看视频| 亚洲国产日韩一区二区| 日韩大片免费观看网站| 婷婷成人精品国产| 欧美+日韩+精品| 热re99久久精品国产66热6| 欧美激情 高清一区二区三区| 欧美变态另类bdsm刘玥| 久久国产精品男人的天堂亚洲 | 精品一品国产午夜福利视频| videosex国产| av福利片在线| 五月天丁香电影| 看免费成人av毛片| 亚洲av福利一区| 亚洲五月色婷婷综合| 人妻一区二区av| 制服诱惑二区| 国产亚洲精品久久久com| 欧美激情 高清一区二区三区| 五月伊人婷婷丁香| 91精品伊人久久大香线蕉| 97超碰精品成人国产| 日本wwww免费看| 日本黄色片子视频| 夫妻性生交免费视频一级片| 亚洲精品乱久久久久久| 亚洲欧美一区二区三区国产| 亚洲在久久综合| 免费不卡的大黄色大毛片视频在线观看| 国产片特级美女逼逼视频| 黄色配什么色好看| 国产免费一区二区三区四区乱码| 日本-黄色视频高清免费观看| 大香蕉97超碰在线| 91久久精品电影网| 久久99热6这里只有精品| 成人手机av| 女的被弄到高潮叫床怎么办| 成人国产麻豆网| 久久精品国产鲁丝片午夜精品| 久久久久久久久大av| 乱码一卡2卡4卡精品| 亚洲综合色网址| 免费人成在线观看视频色| 国产精品三级大全| 99国产精品免费福利视频| av在线播放精品| 女性生殖器流出的白浆| 国语对白做爰xxxⅹ性视频网站| 大陆偷拍与自拍| 国产精品嫩草影院av在线观看| 我的老师免费观看完整版| 97精品久久久久久久久久精品| 美女大奶头黄色视频| 水蜜桃什么品种好| 久久久午夜欧美精品| 三级国产精品片| 国产 一区精品| 熟妇人妻不卡中文字幕| 桃花免费在线播放| 日韩欧美一区视频在线观看| 久久ye,这里只有精品| 国产探花极品一区二区| 丝袜美足系列| 国产免费福利视频在线观看| 国产成人aa在线观看| av在线观看视频网站免费| 中文字幕亚洲精品专区| 国产国语露脸激情在线看| 国产日韩欧美亚洲二区| 欧美国产精品一级二级三级| 国产精品不卡视频一区二区| 少妇猛男粗大的猛烈进出视频| 大陆偷拍与自拍| 久久久久网色| 男的添女的下面高潮视频| 亚洲人与动物交配视频| 国产一区二区在线观看日韩| 亚洲国产精品国产精品| 色哟哟·www| 免费看不卡的av| 视频中文字幕在线观看| 美女脱内裤让男人舔精品视频| 能在线免费看毛片的网站| 人成视频在线观看免费观看| 国产爽快片一区二区三区| 国产精品国产三级专区第一集| 国产亚洲一区二区精品| 人人妻人人澡人人看| 另类精品久久| 桃花免费在线播放| h视频一区二区三区| 综合色丁香网| 国产高清不卡午夜福利| 九色亚洲精品在线播放| 哪个播放器可以免费观看大片| 十八禁高潮呻吟视频| 国产精品免费大片| 国产毛片在线视频| 黄色毛片三级朝国网站| 精品国产露脸久久av麻豆| 99久国产av精品国产电影| 18在线观看网站| 成人二区视频| 极品人妻少妇av视频| 亚洲丝袜综合中文字幕| 午夜免费鲁丝| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美一区二区三区国产| 少妇丰满av| av在线老鸭窝| 久久久久人妻精品一区果冻| 日韩熟女老妇一区二区性免费视频| 婷婷色综合www| freevideosex欧美| 蜜桃久久精品国产亚洲av| 中国三级夫妇交换| 国产精品国产av在线观看| 国产成人av激情在线播放 | 最近中文字幕高清免费大全6| 精品酒店卫生间| 一级毛片我不卡| 久久久久久久久久久久大奶| 女的被弄到高潮叫床怎么办| 亚洲精品色激情综合| 亚洲精华国产精华液的使用体验| 午夜福利视频精品| 边亲边吃奶的免费视频| 久久精品久久久久久久性| 日日撸夜夜添| 久久这里有精品视频免费| 日本爱情动作片www.在线观看| 你懂的网址亚洲精品在线观看| 欧美精品亚洲一区二区| 又大又黄又爽视频免费| 久久久a久久爽久久v久久| 亚洲av电影在线观看一区二区三区| 久久久久国产网址| 大香蕉97超碰在线| 亚洲熟女精品中文字幕| tube8黄色片| 亚洲欧洲精品一区二区精品久久久 | 插逼视频在线观看| 亚洲欧美一区二区三区国产| 草草在线视频免费看| 国产成人精品一,二区| 丰满迷人的少妇在线观看| 久久久久国产网址| 国产精品久久久久久精品电影小说| av播播在线观看一区| 亚洲av在线观看美女高潮| 男女国产视频网站| 国产免费现黄频在线看| 少妇 在线观看| 精品久久久久久电影网| 在线观看免费视频网站a站| 啦啦啦视频在线资源免费观看| 寂寞人妻少妇视频99o| 在线观看www视频免费| 亚洲无线观看免费| 亚洲精品国产av成人精品| 国产精品国产三级专区第一集| 高清在线视频一区二区三区| 一本—道久久a久久精品蜜桃钙片| 久久毛片免费看一区二区三区| 国产高清有码在线观看视频| 寂寞人妻少妇视频99o| 丝袜在线中文字幕| 在线观看三级黄色| 欧美激情极品国产一区二区三区 | 99九九在线精品视频| 国产一区有黄有色的免费视频| 日韩熟女老妇一区二区性免费视频| 欧美性感艳星| 美女中出高潮动态图| 久久鲁丝午夜福利片| 亚州av有码| 日韩视频在线欧美| 亚洲av福利一区| 久久狼人影院| 日本wwww免费看| 国产精品麻豆人妻色哟哟久久| 成人国语在线视频| 精品亚洲成国产av| 我的女老师完整版在线观看| 精品一区二区三卡| 亚洲内射少妇av| 亚洲精品久久成人aⅴ小说 | 亚洲精品av麻豆狂野| 成人免费观看视频高清| 大片免费播放器 马上看| 亚洲av二区三区四区| 久久国产亚洲av麻豆专区| 在线天堂最新版资源| 国产亚洲精品久久久com| 亚洲国产精品一区二区三区在线| 国产片特级美女逼逼视频| 亚洲精品亚洲一区二区| 国产视频内射| 一级,二级,三级黄色视频| 蜜桃国产av成人99| 少妇被粗大猛烈的视频| 内地一区二区视频在线| 久久精品久久久久久噜噜老黄| 日本91视频免费播放| 黑人高潮一二区| 成人综合一区亚洲| a级毛片免费高清观看在线播放| 精品人妻熟女毛片av久久网站| 色婷婷久久久亚洲欧美| 欧美bdsm另类| 亚洲丝袜综合中文字幕| 高清不卡的av网站| 91午夜精品亚洲一区二区三区| 中文字幕亚洲精品专区| 日韩强制内射视频| 99视频精品全部免费 在线| 丝袜在线中文字幕| 午夜久久久在线观看| 人体艺术视频欧美日本| 亚洲精品第二区| 中国三级夫妇交换| 不卡视频在线观看欧美| 国产精品一区www在线观看| 九九久久精品国产亚洲av麻豆| 国产精品人妻久久久久久| 亚洲四区av| 国产欧美日韩综合在线一区二区| 亚洲精品一二三| 国产在视频线精品| 一区在线观看完整版| h视频一区二区三区| 好男人视频免费观看在线| 国产黄频视频在线观看| h视频一区二区三区| 亚洲婷婷狠狠爱综合网| 男人操女人黄网站| h视频一区二区三区| 日韩欧美精品免费久久| 在线观看免费日韩欧美大片 |