#### Abstract

This paper presents a modelling of an unmanned airship. We are studying a quadrotor flying wing. The modelling of this airship includes an aerodynamic study. A special focus is done on the computation of the added masses. Considering that the velocity potential of the air surrounding the airship obeys the Laplace's equation, the added masses matrix will be determined by means of the velocity potential flow theory. Typically, when the shape of the careen is quite different from that of an ellipsoid, designers in preprocessing prefer to avoid complications arising from mathematical analysis of the velocity potential. They use either complete numerical studies, or geometric approximation methods, although these methods can give relatively large differences compared to experimental measurements performed on the airship at the time of its completion. We tried to develop here as far as possible the mathematical analysis of the velocity potential flow of this unconventional shape using certain assumptions. The shape of the careen is assumed to be an elliptic cone. To retrieve the velocity potential shapes, we use the spheroconal coordinates. This leads to the Lamé's equations. The whole system of equations governing the interaction air-structure, including the boundary conditions, is solved in an analytical setting.

#### 1. Introduction

The rapid expansion of the capabilities of airships in the last decade and the growing of the range of missions they designed to support lead to the design of complex shapes of the careens.

Exhaustive studies in that topic were presented by [1, 2]. Traditionally, ellipsoidal shapes are used for airships [3–5]. However, and in order to optimize their performances, different original shapes have been tested in the last years. This is due to the advances in aerodynamics, control theory, and embedded electronics. The airship studied here departs with the traditional shapes. The MC500 is a flying wing (Figure 1), developed by the French network DIRISOFT.

The MC500 is an experimental prototype for a set of great innovating airships. A precise dynamics model is needed for this kind of unmanned airships including the air-structure interaction. This will enable the elaboration of convenient algorithms of control, stabilization, or navigation of these flying objects.

The aerodynamic forces have a large influence on the dynamic behaviour of the airships. Among the aerodynamic solicitations, the added masses phenomenon is one of the most important. In fact, when hovering or in case of low speed displacement, the lift and drag forces and torques could be neglected.

The added masses phenomenon is well known for airships and similarly for submarines. When an airship moves in an incompressible and infinite inviscid fluid, the kinetic energy of the fluid produces an effect equivalent to an important increase of the mass and the inertia moments of the body. This contribution may be of the same magnitude as the terms of mass or inertia of the airship. Apart from ellipsoidal shapes [6] or thin rectangular plates [7] where the theory is well established for many decades, the analysis of this problem for unconventional shapes in preprocessing is usually done by approximate methods. We can quote the geometric method, consisting of an extension of a 2D analysis and requiring the introduction of empirical parameters [8, 9], or the fully numerical methods consisting of modeling both the airship and the surrounding air [10, 11].

As opposed to other treatments of this problem, the derivation proposed here is based on the velocity potential flow theory [6, 7, 12]. We tried to pursue the analytical study to an advanced stage with some assumptions. The shape of the careen of the MC500 is considered as a cone with elliptic section. When we consider that the velocity potential of the air obeys the Laplace’s equation, the added masses matrix could be determined by solving this equation, using the spheroconal coordinates. Solving the Laplace’s equation for such configuration leads to the Lamé’s equations. We have dedicated a section of this paper for the determination of the ellipsoidal harmonic series issued from these equations with the specific boundary conditions. Such series developed for the first time by Lamé [13] were improved particularly by Liouville and Sturm in their famous theory and by Hermite [14] and were applied in different fields. Byerly [15] and later Hobson [16] gave an important comprehensive literature about this topic.

Significant recent works, such as the works of Boersma and Jansen [17] in electromagnetic field, and those of Garmier and Barriot [18] in astrophysics could be quoted. In this work, we are trying to apply this theory for the air-structure interaction.

#### 2. Dynamic Model

##### 2.1. Kinematics

The airship MC500 is assumed to be a rigid flying object.

We use two reference frames. First an earth-fixed frame assumed to be Galilean. Then a local reference frame fixed to the airship. Its axes are selected as follows:** and are the transverse axis of the airship, : the normal axis directed downwards.

To describe the orientation of the airship, we use a parameterization in yaw, pitch, and roll. The configuration of the object is described by means of three rotations defined by three angles of orientation, that is, the yaw , pitch , and roll .

