System and method for finite element simulation of helicopter turbulence

Abstract

The present invention provides a turbulence model that has been developed for blade-element helicopter simulation. This model uses an innovative temporal and geometrical distribution algorithm that preserves the statistical characteristics of the turbulence spectra over the rotor disc, while providing velocity components in real time to each of five blade-element stations along each of four blades, for a total of twenty blade-element stations. The simulator system includes a software implementation of flight dynamics that adheres to the guidelines for turbulence set forth in military specifications. One of the features of the present simulator system is that it applies simulated turbulence to the rotor blades of the helicopter, rather than to its center of gravity. The simulator system accurately models the rotor penetration into a gust field. It includes time correlation between the front and rear of the main rotor, as well as between the side forces felt at the center of gravity and at the tail rotor. It also includes features for added realism, such as patchy turbulence and vertical gusts in to which the rotor disc penetrates. These features are realized by a unique real time implementation of the turbulence filters. The new simulator system uses two arrays one on either side of the main rotor to record the turbulence field and to produce time-correlation from the front to the rear of the rotor disc. The use of Gaussian Interpolation between the two arrays maintains the statistical properties of the turbulence across the rotor disc.The present simulator system and method may be used in future and existing real-time helicopter simulations with minimal increase in computational workload.

Description

ORIGIN OF THE INVENTION

The invention described herein was made in the performance of work under a NASA contract and is subject to the provision of Section 305 of the National Aeronautics and Space Act of 1958, Public Law 85568 (72 Stat. 435; 42 U.S.C. 2457).

BACKGROUND OF THE INVENTION

1. Technical Field of the Invention

The present invention is generally directed at aircraft simulation systems. More specifically, this invention relates to a new system and method for creating turbulence, including progressive gust penetration by the rotor, for use in real-time piloted helicopter simulations.

2. Description of the Prior Art

The modern helicopter is no longer a vehicle used for simple missions where only its hovering capability is required. The helicopter, and rotorcraft in general, are increasingly faced with complex missions which push the aircraft to its design limits. One emerging requirement which challenges helicopters is stabilized flight through moderate or severe atmospheric turbulence to accomplish high workload mission tasks.

Rotorcraft, both civilian and military, now compete with many other forms of transportation and have often shown great reliability and productivity in mission performance. These rotorcraft may be required to perform missions such as hovering over, and landing on moving platforms in high turbulence conditions. In light of the dynamic flight conditions which are encountered routinely, it is essential to have the capability to accurately predict rotorcraft performance in a turbulent atmosphere.

It is therefore clear that there is a great and still unsatisfied need for real-time rotorcraft atmospheric turbulence models in flight simulators. One of the most common grievances that many pilots have is that while they do feel turbulence acting on the helicopter in the flight simulator, the motion is not similar to real atmospheric turbulence.

One of the most severe limitations of the current rotorcraft atmospheric turbulence models is the lack of understanding of the fundamental physics underlying the rotorcraft response to atmospheric turbulence, which is different from that of the fixed wing aircraft. There is a fundamental difference between the ways in which rotary wing aircraft, such as helicopters, and fixed wing aircraft, such as airplanes, experience atmospheric turbulence. The difference is primarily due to blade rotational velocity. Whereas a wing element of an airplane has pure translational motion as the airplane cuts through a turbulence field, a blade element of a helicopter has translational as well as rotational motion. In an airplane both the turbulence and the response to turbulence are stationary random processes which require conventional body-fixed sampling, By comparison, the turbulence as well as the response to turbulence in a helicopter are cyclostationary and require blade-fixed sampling.

The rotational velocity introduces appreciable spatial distribution of turbulence velocities over the rotor disc. For conventional helicopters operating at attitudes of 1000 feet or below, the rotational velocity effects are not negligible. Research efforts have been made to provide methods with an analytical basis to predict and simulate turbulence and response statistics. One such simulation method is based on representing turbulence sample functions in terms of multi-variable sinusoids with random phases and provides second-order statistics of covariances and spectral densities.

Another turbulence simulation method is described in Costello, M. F., "A Theory for the Analysis of Rotorcraft Operating in Atmospheric Turbulence," Proceedings of the 46th Annual National Forum of the American Helicopter Society, Washington, D. C., May 1990, pp. 1003-1015. This article describes turbulence over a rotor blade that is approximated by a series of radial shape functions. With this expansion, a stochastic state space model is formed where the system dynamics matrix and the control matrix are periodic with a period equal to the rotor rotational speed and input to the plant is generated by independent white noise sources.

The following articles reflect some of attempted methods for improving turbulence simulation models, all of which are incorporated herein by reference:

Howlett, J. J., "UH-60A Black Hawk Engineering Simulation Program: Vol. I -Mathematical Model," NASA CR-1 66309, Dec. 1981.

Prasad, J. V. R., Riaz, J., Gaonkar, G. H., Yingyi, D., "Real Time Implementation Aspects of a Rotorcraft Turbulence Simulation Method," 49th Annual Forum of the American Helicopter Society, St. Louis, Mo., May 1993, p. 459.

Dahl, H. J. and Faulkner, A. J., "Helicopter Simulation in Atmospheric Turbulence," Vertica, Vol.3, 1979, pp. 65-78.

Judd, M. and Newman, S. J., "An Analysis of Helicopter Rotor Response due to Gust and Turbulence," Vertica, Vol. 1, 1977, pp. 179188.

