Helicopter Flight Model Config (XML): Difference between revisions

From Bohemia Interactive Community
Jump to navigation Jump to search
(Created page with "Category:Take On Helicopters | Category:Take On Helicopters: Editing Text copied for RotorLib programmers guide by [http://www.rtdynamics.com RTDynamics'] == What is RotorLib FDM...")
 
m (Text replacement - "[[Image:" to "[[File:")
 
(31 intermediate revisions by 6 users not shown)
Line 1: Line 1:
Category:Take On Helicopters | Category:Take On Helicopters: Editing
{{Messagebox/BI
Text copied for RotorLib programmers guide by [http://www.rtdynamics.com RTDynamics']
|This documentation was reproduced from the official RotorLib Programmer's Guide with permission from {{Link|http://www.rtdynamics.com|RTDynamics}}.
== What is RotorLib FDM ==
|This version of the documentation (4.0) is based on a version of the library that was used in TKOH patch 1.04.
RotorLib FDM is RTDynamics' rotary wing aircraft flight dynamics model solution. It includes aerodynamic models, ground contact model, flight instrument models and flight control system simulation.
[[Arma 3 Helicopters]] DLC introduced new version (4.1.1) that contains simulation of vortex ring state. See [[Arma 3: Advanced Helicopter Flight Model]] for more information.}}


RotorLib FDM is designed for pilot-in-the-loop applications such as flight training simulators, mission rehearsal systems and combat training simulators.


== RotorLib Flight Dynamics Model ==
== What is RotorLib FDM? ==


=== GenHeli500 model ===
RotorLib FDM is RTDynamics' rotary wing aircraft flight dynamics model solution. It includes aerodynamic models, ground contact model, flight instrument models and flight control system simulation.
The generic RotorLib FDM helicopter flight dynamics model is implemented in ''RotorLibFDM::GenHeli500'' class. This class uses several mathematical models as components to create a generic helicopter model. GenHeli500 consists of:


* one main rotor and one tail rotor (both based on ''FGRotorMomentumTheory'' model)
RotorLib FDM is designed for pilot-in-the-loop applications such as flight training simulators, mission rehearsal systems and combat training simulators.


* vertical and horizontal stabilizers ( based on ''FGHelicopterStabilizerWing'' )


* fuselage model ( ''FGHeliFuselage100'' )
== RotorLib Flight Dynamics Model ==


* a user defined number of ground contacts which can be configured as steerable or non-steerable gears and body contacts .
=== GenHeli600 model ===


* extensible engine model
The generic RotorLib FDM helicopter flight dynamics model is implemented in ''RotorLibFDM::GenHeli600 class''. This class uses several mathematical models as components to create a generic helicopter model. GenHeli600 consists of:


* rotor models
* vertical and horizontal stabilizers & wings
* fuselage model
* ground contacts (skids, wheels, inflatable buoys)
* engine model
* drive train model
* drive train model
* flight control system
* flight instruments
Forces and moments computed by these models are used by a real-time rigid body dynamics model to update the rigid body state of the helicopter. Below you can find an overview of the flight dynamics model.


* extensible flight control system
[[File:RLIBDOC01.jpg|thumb|1000px|center]]


* flight instruments
=== GenHeli600 configuration ===


Forces and moments computed by models above are used by a real-time rigid body dynamics model to update the rigid body state of the helicopter. Below you can find an overview of the rigid body dynamics used in RotorLib.
GenHeli600 flight dynamics model is configured through XML configuration files. These configuration files are given to the constructor of the GenHeli600 class.


[[Image:RLIBDOC01.jpg]]
Configuration files start with an FDM type declaration and file version definition. Configuration parameters are defined between the opening and closing tags of this type declaration. The simplified layout of the configuration is provided below.


=== GenHeli500 configuration ===
<syntaxhighlight lang="xml">
GenHeli500 flight dynamics model is configured through XML configuration files. These configuration files are given to the constructor of the ''GenHeli500'' class.
<GenHeli600 version="1.0">
<!-- Physical parameters of the Main Body -->
<MomentsOfInertia Ixx Iyy Izz Iyz /><!-- Kg*m^2 -->
<Mass mass /><!-- Kg -->
<CenterOfGravity x y z /><!-- meters, relative to reference datum -->


Configuration files start with an FDM type declaration and file version definition. Configuration parameters are defined between the opening and closing tags of this type declaration. The simplified layout of the configuration is provided below.  
<Fuselage>
<!-- GenHeli600 fuselage model is implemented using RotorLibFDM::FGHeliFuselage100 model. See Fuselage section below for details. -->
</Fuselage>


<GenHeli500 version="1.0">          
<Engines>
                     
<Engine type>
  <MomentsOfInertia Ixx Iyy Izz Iyz />                          
<!-- Configuration for the chosen engine models must be provided inside the <Engine> </Engine> tags -->
 
</Engine>
  <Mass mass/>                                              
.
 
.
  <Fuselage>        
.
  --> Fuselage parameters <--                         
<Engine type>
  </Fuselage>
<!-- The number of engines that can be added to a helicopter is not limited. -->
</Engine>
</Engines>


  <Engines>
<DriveTrains>
    <Engine type>
<DriveTrain type>
      --> Engine parameters <--
<!-- Configuration for the chosen engine model must be provided inside the <DriveTrain> </DriveTrain> tags -->
    </Engine>
</DriveTrain>
  </Engines>
.
 
.
  <ControlSystem>
.
    <FCSComponent type >
<DriveTrain type>
      --> Flight contol system component parameters <--
<!-- The number of drive trains that can be added to a helicopter is not limited. -->
    </FCSComponent>
</DriveTrain>
    .
</DriveTrains>
    .  
    .
    <FCSComponent type >  
      --> Flight contol system component parameters <--
    </FCSComponent>
  </ControlSystem>
 
  <MainRotor>
    --> Main rotor parameters <--
  </MainRotor>


  <TailRotor>
<ControlSystem>
    --> Tail rotor parameters <--
<CSComponent type >
  </TailRotor>
<!-- Flight control system configuration. In this section a multitude of flight control system components can be defined. -->
 
<!-- Each component should be defined in a separate <CSComponent> section. -->
  <HorizontalStabilizerRight> 
<!-- type parameter in CSComponent defines which component type (e.g. FCSRotationDamper ) will be created. -->
    --> Right horizontal stabilizer parameters <--
</CSComponent>
  </HorizontalStabilizerRight>
.
 
.
  <HorizontalStabilizerLeft>
.
    --> Left horizontal stabilizer parameters <--  
<CSComponent type >
  </HorizontalStabilizerLeft>
<!-- Inside the <CSComponent> </CSComponent> configuration parameters of the chosen component model has to be provided. -->
 
</CSComponent>
  <VerticalStabilizer>
</ControlSystem>
    --> Vertical stabilizer parameters <--
  </VerticalStabilizer> 
   
  <GroundContacts>
    <GeneralConfiguration
                            gearExtensionTime gearDragCoefficient
                            gearReferenceArea maxSteerableGearDeflection
                            maxSteerableGearDeflection
                            steerableGearSpeedCorrectionFactor/>  


    <GroundContact type>
<Rotors>
      --> Ground contact parameters <--
<Rotor type>
    </GroundContact>  
<!-- Type of the rotor model (e.g. FGRotorBladeElementTheory) must be defined with the type attribute in <Rotor> tag. -->
      .
<!-- Configuration of the chosen model resides inside <Rotor></Rotor> tags. -->
      .
</Rotor>
      .
.
    <GroundContact type>
.
      --> Ground contact parameters <--
.
    </GroundContact>  
<Rotor type>
<!-- Number of rotors that can be added is not limited. See Rotor Models section for information about rotor models. -->
</Rotor>
</Rotors>


  </GroundContacts>
<!-- N stabilizers. Name of the stabilizers are to be given by the developers. -->
<StabilizersAndWings>
<StabilizerName_1>
<!-- Each stabilizer is defined within opening and closing tags which are named by the developer. -->
<!-- For example, if it is desired that the first stabilizer is called HorizontalStabilizerRight, -->
<!-- than the stabilizer configuration is made within <HorizontalStabilizerRight> </HorizontalStabilizerRight> tags.-->
<!-- Stabilizers and wings use RTD::FGHelicopterStabilizerWing wing model which supports rotor downwash effects. -->
<!-- See Stabilizers section for details.-->
</StabilizerName_1>
.
.
.
<StabilizerName_n>
<!-- The number of stabilizers and wings is not limited. -->
</StabilizerName_n>
</StabilizersAndWings>
<GroundContacts>
<GeneralConfiguration
gearExtensionTime gearDragCoefficient
gearReferenceArea maxSteerableGearDeflection
maxSteerableGearDeflection
steerableGearSpeedCorrectionFactor/ >
<GroundContact type>
<!-- See Ground Contacts section for details -->
</GroundContact>
.
.
.
<GroundContact type>
<!-- The number of Ground contacts is not limited. -->
</GroundContact>
</GroundContacts>
</GenHeli600>
</syntaxhighlight>


</GenHeli500>
=== Rotor Model ===
RTDynamics SDK provides two rotor models which are implemented in ''FGRotorBladeElementTheory'' and ''FGRotorMomentumTheory'' classes. These models simulate aerodynamic and mechanical behavior of a rotor. Both of the models are suitable as main and tail rotors, however, a good tradeoff between performance and fidelity is obtained when the blade element model is used for the main rotor and the momentum theory based model is used for the tail rotor. Inputs and outputs of the model are summarized below:


* '''MomentsOfInertia''' Moments of inertia of the helicopter body are defined here. Available parameters are Ixx, Iyy, Izz and Iyz (in Kg*m^2)
[[File:RLIBDOC02.jpg|thumb|800px|center|Figure 6.1 Rotor model inputs and outputs]]
* '''Mass''' mass parameter defines the mass of the helicopter in Kg.
* '''Fuselage''' GenHeli500 fuselage model is implemented using RotorLibFDM::FGHeliFuselage100 model. See Fuselage section below for details.
* '''Engines''' In this section, engines are configured. Configuration for the chosen engine model must be provided inside the <Engine> </Engine> tags.
* '''ControlSystem''' Automatic flight control system configuration. In this section a multitude of flight control system components can be defined. Each component should be defined in a separate <FCSComponent> section. type parameter in FCSComponent defines which component type (e.g. FCSRotationDamper ) will be created. Inside the <FCSComponent> </FCSComponent> configuration parameters of the chosen component model has to be provided.
* '''MainRotor and TailRotor''' In these sections main and tail rotors are configured. Rotors use RTD::FGRotorMomentumTheory rotor model. See Rotor model section for details.
* '''HorizontalStabilizerRight , HorizontalStabilizerLeft and VerticalStabilizer''' Stabilizers use RTD::FGHelicopterStabilizerWing wing model. See Stabilizers section for details.
* '''GroundContacts''' In this section ground contacts of the helicopter are configured. See Ground Contacts section for details.


=== Rotor Model ===
Following figure summarizes some of the conventions used in the models. The figure is prepared for a CCW rotor, but except the azimuth (psi) the definitions remain the same. Azimuth angle increases always in the direction of rotation and differs for CCW and CW rotors.
''GenHeli500'' uses ''FGRotorMomentumTheory'' class as rotor model. This class simulates aerodynamic and mechanical behavior of a rotor. It is used for both main and tail rotors. ''FGRotorMomentumTheory'' uses momentum theory and certain blade element theory elements to compute the forces and moments. It also outputs values, such as induced velocity or wake skew angle, that are used as input by other mathematical models, such as fuselage model. Inputs and outputs of the model are summarized below:


'''Rotor model inputs and outputs'''
[[File:RLIBDOC04.jpg|thumb|800px|center|Figure 6.3 Flapping geometry and equation]]


[[Image:RLIBDOC02.jpg]]
Rotor models compute forces (thrust, drag and side force) and moments (pitching, rolling and yawing) at the rotor hub and apply them to the aircraft at its center of gravity. The direction convention for forces and moments is independent of rotor rotation direction.


A helicopter rotor produces thrust by setting the air around it in motion. As a result, a rotor operates in the airflow which is induced by itself. Because of the induced airflow, rotor blades operate in a downward air stream which reduces the effective angle of attack which in turn results in reduced lift. In momentum theory of rotors, the law of conservation of momentum is used on the airflow and an equation that relates T (thrust produced by the rotor) to Vi (speed of the air stream at the rotor disc) is obtained. A second equation that relates T and Vi can be written with blade element theory. Resulting equations are not linear and they are not solvable with algebraic methods hence an iterative numerical solution is implemented.
[[File:RLIBDOC03.jpg|thumb|800px|center|Figure 6.2 Main rotor forces]]


Simulation update of the rotor model begins by solving above equations and finding T (Thrust) and Vi (induced velocity). Once we know T and Vi, H (Drag) and Y(Side force) are computed.]
Rotor models compute induced velocity and rotor wake skew angle which are used by other simulation model components to simulate interference effects, e.g. downwash on fuselage. Induced velocity is positive downwards. Rotor wake skew angle is defined zero downwards and increases bacwards (see figure).


'''Main rotor forces'''
[[File:RLIBDOC05.jpg|thumb|800px|center|Figure 6.4 Rotor wake skew angle]]


[[Image:RLIBDOC03.jpg]]
==== FGRotorBladeElementTheory configuration format ====


The next step in rotor update is updating flapping dynamics. A fully articulated helicopter rotor has flapping, lead-lag and feathering hinges. These hinges provide a certain amount of independent motion capability to the blades. Blades can lead and lag the rotation of the hub on the disc plane. Moreover, because of the flapping hinge, they can move below or above the disc plane. These extra degrees of freedom have effects on power consumption, forces, moments and handling quality. ''FGRotorMomentumTheory'' simulates first order flapping dynamics: It computes a1,b1and a0 (coning angle). Higher order terms in the flapping equation and the lead-lag motion is ignored since their effect on the rotor is less important compared to the first order flapping terms. See following figure for the details of the flapping equation and geometry.
To configure a blade element theory model based rotor, a <Rotor> section with type FGRotorBladeElementTheory must be defined in the configuration file. Configuration parameters of the rotor must be in this section.


'''Flapping geometry and equation'''
The layout of the configuration format is as follows.


[[Image:RLIBDOC04.jpg]]
<syntaxhighlight lang="xml">
<Rotor type="FGRotorBladeElementTheory">


After flapping dynamics, power is computed. Power is computed as a sum of following components:
<Geometry>
<Position x y z /> <!-- RTDynamics coordinate system ( Y forward, X rightward, Z up) relative to the vehicle body coordinate system. -->
<Tilt lateralTilt longitudinalTilt /> <!-- degrees, positive values denote a right/back tilt -->
<Blades radius chord count hingeOffset flappingHingeLowerLimit flappingHingeHigherLimit />
<Rotation clockwise /> <!-- true/false -->
<Delta3 angle /> <!-- pitch-flap coupling in degrees. see "blade azimuth, feathering and flapping diagram" -->
</Geometry>


