• <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

    亚洲国产中文字幕在线视频| 欧美日本视频| 久久久精品欧美日韩精品| 在线十欧美十亚洲十日本专区| 桃红色精品国产亚洲av| 欧美+亚洲+日韩+国产| 精品久久久久久久人妻蜜臀av| 国模一区二区三区四区视频 | 999久久久精品免费观看国产| 中文在线观看免费www的网站| 美女扒开内裤让男人捅视频| 国产单亲对白刺激| xxxwww97欧美| 老司机午夜十八禁免费视频| 99久久精品国产亚洲精品| 午夜激情福利司机影院| 十八禁人妻一区二区| 这个男人来自地球电影免费观看| 亚洲av第一区精品v没综合| 精品免费久久久久久久清纯| 香蕉丝袜av| 99国产精品一区二区蜜桃av| 18禁美女被吸乳视频| 免费观看人在逋| 亚洲欧美一区二区三区黑人| а√天堂www在线а√下载| 18禁美女被吸乳视频| av视频在线观看入口| 天天躁日日操中文字幕| 高清毛片免费观看视频网站| 一本久久中文字幕| 一本久久中文字幕| 黄频高清免费视频| 精品久久久久久久毛片微露脸| 国产一区二区在线av高清观看| 亚洲中文av在线| 免费看美女性在线毛片视频| 国产黄色小视频在线观看| 日本a在线网址| 一夜夜www| 成年女人看的毛片在线观看| 久久久久久九九精品二区国产| 久久精品人妻少妇| 国产精品影院久久| 国产精品香港三级国产av潘金莲| 村上凉子中文字幕在线| 国产精品久久久久久亚洲av鲁大| 午夜免费成人在线视频| 午夜福利18| 亚洲 国产 在线| 18禁观看日本| 国产三级中文精品| 色在线成人网| 日韩成人在线观看一区二区三区| 99久久精品热视频| 九色成人免费人妻av| 长腿黑丝高跟| 午夜亚洲福利在线播放| 色播亚洲综合网| a级毛片在线看网站| 亚洲人成电影免费在线| 国产亚洲欧美98| 国产激情久久老熟女| 级片在线观看| 国产 一区 欧美 日韩| 国产在线精品亚洲第一网站| 国产aⅴ精品一区二区三区波| 淫妇啪啪啪对白视频| 99久久精品一区二区三区| 久久久久久久久中文| 两性午夜刺激爽爽歪歪视频在线观看| 欧美精品啪啪一区二区三区| 久久久久国产精品人妻aⅴ院| 一级黄色大片毛片| 亚洲 欧美 日韩 在线 免费| 身体一侧抽搐| 嫩草影视91久久| 亚洲人成网站在线播放欧美日韩| 老司机在亚洲福利影院| 1000部很黄的大片| 国产三级中文精品| 久久国产乱子伦精品免费另类| 此物有八面人人有两片| 欧美在线黄色| 美女午夜性视频免费| 婷婷六月久久综合丁香| 国产91精品成人一区二区三区| 国产精品一区二区三区四区久久| 综合色av麻豆| 偷拍熟女少妇极品色| 91在线精品国自产拍蜜月 | 国产av不卡久久| 久久久国产精品麻豆| av在线天堂中文字幕| 啪啪无遮挡十八禁网站| 1024手机看黄色片| 99久国产av精品| 亚洲欧洲精品一区二区精品久久久| 免费观看精品视频网站| 亚洲av片天天在线观看| 国产亚洲精品av在线| 黑人巨大精品欧美一区二区mp4| 亚洲精品在线美女| 亚洲精品在线观看二区| 国产精品久久久久久精品电影| 天堂√8在线中文| 亚洲国产欧美人成| 亚洲九九香蕉| 久久天躁狠狠躁夜夜2o2o| 精品欧美国产一区二区三| 亚洲国产精品久久男人天堂| 精品国产乱子伦一区二区三区| 搡老熟女国产l中国老女人| 又黄又爽又免费观看的视频| 精品不卡国产一区二区三区| 搡老岳熟女国产| 国产精品日韩av在线免费观看| 亚洲国产欧美一区二区综合| 久久久色成人| 色av中文字幕| 丁香六月欧美| 国产精品国产高清国产av| 一区二区三区高清视频在线| 国产三级在线视频| 亚洲,欧美精品.| 亚洲av电影在线进入| 欧美黑人巨大hd| 国产精品美女特级片免费视频播放器 | 男人的好看免费观看在线视频| 可以在线观看的亚洲视频| 婷婷丁香在线五月| 中文在线观看免费www的网站| 久久久久久久久免费视频了| 中文资源天堂在线| 又紧又爽又黄一区二区| 国产精品综合久久久久久久免费| 欧美日韩黄片免| 国产精品久久久久久亚洲av鲁大| 精华霜和精华液先用哪个| 在线永久观看黄色视频| 床上黄色一级片| 国产伦人伦偷精品视频| 嫁个100分男人电影在线观看| 精品无人区乱码1区二区| 中文在线观看免费www的网站| 99re在线观看精品视频| 欧美黑人巨大hd| 亚洲欧洲精品一区二区精品久久久| 一级黄色大片毛片| 色综合站精品国产| 欧美一级毛片孕妇| 欧美日韩黄片免| 99精品欧美一区二区三区四区| 国产麻豆成人av免费视频| 特级一级黄色大片| 国产主播在线观看一区二区| 午夜激情欧美在线| 日本成人三级电影网站| 综合色av麻豆| 亚洲欧洲精品一区二区精品久久久| 99久久久亚洲精品蜜臀av| 成在线人永久免费视频| 欧美另类亚洲清纯唯美| 亚洲精品一区av在线观看| 亚洲成av人片在线播放无| 嫩草影院入口| 国产亚洲av嫩草精品影院| 国产成人aa在线观看| 老司机午夜福利在线观看视频| 性欧美人与动物交配| 国产激情欧美一区二区| 神马国产精品三级电影在线观看| 国产高清视频在线观看网站| 毛片女人毛片| 亚洲国产欧美网| 成人特级黄色片久久久久久久| 在线永久观看黄色视频| 久久中文字幕一级| 蜜桃久久精品国产亚洲av| 国产精品国产高清国产av| 中文字幕人成人乱码亚洲影| 男插女下体视频免费在线播放| 人妻丰满熟妇av一区二区三区| 我的老师免费观看完整版| 亚洲中文av在线| 黄色片一级片一级黄色片| 日韩大尺度精品在线看网址| 大型黄色视频在线免费观看| 夜夜躁狠狠躁天天躁| 欧美激情久久久久久爽电影| 国产精品久久久久久精品电影| 成人鲁丝片一二三区免费| 免费大片18禁| 成人特级av手机在线观看| 亚洲国产看品久久| 国产高潮美女av| 又黄又粗又硬又大视频| 级片在线观看| 性欧美人与动物交配| 日本黄色片子视频| 三级毛片av免费| 亚洲精品中文字幕一二三四区| 综合色av麻豆| 国产激情欧美一区二区| 亚洲aⅴ乱码一区二区在线播放| 成人午夜高清在线视频| 香蕉丝袜av| 在线观看美女被高潮喷水网站 | 2021天堂中文幕一二区在线观| 身体一侧抽搐| 999久久久精品免费观看国产| 国产精品一区二区精品视频观看| 熟妇人妻久久中文字幕3abv| 亚洲午夜理论影院| 在线观看舔阴道视频| 国产久久久一区二区三区| 两人在一起打扑克的视频| 亚洲国产色片| 精品欧美国产一区二区三| 午夜福利视频1000在线观看| 欧美成人性av电影在线观看| 岛国在线免费视频观看| 国产成人精品久久二区二区91| 黄色视频,在线免费观看| 成人特级av手机在线观看| 免费观看人在逋| 国产黄色小视频在线观看| 精品国产乱码久久久久久男人| 国产精品女同一区二区软件 | 亚洲自拍偷在线| 人人妻人人澡欧美一区二区| 在线视频色国产色| 国产精品一及| 深夜精品福利| 欧美一区二区国产精品久久精品| 国产91精品成人一区二区三区| 97超视频在线观看视频| 中文字幕人妻丝袜一区二区| 一二三四社区在线视频社区8| 亚洲av美国av| АⅤ资源中文在线天堂| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 国产黄a三级三级三级人| 国产精品日韩av在线免费观看| or卡值多少钱| 欧美在线一区亚洲| 欧美zozozo另类| 性色av乱码一区二区三区2| 日日干狠狠操夜夜爽| 搡老妇女老女人老熟妇| 在线观看午夜福利视频| 99在线视频只有这里精品首页| 欧美成人性av电影在线观看| 这个男人来自地球电影免费观看| 男插女下体视频免费在线播放| 丝袜人妻中文字幕| 亚洲色图 男人天堂 中文字幕| xxxwww97欧美| 欧美极品一区二区三区四区| 国产aⅴ精品一区二区三区波| 欧美成人性av电影在线观看| 国产高清有码在线观看视频| 91av网站免费观看| 夜夜爽天天搞| 久久性视频一级片| 久久久久性生活片| 午夜福利成人在线免费观看| 在线播放国产精品三级| 国产在线精品亚洲第一网站| 国产精品野战在线观看| 成人精品一区二区免费| 国产熟女xx| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩黄片免| 我的老师免费观看完整版| av在线蜜桃| 国产精品久久久av美女十八| 午夜a级毛片| 美女黄网站色视频| 午夜福利欧美成人| 久久久久久久久久黄片| 最新中文字幕久久久久 | cao死你这个sao货| 一边摸一边抽搐一进一小说| 男人和女人高潮做爰伦理| 久久久久性生活片| 天堂网av新在线| 国产免费男女视频| 动漫黄色视频在线观看| 久久精品人妻少妇| 97人妻精品一区二区三区麻豆| avwww免费| 麻豆国产97在线/欧美| 国产精品自产拍在线观看55亚洲| 操出白浆在线播放| 我要搜黄色片| 国产黄a三级三级三级人| 一夜夜www| 成年人黄色毛片网站| 日本撒尿小便嘘嘘汇集6| 欧美国产日韩亚洲一区| 久久久久久久久免费视频了| 国产1区2区3区精品| 免费无遮挡裸体视频| 久久久久久九九精品二区国产| 男人舔女人下体高潮全视频| 国产午夜精品久久久久久| 精品免费久久久久久久清纯| 亚洲欧美日韩高清专用| 成人亚洲精品av一区二区| 91老司机精品| 国产亚洲精品一区二区www| 免费av不卡在线播放| ponron亚洲| 99久久精品热视频| 久久久国产精品麻豆| av视频在线观看入口| 亚洲五月天丁香| 久久久久性生活片| 亚洲 欧美一区二区三区| 国产精品99久久99久久久不卡| 免费看a级黄色片| АⅤ资源中文在线天堂| 麻豆av在线久日| 欧美一级毛片孕妇| 一区二区三区国产精品乱码| 国产精品亚洲一级av第二区| 日韩欧美三级三区| 观看美女的网站| 一进一出抽搐动态| 怎么达到女性高潮| 成人特级黄色片久久久久久久| 巨乳人妻的诱惑在线观看| a级毛片a级免费在线| 18美女黄网站色大片免费观看| 成熟少妇高潮喷水视频| 亚洲国产日韩欧美精品在线观看 | av天堂中文字幕网| 国产三级在线视频| 久久久国产精品麻豆| 欧美乱妇无乱码| 高清在线国产一区| 国内精品一区二区在线观看| 2021天堂中文幕一二区在线观| 成年女人永久免费观看视频| a在线观看视频网站| www.熟女人妻精品国产| 精品国内亚洲2022精品成人| 亚洲午夜精品一区,二区,三区| 最好的美女福利视频网| 成年女人永久免费观看视频| 精品一区二区三区视频在线 | 99久国产av精品| 欧美乱色亚洲激情| 免费看a级黄色片| 极品教师在线免费播放| 99精品久久久久人妻精品| 免费观看精品视频网站| 午夜福利在线观看免费完整高清在 | 亚洲人成电影免费在线| 精品国内亚洲2022精品成人| 国产亚洲精品久久久久久毛片| 久久精品国产清高在天天线| 99久久无色码亚洲精品果冻| 少妇的丰满在线观看| 又黄又粗又硬又大视频| 亚洲美女黄片视频| 日韩欧美免费精品| 在线看三级毛片| 亚洲av成人一区二区三| 好男人电影高清在线观看| 变态另类丝袜制服| 成人性生交大片免费视频hd| 免费在线观看日本一区| 亚洲熟女毛片儿| 久久精品国产亚洲av香蕉五月| av视频在线观看入口| 成年免费大片在线观看| 国产av麻豆久久久久久久| 国产爱豆传媒在线观看| 国产精品日韩av在线免费观看| 国产成人av教育| 91在线观看av| 亚洲成人精品中文字幕电影| 男女之事视频高清在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲国产日韩欧美精品在线观看 | 又粗又爽又猛毛片免费看| 九九热线精品视视频播放| 精品一区二区三区四区五区乱码| 少妇的逼水好多| 国产爱豆传媒在线观看| 精品久久久久久久毛片微露脸| 欧美色欧美亚洲另类二区| 国产伦在线观看视频一区| 亚洲国产精品sss在线观看| 亚洲av第一区精品v没综合| 老司机福利观看| 老熟妇乱子伦视频在线观看| 亚洲av美国av| 欧美丝袜亚洲另类 | 久久久国产成人免费| 午夜久久久久精精品| 日本熟妇午夜| 久久久成人免费电影| 午夜福利在线观看吧| 久久久久久久久免费视频了| 日本免费一区二区三区高清不卡| 亚洲色图av天堂| 最新中文字幕久久久久 | 99久国产av精品| 欧美激情久久久久久爽电影| 国产伦精品一区二区三区视频9 | 午夜影院日韩av| 91在线精品国自产拍蜜月 | 99久国产av精品| 免费在线观看视频国产中文字幕亚洲| 精品国产三级普通话版| 午夜激情福利司机影院| 国产探花在线观看一区二区| 亚洲欧美日韩无卡精品| 一本久久中文字幕| 国产精品野战在线观看| 美女大奶头视频| 国产美女午夜福利| 成人特级黄色片久久久久久久| a级毛片在线看网站| 男女那种视频在线观看| 久久婷婷人人爽人人干人人爱| av女优亚洲男人天堂 | 日本一二三区视频观看| 日韩精品中文字幕看吧| 中文字幕熟女人妻在线| 啦啦啦韩国在线观看视频| 精品国产美女av久久久久小说| 老司机午夜十八禁免费视频| 亚洲中文字幕日韩| 亚洲国产欧美一区二区综合| 国产精品国产高清国产av| 十八禁网站免费在线| www.精华液| 最新中文字幕久久久久 | 欧美日本视频| 亚洲av免费在线观看| 欧美极品一区二区三区四区| 91av网一区二区| 国产亚洲av嫩草精品影院| 免费看美女性在线毛片视频| 国产激情久久老熟女| 免费看a级黄色片| 18禁黄网站禁片午夜丰满| 9191精品国产免费久久| 亚洲国产精品合色在线| 俄罗斯特黄特色一大片| 在线国产一区二区在线| 老汉色av国产亚洲站长工具| 国产人伦9x9x在线观看| 五月伊人婷婷丁香| 精品一区二区三区视频在线 | 免费在线观看日本一区| 99热只有精品国产| 国产97色在线日韩免费| 色哟哟哟哟哟哟| 久久性视频一级片| av福利片在线观看| 免费高清视频大片| www.熟女人妻精品国产| 国产不卡一卡二| 99久久综合精品五月天人人| 免费在线观看影片大全网站| 最近最新中文字幕大全免费视频| 欧美日韩国产亚洲二区| 伊人久久大香线蕉亚洲五| 一区二区三区激情视频| 国产一区二区在线av高清观看| 欧美极品一区二区三区四区| 夜夜躁狠狠躁天天躁| 亚洲午夜理论影院| 欧美黄色片欧美黄色片| 成年女人永久免费观看视频| 亚洲美女视频黄频| 嫩草影院入口| 国产毛片a区久久久久| 久久精品aⅴ一区二区三区四区| 久久中文字幕一级| 国产精品综合久久久久久久免费| 12—13女人毛片做爰片一| 久久婷婷人人爽人人干人人爱| 亚洲aⅴ乱码一区二区在线播放| 中文资源天堂在线| 亚洲午夜理论影院| 天天添夜夜摸| 午夜精品一区二区三区免费看| 亚洲av美国av| 一卡2卡三卡四卡精品乱码亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 麻豆国产av国片精品| 麻豆久久精品国产亚洲av| 亚洲av日韩精品久久久久久密| 九色国产91popny在线| 久久精品91蜜桃| 母亲3免费完整高清在线观看| 欧洲精品卡2卡3卡4卡5卡区| 午夜a级毛片| 毛片女人毛片| 精品福利观看| 亚洲无线观看免费| 91字幕亚洲| 日本黄色视频三级网站网址| 欧美一级a爱片免费观看看| 午夜免费成人在线视频| 色精品久久人妻99蜜桃| 国产毛片a区久久久久| www.999成人在线观看| 日本一本二区三区精品| 亚洲真实伦在线观看| 午夜福利视频1000在线观看| 亚洲在线自拍视频| 欧美极品一区二区三区四区| 亚洲最大成人中文| 在线观看一区二区三区| 天天躁日日操中文字幕| 国产伦精品一区二区三区四那| 日韩国内少妇激情av| 91麻豆精品激情在线观看国产| 校园春色视频在线观看| 精华霜和精华液先用哪个| 国产精品久久电影中文字幕| 亚洲 欧美 日韩 在线 免费| 成人av一区二区三区在线看| 色吧在线观看| 亚洲精品美女久久av网站| 亚洲一区二区三区不卡视频| 久久天躁狠狠躁夜夜2o2o| 国产精品一区二区三区四区免费观看 | 脱女人内裤的视频| 变态另类丝袜制服| 国产av麻豆久久久久久久| 51午夜福利影视在线观看| 国产精品野战在线观看| av视频在线观看入口| 亚洲欧美精品综合一区二区三区| 午夜福利在线观看吧| 两个人的视频大全免费| 免费观看的影片在线观看| 亚洲精品一区av在线观看| 伦理电影免费视频| 国产私拍福利视频在线观看| 成人一区二区视频在线观看| 亚洲av免费在线观看| 在线观看舔阴道视频| 午夜福利在线观看吧| 亚洲熟女毛片儿| xxxwww97欧美| 在线观看一区二区三区| www.精华液| 久久香蕉国产精品| 欧美最黄视频在线播放免费| 两个人视频免费观看高清| 国产极品精品免费视频能看的| 国产亚洲精品久久久com| 免费看a级黄色片| 亚洲avbb在线观看| 九九久久精品国产亚洲av麻豆 | 老司机深夜福利视频在线观看| 中国美女看黄片| 女同久久另类99精品国产91| 午夜激情欧美在线| avwww免费| 国产成人啪精品午夜网站| 综合色av麻豆| 成人一区二区视频在线观看| 一个人观看的视频www高清免费观看 | 免费无遮挡裸体视频| 国产精品自产拍在线观看55亚洲| 久久国产精品影院| 97超级碰碰碰精品色视频在线观看| 波多野结衣高清无吗| 99久久成人亚洲精品观看| 五月玫瑰六月丁香| 国产1区2区3区精品| 亚洲欧美日韩卡通动漫| 特级一级黄色大片| 男人的好看免费观看在线视频| 亚洲精品一区av在线观看| 精品一区二区三区四区五区乱码| 日本 av在线| 少妇的逼水好多| 国产一区二区激情短视频| 亚洲色图av天堂| 狂野欧美白嫩少妇大欣赏| 神马国产精品三级电影在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产精品日韩av在线免费观看| 窝窝影院91人妻| 男女视频在线观看网站免费| 97碰自拍视频| bbb黄色大片| 久久久成人免费电影| 亚洲午夜理论影院| 高清毛片免费观看视频网站| 久久国产精品影院| 别揉我奶头~嗯~啊~动态视频| 亚洲中文字幕日韩| 国产日本99.免费观看| 婷婷六月久久综合丁香| 国产在线精品亚洲第一网站| 国产黄色小视频在线观看| 麻豆成人av在线观看| 母亲3免费完整高清在线观看| 亚洲熟女毛片儿| 黄片大片在线免费观看| 少妇人妻一区二区三区视频| 国产精品久久电影中文字幕|