Abstract

Articulating landing gear that use closed-loop feedback control are proven to expand the landing capabilities of rotorcraft on sloped and rough terrain. These systems are commonly referred to as robotic landing gear (RLG). Modern RLG systems have limitations for landing on dynamic platforms because their controllers do not incorporate fuselage roll and roll rate feedback. This work presents a proven crashworthy cable-driven RLG system for the commercial S-100 Camcopter that expands static landing zone limits by a factor of three and enables dynamic platform landings in rough sea state (SS) conditions. A new roll and foot-force feedback fused control algorithm is developed to enable ship deck landings of an RLG equipped S-100 without the need for deck lock or advanced vision-based landing systems. Multibody dynamic simulations of the aircraft, landing gear, and new control system show the benefits of this combined roll and force feedback approach. Results include experimental dynamic landings on platforms rolling under sinusoidal motion and simulated SS conditions. The experiments demonstrate, in a limited fashion, the usability of the RLG through ground experimentation, and the results are compared to simulations. Additional simulations of landings of the S-100 with rigid and active landing gear with more challenging landing conditions than experimentally tested are presented. Such results aid in understanding how RLG with this new roll and contact force fused controller prevent dynamic rollover.

1 Introduction

Throughout their history, rotorcraft have proven to be a vital tool for accessing landscapes that traditional fixed-wing aircraft are unable to reach due to their smaller landing zone footprint. Pilots of rotorcraft flight systems often face the complexity of landing on uneven, rugged, or moving surfaces in order to complete their mission. These complexities lead to operational limits for the landing conditions; pilots or flight control systems are allowed to operate in. For example, rotorcraft acceptable ground slope limits range from 6 deg to 15 deg, depending on their design (critical roll angle for dynamic rollover) and environmental conditions [1]. Additionally, pilot workload and safety risks increase for landings on mobile surfaces, such as maritime ship decks, because of the periodic, coupled landing surface dynamics [25]. Rotorcraft are critical to civil and military missions in maritime environments [68]. As such, there is a broad area of technological developments focused on improving mission safety.

The first set of technologies rely on modification to the maritime vessel itself. Deck lock mechanisms exist that lock the rotorcraft to a point on a landing deck in rough seas, see Fig. 1. The purpose of these tools is to ensure the rotorcraft does not bounce or slide during a landing. A drawback of these solutions is that they prevent the natural sliding of the rotorcraft and can cause damage to rotor mechanical components because of high lateral loads. Motion compensation ship decks have been conceptualized, but not commercially realized. These decks measure the maritime vessels' dynamic motion and adjust the landing deck's position and orientation accordingly [10]. Given the size and cost of this system, it has not been commercially viable to date.

Fig. 1
S-100 unmanned rotorcraft using a deck lock mechanism [9]
Fig. 1
S-100 unmanned rotorcraft using a deck lock mechanism [9]
Close modal

Improved training simulations, sensors suites, and flight controllers have all been explored to reduce workload and maximize safety margins associated with shipboard rotorcraft landings as well. Reber and Bernard developed a game-based learning, high-fidelity training simulation that removed the need for highly skilled instructor oversight [11]. Their research showed that trainees of their program operated with improved safety margins during normal and extreme ship deck operations when compared to a control group, showing the importance of human factors and training. Xu et al. [12] and Saripalli et al. [13] approached the problem through a combination of sensor suite and control algorithm development to implement computer vision driven landing guidance. Using a reference object of known size, shape, and orientation on the landing surface, they showed preliminary results that IR and color spectrum imagery can be used to land on a moving platform. This body of research did not consider or attempt to verify operational capability during high sea state (SS) conditions, limiting its current use in the most challenging conditions. Numerous researchers have approached the problem through unique deep reinforcement learning [14], nature inspired controllers [15], or optimal control techniques to path plan and land on a moving platform [1618]. These research topics showed increases in safety margins or reduced workload during approach, but none of them focus directly on dynamic rollover prevention or impact attenuation during the landing event. Robotic landing gear (RLG) for rotorcraft were developed as a solution to these challenges.

Robotic legs allow a given system to maneuver through irregular and rugged terrains, and they have been applied to a number of ground applications for quadrupedal and bipedal locomotion [1922]. Modern research and development into RLG for rotorcraft was reintroduced by Manivannan et al. [23], but the concepts and idea of RLG for rotorcraft have existed since the 1950s [24,25]. Manivannan's concept and theoretical design provided a replacement to skid gear with three or four articulating legs that enabled rotorcraft landings on sloped or rugged static surfaces. This RLG design concept was shown by Kiefer et al. [26] to reduce peak loads on rotorcraft occupants through active control of the legs so they act as adaptive shock absorbers. An experimental flight vehicle based on Manivannan's work demonstrated, in a limited fashion, the viability and capability of RLG for rotorcraft on a 120 kg unmanned helicopter using a four-legged design [27]. A next generation crashworthy RLG designed by Di Leo et al. [28,29] implements a cable-driven, four-bar linkage mechanism. This mechanism provides single degree-of-freedom motion while housing all actuation electronics inside the fuselage and minimizes weight through the use of composite materials. León et al. presented an integration and flight tests of this next generation design on an S-100 Camcopter with a maximum takeoff weight of 200 kg, see Fig. 2 [30]. Similar concepts designed and experimentally verified by Stolz et al. [31] using a 78 kg helicopter and Huang et al. [32] on a 20 kg octocopter are additional examples of articulating landing gear concepts. The works by Manivannan et al., Stolz et al., and León et al. only used force feedback to control leg orientation during the landing, while Huang's work used infrared sensors on the bottom of each landing leg. The use of force feedback or foot-to-ground distance measurements is not sufficient for shipboard rotorcraft landings because the feedback saturates and is no longer meaningful once each leg makes contact and the rotorcraft drops collective.

Fig. 2
S-100 Camcopter with cable-driven, four-bar RLG prior to flight tests
Fig. 2
S-100 Camcopter with cable-driven, four-bar RLG prior to flight tests
Close modal

For any system where the landing platform moves before and after all legs make contact, there is a need for closed-loop feedback channel(s) in addition to force/foot–ground distance information. Kim and Costello [33] developed one solution with a theoretical study on virtual model control for four-legged RLG landings on naval ship decks. This controller uses feedback from the rotorcraft orientation states to optimally compute the minimum actuator torques necessary to maintain a level fuselage. Their simulations found that such a controller in conjunction with four-legged RLG grasping mechanisms on each leg will absorb impact and conform to ship decks in rough seas with mean wave heights from 4 to 6 m and modal wave periods of 9–17 s, which is also referred to as sea state 6. Kim and Costello's work is promising but it did not provide experimental validation, and it may suffer from complications due to noise and bias in measurements or actuator response of the physical system [34]. Even with this broad scope of research, there is not at present an experimentally validated, robust RLG system for ship deck landings.

This work presents the design and experimental ground tests of roll and force feedback controlled cable-driven, four-bar link mechanism RLG for dynamic platform landings with an emphasis on the ship deck landing problem. Integration of the controller for this work is on the cable-driven, four-bar link RLG system previously used on the S-100 Camcopter during flight tests as presented by León et al. [30]. Section 2 is an overview of the system level design, including RLG kinematics, state estimation, and control for the S-100 platform. Section 3 details the multibody dynamics simulation tool used to guide the design and performance characterization of the S-100 RLG platform in a variety of static and mobile landing deck scenarios. Finally, experimental setup and results are presented and compared to simulations in Sec. 4. Section 5 concludes the paper.

2 System Design

This section details the subsystems pertinent to the cable-driven, four-bar RLG system as integrated with the S-100 Camcopter. These subsystems include the structural mechanism, sensing and state estimation, and multisensor feedback controller.

2.1 Cable-Driven, Four-Bar Linkage Robotic Landing Gear.

Figure 3(a) illustrates a wireframe diagram of the S-100 cable-driven, four-bar RLG in the neutral configuration. A continuous cable (marked in red), gas spring, and actuation/spool are the three primary components that enable motion. The gas spring maintains tension in the cable by applying an upward force (i.e., retracts the legs upward). The cable acts as the opposing force to move legs downward, and it is routed through a spool such that the spool's rotation changes the effective ground angle absorbed by the legs, see Fig. 3(b). This cable doubles back around the crank-leg pin joint and terminates within the fuselage on a shock absorber. The total cable length determines the neutral stance and the maximum achievable ground angle, as studied by Di Leo et al. [29].

Fig. 3
(a) Concept diagram of the cable-driven, four-bar link mechanism RLG and (b) differential motion of the RLG concept via rotation of the cable's center
Fig. 3
(a) Concept diagram of the cable-driven, four-bar link mechanism RLG and (b) differential motion of the RLG concept via rotation of the cable's center
Close modal

The authors chose to implement single degree-of-freedom differential motion based on (1) the normal flight operations of rotorcraft on sloped surfaces, (2) the dynamics of a rolling and pitching ship deck, and (3) the specific configuration of the S-100 landing gear. Normal sloped surface operation suggested by the Federal Aviation Administration (FAA) to prevent dynamic rollover [1] recommends landing on a slope aligned with the roll axis of the rotorcraft.

In addition to this recommendation by the FAA, representative ship deck motion up through sea state 6 does not require rotorcraft pitch control which would require a symmetric retraction/extension of the legs. Illustrated in Fig. 4(a) is a representative maritime vessel, the DDG-54 Arleigh-Burke class destroyer with a length of 154 m, a width of 18.8 m, and a draft of 7 m. These physical properties were used to generate the sea state 6 roll and pitch data in Fig. 4(b) for a period of 1 min. This simulation, described further in Sec. 3.2, used a modal wave period of 17 s and a significant wave height of 6 m with the ship moving at 10 kn. This simulation produced roll on the order of ±12 deg on the landing deck and a pitch of approximately ±5 deg. Even for passive landing gear, ±5 deg of pitch is well within most rotorcraft operational limits, but the roll is above ±10 deg repeatedly. Between the FAA recommendation, and ship deck dynamics, any rotorcraft operating on sloped or mobile surfaces would see a large expansion in mission capabilities with a single degree-of-freedom aligned with the aircraft's roll axis. In addition to these general observations, there are two S-100 specific factors that limit the need for pitch-axis control: (1) the S-100 has a shock-absorber tail and (2) the forward landing leg contact points are lower than the rear leg of the S-100 (i.e., the forward legs will make contact first for small pitch angles of the landing surface).

Fig. 4
(a) Arleigh-Burke class naval vessel [35] and (b) its roll and pitch motion from a sea state 6 simulated motion
Fig. 4
(a) Arleigh-Burke class naval vessel [35] and (b) its roll and pitch motion from a sea state 6 simulated motion
Close modal

Each of the provided considerations shows that there are clear advantages to having RLG for control along the rotorcraft's roll axis. Therefore, the authors chose to focus on actuation along the roll axis and left actuation along the pitch axis to future work, since pitch control may be valuable for other rotorcraft. A more detailed analysis of the design constraints and crashworthiness (up to 5 m/s) of the S-100 design of a cable-driven, four-bar RLG (see Fig. 2) was initially presented by Di Leo et al. [28,29] and expanded upon by León et al. [30]. The same RLG flight tested by León et al. is the system used throughout the remainder of this work.