1. Induced power : Power consumption due to induced air flow.
<MechanicalFeatures Ir Ib maxBrakeTorque />
<!-- Ir : Moment of inertia (rotational inertia) of the rotor about the shaft axis in kg*m^2. -->
<!-- Ib : Moment of inertia (rotational inertia) of a blade about the flapping hinge in kg*m^2. -->
<!-- maxBrakeTorque : Maximum torque that is applied by the rotor brake in N*m. -->


2. Parasite power : Power consumption due to the drag of fuselage, rotor hub, external stores etc.
<AerodynamicFeatures>


3. Climbing power : Power consumption/gain due to altitude change (potential energy change).
<InducedVelocity thrustGain rollMomentGain pitchMomentGain groundEffectFactor />
<!-- unitless, empirical correction factors, default 1 = no effect -->
<!-- thrustFactor: strength of the induced velocity depending on the rotor thrust -->
<!-- groundEffectFactor : strength of the ground effect -->


4. Profile power : Power consumption due to blade profile drag.
<DesignRPM rpm /><!-- If no engine is connected to the rotor, it will rotate with this speed -->


By computing the power consumption of the rotor we also get information about the torque acting on the rotor shaft ( P = angular speed x torque ). This torque is the load the rotor is putting on the engine through the drive train, hence it is an input to the drive train model. Moreover, it also acts on the helicopter fuselage and tries to rotate it in the opposite direction of rotor and hence has effects on flight handling.
<Controls swashPlatePhaseAngle /> <!-- degrees, see "blade azimuth, feathering and flapping diagram" -->


After power, moments at the CoG of the helicopter and rotor wake skew angle are computed. Rotor wake skew angle, along with Vi, is used by other components (fuselage, stabilizers and tail rotor) as input.
<Blades twist><!-- degrees, usually a negativ value -->


'''Rotor wake skew angle'''
<LiftCoefficient>


[[Image:RLIBDOC05.jpg]]
Mach_1 Mach_2 ... Mach_m
AoA_1 Cl_1_1 Cl_1_2 ... Cl_1_m
AoA_2 Cl_2_1 Cl_2_2 ... Cl_2_m


If the model is used to simulate a tail rotor, interference with main rotor and fuselage is taken into account in all above steps.
<!-- . . . -->


==== FGRotorMomentumTheory configuration format ====
AoA_n Cl_n_1 Cl_n_2 ... Cl_n_m
To configure the main rotor, a <MainRotor> section must be defined in the configuration file. Configuration parameters of the rotor has to be in this section. For the tail rotor, <TailRotor> section must be used.


Simplified layout of the model configuration is provided below.
</LiftCoefficient>


<MainRotor>  <-- For tail rotor use <TailRotor>


  <Geometry>
<DragCoefficient>
    <Position x y z />
    <Tilt lateralTilt longitudinalTilt />
    <Blades radius chord count hingeOffset />     
  </Geometry>


  <MechanicalFeatures Ir gearRatio Ib  />
Mach_1 Mach_2 ... Mach_m
AoA_1 Cd_1_1 Cd_1_2 ... Cd_1_m
AoA_2 Cd_2_1 Cd_2_2 ... Cd_2_m
.
.
.
AoA_n Cd_n_1 Cd_n_2 ... Cd_n_m


  <AerodynamicFeatures>
</DragCoefficient>
    <Blades liftCurveSlope zeroLiftDragCoefficient twist />
    <Controls A1UL A1LL B1UL B1LL thetaLL thetaUL />
    <Power inducedPowerCorrectionFactor />     
    <GroundEffect magnitude />
    <DesignRPM rpm/>
  </AerodynamicFeatures>


<MainRotor>
<!-- 2D look up tables describe the lift and drag coefficients for various Mach number and angle of attack combinations. -->
===== Geometry =====
<!-- Lift and drag coefficients are unitless aerodynamic parameters. -->
<!-- Mach numbers must start from zero and go until the largest expected maximum Mach number for the advancing blade's tip. -->
<!-- Angle of attack range must be [-180,180] to ensure that the tables provide a look up value in all flight conditions. -->


This section defines the rotor geometry. All vectors are in RTDynamics coordinate system ( Y forward, X rightward, Z up). They are defined relative to the helicopter body coordinate system.
</Blades>


* Position : position of the rotor hub relative to the center of gravity of the helicopter. In meters.
</AerodynamicFeatures>
* Tilt : orientation of the rotor.
** lateralTilt : Lateral tilt of the rotor in degrees. Positive values denote a rightward tilt.
** longitudinalTilt : Longitudinal tilt of the rotor in degrees. Positive values denote a backward tilt. This parameter is negative for most of the helicopters.
* Blades : Blade parameters.
** count : Number of blades.
** radius : Radius of the blades in meters.
** chord : Chord (width) of the blades in meters.
** hingeOffset : Distance of the flapping hinge from rotation center in meters.


===== MechanicalFeatures : Mechanical features of the rotor. =====
</Rotor>
</syntaxhighlight>


* Ir : Moment of inertia (rotational inertia) of the rotor about the axis of rotation. In kg*m^2.
==== FGRotorMomentumTheory configuration format ====
To configure a momentum theory model based rotor, a <Rotor> section with type ''FGRotorMomentumTheory'' must be defined in the configuration file. Configuration parameters of the rotor must be in this section.


* Ib : Moment of inertia (rotational inertia) of a blade about the flapping hinge. In kg*m^2.
A simplified layout of the model configuration is provided below.


* gearRatio : Gear ratio between the engine and the rotor. rotor RPM = engine RPM / gearRatio.
<syntaxhighlight lang="xml">
<Rotor type="FGRotorMomentumTheory">
<Geometry>
<Position x y z /><!-- RTDynamics coordinate system ( Y forward, X rightward, Z up) relative to the vehicle body coordinate system. -->
<Tilt lateralTilt longitudinalTilt /> <!-- degrees, positive values denote a right/back tilt -->
<Blades radius chord count hingeOffset flappingHingeLowerLimit flappingHingeHigherLimit />
<Rotation clockwise /> <!-- true/false -->
<Delta3 angle /> <!-- pitch-flap coupling in degrees. see "blade azimuth, feathering and flapping diagram" -->
</Geometry>


===== AerodynamicFeatures =====
<MechanicalFeatures Ir Ib maxBrakeTorque />
<!-- Ir : Moment of inertia (rotational inertia) of the rotor about the shaft axis in kg*m^2. -->
<!-- Ib : Moment of inertia (rotational inertia) of a blade about the flapping hinge in kg*m^2. -->
<!-- maxBrakeTorque : Maximum torque that is applied by the rotor brake in N*m. -->


* Blades
<AerodynamicFeatures>
** liftCurveSlope : Slope of the blade lift curve. X axis of the curve is assumed to be in degrees. This parameter does not have a unit. Example: If a lift curve has Clmax = 1.5 at 15 degrees , liftCurveSlope will be 1.5/15 = 0.1.
<Blades liftCurveSlope zeroLiftDragCoefficient twist />
** zeroLiftDragCoefficient : Drag coefficient of the blades at zero lift. This parameter does not have a unit.
<!-- liftCurveSlope : unitless, Slope of the blade lift curve. X axis of the curve is assumed to be in degrees. -->
** twist : Blade twist in degrees.
<!-- Example: If a lift curve has Clmax = 1.5 at 15 degrees , liftCurveSlope will be 1.5/15 = 0.1. -->
* Controls
<!-- zeroLiftDragCoefficient : unitless, Drag coefficient of the blades at zero lift. -->
** A1UL : Upper limit of lateral swash plate tilt in degrees.
<!-- twist : Blade twist in degrees. -->
** A1LL : Lower limit of lateral swash plate tilt in degrees.
<Power inducedPowerCorrectionFactor /> <!-- an empirical correction factor for the computation of the induced power -->
** B1UL : Upper limit of longitudinal swash plate tilt in degrees.
<GroundEffect magnitude /> <!-- unitless, empirical correction factor -->
** B1LL : Lower limit of longitudinal swash plate tilt in degrees.
<DesignRPM rpm /> <!-- If no engine is connected to the rotor, it will rotate with this speed -->
** thetaLL : Minimum pitch of the blades due to collective input. When collective input is 0 blade pitch will be thetaLL degrees.
</AerodynamicFeatures>
** thetaUL : Maximum pitch of the blades due to collective input. When collective input is 1 blade pitch will be thetaUL degrees.
<MainRotor>
* Power
</syntaxhighlight>
inducedPowerCorrectionFactor : Due to various reasons, momentum theory requires an empirical correction factor for the computation of the induced power. This parameter is used for this purpose.
* DesignRPM
rpm : Design RPM of the rotor. If no engine is connected to the helicopter, rotor will rotate with this speed. In revolutions per minute.
* GroundEffect
magnitude : An empirical parameter which defines the magnitude of the ground effect. This parameter does not have a unit.


=== Fuselage ===
=== Fuselage ===
Fuselage model is implemented in ''FGHeliFuselage100'' class. Inputs and outputs of the model are summarized in the following diagram:
Fuselage model is implemented in ''FGHeliFuselage100'' class. Inputs and outputs of the model are summarized in the following diagram:


'''''FGHeliFuselage100'' model inputs and outputs.'''
[[File:RLIBDOC06.jpg|thumb|center|800px|Figure 6.5 FGHeliFuselage100 model inputs and outputs]]
 
[[Image:RLIBDOC06.jpg]]


All forces and moments are computed at the center of gravity.
All forces and moments are computed at the center of gravity.
Line 232: Line 266:
During simulation update, fuselage model initially computes the effective relative wind (in other words, air stream velocity) acting on the fuselage. To compute the effective relative wind helicopter velocity, wind velocity and the main rotor downwash is used. Downwash is characterized by rotor wake skew angle and induced velocity, which are outputs of the rotor model (see Rotor model section for details). Unfortunately, these values can not be used directly to compute the rotor's effect on the fuselage because they are only reference values and the actual induced velocity changes depending on the position on the helicopter. Different parts of the fuselage receive different amount of induced air stream. Therefore there are two curve sets in the configuration file ( ''MainRotorDownwash and MainRotorInplanewash'') which define downwash and inplanewash coefficients depending on the rotor wake skew angle. Using induced velocity output of the main rotor and these curves (as look-up tables indexed by rotor wake angle) final air stream velocity vector is computed.
During simulation update, fuselage model initially computes the effective relative wind (in other words, air stream velocity) acting on the fuselage. To compute the effective relative wind helicopter velocity, wind velocity and the main rotor downwash is used. Downwash is characterized by rotor wake skew angle and induced velocity, which are outputs of the rotor model (see Rotor model section for details). Unfortunately, these values can not be used directly to compute the rotor's effect on the fuselage because they are only reference values and the actual induced velocity changes depending on the position on the helicopter. Different parts of the fuselage receive different amount of induced air stream. Therefore there are two curve sets in the configuration file ( ''MainRotorDownwash and MainRotorInplanewash'') which define downwash and inplanewash coefficients depending on the rotor wake skew angle. Using induced velocity output of the main rotor and these curves (as look-up tables indexed by rotor wake angle) final air stream velocity vector is computed.


'''Velocity of air stream at the fuselage CoP (center of pressure).'''
[[File:RLIBDOC07.jpg|thumb|center|800px|Figure 6.6 Velocity of air stream at the fuselage CoP (center of pressure)]]
 
[[Image:RLIBDOC07.jpg]]


Both the magnitude and direction of the incoming air stream affect the resulting forces and moments so the angle of attack and side slip angles for the effective relative wind are computed. Configuration parameters of the fuselage includes two curves for each force and moment. One of them bases on angle of attack and the other on side slip angle. These are used as look-up tables for computing the forces and the moments.
Both the magnitude and direction of the incoming air stream affect the resulting forces and moments so the angle of attack and side slip angles for the effective relative wind are computed. Configuration parameters of the fuselage includes two curves for each force and moment. One of them bases on angle of attack and the other on side slip angle. These are used as look-up tables for computing the forces and the moments.


'''Forces acting on the fuselage.'''
[[File:RLIBDOC08.jpg|thumb|center|800px|Figure 6.7 Forces acting on the fuselage]]


[[Image:RLIBDOC08.jpg]]


==== ''FGHeliFuselage100'' configuration format ====
==== FGHeliFuselage100 configuration format ====
Simplified layout of the model configuration is provided below.


Simplified layout of the model configuration is provided below.  
<syntaxhighlight lang="xml">
<Fuselage CoPX CoPY CoPZ>
<MainRotorInplanewash>
<Low_a1 a1>
WakeSkewAngle_1 Cinplanewash_1
WakeSkewAngle_2 Cinplanewash_2
.
.
.
WakeSkewAngle_n Cinplanewash_n
</Low_a1>


<Fuselage CoPX CoPY CoPZ>
<Medium_a1 a1>
<!-- As in Low_a1 -->
</Medium_a1>


  <MainRotorInplanewash>
<High_a1 a1>
      <Low_a1 a1>
<!-- As in Low_a1 -->
        WakeSkewAngle_1 Cinplanewash_1
</High_a1>
        WakeSkewAngle_2 Cinplanewash_2
</MainRotorInplanewash>
        .
        .
        .
        WakeSkewAngle_n Cinplanewash_n
      </Low_a1>   
     
      <Medium_a1 a1>
        --> As in Low_a1 <--
      </Medium_a1>   
     
      <High_a1 a1>
        --> As in Low_a1 <--
      </High_a1>                    
    </MainRotorInplanewash>
 
    <MainRotorDownwash>
      <Low_a1 a1>
        WakeSkewAngle_1 Cdownwash_1
        WakeSkewAngle_2 Cdownwash_2
        .
        .
        .
        WakeSkewAngle_n Cdownwash_n
      </Low_a1>   
     
      <Medium_a1 a1>
        --> As in Low_a1 <--
      </Medium_a1>   
     
      <High_a1 a1>
        --> As in Low_a1 <--
      </High_a1> 
    </MainRotorDownwash> 
 
    <LiftDueToAoA>
      AngleOfAttack_1  Cliftaoa_1
      AngleOfAttack_2  Cliftaoa_2
        .
        .
        .
      AngleOfAttack_n  Cliftaoa_n
    </LiftDueToAoA> 


    <LiftDueToSideslip>
<MainRotorDownwash>
      SideSlip_1  Cliftsideslip_1
<Low_a1 a1>
      SideSlip_2  Cliftsideslip_2
WakeSkewAngle_1 Cdownwash_1
        .
WakeSkewAngle_2 Cdownwash_2
        .
.
        .
.
      SideSlip_n  Cliftsideslip_n
.
    </LiftDueToSideslip>
WakeSkewAngle_n Cdownwash_n
</Low_a1>


    <SideForceDueToAoA>
<Medium_a1 a1>
      --> Similar to LiftDueToAoA <--