The whole transformation between the local reference frame and the global frame is given by Goldstein [19]: such as . We denote by ; .

Using the rotation matrix , the expression of the linear speed in the reference frame is given by The angular speed is expressed as follows: with It is noticed that the parameterization by the Euler angles has a singularity in .

This parameterization is acceptable because it is not possible for an airship to reach this singular orientation of 90 degrees pitching angle.

##### 2.2. Dynamics

As mentioned previously, a specificity of the lighter than air vehicles is illustrated by the added masses phenomenon. When a big and light object moves in the air, the kinetic energy of the particles of air produces an effect equivalent to an important growing of the mass and inertia of the body. As the airship displays a very large volume, its added masses and inertias become significant.

The basis of the analysis of the motion of a rigid body immersed in a perfect fluid is established by Lamb [6]. He proves that the kinetic energy of the body with the surrounding fluid can be expressed as a quadratic function of the six components of the translation and rotation velocity as follows: where is the added mass matrix due to the motion of the surrounding air and is the mass matrix of the body.

For an airship fully immersed in the air, the added mass matrix is assumed to be positive and definite. As the added kinetic energy , it will be discussed in the next section. The whole mass matrix is assumed symmetric block-diagonal matrix: For the computation of the whole dynamics model, we choose to use the Kirchoff’s equation [20]: where is the vectorial product, and are, respectively, the external forces and torques, including the rotors effects, the weight (mg), the buoyancy , and the aerodynamic lift () and drag ().

The dynamical system of the airship becomes [21–23]

*Description of the Rotors*

The MC500 has four electric engines driving rotors. Each rotor has two parallel contrarotating propellers to avoid any aerodynamic torque (Figure 2). The rotor can swivel in two directions. A rotation of angle around the axis () and a rotation of angle around an axis normal to and initially coinciding with the axis (). A fictive axis completes the rotor frame.

Let us denote by the position of the rotor . We can then define a rotation matrix between the frame and the local frame such as If we suppose that the intensity of the thrust force of the rotor is , this force could be introduced in the second member of the dynamic equation as where is a unitary vector along the axis.

The torque produced by this rotor in the centre of inertia is .

*Weight and Buoyancy*

An important characteristic of the airships is the buoyancy . This force represents a natural static lift, corresponding roughly to 1 Kg for each m^{3} of helium involved in the careen. We suppose here that this force is applied in the centre of buoyancy different from the centre of inertia :
where is the volume of the careen, is the density of the air, and the gravity.

Let us note and the force and the moment due to the weight and buoyancy.

We have

*Aerodynamic Forces *

Such as other flying objects, the airships are subjected to aerodynamic forces. The resultant of these forces could be divided into two component forces, one parallel to the direction of the relative wind and opposite to the motion, called Drag, and the other perpendicular to the relative wind, called Lift. The MC500 is designed with an original shape oriented to a best optimization of the ratio lift upon drag forces.

However, and as first study, we try to evaluate the behaviour of the airship in the case of low velocity or while hovering. In these cases, the effect of these forces could be neglected.

The global dynamic system could be expressed in a compact form as follows [21, 24]: with such as

#### 3. Computation of the Added Masses

##### 3.1. Spheroconal Coordinates

Much of the current airships developed and studied in the literature are considered as ellipsoidal. This particular shape has a wide popularity in this field due to the existence of a large and reliable knowledge and experimentations for both ellipsoidal airships and their alter ego submarines. However for our flying wing airship, we should use a more convenient approximation for the aerodynamic study.

The MC500 airship is a collection of aerodynamic profiles, with symmetry around the - axis. Each transverse section of the careen parallel to the plan - gives roughly an ellipse. This motivated us to model as first assumption the airship as an elliptic cone (Figure 3).

To take into account the interaction of the airship with the surrounding fluid, a model of the flow is needed.

A variation of this study may be performed by calculating the fluid pressure around the airship through the Bernoulli equation [25].

Here, we use the potential flow theory with the following assumptions.(a)The air can be considered as an ideal fluid with irrotational flow and uniform density , that is, an incompressible fluid.(b)A velocity potential exists and satisfies the Laplace’s equation throughout the fluid domain: and satisfies the nonlinear free surface condition, body boundary condition, and initial conditions.

