I Introduction
The key hypothesis this work tries to inspect is: ”The feasibility of performing flip turns, a well known attribute of a batlike landing maneuver, through the manipulation of inertial dynamics while excluding aerodynamic forces.” Bats (and birds) possess no energyhungry motors widely used in flying robots for thrust vectoring yet they are more capable than any of these systems when agility and, of course, energy efficiency of flight are concerned.
Flying vertebrates apply the combination of inertial dynamics and aerodynamics manipulations to showcase extremely agile maneuvers. Unlike rotary and fixedwing systems wherein aerodynamic surface (i.e., ailerons, rudders, propellers, etc.) come with the sole purpose of aerodynamic force adjustment, the wings (also called appendages) in birds and bats possess more sophisticated roles. It is known that birds perform zeroangularmomentum turns by making differential adjustments (e.g., collapsing armwings) in the inertial forces led by one wing versus the other. Or, bats apply a similar mechanism to perform sharp banking turns [11].
Among these maneuvers, landing (or perching), which flying vertebrates do it in one way or another for a variety of reasons (e.g., transition to walking, resting on a perch, hanging from the ceiling of a cave, etc.), is an interesting maneuver to take inspiration from for aerial robot designs. Perching birds rotate their wings so that the aerodynamic drag is increased by creating a highpressure region inside of the wings and a lowpressure region behind the wings. This brings the wings to a stalled condition at which point the generated lift is equal to zero and the animal falls naturally while employing the legs as a landing gear.
Bats do it in a radically different way. After the selfcreated stalled condition, they manifest an acrobatic heelsabovehead maneuver that involves catapulting the lower body in a similar way that a free style swimmer flip turns.
Perching insects and birds have been the source of inspiration and biomimicry of them has led to interesting robot designs in recent years. Remarkably, the biomimicry of batlike landing is overlooked mainly because not only the aerodynamics adjustments are involved but also unique design provisions are required to allow for the manipulation of inertial dynamics.
A part from the ordeals associated with the design and control of a robot that can land similar to a biological bat a number of unique applications can be identified for these systems. For instance, in scenarios wherein maintaining a high vantage point for extended times is demanded (surveillance and reconnaissance) and limited power budget does not allow to hover for extended time periods hanging from elevated structures such as steel frames in buildings can allow these systems not only safely accomplish their missions but also harvest energy within the time period the system is the most vulnerable.
In addition to the exciting applications, a batstyle landing maneuver is extremely rich in dynamics and control and its characteristics are overlooked. Much of attention has been paid to simpler dynamics such as hovering and straight flight. While mathematical models of insectstyle, rotary and fixedwing robots of varying size and complexity are relatively well developed, models of airborne, fluidicbased vertebrates locomotion remain largely open due to the complex body articulation involved in their flight.
The mainstream school of thought inspired by insect flight has conceptualized wing as a massless, rigid structure, which is nearly planar and translates – as a whole or in twothree rigid parts – through space. In this view, wings possess no inertial effect, are fast that yield twotimescale dynamics, permit quasistatic external force descriptions, and tractable dynamical system. Unfortunately, these paradigms fail to provide insight into airborne, vertebrate locomotion and an ingredient of a more complete and biologically realistic model is missing, that is, the manipulation of inertial dynamics. The manipulation of inertial dynamics is an underappreciated aspect in existing paradigms.
The objective of this work is to manifest closedloop aerial body reorientation and preparation for landing through the manipulation of inertial dynamics using a tiny ballistic robot whose characteristics are carefully scaled to match that of a small bat. It is shown that despite a number of prohibitive restrictions, and lack of multithruster designs typically found in quadrotors, extremely fast body reorientations and preparations for landing through inertial dynamics manipulations not only is possible but also is an effective biologically meaningful solution.
A brief overview of the system will be presented followed by a hybridmodel description of the landing maneuver. A closedloop controller will be designed and the simulation and experimental results will be reported at the end followed by final remarks and conclusion.
Ii A hybridmodel description of upsidedown perching
A little reflection reveals, however, that a perching bat dynamics like other complex animal behavior, which emerges from complex interactions between neural and sensorymotor systems, can be described with tractable mathematical models. The legitimate question that is posed here is: ”Can a hybrid model capture the dynamics?” What encourages me to pose the question directly stems from inspecting highspeed images of the landing maneuver of bats and that in a time envelope not exceeding about one fifth of a second before landing completion (touchdown) the appendages are fully retracted strongly suggesting that the aerodynamic forces are less likely to be as important as other forces. Based on a similar observation, [3] conducted a simulationbased research and concluded that the inertial dynamics manipulation likely are the dominant players in a landing bat.
It is assumed for two modes: mode (1) with aerodynamics dominance (inertial forces are not negligible but are less effective) and mode (2) within which the inertial dynamics are dominant. In a perching bat, a transition from one mode to another triggered by propreoceptive sensing and neural reflexes occur at a remarkable high speed.
This hybrid model will allow the individual examination of each mode. Fluidstructure interactions are hard to model, however, the aerial flip turns (backspins) are mainly the result of differential inertial forces and are rigorously explainable by mathematical models. To describe mode (1), the current schemes can take a number of experimentally validated and successful forms [10]. Therefore, the focus will be mode (2) and that what would be the best strategy to prepare for landing on steel frames upsidedown.
One may raise the concern that how this transition – i.e., mode (1) to (2) – will occur in the first place, or, how the aerodynamic forces could vanish quickly? I will draw inspiration from biological examples and that there are a variety of deployable morphofunctional systems to consider [14]. In bats, the transition occur similar to birds by facing the inside of the wings towards incoming flow to initiate stalled conditions at which point the lift force vanishes followed by swiftly collapsing the wings at the onset of the landing maneuver, which almost likely completely zeros the effects of delayed stall and aerodynamic resistance during the backspin motion.
Without such a mediolateral movement in wings, at high angle of attack maneuvers the aerodynamic forces will be roughly normal to the wing surface at all times resulting in large pressure forces that dominate the shear viscose forces acting parallel to the wings. Because fixedwing systems have no retraction mechanism in their aerodynamic control surfaces they cannot easily backspin. Rotorywing systems to be able to preform backspin, they require thruster units that can generate positive and negative aerodynamic force relative to the body and that can be achieved through a number of costly approaches. Another issue concerning fluidstructure interaction can arise here. The ”ceiling effects” are unknown for us but I believe this phenomena can cause instability in a similar way that ground effects can lead to control issues and instability for nearground hovering rotorywing systems.
An informal comparison of a landing bat and a perching rotorcraft is meaningfully relatable to the comparison between the performance of a pitcher and a robotic arm in throwing the baseball in a baseball game. One employs natural dynamics to perfection to pitch a ball the other is carefully restricted to its predefined joint trajectories completely suppressing the natural dynamics. Of course, unlike the fullyactuated manipulator, a quadrotor is underactuated and its internal dynamics is invariant of the supervisory controller which allows for a limited contribution from the natural dynamics.
Iii Brief Overview of ”Harpoon” System
The hypothesis explained earlier is tested with the help of a selfsustained ballistic object called ”Harpoon” (shown in Fig. 1). This robot is considered as the landing gear in a morphofunctional machine. Here, the objective is to reorient the robot towards an imaginary ceiling surface and extend (or shoot) the landinggear towards the ceiling. This scenario is designed after carefully observing the highspeed images of bat landing maneuvers and once successfully reconstructed it will capture the dynamics once inspected by [3, 12, 2, 8].
Harpoon platform is lasercut from a thin carbon fiber plate and hosts: two actuators, two motor controller, a receiver and a processor. The overall design allows for a limited selfsustained sensing, actuation, computation and communication between the robot and master computer. A brushless motor adjusts the angle of a bob (appendage) with respect to the body. This bob is designed to precisely capture the body to appendage ratio found in bats. A trigger mechanism releases the landing gear once it is faced towards the surface. Since the projectile is charged with massive elastic energy beyond the output performance of the trigger actuator, a mechanism with 1:40 mechanical advantage is designed and fabricated at submillimeterscale resolution as shown in Fig. 1.
The orientation angles including the Euler parameters roll, pitch and yaw are estimated at 200 Hz using an OptiTrack system with six cameras in a 8 m
^{3} space and are fed back to a discrete controller with known asymptotic stability properties. Then, the computed control inputs are fed back to the robot in a wireless communication bus in Ethernet frames. This architecture allows for a robust communication for fast sensing and actuation. The computations take place offline and it allows for an expensive realtime processing otherwise impossible with the limited onboard processing capabilities.Additionally, a computerautomated mechanism, called releaser, ensures the robot is released with zeroangular momentum, therefore, the ballistic motion will be of zeroangularmomentum nature [1]. This condition may put restrictions on the boundaries of the trajectories from mode (1), however, a change in the mechanism design can easily provide nonzeroangularmomentum turns, which can easily allow for the study of nonzeroangularmomentum backspins. To completely suppress the aerodynamic effects, wings are excluded from the model and the host morphofunctional machine will be unveiled in the future works.
Iv Manipulation of inertial dynamics
The system dynamics and the equations of motion are described in a coordinatefree fashion, i.e., a system on manifold without considering any local coordinate charts. The kinetic energy is considered as a Riemannian metric and in writing the discrete Hamilton’s principle associated Riemannian interconnections are considered. The discretetime model on manifold with bodyfixed forces and invariant kinetic energies is erected based on the Lie group variational integrator evolving on which unlike a rigid body results in a nontrivial dynamical model.
The numerical integrator obtained from the discrete variational principle exhibits excellent geometric conservation properties [6] and because of this it is used here. This means that the robot attitude automatically evolves on the rotation group embedded in the space of special orthogonal matrices whereas the required angular velocity of appendage is computed at the level of the Lie algebra and the matrix exponential are employed to update the reference solutions.
Last, an almostglobally stabilizing, geometric discretetime controller based on [13] will be applied to regulate the attitude in preparation for landing.
Iva Notation and Assumptions
It is assumed that a body coordinate frame is fixated to the robot and it position, orientation and angular velocity are denoted by , and , respectively. Another body coordinate frame is attached to the appendage where the orientation and angular velocities are denoted by and , respectively. Therefore, the configuration space (Cspace) is . The physical properties including the inertia matrix and mass are denoted by and . A small change in is denoted by and denotes variations on . The Lagrangian functional, kinetic and potential energy are denoted by , and , respectively. Other widely used notations such as for wedge operator, for a discrete value of at the kth sample, for the trace of a matrix, as the transpose operator, as the eigenvalues of a matrix and as the variation operator are adopted. The control input is denoted by
, the identity matrix is denoted by
, is the zaxis unit vector and are the distance between the joint and the center of mass in body and appendage.IvB MoserVeselov Description of the Discretetime, Constrained Dynamics
Sensing and actuation occur in discretetime, therefore, the discretetime version of the dynamical system on is considered and that arises the constrained dynamics once investigated by [9]. Indirect methods based on Rodrigues’ approach to approximate matrix exponential and firstorder approximations of variations on Lie groups are considered as well which will be explained in this section.
IvB1 Schurform update law to and
The discretized version of Hamilton’s principle determines the equations of motion from a variational problem. The Lie group variational dynamics for the system of body and appendage when external forces are absent are given by:
(1)  
(2)  
(3)  
(4) 
where allows for the inertial dynamics contributions from the appendages. In Eqs. 1, 2, 3 and 4, the solution to the discretetime dynamical system can be resolved first by obtaining an answer to , which involves solving the continuoustime algebraic Riccati equation (ARE) given by Eq. 1. However, the solutions to the ARE are not unique and for these solution to belong to special orthogonal solutions () the convex constraint given by Eq. 4 must be satisfied [9]. In general, a change in the orientation of the body cannot be uniquely mapped to as are minuscule changes in . To obtain , the steps from [4] are applied. The Hamiltonian matrix is given by
(5) 
and after synthesizing the Schur form of the Hamiltonian Matrix , is partitioned as . Hence, the special orthogonal solutions to the ARE are given by
(6) 
IvB2 Rodrigues’ approximation to the variations in and
In addition to the Schurform solutions, another approach is considered. In this way, special orthogonal solutions to Eqs. 14 are iteratively resolved by expressing as an exponential form of where belongs to . The Rodrigues’ formula for the changes in the body orientation gives
and the ARE given in Eq. 1 is rewritten into the equivalent vector form given by
Rodrigues’ approach was considered in addition to the Schur method to provision for the numeric difficulties endured when attempted to use Eqs. 14 in a realtime reference trajectory governor. The Eqs. 14 yield the symplectic geometry of the problem in hand and is followed by computational issues about solutions in finding invariant subspaces of a Hamiltonian matrics. Other methods that could take into account the particular structure of Hamiltonian matrices are not less costly. Worthy of noting is that it was noticed – in this experimentally motivated work – extra care must be paid to as not all of the sampling rates (time intervals denoted by ) can guarantee the existence and uniqueness of the solutions in Eqs. 14. Intuitively, this makes sense as a relatively large time interval leads to an approximations of that violates the topological structure of . The following equation can help observe this problem
(7) 
Obviously, cannot take any arbitrary number. This equation is the direct result of defining such that . Using the kinematic relationship and considering an approximation for Eq. 7 is obtained.
Almost all of these methods resolve unique positive definite solutions for ARE when sample intervals are small, however, showing how the size of the sample intervals can violate the Hamiltonian matrix with no imaginary eigenvalues is well beyond the scope of this work.
Although the sufficient and necessary convex quadratic constraint in Eq. 4 is a strong condition that secures the existence of such a solution, its dependence to the sampling time intervals can be limiting.
When is positive definite then Eq. 1 has a unique positive semidefinite solution , i.e., has positive real parts. Consequently,
is the special orthogonal matrix being sought. When
is positive definite an invertible matrix
exists such that . Consequently, and will be stabilizable and detectable, respectively, which results in with positive real parts. Since are the form they have negative real part. Therefore, and are stabilizable, that is, is stabilizable and is detectable.IvB3 Baseline coordinatefree system on manifold
Based on Hamilton’s Principle, the action integral is given from the Lagrangian and its variations are zero. Because the configuration group is the Lie group, this variation should be consistent with this geometry, therefore, the varied in continuoustime mode is given by
and in discretetime mode is given by
Consequently the variations of the action integral are obtained. In driving the variations of continuoustime action integral, is used. Note that the continuous and discretetime variations vanish at the time boundaries, that is, and for the continuous system and and for the discrete system. From Hamilton’s principle, the above continuous and discretetime variations in the action integral should be zero for all variation of Lie algebra. Hence, the continuous time system on manifold is given by:
(8) 
where and are the states and control inputs, respectively; and are obtained after sorting the varaitions and applying some basic algebra. The symmetric matrix is given by
(9) 
and is given by
(10) 
where , , , and for . The discreteization of this continuous dynamics (Eq.14) is achieved by applying the discretetime variation explained above and applying Hamilton’s principle. Next, I will briefly describe the closedloop feedback used to regulate the body orientation.
IvC Discrete feedback control
A scalar measure of the attitude error is given by the rotation angle about the eigenaxis needed to rotate the robot form to the desired attitude and is given by . This metric is applied to indicate the time of launching the landing gear.
The controller applied here does not require any knowledge in regard to the mass and inertia of the robot body and appendage. A major limitation in designing the closedloop feedback for this constrained dynamics is that not all of the arbitrary reference trajectories can be enforced by the tracking controller and the tracking problem is subject to the constraints given by Eq. 4.
The issue of constraint satisfaction in nonlinear systems affine and nonaffineincontrol and its separation from the issue of closedloop system design have been addressed very extensively. Almost all of these works do so by introducing a reference governor, an auxiliary nonlinear device that operates between the reference command and the input to the closedloop system. A similar approach to [5] is adopted here.
The control input with known stability properties [13] is given by
(11) 
where is the feedback gain and can be calculated by the following relation,
(12) 
where is a symmetric positivedefinite matrix and is the reference input to the body orientation and is resolved by finding in a nonlinear optimizer subject to the following constraints:
(13)  
(14)  
(15)  
(16) 
Note that is the desired body orientation before considering the convex constraint given by Eq. 4.
V Numerical & Experimental Results
The appendage inertial model is scaled to account for 16% of the total body mass according to [8, 7]. A total estimated body weight of 30 gr is considered, which is obtained after assuming a density of according to [7], while equally distanced bobs from the body joint (shown in Fig. 2) are considered. A simple body geometry is considered for a seamless inertia matching.
When only the dynamics restricted to the longitudinal plan is considered the Cspace is the twosphere , where is a unit vector (Note that the body and appendage angles are embedded in ). The tangent space at , namely , is identified with where represents the angular velocity.
A fixed desired body orientation is constructed according to the Euler convention by ramping the pitch angles to 120 and 180 deg and updated by the reference governor subject to the constraints explained above (shown in Fig. 3(a) and 3(b)). Samples (purple dots) are considered by the optimizer at the neighborhood of the desired angles (black dots) before resolving the actual trajectories. Note that due to issues in the cluttered raster image some of the purple dots are plotted. In Fig. 3(c), the orientations of a 3Drigid body representing the body (note that the appendage is not shown) are shown along the time axis for a simultaneous tracking of roll and pitch angels 20 and 60 deg, respectively. In Fig. 3(d), the simulated Euler angles are shown.
In experiment, only the longitudinal motion is considered. The framework was tested in a series of experiments devised with an offboard computer, motion capture system and highspeed imaging equipment. In Fig. 3, the snapshots from one experiment are shown. The computerautomated releaser mechanism releases the robot with zeroangular momentum in a free fall configuration while the motion capture system captures the body and appendage orientations. The commanded (shown in Fig. 4) are resolved and transmitted to the robot for the desired body orientation angle 180 deg.
Vi Conclusion
The ability of animals to reorient their bodies midair without the benefit of net external torques is known yet remarkably has been overlooked in aerial drone design. Instead, current designs consider thrust vectoring for agile aerial robotics, which is known to be limited based on aerodynamics laws. This work used a tiny robot called Harpoon with prohibitive restrictions associated with size and payload not allowing for multithruster designs typically found in quadrotors and demonstrated that extremely fast body reorientation and preparation for upsidedown landing unique to bats is possible through the closedloop manipulations of inertial dynamics in the robot. To prepare for landing, a rubberbandpropelled landing gear design and associated trigger mechanism was proposed. The closedloop manipulations of inertial dynamics took place based on a symplectic description of the dynamical system (body and appendage), which is known to exhibit an excellent geometric conservation properties.
References
 [1] A distributed control model for the airrighting of a cat. pp. 9 (en). Cited by: §III.
 [2] (201303) Glide performance and aerodynamics of nonequilibrium glides in northern flying squirrels (Glaucomys sabrinus). Journal of The Royal Society Interface 10 (80), pp. 20120794. External Links: Link, Document Cited by: §III.
 [3] (2010) Fruit flies modulate passive wing pitching to generate inflight turns. Physical review letters 104 (14), pp. 148101. Cited by: §II, §III.
 [4] (1986) A symplectic qr like algorithm for the solution of the real algebraic riccati equation. IEEE Transactions on Automatic Control 31 (12), pp. 1104–1113. Cited by: §IVB1.
 [5] (2002) Nonlinear tracking control in the presence of state and control constraints: a generalized reference governor. Automatica 38 (12), pp. 2063–2073. Cited by: §IVC.