<!-- As in Low_a1 -->
    </SideForceDueToAoA>
</Medium_a1>


    <SideForceDueToSideslip>
<High_a1 a1>
      --> Similar to LiftDueToSideslip <--
<!-- As in Low_a1 -->
    </SideForceDueToSideslip> 
</High_a1>
   
</MainRotorDownwash>
    <DragDueToAoA>   
      --> Similar to LiftDueToAoA <--
    </DragDueToAoA>
   
    <DragDueToSideslip>
      --> Similar to LiftDueToSideslip <--
    </DragDueToSideslip>


    <PitchingMomentDueToAoA>
<LiftDueToAoA>
      --> Similar to LiftDueToAoA <--
AngleOfAttack_1 Cliftaoa_1
    </PitchingMomentDueToAoA>
AngleOfAttack_2 Cliftaoa_2
.
.
.
AngleOfAttack_n Cliftaoa_n
</LiftDueToAoA>


    <PitchingMomentDueToSideslip>
<LiftDueToSideslip>
      --> Similar to LiftDueToSideslip <--
SideSlip_1 Cliftsideslip_1
    </PitchingMomentDueToSideslip>
SideSlip_2 Cliftsideslip_2
.
.
.
SideSlip_n Cliftsideslip_n
</LiftDueToSideslip>


    <RollingMomentDueToAoA>
<SideForceDueToAoA>
      --> Similar to LiftDueToAoA <--
<!-- Similar to LiftDueToAoA -->
    </RollingMomentDueToAoA>
</SideForceDueToAoA>


    <RollingMomentDueToSideslip>
<SideForceDueToSideslip>
      --> Similar to LiftDueToSideslip <--
<!-- Similar to LiftDueToSideslip -->
    </RollingMomentDueToSideslip>
</SideForceDueToSideslip>


    <YawingMomentDueToAoA>  
<DragDueToAoA>
      --> Similar to LiftDueToAoA <--
<!-- Similar to LiftDueToAoA -->
    </YawingMomentDueToAoA>
</DragDueToAoA>
   
    <YawingMomentDueToSideslip>   
      --> Similar to LiftDueToSideslip <--
    </YawingMomentDueToSideslip>


  </Fuselage>
<DragDueToSideslip>
<!-- Similar to LiftDueToSideslip -->
</DragDueToSideslip>


Fuselage
<PitchingMomentDueToAoA>
<!-- Similar to LiftDueToAoA -->
</PitchingMomentDueToAoA>


CoPX CoPY CoPZ : Position of center of pressure relative to the fuselage center of gravity. These parameters are in meters.
<PitchingMomentDueToSideslip>
<!-- Similar to LiftDueToSideslip -->
</PitchingMomentDueToSideslip>


MainRotorDownwash and MainRotorInplanewash
<RollingMomentDueToAoA>
<!-- Similar to LiftDueToAoA -->
</RollingMomentDueToAoA>


Each of these sections contain three curves which define the relationship between the induced velocity acting on the fuselage and the induced velocity of the main rotor. MainRotorDownwash is used to compute Cdownwash and MainRotorInplanewash is used to compute Cinplanewash. These are used to compute the rotorwash acting on the fuselage as follows:
<RollingMomentDueToSideslip>
<!-- Similar to LiftDueToSideslip -->
</RollingMomentDueToSideslip>
 
<YawingMomentDueToAoA>
<!-- Similar to LiftDueToAoA -->
</YawingMomentDueToAoA>


Downwash (positive downwards) = Main Rotor Induced Velocity * Cdownwash
<YawingMomentDueToSideslip>
<!-- Similar to LiftDueToSideslip -->
</YawingMomentDueToSideslip>
</Fuselage>
</syntaxhighlight>


Inplanewash (positive backwards) = Main Rotor Induced Velocity * Cinplanewash
===== Fuselage =====
CoPX CoPY CoPZ: Position of center of pressure relative to the fuselage center of gravity. These parameters are in meters.
 
====== MainRotorDownwash and MainRotorInplanewash ======
Each of these sections contain three curves which define the relationship between the induced velocity acting on the fuselage and the induced velocity of the main rotor. MainRotorDownwash is used to compute Cdownwash and MainRotorInplanewash is used to compute Cinplanewash. These are used to compute the rotorwash acting on the fuselage as follows:
* Downwash (positive downwards) = Main Rotor Induced Velocity * Cdownwash
* Inplanewash (positive backwards) = Main Rotor Induced Velocity * Cinplanewash


MainRotorDownwash and MainRotorInplanewash can contain three curves for three different a1 (main rotor longitudinal flapping in degrees) values. Depending on current a1, two nearest curves are selected and the coefficients for the current wake skew angle are looked up. Final coefficient is found by linear interpolation between the look-up values.
MainRotorDownwash and MainRotorInplanewash can contain three curves for three different a1 (main rotor longitudinal flapping in degrees) values. Depending on current a1, two nearest curves are selected and the coefficients for the current wake skew angle are looked up. Final coefficient is found by linear interpolation between the look-up values.


Force curves ( Lift, Drag and Side force due to angle of attack and side slip )
====== Force curves (Lift, Drag and Side force due to angle of attack and side slip) ======
 
These curves contain lift, drag and side force coefficients depending on angle of attack and side slip. Coefficients due to angle of attack and side slip are found by using these curves as look-up tables. An example computation is provided below:
These curves contain lift, drag and side force coefficients depending on angle of attack and side slip. Coefficients due to angle of attack and side slip are found by using these curves as look-up tables. An example computation is provided below:
* Clift = Cliftaoa + Cliftsideslip
* Lift = dynamic pressure * Clift


Clift = Cliftaoa + Cliftsideslip
====== Moment curves (Pitching, yawing and rolling moments due to angle of attack and side slip) ======
 
Lift = dynamic pressure * Clift
 
Moment curves (Pitching, yawing and rolling moments due to angle of attack and side slip)
 
These curves contain pitching moment, yawing moment and rolling moment coefficients depending on angle of attack and side slip. Coefficients due to angle of attack and side slip are found by using these curves as look-up tables. An example computation is provided below:
These curves contain pitching moment, yawing moment and rolling moment coefficients depending on angle of attack and side slip. Coefficients due to angle of attack and side slip are found by using these curves as look-up tables. An example computation is provided below:
 
* Cpitch = Cpitchaoa + Cpitchsideslip
Cpitch = Cpitchaoa + Cpitchsideslip
* PitchingMoment = dynamics pressure * Cpitch
 
PitchingMoment = dynamics pressure * Cpitch


=== Stabilizers ===
=== Stabilizers ===
Stabilizers are simulated with F''GHelicopterStabilizerWing'' class. It is derived from ''FGWing'' and inherits most of the functionality from it. It extends ''FGWing'' by adding downwash and dynamic pressure change effects.


'''''FGHelicopterStabilizerWing'' model inputs and outputs.'''
Stabilizers are simulated with ''FGHelicopterStabilizerWing'' class. It is derived from ''FGWing'' and inherits most of the functionality from it. It extends ''FGWing'' by adding downwash and dynamic pressure change effects.


[[Image:RLIBDOC09.jpg]]
[[File:RLIBDOC09.jpg|thumb|center|800px|Figure 6.8 FGHelicopterStabilizerWing model inputs and outputs]]


==== ''FGHelicopterStabilizerWing'' configuration ====
==== ''FGHelicopterStabilizerWing'' configuration ====
GenHeli500 allows definition of one vertical and two horizontal stabilizers. The definitions of the stabilizers are made in <VerticalStabilizer>, <HorizontalStabilizerRight> and <HorizontalStabilizerLeft> sections in the configuration file. Simplified layout of the configuration is provided below.
GenHeli600 allows definition of one vertical and two horizontal stabilizers. The definitions of the stabilizers are made in <VerticalStabilizer>, <HorizontalStabilizerRight> and <HorizontalStabilizerLeft> sections in the configuration file. Simplified layout of the configuration is provided below.
 
<HorizontalStabilizerRight>
 
  <Geometry>
    <Position x y z />
    <Up x y z />
    <Front x y z />
    <SurfaceArea area />     
  </Geometry>
  <AerodynamicFeatures>
    <DynamicPressureRatio>
      <AoA>
        AoA_1 ratio_1
        AoA_2 ratio_2
        .
        .
        .
        AoA_n ratio_n
      </AoA>


      <Sideslip>
<syntaxhighlight lang="xml">
        Sideslip_1 ratio_1
<HorizontalStabilizerRight>
        Sideslip_2 ratio_2
<Geometry>
        .
<Position x y z /><!-- aerodynamic center of the wing with respect to the reference datum in meters -->
        .
<Up x y z /><!-- vector of the wing UP direction in vehicle body coordinates -->
        .
<Front x y z /><!-- vector of the wing FRONT direction in vehicle body coordinates -->
        Sideslip_n ratio_n
<SurfaceArea area /><!-- m2 -->
      </Sideslip>
</Geometry>
    </DynamicPressureRatio>
    <MainRotorDownwash>
      <Low_a1 a1>
        WakeSkew_1 DownwashFactor1
        WakeSkew_2 DownwashFactor2
        .
        .
        .
        WakeSkew_n DownwashFactorn
      </Low_a1>


      <Medium_a1 a1>
<AerodynamicFeatures>
        WakeSkew_1 DownwashFactor_1
<DynamicPressureRatio>
        WakeSkew_2 DownwashFactor_2
<AoA>
        .
AoA_1 ratio_1
        .
AoA_2 ratio_2
        .
.
        WakeSkew_n DownwashFactor_n
.
      </Medium_a1>
.
AoA_n ratio_n
</AoA>


      <High_a1 a1>
<Sideslip>
        WakeSkew_1 DownwashFactor_1
Sideslip_1 ratio_1
        WakeSkew_2 DownwashFactor_2
Sideslip_2 ratio_2
        .
.
        .
.
        .
.
        WakeSkew_n DownwashFactor_n
Sideslip_n ratio_n
      </High_a1>
</Sideslip>
    </MainRotorDownwash>
</DynamicPressureRatio>
     
    <MainRotorInplanewash>
      --> As in MainRotorDownwash <--
    </MainRotorInplanewash>
               
    <LiftCurve>
      AoA_1 Cl_1
      AoA_2 Cl_2
      .
      .
      .
      AoA_n Cl_n
    </LiftCurve>
   
    <DragCurve >
      AoA_1 Cd_1
      AoA_2 Cd_2
      .
      .
      .
      AoA_n Cd_n
    </DragCurve>   


  </AerodynamicFeatures>  
<MainRotorDownwash>
<Low_a1 a1>
WakeSkew_1 DownwashFactor1
WakeSkew_2 DownwashFactor2
.
.
.
WakeSkew_n DownwashFactorn
</Low_a1>


<HorizontalStabilizerRight>
<Medium_a1 a1>
WakeSkew_1 DownwashFactor_1
WakeSkew_2 DownwashFactor_2
.
.
.
WakeSkew_n DownwashFactor_n
</Medium_a1>


Geometry
<High_a1 a1>
WakeSkew_1 DownwashFactor_1
WakeSkew_2 DownwashFactor_2
.
.
.
WakeSkew_n DownwashFactor_n
</High_a1>
</MainRotorDownwash>


Position : Position of the aerodynamic center of the wing with respect to the aircraft CoG in meters. In RTDynamics coordinate system.
<MainRotorInplanewash>
<!-- As in MainRotorDownwash -->
</MainRotorInplanewash>


Up : A vector which determines the up direction of the wing in aircraft body coordinate system.
<LiftCurve>
AoA_1 Cl_1
AoA_2 Cl_2
.
.
.
AoA_n Cl_n
</LiftCurve>


Front : A vector which determines the direction of the front in aircraft body coordinate system.
<DragCurve >
AoA_1 Cd_1
AoA_2 Cd_2
.
.
.
AoA_n Cd_n
</DragCurve>
</AerodynamicFeatures>
<HorizontalStabilizerRight>
</syntaxhighlight>


SurfaceArea
'''DynamicPressureRatio'''
 
area : Wing surface area in square meters.
 
AerodynamicFeatures
 
DynamicPressureRatio


It is possible that due to interference with the fuselage, dynamic pressure observed at the stabilizers are different than the one which is computed using air density and the velocity. This effect is simulated using curves which contain the dynamic pressure ratio vs. fuselage angle of attack and side slip angles. Following equation is used:
It is possible that due to interference with the fuselage, dynamic pressure observed at the stabilizers are different than the one which is computed using air density and the velocity. This effect is simulated using curves which contain the dynamic pressure ratio vs. fuselage angle of attack and side slip angles. Following equation is used:


Stabilizer Dynamic Pressure = Dynamic Pressure * RatioAoA * RatioSideslip
'''Stabilizer Dynamic Pressure = Dynamic Pressure * RatioAoA * RatioSideslip'''


Here RatioAoA is looked up from <AoA> section (angles in degrees, ratio unitless), and RatioSideslip from <sideslip> (angles in degrees, ratio unitless) sections.
Here RatioAoA is looked up from <AoA> section (angles in degrees, ratio unitless), and RatioSideslip from <sideslip> (angles in degrees, ratio unitless) sections.


MainRotorDownwash and MainRotorInplanewash
'''MainRotorDownwash and MainRotorInplanewash'''


Each of these sections contain three curves which define the relationship between the induced velocity acting on the stabilizer and the induced velocity of the main rotor. MainRotorDownwash is used to compute Cdownwash and MainRotorInplanewash is used to compute Cinplanewash. These are used to compute the rotorwash acting on the stabilizer as follows:
Each of these sections contain three curves which define the relationship between the induced velocity acting on the stabilizer and the induced velocity of the main rotor. MainRotorDownwash is used to compute Cdownwash and MainRotorInplanewash is used to compute Cinplanewash. These are used to compute the rotorwash acting on the stabilizer as follows:


Downwash (positive downwards) = Main Rotor Induced Velocity * Cdownwash
'''Downwash (positive downwards) = Main Rotor Induced Velocity * Cdownwash'''


Inplanewash (positive backwards) = Main Rotor Induced Velocity * Cinplanewash
'''Inplanewash (positive backwards) = Main Rotor Induced Velocity * Cinplanewash'''


MainRotorDownwash and MainRotorInplanewash can contain three curves for three different a1 (main rotor longitudinal flapping in degrees) values. Depending on current a1, two nearest curves are selected and the coefficients for the current wake skew angle are looked up. Final coefficient is found by linear interpolation between the look-up values.
MainRotorDownwash and MainRotorInplanewash can contain three curves for three different a1 (main rotor longitudinal flapping in degrees) values. Depending on current a1, two nearest curves are selected and the coefficients for the current wake skew angle are looked up. Final coefficient is found by linear interpolation between the look-up values.


LiftCurve and DragCurve
'''LiftCurve and DragCurve'''