Dang, Y., G. Gaonkar, Prasad, J., Zhang, H., "Parallel Computing of Helicopter Response to Turbulence Toward Real-Time Implementation," 50th Annual Forum of the American Helicopter Society, Washington, D.C., Vol. 2, May 1994, pp. 869-882.

However, the foregoing conventional methods have not proven to be totally satisfactory in real-time simulations since they either take excessive execution time or they do not yield realistic turbulence according to pilot opinion.

SUMMARY OF THE INVENTION

The present invention provides a unique simulator system and method for use in future and existing real-time helicopter simulations with minimal increase in computational workload.

The present invention further provides a turbulence model developed for blade-element helicopter simulation. This model uses an innovative temporal and geometrical distribution algorithm that preserves the statistical characteristics of the turbulence spectra over the rotor disc, while providing velocity components in real time to each of five blade-element stations along each of four blades, for a total of twenty blade-element stations.

An initial investigation of the simulator system was performed using a piloted, motion-based simulation of the Sikorsky UH60A Black Hawk helicopter. Although only the vertical component of stochastic turbulence was used in this initial investigation, vertical turbulence components induced vehicle responses in all translational and rotational degrees of freedom of the helicopter.

The single-degree-of-freedom configuration of the present simulator system was compared to a conventional full 6-degrees-of-freedom baseline configuration, where translational velocity inputs are superimposed at the vehicle center of gravity, and rotational velocity inputs are created from filters that approximate the immersion rate into the turbulent field.

The simulator system includes a software implementation of flight dynamics that adheres to the guidelines for turbulence set forth in military specifications. One of the features of the present simulator system is that it applies simulated turbulence to the rotor blades of the helicopter, rather than to its center of gravity. The simulator system accurately models the rotor penetration into a gust field. It includes time correlation between the front and rear of the main rotor, as well as between the side forces felt at the center of gravity and at the tail rotor. It also includes features for added realism, such as patchy turbulence and vertical gusts into which the rotor disc penetrates.

Briefly, the foregoing and other previously unattainable features and advantages of the present invention are realized by a unique real time implementation of the turbulence filters. The new simulator system uses two arrays, one on either side of the main rotor, to record the turbulence field and to produce time-correlation from the front to the rear of the rotor disc. The use of "Gaussian Interpolation" between the two arrays maintains the statistical properties of the turbulence across the rotor disc.

BRIEF DESCRIPTION OF THE DRAWINGS

The above and other features of the present invention and the manner of attaining them will become apparent, and the invention itself will be best understood, by reference to the following description and the accompanying drawings, wherein:

FIG. 1 is an isometric view of a helicopter cockpit module for use in conjunction with the simulator system and method according to the present invention;

FIG. 2 is a top view of a simulated rotor disc forming part of the simulator system and method of the present invention, showing the rotor velocities and the added turbulence velocities;

FIGS. 3A, 3B, 3C and 3D represent a functional block diagram of the present simulator system titled "TURBR";

FIG. 4 represents a block diagram of the Patches of Turbulence circuit used in the present simulator system; and

FIG. 5 represents a block diagram of the Vertical Gusts of Turbulence circuit used in the present simulator system.

Similar numerals refer to similar elements in the drawing. It should be understood that the sizes of the different components in the Figures are not in exact proportion, and are shown for visual clarity and for the purpose of explanation.

DETAILED DESCRIPTION OF THE INVENTION

As used herein, turbulent fluid motion is an irregular condition of flow in which the various quantities show a random variation with time and space coordinates, so that statistical distinct average values can be discerned. According to this definition, the flow has to satisfy the condition of irregularity which highlights an important feature of turbulence. Because of the irregular nature of any turbulent flow, it is very difficult to describe it by deterministic functional relationships. On the other hand, turbulence velocities are generally random in the sense that it is possible to describe it by the laws of probability, hence a statistical approach is often adopted to describe atmospheric turbulence velocities.

The flight simulator used for performing the piloted simulation according to the present invention is NASA Ames Research Center's Vertical Motion Simulator (VMS). This simulator is a general purpose research tool that retains a high degree of flexibility to enable tailoring for a wide range of fixed and rotary wing applications. The simulator can be configured to meet the needs of a particular task by selecting hardware and software options. In a specific experiment related to the present invention, the simulator was configured for the purpose of evaluating the influence of stochastic turbulence created at the blade-element stations of a helicopter rotor simulator system.

Some key facilities of the simulator system are comprised of a motion simulator, a visual simulator, pilot controls, primary flight instrumentation, and a software or processor implementation of flight dynamics. These components are briefly described below.

The motion simulator system generally provides a cockpit module 10 with motion in all six axes (roll, pitch, yaw, heave, lateral and longitudinal), as shown in FIG. 1. The motion simulator system is capable of high accelerations, velocities and displacements, and notably, the performance described herein can be achieved simultaneously in all six degrees of freedom.

The visual simulator system generally provides a visual representation of the flight environment through monitors mounted in the cockpit module 10. The cockpit module 10, pilot controls and primary flight instrumentation are similar to those of a GENHEL UH60A Black Hawk helicopter. The software implementation of flight dynamics including all the requisite motion and visual communication software, has been executed on a single processor, with a cycle time of about 12 milliseconds.