Kinematic relationships are useful to understand the S-100 RLG's ideal performance limitations and how that relates to ship deck dynamics. Specifically, the kinematic relationship between the actuator speed and the rate at which the legs absorb ground angle is of interest. To start, we define the key parameters used in deriving this relationship. Figures 3(a) and 3(b) defined the ground angle, γ, the leg angle relative to the fuselage, θ, distance between the follower-fuselage pin and the crank-leg pins, s, crank length, L1, and the distance between the crank and follower components, L2. Subscripts of L or R refer to the measurements on the left or right half of the RLG assembly, respectively. The first kinematic relationship we consider is the mapping of leg angle relative to the fuselage θ to the distance s, which is the distance between the follower-fuselage and crank-leg pin joints. This relationship for the left and right side of the assembly is formulated as
(1)
(2)
by applying the law of cosines to the four-bar link mechanism. These two equations may be rearranged and generalized as
(3)
Equation (3) now maps the leg angle relative to the fuselage (as defined in Fig. 3(b)) to the distance between the follower-fuselage pin and the crank-leg pin of the four-bar linkage. This is a useful kinematic relation because the angles θL/R are simple to measure with sensors. The time derivative of Eq. (3) is
(4)
where s˙ is the time rate of change of s. Since the cable is continuous through the spool, s˙=s˙L=s˙R. The rotation rate at which the spool pulls-in or releases cable is β˙. Cable speed is denoted by s˙ and is related to spool rotation rate by s˙=β˙rspool/2. Note that s˙ is a function of β˙/2 because the actuation cable routes around the crank-leg joint and doubles back to the fuselage where it is secured to a shock absorber. Incorporating this relationship into Eq. (4) yields
(5)

Equation (5) now relates the rotation rate of the legs relative to the fuselage and the actuator's rotation speed.

Now we relate the legs' rotational velocity θ˙L/R to the maximum ground angle rate of change γ˙. The ground angle γ is computed as
(6)
(7)
The time derivate of γ is
(8)
(9)

These relationships give useful insight when input into the relationship for ground angle absorption rate, γ˙, where B is the horizontal distance between the legs, and Δh˙ is the descent velocity that the legs are able to track. B is presumed constant. This assumption is acceptable given that B changes only by a small amount, 3%, over the gear's entire range of motion. Equations (5)(9) and actuator performance (β˙) are combined to yield kinematic limits for the S-100 RLG such as maximum achievable ground angle γmax and its associated rate γ˙max. Table 1 lists the necessary supporting specifications to compute the maximum ground roll angle and ground roll angle absorption rate. The S-100 RLG's unique structure constrains the motion such that the maximum roll angle absorption rate, γ˙max, changes as a function of ground roll angle, γmax.

Table 1

S-100 RLG kinematic parameter values

VariableValueVariableValue
L10.305 mL20.102 m
θ080 degs00.304 m
θmax120 degθmin44 deg
B1.27 mrspool0.0168 m
β˙max14 rad/s
VariableValueVariableValue
L10.305 mL20.102 m
θ080 degs00.304 m
θmax120 degθmin44 deg
B1.27 mrspool0.0168 m
β˙max14 rad/s

Figure 5 shows this limit boundary up to the maximum achievable ground roll angle. The RLG configuration are able to conform up to 17 deg with a maximum leg velocity of 0.69 m/s. Maximum leg velocity, which stems from β˙max, is fundamentally limited by the S-100's onboard available power and the original design constraints from Refs. [2830]. Actuator performance was experimentally measured, so it inherently incorporates cable friction and gas damper force and viscous damping. This maximum leg velocity does not change substantially across the system's range of motion, but the geometry causes a change in the maximum ground angle absorption rate, see Fig. 5. Angular roll rates consistent with the kinematic equations and computed results were achieved by León et al. [30] during their ground and flight test operations with this system. The x's in the plot show a sample of the performance achieved by the S-100 RLG during previous ground and flight tests. León et al. also showed that landings at 0.5 m/s were possible with minimal perturbations on roll angle on slopes up to 15 deg (3× the S-100's safe operational roll angle), but landings faster than 0.5 m/s inevitably have some notable roll perturbations because the legs are physically limited in acceleration and maximum speed.

Fig. 5
Operational limits of the S-100 RLG compared to relevant experimental data from sea state 6 simulations
Fig. 5
Operational limits of the S-100 RLG compared to relevant experimental data from sea state 6 simulations
Close modal

Overlaid on Fig. 5 is a region of roll angles and angular rates from a sea state 6 simulation with a significant margin between the peak performance limit of the RLG compared to expected landing deck conditions. Based on the kinematic limits alone, the design of the S-100 RLG from León et al. appears have adequate performance for dynamic platform landing on the deck of the DDG-54, a modern maritime vessel. Note that this overlay only considers the rolling rate of the deck. It does not incorporate the combined ranges of aircraft descent rate and deck heave velocity. The following section detail the control algorithm used to conduct dynamic platform landings.

2.2 Feedback Controller.

Previous research on experimental testing of RLG systems implemented force feedback control using ground contact sensors with no fuselage orientation feedback during or after the landing [23,27]. The primary assumptions of those systems are (1) that the closed-loop control will react faster than the aircraft descends and (2) the aircraft holds a steady descent rate throughout the landing. Both are idealizations of the problem. The first assumption means that the legs will not cause the aircraft to roll or pitch throughout the landing. This assumption is acceptable for low descent rates while landing on static surfaces with ideal force sensors and high update rate feedback control because the closed-loop system reacts faster than the aircraft descends. This is not a valid assumption during rapid landings or on helicopters with unique landing gear configurations (see Ref. [30]). The second assumption is true for static landing surfaces, but it cannot be true for landings on dynamic platforms. The orientation and position of the landing surface results in variable velocity between the aircraft and landing surface. For this reason, the second assumption used on force feedback controlled landing gear is not valid for landing conditions considered in this work. The lack of closed-loop roll and/or pitch feedback is a limitation of contact feedback RLG controllers.

This work fuses roll feedback with previously tested foot-force feedback loop to eliminate the assumptions of previous works. The roll and force feedback control designed here follows three operational stages illustrated in Fig. 6. The first stage has the aircraft descend with no contact on either foot. Regardless of the roll angle during the descent phase, the legs do not actuate. Once one leg makes contact with the landing surface, the second phase initiates. The RLG controller incorporates roll angle and force feedback to minimize dispersion in inertial roll angle until both legs make contact. The force sensor measurement range on both legs limits the effectiveness of force data after landing; thus, only roll angle feedback is used to reach a safe roll angle in the final stage. Nominally, this last phase ends when the aircraft reaches a minimal roll angle within an acceptable deadband, ϕDB, on a static surface or when the frame is at a minimal roll angle and the legs are in their neutral position. The included roll feedback allows the swift detection of whether or not the landing surface is in motion.

Fig. 6
Diagram illustrating three phases of control using RLG with available roll and force feedback
Fig. 6
Diagram illustrating three phases of control using RLG with available roll and force feedback
Close modal

Figure 7 shows the high-level block diagram of the combined roll and force feedback controller used in this work. There are two control loops using three feedback signals. The outer loop control combines foot-force data and fuselage roll data to command a desired leg angle, then the inner loop control monitors and changes the motor's speed, β˙, to achieve the desired leg angle. The leg-angle estimator and foot-force feedback estimators were previously presented by León et al. [30], and the roll estimate used for loop closure is provided by an attitude and heading reference system sensor (AHRS). The RLG control starts with a conditional proportional–derivative (PD) controller that uses desired and estimated foot forces in addition to the desired and estimated roll angles as inputs. A lower-level outline of the conditional PD control operates and calculations are outlined in Algorithm 1.

Fig. 7
Roll sensor fused RLG control algorithm block diagram
Fig. 7
Roll sensor fused RLG control algorithm block diagram
Close modal
Algorithm 1

RLG conditional PD control

1: Control interrupt driven by F̂R/L>Fd begins loop
2: repeat
3:   ifF̂L>Fd and F̂R<Fdthen
4:    eFk=FdF̂Lk
5:   else ifF̂R>Fd and F̂L<Fdthen
6:    eFk=F̂RkFd
7:   else [F̂R>Fd and F̂L>Fd]
8:    eFk=0
9:   end if
10:
11:   CHECK F̂L and F̂R FOR BOUNCE
12:
13:   if|ϕ̂|>ϕDBthen
14:    eϕk=ϕdϕ̂
15:   else [|ϕ̂|<ϕDB]
16:    eϕk=0
17:   end if
18:
19:   δθkF=KPeFk+KD(αFe˙Fk+(1αF)e˙Fk1)
20:   δθkϕ=KPeϕk+KD(αϕe˙ϕk+(1αϕ)e˙ϕk1)
21:   δθk=δθkϕ+δθkF
22:   θck=θ̂Lk+δθk
Ensure: θck<θmax and θck>θmin
23:
24:   eθk=θckθ̂Lkδθk
25:   β˙k=β˙0+KPeθk+KD(αθe˙θk+(1αθ)e˙θk1)
Ensure:  |β˙k|<β˙max
26: until|ϕ̂|<ϕDB and θ̂L/Rkθ0
1: Control interrupt driven by F̂R/L>Fd begins loop
2: repeat
3:   ifF̂L>Fd and F̂R<Fdthen
4:    eFk=FdF̂Lk
5:   else ifF̂R>Fd and F̂L<Fdthen
6:    eFk=F̂RkFd
7:   else [F̂R>Fd and F̂L>Fd]
8:    eFk=0
9:   end if
10:
11:   CHECK F̂L and F̂R FOR BOUNCE
12:
13:   if|ϕ̂|>ϕDBthen
14:    eϕk=ϕdϕ̂
15:   else [|ϕ̂|<ϕDB]
16:    eϕk=0
17:   end if
18:
19:   δθkF=KPeFk+KD(αFe˙Fk+(1αF)e˙Fk1)
20:   δθkϕ=KPeϕk+KD(αϕe˙ϕk+(1αϕ)e˙ϕk1)
21:   δθk=δθkϕ+δθkF
22:   θck=θ̂Lk+δθk
Ensure: θck<θmax and θck>θmin
23:
24:   eθk=θckθ̂Lkδθk
25:   β˙k=β˙0+KPeθk+KD(αθe˙θk+(1αθ)e˙θk1)
Ensure:  |β˙k|<β˙max
26: until|ϕ̂|<ϕDB and θ̂L/Rkθ0

Algorithm 1 generally mimics the procedures of the flight software, but it does not include specific logic of health monitor checks, control value clamping, and additional sanity checks on control outputs that exist in the flight software. In this algorithm, eFk,eθk, and eϕk are the errors in foot force, leg angle, and fuselage roll angle, respectively. The desired force for the legs to keep in ground contact is Fd, and the desired fuselage roll angle is ϕd. The desired roll angle, ϕd, is always 0deg, so the controller tries to maintain a level fuselage.