These curves contain the coefficient of lift and drag vs. angle of attack. Coefficients are unitless and angle of attack is in degrees.
These curves contain the coefficient of lift and drag vs. angle of attack. Coefficients are unitless and angle of attack is in degrees.


=== Engine and Drive Train ===
=== Engine and Drive Train ===
RotorLib FDM has an extensible engine model architecture. Engine models are classes that are derived from ''RTD::EngineModelBase''. An existing or user-developed engine model can be attached to the helicopter by defining it in the XML configuration file. To create an engine, an <Engine> section should be defined and engine type has to be set. The configuration parameters of the engine should be between opening and closing tags of the section. Here is an example which creates and configures a ConstantRPMEngine.
RotorLib FDM has an extensible engine model architecture. Engine models are classes that are derived from ''RTD::EngineModelBase''. An existing or user-developed engine model can be attached to the helicopter by defining it in the XML configuration file. To create an engine, an <Engine> section should be defined and engine type has to be set. The configuration parameters of the engine should be between opening and closing tags of the section. Here is an example which creates and configures a ConstantRPMEngine.


<Engine type="ConstantRPMEngine">
<syntaxhighlight lang="xml">
  <Performance maxPower="3000" />
<Engine type="ConstantRPMEngine">
  <ControllerParams P="0.04" I="1.0" D="0.0"  
<Performance maxPower="3000" />
                      offset="200" targetRPM="20900"/>
<ControllerParams P="0.04" I="1.0" D="0.0" offset="200" targetRPM="20900" />
</Engine>
</Engine>
</syntaxhighlight>


Following figure provides an overview of RotorLib's engine and transmission model. The central component of the model is the transmission which receives the torques as input and computes the final RPM's of the rotors and the engine shaft. The resistance torques applied by the rotors to the transmission are the torques generated by the air resistance on the rotor blades. The engine torque shown in the figure is the torque produced by the engine to rotate the rotors. Transmission uses these torques and the information about rotor gear ratios and rotor rotational inertias to compute angular accelerations and finally the final RPM's for the current frame.
Following figure provides an overview of RotorLib's engine and transmission model. The central component of the model is the transmission which receives the torques as input and computes the final RPM's of the rotors and the engine shaft. The resistance torques applied by the rotors to the transmission are the torques generated by the air resistance on the rotor blades. The engine torque shown in the figure is the torque produced by the engine to rotate the rotors. Transmission uses these torques and the information about rotor gear ratios and rotor rotational inertias to compute angular accelerations and finally the final RPM's for the current frame.


'''RotorLib engine and transmission model overview'''
[[File:RLIBDOC10.jpg|thumb|center|800px|Figure 6.9 RotorLib engine and transmission model overview]]
 
[[Image:RLIBDOC10.jpg]]


It is important to note that the engine does not update its own RPM. Engines and rotors provide only the torques acting on the system. Angular accelerations and hence the final RPM's are computed by the transmission model by solving an equation system where solutions for all variables are found simultaneously.
It is important to note that the engine does not update its own RPM. Engines and rotors provide only the torques acting on the system. Angular accelerations and hence the final RPM's are computed by the transmission model by solving an equation system where solutions for all variables are found simultaneously.


==== ConstantRPMEngine ====
==== ConstantRPMEngine ====
ConstantRPMEngine is a helicopter engine model provided with RotorLib library. This engine model has an internal control system that tries to keep the rotor RPM at a given rate. The controller determines the torque output for this purpose. It is possible to define the maximum power output of a ConstantRPMEngine. This power will be available at the sea level. ConstantRPMEngine computes the actual available power depending on density altitude. Configuration format of the model is below.
ConstantRPMEngine engine model has an internal PID controller that tries to maintain the rotor RPM at a given rate. The input of the controller is RPM error ( targetRPM - currentShaftRPM ) and its output is shaft torque. The torque that is requested by the controller is not directly applied to the shaft. Final torque output is computed as follows:
 
'''Pmax = emergencyPower * (airDensity / 1.22406) * HP_to_Watt'''
 
'''Tmax = min( maxTorque , Pmax / (currentShaftRPM * RPM_to_radsec) )'''


<Engine type="ConstantRPMEngine">
'''Toutput = min( Tcontroller - currentShaftRPM2 * rotationResistance , Tmax )'''


  <Performance maxPower />
Configuration format of the model is as follows:
  <ControllerParams P I D offset targetRPM/>


<Engine>
<syntaxhighlight lang="xml">
<Engine type="ConstantRPMEngine">
<General inertia throttleAlwaysFull />
<Performance emergencyPower takeOffPower maximumContinuousPower maxTorque rotationResistance />
<ControllerParams P I D offset targetRPM />
<Engine>
</syntaxhighlight>


Performance
*'''inertia''' : Rotational inertia of the engine in kg*m2.
*'''throttleAlwaysFull''' : When set to true, which is the default, throttle input is ignored and always 100% throttle is assumed.
*'''maximumContinuousPower , takeOffPower , emergencyPower''' : Available engine power for maximum continuous and takeoff ratings in HP. These parameters are not used in the mathematical model but only as reference. TKOH engine has torque limits and damage defined in *.cpp config files.
* '''maxTorque''' : Maximum torque the engine can apply to the shaft. This parameter must be defined in N*m.
* '''rotationResistance''' : Engines have a small resistance to the rotation of the shaft due to internal friction. This resistance reduces the efficiency of the engines. In ConstantRPMengine, resistance is modeled as a torque that is proportional to the square of the RPM. rotationResistance parameter is used to adjust the strength of the resistance and has the unit N.m/ RPM2.
*A '''PID controller''' is used to regulate the engine torque output.
** output = P * error + clamp( I * errorIntegral , integralLimitLow , integralLimitHigh ) + D * errorDerivative + offset
** P, I, D, offset, integralLimitLow and integralLimitHigh are the configuration parameters. error is the input of the PID (in this case the RPM error). errorIntegral and errorDerivative are automatically computed. output is the torque, in N*m, that will be applied to the output shaft.


Performance characteristics of the engine.
=== Control Systems ===
RTDynamics SDK has an extensible control system architecture. Existing and/or user developed control system components can be attached to the vehicles and configured using XML configuration files. Control system components are created and configured in <ControlSystem> section of the configuration file. In this section, for each component ( in other words, for each stage ) a <CSomponent> subsection must be defined and the component type must be set using the ''type'' attribute. Following example defines a single stage of type ''FCSRotationDamper''


maxPower : Maximum power available at sea level. As the density altitude increases available power decreases.
<syntaxhighlight lang="xml">
<ControlSystem>
<CSomponent type="FCSRotationDamper">
<YawController P="2" D="0" />
<PitchController P="2.0" D="0" />
<RollController P="2" D="0.2" />
</CSomponent>
</ControlSystem>
</syntaxhighlight>


ControllerParams
If multiple components are defined in the configuration, they are connected to each other in a cascade, i.e. each stage receives its inputs from the previous stage. Only the first stage receives its inputs directly from the application. A disabled stage passes its inputs to the next stage without any change. By enabling/disabling the stages, different automatic control system configurations can be obtained.


A PID controller is used to regulate the engine torque output. These are parameters of this controller.
[[File:RLIBDOC11.jpg|thumb|center|800px|Figure 6.10 Overview of control system model]]


P : Proportional coefficient of the controller.
Each stage is a control system component. Input and output of each stage is of type ''RTD::InoutData''
 
Control system components are not limited to operate only on their inputs and outputs. ''RTD::InputData'' is only the default input/output format. The information flow to and from a stage can be extended in two ways:
 
* Control system components can read from and write to all simulation variables of a vehicle. E.g. ''GenericControlSystem'' component works this way.


I : Integral coefficient of the controller.
* Components can call each other's class methods to implement master/slave relationships. E.g. RotorLib FDM autopilot implementation works this way.


D : Derivative coefficient of the controller.
==== RTD::GenericControlSystem ====
GenericControlSystem component uses Generic Mathematical Model module.


offset : "Zero point" of the controller. When error, error derivative and error sum is zero resulting value will be offset.
'''MathExpression'''


targetRPM : The engine shaft RPM the controller will try to reach.
A MathExpression operation contains multiple mathematical expressions separated with comma. These mathematical expressions are compiled into a byte code and can be computed very fast in run-time.


=== Automatic flight control system ===
Mathematical expressions consist of standard algebraic operators, parantheses, a conditional operator, a predefined set of functions and the assignment operator ( "=" ). The aliases that are defined in the Variables section can be used to refer to simulation variables in the expressions. More information about the syntax and available functions is provided in the MathExpression Syntax Appendix.
RotorLib FDM has an extensible automatic flight control system architecture. Existing and/or user developed control system components can be attached to the helicopter and configured using XML configuration files. Control system components are created and configured in <ControlSystem> section of the configuration file. In this section, for each component ( in other words, for each stage ) a <FSComponent> subsection has to be defined and component type has to be set using type attribute. Following example defines a single stage of type FCSRotationDamper:


<ControlSystem>
In the following example average engine RPM and pitch rate in deg/sec are computed. The results are written to two newly created simulation variables.
  <FSComponent type="FCSRotationDamper">
  <YawController P="2" D="0"/>
  <PitchController P="2.0" D="0"/>
  <RollController P="2" D="0.2"/>
  </FSComponent>
</ControlSystem>