Finally we suppose that the velocity of the air is null far from the airship: (c) The velocity of the fluid obeys the equation:

In his study [6], Lamb proves that the kinetic energy of the fluid surrounding a moving rigid body can be expressed as a quadratic function of the six components of the translation and rotation velocity and can also be expressed as function of the velocity potential of this fluid by the following relation:
is an outward normal vector. can be splitted on
where appears as a collection of six spatial potential Φ_{i} functions of , , that verify the Laplace’s equation. To extract the terms of the added mass matrix, we should define the spatial potential of the moving fluid. In some cases it could be determined using geometric means [8, 26]. Here we choose to use the potential flow theory. We can then express the components of the added mass matrix in function of the velocity potential flow of the air surrounding the airship such as
For solving the Laplace’s equation (3.1) and according to the configuration of the careen, we use the spheroconal coordinates [16]. This assumption seems acceptable in that stage.

A first parameterisation by the coordinates () is used, for example, by Boersma and Jansen [17], Kraus and Levine [27] (Figure 4), where is the distance of a point to the origin, is the azimuthal angle (), and is the longitudinal angle. represents the external surface of the elliptic cone.

This parameterisation has the advantage to be physically significant. However, in our case, it leads to intractable calculations. We prefer to use the parameterisation given by Hobson [16]. The equation of the surface of an elliptic cone is given by G. A. Korn and T. M. Korn [28] such as , , and are parameters that define the elliptic cone.

With is always the distance of a point to the origin. We define then the Cartesian coordinates (, , ) as By analogy with the azimuthal angle in the previous parameterisation, we can introduce an angle defined by with .

The surface of the elliptic cone is now defined by . By varying this parameter, one can see its influence on the shape of the cone in Figure 5.

To express the components of the fluid velocity in the conical elliptic reference frame, we define three unitary vectors , , . With are the scale factors. One can easily verify that these three angles form an orthonormal basis, and the velocity of the fluid could be expressed as The surface of the elliptic cone is defined by ; appears as the normal vector to the surface in a given point.

##### 3.2. Laplace’s Equation

With the spheroconal coordinates, the Laplace’s equation (3.1) becomes [16] with To solve the Laplace’s equation (3.13), we use the traditional separation of variables: Replacing (3.16) in (3.12) gives This equation can be splitted on three separated differential equations: And by replacing , , and by their expressions in , one can obtain The solutions of these equations are ellipsoidal harmonics or Lamé’s functions. The Lamé’s equations are usually written in a general form as For a given integer, we can find such that the particular solution of (3.19) can be written as [18] Those solutions are the Lamé’s functions of first kind with where depends on the integer , as shown in Table 1: is called the principal product.

There are four different Lamé’s functions that differ from their characteristics. We present them in Table 1.

is an eigenvalue. We remark that corresponds to a value of the parameter in (3.19), which determines as many equations as there are values of .

For a given value of , the unknowns are the coefficients (3.23) and the parameters . We consider as a first step the functions ; the other functions could be defined similarly.

If we introduce the expressions (3.20) and (3.21) in (3.19), we obtain the following recurrence’s relation: or Knowing that , the iteration will stop at rank . And if we introduce the conditions: the whole recurrence’s relations (3.28) can be written in a matrix form such as or is a square matrix of dimension . The vector is an eigenvector associated to the eigenvalue .

There are basically eigenvalues and eigenvectors. In fact we show that it is possible to find a diagonal matrix and a symmetric matrix such as with

The matrix is diagonalizable such as , where is diagonal and is an orthogonal matrix. Accordingly, is diagonalizable and admits separate eigenvalues associated to eigenvectors , so functions for .

has the same eigenvalues than . Knowing that is symmetric, its eigenvalues and eigenvectors are obtained conventionally using the algorithm QR. Then the eigenvectors of are given by .

The computation of the three other functions , , and is in the same manner.

Each of the functions , , , and admits exactly zeros in the interval .

For , the value of corresponding to is equal to corresponding to , and the values of corresponding to and are identical.