The initial conditional PD algorithm is interrupt driven upon initial contact on either foot, so asynchronous inputs do not impede performance. The forces measured by the left or right foot that initiates control are FL or FR, respectively. The loop runs continuously after a force sensor crosses the predetermined desired force Fd. Force sensors designed and characterized in Ref. [36] were configured in an array on the bottom of each leg. The force threshold and how it was chosen for the S-100 implementation were discussed by León et al. [30] in more detail, but is suffices for this work to say the threshold is 5 N. Line 22 of the algorithm denotes that the primary leg measurement used to track leg-angle control is the left leg. Therefore, the force error force computations (lines 4 and 6) are different depending on whether the left or right food makes contact first, so that the leg-angle control computations drive the legs in the proper direction. Note that the system measures the left and right leg angles, and it can dynamically switch between the left and right leg as the primary leg-angle measurement for feedback control based on sensor health or loss of communication.

The derivative control equations implement a first-order, discrete low-pass filter on the derivative control to limit noise from the incoming signals [37]. Each channel of control (roll and foot force) has their own unique low-pass filter parameter, α, which is chosen based on the frequency characteristics of each input channel. There are more elegant solutions for the filtering of this data, but the authors found the first-order, low-pass filter to be effective for this application.

Notable from this algorithm is that the controller uses force feedback until the aircraft exceeds the allowable deadband along the roll axis. For this reason, it is advantageous to minimize the deadband of the system. Methods to minimize deadband include high dynamic accuracy roll measurements and/or actuation methods that can run continuously with little or no regard for power consumption. For this system, a VectorNAV VN-100 was used for roll and pitch feedback. The company lists the sensor's static roll accuracy to 0.5 deg and dynamic roll accuracy to 1.0 deg, so it is suitable for testbed and experimental use but a higher accuracy sensor would be ideal for a robust, fielded system. This dynamic roll accuracy was used as the roll deadband, ϕDB, value for the experiments in work.

While the leg-angle control does not have a feed forward term, the motor control does. This feed forward term, β˙0, is the required minimum spool rotation rate command to overcome static friction of the system. This algorithm succinctly fuses the roll angle feedback, and with a minor change to the δθkϕ equation it also can fuse roll rate, p, feedback. This minor change takes on the form
(10)

and the use of this alternate computation depends on the quality of roll rate feedback from an AHRS.

3 Simulation and Experimental Setup

This section elucidates the simulation tools and their integration as a performance prediction model for the S-100 RLG. First, a multibody dynamics tool is described that has a proven track record for RLG development. Then, a ship deck motion simulation tool and standardized sea state conditions that are commonly used are described.

3.1 Rigid Body Dynamics Simulation Tool.

This work implements a modified multibody dynamics tool similar to those used by Leylek et al. [38], Gross et al. [39], and Kim and Costello [33]. This same tool was used to size the S-100 RLG for crashworthiness in the work of Di Leo et al. [29]. In this work, modifications are made to characterize the control and sensing aspects of the dynamic platform landing sensing, state estimation, and control. All of the modifications and improvements rely on experimental feedback from earlier works [29,30,36] to produce a useful simulation tool for analyzing landing scenarios outside the experimental results presented in this work. The simulation tool starts with a general form of multibody dynamic equations of motion
(11)
In this relation, X¯ is the state vector containing states from any body or contact model, F¯ is a vector containing the nominal Newton–Euler equations of motion for a six degrees-of-freedom (6DOF) rigid body without constraints, and U¯ is a vector of constraint forces and moments that are mapped through the linear operator [G] appropriately to each body. The form factor of this equation presents an opportunity to apply a constraint stabilization controller to compute U¯. The set of translation and rotation constraints, E¯, are derived for the joints and body interactions of the given dynamic system then reconfigured such that
(12)
The first and second derivatives of this constraint relation yield
(13)
(14)
(15)
(16)
Feedback linearization is then applied to the dynamics of the constraint equation, Eq. (14), to compute the reaction forces and moments in U¯ as
(17)

Proper selection of the damping ratio, ζ, and the natural frequency, ωn, provide a stable constraint stabilization controller. The optimal values of ζ and ωn depend on the simulation time-step and dynamics or vibration frequencies of the system, but a wide range of values maintain constraints and system stability.

This tool has a number of subsystem models that can be included or neglected in the force and moment calculations of any body in a modular manner. Figure 8 illustrates how the subsystem models are defined, then incorporated into the computation of F¯. This figure also calls out the three models used for this work, which includes gravity, ground contact, and a first-order rotor thrust model.

Fig. 8
Diagram illustrating where subsystem models and user-defined sensing and control are incorporated into the multibody dynamics computation loop
Fig. 8
Diagram illustrating where subsystem models and user-defined sensing and control are incorporated into the multibody dynamics computation loop
Close modal
Ground contact is modeled as a standard linear solid (SLS) with friction. This model mimics ideal viscoelastic behavior with a single relaxation time of the landing surface material [40,41]. The structure of the contact model, as seen in Fig. 9, requires the force derivative at any given time to be computed. This force derivative is computed as
(18)

for each leg that is in contact. A coefficient of friction, static and kinetic, is applied to the normal force computation for tangential sliding force. While not a perfect assumption of contact, it is a suitable and fast model of contact for RLG contact and feedback control. In this work, we are interested in dynamic landing surfaces as well. A simple extension to the SLS model is used for dynamic ground contact where δ and δ˙ are computed from the relative velocity of a foot and the ground at any given point. The nominal values of the spring constants and damping coefficients are provided in Table 2. These ground contact parameter values were tuned by experimentation with drop tests of the experimental landing gear with an S-100 airframe surrogate in Refs. [28] and [29].

Fig. 9
Diagram of the SLS computation model shown with a force normal to the plane of contact for the case of a static landing surface
Fig. 9
Diagram of the SLS computation model shown with a force normal to the plane of contact for the case of a static landing surface
Close modal
Table 2

List of SLS ground contact parameters used during simulations

ParameterNominal value
kv100,000 N/m
k50,000 N/m
c5000 N s/m
ParameterNominal value
kv100,000 N/m
k50,000 N/m
c5000 N s/m
A first-order rotor thrust model is pertinent to this work because the drop in collective of the S-100's rotors leads to a transient period of load reduction in which dynamic rollover is possible. Rotor thrust has not been modeled for dynamic simulations of the S-100 integrated system prior to this work. This transient period of load reduction when the collective drop is initiated by a pilot or flight control system puts the rotorcraft at risk of dynamic rollover on the moving landing surface. The model assumes that the rotor collective and cyclic react as first-order models on the commanded inputs. Prior work by Kim and Costello [33] detailed this rotor model in the context of a pilot control model. With uncertainty in how pilot control inputs are mapped or filtered on the S-100, the pilot model was neglected here, and a constant rotor time constant of 0.5 s was used. The complete rotor model then follows as
(19)
(20)
(21)

where T is the rotor thrust, ηlong is the thrust angle along the rotorcraft longitudinal plane, and ηlat is the thrust angle along the rotorcraft lateral plane. The angles ηlong and ηlat stem from cyclic input to the rotor, but are not actively controlled within this work's simulations. The parameters τT and τη are the time constants for the first-order responses of the rotor thrust and thrust angles. Finally, Tc,ηlongc,andηlatc are the desired rotor trust, desired thrust angle along the rotorcraft longitudinal plane, and the desired thrust angle along the rotorcraft lateral plane, respectively.

There are 15 rigid bodies used for dynamic simulation and constraint stabilization using the aforementioned framework, see Fig. 10. While most bodies are direct counterparts to their physical counterparts, such as the fuselage, BD1, and rotor, BD15, other parts are decomposed into multiple rigid bodies in order to introduce elastic joints. These are necessary for understanding all coupled dynamic interactions such as pitch rate to roll inertial coupling of the rotor to the system [4244]. For example, the cable on the left leg assembly (right side of the image) is represented by two bodies, BD10 and BD11. These bodies have variable length to mimic changes in rope length resultant of spool rotations. More important than the ability to vary geometry, these bodies are connected by an elastic joint, denoted by a black dot in Fig. 10. This allows us to capture some of the elastic behavior of the real bodies. Similarly, the carbon fiber reinforced legs are represented by BD4 and BD13 as well as BD7 and BD14 with elastic joints in between. The spring stiffness and damping of each elastic joint (marked by black dots) were experimentally derived from drop tests of the full-scale S-100 RLG [28,29]. The rear leg of the aircraft, BD12, also connects to the fuselage via an elastic joint, though not shown in the figure. Joint elastic parameters for the rear leg were directly provided by the aircraft supplier, but were not experimentally verified. Finally, BD15 rotates with the S-100's rotor nominal frequency, and all thrust model outputs are applied to this body's center of gravity.

Fig. 10
Diagram of the 15 connected bodies used to simulate the S-100 and RLG
Fig. 10
Diagram of the 15 connected bodies used to simulate the S-100 and RLG
Close modal

At this time, aerodynamics and ship wake are not modeled in these simulations. The simulations considered in this work assume the rotorcraft approach is complete, and the simulations initiate just prior to first contact of the gear during landing. Other works have studied aircraft (, namely, rotor) feedback [45,46] and path planning [17] through the ship's wake field. This work focuses on the landing problem when a rotorcraft may enter a dangerous dynamic rollover scenario caused by aggressive roll and roll rate of a ship deck.

The authors used the simulation tool to characterize roll dynamics of S-100 RLG in a rigid, locked configuration. This configuration is the comparison baseline for improvements by the active gear. Table 3 lists the variables for the Monte Carlo. The variable μf is the friction with the landing surface, u is the aircraft's body forward velocity, w is the aircraft's body descent velocity, and θAC is the aircraft's pitch angle. Each range of values considered is within the aircraft's expected landing envelope, and values simulated were generated using a uniform distribution.

Table 3

List of Monte Carlo variables for locked S-100 RLG dynamic rollover characterization

ParameterNominalRange simulated
τT0.75 s[0.25, 1.25] s
μf0.7[0.25, 0.75]
ϕ−2 deg[−3, 0] deg
θAC0 deg[−1, 1] deg
u0 m/s[−0.1, 0.1] m/s
w0.5–0.6 m/s[0, 1] m/s
γ =±5 deg[−15, 15] deg
ParameterNominalRange simulated
τT0.75 s[0.25, 1.25] s
μf0.7[0.25, 0.75]
ϕ−2 deg[−3, 0] deg
θAC0 deg[−1, 1] deg
u0 m/s[−0.1, 0.1] m/s
w0.5–0.6 m/s[0, 1] m/s
γ =±5 deg[−15, 15] deg