More info about PID controler concept at Wiki:[http://en.wikipedia.org/wiki/PID_controller]
<syntaxhighlight lang="xml">
<Variables>
<Variable alias="rpm0" name="Engine_0/rpm" />
<Variable alias="rpm1" name="Engine_1/rpm" />
<Variable alias="pitchRateRadSec" member="X" name="angularVelocityBCS" />
<Variable alias="pitchRateDegSec" dataType="REAL" name="pitchRateDegSec" />
<Variable alias="averageRPM" dataType="REAL" name="averageRPM" />
</Variables>


If multiple components are defined in the configuration, they are connected to each other in a cascade, i.e. each stage receives its inputs from the previous stage. Only the first stage receives its inputs directly from the pilot. A disabled stage passes its inputs to the next stage without any change. By enabling/disabling the stages, different automatic control system configurations can be obtained. If all stages are disabled, pilot inputs directly control the actuators.
<MathExpression>
<Expression>
averageRPM = (rpm0+rpm1)/2.0,
pitchRateDegSec = pitchRateRadSec * 180.0 / 3.1415
</Expression>
</MathExpression>
</syntaxhighlight>


'''Overview of control system model.'''
The first three variable declerations refer to existing simulation variables that are created when the vehicle is created and initialized. The next two, pitchRateDegSec and averageRPM , are created by the generic model. We compute new values and assign to these variables in the Expression section.


[[Image:RLIBDOC11.jpg]]
'''RTD::FCSRotationDamper'''


==== RTD::FCSRotationDamper ====
This control system component is used to dampen the rotational motion of the helicopter in all three axes. It can be used for both airplane and helicopters. FCSRotationDamper contains, for each rotational degree of freedom, a PD controller.
This control system component is used to dampen the rotational motion of the helicopter in all three axes. It can be used for both airplane and helicopters. FCSRotationDamper contains, for each rotational degree of freedom, a PD controller.


<FSComponent type="FCSRotationDamper">
<syntaxhighlight lang="xml">
  <YawController P D/>
<FCSComponent type="FCSRotationDamper">
  <PitchController P D/>
<YawController P D />
  <RollController P D/>
<PitchController P D />
<FSComponent type="FCSRotationDamper">
<RollController P D />
</FCSComponent>
</syntaxhighlight>


YawController
YawController
Line 606: Line 647:
D : Derivative coefficient of the controller.
D : Derivative coefficient of the controller.


=== Ground contacts ===
More info about '''PID controller''' concept on Wikipedia: http://en.wikipedia.org/wiki/PID_controller
Interaction with terrain is simulated by ground contacts, classes that are derived from RTD::FGGroundContactBase. Highest fidelity ground contact implementation is RTD::FGGroundContact200 class in the current version. Initially gears are extended.
 
Ground contacts are defined in <GroundContacts> section of the configuration file. Each new ground contact has to be defined in its own section with name <GroundContact>. Also the type of the ground contact has to be declared here using the type attribute. In <GroundContacts> section <GeneralConfiguration> subsection is used to define some parameters that are common to all ground contacts.
 
<GroundContacts>
  <GeneralConfiguration
                          gearExtensionTime gearDragCoefficient
                          gearReferenceArea maxSteerableGearDeflection
                          steerableGearSpeedCorrectionFactor/>
 
  <GroundContact type>
    -> Parameters Related to chosen ground contact model <-
  </GroundContact>
</GroundContacts>
 
 
GeneralConfiguration
 
gearExtensionTime : Time needed to extend/retract landing gears in seconds.
 
gearDragCoefficient : Coefficient of drag for landing gears.
 
gearReferenceArea : Reference area of the landing gears to be used in air drag computation.
 
maxSteerableGearDeflection : Maximum deflection (in degrees) of the steerable gears.
 
steerableGearSpeedCorrectionFactor : As the speed of the aircraft increases, deflections of steerable gears are automatically reduced to keep the aircraft controllable. This parameter defines the strength of this correction.
 
GroundContact
 
A multitude of <GroundContact> </GroundContact> sections can be defined. The type parameter defines which ground contact model (E.g. FGGroundContact200 ) has to be used. Inside the opening and closing tags, configuration parameters of the chosen model has to be provided.
 
 
==== FGGroundContact200 configuration ====
FGGroundContact200 implements a multibody dynamics based ground contact model. It supports inclined and moving surfaces (e.g. for ship deck landing).
 
<GroundContact type="FGGroundContact200">
 
  <Features isBrakable isRetractable isSteerable />
  <Suspension>
    <MountStation x y z/>
    <SliderAxis x y z />
    <SuspensionLength length />
    <SpringParameters springConstant
                                  damperConstant
                                  maxCompression/>
  </Suspension>
  <Wheel mass radius>
    <RollAxis x y z />
    <Friction slidingFriction resistanceTorqueNoBrake
                          resistanceTorqueFullBrake />
  </Wheel>
  <SteeringController maxTorque maxRotSpeed
                                  controllerGain inverseRotation
                                  keepCenter />
  <Damage thresholdLoad />
 
<GroundContact>
 
 
Features
 
General features of the ground contact.
 
isBrakable : Is the ground contact a brakable landing gear? (true or false)
 
isRetractable : Is the ground contact a retractable landing gear? (true or false)
 
isSteerable : Is the ground contact a steerable landing gear? (true or false)
 
Suspension
 
Parameters of the suspension.
 
MountStation : Position at which the landing gear is mounted onto the aircraft body. Relative to center of gravity and in meters.
 
SliderAxis : This is a unit direction vector which determines along which direction suspension slider can move. It is defined relative to aircraft body coordinates in RTDDynamics coordinate system convention.
 
SuspensionLength : Length of the suspension, i.e. the distance from mount station to middle point of wheel when suspensions are not compressed or decompressed. In meters.
 
SpringParameters
 
springConstant : Spring constant of the suspension spring. In N/m.
 
damperConstant: Damper constant of the suspension spring. In N*sec/meter.


maxCompression : Maximum amount the suspension spring can be compressed. This value should be less than or equal to suspension length. In meters.
==== Helicopter Control System Modeling ====
Helicopter flight control is accomplished by controlling the rotor swash plates, which, by axial movements or tilting, can change the blade pitch collectively (same amount of pitch at every blade azimuth) or cyclically (periodically changing amounts of pitch depending on current azimuth of the blade). Aerodynamic effects, that are produced by blade pitch control, can be used to control the thrust and rotational moments that are produced by a rotor.


Wheel
Pilots control the swash plates, and hence the helicopters, with cyclic stick, collective lever and anti-torque pedals. Flight control system modeling is about modeling the way these control devices are connected to the swash plates. Every helicopter has a different mechanical mixing method, different actuators, and different stability augmentation methods. Because of this, GenHeli600 does not connect pilot inputs to the swash plates automatically. Instead, GenericControlSystem component must be used to connect control system outputs to the swash plates.


Parameters of the wheel.
In the following example a simple flight control system is modeled. The first control system stage, which receives the pilot inputs, is a FCSRotationDamper. FCSRotationDamper modifies the inputs to reduce the angular rates. For example, if the helicopter has a positive roll rate (rolling right), FCSRotationDamper will reduce the RTD::PilotInputs::m_rollPitchControl[0] slightly to dampen this roll motion. The output of the stage will include the reduced value.


mass : Mass in Kg.
The second stage of the control system is a GenericControlSystem. Input to this component is the output of the first stage. In this stage, m_rollPitchControl input is used to tilt the main rotor swash plate and hence the rotor itself to produce pitching and rolling moments. m_collective is used to move the main rotor swash plate vertically for controlling the rotor thrust. Finally, m_pedals input is used to move the tail rotor swash plate axially to provide tail rotor thrust control (i.e yaw control).


radius : Radius of the wheel in meters.
'''''Note'''
In reality, most of the helicopters have a more complex mechanical mixer unit. For example, collective lever is usually connected to the tail rotor so that when collective lever is moved, some amount of yaw moment compensation is provided automatically. There are various similar mixing methods that are used to alleviate control coupling problems.''


rotationalInertia : Rotational inertia of the wheel in Kg*m^2.
We have to change the simulation variables that represent swash plate state to connect the control inputs to the swash plates as described above. These variables are listed below. Note that, all variables are in radians :


RollAxis : A direction vector that defines the axis of wheel rotation. It is defined relative to aircraft body coordinates.
* MainRotor/theta0 : Main rotor collective pitch.
* MainRotor/A1 : Main rotor cyclic blade pitch lateral component
* MainRotor/B1: Main rotor cyclic blade pitch longitudinal component
* TailRotor/theta0 : Tail rotor collective pitch
* TailRotor/A1 : Tail rotor cyclic blade pitch lateral component
* TailRotor/B1: Tail rotor cyclic blade pitch longitudinal component


Friction
[[File:RLIBDOC22.jpg|thumb|center|800px|Helicopter control system example]]


slidingFriction : Coefficient of friction between the tire and the ground.
In the example control system, following control ranges are used:


resistanceTorqueNoBrake and resistanceTorqueFullBrake: An effective resistance torque from linear interpolation of resistanceTorqueNoBrake and resistanceTorqueFullBrake depending on the amount of applied wheel brake is applied to the wheel continuously to simulate rolling friction. The resistance torque acts in a direction that opposes the current rotation of the wheel. In N*m.
* Collective [6.5, 20] deg
* Lateral cyclic [-13 , 9] deg
* Longitudinal cyclic [-12 , 16] deg
* Pedals [30,-30] deg


SteeringController
You can find the configuration of the described system below. Note that, ''latCyclic , lonCyclic'' and'' pedals'' are in [-1,1] range. This is why they are mapped into [0,1] range inside the linearInterp function.


maxTorque : The steering controller rotates the wheel left or right depending on the pilot's input for steering on the ground. This parameter defines how much torque can be applied to the wheel for this purpose. This parameter is in N*m.
<syntaxhighlight lang="xml">
<ControlSystem>
<FCSComponent type="FCSRotationDamper">
<YawController P="2" D="0" />
<PitchController P="2.0" D="0" />
<RollController P="0.4" D="0.01" />
</FCSComponent>


maxRotSpeed : This parameter defines how fast the wheel can rotate around the vertical axis. In Deg/sec.
<CSComponent type="GenericControlSystem">
<Variables>
<Variable name="CSComponent_1/ControlOutput/pitchRollControl" member="roll" alias="latCyclic" />
<Variable name="CSComponent_1/ControlOutput/pitchRollControl" member="pitch" alias="lonCyclic" />
<Variable name="CSComponent_1/ControlOutput/pedals" alias="pedals" />
<Variable name="CSComponent_1/ControlOutput/collective" alias="collective" />


controllerGain : The torque applied to the wheel to rotate it is proportional to the difference of pilot's desired wheel direction and current direction multiplied by \e m_steeringControllerGain (the result is clamped at \e m_steeringMaxTorque ). This parameter is unitless.
<Variable name="MainRotor/A1" alias="A1m" />
<Variable name="MainRotor/B1" alias="B1m" />
<Variable name="MainRotor/theta0" alias="Th0m" />


inverseRotation : If this parameter is set to true, the wheel will rotate in the opposite direction of the pilot input (if pilot input is to the right wheel will rotate to the left etc. ). This is necessary if the steering wheel is at the back instead of front.
<Variable name="TailRotor/A1" alias="A1t" />
<Variable name="TailRotor/B1" alias="B1t" />
<Variable name="TailRotor/theta0" alias="Th0t" />
</Variables>


keepCenter : If this parameter is set to true, the wheel will not be affected by the pilot's control input but instead will try to keep the wheel in center always. This is useful to model free rotating wheels (a wheel that can rotate around vertical axis but is not actively controlled by actuators ).
<MathExpression>
<Expression>
A1m = linearInterp( -13 , 9 , (latCyclic + 1)/2 ),
B1m = linearInterp( -12 , 16 , (lonCyclic + 1)/2 ),
Th0m= collective * 13.5 + 6.5,


Damage
A1t = 0.0,
B1t = 0.0,
Th0t = linearInterp( 30 , -30 , (pedals + 1)/2 ),


thresholdLoad : During landing the forces acting on the wheel are transferred to the fuselage. Beyond a certain threshold the transferred forces result in damage in the aircraft. The load defined in thresholdLoad (in N) is the threshold used for this purpose.
A1m = degToRad(A1m),
B1m = degToRad(B1m),
Th0m = degToRad(Th0m),
Th0t = degToRad(Th0t)
</Expression>
</MathExpression>
</CSComponent>
</ControlSystem>
</syntaxhighlight>


==== RotorLibFDM::FCSHelicopterAutoTrim ====
==== RotorLibFDM::FCSHelicopterAutoTrim ====
Line 734: Line 730:
FCSHelicopterAutoTrim is suitable for simulators where simple input devices, such as a joystick, is used. In case real helicopter input devices with trim lock controls are used, FCSHelicopterAutoTrim should not be used.
FCSHelicopterAutoTrim is suitable for simulators where simple input devices, such as a joystick, is used. In case real helicopter input devices with trim lock controls are used, FCSHelicopterAutoTrim should not be used.


<FSComponent type="FCSHelicopterAutoTrim">
<syntaxhighlight lang="xml">
  <YawRateAndSideslipController printDebugInfo>
<FCSComponent type="FCSHelicopterAutoTrim">
    <YawRatePID  P I D offset integralMin integralMax />
<YawRateAndSideslipController printDebugInfo>
    <SideslipPID P I D offset integralMin integralMax />
<YawRatePID  P I D offset integralMin integralMax />
    <Transition speedLow speedHigh />
<SideslipPID P I D offset integralMin integralMax />
    <Maximums maxYawRate maxSideslip />
<Transition speedLow speedHigh />
    <InputMapping nonlinearity />
<Maximums maxYawRate maxSideslip />
  </YawRateAndSideslipController>
<InputMapping nonlinearity />
  <PitchRateController printDebugInfo >  
</YawRateAndSideslipController>
    <PID P I D offset integralMin integralMax />
<PitchRateController printDebugInfo>
    <MaxPitchRate maxPitchRate />
<PID P I D offset integralMin integralMax />
  </PitchRateController>
<MaxPitchRate maxPitchRate />
  <RollRateController printDebugInfo >
</PitchRateController>
    <PID P I D offset integralMin integralMax />
<RollRateController printDebugInfo>
    <MaxRollRate maxRollRate />
<PID P I D offset integralMin integralMax />
  </RollRateController>
<MaxRollRate maxRollRate />
<FSComponent type="FCSHelicopterAutoTrim">
</RollRateController>
</FCSComponent>
</syntaxhighlight>


YawRateAndSideslipController
YawRateAndSideslipController
Line 823: Line 821:


nonlinearity : In some cases, it might be better to pass the pilot input to the control system non-linearly. For example, in a desktop simulator, where the input device is a common joystick, we might want to have a better precision at the zero point. We can define an exponential relationship between pilot's input and the control system input using nonlinearity parameter. The actual input that is given to the control system is computed as follows: input to controller = (|pilot input| ^ nonlinearity) * sign(pilot input).
nonlinearity : In some cases, it might be better to pass the pilot input to the control system non-linearly. For example, in a desktop simulator, where the input device is a common joystick, we might want to have a better precision at the zero point. We can define an exponential relationship between pilot's input and the control system input using nonlinearity parameter. The actual input that is given to the control system is computed as follows: input to controller = (|pilot input| ^ nonlinearity) * sign(pilot input).
=== Ground contacts ===
Interaction with terrain is simulated by ground contacts, classes that are derived from RTD::FGGroundContactBase. Highest fidelity ground contact implementation is RTD::FGGroundContact200 class in the current version. Initially gears are extended.
Ground contacts are defined in <GroundContacts> section of the configuration file. Each new ground contact has to be defined in its own section with name <GroundContact>. Also the type of the ground contact has to be declared here using the type attribute. In <GroundContacts> section <GeneralConfiguration> subsection is used to define some parameters that are common to all ground contacts.
<syntaxhighlight lang="xml">
<GroundContacts>
<GeneralConfiguration
gearExtensionTime gearDragCoefficient
gearReferenceArea maxSteerableGearDeflection
steerableGearSpeedCorrectionFactor />
<GroundContact type>
<!-- Parameters Related to chosen ground contact model -->
</GroundContact>
</GroundContacts>
</syntaxhighlight>
===== GeneralConfiguration =====
* gearExtensionTime: Time needed to extend/retract landing gears in seconds.
* gearDragCoefficient: Coefficient of drag for landing gears.
* gearReferenceArea: Reference area of the landing gears to be used in air drag computation.
* maxSteerableGearDeflection: Maximum deflection (in degrees) of the steerable gears.
* steerableGearSpeedCorrectionFactor: As the speed of the aircraft increases, deflections of steerable gears are automatically reduced to keep the aircraft controllable. This parameter defines the strength of this correction.
===== GroundContact =====
A multitude of <GroundContact> </GroundContact> sections can be defined. The type parameter defines which ground contact model (E.g. FGGroundContact200) has to be used. Inside the opening and closing tags, configuration parameters of the chosen model has to be provided.
==== FGGroundContact200 configuration ====
FGGroundContact200 implements a multibody dynamics based ground contact model. It supports inclined and moving surfaces (e.g. for ship deck landing).
<syntaxhighlight lang="xml">
<GroundContact type="FGGroundContact200">
<Features isBrakable isRetractable isSteerable />
<Suspension>
<MountStation x y z />
<SliderAxis x y z />
<SuspensionLength length />
<SpringParameters
springConstant
damperConstant
maxCompression />
</Suspension>
<Wheel mass radius>
<RollAxis x y z />
<Friction
slidingFriction
resistanceTorqueNoBrake
resistanceTorqueFullBrake />
</Wheel>
<SteeringController
maxTorque
maxRotSpeed
controllerGain
inverseRotation
keepCenter />
<Damage thresholdLoad />
</GroundContact>
</syntaxhighlight>
===== Features =====
General features of the ground contact.
* isBrakable: Is the ground contact a brakable landing gear? (true or false)
* isRetractable: Is the ground contact a retractable landing gear? (true or false)
* isSteerable: Is the ground contact a steerable landing gear? (true or false)
===== Suspension =====
Parameters of the suspension.
* MountStation: Position at which the landing gear is mounted onto the aircraft body. Relative to center of gravity and in meters.
* SliderAxis: This is a unit direction vector which determines along which direction suspension slider can move. It is defined relative to aircraft body coordinates in RTDDynamics coordinate system convention.
* SuspensionLength: Length of the suspension, i.e. the distance from mount station to middle point of wheel when suspensions are not compressed or decompressed. In meters.
* SpringParameters
** springConstant: Spring constant of the suspension spring. In N/m.
** damperConstant: Damper constant of the suspension spring. In N*sec/meter.
** maxCompression: Maximum amount the suspension spring can be compressed. This value should be less than or equal to suspension length. In meters.
===== Wheel =====
Parameters of the wheel.
* mass: Mass in Kg.
* radius: Radius of the wheel in meters.
* rotationalInertia: Rotational inertia of the wheel in Kg*m^2.
** RollAxis: A direction vector that defines the axis of wheel rotation. It is defined relative to aircraft body coordinates.
** Friction
*** slidingFriction: Coefficient of friction between the tire and the ground.
*** resistanceTorqueNoBrake and resistanceTorqueFullBrake: An effective resistance torque from linear interpolation of resistanceTorqueNoBrake and resistanceTorqueFullBrake depending on the amount of applied wheel brake is applied to the wheel continuously to simulate rolling friction. The resistance torque acts in a direction that opposes the current rotation of the wheel. In N*m.
===== SteeringController =====
* maxTorque: The steering controller rotates the wheel left or right depending on the pilot's input for steering on the ground. This parameter defines how much torque can be applied to the wheel for this purpose. This parameter is in N*m.
* maxRotSpeed: This parameter defines how fast the wheel can rotate around the vertical axis. In Deg/sec.
* controllerGain: The torque applied to the wheel to rotate it is proportional to the difference of pilot's desired wheel direction and current direction multiplied by \e m_steeringControllerGain (the result is clamped at \e m_steeringMaxTorque ). This parameter is unitless.
* inverseRotation: If this parameter is set to true, the wheel will rotate in the opposite direction of the pilot input (if pilot input is to the right wheel will rotate to the left etc. ). This is necessary if the steering wheel is at the back instead of front.
* keepCenter: If this parameter is set to true, the wheel will not be affected by the pilot's control input but instead will try to keep the wheel in center always. This is useful to model free rotating wheels (a wheel that can rotate around vertical axis but is not actively controlled by actuators ).
===== Damage =====
* thresholdLoad: During landing the forces acting on the wheel are transferred to the fuselage. Beyond a certain threshold the transferred forces result in damage in the aircraft. The load defined in thresholdLoad (in N) is the threshold used for this purpose.
{{GameCategory|tkoh|Editing}}
{{GameCategory|arma3|Vehicle Configuration}}

Latest revision as of 23:10, 20 November 2023

bi symbol white.png
This documentation was reproduced from the official RotorLib Programmer's Guide with permission from RTDynamics.
This version of the documentation (4.0) is based on a version of the library that was used in TKOH patch 1.04. Arma 3 Helicopters DLC introduced new version (4.1.1) that contains simulation of vortex ring state. See Arma 3: Advanced Helicopter Flight Model for more information.


What is RotorLib FDM?

RotorLib FDM is RTDynamics' rotary wing aircraft flight dynamics model solution. It includes aerodynamic models, ground contact model, flight instrument models and flight control system simulation.

RotorLib FDM is designed for pilot-in-the-loop applications such as flight training simulators, mission rehearsal systems and combat training simulators.


RotorLib Flight Dynamics Model

GenHeli600 model

The generic RotorLib FDM helicopter flight dynamics model is implemented in RotorLibFDM::GenHeli600 class. This class uses several mathematical models as components to create a generic helicopter model. GenHeli600 consists of:

  • rotor models
  • vertical and horizontal stabilizers & wings
  • fuselage model
  • ground contacts (skids, wheels, inflatable buoys)
  • engine model
  • drive train model
  • flight control system
  • flight instruments

Forces and moments computed by these models are used by a real-time rigid body dynamics model to update the rigid body state of the helicopter. Below you can find an overview of the flight dynamics model.

RLIBDOC01.jpg

GenHeli600 configuration

GenHeli600 flight dynamics model is configured through XML configuration files. These configuration files are given to the constructor of the GenHeli600 class.

Configuration files start with an FDM type declaration and file version definition. Configuration parameters are defined between the opening and closing tags of this type declaration. The simplified layout of the configuration is provided below. 

<GenHeli600 version="1.0">
	<!-- Physical parameters of the Main Body -->
	<MomentsOfInertia Ixx Iyy Izz Iyz /><!-- Kg*m^2 -->
	<Mass mass /><!-- Kg -->
	<CenterOfGravity x y z /><!-- meters, relative to reference datum -->

	<Fuselage>
		<!-- GenHeli600 fuselage model is implemented using RotorLibFDM::FGHeliFuselage100 model. See Fuselage section below for details. -->
	</Fuselage>

	<Engines>
		<Engine type>
			<!-- Configuration for the chosen engine models must be provided inside the <Engine> </Engine> tags -->
		</Engine>
		.
		.
		.
		<Engine type>
			<!-- The number of engines that can be added to a helicopter is not limited. -->
		</Engine>
	</Engines>

	<DriveTrains>
		<DriveTrain type>
			<!-- Configuration for the chosen engine model must be provided inside the <DriveTrain> </DriveTrain> tags -->
		</DriveTrain>
		.
		.
		.
		<DriveTrain type>
			<!-- The number of drive trains that can be added to a helicopter is not limited. -->
		</DriveTrain>
	</DriveTrains>

	<ControlSystem>
		<CSComponent type >
			<!-- Flight control system configuration. In this section a multitude of flight control system components can be defined. -->
			<!-- Each component should be defined in a separate <CSComponent> section. -->
			<!-- type parameter in CSComponent defines which component type (e.g. FCSRotationDamper ) will be created. -->
		</CSComponent>
		.
		.
		.
		<CSComponent type >
			<!-- Inside the <CSComponent> </CSComponent> configuration parameters of the chosen component model has to be provided. -->
		</CSComponent>
	</ControlSystem>

	<Rotors>
		<Rotor type>
			<!-- Type of the rotor model (e.g. FGRotorBladeElementTheory) must be defined with the type attribute in <Rotor> tag. -->
			<!-- Configuration of the chosen model resides inside <Rotor></Rotor> tags. -->
		</Rotor>
		.
		.
		.
		<Rotor type>
			<!-- Number of rotors that can be added is not limited. See Rotor Models section for information about rotor models. -->
		</Rotor>
	</Rotors>

	<!-- N stabilizers. Name of the stabilizers are to be given by the developers. -->
	<StabilizersAndWings>
		<StabilizerName_1>
			<!-- Each stabilizer is defined within opening and closing tags which are named by the developer. -->
			<!-- For example, if it is desired that the first stabilizer is called HorizontalStabilizerRight, -->
			<!-- than the stabilizer configuration is made within <HorizontalStabilizerRight> </HorizontalStabilizerRight> tags.-->
			<!-- Stabilizers and wings use RTD::FGHelicopterStabilizerWing wing model which supports rotor downwash effects. -->
			<!-- See Stabilizers section for details.-->
		</StabilizerName_1>
		.
		.
		.
		<StabilizerName_n>
			<!-- The number of stabilizers and wings is not limited. -->
		</StabilizerName_n>
	</StabilizersAndWings>
	<GroundContacts>
		<GeneralConfiguration
			gearExtensionTime gearDragCoefficient
			gearReferenceArea maxSteerableGearDeflection
			maxSteerableGearDeflection
			steerableGearSpeedCorrectionFactor/ >
		<GroundContact type>
			<!-- See Ground Contacts section for details -->
		</GroundContact>
		.
		.
		.
		<GroundContact type>
			<!-- The number of Ground contacts is not limited. -->
		</GroundContact>
	</GroundContacts>
</GenHeli600>

Rotor Model

RTDynamics SDK provides two rotor models which are implemented in FGRotorBladeElementTheory and FGRotorMomentumTheory classes. These models simulate aerodynamic and mechanical behavior of a rotor. Both of the models are suitable as main and tail rotors, however, a good tradeoff between performance and fidelity is obtained when the blade element model is used for the main rotor and the momentum theory based model is used for the tail rotor. Inputs and outputs of the model are summarized below:

Figure 6.1 Rotor model inputs and outputs

Following figure summarizes some of the conventions used in the models. The figure is prepared for a CCW rotor, but except the azimuth (psi) the definitions remain the same. Azimuth angle increases always in the direction of rotation and differs for CCW and CW rotors.

Figure 6.3 Flapping geometry and equation

Rotor models compute forces (thrust, drag and side force) and moments (pitching, rolling and yawing) at the rotor hub and apply them to the aircraft at its center of gravity. The direction convention for forces and moments is independent of rotor rotation direction.

Figure 6.2 Main rotor forces

Rotor models compute induced velocity and rotor wake skew angle which are used by other simulation model components to simulate interference effects, e.g. downwash on fuselage. Induced velocity is positive downwards. Rotor wake skew angle is defined zero downwards and increases bacwards (see figure).

Figure 6.4 Rotor wake skew angle

FGRotorBladeElementTheory configuration format

To configure a blade element theory model based rotor, a <Rotor> section with type FGRotorBladeElementTheory must be defined in the configuration file. Configuration parameters of the rotor must be in this section.

The layout of the configuration format is as follows. 

<Rotor type="FGRotorBladeElementTheory">

	<Geometry>
		<Position x y z /> <!-- RTDynamics coordinate system ( Y forward, X rightward, Z up) relative to the vehicle body coordinate system. -->
		<Tilt lateralTilt longitudinalTilt /> <!-- degrees, positive values denote a right/back tilt -->
		<Blades radius chord count hingeOffset flappingHingeLowerLimit flappingHingeHigherLimit />
		<Rotation clockwise /> <!-- true/false -->
		<Delta3 angle /> <!-- pitch-flap coupling in degrees. see "blade azimuth, feathering and flapping diagram" -->
	</Geometry>

	<MechanicalFeatures Ir Ib maxBrakeTorque />
	<!-- Ir : Moment of inertia (rotational inertia) of the rotor about the shaft axis in kg*m^2. -->
	<!-- Ib : Moment of inertia (rotational inertia) of a blade about the flapping hinge in kg*m^2. -->
	<!-- maxBrakeTorque : Maximum torque that is applied by the rotor brake in N*m. -->

	<AerodynamicFeatures>

		<InducedVelocity thrustGain rollMomentGain pitchMomentGain groundEffectFactor />
		<!-- unitless, empirical correction factors, default 1 = no effect -->
		<!-- thrustFactor: strength of the induced velocity depending on the rotor thrust -->
		<!-- groundEffectFactor : strength of the ground effect -->

		<DesignRPM rpm /><!-- If no engine is connected to the rotor, it will rotate with this speed -->

		<Controls swashPlatePhaseAngle /> <!-- degrees, see "blade azimuth, feathering and flapping diagram" -->

		<Blades twist><!-- degrees, usually a negativ value -->

		<LiftCoefficient>

						Mach_1		Mach_2	...	Mach_m
			AoA_1		Cl_1_1		Cl_1_2	...	Cl_1_m
			AoA_2		Cl_2_1		Cl_2_2	...	Cl_2_m

			<!-- . . . -->

			AoA_n		Cl_n_1		Cl_n_2	...	Cl_n_m

		</LiftCoefficient>


		<DragCoefficient>

						Mach_1		Mach_2	...	Mach_m
			AoA_1		Cd_1_1		Cd_1_2	...	Cd_1_m
			AoA_2		Cd_2_1		Cd_2_2	...	Cd_2_m
			.
			.
			.
			AoA_n		Cd_n_1		Cd_n_2	...	Cd_n_m

		</DragCoefficient>

		<!-- 2D look up tables describe the lift and drag coefficients for various Mach number and angle of attack combinations. -->
		<!-- Lift and drag coefficients are unitless aerodynamic parameters. -->
		<!-- Mach numbers must start from zero and go until the largest expected maximum Mach number for the advancing blade's tip. -->
		<!-- Angle of attack range must be [-180,180] to ensure that the tables provide a look up value in all flight conditions. -->

		</Blades>

	</AerodynamicFeatures>

</Rotor>

FGRotorMomentumTheory configuration format

To configure a momentum theory model based rotor, a <Rotor> section with type FGRotorMomentumTheory must be defined in the configuration file. Configuration parameters of the rotor must be in this section.

A simplified layout of the model configuration is provided below. 

<Rotor type="FGRotorMomentumTheory">
<Geometry>
	<Position x y z /><!-- RTDynamics coordinate system ( Y forward, X rightward, Z up) relative to the vehicle body coordinate system. -->
	<Tilt lateralTilt longitudinalTilt /> <!-- degrees, positive values denote a right/back tilt -->
	<Blades radius chord count hingeOffset flappingHingeLowerLimit flappingHingeHigherLimit />
	<Rotation clockwise /> <!-- true/false -->
	<Delta3 angle /> <!-- pitch-flap coupling in degrees. see "blade azimuth, feathering and flapping diagram" -->
</Geometry>

<MechanicalFeatures Ir Ib maxBrakeTorque />
<!-- Ir : Moment of inertia (rotational inertia) of the rotor about the shaft axis in kg*m^2. -->
<!-- Ib : Moment of inertia (rotational inertia) of a blade about the flapping hinge in kg*m^2. -->
<!-- maxBrakeTorque : Maximum torque that is applied by the rotor brake in N*m. -->

<AerodynamicFeatures>
	<Blades liftCurveSlope zeroLiftDragCoefficient twist />
	<!-- liftCurveSlope : unitless, Slope of the blade lift curve. X axis of the curve is assumed to be in degrees. -->
	<!-- Example: If a lift curve has Clmax = 1.5 at 15 degrees , liftCurveSlope will be 1.5/15 = 0.1. -->
	<!-- zeroLiftDragCoefficient : unitless, Drag coefficient of the blades at zero lift. -->
	<!-- twist : Blade twist in degrees. -->
	<Power inducedPowerCorrectionFactor /> <!-- an empirical correction factor for the computation of the induced power -->
	<GroundEffect magnitude /> <!-- unitless, empirical correction factor -->
	<DesignRPM rpm /> <!-- If no engine is connected to the rotor, it will rotate with this speed -->
</AerodynamicFeatures>
<MainRotor>

Fuselage

Fuselage model is implemented in FGHeliFuselage100 class. Inputs and outputs of the model are summarized in the following diagram:

Figure 6.5 FGHeliFuselage100 model inputs and outputs

All forces and moments are computed at the center of gravity.

During simulation update, fuselage model initially computes the effective relative wind (in other words, air stream velocity) acting on the fuselage. To compute the effective relative wind helicopter velocity, wind velocity and the main rotor downwash is used. Downwash is characterized by rotor wake skew angle and induced velocity, which are outputs of the rotor model (see Rotor model section for details). Unfortunately, these values can not be used directly to compute the rotor's effect on the fuselage because they are only reference values and the actual induced velocity changes depending on the position on the helicopter. Different parts of the fuselage receive different amount of induced air stream. Therefore there are two curve sets in the configuration file ( MainRotorDownwash and MainRotorInplanewash) which define downwash and inplanewash coefficients depending on the rotor wake skew angle. Using induced velocity output of the main rotor and these curves (as look-up tables indexed by rotor wake angle) final air stream velocity vector is computed.

Figure 6.6 Velocity of air stream at the fuselage CoP (center of pressure)

Both the magnitude and direction of the incoming air stream affect the resulting forces and moments so the angle of attack and side slip angles for the effective relative wind are computed. Configuration parameters of the fuselage includes two curves for each force and moment. One of them bases on angle of attack and the other on side slip angle. These are used as look-up tables for computing the forces and the moments.

Figure 6.7 Forces acting on the fuselage


FGHeliFuselage100 configuration format

Simplified layout of the model configuration is provided below. 

<Fuselage CoPX CoPY CoPZ>
	<MainRotorInplanewash>
		<Low_a1 a1>
			WakeSkewAngle_1 Cinplanewash_1
			WakeSkewAngle_2 Cinplanewash_2
			.
			.
			.
			WakeSkewAngle_n Cinplanewash_n
		</Low_a1>

		<Medium_a1 a1>
			<!-- As in Low_a1 -->
		</Medium_a1>

		<High_a1 a1>
			<!-- As in Low_a1 -->
		</High_a1>
	</MainRotorInplanewash>

	<MainRotorDownwash>
		<Low_a1 a1>
			WakeSkewAngle_1 Cdownwash_1
			WakeSkewAngle_2 Cdownwash_2
			.
			.
			.
			WakeSkewAngle_n Cdownwash_n
		</Low_a1>

		<Medium_a1 a1>
			<!-- As in Low_a1 -->
		</Medium_a1>

		<High_a1 a1>
			<!-- As in Low_a1 -->
		</High_a1>
	</MainRotorDownwash>

	<LiftDueToAoA>
		AngleOfAttack_1		Cliftaoa_1
		AngleOfAttack_2		Cliftaoa_2
		.
		.
		.
		AngleOfAttack_n		Cliftaoa_n
	</LiftDueToAoA>

	<LiftDueToSideslip>
		SideSlip_1		Cliftsideslip_1
		SideSlip_2		Cliftsideslip_2
		.
		.
		.
		SideSlip_n		Cliftsideslip_n
	</LiftDueToSideslip>

	<SideForceDueToAoA>
		<!-- Similar to LiftDueToAoA -->
	</SideForceDueToAoA>

	<SideForceDueToSideslip>
		<!-- Similar to LiftDueToSideslip -->
	</SideForceDueToSideslip>

	<DragDueToAoA>
		<!-- Similar to LiftDueToAoA -->
	</DragDueToAoA>

	<DragDueToSideslip>
		<!-- Similar to LiftDueToSideslip -->
	</DragDueToSideslip>

	<PitchingMomentDueToAoA>
		<!-- Similar to LiftDueToAoA -->
	</PitchingMomentDueToAoA>

	<PitchingMomentDueToSideslip>
		<!-- Similar to LiftDueToSideslip -->
	</PitchingMomentDueToSideslip>

	<RollingMomentDueToAoA>
		<!-- Similar to LiftDueToAoA -->
	</RollingMomentDueToAoA>

	<RollingMomentDueToSideslip>
		<!-- Similar to LiftDueToSideslip -->
	</RollingMomentDueToSideslip>

	<YawingMomentDueToAoA>
		<!-- Similar to LiftDueToAoA -->
	</YawingMomentDueToAoA>

	<YawingMomentDueToSideslip>
		<!-- Similar to LiftDueToSideslip -->
	</YawingMomentDueToSideslip>
</Fuselage>
Fuselage

CoPX CoPY CoPZ: Position of center of pressure relative to the fuselage center of gravity. These parameters are in meters.

MainRotorDownwash and MainRotorInplanewash

Each of these sections contain three curves which define the relationship between the induced velocity acting on the fuselage and the induced velocity of the main rotor. MainRotorDownwash is used to compute Cdownwash and MainRotorInplanewash is used to compute Cinplanewash. These are used to compute the rotorwash acting on the fuselage as follows:

  • Downwash (positive downwards) = Main Rotor Induced Velocity * Cdownwash
  • Inplanewash (positive backwards) = Main Rotor Induced Velocity * Cinplanewash

MainRotorDownwash and MainRotorInplanewash can contain three curves for three different a1 (main rotor longitudinal flapping in degrees) values. Depending on current a1, two nearest curves are selected and the coefficients for the current wake skew angle are looked up. Final coefficient is found by linear interpolation between the look-up values.

Force curves (Lift, Drag and Side force due to angle of attack and side slip)

These curves contain lift, drag and side force coefficients depending on angle of attack and side slip. Coefficients due to angle of attack and side slip are found by using these curves as look-up tables. An example computation is provided below:

  • Clift = Cliftaoa + Cliftsideslip
  • Lift = dynamic pressure * Clift
Moment curves (Pitching, yawing and rolling moments due to angle of attack and side slip)

These curves contain pitching moment, yawing moment and rolling moment coefficients depending on angle of attack and side slip. Coefficients due to angle of attack and side slip are found by using these curves as look-up tables. An example computation is provided below:

  • Cpitch = Cpitchaoa + Cpitchsideslip
  • PitchingMoment = dynamics pressure * Cpitch

Stabilizers

Stabilizers are simulated with FGHelicopterStabilizerWing class. It is derived from FGWing and inherits most of the functionality from it. It extends FGWing by adding downwash and dynamic pressure change effects.

Figure 6.8 FGHelicopterStabilizerWing model inputs and outputs

FGHelicopterStabilizerWing configuration

GenHeli600 allows definition of one vertical and two horizontal stabilizers. The definitions of the stabilizers are made in <VerticalStabilizer>, <HorizontalStabilizerRight> and <HorizontalStabilizerLeft> sections in the configuration file. Simplified layout of the configuration is provided below. 

<HorizontalStabilizerRight>
	<Geometry>
		<Position x y z /><!-- aerodynamic center of the wing with respect to the reference datum in meters -->
		<Up x y z /><!-- vector of the wing UP direction in vehicle body coordinates -->
		<Front x y z /><!-- vector of the wing FRONT direction in vehicle body coordinates -->
		<SurfaceArea area /><!-- m2 -->
	</Geometry>

	<AerodynamicFeatures>
		<DynamicPressureRatio>
			<AoA>
				AoA_1 ratio_1
				AoA_2 ratio_2
				.
				.
				.
				AoA_n ratio_n
			</AoA>

			<Sideslip>
				Sideslip_1 ratio_1
				Sideslip_2 ratio_2
				.
				.
				.
				Sideslip_n ratio_n
			</Sideslip>
		</DynamicPressureRatio>

		<MainRotorDownwash>
			<Low_a1 a1>
				WakeSkew_1 DownwashFactor1
				WakeSkew_2 DownwashFactor2
				.
				.
				.
				WakeSkew_n DownwashFactorn
			</Low_a1>

			<Medium_a1 a1>
				WakeSkew_1 DownwashFactor_1
				WakeSkew_2 DownwashFactor_2
				.
				.
				.
				WakeSkew_n DownwashFactor_n
			</Medium_a1>

			<High_a1 a1>
				WakeSkew_1 DownwashFactor_1
				WakeSkew_2 DownwashFactor_2
				.
				.
				.
				WakeSkew_n DownwashFactor_n
			</High_a1>
		</MainRotorDownwash>

		<MainRotorInplanewash>
			<!-- As in MainRotorDownwash -->
		</MainRotorInplanewash>

		<LiftCurve>
			AoA_1 Cl_1
			AoA_2 Cl_2
			.
			.
			.
			AoA_n Cl_n
		</LiftCurve>

		<DragCurve >
			AoA_1 Cd_1
			AoA_2 Cd_2
			.
			.
			.
			AoA_n Cd_n
		</DragCurve>
	</AerodynamicFeatures>
<HorizontalStabilizerRight>

DynamicPressureRatio

It is possible that due to interference with the fuselage, dynamic pressure observed at the stabilizers are different than the one which is computed using air density and the velocity. This effect is simulated using curves which contain the dynamic pressure ratio vs. fuselage angle of attack and side slip angles. Following equation is used:

Stabilizer Dynamic Pressure = Dynamic Pressure * RatioAoA * RatioSideslip

Here RatioAoA is looked up from <AoA> section (angles in degrees, ratio unitless), and RatioSideslip from <sideslip> (angles in degrees, ratio unitless) sections.

MainRotorDownwash and MainRotorInplanewash

Each of these sections contain three curves which define the relationship between the induced velocity acting on the stabilizer and the induced velocity of the main rotor. MainRotorDownwash is used to compute Cdownwash and MainRotorInplanewash is used to compute Cinplanewash. These are used to compute the rotorwash acting on the stabilizer as follows:

Downwash (positive downwards) = Main Rotor Induced Velocity * Cdownwash

Inplanewash (positive backwards) = Main Rotor Induced Velocity * Cinplanewash

MainRotorDownwash and MainRotorInplanewash can contain three curves for three different a1 (main rotor longitudinal flapping in degrees) values. Depending on current a1, two nearest curves are selected and the coefficients for the current wake skew angle are looked up. Final coefficient is found by linear interpolation between the look-up values.

LiftCurve and DragCurve

These curves contain the coefficient of lift and drag vs. angle of attack. Coefficients are unitless and angle of attack is in degrees.

Engine and Drive Train

RotorLib FDM has an extensible engine model architecture. Engine models are classes that are derived from RTD::EngineModelBase. An existing or user-developed engine model can be attached to the helicopter by defining it in the XML configuration file. To create an engine, an <Engine> section should be defined and engine type has to be set. The configuration parameters of the engine should be between opening and closing tags of the section. Here is an example which creates and configures a ConstantRPMEngine. 

<Engine type="ConstantRPMEngine">
	<Performance maxPower="3000" />
	<ControllerParams P="0.04" I="1.0" D="0.0" offset="200" targetRPM="20900" />
</Engine>

Following figure provides an overview of RotorLib's engine and transmission model. The central component of the model is the transmission which receives the torques as input and computes the final RPM's of the rotors and the engine shaft. The resistance torques applied by the rotors to the transmission are the torques generated by the air resistance on the rotor blades. The engine torque shown in the figure is the torque produced by the engine to rotate the rotors. Transmission uses these torques and the information about rotor gear ratios and rotor rotational inertias to compute angular accelerations and finally the final RPM's for the current frame.

Figure 6.9 RotorLib engine and transmission model overview

It is important to note that the engine does not update its own RPM. Engines and rotors provide only the torques acting on the system. Angular accelerations and hence the final RPM's are computed by the transmission model by solving an equation system where solutions for all variables are found simultaneously.

ConstantRPMEngine

ConstantRPMEngine engine model has an internal PID controller that tries to maintain the rotor RPM at a given rate. The input of the controller is RPM error ( targetRPM - currentShaftRPM ) and its output is shaft torque. The torque that is requested by the controller is not directly applied to the shaft. Final torque output is computed as follows:

Pmax = emergencyPower * (airDensity / 1.22406) * HP_to_Watt

Tmax = min( maxTorque , Pmax / (currentShaftRPM * RPM_to_radsec) )

Toutput = min( Tcontroller - currentShaftRPM2 * rotationResistance , Tmax )

Configuration format of the model is as follows: 

<Engine type="ConstantRPMEngine">
	<General inertia throttleAlwaysFull />
	<Performance emergencyPower takeOffPower maximumContinuousPower maxTorque rotationResistance />
	<ControllerParams P I D offset targetRPM />
<Engine>
  • inertia : Rotational inertia of the engine in kg*m2.
  • throttleAlwaysFull : When set to true, which is the default, throttle input is ignored and always 100% throttle is assumed.
  • maximumContinuousPower , takeOffPower , emergencyPower : Available engine power for maximum continuous and takeoff ratings in HP. These parameters are not used in the mathematical model but only as reference. TKOH engine has torque limits and damage defined in *.cpp config files.
  • maxTorque : Maximum torque the engine can apply to the shaft. This parameter must be defined in N*m.
  • rotationResistance : Engines have a small resistance to the rotation of the shaft due to internal friction. This resistance reduces the efficiency of the engines. In ConstantRPMengine, resistance is modeled as a torque that is proportional to the square of the RPM. rotationResistance parameter is used to adjust the strength of the resistance and has the unit N.m/ RPM2.
  • A PID controller is used to regulate the engine torque output.
    • output = P * error + clamp( I * errorIntegral , integralLimitLow , integralLimitHigh ) + D * errorDerivative + offset
    • P, I, D, offset, integralLimitLow and integralLimitHigh are the configuration parameters. error is the input of the PID (in this case the RPM error). errorIntegral and errorDerivative are automatically computed. output is the torque, in N*m, that will be applied to the output shaft.

Control Systems

RTDynamics SDK has an extensible control system architecture. Existing and/or user developed control system components can be attached to the vehicles and configured using XML configuration files. Control system components are created and configured in <ControlSystem> section of the configuration file. In this section, for each component ( in other words, for each stage ) a <CSomponent> subsection must be defined and the component type must be set using the type attribute. Following example defines a single stage of type FCSRotationDamper 

<ControlSystem>
	<CSomponent type="FCSRotationDamper">
		<YawController P="2" D="0" />
		<PitchController P="2.0" D="0" />
		<RollController P="2" D="0.2" />
	</CSomponent>
</ControlSystem>

If multiple components are defined in the configuration, they are connected to each other in a cascade, i.e. each stage receives its inputs from the previous stage. Only the first stage receives its inputs directly from the application. A disabled stage passes its inputs to the next stage without any change. By enabling/disabling the stages, different automatic control system configurations can be obtained.

Figure 6.10 Overview of control system model

Each stage is a control system component. Input and output of each stage is of type RTD::InoutData

Control system components are not limited to operate only on their inputs and outputs. RTD::InputData is only the default input/output format. The information flow to and from a stage can be extended in two ways:

  • Control system components can read from and write to all simulation variables of a vehicle. E.g. GenericControlSystem component works this way.
  • Components can call each other's class methods to implement master/slave relationships. E.g. RotorLib FDM autopilot implementation works this way.

RTD::GenericControlSystem

GenericControlSystem component uses Generic Mathematical Model module.

MathExpression

A MathExpression operation contains multiple mathematical expressions separated with comma. These mathematical expressions are compiled into a byte code and can be computed very fast in run-time.

Mathematical expressions consist of standard algebraic operators, parantheses, a conditional operator, a predefined set of functions and the assignment operator ( "=" ). The aliases that are defined in the Variables section can be used to refer to simulation variables in the expressions. More information about the syntax and available functions is provided in the MathExpression Syntax Appendix.

In the following example average engine RPM and pitch rate in deg/sec are computed. The results are written to two newly created simulation variables. 

<Variables>
	<Variable alias="rpm0" name="Engine_0/rpm" />
	<Variable alias="rpm1" name="Engine_1/rpm" />
	<Variable alias="pitchRateRadSec" member="X" name="angularVelocityBCS" />
	<Variable alias="pitchRateDegSec" dataType="REAL" name="pitchRateDegSec" />
	<Variable alias="averageRPM" dataType="REAL" name="averageRPM" />
</Variables>

<MathExpression>
	<Expression>
		averageRPM = (rpm0+rpm1)/2.0,
		pitchRateDegSec = pitchRateRadSec * 180.0 / 3.1415
	</Expression>
</MathExpression>

The first three variable declerations refer to existing simulation variables that are created when the vehicle is created and initialized. The next two, pitchRateDegSec and averageRPM , are created by the generic model. We compute new values and assign to these variables in the Expression section.

RTD::FCSRotationDamper

This control system component is used to dampen the rotational motion of the helicopter in all three axes. It can be used for both airplane and helicopters. FCSRotationDamper contains, for each rotational degree of freedom, a PD controller. 

<FCSComponent type="FCSRotationDamper">
	<YawController P D />
	<PitchController P D />
	<RollController P D />
</FCSComponent>

YawController

Configures the PD controller that damps yawing motion.

P : Proportional coefficient of the controller.

D : Derivative coefficient of the controller.

PitchController

Configures the PD controller that damps pitching motion.

P : Proportional coefficient of the controller.

D : Derivative coefficient of the controller.

RollController

Configures the PD controller that damps rolling motion.

P : Proportional coefficient of the controller.

D : Derivative coefficient of the controller.

More info about PID controller concept on Wikipedia: http://en.wikipedia.org/wiki/PID_controller

Helicopter Control System Modeling

Helicopter flight control is accomplished by controlling the rotor swash plates, which, by axial movements or tilting, can change the blade pitch collectively (same amount of pitch at every blade azimuth) or cyclically (periodically changing amounts of pitch depending on current azimuth of the blade). Aerodynamic effects, that are produced by blade pitch control, can be used to control the thrust and rotational moments that are produced by a rotor.

Pilots control the swash plates, and hence the helicopters, with cyclic stick, collective lever and anti-torque pedals. Flight control system modeling is about modeling the way these control devices are connected to the swash plates. Every helicopter has a different mechanical mixing method, different actuators, and different stability augmentation methods. Because of this, GenHeli600 does not connect pilot inputs to the swash plates automatically. Instead, GenericControlSystem component must be used to connect control system outputs to the swash plates.

In the following example a simple flight control system is modeled. The first control system stage, which receives the pilot inputs, is a FCSRotationDamper. FCSRotationDamper modifies the inputs to reduce the angular rates. For example, if the helicopter has a positive roll rate (rolling right), FCSRotationDamper will reduce the RTD::PilotInputs::m_rollPitchControl[0] slightly to dampen this roll motion. The output of the stage will include the reduced value.

The second stage of the control system is a GenericControlSystem. Input to this component is the output of the first stage. In this stage, m_rollPitchControl input is used to tilt the main rotor swash plate and hence the rotor itself to produce pitching and rolling moments. m_collective is used to move the main rotor swash plate vertically for controlling the rotor thrust. Finally, m_pedals input is used to move the tail rotor swash plate axially to provide tail rotor thrust control (i.e yaw control).

Note In reality, most of the helicopters have a more complex mechanical mixer unit. For example, collective lever is usually connected to the tail rotor so that when collective lever is moved, some amount of yaw moment compensation is provided automatically. There are various similar mixing methods that are used to alleviate control coupling problems.

We have to change the simulation variables that represent swash plate state to connect the control inputs to the swash plates as described above. These variables are listed below. Note that, all variables are in radians :

  • MainRotor/theta0 : Main rotor collective pitch.
  • MainRotor/A1 : Main rotor cyclic blade pitch lateral component
  • MainRotor/B1: Main rotor cyclic blade pitch longitudinal component
  • TailRotor/theta0 : Tail rotor collective pitch
  • TailRotor/A1 : Tail rotor cyclic blade pitch lateral component
  • TailRotor/B1: Tail rotor cyclic blade pitch longitudinal component
Helicopter control system example

In the example control system, following control ranges are used:

  • Collective [6.5, 20] deg
  • Lateral cyclic [-13 , 9] deg
  • Longitudinal cyclic [-12 , 16] deg
  • Pedals [30,-30] deg

You can find the configuration of the described system below. Note that, latCyclic , lonCyclic and pedals are in [-1,1] range. This is why they are mapped into [0,1] range inside the linearInterp function. 

<ControlSystem>
	<FCSComponent type="FCSRotationDamper">
		<YawController P="2" D="0" />
		<PitchController P="2.0" D="0" />
		<RollController P="0.4" D="0.01" />
	</FCSComponent>

	<CSComponent type="GenericControlSystem">
		<Variables>
			<Variable name="CSComponent_1/ControlOutput/pitchRollControl" member="roll" alias="latCyclic" />
			<Variable name="CSComponent_1/ControlOutput/pitchRollControl" member="pitch" alias="lonCyclic" />
			<Variable name="CSComponent_1/ControlOutput/pedals" alias="pedals" />
			<Variable name="CSComponent_1/ControlOutput/collective" alias="collective" />

			<Variable name="MainRotor/A1" alias="A1m" />
			<Variable name="MainRotor/B1" alias="B1m" />
			<Variable name="MainRotor/theta0" alias="Th0m" />

			<Variable name="TailRotor/A1" alias="A1t" />
			<Variable name="TailRotor/B1" alias="B1t" />
			<Variable name="TailRotor/theta0" alias="Th0t" />
		</Variables>

		<MathExpression>
			<Expression>
				A1m = linearInterp( -13 , 9 , (latCyclic + 1)/2 ),
				B1m = linearInterp( -12 , 16 , (lonCyclic + 1)/2 ),
				Th0m= collective * 13.5 + 6.5,

				A1t = 0.0,
				B1t = 0.0,
				Th0t = linearInterp( 30 , -30 , (pedals + 1)/2 ),

				A1m = degToRad(A1m),
				B1m = degToRad(B1m),
				Th0m = degToRad(Th0m),
				Th0t = degToRad(Th0t)
			</Expression>
		</MathExpression>
	</CSComponent>
</ControlSystem>

RotorLibFDM::FCSHelicopterAutoTrim

This control system component is used to trim the helicopter automatically and hence remove the coupling between controls (i.e. collective / pedal coupling). It translates pilot inputs to requested heading, pitch and roll rates or to requested side slip and controls the actuators to achieve the desired state. The final result is a simpler flight.

FCSHelicopterAutoTrim is suitable for simulators where simple input devices, such as a joystick, is used. In case real helicopter input devices with trim lock controls are used, FCSHelicopterAutoTrim should not be used. 

<FCSComponent type="FCSHelicopterAutoTrim">
	<YawRateAndSideslipController printDebugInfo>
		<YawRatePID  P I D offset integralMin integralMax />
		<SideslipPID P I D offset integralMin integralMax />
		<Transition speedLow speedHigh />
		<Maximums maxYawRate maxSideslip />
		<InputMapping nonlinearity />
	</YawRateAndSideslipController>
	<PitchRateController printDebugInfo>
		<PID P I D offset integralMin integralMax />
		<MaxPitchRate maxPitchRate />
	</PitchRateController>
	<RollRateController printDebugInfo>
		<PID P I D offset integralMin integralMax />
		<MaxRollRate maxRollRate />
	</RollRateController>
</FCSComponent>

YawRateAndSideslipController

This controller controls the pedals to reach and maintain a certain yaw rate or side slip. Yaw rate control mode is active at low speeds, whereas side slip control mode is activated at high speeds. The mode change occurs smoothly within a user defined speed range (see Transition below). Within this transition region, both yaw rate and side slip control modes are active but their contribution to the final output varies depending on speed.

printDebugInfo : A boolean parameter (possible values "true" and "false" ) to enable/disable debug output to the console from the controller. This is useful for the parametrization of the controller.

YawratePID: Parameters of the PID controller that is used to control yaw rate. Error input to this PID is desired yaw rate - current yaw rate. Desired yaw rate is determined by the pedal input, i.e. pilot's pedal input is translated to a desired yaw rate by the controller.

P : Proportional coefficient of the controller.

I : Integral coefficient of the controller.

D : Derivative coefficient of the controller.

offset : "Zero point" of the controller. When error, error derivative and error sum is zero resulting value will be offset.

IntegralMin and integralMax : Minimum/Maximum value the error sum can reach. If the sum of errors go below/above this value it is clamped to this value.

SideslipPID: Parameters of the PID controller that is used to control side slip. Error input to this PID is desired side slip - current side slip. Desired side slip is determined by the pedal input, i.e. pilot's pedal input is translated to a desired side slip by the controller.

P : Proportional coefficient of the controller.

I : Integral coefficient of the controller.

D : Derivative coefficient of the controller.

offset : "Zero point" of the controller. When error, error derivative and error sum is zero resulting value will be offset.

IntegralMin and integralMax : Minimum/Maximum value the error sum can reach. If the sum of errors go below/above this value it is clamped to this value.

Transition

YawRateAndSideslipController is designed so that it interprets the pilot pedal input as desired yaw rate at low speeds and as desired side slip at high speeds. A transition region based on speed can be defined to enable a smooth transition between these modes.

speedLow and speedHigh : Below speedLow, the pedal input will be interpreted as desired yaw rate and after speedHigh as desired side slip. Between these two speeds, the outputs of the Yaw Rate and Side slip subcontrollers will be linearly mixed.

Maximums

maxYawRate : Maximum yaw rate that can be requested from the controller. This parameter is in deg/s and must be positive always.

maxSideslip : Maximum side slip angle that can be requested from the controller. This parameter is in degrees and must be positive always.

InputMapping

nonlinearity : In some cases, it might be better to pass the pilot input to the control system non-linearly. For example, in a desktop simulator, where the input device is a common joystick, we might want to have a better precision at the zero point. We can define an exponential relationship between pilot's input and the control system input using nonlinearity parameter. The pedal input that is given to the control system is computed as follows: pedal input to controller = (|pilot pedal input| ^ nonlinearity) * sign(pilot pedal input).

PitchRateController, RollRateController

Controllers for pitch and roll rates.

printDebugInfo : A boolean parameter (possible values "true" and "false" ) to enable/disable debug output to the console from the controller. This is useful for the parametrization of the controller.

PID : Parameters of the PID controller that is used to regulate the rotational motion.

P : Proportional coefficient of the controller.

I : Integral coefficient of the controller.

D : Derivative coefficient of the controller.

offset : "Zero point" of the controller. When error, error derivative and error sum is zero resulting value will be offset.

IntegralMin and integralMax : Minimum/Maximum value the error sum can reach. If the sum of errors go below/above this value it is clamped to this value.

MaxPitchRate and MaxRollRate

These are maximum angular speeds the controllers are trying to reach. E.g. if control stick is pulled 100% right, RollRateController will try to reach a roll rate of maxRollRate. The values are in deg/sec.

InputMapping

nonlinearity : In some cases, it might be better to pass the pilot input to the control system non-linearly. For example, in a desktop simulator, where the input device is a common joystick, we might want to have a better precision at the zero point. We can define an exponential relationship between pilot's input and the control system input using nonlinearity parameter. The actual input that is given to the control system is computed as follows: input to controller = (|pilot input| ^ nonlinearity) * sign(pilot input).

Ground contacts

Interaction with terrain is simulated by ground contacts, classes that are derived from RTD::FGGroundContactBase. Highest fidelity ground contact implementation is RTD::FGGroundContact200 class in the current version. Initially gears are extended.

Ground contacts are defined in <GroundContacts> section of the configuration file. Each new ground contact has to be defined in its own section with name <GroundContact>. Also the type of the ground contact has to be declared here using the type attribute. In <GroundContacts> section <GeneralConfiguration> subsection is used to define some parameters that are common to all ground contacts. 

<GroundContacts>
	<GeneralConfiguration
		gearExtensionTime gearDragCoefficient
		gearReferenceArea maxSteerableGearDeflection
		steerableGearSpeedCorrectionFactor />
		<GroundContact type>
		<!-- Parameters Related to chosen ground contact model -->
	</GroundContact>
</GroundContacts>
GeneralConfiguration
  • gearExtensionTime: Time needed to extend/retract landing gears in seconds.
  • gearDragCoefficient: Coefficient of drag for landing gears.
  • gearReferenceArea: Reference area of the landing gears to be used in air drag computation.
  • maxSteerableGearDeflection: Maximum deflection (in degrees) of the steerable gears.
  • steerableGearSpeedCorrectionFactor: As the speed of the aircraft increases, deflections of steerable gears are automatically reduced to keep the aircraft controllable. This parameter defines the strength of this correction.
GroundContact

A multitude of <GroundContact> </GroundContact> sections can be defined. The type parameter defines which ground contact model (E.g. FGGroundContact200) has to be used. Inside the opening and closing tags, configuration parameters of the chosen model has to be provided.

FGGroundContact200 configuration

FGGroundContact200 implements a multibody dynamics based ground contact model. It supports inclined and moving surfaces (e.g. for ship deck landing). 

<GroundContact type="FGGroundContact200">
	<Features isBrakable isRetractable isSteerable />
	<Suspension>
		<MountStation x y z />
		<SliderAxis x y z />
		<SuspensionLength length />
		<SpringParameters
			springConstant
			damperConstant
			maxCompression />
	</Suspension>
	<Wheel mass radius>
		<RollAxis x y z />
		<Friction
			slidingFriction
			resistanceTorqueNoBrake
			resistanceTorqueFullBrake />
	</Wheel>
	<SteeringController
		maxTorque
		maxRotSpeed
		controllerGain
		inverseRotation
		keepCenter />
	<Damage thresholdLoad />
</GroundContact>
Features

General features of the ground contact.

  • isBrakable: Is the ground contact a brakable landing gear? (true or false)
  • isRetractable: Is the ground contact a retractable landing gear? (true or false)
  • isSteerable: Is the ground contact a steerable landing gear? (true or false)
Suspension

Parameters of the suspension.

  • MountStation: Position at which the landing gear is mounted onto the aircraft body. Relative to center of gravity and in meters.
  • SliderAxis: This is a unit direction vector which determines along which direction suspension slider can move. It is defined relative to aircraft body coordinates in RTDDynamics coordinate system convention.
  • SuspensionLength: Length of the suspension, i.e. the distance from mount station to middle point of wheel when suspensions are not compressed or decompressed. In meters.
  • SpringParameters
    • springConstant: Spring constant of the suspension spring. In N/m.
    • damperConstant: Damper constant of the suspension spring. In N*sec/meter.
    • maxCompression: Maximum amount the suspension spring can be compressed. This value should be less than or equal to suspension length. In meters.
Wheel

Parameters of the wheel.

  • mass: Mass in Kg.
  • radius: Radius of the wheel in meters.
  • rotationalInertia: Rotational inertia of the wheel in Kg*m^2.
    • RollAxis: A direction vector that defines the axis of wheel rotation. It is defined relative to aircraft body coordinates.
    • Friction
      • slidingFriction: Coefficient of friction between the tire and the ground.
      • resistanceTorqueNoBrake and resistanceTorqueFullBrake: An effective resistance torque from linear interpolation of resistanceTorqueNoBrake and resistanceTorqueFullBrake depending on the amount of applied wheel brake is applied to the wheel continuously to simulate rolling friction. The resistance torque acts in a direction that opposes the current rotation of the wheel. In N*m.
SteeringController
  • maxTorque: The steering controller rotates the wheel left or right depending on the pilot's input for steering on the ground. This parameter defines how much torque can be applied to the wheel for this purpose. This parameter is in N*m.
  • maxRotSpeed: This parameter defines how fast the wheel can rotate around the vertical axis. In Deg/sec.
  • controllerGain: The torque applied to the wheel to rotate it is proportional to the difference of pilot's desired wheel direction and current direction multiplied by \e m_steeringControllerGain (the result is clamped at \e m_steeringMaxTorque ). This parameter is unitless.
  • inverseRotation: If this parameter is set to true, the wheel will rotate in the opposite direction of the pilot input (if pilot input is to the right wheel will rotate to the left etc. ). This is necessary if the steering wheel is at the back instead of front.
  • keepCenter: If this parameter is set to true, the wheel will not be affected by the pilot's control input but instead will try to keep the wheel in center always. This is useful to model free rotating wheels (a wheel that can rotate around vertical axis but is not actively controlled by actuators ).
Damage
  • thresholdLoad: During landing the forces acting on the wheel are transferred to the fuselage. Beyond a certain threshold the transferred forces result in damage in the aircraft. The load defined in thresholdLoad (in N) is the threshold used for this purpose.
Retrieved from "https://community.bistudio.com/wiki?title=Helicopter_Flight_Model_Config_(XML)&oldid=349865"