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

    Fault diagnosis of an intelligent hydraulic pump based on a nonlinear unknown input observer

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

    Zhonghai MA,Shaoping WANG,Jian SHI,Tongyang LI,Xingjian WANG

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

    1.Introduction

    A hydraulic pump is the key component in an aircraft hydraulic system,on whose operating performance the quality of the pump has an important impact.In order to meet those requirements of modern aircraft,such as high speed,high mobility,low weight,and high load capacity,the intelligent hydraulic pump system(IHPS)has been paid more and more attention to for improving the efficiency of an aircraft hydraulic system.An aircraft hydraulic system with an IHPS has become an important future development direction,1–3and the requirements for reliability and security also grow significantly,4which implies that an efficient fault diagnosis will improve the reliability of the IHPS significantly.

    An IHPS is generally defined as a kind of pump source system whose output can be easily controlled by the virtue of an intelligent controller to meet those requirements of an actual aircraft hydraulic system.To realize feedback control of the output,necessary sensors are installed on an IHPS,including pressure sensors,displacement sensors,and temperature sensors.5Based on an IHPS’s state parameters and its actual operating condition,the controller adjusts the displacement of the pump in accordance with the pressure signal,and then achieves optimal matching with the load.The advantage of an IHPS is that,under the conditions of different aircraft flight statuses and commands from the flight control computer,the output of the IHPS can match the load reasonably.Thereby,an IHPS can reduce the invalid power and heat of an aircraft hydraulic system to improve the aircraft’s work efficiency.6Up to now,a lot of research on the structure design of an IHPS has been done.7–9Ma et al.proposed a new type of axial piston variable displacement pump,7which was an improvement on the basis of an A4V pump and used a servo valve to control the variable mechanism directly.This improved pump can sense the load correctly and timely based on the output pressure command of the intelligent pump controller,and can provide necessary information for monitoring the pump operating condition and fault diagnosis.8However,this intelligent pump is more complex than a traditional hydraulic pump in structure.Comparing to a traditional hydraulic pump,an IHPS installs many additional sensors to attain pump state signals.10Furthermore,an IHPS integrates a servo valve to control the pump variable mechanism,so as to realize servo control.9The complex structure has made the structure design and signal processing of the IHPS a tough job,11,12and the failure probability of the IHPS also increases with a dynamic change of the operating condition.Some previous studies have been mainly focused on the fault diagnosis of a common piston pump.Lu et al.presented a fault diagnosis method of piston pumps using a pump discharge pressure signal.13Wang and Hu proposed a vibration-signals-based fault diagnosis method using the fuzzy technique.14Some artificial neural network(ANN)models have also been widely applied to pump fault diagnosis considering their nonlinearity,adaptability,and robustness.15Lu et al.presented a diagnosis method for hydraulic pumps based on the chaotic parallel radial basis function(RBF)network.16

    These above methods have been proven to be effective for piston-pump fault diagnosis to a certain extent.However,many limitations still exist obviously in these methods,which becomes pretty severe especially when we apply these methods to IHPS fault diagnose.Due to structural complexity and flexible operation conditions,the mechanism of IHPS failure is more complicated,and the fault detection and diagnosis for an IHPS is more difficult.Comparing to common variable pumps,the working modes of an IHPS switch more frequently.According to a real aircraft flight profile,an IHPS would switch its operation within four kinds of different working modes,that is,constant pressure mode,constant flow mode,constant power mode,and load sensitive mode.17The outputs of flow and pressure in each mode are totally different,which will obviously increase the fatigue and aging of the reset spring,as well as the wear of the variable cylinder of the pump.In such situations,the swashplate angle cannot be adjusted with a command signal,and a flow regulation of the IHPS cannot be performed.Furthermore,the leakage due to the wear of three friction pairs(swashplate and slipper,valve plate and cylinder block,and piston and cylinder bore)becomes more complex,in which the aging of sealing elements and the failure of the connection between the shaft and the roller bearing should be considered.18When it comes to the two aforementioned common faults of an IHPS,traditional fault diagnosis methods are not feasible for the analysis of intelligent pumps.8,19On one hand,due to long life cycles and actual complex flight conditions,it is difficult to obtain a large quantity of faulty pump samples and fault data.Therefore,the expert knowledge of pump failure is not perfect,and it brings difficulties to an accurate fault diagnosis based on artificial intelligence which needs long-term practical experience and a large quantity of fault information.20–22On the other hand,because of its severe working environment,the slowly varying fault value of leakage and aging may be overwhelmed by environmental noise and cannot be separated from the environmental interference.The dynamic change of various conditions leads to complex fault characteristic signals,which are difficult to analyze with a signal-based diagnosis method using a fixed threshold.23–25Therefore,it is urgent to find an effective fault diagnosis method which is suitable for the complex technical characteristics of an IHPS to ensure its reliable operation.

    In this paper,to achieve higher reliability and security and improve the performance of an IHPS,a fault diagnosis method based on a nonlinear unknown input observer(NUIO)is proposed to realize IHPS fault detection.Different from factors of a full-order Luenberger-type unknown input observer presented by Chen et al.and Duan et al.26,27nonlinear factors of the IHPS are considered.By comparing the measured signals of the object to be diagnosed and the prior information of the system expressed by a model,a residual error is generated by utilizing a NUIO,and then the residual error is analyzed and processed to realize fault diagnosis.28Based on analysis of IHPS working mechanism and typical failure modes,an accurate nonlinear mathematical model of the IHPS is established.Then,a dynamic fault diagnosis method based on dynamic characteristics and variability of the signals is proposed.29,30Finally,based on the MATLAB simulation platform,the simulation analysis of two typical failure modes of the IHPS is carried out,and the results verify the validity and accuracy of this fault diagnosis method.

    The remaining part of this paper is organized as follows.In Section 2,a mathematical model of an IHPS is established considering the working mechanism and failure modes.Section 3 presents the IHPS NUIO on the basis of output pressure and swashplate angle signals.Section 4 carries out model based fault diagnose based on the NUIO for typical failure modes of the intelligent pump according to detailed parameters.Section 5 gives conclusions.

    2.Mechanism analysis of an intelligent pump

    2.1.Structure and operational mechanism of an intelligent pump

    A conventional axial piston pump,as shown in Fig.1,is widely used to provide constant high-pressure oil to manipulate an outside load,which leads to power dissipation under a smallload demand.Therefore,an intelligent pump appears,which consists of a reciprocating cylinder swashplate,an embedded control circuit,and a servo variable mechanism.The variable displacement pump of the IHPS is controlled by a microcomputer through a feedback signal.The control circuit accepts the command from the flight computer according to the flight profile and controls the angle of the swashplate to get appropriate flow.This kind of design can make the intelligent pump provide a flow dynamically according to the load demand in order to save energy.In addition,the intelligent pump can also give an appropriate failure treatment when a failure happens.The core of the IHPS is a variable displacement pump controlled by a microcomputer,which makes it adaptable to various load requirements,capable of dealing with fault states,and have good transient response ability.This adaptive adjustment of the swashplate overcomes the throttling loss of a traditional axial piston pump under constant pressure,and reduces the invalid power and heat of an aircraft hydraulic system.31In order to figure out the dynamic characteristics of the intelligent pump,it is necessary to establish its mathematical model and analyze its dynamic performance.

    The theoretical flow of the intelligent pumpQtis related to the area of a piston,the number of pistons,the swashplate angle,and the speed of the pump as follows:

    Fig.1 A certain type of intelligent pump.

    wheredis the diameter of a piston,Zis the number of pistons,nis the speed of the pump,γ is the angle of the swashplate of the pump,andris the radius of the piston distribution circle.The theoretical flowQtis nonlinear with the angle of the swashplate γ.

    Based on Eq.(1),the theoretical angular displacementqt(γ)with respect to the angle of the swashplate γ can be expressed as

    In order to adjust the output flow of the intelligent pump dynamically,the variable mechanism controlled by the embedded microcomputer is shown in Fig.2.As shown in Fig.2,the motor provides a torqueMto drive the hydraulic pump,and the output pressure of the systemPsand the load pressure of the systemPLcan be measured by pressure sensors respectively,while the temperatures of the system are also obtained by three temperature sensors.The displacement of the swashplate is measured by a LVDT(linear variable differential transformer)sensor to calculate the swashplate angle γ.The temperatures of the output,return oil,and shell of the pump are measured by temperature sensors.Its detailed modeling and analysis is the basis and precondition of further research on the pump,even the IHPS.

    2.1.1.Flow equation of swashplate mechanism

    Assuming that the pressure of the pump source with the variable mechanismPVSis a constant,the swashplate variable mechanism could be regarded as an independent hydraulic cylinder system controlled by a load valve.32When the spool displacementxs>0,the flow equation of the servo valve can be described as

    whereQVLis the load flow of the rod cavity within the variable mechanism,andQVL1is the load flow of the non-rod cavity within the variable mechanism.KVQandKVQ1are the flow amplification coefficients in the variable mechanism of the swashplate,which can be described respectively as follows:

    When the spool displacementxs<0,the flow equation of the servo valve is shown as

    2.1.2.Flow equation of variable cylinder

    whereAVis the piston area of the variable cylinder,lis the distance between the axis of the cylinder and the rotation center of the swashplate,VVis the total volume of the variable cylinder,Eyis the elastic modulus of the equivalent volume of hydraulic oil,CVLis the total leakage coefficient of the variable cylinder,andCVipandCVepare the internal and external leakage coefficients of the variable cylinder,respectively.

    Fig.2 Variable mechanism of intelligent pump.

    2.1.3.Torque equation of variable cylinder

    whereMLis the output torque of the variable cylinder,MVIis the total inertia torque of the swashplate variable mechanism and the load on the rotary shaft,MVBis the total viscous damping torque of the swashplate variable mechanism,MVKis the elastic torque of the spring,andMLfis the external disturbance torque of the swashplate variable mechanism.Detailed calculations of the above torques can be described as follows:

    whereMSIandMCIare the inertia torques of the swashplate and the variable cylinder with respect to the rotating shaft of the swashplate,respectively,ISandICare the inertia of the swashplate and the variable cylinder with respect to the rotating shaft of the swashplate,MSBis the viscous damping torque of the swashplate,MCBis the viscous damping torque of the variable cylinder,BSandBCare the viscous damping coefficients of the swashplate and the variable cylinder,respectively,Ris the equivalent radius of the swashplate bearing,andKVis the load elastic stiffness.

    Based on the analysis and calculation of each torque,the dynamics of hydraulic servo systems are strongly nonlinear and uncertain,due to fluid compressibility,friction,and aging.Therefore,Eq.(11)can be simplified as follows:

    whereIVis the total inertia of the swashplate variable mechanism,BVis the total viscous damping coefficient,Vpis the pressure coefficient of the disturbance torque,andMfis some uncertain torque caused by friction and other disturbances which are nonlinear.

    The spool displacementxscan be expressed as follows:

    whereKcis the gain of the controller and the sensor,which converts a pressure signal into a displacement signal that can be the output of the servo valve.Puis the control input signal of pressure calculated by the IHPS.

    2.1.4.Mathematical model of intelligent pump

    The relationship betweenQtand the output pressure of the pumpPscan be shown as

    where|QL|is the load flow of the pump,Cplis its total pressure leakage coefficient,andVsis the volume of the output.

    According to Eq.(2),(4)and(5),the theoretical flowQtis nonlinear with the angle of the swashplate γ,QVLis nonlinear with the statePVL,which also includes the nonlinear termMfshown in Eq.(18).The nonlinear terms in the system make it more complex,so the traditional transfer function method cannot be used to solve multivariable systems and detect failure situations from inner states.In order to take into account nonlinear factors and observe internal states,the state space method is used to describe the IHPS comprehensively,in which inner states are obtained via an observer.

    According to Eqs.(1),(3),(9),(18),and(22),its state space can be expressed as follows:

    where x(t)∈ R4×1is the state vector,y(t)∈ R2×1is the vector of measurable signals,υ is the vector of measurable signals,u is the control input vector which is assumed to be bounded,and durepresents the unknown input vector.A is the system matrix,B the control matrix,C the output matrix,and Bdthe output matrix,which are matrices with appropriate dimensions.Without losing generality,matrix Bdis assumed to be of full-column rank,whileg(υ,y,u)is the control input function.In the IHPS,the swashplate angle γ and the actual output pressurePscan be obtained by the LVDT and pressure sensors installed on the pump,respectively.An observer can be designed based on the input and signals measured by sensors.

    The matrixes with nonlinear parameters,such asg(υ,y,u)andqt(γ),can be obtained as follows:

    In order to carry out fault diagnosis of the intelligent pump with a high precision,it is necessary to observe its inner states which are sensitive to typical failures without any knowledge of unknown inputs.

    2.2.Main failure modes and influencing factors of IHPS

    Previous research has shown that the three main failure modes of an axial piston pump are embodied in leakage,fatigue damage,and aging.18The most important failure mode is the wear of three friction pairs(swashplate and slipper,valve plate and cylinder block,and piston and cylinder bore)33and the bearing and servo of the variable mechanism.When the wear and tear of friction pairs occur time to time,the pump will fail eventually because of leakage.Most of the components in the pump perform rotation and reciprocating motions,and under alternating stress,their fatigue is unnegligible,24such as the fatigue of the spring of the variable mechanism and stress raiser.Aging associated with sealing elements also increases the leakage of the pump.Based on the embedded controller,the failure mode of the IHPS also consists of short circuit,open circuit,and parameter shift.

    According to the analysis above,the most important failure effect is the leakage of the intelligent pump.Herein,Cpldenotes the total leakage coefficient of the intelligent pump.The leakage coefficient increases with the failure development of the intelligent pump under cyclic applied loads.KVindicates the elastic stiffness of the variable mechanism of the swashplate,which gradually decreases with an increase of the fault degree of the variable mechanism.

    3.Nonlinear unknown input observer(NUIO)design

    Due to dynamic operational conditions,the signal-based diagnosis method can only be a rough judgment of a failure and some sensors are not sensitive to early failures.The modelbased diagnosis method can be adopted to analyze the change of a dynamic system using deep knowledge.In addition,according to the analysis of the IHPS,nonlinear factors such as the control input functiong(υ,y,u)and the theoretical angular displacementqt(γ)can be obtained through establishing a nonlinear observer.Hence,this paper uses a NUIO to realize the fault diagnosis of the IHPS.Generally,fault diagnosis based on a mathematical model includes two stages:residual error generation and residual error evaluation.There are three methods to estimate a residual error,viz.parameter estimation method,state estimation method,and equivalent space method.A residual error is generated by comparing the measured information of an object to be diagnosed with the prior information of the system expressed by a model.If the residual error is greater than a threshold,the system is considered to be faulty.Since the intelligent pump is very close to an engine,its input or interference might be mixed with unknown factors such as noise and uncertainty.However,the distribution and bound of its input could be specified according to history and prior experience,which leads to the use of the NUIO.Based on the information given by the distribution matrix,the unknown input(interference)can be extracted from the residual error.Because the input and output information required for model-based fault diagnosis is only related to an open-loop system,it is not necessary to consider a controller in the design of a fault diagnosis scheme,which is consistent with the principle of separation of control theory.As long as the input of a regulator is used,it is the same for the fault diagnosis no matter whether the system is in closed or open loop.In order to simplify the model description,the uncertainty of the IHPS is summarized as the additive disturbance term in dynamic nonlinear Eq.(23)in this paper.In order to describe the intelligent pump comprehensively,the mathematical model can be described as follows:

    Therefore,the corresponding observer for the IHPS can be designed as

    Assuming that the estimation error is e(t)=^x(t)-x(t),its differential equation can be described as

    then e(t)goes to zero asymptotically and it is invariant with respect to the unknown input du(t).The notation I indicates the identity matrix.

    Eqs.(33)and(34)are equivalent to Eq.(35)as follows:

    All possible solutions of Eq.(35)will be

    Fig.3 Nonlinear unknown input observer(NUIO)structure.

    where Y is an arbitrary compatible matrix and X+= (XTX)-1XT.

    The state estimation error will be

    If the eigenvalues of Njare stable,then e(t)will approach zero asymptotically,i.e.,^x(t)→ x(t),which means that the observer in Eq.(28)is a NUIO for the IHPS according to the definition of unknown input observer.The design of this NUIO involves solving Eqs.(30)–(34)and guaranteeing that all eigenvalues of the system matrix Njare stable.The necessary and sufficient condition of the existence of a NUIO will be shown as follows26:

    According to the analysis of the NUIO,the key step is to design matrix K(ρ).Once it is available,other parameters of the NUIO can be obtained from Eqs.(30)–(34).The only constraint on K(ρ)is to make the dynamic characteristic matrix N(ρ)stable.Therefore,in the process of calculation,it is necessary to find a suitable value,which leads the state estimation error to zero asymptotically.The flowchart of the process within a specific NUIO is shown in Fig.4.

    In the design of a NUIO,the most important part is that,when matrix CBdis full-column rank,matrix N(ρ)could be kept stable by selecting a proper matrix K(ρ).In the IHPS,based on the characteristics of the dynamic system,the poles of matrix N(ρ)are configured to ensure the stability of the observer.

    4.Fault diagnosis of IHPS based on NUIO

    Based on the NUIO,a simulation model of the IHPS is established by using the Simulink toolbox,and the fault diagnosis system with a residual error is proposed.The simulation parameters of the IHPS are shown in Table 1.

    The concrete parameters of the IHPS can be obtained through calculation,and the values of N(ρ),E,L(ρ),and M are as follows:

    Fig.4 Flowchart of process of a nonlinear unknown input observer.

    where the eigenvalues of each matrix N(ρ)are negative,which satisfies the requirement of the observer system.

    Based on the MATLAB platform,a simulation model of the intelligent pump is established,and an observer is designed according to the method in Section 3.Simulation will be carried out for two kinds of typical IHPS faults:one is the leakage of the intelligent pump,and the other is the fatigue damage of the variable mechanism.

    When the intelligent pump operates normally,its output pressure varies from 21 MPa to 28 MPa.The input pressure instructions are set according to the pressure requirements under four phases of a flight profile shown in Table 2.The specific output pressure is shown in Fig.5.Obviously,the output pressure of the intelligent pumpPsfollows the command signalPuwell in normal operation.The residual errors of the swashplate angle and the output pressure obtained by the observer are almost equal to zero as shown in Fig.6(a)and(b),which means that the observer can approach a normal intelligent pump.

    Suppose that the IHPS is required to track a sine wave signal in the presence of system faults on different severity levels.The tracking trajectories and residual errors of the IHPS,for an input signalPu=(25+sint)×106under a normal condition,are shown in Fig.7.The residual errors of the swashplate angle and output pressure of the IHPS are almost equal to zero,as shown in Fig.8(a)and(b).

    In order to verify the fault diagnosis ability of the NUIO,the following failures are injected into the system:

    (1)The failure of the variable mechanism(simply expressed as fault I)is injected when the operational time>15 s by adjusting the spring stiffness of the variable mechanism to 50%of its original value.

    (2)The leakage failure of the intelligent pump(simply expressed as fault II)is injected when the operational time>15 s through increasing the total leakage coefficientCplto 105%of its original value.

    (3)The spring stiffness is adjusted to 50%,25%,10%,or 5%of its original value when the operational time>10 s.

    (4)The total leakage coefficient is adjusted to 105%,110%,115%,or 120%of its original value when the operational time>10 s.

    The simulation results are shown in Fig.9.

    In Fig.9(a)and(b),the residual errors are almost zero before the failure is injected.When the operational time=15 s,fault I is injected and the amplitude of the residual error of the swashplate angle is lower than zero,which means that a decrease of the spring stiffness directly leads to the failure of the variable mechanism.To reach torque balance with the output pressure set by the pressure signals command,the actual swashplate angle should increase with its new spring stiffness,which means that its angle amplitude should be higher than that of the ideal model.Fig.9(b)shows the increase of the residual error of the output pressure in 5 s,which means that the output pressure of the actual intelligent pump is lower than that of the ideal model.The result conforms to the fact that the variable mechanism becomes unable to produce high pressure when a failure happens.It can be said that the residual error is sensitive to fault I and can be used to detect this kind of fault.

    Table 1 Simulation parameters of the IHPS.

    In Fig.10(a),when the operational time=15 s,the residual error of the swashplate angle is larger than zero,that is,when fault II occurs,the actual value will be less than the observed one.The reason can be explained that with an increase of leakage in the intelligent pump,the torque used to balance the spring of the variable mechanism decreases.In order to reach a new balance state of torque,the actual value of the swashplate angle will be less than that observed.Meanwhile,as shown in Fig.10(b),the actual value is higher than the observed value,which can be explained as follows:in order to make up for the pressure loss caused by leakage,the actual output pressure should be higher than that calculated by the ideal model.Thus,the residual error of the output pressure is less than zero.As shown in Fig.10(a)and(b),the residual errors can also be used to detect fault II in real time.

    Table 2 Four operation modes of the intelligent pump.

    Fig.5 Output pressure Psfollowing the command signal Puin a normal condition.

    Fig.6 Residual errors of swashplate angle and output pressure of normal operation.

    Fig.7 Output pressure Psfollowing the command signal Puin a normal condition(Pu=(25+sin t)×106).

    Fig.8 Residual errors of swashplate angle and output pressure of normal operation(Pu=(25+sin t)×106).

    Fig.9 Residual errors of swashplate angle and output pressure of fault I.

    This simulation is designed to show the relationship between the severity of fault and the different value of spring stiffness.It is obvious in Fig.11(a)and(b)that the amplitudes of the residual errors of the swashplate angle and output pressure are constantly increasing,which means that,the higher the fault degree is,the greater the values of the residual errors will be.Therefore,the threshold of residual errors can be found via this method.The minimum acceptable level of fault degree can be set according to the simulation results.The threshold can reflect the early failure of the system quickly and its value must be within the range of the sensor’s resolution.The system is considered to be in the state of failure when the residual error exceeds the set threshold.

    Similar to the results obtained in simulation(3),as shown in Fig.12(a)and(b),the amplitudes of the residual errors of the swashplate angle and output pressure are constantly increasing respectively,but the directions of the residual errors are different from that of the variable mechanism failure.The more serious leakage failure happens,the greater the residual error values will be in different directions compared to another failure mode.

    Therefore,in terms of these two kinds of typical failure mode in the IHPS,the NUIO can accurately detect and isolate a failure.As shown in Fig.13,the failure modes distribute in different regions according to the simulation results,so a simple directional residual error can be used for fault isolation easily.At the same time,the threshold of the residual error can be found through this method.Some strategies can be used to determine the threshold and failure degree with an online method using the NUIO.

    Fig.10 Residual errors of swashplate angle and output pressure in fault II.

    Fig.11 Residual errors of swashplate angle and output pressure in different degrees of fault I.

    Fig.12 Residual errors of swashplate angle and output pressure in different degree of fault II.

    Fig.13 Failure distribution in different degrees of failure.

    Through the simulation above,the state observation of the dynamic IHPS is completed with the NUIO.The fault state of the system can be monitored with the residual error signal,and the fault detection of the system can be achieved by setting the threshold of failure in advance.

    5.Conclusions

    In this paper,model-based fault diagnosis is adopted to diagnose a real-time IHPS.Under an actual flight profile,with the output pressure and swashplate angle signals,the influences of nonlinear factors are taken into account and the dynamic performance is also embodied at the same time.Results show that the residual errors obtained from the NUIO can be used to diagnose the fault within the IHPS.It can be seen that the residual errors of the output pressure and swashplate angle signals could not only diagnose the health status of the IHPS accurately,but also identify the specific type of fault.The method proposed in this paper is an effective fault diagnosis method suitable for real-time systems including the IHPS as well as other hydraulic systems.

    Acknowledgements

    This study was co-supported by the National Natural Science Foundation of China(Nos.51620105010,51575019 and 51675019),National Basic Research Program of China(No.2014CB046400),and 111 Program of China.

    1.Rydberg KE.Energy efficient hydraulic systems and regenerative capabilities,Proceedings of the ninth scandinavian international conference on fluid power;2005 June 1–3;Sweden:Linko¨ping Universitet;2005.

    2.Chen B,Wang ZL,Qiu LH.Main developmental trend of aircraft hydraulic systems.Acta Aeronautica et Astronautica Sinica1998;19(7):S1–6[Chinese].

    3.Wang ZL,Chen B,Qiu LH.Trends of future aircraft hydraulic system.Hydraul Pnenmatics Seals2000;1:14–8[Chinese].

    4.Liu Y,Zuo MJ,Li YF,Huang HZ.Dynamic reliability assessment for multi-state systems utilizing system-level inspection data.IEEE Trans Reliab2015;64(4):1287–99.

    5.Wu ZM,Shi QL,Bai PF,Zhong F.Simulation analysis of load sensing system based on amesim.MachToolHydraul2013;44:38–41[Chinese].

    6.Mitchell JP.Load sensing hydraulic control system for variable displacement pump.United States patent US 20016216456;2001 Apr 17.

    7.Ma JG,Wang SF Wang ZL.Study of intelligent pump scheme.Chin Hydraul Pneumatics2002;11:6–8[Chinese].

    8.Zhou RX,Lin TQ,Han JD,Yan DC.Fault diagnosis of airplane hydraulic pump.Proceedings of the 4th world congress on intelligent control and automation.2002 June 10–14;Shanghai,China;2002.

    9.Lee JG,Kim OH.Development of a new hydraulic servo cylinder with mechanical feedback.Control Eng Practice1999;7(3):327–34.

    10.Wang XJ,Lin SR,Wang SP,He ZM,Zhang C.Remaining useful life prediction based on the wiener process for an aviation axial piston pump.Chin J Aeronaut2016;29(3):779–88.

    11.Li YH,Wang ZL.Development of airborne intelligent power supply system.J Beijing Univ Aeronaut Astronaut2004;30(6):493–7[Chinese].

    12.Wang SF,Ma JG,Wang ZL.Key technique for the research of airborne intelligent power supply system.Mach Tool Hydraul2003;4:85–7[Chinese].

    13.Lu CQ,Wang SP,Zhang C.Fault diagnosis of hydraulic piston pumps based on a two-step emd method and fuzzy c-means clustering.Proc Inst Mech Eng Part C:J Mech Eng Sci2016;230(16):2913–28.

    14.Wang JP,Hu HT.Vibration-based fault diagnosis of pump using fuzzy technique.Measurement2006;39(2):176–85.

    15.Gao Y,Zhang Q,Kong X.Wavelet–based pressure analysis for hydraulic pump health diagnosis.TransASAE2003;46(4):969–76.

    16.Lu C,Ma N,Wang ZP.Fault detection for hydraulic pump based on chaotic parallel rbf network.EURASIP J Adv Signal Process2011;2011(49):1–10.

    17.Piao XK.Research about pressure selection of hydraulic system in civil aircraft.Fluid Power Trans Control2011;6:22–4[Chinese].

    18.Ma J.Typical failure modes and accelerated lifetime test methods for aircraft hydraulic pump.Chin Hydraul Pneumatics2015;7:1–6[Chinese].

    19.Willsky AS.A survey of design methods for failure detection in dynamic systems.Automatica1976;12(6):601–11.

    20.Isermann R,Freyermuth B.Process fault diagnosis based on process model knowledge.Adv Inform Process Autom Control1990;31(4):21–34.

    21.Isermann R.Fault diagnosis of machines via parameter estimation and knowledge processing—tutorial paper.Automatica1993;29(4):815–35.

    22.Liu Y,Lin P,Li YF,Huang HZ.Bayesian reliability and performance assessment for multi-state systems.IEEE Trans Reliab2015;64(1):394–409.

    23.Palazzolo JJ,Scheunemann LD,Hartin JR.Leakage fault detection method for axial-piston variable displacement pumps.Proceedings of IEEE aerospace conference;2008 March 1–8;Montana,Big Sky:IEEE;2008.

    24.Liu S,Liu X,Han X,Cui X.Fault diagnosis of pump valve spring based on improved singularity analysis.J Vibroeng2014;16(2):704–12.

    25.Liu Y,Shi Y,Zhou Q,Xiu R.A sequential sampling strategy to improve the global fidelity of metamodels in multi-level system design.Struct Multidiscip Opt2016;53(6):1–19.

    26.Chen J,Patton RJ,Zhang HY.Design of unknown input observers and robust fault detection filters.Int J Control1996;63(1):85–105.

    27.Duan GR,Patton RJ.Robust fault detection using luenbergertype unknown input observers-a parametric approach.Int J Syst Sci2001;32(4):533–40.

    28.Isermann R,Balle′P.Trends in the application of model-based fault detection and diagnosis of technical processes.Control Eng Practice1997;5(5):709–19.

    29.Isermann R.Model base fault detection and diagnosis methods.Proceedings of the american control conference1995 June 21–23;Washington,Seattle:IEEE;2002.

    30.Isermann R.Model-based fault-detection and diagnosis–status and applications.Ann Rev Control2005;29(1):71–85.

    31.Mileti JA,Lawhead PM.Controlled pressure pumps for more efficient hydraulic systems.SAE technical paper;1986.Report No.:861844.

    32.Ma X,Chen B,Gan XH,Sun ZJ.The modeling and simulation of electro-hydraulic proportional valve-controlled cylinder position servo system.Mach Des Manuf2008;4:43–5.

    33.Huang XK,Zhang C,Wang SP,Liang XG.An accelerated life test model of hydraulic piston pumps based on dependency analysis:Model development.Proceedings of the 10th industrial electronics and applications(ICIEA)conference;2015 June 15–17;Auckland,New Zealand:IEEE;2015.

    人妻丰满熟妇av一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成av人片在线播放无| 伦理电影大哥的女人| 亚洲欧美日韩无卡精品| 最近视频中文字幕2019在线8| 免费观看的影片在线观看| 亚洲男人的天堂狠狠| 午夜视频国产福利| 老熟妇仑乱视频hdxx| 久久久久久久久久黄片| 青草久久国产| 啦啦啦观看免费观看视频高清| 99视频精品全部免费 在线| 国产不卡一卡二| 极品教师在线免费播放| 久久精品国产99精品国产亚洲性色| 少妇的逼水好多| 99久久久亚洲精品蜜臀av| 日韩欧美三级三区| 日韩亚洲欧美综合| 免费观看精品视频网站| 久久久久久九九精品二区国产| 深夜精品福利| 俺也久久电影网| 女人被狂操c到高潮| 怎么达到女性高潮| 亚洲中文字幕一区二区三区有码在线看| av国产免费在线观看| 欧美性猛交╳xxx乱大交人| 日韩 亚洲 欧美在线| 啦啦啦观看免费观看视频高清| 国产一级毛片七仙女欲春2| 香蕉av资源在线| 性插视频无遮挡在线免费观看| 亚洲欧美精品综合久久99| 亚洲,欧美,日韩| 国产三级在线视频| 中文字幕熟女人妻在线| 亚洲精品456在线播放app | 精品一区二区三区视频在线观看免费| 长腿黑丝高跟| 中文字幕精品亚洲无线码一区| 成人一区二区视频在线观看| 国产三级中文精品| 欧美日韩乱码在线| 99久国产av精品| 色5月婷婷丁香| 午夜福利在线在线| 免费大片18禁| 欧美极品一区二区三区四区| av天堂中文字幕网| 国产精品一及| 给我免费播放毛片高清在线观看| 国产精品野战在线观看| 久久人人精品亚洲av| 偷拍熟女少妇极品色| 夜夜夜夜夜久久久久| 91在线观看av| 3wmmmm亚洲av在线观看| 99久久无色码亚洲精品果冻| 观看免费一级毛片| 日本免费一区二区三区高清不卡| 一级毛片久久久久久久久女| 免费黄网站久久成人精品 | 好男人在线观看高清免费视频| 国产 一区 欧美 日韩| 中文亚洲av片在线观看爽| www.色视频.com| 97热精品久久久久久| 精品久久久久久久久久久久久| 国产精品1区2区在线观看.| 热99在线观看视频| 99riav亚洲国产免费| 精品人妻一区二区三区麻豆 | 九色国产91popny在线| 日韩中文字幕欧美一区二区| 哪里可以看免费的av片| 亚洲人成伊人成综合网2020| 久久久久国产精品人妻aⅴ院| 欧美bdsm另类| 美女cb高潮喷水在线观看| 国产高清三级在线| 久久欧美精品欧美久久欧美| 亚洲av熟女| 蜜桃久久精品国产亚洲av| .国产精品久久| 国产私拍福利视频在线观看| 三级毛片av免费| 午夜免费男女啪啪视频观看 | 我的老师免费观看完整版| 变态另类丝袜制服| 99热6这里只有精品| 国产亚洲精品av在线| 日韩成人在线观看一区二区三区| 成人无遮挡网站| 又黄又爽又刺激的免费视频.| 国产精品乱码一区二三区的特点| 国产高潮美女av| 婷婷色综合大香蕉| 久久久久久大精品| 亚洲五月婷婷丁香| 亚洲avbb在线观看| 美女被艹到高潮喷水动态| 天天一区二区日本电影三级| 精品一区二区三区视频在线观看免费| 丁香欧美五月| 天堂√8在线中文| 免费人成视频x8x8入口观看| 别揉我奶头~嗯~啊~动态视频| 色综合欧美亚洲国产小说| 免费观看精品视频网站| 亚洲无线观看免费| 一区二区三区激情视频| 国内揄拍国产精品人妻在线| 久久午夜福利片| 久久国产精品人妻蜜桃| 国产男靠女视频免费网站| 日本a在线网址| 一卡2卡三卡四卡精品乱码亚洲| 黄色日韩在线| 国产精品,欧美在线| 亚洲欧美日韩高清专用| 可以在线观看毛片的网站| 亚洲精品一区av在线观看| 最近中文字幕高清免费大全6 | 天堂网av新在线| 成人av一区二区三区在线看| 欧美性感艳星| 久久人妻av系列| 在线观看av片永久免费下载| 一个人免费在线观看电影| 欧美另类亚洲清纯唯美| 动漫黄色视频在线观看| 免费大片18禁| 男女那种视频在线观看| 精品久久久久久,| 免费搜索国产男女视频| 久久久久久久亚洲中文字幕 | 日韩欧美在线二视频| 国产精品久久久久久久电影| 国产三级黄色录像| 国产精品久久电影中文字幕| 一区二区三区激情视频| 超碰av人人做人人爽久久| 亚洲av电影不卡..在线观看| 12—13女人毛片做爰片一| a级毛片免费高清观看在线播放| 最新在线观看一区二区三区| 亚洲第一区二区三区不卡| 亚洲无线在线观看| 国内久久婷婷六月综合欲色啪| 久久精品91蜜桃| 国产一区二区激情短视频| 久久精品91蜜桃| 免费高清视频大片| 中出人妻视频一区二区| 99视频精品全部免费 在线| 精品久久国产蜜桃| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕av成人在线电影| 免费大片18禁| 级片在线观看| 久久国产乱子伦精品免费另类| 色精品久久人妻99蜜桃| 在线国产一区二区在线| 五月伊人婷婷丁香| 99在线视频只有这里精品首页| 亚洲无线观看免费| 97超视频在线观看视频| 性欧美人与动物交配| 日本成人三级电影网站| 深爱激情五月婷婷| 国产色婷婷99| 午夜a级毛片| 男人和女人高潮做爰伦理| 亚洲成av人片在线播放无| 国产主播在线观看一区二区| 国内精品久久久久久久电影| 亚洲人成网站在线播放欧美日韩| 免费观看人在逋| 天堂影院成人在线观看| 欧美又色又爽又黄视频| 欧美黑人巨大hd| 99久久精品一区二区三区| 精华霜和精华液先用哪个| 男插女下体视频免费在线播放| 午夜免费男女啪啪视频观看 | 色在线成人网| 18禁裸乳无遮挡免费网站照片| 免费一级毛片在线播放高清视频| www.色视频.com| bbb黄色大片| 欧美绝顶高潮抽搐喷水| 日韩欧美三级三区| 国内精品久久久久久久电影| 在线播放无遮挡| 亚洲五月天丁香| 国产人妻一区二区三区在| 国产大屁股一区二区在线视频| 99视频精品全部免费 在线| 男女做爰动态图高潮gif福利片| 嫩草影院入口| 国产又黄又爽又无遮挡在线| 不卡一级毛片| 日韩欧美三级三区| 一区二区三区免费毛片| 观看免费一级毛片| 亚洲成人久久性| or卡值多少钱| 精品日产1卡2卡| 国产精品98久久久久久宅男小说| 九九久久精品国产亚洲av麻豆| 最后的刺客免费高清国语| 国产aⅴ精品一区二区三区波| 身体一侧抽搐| 亚洲av成人av| 亚洲成人久久爱视频| 国产aⅴ精品一区二区三区波| 伦理电影大哥的女人| 中文字幕人成人乱码亚洲影| 日本a在线网址| 成人亚洲精品av一区二区| 国产免费男女视频| 人人妻,人人澡人人爽秒播| 亚洲av电影在线进入| 99精品久久久久人妻精品| 亚洲在线自拍视频| 免费看光身美女| 久久精品人妻少妇| 99在线人妻在线中文字幕| 禁无遮挡网站| 男女那种视频在线观看| 全区人妻精品视频| 国产私拍福利视频在线观看| 一区二区三区免费毛片| 高清在线国产一区| 国产精品98久久久久久宅男小说| 精品久久久久久久久亚洲 | 国产一区二区在线av高清观看| 亚洲国产色片| 中文字幕av在线有码专区| 国产av在哪里看| 少妇裸体淫交视频免费看高清| 亚洲国产色片| 久久久精品大字幕| 国产视频一区二区在线看| 国产精品av视频在线免费观看| 中亚洲国语对白在线视频| 高清在线国产一区| 国产精品亚洲一级av第二区| 国产91精品成人一区二区三区| 嫁个100分男人电影在线观看| 久久国产精品人妻蜜桃| 日本黄色视频三级网站网址| 日本与韩国留学比较| 久久久精品欧美日韩精品| 亚洲最大成人手机在线| 国内精品久久久久精免费| 成人av一区二区三区在线看| 国产成人福利小说| 1024手机看黄色片| av专区在线播放| 国产精品影院久久| 乱人视频在线观看| 99久国产av精品| 在线a可以看的网站| 男女床上黄色一级片免费看| 久99久视频精品免费| 亚洲无线观看免费| 久久精品夜夜夜夜夜久久蜜豆| 深爱激情五月婷婷| 欧美3d第一页| 国产精品亚洲美女久久久| 国产一区二区在线观看日韩| 亚洲aⅴ乱码一区二区在线播放| 国产精品永久免费网站| 国产伦在线观看视频一区| 久久久久亚洲av毛片大全| 舔av片在线| 国产成人欧美在线观看| 日韩精品中文字幕看吧| 他把我摸到了高潮在线观看| 禁无遮挡网站| 国产欧美日韩精品一区二区| 此物有八面人人有两片| 一夜夜www| 日韩欧美国产在线观看| 成人毛片a级毛片在线播放| 中文字幕熟女人妻在线| 亚洲欧美精品综合久久99| 中文亚洲av片在线观看爽| 欧美+亚洲+日韩+国产| 天堂√8在线中文| 欧美成狂野欧美在线观看| 国产欧美日韩精品一区二区| 亚洲欧美日韩无卡精品| 99国产综合亚洲精品| 一卡2卡三卡四卡精品乱码亚洲| 国产视频一区二区在线看| 亚洲国产欧洲综合997久久,| 日韩亚洲欧美综合| 9191精品国产免费久久| 一a级毛片在线观看| 欧美一区二区国产精品久久精品| 国产成人aa在线观看| 99热精品在线国产| 午夜免费成人在线视频| 嫩草影院新地址| 乱码一卡2卡4卡精品| 国产精品综合久久久久久久免费| 国产精品女同一区二区软件 | 在线播放国产精品三级| 免费在线观看影片大全网站| 美女xxoo啪啪120秒动态图 | 精品国产三级普通话版| 欧美zozozo另类| 老司机午夜十八禁免费视频| 激情在线观看视频在线高清| 日韩精品青青久久久久久| 亚洲av第一区精品v没综合| 91九色精品人成在线观看| 国产三级黄色录像| 精品久久久久久久人妻蜜臀av| 长腿黑丝高跟| 美女被艹到高潮喷水动态| 国产高清视频在线观看网站| 国产在线精品亚洲第一网站| 欧美黄色片欧美黄色片| 又黄又爽又刺激的免费视频.| 亚洲综合色惰| 又爽又黄a免费视频| 国内揄拍国产精品人妻在线| 久久伊人香网站| 国产不卡一卡二| 久久草成人影院| 日本与韩国留学比较| x7x7x7水蜜桃| 久久人人精品亚洲av| 很黄的视频免费| 亚洲 国产 在线| 久99久视频精品免费| 麻豆久久精品国产亚洲av| 国产精华一区二区三区| 永久网站在线| 久久午夜亚洲精品久久| 日韩欧美 国产精品| 欧美成狂野欧美在线观看| 夜夜看夜夜爽夜夜摸| 97热精品久久久久久| 亚洲精品一卡2卡三卡4卡5卡| 欧美在线一区亚洲| 三级男女做爰猛烈吃奶摸视频| 色播亚洲综合网| 久久久久国产精品人妻aⅴ院| 天美传媒精品一区二区| 婷婷亚洲欧美| 丝袜美腿在线中文| 久久国产精品影院| 久久性视频一级片| 国产免费一级a男人的天堂| 亚洲av免费高清在线观看| 动漫黄色视频在线观看| 99在线视频只有这里精品首页| 日本黄大片高清| 国产精华一区二区三区| 12—13女人毛片做爰片一| 久9热在线精品视频| 国内精品美女久久久久久| 国产亚洲av嫩草精品影院| 亚洲精品在线美女| 成人av在线播放网站| 精品一区二区三区视频在线观看免费| 欧美zozozo另类| 中文字幕av在线有码专区| av欧美777| 精品一区二区三区视频在线| 99国产精品一区二区三区| 日本 av在线| 日韩有码中文字幕| 婷婷亚洲欧美| 国产精品自产拍在线观看55亚洲| 亚洲av成人精品一区久久| 一级a爱片免费观看的视频| 欧美乱妇无乱码| 欧美日韩国产亚洲二区| 亚洲片人在线观看| 亚洲人成电影免费在线| 内射极品少妇av片p| 国产成年人精品一区二区| 又爽又黄无遮挡网站| 悠悠久久av| 欧美成人a在线观看| 亚洲国产日韩欧美精品在线观看| 亚洲,欧美精品.| 特大巨黑吊av在线直播| 国产高潮美女av| 精品久久久久久久末码| 91狼人影院| 亚洲自偷自拍三级| 久久伊人香网站| 久久国产乱子伦精品免费另类| 国产淫片久久久久久久久 | 亚洲美女视频黄频| 国内精品久久久久久久电影| 精品国产亚洲在线| 日韩av在线大香蕉| 亚洲欧美日韩东京热| 女人十人毛片免费观看3o分钟| 国产乱人伦免费视频| 性色avwww在线观看| 国产一级毛片七仙女欲春2| 午夜精品久久久久久毛片777| 好男人电影高清在线观看| 99久久精品国产亚洲精品| 成人三级黄色视频| 日韩欧美在线二视频| 宅男免费午夜| 99久久精品一区二区三区| 亚洲狠狠婷婷综合久久图片| 亚洲欧美日韩东京热| 日本 欧美在线| 女同久久另类99精品国产91| 亚洲国产色片| 人人妻人人澡欧美一区二区| 亚洲人成网站在线播| 黄色女人牲交| 99精品在免费线老司机午夜| 噜噜噜噜噜久久久久久91| 久久久久久久久中文| 免费大片18禁| 免费av不卡在线播放| 青草久久国产| 日韩欧美国产一区二区入口| 老司机福利观看| 免费av毛片视频| 亚洲人成网站高清观看| 亚洲av电影在线进入| 日韩国内少妇激情av| 最近最新免费中文字幕在线| 国产一区二区三区视频了| 精品人妻熟女av久视频| 美女大奶头视频| 国产精品永久免费网站| 久久久久久久久大av| 日本一本二区三区精品| 90打野战视频偷拍视频| 一个人看视频在线观看www免费| 男女视频在线观看网站免费| 永久网站在线| 99热这里只有是精品在线观看 | 欧美一级a爱片免费观看看| 成年女人永久免费观看视频| 日韩国内少妇激情av| 中国美女看黄片| 12—13女人毛片做爰片一| 99久久精品国产亚洲精品| 国产不卡一卡二| 国产色爽女视频免费观看| 国产黄a三级三级三级人| 亚洲国产日韩欧美精品在线观看| 夜夜爽天天搞| 国产精品一区二区性色av| 亚洲欧美激情综合另类| 国产主播在线观看一区二区| 精品国内亚洲2022精品成人| 午夜亚洲福利在线播放| 国产亚洲精品久久久久久毛片| 一级作爱视频免费观看| 亚洲欧美日韩无卡精品| 亚洲不卡免费看| 日本 av在线| 欧美高清性xxxxhd video| 国内揄拍国产精品人妻在线| 国产单亲对白刺激| 啦啦啦韩国在线观看视频| 日韩欧美精品v在线| 亚洲欧美激情综合另类| 嫩草影视91久久| 无人区码免费观看不卡| 国内少妇人妻偷人精品xxx网站| 免费观看的影片在线观看| 日本成人三级电影网站| www.www免费av| 欧美日韩综合久久久久久 | 中文字幕熟女人妻在线| 精品久久久久久久人妻蜜臀av| 国产高清视频在线播放一区| av在线天堂中文字幕| 欧美成狂野欧美在线观看| 三级国产精品欧美在线观看| 简卡轻食公司| 在线播放无遮挡| 亚洲最大成人av| netflix在线观看网站| 欧美一区二区国产精品久久精品| 少妇的逼水好多| 女生性感内裤真人,穿戴方法视频| 看片在线看免费视频| 日韩欧美一区二区三区在线观看| 久99久视频精品免费| 亚洲真实伦在线观看| 老司机福利观看| 久久精品夜夜夜夜夜久久蜜豆| 99久久成人亚洲精品观看| 久久精品国产亚洲av香蕉五月| 一本精品99久久精品77| 国产午夜精品久久久久久一区二区三区 | 夜夜爽天天搞| h日本视频在线播放| 丝袜美腿在线中文| 亚洲av五月六月丁香网| 久久久久久久精品吃奶| 激情在线观看视频在线高清| 丝袜美腿在线中文| 网址你懂的国产日韩在线| 性插视频无遮挡在线免费观看| 国语自产精品视频在线第100页| 欧美激情在线99| 国产精品99久久久久久久久| 日韩欧美免费精品| 内射极品少妇av片p| 最近视频中文字幕2019在线8| 小蜜桃在线观看免费完整版高清| 午夜福利免费观看在线| 中文字幕av在线有码专区| 麻豆国产97在线/欧美| 脱女人内裤的视频| АⅤ资源中文在线天堂| 国产激情偷乱视频一区二区| 90打野战视频偷拍视频| 在线观看舔阴道视频| 我要搜黄色片| 欧美日韩瑟瑟在线播放| 亚洲欧美日韩卡通动漫| 三级男女做爰猛烈吃奶摸视频| 禁无遮挡网站| 91午夜精品亚洲一区二区三区 | 日韩大尺度精品在线看网址| АⅤ资源中文在线天堂| 亚洲成av人片在线播放无| 51午夜福利影视在线观看| 国产精品乱码一区二三区的特点| 精品免费久久久久久久清纯| 成人高潮视频无遮挡免费网站| 一个人免费在线观看电影| 欧美+日韩+精品| 69av精品久久久久久| av福利片在线观看| 男插女下体视频免费在线播放| av在线蜜桃| 国产成人av教育| 久久人人精品亚洲av| 校园春色视频在线观看| 高清日韩中文字幕在线| 亚洲在线观看片| 听说在线观看完整版免费高清| 欧美日韩瑟瑟在线播放| 亚洲av熟女| 免费av不卡在线播放| 在现免费观看毛片| 神马国产精品三级电影在线观看| 色视频www国产| 丁香欧美五月| 亚洲第一电影网av| 国内揄拍国产精品人妻在线| 国产aⅴ精品一区二区三区波| 天堂√8在线中文| 国产精品不卡视频一区二区 | 天天一区二区日本电影三级| 99国产精品一区二区蜜桃av| 国产精品乱码一区二三区的特点| 中国美女看黄片| 91字幕亚洲| 成熟少妇高潮喷水视频| 丰满人妻一区二区三区视频av| 欧美最黄视频在线播放免费| 综合色av麻豆| 少妇熟女aⅴ在线视频| 一区二区三区免费毛片| 免费在线观看影片大全网站| 欧美色欧美亚洲另类二区| 在线观看一区二区三区| 十八禁国产超污无遮挡网站| 琪琪午夜伦伦电影理论片6080| 熟妇人妻久久中文字幕3abv| 嫩草影院入口| 亚洲成人久久性| 国产成人欧美在线观看| 色综合亚洲欧美另类图片| 欧美最黄视频在线播放免费| 18禁在线播放成人免费| 真人做人爱边吃奶动态| 日本a在线网址| 欧美+亚洲+日韩+国产| 观看免费一级毛片| 欧美激情在线99| 亚洲专区国产一区二区| 性欧美人与动物交配| 黄片小视频在线播放| 成人精品一区二区免费| 欧美成人性av电影在线观看| 小说图片视频综合网站| 国产亚洲精品久久久com| 日本免费一区二区三区高清不卡| 又爽又黄无遮挡网站| 精华霜和精华液先用哪个| av在线老鸭窝| 欧美一区二区亚洲| 亚洲人成网站在线播放欧美日韩| 亚洲成人久久爱视频| 亚洲片人在线观看| 精品国内亚洲2022精品成人| 一区二区三区四区激情视频 | 黄片小视频在线播放|