[6]
(2006)
Geometric numerical integration: structurepreserving algorithms for ordinary differential equations
. Vol. 31, Springer Science & Business Media. Cited by: §IV.  [7] (201105) Wholebody kinematics of a fruit bat reveal the influence of wing inertia on body accelerations. Journal of Experimental Biology 214 (9), pp. 1546–1553 (en). External Links: ISSN 00220949, 14779145, Link, Document Cited by: §V.
 [8] (2008) Kinematics of slow turn maneuvering in the fruit bat Cynopterus brachyotis.. The Journal of experimental biology 211 (Pt), pp. 3478–3489. External Links: Document Cited by: §III, §V.
 [9] (1991) Discrete versions of some classical integrable systems and factorization of matrix polynomials. Communications in Mathematical Physics 139 (2), pp. 217–243. Cited by: §IVB1, §IVB.
 [10] (2017) A biomimetic robotic platform to study flight specializations of bats. Science Robotics 2 (3), pp. Art–No. Cited by: §II.
 [11] (2012) Upstroke wing flexion and the inertial cost of bat flight. Proceedings of the Royal Society B: Biological Sciences 279 (1740), pp. 2945–2950. Cited by: §I.
 [12] (2009) Bats go headunderheels: the biomechanics of landing on a ceiling.. The Journal of experimental biology 212 (Pt), pp. 945–953. External Links: Document Cited by: §III.
 [13] (2009) Inertiafree spacecraft attitude tracking with disturbance rejection and almost global stabilization. Journal of guidance, control, and dynamics 32 (4), pp. 1167–1178. Cited by: §IVC, §IV.
 [14] (2003) Deployable Structures in Biology. In Morphofunctional Machines: The New Species, F. Hara and R. Pfeifer (Eds.), pp. 23–40 (en). External Links: ISBN 9784431678694 Cited by: §II.
Comments
There are no comments yet.