The simulations all initiate with one foot in contact with the ground surface, and the commanded thrust changes to the minimum value. This is designed to mimic how pilots land the aircraft, where they initiate contact then lower collective to the minimum value. Figure 11 shows cumulative distribution functions (CDFs) of the peak fuselage roll angle absolute value of eight landing surfaces, γ=±5deg,7deg,10deg,and15deg, from this study. The x-axis is limited to 20 deg to show the fine detail for each slope, any portion of the CDF beyond 20 deg represents dynamic rollover. The results of the study show that the S-100 RLG configuration with the gear inactive, or locked, shows no risk with landings at γ=5deg. The risk for dynamic rollover first appears at γ=7deg. The simulations suggest that less than 5% of landings in the nominal envelope on a 7 deg slope result in dynamic rollover. By γ=10deg, the risk increases with more than 40% of landings in an otherwise nominal landing envelope result in dynamic rollover, while γ=15deg results in more than 85% of landings with dynamic rollover. These results show that the S-100 would benefit from active gear on static surfaces at low as 7 deg. The authors note that these CDF plots will shift downward on dynamic platforms with similar γ because of the added interface dynamics between the aircraft and landing surface.

Fig. 11
Peak roll angle cumulative distribution of Monte Carlo simulations of the locked S-100 RLG landing on static landing surfaces within nominal operation envelope
Fig. 11
Peak roll angle cumulative distribution of Monte Carlo simulations of the locked S-100 RLG landing on static landing surfaces within nominal operation envelope
Close modal

3.2 Dynamic Platform Simulation Tools.

Three types of landing surface dynamic motion programs were simulated and experimented with in this work. The first is a static surface. This is the most basic landing surface for RLG to interact with and has been extensively simulated and tested (see Fig. 11 and Refs. [29] and [30]). This section describes the sinusoidal and maritime ship deck landing platforms and how their motion is generated for simulation and experiments.

3.2.1 Sinusoidal Platform Dynamics.

The S-100 RLG has significant controllability along the aircraft roll axis, so all platforms undergo roll motion in this work. Roll motion simulations and experiments provide data of the system's performance decoupled from other motion (pitch rate and translational velocities). Roll sinusoidal motion is described by
(22)
where A is the desired sine wave amplitude, f is the frequency, t is the current simulation/experiment time, and t0 is the time shift required to start with a desired roll angle. This assumes the oscillation always occurs around a 0 deg center point. In the simulation tool, a user defines the desired amplitude, A, initial roll angle, ϕ0, and the maximum angular rate desired, ϕ˙r. The tool determines t0 as the first time where the sinusoid is the desired roll angle from the amplitude and initial roll angle inputs. Additionally, the frequency is computed as
(23)

A sine wave is generated for the simulation using the resultant frequency, f, and time shift, t0, and input into the multibody dynamic simulation tool.

3.2.2 Ship Deck Dynamics.

A ship deck simulation tool provides realistic landing deck 6DOF motion, given the vessel's geometry, mass properties, and sea state parameters. This simulation tool is crucial to understanding the benefits of a roll-fused RLG controller because rotorcraft often must land on dynamic maritime ship decks. The ship motion simulation consists of two parts, the Ship Motion Program (SMP) [47] and the Simulation Time History (STH) [48]. Given a ship with specified dimensions and hull characteristics, the SMP creates 6DOF response transfer functions for regular waves. The STH uses the transfer functions from the SMP to predict the ship motion due to ocean waves. This ship motion simulation generates deck motion as a function of wave motion parameters.

The origin of the ship is located at the intersection of the ship's forward perpendicular, centerline, and baseline. This point is easily identifiable and does not vary with changes in ship draft or trim. The direction of the x, y, and z coordinates from the ship's origin is depicted in Fig. 12. Inputs to the simulation are the hull characteristics (ship length, draft, and trim). The ship configuration used for the simulation shown here is the same as defined in Sec. 2. The SMP generates the ship's 6DOF response transfer functions for regular waves. It also provides absolute and relative 6DOF information (surge, sway, heave, roll, pitch, and yaw) of the ship at the origin and other desired points on the ship. The STH provides realistic, random wave time histories of 6DOF ship responses using the regular wave response transfer functions obtained from the SMP. Random waves are simulated using a two parameter Bretschneider wave spectral model to define the frequency content of the random sea waves [49]. The two parameters used are the significant wave height and the modal (peak-to-peak) wave period. Strip theory is used to obtain the response for a ship advancing at constant forward speed with arbitrary heading in regular sinusoidal waves. The Access Time History Computer Program then uses ship origin responses from the STH to generate response time histories at other locations on the ship. The inputs to the STH are the transfer functions for a particular ship generated by the SMP, sea state (significant wave height and mean wave frequency) and the steaming condition (ship heading relative to predominant wave direction, ship speed). The outputs from the STH used in this work are 6DOF time histories for motion of the helicopter landing surface of the ship.

Fig. 12
SMP and STH simulation tool center of gravity location and motion definitions
Fig. 12
SMP and STH simulation tool center of gravity location and motion definitions
Close modal

Nominal ranges of wave height and wave modal period are correlated to standard SS codes. Table 4 lists the sea state parameters from SS 3 through SS 6. These values come from the nominal conditions list by Bales [50] as measured in the North Pacific. Simulations in this work use the largest wave height and average wave modal period. An example of roll and pitch data using SS 6 wave conditions with STH was shown in Fig. 4(b). Sections 4.2.14.2.3 use a combination of static and dynamic landing surfaces to better understand the theoretical limits of the S-100 RLG.

Table 4

Sea state parameters used for simulations and ground tests

Sea stateWave height (m)Modal period (s)
30.50–1.255.20–15.5 (7.50 mean)
52.50–4.007.20–16.4 (9.70 mean)
64.00–6.009.30–16.5 (13.8 mean)
Sea stateWave height (m)Modal period (s)
30.50–1.255.20–15.5 (7.50 mean)
52.50–4.007.20–16.4 (9.70 mean)
64.00–6.009.30–16.5 (13.8 mean)

3.3 Experimental Test Setup.

Experiments of the integrated system were performed at Boeing's Mesa, AZ facility with a Sarnicola Simulation Systems hexad series 6DOF hydraulic motion table. This S-100 RLG system uses the same drivetrain as designed and flight tested by León et al. [30]. Figure 13 illustrates the integrated S-100 RLG on the motion table. Denoted by arrows are the shore power cable that was always connected during RLG operation on the aircraft and a crane connection point.

Fig. 13
S-100 RLG on the Sarnicola motion table for experimental tests of roll-fused RLG control
Fig. 13
S-100 RLG on the Sarnicola motion table for experimental tests of roll-fused RLG control
Close modal

The crane used for experiments lowers the system at a constant 0.18 m/s for these experiments. The orientation of the aircraft during crane descent is approximately level, unlike the flight system that trims at −2 deg. The Sarnicola motion table can maintain angular rates in excess of 30 deg/s, but has hard limits at ±14 in. of heave (vertical) motion [51]. This limitation meant experiments could not fully mimic ship deck surge–sway–heave conditions, but experiments were conducted up to the acceptable limits of the table for translation motion. Operation of the configuration was limited to 10 deg/s in roll, pitch, and yaw rate, and SS 5 with the restriction to surge, sway, and heave to 10% of maximum values. The experimental setup as provided does not mimic the full range of landing speeds, and it does not include dynamic interactions from spinning rotor blades, so this experimental system will perform better than the flight system. Nevertheless, the experiments prove necessary to compare with the results of the simulation tools, so that the entire landing envelope may be further explored with the simulation tools.

4 Results and Discussion

Two categories of dynamic platform landings were tested as a part of this work; sinusoidal and sea state. Sinusoidal table motion required a wave frequency and amplitude, but provides a consistent and decoupled input for the RLG to track. These tests were completed to better understand the hardware's physical limitations before experimental tests on sea state platform landings.

4.1 Sinusoidal Platform Experiments.

Sinusoidal landing experiments with roll amplitude of ±10 deg and peak roll rates between 3 and 5 deg/s were conducted. These experiments were conducted to measure pure roll rate response of the RLG decoupled from pitch motion and translational velocities of the landing surface because the roll is axis where this RLG was designed to mitigate dispersions. These experiments began with the S-100 on the table, then the table activated. With the aircraft on the table, the force feedback portion of control is effectively turned off from the beginning. Numerous experiments were conducted prior to the examples presented in this section to converge on a smooth and fast response of the RLG. Figure 14 shows still images of the experiment (a), still images of the simulations (b), and experimental data (c) and (d), of a motion table commanded ±10 deg, 3 deg/s peak roll rate sinusoid as an example of these tests. The images in Figs. 14(a) and 14(b) illustrate the RLG and fuselage undergoing roll feedback control and maintaining the fuselage within or at the bounds of ϕDB for the experiment and simulation, respectively. As expected from simulations and kinematic relations, the feedback controller and drivetrain are consistently able to absorb roll dynamics of the landing platform. Experiments were run for less than 30 s to prevent any motor overheating issues.

Fig. 14
(a) S-100 RLG response to a ±10 deg sinusoid with peak roll rates of ±3 deg/s where the aircraft starts the experiment on the table, (b) still images of a simulation under the same conditions, and (c) and (d) compare the experimental roll angle ϕ, landing platform angle γ̂, and RLG absorbed ground angle γRLG to the simulated performance of the system. S-100 images are Copyright 2018 by Boeing: Approved for public release.
Fig. 14
(a) S-100 RLG response to a ±10 deg sinusoid with peak roll rates of ±3 deg/s where the aircraft starts the experiment on the table, (b) still images of a simulation under the same conditions, and (c) and (d) compare the experimental roll angle ϕ, landing platform angle γ̂, and RLG absorbed ground angle γRLG to the simulated performance of the system. S-100 images are Copyright 2018 by Boeing: Approved for public release.
Close modal

The plot in Fig. 14(c) shows fuselage, RLG, and estimated table angle data from this experiment. Roll controlled RLG maintains the fuselage roll angle (green data points) within or at the ϕDB (green shaded region) boundary for the entirety of the table's dynamic motion (red dashed line). Figure 14(d) shows the simulated response of the S-100 subjected to the similar landing conditions. The simulated sinusoid was computed using the method described in Sec. 3.2.1. All sinusoid simulations for the remainder of this work were precomputed using the peak rate and amplitude of the associated experimental test. These example simulations may not directly match the dynamics of the table, since there is inertia and lag in the table's motion and the virtual table moves without lag or inertia. The motion table had a peak roll rate of 3.5 deg/s, while the simulation was ideal at 3.0 deg/s. The authors computed the peak angular rate of the RLG using measurements from all angular encoders and computing a filtered, moving average derivate of the measurements. During the experiment, the RLG achieved angular rates between 4.5 and 4.8 deg/s. Figure 14(c) shows the table reached peaks of ±11 deg based on the onboard computed estimate of ground angle, γ̂, rather than the commanded ±10 deg. This was accounted for in the matching simulation. Regardless of the motion platform's performance, the roll feedback controller operated nominally for the entire sinusoid experiment. Furthermore, the simulation tool accurately represented the response of the experimental system.