The experimental simulations conducted at NASA Ames Research Center met with very favorable pilot opinion. Some of these simulations were implemented into the GENHEL UH60A mathematical model described in Howlett, J. J., "UH-60A Black Hawk Engineering Simulation Program: Vol. 1- Mathematical Model," NASA CR-166309, Dec. 1981. An heuristic model of turbulence was developed that preserves the Gaussian statistical characteristics of turbulence filters over the rotor disc. This model avoids the large computational and storage requirements of conventional investigations into rotorcraft turbulence models, such as those described in Prasad, J. V. R., Riaz, J., Gaonkar, G. H., Yingyi, D., "Real Time Implementation Aspects of a Rotorcraft Turbulence Simulation Method," 49th Annual Forum of the American Helicopter Society, St. Louis, Mo., May 1993, p. 459.

The present simulator system is designed to create and propagate all three translational components of turbulence in real time. The turbulence velocity components are superimposed at each blade element location as a function of the helicopter velocity, with distribution according to the blade element instantaneous geometry.

The entire stochastic rotor element turbulence model will be described herein. Examples from this model are presented for the helicopter in low speed flight (i.e., 10 knots). In this study the Dryden spectral form is selected due to its computational efficiency, although a discrete model of a curve-fit approximation to the von Karman form may alternatively be used. Spectral analysis of the simulation shows that the rotor-to-body attenuation is so severe that the original gust spectra appear to be of secondary consideration.

Conventional military specifications ("MIL SPEC") angular turbulence filters are a substitute for a finite-element distribution of translational turbulence in the air mass over the airframe. In the present simulator system only the rotor disc may be considered as an "airframe", and the translational turbulence in the air mass is distributed over the rotor disc to each blade element. This distribution produces rotor moments and forces. The MIL SPEC translational turbulence filters may be used in the present simulator model.

Referring to FIG. 2, it is a diagrammatic representation of the rotor geometry, illustrating a rotor disc 12 in top plan view. The helicopter horizontal velocity component in the forward direction is represented by V.sub.H, and the added turbulence velocity components are represented by .DELTA.u.sub.L, .DELTA.v.sub.L, .DELTA.w.sub.L, .DELTA.u.sub.R, .DELTA.v.sub.R, .DELTA.w.sub.R. Since the rotor disc 12 encounters turbulence at its leading edge 14, the line of turbulence which translates along the rotor disc 12 based on the velocity of the helicopter, is represented by the onset line 16. The onset line 16 is perpendicular to the vehicle forward flight path (i.e., the direction of V.sub.H,), and tangential to the leading edge of the rotor disc 12. The vertical added turbulence velocity components .DELTA.w.sub.L and .DELTA.w.sub.R and the horizontal added turbulence velocity components .DELTA.u.sub.L, .DELTA.v.sub.L, .DELTA.U.sub.R, and .DELTA.v.sub.R, which will be described later in greater detail, are represented by corresponding arrows at the left and right ends of the onset line 16.

From Taylor's hypothesis, the turbulent flow is stationary and homogeneous in the longitudinal direction. Turbulence filter outputs are statistically valid at either fixed spatial locations or at points translating with the helicopter, as is the case in the conventional body-fixed formulation. The outputs of the turbulence filters are applicable at points on the onset line 16 at the discrete time points. Interpolation perpendicular to the onset line 16 requires a transport delay that is proportional to the distance along the flight path to the blade element location, and inversely proportional to the helicopter aerodynamic velocity. The turbulence is assumed uncorrelated at each side of the rotor disc 12, and homogeneous in the lateral dimension. The statistical properties of the turbulence are then invariant over the entire rotor disc 12 by using "Gaussian Interpolation" between values applicable on each side of the rotor disc 12.

In a particular simulation using the Black Hawk helicopter model, the rotor system has four blades (N=4). For each of the blades the number of segments used in the simulation was five (M=5). The blade and segment indices are defined, as follows:

1<n<N (blades) 1<m<M (segments)

The finite element analysis is performed using in part the equal-annuli algorithm described in Howlett, J. J., "UH-60A Black Hawk Engineering Simulation Program: Vol. I -Mathematical Model," NASA CR-1 66309, Dec. 1981, to generate the radii (r.sub.m) to the blade stations. r.sub.m may be computed as follows: ##EQU1## where the rotor radius is R=268 ft, the hinge offset is e=1.2 ft, and the spar length is e'=2.2 ft.

The azimuth angle of the n.sup.th blade with respect to the helicopter aft centerline is given by the following equation: ##EQU2## where .OMEGA.=27 rad/sec is the nominal rotor RPM.

In the present model, rather than defining the translational disturbance velocities at the helicopter center of gravity, they are defined at two onset points 20, 21 located along the onset line 16. Using Gaussian inputs, three independent turbulence velocity histories are computed at each of these two points 20, 21 in real time. These histories are updated each cycle time, and the pertinent values from these histories are distributed to the individual segment locations by temporal and geometrical algorithms.

The helicopter horizontal velocity component in the forward direction (v.sub.H) may be defined approximately by the following equation, where u.sub.b is the helicopter longitudinal aerodynamic velocity and v.sub.b is the lateral aerodynamic velocity: ##EQU3##

Using the velocity v.sub.H, the sideslip angle ".beta." and the perpendicular onset line 16 shown in FIG. 2, the aerodynamic azimuth angle (wind axis) .PSI..sub.n for the n.sup.th blade would then be defined by the following equation:

.PSI..sub.n =.PSI..sub.n =.beta. where .PSI..sub.n is the n.sup.th blade geometrical azimuth angle from aft.

This geometry is sufficient to establish the longitudinal distribution. The left and right Dryden filters receive Gaussian inputs, so the outputs of these linear filters are also Gaussian. Furthermore, a transport delay in these outputs does not alter the statistical characteristics. The outputs of the filters are stored in tables, which are updated each cycle time. The value applicable for a given element is determined by selecting a value from the tables. Given that the table size is established (i.e., K.sub.M =500) for storing time histories of the outputs of the filters, the minimum aerodynamic velocity V.sub.min in that may be accommodated is determined by distributing the table over the rotor diameter: ##EQU4##

If a smaller minimum velocity is required, a larger table size must be selected. For the simulation cycle time .DELTA.t=0.012 sec, the minimum velocity for K.sub.M =500 is 8.944 ft/sec (5.3 knots). The rotor-plane aerodynamic velocity used in the turbulence filters is restricted to this minimum. ##EQU5##

For this simulation, independent of vehicle velocity, exactly K.sub.M .DELTA.t=6 seconds of data is stored in the tables. Time histories of the outputs of six turbulence filters are contained in six separate tables, each of length K.sub.M. At the minimum velocity the individual cells correspond to a length of 0.10732 ft (such that 500 cells of the tables span the rotor diameter), whereas at a velocity of 100 ft/sec (59.2 knots), for example, the cells are each 1.2 ft long such that only 45 cells are needed to span the rotor diameter.

In order to select the velocities to be used from the tables, for any blade element, an integer index in the table for the blade element turbulence (k.sub.m, n) is computed to determine the blade element perpendicular distance from the onset line. Defining,

r.sub.m,n =R+r.sub.m cos .PSI..sub.n the delay index (k.sub.m, n) is defined by the following equation: ##EQU6## where the upper brackets denote the integer ceiling operation. The tables extend beyond the rotor disc for higher vehicle velocities than the minimum, such that stationary outputs are available for table size increases caused by decreases in velocity.

The Dryden model will now be described in more detail. It is known that the low-altitude vertical turbulence scale length (L.sub.w) is given as a piecewise continuous function of altitude (h), as follows: ##EQU7## and the in-plane scale lengths L.sub.u and L.sub.v are given by the following functions: ##EQU8## where f.sub.uv is defined by:

f.sub.uv =0.177+0.000823h The horizontal turbulence RMS intensities .sigma..sub.u,.sigma..sub.v are a function of the selected vertical turbulence intensity .sigma..sub.w, .sigma..sub.u =.sigma..sub.v =.sigma..sub.w f.sub.uv.sup.-0.4

The results shown herein are for .sigma..sub.w =5 ft/sec, h=200 ft, and v.sub.H =16.9 ft/sec (10 knots).

Using Taylor's frozen field hypothesis, the Dryden form of the linear turbulence transfer functions are given by the following equations: ##EQU9##

The angular transfer functions are developed from partial differential equations relating these translational disturbances to the immersion rate of the aircraft into the turbulent field. These equations use a parameter "b", which is the wing span length of the helicopter. In an exemplary simulation "b" is approximately 26 ft. ##EQU10##

These angular disturbances are not used in the illustrated simulation model because the geometry is considered in the points of application of the translational turbulence excitations. The asymmetric turbulence velocities produce moments, which, in turn, produce angular activity. The angular disturbence equations are meant to approximate the body dynamics. The new turbulence implementation uses the body's own interactions with the rotor to produce angular disturbances.

Continuous transfer functions involving random inputs are typically converted to discrete form using the zero-order hold formulation, where the input is assumed constant over each computer cycle .DELTA.t. The z-transforms of Laplace functions f.sub.i (s) then become, ##EQU11## which may then be converted to difference equations for discrete implementation.

Using this technique on the Dryden transfer functions produces the following coefficients: ##EQU12##

These coefficients are used in the following six difference equations, each driven by an independent Gaussian noise source n.sub.i, with a zero mean value and a unity standard deviation.

.DELTA.u.sub.L (k)=f.sub.1 .DELTA.u.sub.L (k-1)+f.sub.2 n.sub.1 (k) .DELTA.u.sub.R (k)=f.sub.1 .DELTA.u.sub.R (k-1)+f.sub.2 n.sub.2 (k) .DELTA.v.sub.L (k)=g.sub.1 .DELTA.v.sub.L (k-1)+g.sub.2 .DELTA.v.sub.L (k-2) +g.sub.3 n.sub.3 (k)+g.sub.4 n.sub.3 (k-1) .DELTA.v.sub.R (k)=g.sub.1 .DELTA.v.sub.R (k-1)+g.sub.2 .DELTA.v.sub.R (k-2) +g.sub.3 n.sub.4 (k)+g.sub.4 n.sub.4 (k-1) .DELTA.w.sub.L (k)=h.sub.1 .DELTA.w.sub.L (k-1)+h.sub.2 .DELTA.w.sub.L (k-2) +h.sub.3 n.sub.5 (k)+h.sub.4 n.sub.5 (k-1) .DELTA.W.sub.R (k)=h.sub.1 .DELTA.W.sub.R (k-1)+h.sub.2 .DELTA.W.sub.R (k-2) +h.sub.3 n.sub.6 (k)+h.sub.4 n.sub.6 (k-1)

