• <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成人不卡在线观看播放网| 亚洲人成电影观看| 亚洲国产欧美一区二区综合| 色精品久久人妻99蜜桃| 午夜精品久久久久久毛片777| 国产片内射在线| 国产精品国产高清国产av| 国产野战对白在线观看| 在线观看一区二区三区| 电影成人av| 色婷婷av一区二区三区视频| 国产精品一区二区免费欧美| 欧美日韩av久久| 侵犯人妻中文字幕一二三四区| 久热爱精品视频在线9| 美女大奶头视频| 午夜福利在线免费观看网站| а√天堂www在线а√下载| 高清毛片免费观看视频网站 | 中文字幕精品免费在线观看视频| 亚洲成人免费av在线播放| av有码第一页| 真人一进一出gif抽搐免费| 中文字幕最新亚洲高清| 91九色精品人成在线观看| xxxhd国产人妻xxx| av在线天堂中文字幕 | 久久久久久免费高清国产稀缺| 9热在线视频观看99| 少妇 在线观看| 91大片在线观看| 亚洲avbb在线观看| 丝袜人妻中文字幕| 国产色视频综合| 久久久国产成人免费| 日本wwww免费看| 18禁黄网站禁片午夜丰满| 免费在线观看亚洲国产| 久久久国产成人精品二区 | 老汉色av国产亚洲站长工具| 免费在线观看黄色视频的| 国产99久久九九免费精品| 亚洲精品在线美女| 亚洲精品国产精品久久久不卡| 亚洲色图 男人天堂 中文字幕| 久久影院123| 欧美乱码精品一区二区三区| 亚洲精品国产色婷婷电影| 国内毛片毛片毛片毛片毛片| 精品国产美女av久久久久小说| 在线免费观看的www视频| 成年人免费黄色播放视频| 久久香蕉国产精品| 99国产综合亚洲精品| 性色av乱码一区二区三区2| 亚洲一码二码三码区别大吗| 亚洲熟妇熟女久久| 99国产综合亚洲精品| 国产精品久久久久久人妻精品电影| av超薄肉色丝袜交足视频| 一a级毛片在线观看| 18禁国产床啪视频网站| 国产成人精品在线电影| 在线av久久热| 欧美日韩精品网址| 国产欧美日韩一区二区三| 最近最新免费中文字幕在线| 丰满迷人的少妇在线观看| 丰满饥渴人妻一区二区三| 欧美老熟妇乱子伦牲交| 日韩有码中文字幕| 中文字幕高清在线视频| 亚洲自拍偷在线| 91av网站免费观看| 久久久久国产一级毛片高清牌| 国产伦一二天堂av在线观看| 啦啦啦免费观看视频1| 国产成人精品无人区| 欧美成人免费av一区二区三区| 欧美性长视频在线观看| 在线观看免费视频网站a站| 免费人成视频x8x8入口观看| 高清在线国产一区| 少妇粗大呻吟视频| 五月开心婷婷网| 18禁黄网站禁片午夜丰满| 久久久久久久久中文| 97碰自拍视频| 日韩欧美免费精品| svipshipincom国产片| 变态另类成人亚洲欧美熟女 | 国产精品日韩av在线免费观看 | 男男h啪啪无遮挡| 韩国精品一区二区三区| 天天添夜夜摸| 日韩一卡2卡3卡4卡2021年| av电影中文网址| 精品久久久久久,| 亚洲精品中文字幕一二三四区| 亚洲欧美精品综合久久99| www日本在线高清视频| 午夜福利欧美成人| 国产色视频综合| 脱女人内裤的视频| 97碰自拍视频| 亚洲黑人精品在线| 成人18禁在线播放| 三上悠亚av全集在线观看| videosex国产| www国产在线视频色| 女性被躁到高潮视频| av欧美777| 麻豆av在线久日| 两个人看的免费小视频| 自线自在国产av| 国内毛片毛片毛片毛片毛片| 熟女少妇亚洲综合色aaa.| 欧美精品亚洲一区二区| 19禁男女啪啪无遮挡网站| 18禁黄网站禁片午夜丰满| 国产亚洲精品一区二区www| 亚洲国产精品sss在线观看 | 女性被躁到高潮视频| 午夜福利一区二区在线看| 国产精品偷伦视频观看了| 久久影院123| 少妇粗大呻吟视频| 国产精品一区二区三区四区久久 | 80岁老熟妇乱子伦牲交| 操出白浆在线播放| 亚洲午夜精品一区,二区,三区| 精品国产美女av久久久久小说| 欧美日韩视频精品一区| 精品国产国语对白av| 一级a爱片免费观看的视频| 色精品久久人妻99蜜桃| 亚洲av片天天在线观看| 久热这里只有精品99| 高清毛片免费观看视频网站 | 在线免费观看的www视频| 国产欧美日韩综合在线一区二区| 久久久久国产精品人妻aⅴ院| 午夜福利一区二区在线看| 欧美日韩亚洲综合一区二区三区_| 国产成人精品在线电影| 国产成年人精品一区二区 | 高清欧美精品videossex| 国产又色又爽无遮挡免费看| 免费看十八禁软件| 热re99久久国产66热| 亚洲欧美日韩无卡精品| 又大又爽又粗| 热re99久久国产66热| 国产真人三级小视频在线观看| 丝袜人妻中文字幕| 无限看片的www在线观看| 国产高清激情床上av| 在线观看免费日韩欧美大片| 午夜91福利影院| 一区二区日韩欧美中文字幕| 人人妻人人澡人人看| 一本大道久久a久久精品| 精品人妻在线不人妻| 中文欧美无线码| 色综合婷婷激情| 美女午夜性视频免费| 国产色视频综合| 咕卡用的链子| 亚洲激情在线av| 女警被强在线播放| 水蜜桃什么品种好| 大陆偷拍与自拍| 熟女少妇亚洲综合色aaa.| 亚洲精品国产区一区二| av网站免费在线观看视频| 黄色怎么调成土黄色| 午夜影院日韩av| 亚洲国产欧美一区二区综合| 欧美乱码精品一区二区三区| 美女高潮喷水抽搐中文字幕| 日本撒尿小便嘘嘘汇集6| 欧美日韩av久久| 日韩三级视频一区二区三区| 多毛熟女@视频| ponron亚洲| 最近最新中文字幕大全免费视频| 真人做人爱边吃奶动态| 久久中文看片网| 国产麻豆69| a级毛片黄视频| 亚洲激情在线av| 亚洲性夜色夜夜综合| 亚洲伊人色综图| 80岁老熟妇乱子伦牲交| ponron亚洲| 日韩视频一区二区在线观看| 可以免费在线观看a视频的电影网站| 老司机深夜福利视频在线观看| 日韩欧美国产一区二区入口| 成在线人永久免费视频| 国产精华一区二区三区| 亚洲一区高清亚洲精品| 中文字幕最新亚洲高清| 在线国产一区二区在线| 99久久人妻综合| 日韩人妻精品一区2区三区| 日日夜夜操网爽| 美女午夜性视频免费| 麻豆国产av国片精品| 不卡av一区二区三区| 男女床上黄色一级片免费看| 亚洲精品国产一区二区精华液| 欧美亚洲日本最大视频资源| 精品免费久久久久久久清纯| 啦啦啦在线免费观看视频4| 男人舔女人下体高潮全视频| 性少妇av在线| 女人被躁到高潮嗷嗷叫费观| 亚洲一区二区三区欧美精品| 99国产精品一区二区三区| 男女床上黄色一级片免费看| 国内久久婷婷六月综合欲色啪| www.自偷自拍.com| 欧美成狂野欧美在线观看| 国产免费男女视频| 女性生殖器流出的白浆| 人人澡人人妻人| 国产野战对白在线观看| 国产国语露脸激情在线看| 变态另类成人亚洲欧美熟女 | 久久精品成人免费网站| 久久精品亚洲熟妇少妇任你| 国产成人精品在线电影| 免费在线观看日本一区| 国产精品久久久人人做人人爽| 人人妻人人澡人人看| 侵犯人妻中文字幕一二三四区| 99精国产麻豆久久婷婷| 国产精品野战在线观看 | 久久 成人 亚洲| 少妇被粗大的猛进出69影院| 国产成人欧美在线观看| ponron亚洲| 亚洲五月天丁香| 看黄色毛片网站| 久久天堂一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 国产欧美日韩一区二区三区在线| 一级片免费观看大全| 亚洲 欧美一区二区三区| 在线观看一区二区三区| 国产单亲对白刺激| av天堂在线播放| 国产色视频综合| 亚洲午夜精品一区,二区,三区| 90打野战视频偷拍视频| 亚洲色图综合在线观看| 国产91精品成人一区二区三区| 最新在线观看一区二区三区| 嫩草影视91久久| 亚洲五月婷婷丁香| 757午夜福利合集在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 欧美在线一区亚洲| 国产伦一二天堂av在线观看| 国产精品日韩av在线免费观看 | 热re99久久精品国产66热6| 高清黄色对白视频在线免费看| 国产精品野战在线观看 | 国产成人精品久久二区二区免费| 亚洲av美国av| 级片在线观看| 一区福利在线观看| 精品久久蜜臀av无| 日韩成人在线观看一区二区三区| 精品久久久久久电影网| 亚洲全国av大片| 这个男人来自地球电影免费观看| 欧美日韩精品网址| 大型av网站在线播放| 丝袜美足系列| 国产一区二区在线av高清观看| 天堂影院成人在线观看| 久热这里只有精品99| 女人被狂操c到高潮| 纯流量卡能插随身wifi吗| 日韩大码丰满熟妇| 在线观看66精品国产| 中文字幕另类日韩欧美亚洲嫩草| 精品国产国语对白av| 99在线人妻在线中文字幕| 在线播放国产精品三级| 日韩欧美一区视频在线观看| 人妻久久中文字幕网| 99国产极品粉嫩在线观看| 成人免费观看视频高清| 国产熟女午夜一区二区三区| 天天添夜夜摸| 国产精品亚洲av一区麻豆| 又黄又爽又免费观看的视频| 国产精品久久久人人做人人爽| 18美女黄网站色大片免费观看| 亚洲精品久久午夜乱码| 亚洲成av片中文字幕在线观看| 色尼玛亚洲综合影院| 巨乳人妻的诱惑在线观看| 老司机午夜十八禁免费视频| 12—13女人毛片做爰片一| 最近最新免费中文字幕在线| 亚洲精品一卡2卡三卡4卡5卡| 免费观看人在逋| 久久国产亚洲av麻豆专区| 在线免费观看的www视频| svipshipincom国产片| 国产又爽黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 乱人伦中国视频| av福利片在线| 女人被躁到高潮嗷嗷叫费观| 在线观看免费日韩欧美大片| 激情视频va一区二区三区| 久久天躁狠狠躁夜夜2o2o| 色婷婷久久久亚洲欧美| 大型av网站在线播放| 欧美 亚洲 国产 日韩一| 国产精品av久久久久免费| 一级毛片精品| 国产精品九九99| 黄色a级毛片大全视频| 他把我摸到了高潮在线观看| 法律面前人人平等表现在哪些方面| 亚洲欧美日韩无卡精品| 久久久久国产一级毛片高清牌| 欧美日本亚洲视频在线播放| 女人被狂操c到高潮| av中文乱码字幕在线| 热99re8久久精品国产| 国产精品综合久久久久久久免费 | 国产黄色免费在线视频| 国产精品98久久久久久宅男小说| 麻豆av在线久日| 午夜激情av网站| 满18在线观看网站| 亚洲午夜精品一区,二区,三区| 亚洲欧美一区二区三区黑人| 女性被躁到高潮视频| 亚洲欧美一区二区三区黑人| 欧美 亚洲 国产 日韩一| 国产有黄有色有爽视频| 一本大道久久a久久精品| 国产有黄有色有爽视频| 免费日韩欧美在线观看| 日日干狠狠操夜夜爽| 欧美乱色亚洲激情| 精品免费久久久久久久清纯| 激情在线观看视频在线高清| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美精品综合一区二区三区| 久久久久国内视频| 女性被躁到高潮视频| 一进一出抽搐动态| 亚洲成国产人片在线观看| 欧美日韩瑟瑟在线播放| 99国产精品免费福利视频| 亚洲第一欧美日韩一区二区三区| 午夜视频精品福利| 国产黄色免费在线视频| 欧美人与性动交α欧美软件| 99国产精品免费福利视频| 婷婷丁香在线五月| 一二三四在线观看免费中文在| 国产欧美日韩一区二区三| 真人一进一出gif抽搐免费| 一区二区三区国产精品乱码| 99re在线观看精品视频| 99国产综合亚洲精品| 国产成人欧美在线观看| 又黄又爽又免费观看的视频| 深夜精品福利| 国产极品粉嫩免费观看在线| 色综合婷婷激情| 精品人妻1区二区| 国产精品久久久久久人妻精品电影| 91大片在线观看| 黄网站色视频无遮挡免费观看| 成人亚洲精品一区在线观看| av片东京热男人的天堂| 欧美日韩av久久| 五月开心婷婷网| 日本三级黄在线观看| 亚洲伊人色综图| 97超级碰碰碰精品色视频在线观看| 欧美日本亚洲视频在线播放| 欧美成人免费av一区二区三区| 国产又色又爽无遮挡免费看| 成人影院久久| 日韩三级视频一区二区三区| 午夜91福利影院| 中文字幕另类日韩欧美亚洲嫩草| 手机成人av网站| 色精品久久人妻99蜜桃| 国产无遮挡羞羞视频在线观看| 亚洲色图 男人天堂 中文字幕| 少妇粗大呻吟视频| av片东京热男人的天堂| 一a级毛片在线观看| 法律面前人人平等表现在哪些方面| 午夜福利一区二区在线看| av天堂久久9| 天堂动漫精品| 国产亚洲精品综合一区在线观看 | 国产亚洲精品综合一区在线观看 | 亚洲中文字幕日韩| 中文字幕人妻熟女乱码| 亚洲少妇的诱惑av| 大码成人一级视频| 婷婷六月久久综合丁香| 久久热在线av| 色综合站精品国产| 国产精品久久久久成人av| 久久久久久久精品吃奶| 女同久久另类99精品国产91| 国产熟女午夜一区二区三区| av中文乱码字幕在线| 久久久久国产一级毛片高清牌| 精品一品国产午夜福利视频| 91老司机精品| 女性生殖器流出的白浆| 不卡av一区二区三区| 国产精品 国内视频| 亚洲中文字幕日韩| 视频区图区小说| 亚洲自拍偷在线| 亚洲精品国产一区二区精华液| 91麻豆精品激情在线观看国产 | 精品国产乱码久久久久久男人| 亚洲国产精品合色在线| 在线观看免费日韩欧美大片| 欧美性长视频在线观看| 中文字幕人妻丝袜制服| 欧美不卡视频在线免费观看 | 国产av在哪里看| 可以免费在线观看a视频的电影网站| 亚洲片人在线观看| 欧美久久黑人一区二区| 国产三级黄色录像| 日韩国内少妇激情av| av网站免费在线观看视频| 日本黄色视频三级网站网址| 热99re8久久精品国产| 成人三级黄色视频| 久久影院123| 国产三级黄色录像| 欧美国产精品va在线观看不卡| 1024视频免费在线观看| 欧美一区二区精品小视频在线| 女人精品久久久久毛片| 亚洲少妇的诱惑av| 精品国产一区二区久久| 无遮挡黄片免费观看| 国产单亲对白刺激| 国产亚洲精品久久久久久毛片| 免费av毛片视频| 深夜精品福利| 黑人猛操日本美女一级片| 亚洲美女黄片视频| 欧美日韩一级在线毛片| 精品一区二区三区视频在线观看免费 | 婷婷六月久久综合丁香| 国产精品亚洲av一区麻豆| 日韩免费av在线播放| 大香蕉久久成人网| 欧美精品亚洲一区二区| 香蕉丝袜av| 午夜日韩欧美国产| 国产亚洲欧美在线一区二区| 久久影院123| 国产高清激情床上av| 十分钟在线观看高清视频www| 国产深夜福利视频在线观看| 午夜精品久久久久久毛片777| 亚洲 欧美 日韩 在线 免费| 国产精品久久视频播放| 9热在线视频观看99| 最近最新中文字幕大全免费视频| 国产视频一区二区在线看| 99热国产这里只有精品6| 一级a爱片免费观看的视频| 丰满的人妻完整版| 亚洲欧洲精品一区二区精品久久久| 十八禁网站免费在线| 香蕉久久夜色| 午夜免费激情av| 人人妻,人人澡人人爽秒播| 校园春色视频在线观看| 亚洲熟妇中文字幕五十中出 | 中国美女看黄片| 午夜精品国产一区二区电影| 自拍欧美九色日韩亚洲蝌蚪91| 91麻豆精品激情在线观看国产 | 麻豆国产av国片精品| xxxhd国产人妻xxx| 久久久国产精品麻豆| 精品国产国语对白av| 国产伦人伦偷精品视频| 久久精品亚洲精品国产色婷小说| 天天添夜夜摸| 看黄色毛片网站| 国产三级黄色录像| 9191精品国产免费久久| 久久久久久大精品| netflix在线观看网站| 少妇裸体淫交视频免费看高清 | 欧美日韩瑟瑟在线播放| 美女高潮喷水抽搐中文字幕| 一级毛片精品| 午夜福利在线免费观看网站| 中文字幕精品免费在线观看视频| 国产成人精品无人区| 水蜜桃什么品种好| 久久久久久久久久久久大奶| av免费在线观看网站| 亚洲五月婷婷丁香| 老熟妇仑乱视频hdxx| 午夜老司机福利片| www.自偷自拍.com| 女人精品久久久久毛片| 天天影视国产精品| 精品久久蜜臀av无| 成人三级黄色视频| 久久久国产一区二区| av欧美777| 国产激情久久老熟女| 久久国产精品男人的天堂亚洲| 欧美成人免费av一区二区三区| 日韩 欧美 亚洲 中文字幕| 熟女少妇亚洲综合色aaa.| 一区福利在线观看| 亚洲欧洲精品一区二区精品久久久| 国产91精品成人一区二区三区| 夜夜躁狠狠躁天天躁| 成年人免费黄色播放视频| 美女福利国产在线| 嫩草影院精品99| 老司机深夜福利视频在线观看| 久久久久久久精品吃奶| 亚洲专区国产一区二区| 久久久久国内视频| 亚洲一区二区三区色噜噜 | 亚洲欧美日韩高清在线视频| 亚洲精品美女久久久久99蜜臀| 亚洲精品av麻豆狂野| 久久香蕉国产精品| 国产熟女xx| 亚洲一码二码三码区别大吗| ponron亚洲| 黑人猛操日本美女一级片| 精品久久久久久久毛片微露脸| 国产亚洲av高清不卡| 最新在线观看一区二区三区| 国产精品一区二区三区四区久久 | 午夜久久久在线观看| 国产精品日韩av在线免费观看 | 18禁国产床啪视频网站| 高清欧美精品videossex| 国产欧美日韩一区二区精品| av欧美777| 亚洲av成人一区二区三| 男人的好看免费观看在线视频 | 国产在线精品亚洲第一网站| 中文字幕人妻丝袜制服| 女人被躁到高潮嗷嗷叫费观| 多毛熟女@视频| 99精国产麻豆久久婷婷| 亚洲一区二区三区色噜噜 | 亚洲成人精品中文字幕电影 | 久久精品国产清高在天天线| 一本综合久久免费| 久久久国产成人免费| 久热这里只有精品99| 韩国av一区二区三区四区| 久久午夜亚洲精品久久| 国产精品香港三级国产av潘金莲| 亚洲中文日韩欧美视频| 国产精华一区二区三区| 久久狼人影院| 亚洲午夜理论影院| 两性夫妻黄色片| 日韩 欧美 亚洲 中文字幕| 丁香六月欧美| www.www免费av| 国产成人系列免费观看| 久久国产亚洲av麻豆专区| 丝袜人妻中文字幕| 天天影视国产精品| 看黄色毛片网站| 国产精品成人在线| 午夜a级毛片| 国产熟女xx| 一级作爱视频免费观看| 另类亚洲欧美激情| 日韩国内少妇激情av| 人人妻人人澡人人看| 免费在线观看完整版高清| 亚洲一区中文字幕在线| 日本vs欧美在线观看视频| 亚洲精品美女久久久久99蜜臀| 国产精品1区2区在线观看.| 女警被强在线播放| 亚洲av美国av| 国产精品国产高清国产av| 性少妇av在线| 亚洲熟妇熟女久久| 在线观看免费视频日本深夜| 狂野欧美激情性xxxx|