The authors conducted additional experiments where the aircraft was lowered onto the platform to close the force and roll feedback loops. Figure 15(a) shows data from all RLG sensors during one of these experiments where the table was commanded a 5 deg/s peak roll rate sinusoid of ±10 deg with the aircraft lowered onto the table at 0.18 m/s. The table only achieves ±8–9 deg but it had peak angular rates of 10 deg/s for reference, which were both accounted for in the comparable simulation shown in Fig. 15. The table angle and rate appreciably lags when compared to the simulation at t = 1–2 s as the table adjusts its control gains to keep up with the added weight of the aircraft. Simulation results with the S-100 landed on a deck undergoing similar motion to the experiment are shown in Fig. 15(b). The differences between the experimental data and simulation data are due to the variable motion rates of the table that the simulation does not account for in this work. Visible in Figs. 15(a) and 15(b) are the bumps in ϕ and γRLG when the landing gear first make contact and begin to conform on the surface. While the aircraft does not undergo landing in the simulation at the same time as the experiment, the motion dynamics from 3 s onward match within 0.5 deg.

Fig. 15
(a) Experimental system roll and ground angle estimates where the S-100 and RLG are lowered at 0.18 m/s onto a platform undergoing a ±10 deg sinusoid with peak angular rate of 10 deg/s and (b) roll and ground angle output from a simulation with similar table dynamics
Fig. 15
(a) Experimental system roll and ground angle estimates where the S-100 and RLG are lowered at 0.18 m/s onto a platform undergoing a ±10 deg sinusoid with peak angular rate of 10 deg/s and (b) roll and ground angle output from a simulation with similar table dynamics
Close modal

These examples of roll sinusoid responsiveness of the S-100 and integrated RLG show the ability of the prototype system to operate with force feedback and roll feedback closed-loop control of the gear from initial contact through the full landing sequence. The RLG prototype conformed on the table with platform roll rates up to 10 deg/s on slopes up to 11 deg. Finally, the simulations most comparable to the experiments show comparable response of the aircraft and the gear.

4.1.1 Sea State Dynamic Platform Experiments.

The final set of experiments conducted on the S-100 were SS 3 and 5 platform landings. Experiments began in a similar fashion to the sinusoidal experiments, with the aircraft on the platform when the experiment began. The authors choose to show two of these experiments in this work, one SS 3 and one SS 5. Both of these experiments used the maximum allowable surge–sway–heave for the table, up to 10% of the maximum value output from STH and SMP as described in Sec. 3.2. Figures 16(a) and 16(b) depict the angle measurements, and (c) and (d) present the angular rates of the gear for these two experiments. The SS 3 data in (a) and (c) show mild waves where RLG is not required for safe aircraft operation due to the low roll angles associated with the platform. This is in stark contrast to the results of the SS 5 experiment. The particular SS 5 test chosen begins with mild waves, and within 10 s the waves begin to increase such that the ship deck peaks near ±10deg and angular rates at and above 5 deg. The operational range of the gear is apparent from these two vastly different experiments. For both systems, the gear maintains the experimental platform within 0.5 deg of the ϕDB bounds. Comparable simulations were executed to further validate the simulation design. Figures 16(e) and 16(f) illustrate the roll angle dynamics of the simulated systems undergoing the same motion table dynamics as the experimental system. The simulations closely match the response of the experimental system, with one minor difference. The SS 3 experimental system overcorrected the roll angle of the fuselage, resulting in a minor difference in roll angle response between the real and simulated systems (see Fig. 16(c) for the spike in γ˙RLG). The set of experiments conducted with the aircraft initialized on the motion table showed enough success to conduct experiments with the S-100 descending onto the platform while in motion.

Fig. 16
Angle measurements and post-processed angular rate estimates of the S-100 roll with the aircraft initialized on the Sarnicola motion table undergoing a sea state 3 (a) and (c) and sea state 5 (b) and (d) experiments, and comparable simulations for sea state 3 and 5 are depicted in (e) and (f), respectively. S-100 images are Copyright 2018 by Boeing: Approved for public release.
Fig. 16
Angle measurements and post-processed angular rate estimates of the S-100 roll with the aircraft initialized on the Sarnicola motion table undergoing a sea state 3 (a) and (c) and sea state 5 (b) and (d) experiments, and comparable simulations for sea state 3 and 5 are depicted in (e) and (f), respectively. S-100 images are Copyright 2018 by Boeing: Approved for public release.
Close modal

Figures 17(a)17(f) show two sea state 5 experimental landings and (g) and (h) depict simulation results with similar conditions as the experiments. Both of the landings descended onto the platform at 0.18 m/s, and operators initiated descent at random times during the platform's dynamic motion. Therefore, the dynamics seen in the experiments look significantly different. Additionally, the dynamic motion programs are limited to 5% of the maximum surge, sway, and heave outputs from SMP and STH simulations for operator and aircraft safety.

Fig. 17
(a) and (b) Angle measurements, (c) and (d) post-processed angular rate estimates, (e) and (f) foot-force measurements, and (g) and (h) simulated aircraft and gear response of the S-100 RLG with the aircraft initialized above the Sarnicola motion table and lowered onto the table, while it undergoes two sea state 5 dynamic motion programs(a) and (b) Angle measurements, (c) and (d) post-processed angular rate estimates, (e) and (f) foot-force measurements, and (g) and (h) simulated aircraft and gear response of the S-100 RLG with the aircraft initialized above the Sarnicola motion table and lowered onto the table, while it undergoes two sea state 5 dynamic motion programs
Fig. 17
(a) and (b) Angle measurements, (c) and (d) post-processed angular rate estimates, (e) and (f) foot-force measurements, and (g) and (h) simulated aircraft and gear response of the S-100 RLG with the aircraft initialized above the Sarnicola motion table and lowered onto the table, while it undergoes two sea state 5 dynamic motion programs(a) and (b) Angle measurements, (c) and (d) post-processed angular rate estimates, (e) and (f) foot-force measurements, and (g) and (h) simulated aircraft and gear response of the S-100 RLG with the aircraft initialized above the Sarnicola motion table and lowered onto the table, while it undergoes two sea state 5 dynamic motion programs
Close modal

Data in the plots of Figs. 17(a), 17(c), and 17(e) present the landing response when the gear make first contact, while the table increases roll angle up to a peak near 10 deg. First contact is made, and the combine impact force of the aircraft and table's heave pushed the force sensor measurement beyond the minimum threshold of force feedback control. Roll and force fused control maintains the aircraft at the boundary of ϕDB until the RLG achieves full contact. RLG tracked angular rates of 9.5 deg/s to match the combined dynamics of the platform and aircraft. The gear makes full contact with both legs and matches the deck motion within 2 s of descent onto the platform. For the next 10 s, the RLG maintains the aircraft within 1.5 deg of ϕDB. The larges deviation in roll occurs at 8 s, at which time the cable friction under the aircraft's weight causes significant differential torque on the drivetrain's spool. The torque on the motor was too large to keep up with the table's motion, so the angular rate trended toward 0 deg/s at this moment. Simulation of this landing showed similar behavior during contact with the table, where the fuselage moves beyond ϕDB before rapid RLG motion starts (see Fig. 17(g)). For the majority of the test, the simulation and experiment match together except for the peak in roll on the experimental system at 8 s. While the simulation includes a motor thermal model for degradation of control, it assumes the motor is in free space with air at 25 °C. Since the experimental drivetrain is enclosed in the fuselage and operated for long periods of time before this experiment, these assumptions cannot mimic the torque (heating) limitation of motor speed at 8 s. The thermal model is accurate for short bursts of operation less than 4 s, as verified from the multiple experiments in this section.

The data in the plots of Figs. 17(b), 17(d), 17(f), and 17(h) present a mild portion of SS 5 conditions. The largest angle magnitude achieved by the motion table during this landing was −6 deg, and the RLG never moved at angular rates above 6 deg/s to track the deck's motion. A simulation chosen for comparison has the same contact timing, but different roll response of the simulated platform. Even with this minor difference, the simulation and experimental system have matching responses to contact and platform motion in the first 4–5 s. The gear never exceeded 0.5 deg outside of ϕDB for either the simulation or experiments. These results further show the simulations as a useful characterization tool for the S-100 RLG system.

In both landing scenarios, the contact force varies depending on contact angle with the deck. The amount of force variation is an unexpected find of the experimental research, and it warrants study of force feedback inclusion after full contact in the future. Both landing scenarios presented in this work support that the experimental prototype has the sensor measurements and drivetrain response to land on a deck in sea state 5 conditions. These experimental results show at a limited level that the experimental prototype S-100 RLG meets the design performance expectations. Experiments and simulations of the same or similar landing conditions match well, so the authors explored additional landing envelope and dynamic platform conditions through simulation.

4.2 Extended Simulation Results.

Additional simulations were conducted on static and dynamic landing platforms for three RLG operational cases: no feedback control (legs locked), force feedback controlled, and roll and force feedback control. No feedback control (i.e., rigid landing gear) was simulated as a part of this section to show a direct comparison between the locked and active landing gear configuration responses with respect to dynamic rollover. This does not represent the dynamic rollover scenario for the nominal static gear on the aircraft. Force feedback control scenarios were simulated to compare to a controller used for previous flight tests [30]. The roll and force feedback controller is not expected to enhance performance during the short time horizon from initial contact to contact on all legs because of the fundamental acceleration (power) limits of the drivetrain as noted in Sec. 2.1. Nevertheless, this roll and force fused controller is expected to prevent dynamic rollover by pushing the rotorcraft's roll angle toward the ideal deadband of ϕDB=±0.1 deg once full contact is achieved. This deadband was chosen to simulate high quality AHRS that are integrated into rotorcraft, not the VN-100 used on the experimental platform. The last consideration for these simulations is the landing velocity. Each simulation presented here considers the case where the aircraft starts at 0 m/s descent, but collective is dropped at first contact, similar to how remote pilots fly the aircraft [30]. Typical aircraft descent rates once both legs have made contact range from 0.55 to 0.6 m/s for this scenario, at the upper end of the RLG's performance specification. Any additional vertical velocity from the platform will push the RLG up to its original design boundary and beyond. However, this method is similar to how pilots land the aircraft, and it accurately represents the flight-tested actuation. The landing simulations presented here represent the approximate performance boundary for the S-100 RLG, not a full range of performance estimates. Furthermore, they are meant to show general trends of landing performance outside the range of the experimental validation in Sec. 4 because of the numerous nonlinearities that may stem from faster, more aggressive landing situations.

4.2.1 Static Landings.

The first landing cases simulated were on static surfaces with the aircraft descending with a roll angle of −2 deg, the nominal approach roll angle of the aircraft [30]. Figure 18(a) illustrates the roll angle results for a 15 deg even slope landing without any feedback control (green line), force feedback controlled RLG (blue line), and force and roll feedback controlled RLG (black dashed line). Four shaded regions are shown in (a) to distinguish risk of dynamic rollover for the given roll angle. These regions were analytically observed as low risk (orange, 5–10% chance) and high risk (red, > 40% chance) for dynamic rollover in Sec. 3.1. As previously noted, the authors have no knowledge of the aircraft's dynamic rollover risk when the original equipment manufacturer rigid landing gear is installed, so this data is a baseline comparison of the active RLG system to the RLG without any active control (i.e., locked in the neutral position). The maximum descent speed achieved during this landing for the force and roll feedback controlled RLG is 0.56 m/s, so it is expected for the locked RLG to diverge and experience dynamic rollover during the short landing period. Since this descent rate is within the performance specification of the gear's drivetrain, the system compensates and maintains a roll angle throughout the landing with no risk of dynamic rollover. Figures 18(b)18(d) give a physical picture of the simulations performed and shown in Fig. 18(a). These images and the simulation outputs clearly show that any feedback control of RLG is beneficial to prevent dynamic rollover on extreme slopes. The benefit of including roll feedback is clearer for landings where the aircraft's inertial cross-coupling of the rotor impacts roll rate.