These computationally efficient equations produce stable outputs for all stable inputs. For the body-fixed formulation these velocities are superimposed at the center of gravity. However, in the present simulator system, these velocities are created at the two onset points 20, 21 and then distributed and superimposed at the various blade element locations. They produce forces and moments because they modify the angle of attack of each blade element.

The foregoing analytical study will now be correlated more specifically to the block diagrams in FIGS. 3A, 3B, 3C, with reference to the following two tables of variables.

TABLE I______________________________________TURBR VARIABLESName Description______________________________________BETA Sideslip AngleDISP Dispersion - standard deviation of vertical turbulenceDT2 Cycle TimeFIRSTELM Value of single input for all segmentsFTR Distance from CG to TR in feetGAINU Scale factor for U-body turbulenceGAINV Scale factor for V-body turbulenceHCG HeightIFIRSTELM Use single input for all segmentsINDX Index for segment arraysITRBGUST Use vertical gustsITRBPAT Use Patchy TurbulenceITURBRW Allow segment inputsITWBD Create vertical body turbulenceIVCG Index for CG for V-arrayIVTR Index for Tail Rotor for V-arrayNx Number of cells in the main rotor arraysPSIR(M) Angle to blade #1 in radiansR, RMR Main Rotor RadiusSEGMID(N) Midpoints of the segments (in fractions or RMR)UB X-velocity (ft/s)UPGMR(M,N) Segment turbulence componentsUTURB U - Body TurbulenceVB Y-velocity (ft/s)VTURB V - Body TurbulenceVYGTR Tail rotor Y-gust componentWTURB W - Body Turbulence______________________________________ TABLE II______________________________________EXTRA SECTIONS(VERTICAL GUSTS, PATCHES OF TURBULENCE)Gust-Name Patch-Name Description______________________________________TIME, A(303) TIME, A(303) Total time since beginning operateDT2, A(168) DT2, A(168) Cycle timeGTLGUS PTLSIG Time since last change of outputGTRNDW PTRNDW Time to wait until next change of outputGLMBD PLMBD Average time to wait between changesGLMBDMX PLMBDMX Max Average timeGLMBDMN PLMBDMN Min Average timeGRMPT PRMPT Time to ramp in the changeGRMPTMX PRMPTMX Max Ramp timeGRMPTMN PRMPTMN Min Ramp timeGVMIN PVMIN Minimum velocityGVMAX PVMAX Maximum velocityGUSTNEW PSIGNEW New output value to ramp towardGOUT POUT Current output value______________________________________

Given the height above ground (h) and vertical standard deviation (.sigma..sub.w), the remaining turbulence specifications can be defined as follows: The function f.sub.h appendix is given by the following equation:

f.sub.h =0.177+0.000823h TABLE III______________________________________TURBULENCE SPECIFICATIONS______________________________________ 37.82 h .ltoreq. 10 hf.sub.h.sup.-1.2 /2 10 < h .ltoreq. 1000L.sub.u, L.sub.v h/2 1000 < h .ltoreq. 175 875 h .gtoreq. 1750 5 h .ltoreq. 10L.sub.w h/2 10 < h .ltoreq. 175 875 h > 1750 (0.18523).sigma..sub.w h .ltoreq. 10.sigma..sub.u, .sigma..sub.v .sigma..sub.w .function..sub.h.sup.-0.4 10 < h < 1000 .sigma..sub.w h .gtoreq. 1000______________________________________

FIG. 3A shows 11 parameters that are input to the simulator system, and generally grouped into four different groups, as illustrated in the following table IV:

TABLE IV______________________________________PARAMETERS GROUPINGGroup Parameter______________________________________Rotor Geometry SEGMID(N), PSIR(M), R, FTR, BETAHelicopter Dynamics UB, VB, HCGSimulator System Variable DT2Turbulence Variable Determined NX, DISPby User______________________________________

FIGS. 3A, 3B, 3C and 3D combined represent a functional block diagram of the present simulator system herein titled "TURBR". For simplicity of description, the circuits shown in FIGS. 3A and 3B will be collectively referred to as the input circuit 100, the circuit shown in FIG. 3C will be referred to as the blade-element turbulence circuit 102, and the circuit shown in FIG. 3D will be referred to as the bodyturbulence circuit 104.

Referring to FIGS. 3A and 3B, the input circuit 100 generates the following output signals that are fed to the blade-element turbulence circuit 102 in FIG. 3C:

AA: corresponds to DT2. AB: corresponds to V.sub.HV, which is calculated as follows: v.sub.HV =(UB.sup.2 +VB.sup.2).sup.1/2 AC: corresponds to SEGMID(N). AD: corresponds to PSIR(M). AE: is calculated as follows: tan.sup.-1 (VB/UB) AF: corresponds to V.sub.DTI, which is calculated as follows: 1/v.sub.DT, where v.sub.DT =(DT2).(V.sub.HV) AG: corresponds to R.

The following output signals in FIG. 3B are fed to the body-turbulence circuit 104 in FIG. 3D:

AH: corresponds to V.sub.HV, and is fed to the Patches of Turbulence circuit 150 (FIG. 4). AJ: corresponds to Lw in the above Table IV - Turbulence Specifications. AK: corresponds to Lu in the above Table IV - Turbulence Specifications. AL: corresponds either to .sigma..sub.W at the output of the Patches of Turbulence circuit 150 in FIG. 4, or to DISP. ITRBPAT is a flag which chooses the "Patches of Turbulence" .sigma..sub.W or a user specified field .sigma..sub.W in the variable display. AM: corresponds to .sigma..sub.u in the above Table IV - Turbulence Specifications. AN: corresponds to .sigma..sub.v in the above Table IV - Turbulence Specifications. AO: corresponds to Lv in the above Table IV - Turbulence Specifications. AP: corresponds to IVCG, which is calculated as follows: IVCG=(Nv-1).times.0.5, where Nv=�(Nx/R)(FTR)!+2 AQ: corresponds to IVTR, which is calculated as follows: �FTR.times.v.sub.DTI .times.cos(BETA)!+IVCG.

A loop 109 performs the following routine: If V.sub.DTN is greater than, or equal to 2R, no recompute is needed. If, on the other hand, V.sub.DTN is less than 2R, the recompute routine 110 is executed. This corresponds to the velocity being below the minimum velocity allowed by the table size, and performs the corresponding equation above relating to Vuv.

With reference to FIG. 3C the input signals AA, AB, AC, AD, AE, AF, AG to the blade-element turbulence circuit 102 are processed as follows:

AA and AB are input to a vertical gust circuit 200 (shown in detail in FIG. 5) whose output may used when a flag (ITRBGUST=1), or is not used when the flag (ITRBGUST=0), and which corresponds to the newest of TRBGUST array.

AD and AE are added to provide the angle of each blade (one blade at a time) relative to the direction of translation V.sub.H. The cosine of the blade angle is calculated and multiplied by r.sub.m (where r.sub.m =AC.times.AG, as described above) and also by AF, to generate INDX. INDX is then simultaneously fed to three arrays W.sub.A, W.sub.R, W.sub.L, where the array W.sub.A holds the outputs of the Vertical Gusts circuit 200, and the arrays W.sub.R and W.sub.L hold the turbulence values found at the onset point. The sine of the blade angle is also calculated and is processed as follows to generate P.sub.m,n :

p.sub.m,n =�{�(Sine of blade angle)*(r.sub.m)!+�AG!}.times.0.5!+0.5.

Input signals CB and CC are respectively fed from the body-turbulence circuit 104 (FIG. 3D) to the arrays W.sub.R and W.sub.L. p.sub.m,n and the outputs of arrays W.sub.R and W.sub.L are fed to a Gaussian Interpolation circuit W.sub.B whose output OW.sub.B is defined by the following equation: ##EQU13##

The outputs of array W.sub.A and circuit W.sub.B are added to provide the total turbulence input UPGMR(M,N) to the blade segment. It should be noted that the circuit inside the rectangle in dashed lines constitutes a loop which is repeated for each blade-segment combination (i.e., twenty times in the present example).

With reference to FIG. 3D, the input signals AH, AJ, AK, AL, AM, AM, AO, AP, AQ to the body-turbulence circuit 104 are processed as follows:

AH, AJ, and a Gaussian random number Z(.sigma..sub.W) are fed into a Dryden turbulence filter W.sub.C for generating signal CC according to the following equation: ##EQU14##

AH, AJ, and the Gaussian random number z(.sigma..sub.w) are fed into a Dryden turbulence filter W.sub.D for generating signal CB according to the following equation: ##EQU15##

Signals AL, AM and AN are fed into a Gaussian Random Numbers calculation filter W.sub.H, one of whose outputs z(.sigma..sub.u) is fed into a Dryden turbulence filter W.sub.E along with signals AH and AK, for generating a signal OW.sub.E according to the following equation: ##EQU16## Signal OW.sub.E is then multiplied by GAINU to generate output signal UTURB.

Another output z(.sigma..sub.v) of the Gaussian Random Numbers calculation filter W.sub.H is fed simultaneously into two arrays W.sub.K and W.sub.M. Signal AP is the index into array W.sub.K, and signal AQ is the index into array W.sub.M. The output of array W.sub.K is fed into a Dryden turbulence filter W.sub.F along with signals AH and AO to generate a signal OW.sub.E according to the following equation: ##EQU17##

Similarly, the output of array W.sub.M is fed into a Dryden turbulence filter W.sub.G along with signals AH and AO to generate a signal OW.sub.G according to the following equation: ##EQU18##

Each of the signals OW.sub.E and OW.sub.G is then multiplied by GAINV, to generate output signals VTURB and VYGTR, respectively.

The two arrays W.sub.K and W.sub.M are useful for correlating the side gusts entering the center of gravity and the side gusts entering the tail rotor. For instance, if the simulation helicopter were flying forward through a sudden side gust, the center of gravity of the helicopter would respond before the tail rotor enters the gust, thus providing a realistic response. In another situation where the helicopter were flying backward through the sudden side gust, the tail rotor will respond to the gust before the center of gravity does. Therefore, the two arrays W.sub.K and W.sub.M may delay the response of either the center of gravity or the tail rotor relative to each other, based on the velocity of the helicopter to provide a realistic response to gust conditions.

