• <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手机在线免费观看| 亚洲av免费高清在线观看| 啦啦啦在线观看免费高清www| av网站免费在线观看视频| 性高湖久久久久久久久免费观看| 亚洲欧美日韩东京热| 日韩一区二区视频免费看| 秋霞在线观看毛片| 女性被躁到高潮视频| 亚洲va在线va天堂va国产| 久久久久久久国产电影| 国产精品一区二区三区四区免费观看| 亚洲真实伦在线观看| 18禁裸乳无遮挡动漫免费视频| 日韩成人av中文字幕在线观看| 欧美xxxx性猛交bbbb| 观看免费一级毛片| 大陆偷拍与自拍| 极品人妻少妇av视频| 日产精品乱码卡一卡2卡三| 欧美高清成人免费视频www| 亚洲av在线观看美女高潮| 99热这里只有精品一区| 香蕉精品网在线| 中文乱码字字幕精品一区二区三区| 日本黄色片子视频| 亚洲国产精品999| 成人毛片a级毛片在线播放| 亚洲国产av新网站| 精品久久久精品久久久| 久久亚洲国产成人精品v| 熟女电影av网| 国产欧美另类精品又又久久亚洲欧美| 2021少妇久久久久久久久久久| av在线app专区| 国产精品一区二区三区四区免费观看| 日韩制服骚丝袜av| 大话2 男鬼变身卡| 99久久综合免费| 久久影院123| 九色成人免费人妻av| 国产成人免费观看mmmm| 菩萨蛮人人尽说江南好唐韦庄| 国产成人精品福利久久| h日本视频在线播放| 国产熟女欧美一区二区| 国产成人一区二区在线| 成人午夜精彩视频在线观看| 少妇被粗大的猛进出69影院 | 精品人妻一区二区三区麻豆| 久热久热在线精品观看| 日日啪夜夜撸| 成人亚洲精品一区在线观看| 内射极品少妇av片p| 亚洲,欧美,日韩| 欧美最新免费一区二区三区| 国产女主播在线喷水免费视频网站| 亚洲一区二区三区欧美精品| 人人妻人人看人人澡| 久久久久久久久久成人| 亚洲国产精品一区二区三区在线| 九九在线视频观看精品| 国产伦在线观看视频一区| 欧美一级a爱片免费观看看| 一级爰片在线观看| 男人添女人高潮全过程视频| 久久久国产精品麻豆| 边亲边吃奶的免费视频| 国产欧美亚洲国产| 中文字幕精品免费在线观看视频 | 亚洲婷婷狠狠爱综合网| .国产精品久久| 少妇人妻久久综合中文| 亚洲怡红院男人天堂| 国产成人91sexporn| 日产精品乱码卡一卡2卡三| 国产成人精品福利久久| 视频区图区小说| 午夜福利在线观看免费完整高清在| 欧美激情国产日韩精品一区| 亚洲av在线观看美女高潮| 亚洲va在线va天堂va国产| 中国三级夫妇交换| 亚洲精品国产av蜜桃| 在线亚洲精品国产二区图片欧美 | 交换朋友夫妻互换小说| 欧美区成人在线视频| 日韩av不卡免费在线播放| 国产一区二区三区av在线| 在线观看国产h片| 精品人妻熟女av久视频| 亚洲人成网站在线观看播放| 亚洲熟女精品中文字幕| 能在线免费看毛片的网站| 97在线人人人人妻| 成人毛片60女人毛片免费| 国产亚洲91精品色在线| 欧美日韩av久久| 国产高清三级在线| 一级二级三级毛片免费看| 亚洲图色成人| 欧美成人午夜免费资源| 国产视频内射| 我要看黄色一级片免费的| 男女国产视频网站| av又黄又爽大尺度在线免费看| 久久久久精品久久久久真实原创| 日韩制服骚丝袜av| 这个男人来自地球电影免费观看 | 久久久久人妻精品一区果冻| 偷拍熟女少妇极品色| 日韩欧美精品免费久久| 久久婷婷青草| av黄色大香蕉| 免费黄频网站在线观看国产| 国产精品久久久久久久久免| 熟妇人妻不卡中文字幕| av不卡在线播放| 香蕉精品网在线| 午夜免费观看性视频| 久久久久久久久久久免费av| 久久精品国产亚洲av涩爱| 久久久国产精品麻豆| 亚洲av日韩在线播放| 成年人免费黄色播放视频 | xxx大片免费视频| 91精品国产九色| 日日撸夜夜添| 久热这里只有精品99| 久热这里只有精品99| 女人久久www免费人成看片| 国产成人精品福利久久| 99国产精品免费福利视频| 日韩在线高清观看一区二区三区| 国产av精品麻豆| 久久免费观看电影| 欧美精品人与动牲交sv欧美| 少妇精品久久久久久久| 91成人精品电影| 精品人妻熟女毛片av久久网站| 国产永久视频网站| 日韩一区二区三区影片| 国产色爽女视频免费观看| 亚洲国产色片| 欧美国产精品一级二级三级 | 又黄又爽又刺激的免费视频.| 亚洲欧美日韩东京热| 永久免费av网站大全| 国产女主播在线喷水免费视频网站| av在线观看视频网站免费| 99久久综合免费| 三级国产精品片| 国产91av在线免费观看| 国产淫语在线视频| 国产精品一区二区性色av| 亚洲国产精品专区欧美| 久久国产亚洲av麻豆专区| 亚洲精品456在线播放app| 免费人成在线观看视频色| .国产精品久久| 五月伊人婷婷丁香| 久久ye,这里只有精品| 少妇的逼好多水| 亚洲综合精品二区| 天堂俺去俺来也www色官网| 精品亚洲成a人片在线观看| 一个人看视频在线观看www免费| 国产男女超爽视频在线观看| 久久精品国产自在天天线| 有码 亚洲区| 男女国产视频网站| 日韩一区二区视频免费看| 中文资源天堂在线| 国产精品国产av在线观看| 精品人妻熟女毛片av久久网站| 国产高清不卡午夜福利| 少妇高潮的动态图| 交换朋友夫妻互换小说| 嘟嘟电影网在线观看| 中文字幕亚洲精品专区| 欧美日韩视频精品一区| 亚洲精品456在线播放app| 蜜桃在线观看..| av在线观看视频网站免费| 日韩成人av中文字幕在线观看| 晚上一个人看的免费电影| 久久久久久久久久久免费av| 国产av一区二区精品久久| kizo精华| 哪个播放器可以免费观看大片| 欧美另类一区| 国产精品99久久99久久久不卡 | 最新的欧美精品一区二区| 草草在线视频免费看| 少妇裸体淫交视频免费看高清| 午夜老司机福利剧场| 天堂俺去俺来也www色官网| 亚洲在久久综合| 欧美人与善性xxx| 亚洲综合色惰| 国产无遮挡羞羞视频在线观看| av有码第一页| 精品久久国产蜜桃| 国产69精品久久久久777片| .国产精品久久| 夫妻午夜视频| 日日啪夜夜爽| 男女边摸边吃奶| 一区二区三区四区激情视频| 亚洲av欧美aⅴ国产| 搡老乐熟女国产| 蜜桃久久精品国产亚洲av| 99re6热这里在线精品视频| 国产高清不卡午夜福利| 黄色欧美视频在线观看| 国产一区二区三区综合在线观看 | 久久久久人妻精品一区果冻| 熟女人妻精品中文字幕| 亚洲三级黄色毛片| 亚洲av中文av极速乱| av一本久久久久| 22中文网久久字幕| 18禁裸乳无遮挡动漫免费视频| 精品久久国产蜜桃| 久久久久网色| 99热这里只有是精品50| 国产有黄有色有爽视频| 一级片'在线观看视频| 男人狂女人下面高潮的视频| 日韩av不卡免费在线播放| 男男h啪啪无遮挡| 99re6热这里在线精品视频| 色吧在线观看| 中文字幕亚洲精品专区| 91久久精品电影网| 久久精品久久精品一区二区三区| av专区在线播放| 五月伊人婷婷丁香| 99热全是精品| 色网站视频免费| a级毛片免费高清观看在线播放| 亚洲国产成人一精品久久久| 91午夜精品亚洲一区二区三区| 亚洲欧美日韩卡通动漫| 国产 精品1| 草草在线视频免费看| 噜噜噜噜噜久久久久久91| 又大又黄又爽视频免费| 久久久久久久久久成人| 亚洲真实伦在线观看| 成人18禁高潮啪啪吃奶动态图 | 精品国产一区二区久久| 国产中年淑女户外野战色| 大陆偷拍与自拍| 9色porny在线观看| 少妇被粗大的猛进出69影院 | 国产 一区精品| 一区二区三区乱码不卡18| 秋霞伦理黄片| 久久免费观看电影| 日韩三级伦理在线观看| 中国国产av一级| 曰老女人黄片| 免费黄频网站在线观看国产| 亚洲精华国产精华液的使用体验| 久久热精品热| 久久ye,这里只有精品| 国产亚洲一区二区精品| 亚洲,欧美,日韩| 欧美bdsm另类| 嫩草影院入口| 在线播放无遮挡| 亚洲国产日韩一区二区| 十八禁网站网址无遮挡 | 黄色欧美视频在线观看| 插逼视频在线观看| 97超视频在线观看视频| 美女cb高潮喷水在线观看| 精品视频人人做人人爽| 交换朋友夫妻互换小说| 天美传媒精品一区二区| 国产亚洲av片在线观看秒播厂| 最新的欧美精品一区二区| 视频中文字幕在线观看| 久久久久久久国产电影| 色94色欧美一区二区| 欧美3d第一页| 国产精品国产av在线观看| 中文欧美无线码| 国产成人a∨麻豆精品| 亚洲精品乱码久久久v下载方式| 少妇的逼水好多| 亚洲国产精品一区二区三区在线| 国产黄片视频在线免费观看| 深夜a级毛片| 九九久久精品国产亚洲av麻豆| 久久 成人 亚洲| 中文字幕免费在线视频6| kizo精华| 亚洲国产精品专区欧美| 丰满少妇做爰视频| 国产免费又黄又爽又色| 在线观看av片永久免费下载| 久久97久久精品| videos熟女内射| 极品少妇高潮喷水抽搐| 在线看a的网站| 免费久久久久久久精品成人欧美视频 | 久久久国产欧美日韩av| 91午夜精品亚洲一区二区三区| 国产视频首页在线观看| 亚洲欧美日韩卡通动漫| 久久ye,这里只有精品| 秋霞在线观看毛片| 大香蕉久久网| 亚洲一级一片aⅴ在线观看| 纵有疾风起免费观看全集完整版| 街头女战士在线观看网站| 亚洲欧美精品专区久久| 七月丁香在线播放| 日本欧美国产在线视频| 麻豆成人av视频| 日韩av免费高清视频| 成人无遮挡网站| 少妇被粗大的猛进出69影院 | 国产真实伦视频高清在线观看| 嫩草影院入口| 国产av一区二区精品久久| 看非洲黑人一级黄片| 老司机影院成人| 国产免费又黄又爽又色| 国产精品国产三级专区第一集| 亚洲,一卡二卡三卡| 少妇精品久久久久久久| 亚洲av国产av综合av卡| 色视频在线一区二区三区| 久久久久久久久久久丰满| 欧美精品一区二区大全| 中文欧美无线码| 午夜免费鲁丝| 免费黄频网站在线观看国产| 国产极品天堂在线| 毛片一级片免费看久久久久| 大又大粗又爽又黄少妇毛片口| 欧美日韩在线观看h| 午夜福利影视在线免费观看| 97超碰精品成人国产| 精品久久久久久电影网| 国产精品一区www在线观看| 精品久久久久久久久亚洲| 国模一区二区三区四区视频| tube8黄色片| 亚洲天堂av无毛| 在现免费观看毛片| 亚洲中文av在线| 日韩不卡一区二区三区视频在线| 国产片特级美女逼逼视频| 十八禁网站网址无遮挡 | 成年av动漫网址| 成人毛片60女人毛片免费| 亚洲色图综合在线观看| 亚洲av成人精品一区久久| 国产黄频视频在线观看| 精品国产国语对白av| 国产成人午夜福利电影在线观看| 黄色毛片三级朝国网站 | 黄片无遮挡物在线观看| 国产在视频线精品| 中国美白少妇内射xxxbb| 亚洲综合色惰| 久久99热6这里只有精品| 欧美日韩综合久久久久久| 亚洲av电影在线观看一区二区三区| 欧美成人午夜免费资源| 伊人亚洲综合成人网| 免费观看无遮挡的男女| 秋霞伦理黄片| 一区二区av电影网| 最近的中文字幕免费完整| 国产在线免费精品| 国产精品人妻久久久影院| 国产69精品久久久久777片| 精品久久久精品久久久| 午夜福利网站1000一区二区三区| 国产男人的电影天堂91| 欧美丝袜亚洲另类| 国产男女超爽视频在线观看| 黄色日韩在线| a级毛片在线看网站| 久久青草综合色| 免费黄网站久久成人精品| 在线观看av片永久免费下载| 久久国产亚洲av麻豆专区| 涩涩av久久男人的天堂| 久久久久久久大尺度免费视频| 国产精品国产三级国产专区5o| 亚洲成人手机| 丁香六月天网| 亚洲真实伦在线观看| 春色校园在线视频观看| 国产伦理片在线播放av一区| 亚洲一级一片aⅴ在线观看| 秋霞在线观看毛片| 中文资源天堂在线| 日韩不卡一区二区三区视频在线| 国产亚洲一区二区精品| 亚洲精品aⅴ在线观看| videossex国产| 国产伦在线观看视频一区| 九色成人免费人妻av| 久久精品国产a三级三级三级| 久久久亚洲精品成人影院| 午夜老司机福利剧场| 一个人免费看片子| 极品教师在线视频| 国产精品一区二区在线观看99| 国产视频首页在线观看| 国产日韩欧美亚洲二区| 精品人妻偷拍中文字幕| 夜夜爽夜夜爽视频| 免费看日本二区| 日韩,欧美,国产一区二区三区| 777米奇影视久久| 国产精品蜜桃在线观看| 国产精品一区二区性色av| 亚洲,一卡二卡三卡| 99热全是精品| 亚洲精品一区蜜桃| 婷婷色麻豆天堂久久| 一级毛片我不卡| 精品酒店卫生间| 亚洲精品国产成人久久av| 精品一区二区三卡| 亚洲欧洲国产日韩| 韩国高清视频一区二区三区| 日本黄大片高清| 性色av一级| 亚洲国产欧美日韩在线播放 | 国产欧美日韩精品一区二区| 这个男人来自地球电影免费观看 | 日日啪夜夜撸| 夜夜骑夜夜射夜夜干| 亚洲久久久国产精品| 午夜激情久久久久久久| 乱人伦中国视频| 婷婷色综合大香蕉| 日韩一区二区三区影片| 永久免费av网站大全| 亚洲国产最新在线播放| 女的被弄到高潮叫床怎么办| 国产日韩欧美在线精品| 99久久精品热视频| av在线老鸭窝| 亚洲精品国产色婷婷电影| 麻豆成人av视频| 中文字幕av电影在线播放| 人体艺术视频欧美日本| 国产成人aa在线观看| 欧美精品一区二区免费开放| 内地一区二区视频在线| 欧美精品国产亚洲| 午夜福利,免费看| 嫩草影院新地址| 三级经典国产精品| 日韩精品免费视频一区二区三区 | 精品酒店卫生间| 精品久久久精品久久久| 国产成人freesex在线| 在线观看免费高清a一片| 免费观看在线日韩| 免费在线观看成人毛片| 日韩精品有码人妻一区| 精品一区二区三卡| 少妇裸体淫交视频免费看高清| 色94色欧美一区二区| 性色av一级| 亚洲精品日韩在线中文字幕| 日韩亚洲欧美综合| 热re99久久精品国产66热6| 日日摸夜夜添夜夜添av毛片| 丰满人妻一区二区三区视频av| 韩国高清视频一区二区三区| 国产有黄有色有爽视频| 国产av国产精品国产| 亚洲国产精品999| 国产国拍精品亚洲av在线观看| 欧美成人精品欧美一级黄| 女性生殖器流出的白浆| 久久久久久久精品精品| 日韩在线高清观看一区二区三区| 欧美少妇被猛烈插入视频| 日韩电影二区| 亚洲自偷自拍三级| 国产精品麻豆人妻色哟哟久久| 国产成人免费观看mmmm| 久久精品国产a三级三级三级| 日日啪夜夜撸| 七月丁香在线播放| 亚洲精品国产av蜜桃| 一区在线观看完整版| 婷婷色av中文字幕| 国产精品偷伦视频观看了| 我要看黄色一级片免费的| 少妇高潮的动态图| 91在线精品国自产拍蜜月| 91精品一卡2卡3卡4卡| 丰满乱子伦码专区| 亚洲精品,欧美精品| 国产亚洲欧美精品永久| 久久久久久人妻| 搡女人真爽免费视频火全软件| 久久久久久久大尺度免费视频| 久久国产乱子免费精品| 麻豆乱淫一区二区| 少妇人妻一区二区三区视频| 一级av片app| av免费观看日本| 自线自在国产av| 女人精品久久久久毛片| 久久久a久久爽久久v久久| 我的女老师完整版在线观看| 国产一区有黄有色的免费视频| 国产亚洲91精品色在线| 久热这里只有精品99| 在线观看国产h片| 亚洲av成人精品一区久久| 一级a做视频免费观看| 中国美白少妇内射xxxbb| 高清不卡的av网站| 黑丝袜美女国产一区| 美女cb高潮喷水在线观看| 夜夜看夜夜爽夜夜摸| 99久久精品国产国产毛片| 国产女主播在线喷水免费视频网站| 国产美女午夜福利| 日韩欧美 国产精品| 欧美精品亚洲一区二区| 日本与韩国留学比较| 乱码一卡2卡4卡精品| 国产一区二区在线观看av| 久久精品国产a三级三级三级| 视频区图区小说| 欧美精品高潮呻吟av久久| 国产精品人妻久久久影院| 国产色婷婷99| 免费久久久久久久精品成人欧美视频 | 亚洲激情五月婷婷啪啪| 曰老女人黄片| 国产女主播在线喷水免费视频网站| 国产精品久久久久成人av| 国产有黄有色有爽视频| 欧美 亚洲 国产 日韩一| 色94色欧美一区二区| 美女大奶头黄色视频| 99九九在线精品视频 | 日韩不卡一区二区三区视频在线| 在线亚洲精品国产二区图片欧美 | 最新的欧美精品一区二区| 人妻少妇偷人精品九色| 成年美女黄网站色视频大全免费 | 黑人巨大精品欧美一区二区蜜桃 | 十分钟在线观看高清视频www | 高清午夜精品一区二区三区| 国产精品久久久久成人av| 精品国产露脸久久av麻豆| 精品一区二区三区视频在线| 日韩av不卡免费在线播放| 最近最新中文字幕免费大全7| 日韩视频在线欧美| 一级毛片电影观看| 久久亚洲国产成人精品v| 中国国产av一级| 久久久国产精品麻豆| 午夜福利影视在线免费观看| 纵有疾风起免费观看全集完整版| 在线免费观看不下载黄p国产| 中文欧美无线码| 在线亚洲精品国产二区图片欧美 | 国产亚洲91精品色在线| av视频免费观看在线观看| 伦理电影免费视频| 国产精品不卡视频一区二区| 日本av手机在线免费观看| 欧美成人午夜免费资源| 亚洲欧美成人综合另类久久久| 精品亚洲乱码少妇综合久久| 久久人妻熟女aⅴ| 亚洲国产欧美在线一区| 自线自在国产av| 少妇人妻 视频| 高清在线视频一区二区三区| 午夜激情久久久久久久| 亚洲国产最新在线播放| 一级毛片黄色毛片免费观看视频| 青青草视频在线视频观看| 少妇人妻精品综合一区二区| 国产黄色免费在线视频| 免费看日本二区| 人人妻人人澡人人看| videos熟女内射| 国产亚洲欧美精品永久| 在线播放无遮挡| 91久久精品国产一区二区成人| 精品国产露脸久久av麻豆| 精品人妻熟女av久视频| 亚洲精品,欧美精品| 久久久久久久久久久丰满| 一级片'在线观看视频| 国产精品国产av在线观看| 日本av手机在线免费观看|