Fig. 18
Static surface landing simulation of a 15 deg slope without active RLG, with force feedback control, and roll and force feedback control with shaded regions for low risk (between –5 deg and 5 deg) and high risk (lower than –10 deg and higher than 10 deg) for dynamic rollover
Fig. 18
Static surface landing simulation of a 15 deg slope without active RLG, with force feedback control, and roll and force feedback control with shaded regions for low risk (between –5 deg and 5 deg) and high risk (lower than –10 deg and higher than 10 deg) for dynamic rollover
Close modal

Figure 19(a) presents the roll angle data for a −10 deg roll angle sloped surface. In this scenario, the pitch–roll rate inertial cross-coupling exacerbates the roll rate increase during landing which the roll feedback controller is able to mitigate. The simulation with force feedback control leaves the aircraft above 5 deg of roll. This is consistent with the findings of León et al. [30] during flight test experiments of the force feedback controller. The force and roll feedback controller is not able to completely mitigate roll angle increase prior to full contact of the legs. This is due to the physical power and speed constraints of the RLG. The added feedback channel (roll) allows the gear to return the aircraft to a safe operating region rapidly after full contact is made. A custom, higher power actuator would be able to mitigate this issue, though. Nevertheless, the roll and force feedback fused control is especially important for limiting the dynamic rollover risk during subsequent takeoff when collective is increased and thrust begins to push the aircraft to roll about contact points. With the fused control system, the aircraft is able to takeoff from a near level configuration unlike the other configurations.

Fig. 19
(a) Roll angle data from a static surface landing simulation of a −10 deg slope without active RLG, with force feedback control, and roll and force feedback control with shaded regions for low risk (between −5 deg and 5 deg) and high risk (lower than −10 deg and higher than 10 deg) for dynamic rollover and (b) images of each simulation to illustrate the data at t = 1.25 s
Fig. 19
(a) Roll angle data from a static surface landing simulation of a −10 deg slope without active RLG, with force feedback control, and roll and force feedback control with shaded regions for low risk (between −5 deg and 5 deg) and high risk (lower than −10 deg and higher than 10 deg) for dynamic rollover and (b) images of each simulation to illustrate the data at t = 1.25 s
Close modal

4.2.2 Roll Sinusoid Platform Landings.

Simulations were conducted landing the S-100 on a moving platform with only sinusoidal roll motion. The sine waves used to generate the deck motion have motion from ±15 deg in the frequency ranges that correspond to waves with maximum angular rates of 3–8 deg/s, respectively. Sinusoid waves of these amplitudes and frequencies are within the sea state 6 simulated ship deck response as shown in Fig. 4(b). The primary benefit of simulations using only roll sinusoidal is the lack of surge, sway, and heave motion which adds to the dynamic interaction. These simulations are designed to be compared to sea state platform simulations to determine the impact of decoupled angular motion versus translation motion on the RLG's overall performance.

Figure 20(a) shows the roll angle data for the three operation modes landing on a platform undergoing sinusoidal motion. The platform's sinusoid had a peak roll rate of ±5 deg/s and was initialized at −10 deg. Figure 20(b) illustrates each operation mode's orientation at t = 1.0 s, when the locked RLG is recovering from a bounce. The aircraft for each simulation has a peak descent speed of 0.61 m/s, near the bound of design specification for the S-100 RLG and well above the experimental tests conducted in Sec. 4.1. Simulation with locked RLG does not roll over, but it has a significant and dangerous bounce. The authors note that the aircraft does not enter dynamic rollover because the platform was moving toward γ=0deg and the risk to enter dynamic is low at that roll angle. Both RLG controllers prevent this large roll angle bounce and keep the aircraft outside the low risk region. Again the force and roll feedback controlled system quickly (<0.5 s) returns the aircraft to a level position even as the deck continues to move. The force feedback controller does not have any sensor feedback to track the deck motion, so the aircraft is a ride along on the platform with a constant offset. This puts the aircraft in an unusual configuration for recovery or take off at a later time. All of these issues are dealt with by the force and roll feedback controlled gear, since it converges toward and maintains the allowable ϕDB.

Fig. 20
(a) Roll angle simulation data landing on a moving platform sinusoid initial condition of −10 deg for locked, force feedback controlled, and force and roll feedback controlled RLG, and (b) renders of the system for each simulation at t = 1.0 s
Fig. 20
(a) Roll angle simulation data landing on a moving platform sinusoid initial condition of −10 deg for locked, force feedback controlled, and force and roll feedback controlled RLG, and (b) renders of the system for each simulation at t = 1.0 s
Close modal

Similar behavior is observed from a landing where the platform's initial condition is 10 deg with positive angular rate. Figure 21(a) shows the roll feedback of the simulations for each operating mode. The authors note that the aircraft has a significant bounce, but does not enter dynamic rollover in these ideal conditions (no wind modeled). In the locked RLG case, there is a 1 deg bias between the aircraft roll and the deck roll angle, which is due to semipermanent stretch of the cable on one side of the RLG assembly. This semipermanent asymmetry is caused by the large impact force on the right leg.

Fig. 21
Moving platform landing (a) roll angle outputs of a deck sinusoid with initial condition of 10 deg for locked, force feedback controlled, and force and roll feedback controlled RLG, and (b) γ˙RLG for the force and roll feedback controlled system
Fig. 21
Moving platform landing (a) roll angle outputs of a deck sinusoid with initial condition of 10 deg for locked, force feedback controlled, and force and roll feedback controlled RLG, and (b) γ˙RLG for the force and roll feedback controlled system
Close modal

Some oscillations are visible on the force and force and roll feedback RLG simulations, which Fig. 21(b) further elucidates for the force and roll feedback RLG simulation. The oscillations are caused by momentary loss of contact. The force feedback control portion of both controllers is not perfectly tuned across a wide range of landing speeds; therefore, the gear can momentarily lose contact. Any momentary loss of contact results in a deceleration of the gear until contact is regained. This prevents the gear from impacting the surface with significant force, further destabilizing roll angle. In this case, the gear loses contact four times during the landing sequence (see Fig. 21(b)). Regardless of this contact bouncing, the RLG is fast enough to bound the roll angle within a safe operating region during the time before full contact is established around 0.4 s. There is then a transient period where the torque on the spool from both cables is too high to achieve significant γ˙RLG. By 0.6 s, the torque stabilizes to a point where the drivetrain can drive the legs at a rate of 7.7 deg/s until the roll is within ϕDB. These sinusoid platform simulations show the force and roll feedback RLG to bound the roll angle to a safe operating region. The force and roll-fused control maintains the roll angle near the desired deadband even with significant pitch–roll inertial cross-coupling caused by the rotor blade and landing gear geometry.

4.2.3 Sea State Platform Landings.

The final set of simulations conducted were on dynamic platforms simulated SS 5 and 6 conditions. These two conditions were chosen since they represent a combined 37% probability of occurrence in the northern hemisphere [50], and the dynamics of the ship decks at these SS pose a risk to crew and aircraft.

Figure 22 depicts the fuselage roll angle output of a 4 s simulation where the aircraft initiates contact on the landing surface while it is near 10 deg. This particular case leads to dynamic rollover of the locked RLG system. There is rolling motion of the fuselage for the active controlled systems as well, but no dynamic rollover or dangerous bouncing behavior. The significant rolling of the fuselage of the actively controlled simulations is because the combined impact velocity of the helicopter decks is 1.3 m/s because of the deck's heave motion. This is almost two times the maximum speed of the active RLG, so they are unable to absorb all of the roll angle. The force and roll feedback RLG returns the fuselage to a safe operating condition within 1 s of initial contact even though the combine system dynamics are well outside of the nominal operating specification.

Fig. 22
Roll angle outputs from a S-100 landing on a sea state 5 simulated dynamic platform
Fig. 22
Roll angle outputs from a S-100 landing on a sea state 5 simulated dynamic platform
Close modal

More detail of the drivetrain's response to this aggressive landing is illustrated in Fig. 23. It is clear the system is accelerating as quickly as possible to keep up with the combined deck and aircraft dynamics until contact is made on all legs and the torque is too high to continue accelerating at the same magnitude. Once the drivetrain slows and the torque stabilizes across the spool, the RLG system begins to drive toward ϕDB at t = 0.7 s. Noise in the data between the full contact and return toward level is from response dynamics of the cable, ground, and elastic joints during the impact.

Fig. 23
RLG ground angle absorption rate, γ˙RLG, outputs from a S-100 landing over the 1 s landing window on a sea state 5 simulated dynamic platform using the fused roll and force feedback controller
Fig. 23
RLG ground angle absorption rate, γ˙RLG, outputs from a S-100 landing over the 1 s landing window on a sea state 5 simulated dynamic platform using the fused roll and force feedback controller
Close modal

Similar responses result from sea state 6 simulated landings. Figure 24(a) depicts the roll response to a sea state 6 landing when contact begins with the deck near −7 deg and 1.28 m/s relative vertical velocity, while Fig. 24(b) shows the response to a sea state 6 landing with contact initiated at 10 deg and a relative velocity between the deck and aircraft being 0.98 m/s. In Fig. 24(a), the difference between the roll angles of the three systems is small because of the large relative velocity and inertial cross-coupling effect on roll rate. This represents an aggressive landing outside the nominal operation bounds of the experimental RLG. The active RLG systems have a minor bounce but do not enter dynamic rollover. The gear's ability to push against the slope prevents the aircraft from a dangerous bounce and rollover. The locked RLG does enter dynamic rollover in the low risk region because of the high landing velocity relative to the nominal value (see Table 3). This ability to combat dynamic rollover is more pronounced for the force and roll feedback system, which pushes the aircraft into the safe operating region well before the deck reaches a safe roll angle.

Fig. 24
RLG ground angle absorption rate, γ˙RLG, outputs from a S-100 landing over the 1 s landing window on a sea state 5 simulated dynamic platform using the fused roll and force feedback controller
Fig. 24
RLG ground angle absorption rate, γ˙RLG, outputs from a S-100 landing over the 1 s landing window on a sea state 5 simulated dynamic platform using the fused roll and force feedback controller
Close modal