In addition, the present model renders a realistic turbulence simulation and provides a correlation to the terrain and other considerations from the physical environment surrounding the helicopter. To this end, the present model is provided with multiple spatial inputs, such as at the two onset points 20, 21, and at two side input locations: the center of gravity and the tail rotor. Turbulence is initiated for all the segments at the two onset points 20, 21. VTURB signal represents the turbulence velocity that is applied to the center of gravity, and VYGTR represents the lateral velocity in the Y direction that is applied to the tail rotor.

Referring now to FIG. 4, it represents a block diagram of the Patches of Turbulence circuit 150. The waiting time Tw is calculated according to the following equation:

Tw=Time-PTLSIG

A software or firmware routine embedded in the circuit 150 inquires whether Tw is greater than the random wait time PTRNDW. If Tw is greater than the random wait time PTRNDW, the software runs a routine 152 followed by a routine 154. If, on the other hand, Tw does not exceed the random wait time PTRNDW, then the software runs routine 154.

The routine 152 will now be described. Vuv is fed into a look up table for generating the corresponding PLMBD, which is then fed into another look-up table for generating output signal PRMPT. Output signals PTRNDW and PSIGNEW are also generated according to the following equations:

PTRNDW=PLMBD.times.In�0.85(U(0,1))+0.1! PSIGNEW=abs(z(DISP* 1.2533)) where U is a uniform random number, and z is a Gaussian random number.

The routine 154 will now be described. The output signal .sigma..sub.w is determined according to the following equation:

.sigma..sub.w =�(POUT-PSIGNEW) (DT2)/(PRMPT)!+POUT

Referring now to FIG. 5, it represents a block diagram of the Vertical Gusts of Turbulence circuit 200. The circuit 200 is similar to the Patches of Turbulence circuit 150 (FIG. 4), with the exception that it includes two routines 202 and 204. The routine 202 is similar to the routine 152 (FIG. 4) except that GSIG is used according to the following equation, to generate an output signal GUSTNEW:

GUSTNEW=z(GSIG), where z is a Gaussian random number.

While specific embodiments of the system and method for finite element simulation of helicopter turbulence have been illustrated and described in accordance with the present invention, modifications and changes of the apparatus, parameters, use and operation will become apparent to those skilled in the art, without departing from the scope of the invention.

Claims

1. A simulator system for use in real-time piloted helicopter simulation, the helicopter having a center of gravity, rotor disc, a tail rotor, a plurality of blades, each of which including a plurality of blade-element stations, the system comprising in combination:

a processor for applying finite element turbulence simulation to the rotor blades of the helicopter, for creating and propagating three translational components of turbulence in real time, and for superposing turbulence velocity components at each blade element station as a function of helicopter velocity, with distribution according to each blade element instantaneous geometry;

an onset line representing a line of turbulence along the rotor disc, said onset line being tangential to the rotor disc, with said turbulence velocity components being applied at two points along said onset line;

said processor periodically updating said turbulence velocity components and distributing the updated turbulence velocity components to the individual blade-segment stations; and

at least one array W.sub.K, W.sub.M for correlating gusts entering the center of gravity and gusts entering the tail rotor, and for delaying the response of either the center of gravity or the tail rotor relative to each other, based on the velocity of the helicopter to provide a realistic response to gust conditions.

2. The simulator system according to claim 1, wherein said processor includes an input circuit, a blade-element turbulence circuit and a body-turbulence circuit.

3. The simulator system according to claim 2, wherein said input circuit generates a plurality of output signals (AA, AB, AC, AD, AE, AF, AG) that are fed to a body-turbulence circuit, and that generally correspond to, and are derived from the helicopter rotor geometry and dynamics parameters, simulator system variables, and turbulence variables.

4. The simulator system according to claim 3, wherein said input circuit further generates a plurality of output signals (AH, AJ, AK, AL, AM, AN, AO, AP, AQ) that are fed to said blade-element turbulence circuit, and that generally correspond to, and are derived from the helicopter rotor geometry and dynamics parameters, simulator system variables, and turbulence variables.

5. The simulator system according to claim 4, wherein said signals AA and AB are input to a vertical gust circuit whose output may used when a flag ITRBGUST=1, or is not used when said flag ITRBGUST =0.

6. The simulator system according to claim 5, wherein said signals AD and AE are added to provide a blade angle of each blade relative to a direction of translation of the helicopter;

wherein the cosine of said blade angle is calculated and multiplied by r.sub.m (where r.sub.m =AC.times.AG) and also by AF, to generate array index INDX; and

wherein INDX is simultaneously fed to an array (W.sub.A) that holds the outputs of said Vertical Gusts circuit, and to two arrays (W.sub.R, W.sub.L,) that hold the turbulence values at said onset line.

7. The simulator system according to claim 6, wherein the sine of said blade angle is calculated and processed as follows to generate P.sub.m,n :

P.sub.m,n =�{�(Sine of blade angle)*(r.sub.m)!+�AG!}.times.0.5!+0.5.

8. The simulator system according to claim 7, wherein signals CB and CC are respectively fed from said body-turbulence circuit to said arrays W.sub.R and W.sub.L ; and

wherein P.sub.m,n and the outputs of said arrays W.sub.R and W.sub.L are fed to a Gaussian Interpolation circuit whose output OW.sub.B is defined by the following equation: ##EQU19##

9. The simulator system according to claim 8, wherein the outputs of said array W.sub.A and said circuit W.sub.B are added to provide a total turbulence input to said blade-segment.