However, the computation shows a significant numerical instability if the value of increases beyond 10. This is due to the computation of the coefficients in the expression (3.21).

A technique was proposed by Dobner and Ritter [29] to stabilize such computation. They proposed to use another expression of the polynomial such as This variation gives a more stable computation.

The product satisfies the potential equation (3.19) within the elliptic cone. However for the external space other kinds of solutions are needed which vanish at infinity.

There is another kind of function solution of (3.19) which vanishes at infinity. These functions are called Lamé’s functions of second kind. They are defined by Hobson [16]: For each value of we have a function . We have then Lamé’s functions of second kind.

Recalling that surface of the elliptic cone is given by the relation , the spatial velocity potential flow Φ, solution of (3.13), is then [16] where are constants to be defined using the boundary conditions.

##### 3.3. Boundary Condition

In addition to the relations (3.1)–(3.3) the velocity potential flow Φ verifies a kinematical boundary condition on the surface of the elliptic cone such as is an outward normal vector to the surface: The relation (3.33) can then be written as with Let us denote . And suppose that and with is the complete elliptic integral .

For a function of , , and on the boundary surface we will have We can then solve the problem for the potential Φ in a given point in the space with the boundary conditions defined on the surface of the airship .

The relations (3.32) and (3.35) give Using (3.37) and (3.38) we obtain and then Φ and by the way .

We can now deduce the different components of the mass matrix .

Finally, the vector of gyroscopic forces can then be expressed in a developed form as This leads to the developed dynamic model:

#### 4. Results and Discussion

In this section we present the computation results of the added mass matrix of the MC500. Since any comparison with a similarly shaped airship is impossible at present, we conducted the calculation by the geometric method that we present here briefly in order to compare the results of our method.

Geometric method is well known and discussed extensively elsewhere [6, 8, 12, 20]. In a 2D analysis the planar coefficients are established for the standard shapes (rectangle, circle, ellipse, etc.). Following the exhaustive study of Brenner [8], we model our flying wing as a truncated cone () with elliptic section.

The airship is divided into a dozen cross-sections to optimize the inclusion of changes in transverse dimensions in the 3D calculation.

The computing of the terms of the added masses matrix can be seen, for example, as where is a 2D added mass coefficient for the forward motion.

According to the large difference of size between the diagonal and off-diagonal terms, we will neglect these last terms, keeping only the term .

##### 4.1. Results

We present here some characteristics of the geometry of the airship: m; m; m; m; m; m; volume m^{3}. Numerical results are presented in Table 2.

##### 4.2. Discussion

The first results described here show that thanks to the application of the velocity potential flow theory to this unconventional airship it is possible to obtain reasonable values of the added masses matrix. To our knowledge, this is the first attempt to compute the added masses using this technique for an elliptic cone-shaped airship.

Some differences can be observed between the two methods to certain terms of the added masses matrix. Experimental studies on the prototype nearing completion will confirm the accuracy of our method.

Although some geometric assumptions are made, it nonetheless demonstrates the capability of this method to compute interesting values of the added masses matrix.

Validation of our technique will allow a good estimation of the added masses matrix in preprocessing phase for this kind of airship. The development of an airship is usually done by iterative techniques (model, stability, propulsion, sizing of the control, rudders, etc.); obtaining a fairly accurate aerodynamic model early design allows better refinement of the final model the airship.

#### 5. Conclusion

In this paper a dynamic model of an unconventional airship is presented. The original shape of the careen induces a necessary reformulation of the dynamic and aerodynamic study of these flying objects. A special focus is put on the computation of the added masses. According with the original shape, an assumption of elliptic cone was made to define the careen. The ellipsoidal harmonics and the Lamé’s equations are revisited for the analysis of the velocity potential flow, according to the constraints of fluid-structure interaction. The first results seem promising.

#### Nomenclature

: | Vector position of the origin expressed in the fixed reference frame |

: | Vector orientation of the pointer in regards to and given by the Euler angles |

: | Vector attitude compared to |

Velocity vector compared to expressed in | |

: | Velocity vector expressed in |

: | Vector of angular velocities expressed in |

: | The velocity vector |

The mass of the airship | |

The identity matrix . |