Each simulation presented shows that the force feedback controller is, at minimum, sufficient for reducing roll instability, but the danger with such a controller is clear in Fig. 24(b). The gear using force feedback absorbs significant ground angle (8 deg) which in turn puts the aircraft at a significant inertial roll angle of 17 deg by the end of the simulation at t = 5 s. At this roll angle, the aircraft is in significant danger of rollover if any collective is applied to the rotor. This is an example of a case where force feedback RLG is worse than a rigid gear configuration, showing the need for continued feedback after landing. The force and roll feedback controlled RLG eliminates this issue by maintaining the aircraft roll near level.

4.2.4 Summary of Extended Simulation Results.

The set of simulations presented here are representative of the conditions; the S-100 RLG would be subjected to in the real world, from static to dynamic simulated ship decks. A rigid landing gear configuration enters dynamic rollover repeatedly even with the ideal conditions for collective (thrust) drop and no simulated wind. Both active RLG configurations prevent roll instability under the same static and mobile platform landing conditions, but the system using only force feedback has a significant probability of the aircraft entering a dangerous condition when landing on a dynamic platform (see Fig. 24(b)). Even with the most aggressive sea state 5 and 6 conditions with platform-aircraft relative impact velocities at 2× the maximum rating of the experimental landing gear, the aircraft with force and roll feedback did not enter dynamic rollover. These simulations show that the force and roll feedback controller is theoretically effective to eliminate roll dynamic instability and ensures safe operating conditions for the S-100 in conditions up to the extremes of sea state 6.

5 Conclusion

This work presented a cable-driven, four-bar link RLG integrated onto the S-100 Camcopter as a method to aid landings on dynamic platforms—namely, maritime ship decks. The experimental S-100 RLG was proven to have ample dynamic response landing on moving platforms using a discrete, real-time roll and foot-force feedback control. Experiments on dynamic motion platform simulating sinusoidal and SS 5 conditions showed the feedback controller and drivetrain are able to maintain a user-defined deadband and maintain safe operating orientation for the S-100. Each rotor-off experiment verified simulation results at a basic level by testing the hardware on a real motion table. The simulation tools complement the experimental findings and further demonstrate that the drivetrain and feedback control chosen for the system are capable of preventing dynamic rollover in aggressive sea state 6 conditions. This work is the first necessary step toward autonomous landings with RLG on static or moving platforms, and commercial applications of this technology are far reaching for rotorcraft and the expanding unmanned vertical takeoff and landing aircraft space. There are a number of future work tasks that could improve the usability of such a four-bar RLG system. A few of these future work tasks are itemized below:

  1. Incorporate aerodynamic and ship wake flow fields into the multibody dynamics simulation to better understand how the wind and gusts impact landing stability after the rotorcraft collective has been dropped. This will aid in the statistical understanding of how RLG performs as compared to locked (rigid) landing gear.

  2. RLG sensing and control systems may be implemented with onboard flight controls and landing approach flight software in order to fully realize the commercial potential of the technology.

  3. The RLG control may be integrated with vision or LIDAR-based landing guidance systems. This would give the landing gear a preemptive image of the landing surface to preconform or prime the actuators prior to landing.

  4. Study the scalability of the four-bar RLG structure for full-size helicopters. The scalability for a cable-driven design is not immediately clear. However, the four-bar link mechanism has distinct advantages that will scale to full-size helicopters.

Acknowledgment

The authors would like to thank DARPA, and in particular, Ashish Bagai, Timothy Chung, Ryan Hofmeister, Natalie Van Osch, Matthew Tallent, and Michael Ol for providing valuable comments and guidance throughout the S-100 program. The authors would like to thank the Boeing-Mesa team, and in particular, Ray Lopez, James Roos, Steve Low, John Waxler, Roy Bergman, James Bonds, and Bill Brady for their assistance during ground and flight testing of the S-100 program. Finally, the authors thank Jacob Wachlin and Kyle Solomon for their guidance and assistance throughout the program. The views, opinions, and/or findings expressed are those of the author(s) and should not be interpreted as representing the official views or policies of the Department of Defense or the U.S. Government.

Funding Data

  • Defense Advanced Research Projects Agency (DARPA) Tactical Technology Office (TTO) (Contract No. FA8650-12-C-7276; Funder ID: 10.13039/100000185).

Nomenclature

B =

distance between the contact points of the left and right legs

c =

standard linear solid contact model parallel damper constant

CF =

linear scaling coefficient between pressure on each foot and the force applied on that foot

F =

standard linear solid contact model force

F¯ =

vector containing six degrees-of-freedom dynamic equations of motion

F̂d =

desired foot contact force during landing

F̂L/R =

force estimate on the left or right foot

[G] =

linear transform to map constraint forces and moments to multibody equations of motion

k =

standard linear solid contact model parallel spring constant

kv =

standard linear solid contact model series spring constant

KD =

derivative control gain

KP =

proportional control gain

L1 =

crank length

L2 =

distance between the crank and follower

P̃ =

normalized pressure measured on a RLG foot

rspool =

RLG drivetrain spool radius

sL/R =

cable length of the left or right half of the RLG assembly

T =

rotor thrust

Tc =

commanded rotor thrust

u =

aircraft body relative forward velocity

U¯ =

vector containing multibody constraint forces and moments

w =

aircraft body relative downward velocity

X¯ =

state vector

α =

low-pass filter time constant

β˙ =

angular speed of the RLG drivetrain spool

γ =

ground angle

γ̂ =

onboard estimate of the ground angle based on roll angle and RLG leg angles

γ˙RLG =

angular rate of the active landing gear

δ =

standard linear solid contact model displacement

Δh =

vertical height between the left and right leg contact points

ζ =

damping ratio

ηlat =

thrust angle relative to the fuselage's lateral plane

ηlong =

thrust angle relative to the fuselage's longitudinal plane

ηlatc =

commanded thrust angle relative to the fuselage's lateral plane

ηlongc =

commanded thrust angle relative to the fuselage's longitudinal plane

θ =

crank orientation angle relative to the fuselage body relative vertical axis

μ =

friction coefficient

τT =

first-order model time constant for rotor thrust

τη =

first-order model time constant for the angle of rotor thrust

ϕ =

roll angle

ϕ̂ =

roll angle estimate

ϕd =

desired fuselage roll angle in feedback controller

ϕDB =

deadband roll angle for roll feedback control

ωn =

natural frequency

References