10. The simulator system according to claim 9, wherein said total turbulence input is determined for each one of said plurality of blade-segments.

11. The simulator system according to claim 10, wherein said signal AH corresponds to a velocity (V.sub.HV) and is fed to a Patches of Turbulence circuit; and

wherein signal AL may correspond either to the output of said a Patches of Turbulence circuit or to an input signal to said Patches of Turbulence circuit.

12. The simulator system according to claim 11, wherein said Patches of Turbulence circuit calculates a waiting time Tw according to the following equation:

Tw=Time-PTLSIG,

where PTLSIC represents the last time;

wherein if Tw is greater than a random wait time PTRNDW, said Patches of Turbulence circuit runs a first routine followed by a second routine, but if Tw does not exceed said random wait time PTRNDW, then said Patches of Turbulence circuit runs said second routine.

13. The simulator system according to claim 12, wherein a rotor-plane aerodynamic velocity Vuv is fed into a look up table for generating a corresponding time delay PLMBD, which is then fed into another look-up table for generating an output signal ramp up time PRMPT.

14. The simulator system according to claim 13, wherein said first routine further generates two output signals random time to wait (PTRNDW) and turbulence magnitude (PSIGNEW) according to the following equations:

PTRNDW=PLMBD.times.In �0.85(U(0,1))+0.1!

PSIGNEW=abs (z(DISP * 1.2533))

where U is a uniform random number, and z is a Gaussian random number.

15. The simulator system according to claim 14, wherein said second routine generates an output signal .sigma..sub.w that represents a turbulence RMS intensity.

16. The simulator system according to claim 15, wherein said Vertical Gusts of Turbulence circuit calculates said ating time Tw, and if Tw is found to be greater than said random wait time PTRNDW, said Vertical Gusts of Turbulence circuit runs a first Vertical Gusts of Turbulence routine followed by a second Vertical Gusts of Turbulence routine, but if Tw does not exceed said random wait time PTRNDW, then said Vertical Gusts of Turbulence circuit runs said second Vertical Gusts of Turbulence routine;

wherein a rotor-plane aerodynamic velocity Vuv is fed into a look up table for generating a corresponding delay time PLMBD, which is then fed into another look-up table for generating an output signal ramp-up time PRMPT;

wherein said first Vertical Gusts of Turbulence routine further generates two output signals time to wait (GTRNDW) and gust magnitude (GUSTNEW) according to the following equations:

GTRNDW=PLMBD.times.In�0.85(U(0,1))+0.1!

GUSTNEW=z(GSIG),

where z is a Gaussian random number, and GSIG is standard derivation of turbulence; and

wherein said second Vertical Gusts of Turbulence routine generates an output signal .sigma..sub.w that represents a turbulence RMS intensity.

17. A method for use in real-time piloted helicopter simulation, the helicopter having a center of gravity, a rotor disc, a rotor tail, and a plurality of blades, each of which including a plurality of blade-element stations, the method comprising the steps of:

applying finite element turbulence simulation to the rotor blades of the helicopter creating and propagating three translational components of turbulence in real time, and superposing turbulence velocity components at each blade element station as a function of helicopter velocity, with distribution according to each blade element instantaneous geometry;

defining an onset line for representing a line of turbulence along the rotor disc, said onset line being tangential to the rotor disc;

applying said turbulence velocity components at two points along said onset line;

periodically updating said turbulence velocity components and distributing the updated turbulence velocity components to the individual blade-segment stations; and

using at least one array W.sub.K, W.sub.M to correlate gusts entering the center of gravity and gusts entering the tail rotor, and for delaying the response of either the center of gravity or the tail rotor relative to each other, based on the velocity of the helicopter to provide a realistic response to gust conditions.

18. The method according to claim 17, wherein an input circuit generates a plurality of output signals (AA, AB, AC, AD, AE, AF, AG) that are fed to a body-turbulence circuit, and that generally correspond to, and are derived from the helicopter rotor geometry and dynamics parameters, simulator system variables, and turbulence variables.

19. The method according to claim 18, wherein said input circuit further generates a plurality of output signals (AH, AJ, AK, AL, AM, AN, AO, AP, AQ) that are fed to said blade-element turbulence circuit, and that generally correspond to, and are derived from the helicopter rotor geometry and dynamics parameters, simulator system variables, and turbulence variables.

20. A software implemented real-time piloted helicopter simulation, wherein the helicopter has a center of gravity, a rotor disc, a rotor tail, and a plurality of blades, each of which including a plurality of blade-element stations, the simulation comprising in combination:

a routine for applying finite element turbulence simulation to the rotor blades of the helicopter, for creating and propagating three translational components of turbulence in real time, and for superposing turbulence velocity components at each blade element station as a function of helicopter velocity, with distribution according to each blade element instantaneous geometry;

a simulated onset line representing a line of turbulence along the rotor disc, said onset line being tangential to the rotor disc, with said turbulence velocity components being applied at two points along said onset line;

said routine periodically updating said turbulence velocity components and distributing the updated turbulence velocity components to the individual blade-segment stations; and

a routine for using at least one array W.sub.K, W.sub.M to correlate gusts entering the center of gravity and gusts entering the tail rotor, and for delaying the response of either the center of gravity or the tail rotor relative to each other, based on the velocity of the helicopter to provide a realistic response to gust conditions.