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

    Review on Key Techniques of Dynamics Modeling for Helicopter Simulator

    2016-02-09 01:52:37,,

    , ,

    Research Center of Flight Simulation and Advanced Engineering Training, College of Civil Aviation,Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P. R. China

    ?

    Review on Key Techniques of Dynamics Modeling for Helicopter Simulator

    GuHongbin*,FuJun,HuJinshuo

    Research Center of Flight Simulation and Advanced Engineering Training, College of Civil Aviation,Nanjing University of Aeronautics and Astronautics, Nanjing 210016, P. R. China

    (Received 22 August 2016; revised 16 November 2016; accepted 20 November 2016)

    We review the state of the art of dynamic modeling methods for helicopter simulators. Major dynamic models, including ARMCOP, GENHEL, FLYRT, BEMAP, F-B412 and UTIAS, are introduced. We address off-axis response problem by focusing on flexible blade model, Pitt-Peters dynamic inflow model, Peters-He finite state inflow model and free wake model which are integrated with the baseline model. With the advances in computing power, efficient free-vortex algorithms and parallel processing, free wake model can be used to simulate in real-time and significantly improves the effectiveness of solving off-axis response problem.

    helicopter; dynamic model; off-aixs response

    0 Introduction

    Flight simulators are indispensable tools for both training and researching. Fig.1 shows the train simulator developed by Sikorsky Aircraft Corporation[1], and Fig.2 the research simulator developed by The University of Liverpool[2]. In the effort of replicating complex flying maneuvers and conditions of the real aircraft, both the simulator mathematical model and the motion and visual systems are becoming increasingly sophisticated. A flight simulator cannot perfectly represent the aircraft in all aspects: The mathematical model of the aircraft has never been fully accurate, and the motion and visual systems have physical limitations that impede the full representation of the sensation of flying always toward perfection. Also, sophisticated simulator mathematical models are not convenient for obtaining physical insight into the mechanism of coupling between the different motion degrees of freedom.

    Regulatory authorities have established standards for the acceptance and certification of simulators. For helicopters, the most widely recognized standard for simulation qualification is Federal Aviation Administration advisory circular AC 120-63[3], which is currently being reworked into joint-aviation-requirement (JAR) STD-1H[4]. The mathematical model (including sophisticated kinematics, structural, aerodynamic, and control dynamic aspects) is vital to flight simulators which must run in real-time. Flight simulator′s fidelity depends on mathematical model′s accuracy.

    Fig.1 UH-60 level D train simulator

    Fig.2 Research simulator at The University of Liverpool

    We concentrate on mathematical model of dynamics, and overview dynamics modeling for the helicopter flight simulator.

    1 Model Overview

    In 1963, the equations of motion for a single-rotor helicopter were presented[5]. The aerodynamic equations were developed through a modified blade element approach, although other alternative techniques were considered. The equations were not constrained to one or a number of given flight conditions but valid for the entire flight regime including forward flight, hover, transition, autorotation, the effects of varying altitude, ground effects, and blade aeroelasticity in twist. They were suitable for analog, digital, and hybrid simulation programs. The equations included dynamic modeling of the main rotor and airframe to provide a mathematical model which will accept control inputs and calculate the resulting vehicle dynamics and orientation. The mathematical model for a single-rotor helicopter[5]was modified and implemented on the Langley Research Program Center real-time simulation (RTS) system[6]. The simulation program, with modifications by using a second-order Adams-Bashforth predictor integration formula with a solution rate of 32 iteration/s to solve the state variables and a modified secant method for trim condition calculation were used in support of two piloted studies of the Sikorsky S-61 helicopter, which was a five-blade, single-rotor, commercial passenger-type helicopter used by New York Airways. It was the early version of helicopter simulator, whose mathematical model was quietly simple when compared to further models.

    ARMCOP model[7]is a mathematical model of a single main rotor helicopter for piloted simulation, which was developed by Ames Research Center. The model has ten degrees of freedom: Six rigid-body, three rotor-flapping, and rotor rotational degrees of freedom. The rotor model assumed rigid blades with rotor forces and moments radially integrated and summed about the azimuth. The fuselage aerodynamic model use a detailed representation over a nominal angle of attack and sideslip range of ±15°, and it used a simplified curve fit at large angles of attack or sideslip. Stabilizing surface aerodynamics were modeled with a lift curve slope between stall limits and a general curve fit for large angles of attack. A generalized stability and control for a simplified engine/governor model, atmospheric turbulence, and a linearized six-degree-of-freedom dynamic model for stability and control analysis. Then for a specific helicopter—UH-60, some revisions were conducted to the ARMCOP model, with the following objectives[8]: (1) Improvement of the fidelity of the UH-60 fuselage aerodynamic model over a wide range of angles of attack and sideslip angles; (2) Modification of the tail rotor aerodynamic model to include the option of canting the tail rotor and modeling its associated aerodynamic effects; (3) Incorporation in the model of the control system for the UH-60 horizontal stabilator with variable incidence and the resultant aerodynamic effects; (4) Incorporation of the UH-60′s pitch bias actuator as part of the stability and control augmentation system.

    Sikorsky General Helicopter Flight Dynamics Simulation, known as GENHEL, provided the army with an engineering simulation for performance and handling quality evaluations. Initially, it was applied to analysis mode, which would eventually be used in the real time pilot-in-the-loop simulation. The model of BLACK HAWK helicopter was a total system, free flight representation based on GENHEL[9], and expanded/improved at Ames[10]. The basic model was a total force, non-linear, large angle representation in six rigid body degrees of freedom. The rotor aerodynamics were developed using a blade element approach. The total rotor forces and moments were developed from a combination of the aerodynamic, mass and inertia loads acting on each simulated blade. In the latter case, which represented the air mass degree of freedom, a uniform downwash was derived from momentum considerations, passed through a first order lag, and then distributed first harmonically as a function of rotor wake skew angle and the aerodynamic hub moment. The fuselage was defined by six component aerodynamic characteristics which were loaded from wind tunnel data which had been extended analytically to large angles. The aerodynamics of the empennage were treated separately from the forward airframe and used the look-up table method. The tail rotor was represented by the linearized closed form Bailey theory solution.

    FLYRT model[11]is a model developed by the McDonnell Douglas Helicopter System′s Fly Real Time for the AH-64A. Unlike GENHEL model ,the FLYRT model used a map-type approach to model its main rotor. This modeling approach was originally conceived in the late 1970s′ to allow real-time operation with the limited computational capabilities available at that time[12]. The problem was that this technique complete depended on and was restricted by this pregenerated map, especially since the code for the generation of more maps was not provided as part of FLYRT. The restriction was deemed surmountable since the computational power of computers enabled the model to run blade-element models in real-time without the need for the generation of a map. A blade-element-type main-rotor module was therefore developed to replace the map-type module in FLYRT. The rotor module held the access to the actual physical parameters of the rotor, so as introduce corrections, enhancements, and new theories as they were developed. The new model was known as the Blade-Element Model for APache (BEMAP)[13]. The basic structure of the new rotor module is based on the main-rotor module of GENHEL. BEMAP was, therefore, a version of FLYRT which were: (1) Equipped with a blade-element main-rotor module; (2) Upgraded with new trim and modified Aquations-of-motion modules; (3) Restructured to improve modularity; (4) Enhanced with flight-data-control-input access capability for flight-data comparisons; and (5) Updated with a versatile plotting option.

    F-B412 is the FLIGHTLAB simulation model of the Bell 412, developed by University of Liverpool[14], features a rigid, articulated, blade- element rotor formulation with flap and lag degrees of freedom. A quasi-steady aerodynamic component was used for the airload calculation on the blade elements. The blade was divided into several segments and each segment was treated as a two-dimensional panel producing aerodynamic loads as nonlinear functions of angle of attack, dynamic pressure and Mach number. The 2D aerofoil table should include lift, drag and pitch moment coefficients as functions of angle of attack and Mach number. There were several modelling options available within FLIGHTLAB for the fuselage aerodynamics, including a panel method and a simple table look-up. For the Bell-412 model, the table look-up option was chosen where the fuselage coefficients were supplied by means of look-up tables as functions of angle of attack and sideslip angle. Empennage aerodynamics tail rotor used the Bailey rotor model. The inflow model was the enhanced Peters-He finite state dynamic inflow model[15]. The F-B412 model produced very good on-axis responses in comparison to the NRC supplied flight test data over a broad speed range (15-120kt). However, off axis responses were mostly unsatisfactory[16].

    The UTIAS helicopter models were loosely based on that in Ref.[9] combined with first principals derivations all programmed in The MathWorks′ Simulink?environment developed by University of Toronto[17-18]. The main rotor models employed a blade-element formulation, The model used the Pitt-Peters dynamic inflow model for determining the total airflow, augmented with dynamic wake distortion which could be disabled for comparisons. The rotor rotational speed was determined using a second-order engine dynamics model. The tail rotor loads were determined using a Bailey rotor model, which contained a number of assumptions such as uniform inflow and a linear lift curve. Aerodynamic forces and moments for the fuselage and empennage were determined using look-up tables. Table 1 represents the differences between several real-time helicopter dynamic models for flight simulators.

    Table 1 Model comparison

    In the early ages, the dynamic model was pure mathematical model representing helicopter′s characteristic. Then GENHEL model stood out as the first engineering model with the look-up table method for some component with flight test data and wind tunnel data. It notably has improved model′s fidelity. It has worked as the baseline model for many other simulator flight dynamic models, such as BEMAP, F-B412, UTIAS model, and also used in a low-cost engineering real-time simulation facility MDRSF at Penn State University[19]. Uniform inflow model in GENHEL tirst developed by Paul and Bernard[20], introduced a term representing the acceleration of a mass of air surrounding the rotor into momentum-disk-type model to explain observations of first-order dynamics of the induced velocity below a hovering rotor. With the advancement in rotor wake study, some new inflow models developed like Pitt-Peters dynamic inflow model[21-22], which used in UTIAS model, and Peters-He finite state dynamic inflow model,which used in F-B412 model. The F-B412 model has been used in HELIFLIGHT simulator at the University of Liverpool[2]. The UTIAS model has been used in UTIAS Flight Research Simulator at the University of Toronto[23]. Helicopter simulator is not only used in industry and commerce, but now also used by universities, who have endeavoured to establish their own research simulator for flight control design, advanced rotorcraft flight dynamics modeling and handling quality analysis.

    2 Off-Axis Response Problem

    An ongoing concern with helicopter simulations is that they often exhibit an incorrect off-axis response to cyclic control inputs when compared with the corresponding flight test data[24-29]. The off-axis response problem is that in a simulator the initial lateral response to a longitudinal cyclic input, and vice versa, notably differ from that observed in the actual aircraft, or even in opposite. To obtain good flying and handling qualities, the off-axis problem must be solved. In fact, FAA Level D rotorcraft flight simulator certification requires that the simulated off-axis response be similar to that of the real rotorcraft[3].

    To understand this perplexing off-axis phenomena, there have been many attempts to explain the primary source of this problem. Most of the theories and methods can be classified as: (1)Flexible main rotor blades[17,30-44]; (2) Augmented Pitt-Peters dynamic inflow model and augmented Peters-He finite state inflow model[16,35-41]; (3) Free wake model[31,42-47].

    2.1 Flexible blade model

    A recent, non-real-time modeling effort by Theodore and Celi[30]has shown that blade flexibility and main rotor wake dynamics can be an important contributor to off-axis handling with a sophisticated non-real-time flight dynamics model. The employed wake model was the Bagai-Leishman free wake model[31], which was capable of modeling the wake geometry changes caused by maneuvers while making no assumptions of the geometry of the wake. The wake was discretized into a number of vortex segments shed from the rotor, which was solved iteratively at each time step and used to determine the local induced velocity at the rotor blade elements. The total number of vortex segments to be tracked was determined by the discretization resolution and the number of rotor revolutions retained in the analysis.

    In the Theodore and Celi reference, a coupled flap-lag-torsion finite element formulation was used to model the blade[32], combined with a modal coordinate transformation to reduce the number of degrees of freedom by retaining the first seven elastic modes. This approach of using a detailed finite element model to predict the flexible modes and then including only the important modes has been used in other applications as well[33], where the full finite element model was not practical to run in real-time. Run times were not given, however the model was not intended to be solved in real-time, and inclusion of the free wake model was described as requiring greater than ten times the overall CPU time compared to the same simulation with a linear inflow model. In addition, a recent paper using this model[34]gave representative run times on recently updated hardware, requiring 12—15 to simulate 2 min of flight time.

    In Ref.[14], using a Ritz expansion approach with constrained elastic modes was computationally compact and efficient and was therefore suitable for use in a real-time simulation. In Fig.3, the effect of including blade flexibility on the dynamic response, and in particular the off-axis response, of two UTIAS helicopter models is examined.

    Fig.3 Bell 206 lateral off-axis response[17]

    2.2 Pitt-Peters dynamic inflow model and Peters-He finite state inflow model

    Krothapalli, et al.[35-40]studied the rotor wake distortion effect during helicopter maneuvering flight by using various analytical tools, and finally obtained a generalized dynamic wake model with quasi-steady wake curvature augmentation in hover. Using the extended generalized wake model coupled with a rigid blade flapping model, the simulated off-axis flapping response to both longitudinal and lateral cyclic control stick inputs have shown fairly well qualitative results, when compared with the Sikorsky bearingless main rotor (SBMR) wind tunnel test data.

    He[41]adopted a full span prescribed vortex wake model consisting of discrete trailing and shed vortex elements to model the rotor wake geometric distortion effect during a pitch and/or roll rate. This rotor geometric wake distortion effect was implemented in a comprehensive rotorcraft simulation program, FLIGHTLAB, and simulated control response of a UH-60 Black Hawk helicopter to a lateral cyclic control stick doublet input in hover was compared with the corresponding flight test data. The comparison showed that the correlation of the off-axis pitch rate response prediction to a lateral cyclic control stick doublet input with flight test data was improved, although the magnitude was under predicted. In Ref.[16], the enhanced Peters-He finite state dynamic inflow model was integrated into the F-B412 model in FLIGHTLAB. It has been found that the pitch/roll off-axis responses in manoeuvres could be improved by including the wake geometry distortion effect in the Peters-He finite state dynamic inflow model, as Fig.4 shows.

    Fig.4 Off-axis response with and without wake distortion effet

    2.3 Free wake model

    Although the dynamic inflow model of Peters et al.[15]has found wide usage, there are known limitations. There are certain flight conditions in which the assumptions of the finite state inflow break down, for example, in certain descending flight conditions or vortex ring state. On the other hand, free-vortex wake models are inherently capable of modeling all important wake distortion effects and well suited for modeling wake-empennage interactions. The basic assumptions of the free-vortex models are perfectly valid for all generalized flight maneuvers, including descending flight conditions and large amplitude maneuvers. It is computationally expensive and therefore, not suitable for real-time applications. However, with advances in computing power, efficient free-vortex algorithms, and parallel processing, it is now feasible to use free wakes in real-time simulation.

    The first major contribution to the understanding of the off-axis response correlation discrepancy problem came from Rosen and Isser[43-44]. Rosen and his colleagues used a complex dynamic model of the rotor wake to show that a geometric distortion of the rotor wake structure occurred under a steady rotor shaft pitch and/or roll motion. This distorted rotor wake geometry induced a lateral inflow gradient perturbation across the rotor disk for a roll motion and a longitudinal inflow gradient perturbation for a pitch motion, which was counter to the kinematic inflow perturbation due to the shaft pitch and/or roll rate. Therefore, they proposed that this inflow change caused by the wake geometric distortion effect could yield a sign reversal in the off-axis flapping response in hover condition. Through including the wake geometric distortion effect during a pitch and/or roll maneuver in their unsteady aerodynamic model (TEMURA), a signi-ficant improvement in the off-axis response correlation with a UH-60 Black Hawk helicopter flight test identification results was achieved.

    Bagai et al. and Bhagwat et al.[31,45-47]used a free vortex method for the rotor wake to study the aerodynamics of a helicopter undergoing a steady pitch and/or roll maneuver. Through solving the rotor wake geometry under a prescribed pitch and/or roll rate using either a five-point central difference algorithm or a time accurate free vortex algorithm[46-47], it has been shown that the rotor pitch and/or roll maneuver was an additional source of distortion to the rotor wake geometry in both hover and forward flight. However, the wake distortion effect strongly depended on flight condition. In hover and low speed forward flight, the distorted rotor wake geometry caused by the pitch and/or roll maneuver induced an inflow gradient across the rotor disk, which was proportional to the rotor pitch and/or roll rate and manifests as a counter effect to the kinematic inflow across the rotor disk. At higher advance ratios, the effect of maneuver induced wake distortion is still existed, but was reduced by about half of that in hover condition.

    Theodore and Celi[48]investigated the effect of rotor wake modeling on the prediction of the off-axis response to pilot cyclic control stick input by incorporating a free wake model for helicopters with both a hingeless and an articulated rotor configuration and compared the simulation results with available flight test data. Although some improvements on the initial off-axis response prediction have been achieved, their results suggesed that deficiencies still existed in the wake model, and that other aerodynamic mechanism associated with the rotor wake dynamics might play a role in affecting the off-axis dynamic response correlation. They suggested that these mechanisms potentially might include the effect of unsteady aerodynamics and the transient effect associated with the development of the wake itself.

    Wachspress and Quackenbush[42]demonstrated the potential for using free wake models in real-time simulation. Then Rosen and Isser[43]showed the CHARM free wake model, integrated with the GENHEL simulation model, finally used in low-cost real-time simulator[19]. The CHARM free-vortex wake model was integrated with the GENHEL flight dynamics simulation of the UH-60A helicopter[49]. Fig.5 shows the schematic of GENHEL and CHARM coupling procedure.

    Fig.5 Schematic of GENHEL/CHARM coupling procedure[48]

    A high fidelity wake model was defined by increasing the spatial and temporal resolution of the wake until a converged response was observed, but this baseline model could not execute in real-time. A parametric study was performed to find the best combination of wake parameters to achieve real-time execution with minimal deviation from the baseline model in terms of the frequency and time responses in the pitch and roll axes. A real-time model was found and showed reasonable agreement with the baseline model as compared to a simple finite-state inflow model (Fig.6).

    Fig.6 Lateral off-axis response of real-time wake model[48]

    A parallel implementation of the free wake model was also investigated. The dynamic model with free wake model for flight simulator cause increased fidelity for real-time simulation and the use of first-principles/physics-based models in real-time flight simulation. Thus, the flight dynamics models could be used more as a predictive tool and not rely as much on empirical models.

    3 Conclusions

    The dynamic model for helicopter flight simulator has becone increasingly complex, from the simple mathematical model with uniform inflow, like ARMCOP model, to component-type complex model with flexible blades and free wake model, like UTIAS model which has improved the off-axis response situation.

    In future, with advances in computing power, efficient algorithms, parallel processing and developing helicopter dynamics theory, more and more complex dynamic model even coupled with computational fluid dynamics (CFD) calculating for rotor wake and aerodynamic interference, full finite element blade model will be developed for flight simulator.

    For training simulator, we will focus on helicopter special situation model such as hazard airflow[1], rortex ring state(VRS)[50], icing[51]and so on. The special situation model coupled with dynamic model will be the next key development.

    Acknowledgements

    This work was supported by the Fundamental Research Funds for the Central Universities; and the Research Innovation Program for College Graduates of Jiangsu Province (No.KYLX_0294).

    [1] ARAGON C R, LONG K R. Airflow hazard visualization for helicopter pilots: Flight simulation study results[C]//American Helicopter Society 61st Annual Forum. Grapevine, TX: AHS, 2005:61(1): 1.

    [2] PADFIELD G D, WHITE M D. Flight simulation in academia-HELIFLIGHT in its first year of operation at the University of Liverpool[J]. Aeronautical Journal, 2003, 107(1075):529-538.

    [3] FAA AC 120-63, Helicopter simulator qualification [S]. [S.l.]: Federal Aviation Administration, 1994.

    [4] JAR STD 1H. Helicopter flight Simulators[S]. The Netherlands: Joint Aviation Authorities, Hoofddorp, 2001.

    [5] TOLER J R, MCINTYRE W. Simulation of helicopter and v/stol aircraft. Volume I: Helicopter analysis report: NAVTRADEVCEN 1205-1[R]. 1963.

    [6] HOUCK J A, GIBSON L H, STEINMETZ G G. A real-time digital computer program for the simulation of a single rotor helicopter:NASA TM X-2874[R]. 1974.

    [7] TALBOT P D, TINLING B E, DECKER W A, et al. A mathematical model of a single main rotor helicopter for piloted simulation:NASA TM-84281[R]. 1982.

    [8] HILBERT K B. A mathematical model of the UH-60 helicopter:NASA TM-85890 [R]. 1984.

    [9] HOWLETT J J. UH-60A black hawk engineering simulation program. Volume 1: Mathematical model:NASA CR-166309[R]. 1981.

    [10]BALLIN M G. Validation of a real-time engineering simulation of the UH-60A helicopter:NASA TM-88360[R]. 1987.

    [11]HARDING J W, BASS S M. Validation of a flight simulation model of the AH-64 apache attack helicopter against flight test data[C]// American Helicopter Society 46th Annual Forum. Washington, D C: AHS,1990: 1309-1325.

    [12] HARRISON J M. An integrated approach to effective analytical support of helicopter design and development[C]// 6th European Rotorcraft and powered Lift Aircraft Forum. Bristol, England: Bristol University, 1980: 1-19.

    [13]MANSUR M H. Development and validation of a blade-element mathematical model for the AH-64A Apache helicopter:NASA TM-108863[R]. 1995.

    [14]GUBBELS A W, CARIGNAN S, ELLIS D K. The NRC bell 412 advanced systems research aircraft—A new facility for airborne simulation[J]. Canadian Aeronautics and Space Journal, 2000, 46(2):106-115.

    [15]PETERS D A, HE C J. Finite state induced flow models. II-three-dimensional rotor disk[J]. Journal of Aircraft, 1995, 32(2):323-333.

    [16]MANIMALA B, WALKER D, PADFIELD G. Rotorcraft simulation modelling and validation for control design and load prediction[C]//31st European Rotorcraft Forum. Florence, Italy: AHS,2005: 13-15.

    [17]HAYCOCK B, GRANT P. A real-time helicopter model with flexible main rotor blades[C]//AIAA Modeling and Simulation Technologies Conference. Portland, Oregon: AIAA, 2011: 6278.

    [18]HAYCOCK B, GRANT P. Dynamic wake distortion in the utias real-time helicopter models[C]//AIAA Modeling and Simulation Technologies Conference. Toronto, Ontario, Canada: AIAA, 2010: 8212.

    [19]HORN J F, BRIDGES D O, LOPES L V, et al. Development of a low-cost, multi-disciplinary rotorcraft simulation facility[J]. Journal of Aerospace Computing, Information, and Communication, 2005, 2(7):267-284.

    [20]CARPENTER P J, FRIDOVICH B. Effect of a rapid blade-pitch increase on the thrust and induced-velocity response of a full-scale helicopter rotor:NACA-TN-3044[R]. 1953.

    [21]GAONKAR G H, PETERS D. Review of dynamic inflow modeling for rotorcraft flight dynamics[J]. Vertica, 1988, 12(3):213-242.

    [22]CHEN R T N. A survey of nonuniform inflow models for rotorcraft flight dynamics and control applications:NASA-TM-102219[R]. 1989.

    [23]SHERARD S, REID L. The effects of simulator motion on handling qualities[C]// AIAA Modeling and Simulation Technologies Conference and Exhibit. Austin, Texas: AIAA,2003:5677.

    [24]ESHOW M M, ORLANDI D, BONAITA G, et al. Results of an A109 simulation validation and handling qualities study:NASA-TM-101062[R]. 1990.

    [25]TAKAHASHI M D. A flight-dynamic helicopter mathematical model with a single flap-lag-torsion main rotor:NASA-TM-102267[R]. 1990.

    [26]HARDING J W, BASS S M. Validation of a flight simulation model of the AH-64 Apache attack helicopter against flight test data[C]//46th American Helicopter Society Annual Forum. Washington D.C.:AHS,1990: 1309-1325.

    [27]CHAIMOVICH M, ROSEN A, RAND O, et al. Investigation of the flight mechanics simulation of a hovering helicopter[C]//48th American Helicopter Society Annual Forum. Washington D.C.:AHS,1992: 1237-1256.

    [28]HE C, LEWIS W. A parametric study of real time mathematical modeling incorporating dynamic wake and elastic blades[C] //48th American Helicopter Society Annual Forum. Washington D.C.:AHS,1992: 1181-1196.

    [29]TISCHLER M B, DRISCOLL J T, CAUFFMAN M G, et al. Study of bearingless main rotor dynamics from frequency-response wind tunnel test data[C]//Aeromechanics Specialists Conference on Aerodynamics, Acoustics and Dynamics. San Francisco, CA: [s.n.],1994:6.

    [30]THEODORE C, CELI R. Helicopter flight dynamic simulation with refined aerodynamics and flexible blade modeling[J]. Journal of Aircraft, 2002, 39(4):577-586.

    [31]BAGAI A, LEISHMAN J G, PARK J. Aerodynamic analysis of a helicopter in steady maneuvering flight using a free-vortex rotor wake model[J]. Journal of the American Helicopter Society, 1999, 44(2):109-120.

    [32]TURNOUR S R, CELI R. Modeling of flexible rotor blades for helicopter flight dynamics applications[J]. Journal of the American Helicopter Society, 1996, 41(1):52-66.

    [33]SCHMIDT D K, RANEY D L. Modeling and simulation of flexible flight vehicles[J]. Journal of Guidance, Control, and Dynamics, 2001, 24(3):539-546.

    [34]ALFRED J S. Rotorcraft brownout mitigation through flight path optimization using a high fidelity rotorcraft simulation model[D].USA:University of Maryland, 2012.

    [35]BAROCELA E B, PETERS D A, KROTHAPALLI K R, et al. The effect of wake distortion on rotor inflow gradients and off-axis coupling[C]// 22nd Atmospheric Flight Mechanics Conference. New Orleans, Louisiana:AIAA, 1997:272-282.

    [36]BAROCELA E B.The effect of wake distortion on dynamic inflow for lifting rotors[D]. St. Louis, Missouri: Washington University, 1997.

    [37]KROTHAPALLI K R, PRASAD J V R, PETERS D A. Improved wake geometry model for a maneuvering rotor[C]// Atmospheric Flight Mechanics Conference of the American Institute of Aeronautics and Astronautics. San Diego, CA:AIAA,1996:29-31.

    [38]KROTHAPALLI K R, PRASAD J V R, PETERS D A. Development of a comprehensive wake theory for lifting rotors[C]// 21st Atmospheric Flight Mechanics Conference. San Diego, California: AIAA, 1996:767-772.

    [39]KROTHAPALLI K R. Helicpoter rotor dynamic inflow modeling for maneuvering flight[D]. Atlanta,Georgia: Georgia Institute of Technology, 1998.

    [40]KROTHAPALLI K R, PRASAD J V R, PETERS D A. Helicopter rotor dynamic inflow modeling for maneuvering flight[J]. Journal of the American Helicopter Society, 2001, 46(2):129-139.

    [41]HE C, LEE C S, CHEN W. Rotorcraft simulation model enhancement to support design, testing and operational analysis[J]. Journal of the American Helicopter Society, 2000, 45(4):284-292.

    [42]WACHSPRESS D A, QUACKENBUSH T R, BOSCHITSCH A H. First-principles free-vortex wake analysis for helicopters and tilt rotors[C]// 59th American Helicopter Society Annual Forum. Phoenix, Arizona:AHS, 2003, 59(2): 1763-1786.

    [43]ROSEN A, ISSER A. A new model of rotor dynamics during pitch and roll of a hovering helicopter[J]. Journal of the American Helicopter Society, 1995, 40(3):17-28.

    [44]ROSEN A, ISSER A. A model of the unsteady aerodynamics of a hovering helicopter rotor that includes variations of the wake geometry[J]. Journal of the American Helicopter Society, 1995, 40(3):6-16.

    [45]PARK J S, LEISHMAN J G. The investigation of unsteady aerodynamics on rotor wake effects in maneuvering flight[J]. Annual Forum Proceedings—American Helicopter Society, 1999, 55(1):467-480.

    [46]BHAGWAT M J, LEISHMAN J G. Rotor aerodynamics during maneuvering flight using a yime-accurate free-vortex wake[J]. Journal of the American Helicopter Society, 2003, 48(3):143-158.

    [47]BHAGWAT M J. Mathematical modeling of the transient dynamics of rotor wakes using a time-accurate free vortex method[D]. College Park, MD: University of Maryland at College Park, 2001.

    [48]THEODORE C, CELI R. Prediction of the off-axis response to cyclic pitch using a maneuvering free wake model[C]// 25th European Rotorcraft Forum. Rome, Italy:AHS, 1999.

    [49]HORN J F, BRIDGES D O, WACHSPRESS D A, et al. Implementation of a free-vortex wake model in real-time simulation of rotorcraft[J]. Journal of Aerospace Computing, Information, and Communication, 2006, 3(3):93-107.

    [50]CHEN C. Development of a simplified inflow model for a helicopter rotor in descent flight[D]. Atlanta, Georgia: Georgia Institute of Technology, 2006.

    [51]AVANZINI G, de MATTEIS G, CISTRIANI L, et al. Development of a simplified helicopter model for piloted training simulation[C]//AIAA Modeling and Simulation Technologies Conference. Portland, Oregon:AIAA, 2011: 6280.

    Prof. Gu Hongbin is the head of Center for Flight Simulation and Advanced Technology for Training at Nanjing University of Aeronautics and Astronautics (NUAA). He received his Ph.D Degree in Aerospace from in 2000. His current research interest is flight simulation and advanced technology for training, such as computer based training and virtual reality.

    Mr.Fu Jun was born in Zhejiang, in 1990. He is currently a Ph.D. candidate in Traffic and Transportation Engineering, College of Civil Aviation Engineering, NUAA. His research interests include helicopter dynamics modeling and atmosphere turbulence modeling.

    Mr.Hu Jinshuo was born in Shandong, in 1988. He is currently a Ph.D. candidate in Traffic and Transportation Engineering, College of Civil Aviation Engineering, NUAA. His research interests include flight control simulation, nonlinear control system design, and applications of advanced control techniques to air vehicles.

    (Executive Editor: Zhang Bei)

    V275 Document code:A Article ID:1005-1120(2016)06-0647-10

    *Corresponding author, E-mail address: ghb@nuaa.edu.cn. How to cite this article: Gu Hongbin, Fu Jun, Hu Jinshuo. Review on key techniques of dynamics modeling for helicopter simulator[J]. Trans. Nanjing Univ. Aero. Astro., 2016,33(6):647-656. http://dx.doi.org/10.16356/j.1005-1120.2016.06.647

    免费在线观看影片大全网站| 国产午夜精品久久久久久一区二区三区 | 少妇熟女aⅴ在线视频| 亚洲一区二区三区色噜噜| 少妇人妻一区二区三区视频| 一个人免费在线观看的高清视频| 久久久久九九精品影院| 99久久成人亚洲精品观看| 亚洲午夜理论影院| 亚洲 欧美 日韩 在线 免费| 男人舔女人下体高潮全视频| or卡值多少钱| 日韩有码中文字幕| 高潮久久久久久久久久久不卡| 成人永久免费在线观看视频| 波多野结衣高清无吗| 长腿黑丝高跟| 亚洲成av人片免费观看| 国产成人a区在线观看| 99久久九九国产精品国产免费| 亚洲真实伦在线观看| www.999成人在线观看| 欧美成人a在线观看| 狂野欧美激情性xxxx| 一区二区三区国产精品乱码| 亚洲av五月六月丁香网| 日韩成人在线观看一区二区三区| а√天堂www在线а√下载| 国产中年淑女户外野战色| 国产高清三级在线| 亚洲精品在线美女| 草草在线视频免费看| 久久久精品欧美日韩精品| 精品一区二区三区视频在线 | 丰满人妻一区二区三区视频av | 9191精品国产免费久久| 嫩草影视91久久| 欧美黄色片欧美黄色片| 一本久久中文字幕| 中文字幕人妻熟人妻熟丝袜美 | 日韩欧美精品免费久久 | 一个人看视频在线观看www免费 | www.熟女人妻精品国产| 韩国av一区二区三区四区| 欧美不卡视频在线免费观看| 日韩国内少妇激情av| 午夜日韩欧美国产| netflix在线观看网站| 九色国产91popny在线| 亚洲天堂国产精品一区在线| 又黄又粗又硬又大视频| 国产免费一级a男人的天堂| 国产成人系列免费观看| 19禁男女啪啪无遮挡网站| 99国产精品一区二区蜜桃av| 日本精品一区二区三区蜜桃| 天天躁日日操中文字幕| 桃色一区二区三区在线观看| 免费在线观看日本一区| 亚洲人成网站在线播放欧美日韩| 国产不卡一卡二| 19禁男女啪啪无遮挡网站| 欧美黑人巨大hd| 在线观看免费午夜福利视频| 最后的刺客免费高清国语| 久久久久久久午夜电影| 免费人成在线观看视频色| 亚洲第一电影网av| 成熟少妇高潮喷水视频| 又黄又粗又硬又大视频| 久久欧美精品欧美久久欧美| 精品国产三级普通话版| 欧美一区二区国产精品久久精品| 最后的刺客免费高清国语| 欧美丝袜亚洲另类 | 国产真实伦视频高清在线观看 | 91在线观看av| 国产一级毛片七仙女欲春2| 午夜影院日韩av| 丰满人妻一区二区三区视频av | 国产av一区在线观看免费| 成人一区二区视频在线观看| av在线蜜桃| 午夜a级毛片| 成人特级av手机在线观看| 亚洲黑人精品在线| 欧美日韩亚洲国产一区二区在线观看| 欧美高清成人免费视频www| 精品国产超薄肉色丝袜足j| 亚洲精品在线观看二区| 国产精品三级大全| 午夜福利视频1000在线观看| 真人一进一出gif抽搐免费| 亚洲国产精品成人综合色| 午夜精品一区二区三区免费看| 1024手机看黄色片| 丝袜美腿在线中文| 午夜亚洲福利在线播放| 国产黄色小视频在线观看| 国产精品爽爽va在线观看网站| 中文字幕人妻丝袜一区二区| 久久久久久久久大av| 亚洲精品在线观看二区| 人人妻,人人澡人人爽秒播| 国语自产精品视频在线第100页| 三级男女做爰猛烈吃奶摸视频| 麻豆久久精品国产亚洲av| 中文字幕人成人乱码亚洲影| 久久6这里有精品| 日本撒尿小便嘘嘘汇集6| 真人做人爱边吃奶动态| 亚洲成人久久性| www.色视频.com| 在线观看美女被高潮喷水网站 | 国产一区二区亚洲精品在线观看| 精品一区二区三区视频在线 | 久久欧美精品欧美久久欧美| 日本精品一区二区三区蜜桃| 最后的刺客免费高清国语| 色av中文字幕| 亚洲午夜理论影院| 不卡一级毛片| 久久中文看片网| 极品教师在线免费播放| 在线免费观看不下载黄p国产 | 亚洲专区中文字幕在线| 两个人的视频大全免费| 精品久久久久久久人妻蜜臀av| 蜜桃亚洲精品一区二区三区| 久久草成人影院| 真人一进一出gif抽搐免费| 真人一进一出gif抽搐免费| 久9热在线精品视频| 久久亚洲真实| 亚洲精品色激情综合| 午夜免费观看网址| 久久中文看片网| АⅤ资源中文在线天堂| 亚洲av电影在线进入| 一个人免费在线观看的高清视频| 最近视频中文字幕2019在线8| 99精品欧美一区二区三区四区| 国产三级黄色录像| 国产精品女同一区二区软件 | 性欧美人与动物交配| 麻豆一二三区av精品| 欧美一区二区精品小视频在线| 久久久久亚洲av毛片大全| 九九久久精品国产亚洲av麻豆| 88av欧美| 90打野战视频偷拍视频| 18禁黄网站禁片午夜丰满| 国产色婷婷99| 精品一区二区三区视频在线观看免费| 他把我摸到了高潮在线观看| 久久久精品欧美日韩精品| 国产麻豆成人av免费视频| 窝窝影院91人妻| 欧美日韩精品网址| 久久精品国产自在天天线| 在线免费观看的www视频| 天美传媒精品一区二区| 中亚洲国语对白在线视频| 18禁美女被吸乳视频| 九九在线视频观看精品| 特级一级黄色大片| 日韩大尺度精品在线看网址| 国内久久婷婷六月综合欲色啪| 亚洲最大成人手机在线| 美女cb高潮喷水在线观看| 在线观看午夜福利视频| 给我免费播放毛片高清在线观看| 男插女下体视频免费在线播放| 女警被强在线播放| 国产淫片久久久久久久久 | 日韩欧美三级三区| 国产免费av片在线观看野外av| 一二三四社区在线视频社区8| 天美传媒精品一区二区| 麻豆成人午夜福利视频| aaaaa片日本免费| 亚洲欧美激情综合另类| 51国产日韩欧美| 真人一进一出gif抽搐免费| 国产亚洲欧美98| 最近最新中文字幕大全电影3| 色综合欧美亚洲国产小说| 熟女人妻精品中文字幕| 国产日本99.免费观看| 国产精品 欧美亚洲| e午夜精品久久久久久久| 免费观看的影片在线观看| 国产色爽女视频免费观看| 亚洲专区国产一区二区| 久久香蕉国产精品| 别揉我奶头~嗯~啊~动态视频| 国产成人欧美在线观看| 中文字幕av在线有码专区| 日本免费一区二区三区高清不卡| 波多野结衣高清作品| 嫁个100分男人电影在线观看| 制服人妻中文乱码| 在线观看av片永久免费下载| 天堂动漫精品| 99国产精品一区二区蜜桃av| 久久人妻av系列| 亚洲中文字幕日韩| 哪里可以看免费的av片| bbb黄色大片| 欧美中文日本在线观看视频| 国产探花在线观看一区二区| 久久久久久九九精品二区国产| 国产成人系列免费观看| 国产黄色小视频在线观看| 色视频www国产| 啦啦啦韩国在线观看视频| 国产精品1区2区在线观看.| 99国产极品粉嫩在线观看| 亚洲精品粉嫩美女一区| 18禁黄网站禁片午夜丰满| 美女被艹到高潮喷水动态| 天天一区二区日本电影三级| 欧美+亚洲+日韩+国产| 在线免费观看的www视频| 欧美日韩一级在线毛片| 中文字幕高清在线视频| 亚洲专区国产一区二区| 免费高清视频大片| 无人区码免费观看不卡| 国产伦一二天堂av在线观看| 国产午夜精品久久久久久一区二区三区 | 99热这里只有是精品50| 免费人成在线观看视频色| av中文乱码字幕在线| 国产精品三级大全| 成人鲁丝片一二三区免费| 在线播放无遮挡| 又黄又粗又硬又大视频| 久久久久久久久大av| 老鸭窝网址在线观看| 国产成人啪精品午夜网站| 婷婷精品国产亚洲av| 国产精品香港三级国产av潘金莲| 午夜福利视频1000在线观看| 9191精品国产免费久久| 男女做爰动态图高潮gif福利片| 69av精品久久久久久| 99热这里只有精品一区| 欧美最新免费一区二区三区 | 欧美高清成人免费视频www| 久久精品人妻少妇| 男女午夜视频在线观看| 俄罗斯特黄特色一大片| 99久久成人亚洲精品观看| 老汉色∧v一级毛片| 日本精品一区二区三区蜜桃| 蜜桃亚洲精品一区二区三区| 国产精品亚洲一级av第二区| 精品国产美女av久久久久小说| 18禁裸乳无遮挡免费网站照片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 舔av片在线| 国产精品一区二区三区四区久久| 女警被强在线播放| 国产精品野战在线观看| 亚洲美女黄片视频| 最近在线观看免费完整版| av专区在线播放| 久久人人精品亚洲av| 日日干狠狠操夜夜爽| 国产精品 欧美亚洲| 精品人妻偷拍中文字幕| 亚洲无线在线观看| 欧美色欧美亚洲另类二区| 女警被强在线播放| 亚洲18禁久久av| 一进一出好大好爽视频| 久久99热这里只有精品18| 亚洲av电影在线进入| 日本 欧美在线| 青草久久国产| 久久草成人影院| 久久久久免费精品人妻一区二区| 久久精品91蜜桃| 久久精品91蜜桃| 悠悠久久av| 亚洲美女视频黄频| www.999成人在线观看| 欧美性感艳星| 高清在线国产一区| 亚洲美女黄片视频| 免费搜索国产男女视频| 可以在线观看毛片的网站| 亚洲国产精品sss在线观看| 99riav亚洲国产免费| 国产精品电影一区二区三区| 成年女人看的毛片在线观看| 成年女人永久免费观看视频| 欧美极品一区二区三区四区| 精品国产超薄肉色丝袜足j| 国产亚洲欧美在线一区二区| 亚洲精品乱码久久久v下载方式 | 在线观看66精品国产| 亚洲美女黄片视频| 俺也久久电影网| 波多野结衣巨乳人妻| 国产一区二区在线观看日韩 | 欧美黑人欧美精品刺激| 国产一区二区亚洲精品在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲成人久久性| 国产精品久久久久久久电影 | 婷婷精品国产亚洲av在线| 午夜免费成人在线视频| 亚洲欧美日韩无卡精品| 久久久久久久久大av| 男女下面进入的视频免费午夜| 久久精品人妻少妇| 欧美黄色片欧美黄色片| 欧美性猛交╳xxx乱大交人| 99在线视频只有这里精品首页| 99久国产av精品| 欧美日本视频| www日本黄色视频网| 久久久久九九精品影院| 波多野结衣巨乳人妻| av欧美777| 给我免费播放毛片高清在线观看| 国产熟女xx| 免费在线观看影片大全网站| 国产一区二区三区在线臀色熟女| 国产综合懂色| 成人国产一区最新在线观看| 欧美成人a在线观看| 精品久久久久久久久久免费视频| 精品久久久久久,| 99国产精品一区二区蜜桃av| xxx96com| 国产真人三级小视频在线观看| 久久精品91无色码中文字幕| 9191精品国产免费久久| 88av欧美| 欧美日韩亚洲国产一区二区在线观看| 村上凉子中文字幕在线| 51午夜福利影视在线观看| 99热这里只有精品一区| 欧美+亚洲+日韩+国产| 两个人视频免费观看高清| 成人鲁丝片一二三区免费| 欧美一级a爱片免费观看看| 欧美日本亚洲视频在线播放| 国产精品久久视频播放| 美女大奶头视频| 51午夜福利影视在线观看| 亚洲在线自拍视频| 不卡一级毛片| 欧美一区二区国产精品久久精品| 真实男女啪啪啪动态图| 日本与韩国留学比较| 免费观看的影片在线观看| 国产精品一区二区免费欧美| 精品国产亚洲在线| 欧美日韩乱码在线| av专区在线播放| 免费看a级黄色片| 啦啦啦免费观看视频1| 色播亚洲综合网| a级一级毛片免费在线观看| 午夜福利18| 极品教师在线免费播放| 国产成人a区在线观看| 国产伦一二天堂av在线观看| 亚洲熟妇熟女久久| svipshipincom国产片| 久久人人精品亚洲av| 亚洲成人久久爱视频| 国产精品98久久久久久宅男小说| 欧美性感艳星| 国产激情偷乱视频一区二区| 亚洲精品色激情综合| 一区福利在线观看| 精品无人区乱码1区二区| 老司机在亚洲福利影院| av中文乱码字幕在线| 一二三四社区在线视频社区8| av片东京热男人的天堂| 人妻夜夜爽99麻豆av| 欧美在线一区亚洲| 亚洲成av人片在线播放无| 岛国视频午夜一区免费看| 一卡2卡三卡四卡精品乱码亚洲| 99久久久亚洲精品蜜臀av| 听说在线观看完整版免费高清| 亚洲国产精品久久男人天堂| 国产精品综合久久久久久久免费| 中亚洲国语对白在线视频| 国产伦精品一区二区三区四那| 免费av毛片视频| 成年女人毛片免费观看观看9| 亚洲人成网站在线播放欧美日韩| 久久精品影院6| 国产精品久久电影中文字幕| 国产野战对白在线观看| 女警被强在线播放| 中文字幕av在线有码专区| 村上凉子中文字幕在线| 精品熟女少妇八av免费久了| 夜夜躁狠狠躁天天躁| 麻豆成人av在线观看| 成人性生交大片免费视频hd| 亚洲欧美日韩高清专用| 国产综合懂色| 99久久久亚洲精品蜜臀av| 国产高清视频在线观看网站| 最近最新免费中文字幕在线| 内地一区二区视频在线| 国产精品一区二区三区四区免费观看 | 成年女人毛片免费观看观看9| 看片在线看免费视频| 在线观看66精品国产| 丰满乱子伦码专区| 精品久久久久久久人妻蜜臀av| 高清日韩中文字幕在线| 国产单亲对白刺激| 婷婷精品国产亚洲av在线| 免费一级毛片在线播放高清视频| 亚洲最大成人中文| 成人精品一区二区免费| 国产三级中文精品| 亚洲av免费在线观看| 午夜精品在线福利| а√天堂www在线а√下载| 99热精品在线国产| 高清日韩中文字幕在线| 亚洲第一电影网av| 中文字幕人妻丝袜一区二区| 蜜桃久久精品国产亚洲av| 国产野战对白在线观看| aaaaa片日本免费| 亚洲 欧美 日韩 在线 免费| a在线观看视频网站| 在线十欧美十亚洲十日本专区| 国产成年人精品一区二区| 真人做人爱边吃奶动态| 中文字幕人妻丝袜一区二区| 国产高清videossex| 午夜影院日韩av| 国产精品av视频在线免费观看| 麻豆国产av国片精品| 日本熟妇午夜| 欧美日韩中文字幕国产精品一区二区三区| 三级毛片av免费| 人妻夜夜爽99麻豆av| 国内精品久久久久久久电影| 又紧又爽又黄一区二区| 日本 欧美在线| 成人特级黄色片久久久久久久| 最近最新免费中文字幕在线| 熟妇人妻久久中文字幕3abv| 日韩欧美三级三区| 久久人妻av系列| 国产精品香港三级国产av潘金莲| 欧美另类亚洲清纯唯美| 又粗又爽又猛毛片免费看| 日日干狠狠操夜夜爽| 免费看光身美女| 国产亚洲欧美在线一区二区| av专区在线播放| 国产毛片a区久久久久| 欧美日韩综合久久久久久 | 全区人妻精品视频| 日日干狠狠操夜夜爽| 亚洲成av人片在线播放无| 国产一区二区亚洲精品在线观看| 男女那种视频在线观看| 久久久久久久久大av| 久久久久久大精品| avwww免费| 欧美极品一区二区三区四区| 97超视频在线观看视频| xxx96com| 九九久久精品国产亚洲av麻豆| 蜜桃亚洲精品一区二区三区| 久久久久亚洲av毛片大全| 99在线视频只有这里精品首页| 色老头精品视频在线观看| 久久久久久国产a免费观看| 啦啦啦观看免费观看视频高清| 亚洲av美国av| 色综合欧美亚洲国产小说| 亚洲精品日韩av片在线观看 | 天堂av国产一区二区熟女人妻| 日本精品一区二区三区蜜桃| 亚洲在线自拍视频| 亚洲精品色激情综合| 在线观看舔阴道视频| 国产高清视频在线观看网站| 久久久精品大字幕| 高清日韩中文字幕在线| 真人做人爱边吃奶动态| 亚洲精品国产精品久久久不卡| 国产精品野战在线观看| 免费人成视频x8x8入口观看| 国内精品一区二区在线观看| 一个人免费在线观看的高清视频| 亚洲中文日韩欧美视频| 色播亚洲综合网| 欧美不卡视频在线免费观看| 色播亚洲综合网| 色老头精品视频在线观看| 亚洲欧美日韩高清专用| 精品一区二区三区人妻视频| 嫁个100分男人电影在线观看| 午夜福利在线观看免费完整高清在 | h日本视频在线播放| 精品久久久久久久末码| 亚洲最大成人中文| 久久久久亚洲av毛片大全| 99热6这里只有精品| 国产99白浆流出| 欧美激情在线99| 小蜜桃在线观看免费完整版高清| 淫秽高清视频在线观看| 成年女人毛片免费观看观看9| 精品一区二区三区视频在线观看免费| 欧美成人一区二区免费高清观看| 无限看片的www在线观看| 国产精品久久视频播放| 欧美不卡视频在线免费观看| 五月伊人婷婷丁香| 51午夜福利影视在线观看| 日韩国内少妇激情av| 免费搜索国产男女视频| 熟女少妇亚洲综合色aaa.| 亚洲中文字幕日韩| 天天添夜夜摸| 我要搜黄色片| 国产精华一区二区三区| 亚洲男人的天堂狠狠| 婷婷精品国产亚洲av| www.999成人在线观看| 亚洲一区高清亚洲精品| 国产欧美日韩精品一区二区| 国内少妇人妻偷人精品xxx网站| 听说在线观看完整版免费高清| 精品日产1卡2卡| 精品一区二区三区视频在线观看免费| 97人妻精品一区二区三区麻豆| 欧美黑人巨大hd| 啦啦啦韩国在线观看视频| 国产免费一级a男人的天堂| 亚洲熟妇中文字幕五十中出| 午夜老司机福利剧场| 制服人妻中文乱码| 丰满人妻一区二区三区视频av | 少妇人妻精品综合一区二区 | а√天堂www在线а√下载| 国产黄片美女视频| 我要搜黄色片| 91麻豆av在线| 午夜影院日韩av| 91麻豆av在线| 最新中文字幕久久久久| 国产精品一区二区免费欧美| 日韩欧美三级三区| 夜夜夜夜夜久久久久| 啦啦啦免费观看视频1| www国产在线视频色| 无限看片的www在线观看| 亚洲av第一区精品v没综合| 日韩国内少妇激情av| 免费看十八禁软件| 天堂影院成人在线观看| 国产欧美日韩精品亚洲av| 精品久久久久久久久久久久久| 欧美三级亚洲精品| 国产三级黄色录像| 欧美乱妇无乱码| 亚洲 欧美 日韩 在线 免费| 国产一区二区三区在线臀色熟女| 窝窝影院91人妻| 久久久国产精品麻豆| 高潮久久久久久久久久久不卡| 99精品欧美一区二区三区四区| 99久国产av精品| 日本五十路高清| 国产亚洲欧美98| 亚洲 欧美 日韩 在线 免费| 欧美日韩福利视频一区二区| 日韩人妻高清精品专区| 一区二区三区免费毛片| 日韩欧美 国产精品| 久久中文看片网| av片东京热男人的天堂| 国内揄拍国产精品人妻在线| 日韩欧美国产在线观看| 色哟哟哟哟哟哟| 国产亚洲av嫩草精品影院| 精华霜和精华液先用哪个| 亚洲天堂国产精品一区在线| 男人舔女人下体高潮全视频| 国内久久婷婷六月综合欲色啪| 免费看a级黄色片| 国产成人av教育| 欧美bdsm另类| 久久精品影院6| 国产亚洲av嫩草精品影院| 亚洲成人久久性| 久久精品国产99精品国产亚洲性色| 女人高潮潮喷娇喘18禁视频| 久久精品影院6| 中文在线观看免费www的网站| 婷婷精品国产亚洲av在线| 欧美日韩乱码在线| 久久久久久久久中文|