1.
FAA,
2012
, “
Advanced Landing Maneuvers
,”
Helicopter Flying Handbook
,
Federal Aviation Administration
, Washington, DC, pp.
10.7
10.9
.
2.
Fang
,
R.
, and
Finlay
,
B.
,
1992
, “
Determination of Limitations for Helicopter Ship-Borne Operations
,”
Aircraft Ship Operations
,
Advisory Group for Aerospace Research and Development (AGARD)
,
Amsterdam, The Netherlands
, pp.
18.1
19.14
.
3.
Ibrahim
,
R. A.
, and
Grace
,
I. M.
,
2010
, “
Modeling of Ship Roll Dynamics and Its Coupling With Heave and Pitch
,”
Math. Probl. Eng.
,
2010
, pp.
1
32
.10.1155/2010/934714
4.
Moon
,
J.
,
Domercant
,
J.
, and
Mavris
,
D.
,
2015
, “
A Simplified Approach to Assessment of Mission Success for Helicopter Landing on a Ship
,”
Int. J. Control, Autom. Syst.
,
13
(
3
), pp.
680
688
.10.1007/s12555-013-0092-y
5.
Forrest
,
J.
,
Owen
,
I.
,
Padfield
,
G.
, and
Hodge
,
S.
,
2012
, “
Ship-Helicopter Operating Limits Prediction Using Piloted Flight Simulation and Time-Accurate Airwakes
,”
J. Aircr.
,
49
(
4
), pp.
1020
1031
.10.2514/1.C031525
6.
2001
, “
Helicopter Operations in the Maritime Environment
,”
Aircr. Eng. Aerosp. Technol.
,
73
(
4
), pp.
2037
2049
.10.1108/aeat.2001.12773dac.001
7.
Varol
,
A. E.
, and
Gunal
,
M. M.
,
2015
, “
Simulating Prevention Operations at Sea Against Maritime Piracy
,”
J. Oper. Res. Soc.
,
66
(
12
), pp.
2037
2049
.10.1057/jors.2015.34
8.
2010
, “
Maritime Force Cornerstone of Relief Operations in Haiti, Navy Says
,”
Sea Technol.
,
51
(
3
), p.
53
.
9.
Schiebel,
2019
, “
Camcopter S-100 UAV
,” Schiebel, Vienna, Austria, accessed Apr. 19, 2019, https://www.airforce-technology.com/projects/camcopters-100uav/
10.
Barge Master,
2014
, “
Launching of the Motion Compensated Helideck
,” Barge Master, Rotterdam, The Netherlands, accessed Apr. 19, 2019, https://www.barge-master.com/item/launching-of-the-motion-compensated-helideck/
11.
Reber
,
E. A.
, and
Bernard
,
B. J.
,
2012
, “
The Sea of Simulation: Improving Naval Shiphandling Training and Readiness Through Game-Based Learning
,” Naval Postgraduate School, Monterey, CA, Technical Report No.
ADA561764
.https://apps.dtic.mil/sti/pdfs/ADA561764.pdf
12.
Xu
,
G.
,
Zhang
,
Y.
,
Ji
,
S.
,
Cheng
,
Y.
, and
Tian
,
Y.
,
2009
, “
Research on Computer Vision-Based for UAV Autonomous Landing on a Ship
,”
Pattern Recognit. Lett.
,
30
(
6
), pp.
600
605
.10.1016/j.patrec.2008.12.011
13.
Saripalli
,
S.
,
Montgomery
,
J.
, and
Sukhatme
,
G.
,
2003
, “
Visually Guided Landing of an Unmanned Aerial Vehicle
,”
IEEE Trans. Rob. Autom.
,
19
(
3
), pp.
371
380
.10.1109/TRA.2003.810239
14.
Polvara
,
R.
,
Sharma
,
S.
,
Wan
,
J.
,
Manning
,
A.
, and
Sutton
,
R.
,
2019
, “
Autonomous Vehicular Landings on the Deck of an Unmanned Surface Vehicle Using Deep Reinforcement Learning
,”
Robotica
,
37
(
11
), pp.
1867
1882
.10.1017/S0263574719000316
15.
Voskuijl
,
M.
,
Padfield
,
G.
,
Walker
,
D.
,
Manimala
,
B.
, and
Gubbels
,
A.
,
2010
, “
Simulation of Automatic Helicopter Deck Landings Using Nature Inspired Flight Control
,”
Aeronaut. J.
,
114
(
1151
), pp.
25
34
.10.1017/S000192400000350X
16.
Sandino
,
L. A.
,
Bejar
,
M.
, and
Ollero
,
A.
,
2011
, “
On the Applicability of Linear Control Techniques for Autonomous Landing of Helicopters on the Deck of a Ship
,”
2011 IEEE International Conference on Mechatronics
, IEEE, Istanbul, Turkey, Apr. 13–15,
pp.
363
368
.10.1109/ICMECH.2011.5971312
17.
Horn
,
J. F.
,
He
,
C.
,
Roark
,
S.
,
Yang
,
J.
,
Tritschler
,
J.
,
Gonzalez
,
G.
, and
Lee
,
D.
,
2016
, “
Autonomous Control Modes and Optimized Path Guidance for Shipboard Landing in High Sea States
,” Office of Naval Research, University Park, PA, Technical Report No. AD1013549.
18.
Tan
,
C. K.
,
Wang
,
J.
,
Paw
,
Y. C.
, and
Liao
,
F.
,
2016
, “
Autonomous Ship Deck Landing of a Quadrotor Using Invariant Ellipsoid Method
,”
IEEE Trans. Aerosp. Electron. Syst.
,
52
(
2
), pp.
891
903
.10.1109/TAES.2015.140850
19.
Raibert
,
M. H.
,
1986
, “
Legged Robots
,”
Commun. ACM
,
29
(
6
), pp.
499
514
.10.1145/5948.5950
20.
Buehler
,
M.
,
Playter
,
R.
, and
Raibert
,
M.
,
2005
, “
Robots Step Outside
,”
International Symposium on Adaptive Motion of Animals and Machines
(
AMAM
), Ilmenau, Germany, pp.
1
4
.https://www.researchgate.net/publication/280776075_Robots_Step_Outside
21.
Ahmadi
,
M.
,
Michalska
,
H.
, and
Buehler
,
M.
,
2007
, “
Control and Stability Analysis of Limit Cycles in a Hopping Robot
,”
IEEE Trans. Rob.
,
23
(
3
), pp.
553
563
.10.1109/TRO.2007.898956
22.
Raibert
,
M.
,
Blankespoor
,
K.
,
Nelson
,
G.
, and
Playter
,
R.
,
2008
, “
Bigdog, the Rough-Terrain Quadruped Robot
,”
IFAC Proc. Vol.
,
41
(
2
), pp.
10822
10825
.10.3182/20080706-5-KR-1001.01833
23.
Manivannan
,
V.
,
Langley
,
J. P.
,
Costello
,
M.
, and
Ruzzene
,
M.
,
2013
, “
Rotorcraft Slope Landings With Articulated Landing Gear
,”
AIAA
Paper No. 2013-5160. 10.2514/6.2013-5160
24.
Sikorsky
,
I. I.
,
1953
, “
Landing Gear for Use on Inclined Surfaces
,” U.S. Patent No. US2630989A.
25.
Mason
,
S.
,
1974
, “
Helicopter Self-Leveling Landing Gear
,” U.S. Patent No. US3857533A.
26.
Kiefer
,
J.
,
Ward
,
M.
, and
Costello
,
M.
,
2016
, “
Rotorcraft Hard Landing Mitigation Using Robotic Landing Gear
,”
ASME J. Dyn. Syst., Meas., Control
,
138
(
3
), p.
031003
.10.1115/1.4032286
27.
DARPA
,
2015
, “
Robotic Landing Gear Could Enable Future Helicopters to Take Off and Land Almost Anywhere
,” Defense Advanced Research Projects Agency (DARPA), Arlington, VA, accessed Apr. 19, 2019, https://www.darpa.mil/news-events/2015-09-10
28.
Di Leo
,
C. V.
,
Leon
,
B.
,
Wachlin
,
J.
,
Kurien
,
M.
,
Rimoli
,
J. J.
, and
Costello
,
M.
,
2018
, “
Cable-Driven Four-Bar Link Robotic Landing Gear Mechanism: Rapid Design and Survivability Testing
,”
AIAA
Paper No. 2018-0491.10.1115/2018-0491
29.
Di Leo
,
C. V.
,
León
,
B.
,
Wachlin
,
J.
,
Kurien
,
M.
,
Krishnan
,
A.
,
Krishnan
,
A.
,
Rimoli
,
J. J.
, and
Costello
,
M.
,
2020
, “
Design of a Crashworthy Cable-Driven Four-Bar Link Robotic Landing Gear System
,”
J. Aircr.
,
57
(
2
), pp.
224
244
.10.2514/1.C035386
30.
León
,
B.
,
Rimoli
,
J.
, and
Leo
,
C. V. D.
,
2019
, “
Ground and Flight Tests of a Cable-Driven Four-Bar Linkage Robotic Landing Gear for Rotorcraft
,”
Vertical Flight Society 75th Annual Forum & Technology Display
, Philadelphia, PA, May 13–16, Paper No. 0294.
31.
Stolz
,
B.
,
Brödermann
,
T.
,
Castiello
,
E.
,
Engelberger
,
G.
,
Erne
,
D.
,
Gasser
,
J.
,
Hayoz
,
E.
,
Müller
,
S.
,
Mühlebach
,
L.
,
Löw
,
T.
,
Scheuer
,
D.
,
Vandeventer
,
L.
,
Bjelonic
,
M.
,
Günther
,
F.
,
Kolvenbach
,
H.
,
Höpflinger
,
M.
, and
Hutter
,
M.
,
2018
, “
An Adaptive Landing Gear for Extending the Operational Range of Helicopters
,”
2018 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS 2018)
, Madrid, Spain, Oct. 1–5, pp.
1757
1763
.10.1109/IROS.2018.8594062
32.
Huang
,
M.
,
Nie
,
H.
,
Zhang
,
M.
,
Wei
,
X.
, and
Yue
,
S.
,
2016
, “
Design of Mission Adaptive Landing Gear for Near Space Travel Lander
,”
J. Vibroeng.
,
18
(
8
), pp.
4949
4963
.10.21595/jve.2016.17078
33.
Kim
,
D.
, and
Costello
,
M.
,
2016
, “
Virtual Model Control of Rotorcraft With Articulated Landing Gear for Shipboard Landing
,”
AIAA
Paper No. 2016-1863.10.1115/2016-1863
34.
Stengel
,
R. F.
,
1994
,
Optimal Control and Estimation
(Dover Books on Advanced Mathematics),
Dover Publications
,
New York
.
35.
Barnes
,
D.
,
2015
, “
USS Curtis Wilbur (DDG 54) Crosses Into 2016
,” Navy Live, East China Sea, At Sea, accessed July 28, 2021, https://www.dvidshub.net/image/1026948/uss-curtis-wilbur-transits-east-china-sea
36.
León
,
B.
,
Rimoli
,
J.
, and
Di Leo
,
C. V.
,
2019
, “
Elastomer Encapsulated Pressure Sensor With Engineered Air Cavity for Force Sensing
,”
IEEE Sens. J.
,
19
(
16
), pp.
6628
6643
.10.1109/JSEN.2019.2912515
37.
Smith
,
J. O.
,
2008
, “
The Simplest Low-Pass Filter
,”
Introduction to Digital Filters: With Audio Applications
,
W3K Publishing
, Stanford, CA, pp.
1
25
.
38.
Leylek
,
E.
,
Ward
,
M.
, and
Costello
,
M.
,
2012
, “
Flight Dynamic Simulation for Multibody Aircraft Configurations
,”
J. Guid., Control, Dyn.
,
35
(
6
), pp.
1828
1842
.10.2514/1.55858
39.
Gross
,
M.
,
Costello
,
M.
, and
Fresconi
,
F.
,
2013
, “
Impact Point Model Predictive Control of a Spin-Stabilized Projectile With Instability Protection
,”
AIAA
Paper No. 2013-4509.10.1115/2013-4509
40.
Rossikhin
,
Y. A.
, and
Shitikova
,
M. V.
,
2011
, “
The Analysis of the Impact Response of a Thin Plate Via Fractional Derivative Standard Linear Solid Model
,”
J. Sound Vib.
,
330
(
9
), pp.
1985
2003
.10.1016/j.jsv.2010.11.010
41.
Carbone
,
G.
, and
Putignano
,
C.
,
2013
, “
A Novel Methodology to Predict Sliding and Rolling Friction of Viscoelastic Materials: Theory and Experiments
,”
J. Mech. Phys. Solids
,
61
(
8
), pp.
1822
1834
.10.1016/j.jmps.2013.03.005
42.
Goraj
,
Z. J.
, and
Cichocka
,
E.
,
2016
, “
Influence of Weak and Strong Gyroscopic Effects on Light Aircraft Dynamics
,”
Aircr. Eng. Aerosp. Technol.
,
88
(
5
), pp.
613
622
.10.1108/AEAT-03-2015-0076
43.
Kendoul
,
F.
,
Fantoni
,
I.
, and
Lozano
,
R.
,
2006
, “
Modeling and Control of a Small Autonomous Aircraft Having Two Tilting Rotors
,”
IEEE Trans. Rob.
,
22
(
6
), pp.
1297
1302
.10.1109/TRO.2006.882956
44.
Rothhaar
,
P. M.
,
Murphy
,
P. C.
,
Bacon
,
B. J.
,
Gregory
,
I. M.
,
Grauer
,
J. A.
,
Busan
,
R. C.
, and
Croom
,
M. A.
,
2014
, “
NASA Langley Distributed Propulsion VTOL Tiltwing Aircraft Testing, Modeling, Simulation, Control, and Flight Test Development
,”
AIAA
Paper No. 2014-2999. 10.2514/6.2014-2999
45.
Syms
,
G.
, and
Zan
,
S. J.
,
1994
, “
Analysis of Rotor Forces in a Ship Airwake
,” Institute for Aerospace Research, National Research Council of Canada, Ottawa, ON, Canada, Technical Report No.
ADA640406
.https://apps.dtic.mil/sti/citations/ADA640406
46.
Snyder
,
M.
,
Kang
,
H.
,
Brownell
,
C.
, and
Burks
,
J.
,
2013
, “
Validation of Ship Air Wake Simulations and Investigation of Ship Air Wake Impact on Rotary Wing Aircraft
,”
Nav. Eng. J.
,
125
(
1
), pp.
69
79
.https://www.usna.edu/MechEngDept/_files/documents/brownell_files/NEJ%202013.pdf
47.
Conrad
,
R. E.
,
2005
, “
SMP95: Standard Ship Motion Program User Manual
,” Ship Hydromechanics Department, Naval Surface Warfare Center: Carderock Division, West Bethesda, MD, Technical Report No. NSWCCD-50-TR-2005/074.
48.
Meyers
,
W. G.
,
Benne
,
C.
, and
Applebee
,
T.
,
1993
, “
Users Manual for the Simulation Time History and Access Time History Programs
,” Ship Hydromechanics Department, Naval Surface Warfare Center: Carderock Division, Bethesda, MD, Technical Report No. ADM000225.
49.
Kinsman
,
B.
,
1965
,
Wind Waves: Their Generation and Propagation on the Ocean Surface
,
Prentice-Hall
, Englewood Cliffs, NJ.
50.
Bales
,
S.
,
1982
, “
Designing Ships to the Natural Environment
,” Association of Scientists and Engineers, Washington, DC, Technical Report No. ADA114147.
51.
Sarnicola Simulation Systems, Inc.,
2007
, “
Hexad
,” Sarnicola Simulation Systems, Inc., Conklin, NY, accessed Apr. 19, 2019, https://sarnicola.com/hexad3000 h.html