Design and Analysis for a Three-Rotational-DOF Flight Simulator of Fighter-Aircraft

Most of researchers focused on traditional six degrees of freedom (DOF) Stewart flight simulator, which can not be adaptive in fighter‐aircraft flight simulator. A three rotational DOF flight simulator of fighter‐aircraft based on dou‐ ble parallel manipulator and hybrid structure is presented. The flight simulator is composed of two identical 3‐RRS (revolute‐revolute‐spherical) spherical parallel manipulators and one cabin, called Twins. The cabin has an additional independent DOF for 360° continuous rotation, so it can be applied as a flight simulator for a fighter‐aircraft to achieve spin maneuvering. Because of the introduction of the hybrid structure and double parallel manipulator of the mechanism, the redundancy exists with respect to both kinematics and actuation. Kinematics is carried out and Jaco‐ bian matrix is established by means of screw theory. The inverse kinematics is given out by the analytical method. 64 groups inverse solutions are showed in a table by permutation. Forward kinematics is solved by an effectively numeri‐ cal method. The forward numerical method is realized based on the analytically inverse kinematics and Jacobian matrix. The numerical examples show that the forward numerical method can be used in real‐time control. The rolling motion is considered in forward kinematics and a numerical example is given out. The proposed flight simulator can spin and there are three rotational DOF with a hybrid structure so that the novel flight simulator can be used in the field of the fighter‐aircraft for pilots to train.


Introduction
Flight simulators are devices in which air crews and pilots can train without the use of an actual aircraft [1]. In flight training, they are used mainly to reduce costs and increase safety. In their most sophisticated form, they simulate an aircraft's vehicular motion, instrumentation and sounds, gravitational forces, radar and electro-optical sensor displays, and out-the-window views. According to the USA's Federal Aviation Administration (FAA) regulations, any device called a Flight Simulator must have at least one motion platform, otherwise it can only be termed a Flight Training Device [2]. Therefore, the vehicular motion platform in flight simulator is one of the most important parts.
Since Stewart's initial use of a 6-DOF parallel manipulator as a flight simulator in 1965 [3], this approach has become standard. Over the past five decades, the Stewart parallel manipulator has been used to make significant contributions to aeronautical research [4][5][6][7]. Even so, there is still the disadvantage of the Stewart-type flight simulator in that its posture rotation range is less than 30°. The limitation of the posture range of a Stewart parallel manipulator hinders its ability to serve as a flight simulator of a fighter-aircraft. The motion of a fighteraircraft involves continuous 360° rolling frequently. In order to achieve continuous 360° rolling, Kim et al. [8] presented an innovative motion base as a flight simulator, based on a 6-DOF parallel mechanism, called Eclipse-II. The Eclipse-II allows continuous 360° rotation in A, B, and C-axes as well as translational motions in X, Y, and Z-axes. However, the rotations of the Eclipse-II parallel

Open Access
Chinese Journal of Mechanical Engineering manipulator are achieved by two circular guide rails, which severely affect the dynamical performance. As a matter of fact, in both the Stewart and Eclipse-II simulators, translational movements are nearly inoperative. Low-cost flight simulator with a reduced-DOF (less than six) platform was proposed by Pouliot and Gosselin [9], which revealed result comparable to those of a 6-DOF Stewart platform. Subsequently, Shui et al. [10] proposed a more advanced and innovative reconfigurable spherical motion generator to enable continuous spherical motion of the flight simulator. The mechanism enables unlimited workspace with respect to 3-DOF spherical motion with rapid, continuous, and precise motion capabilities. However, it is actuated by an electromagnetic motor with a highly complicated configuration that causes the volume of manipulator to be huge. So far, no parallel manipulator can be used as a flight simulator to accomplish a 360° continuous rotation in practice.
Although it is not necessary to be able to continuously rotate in all orientations, a flight simulator of a fighteraircraft needs to be able to roll in an additional 360° to perform most motions of a fighter-aircraft. Therefore, the purpose of this study is to achieve continuous 360° rotation of a fighter-aircraft flight simulator.
Many researchers have focused on hybrid platforms such as the famous Tricept by Neumann [11] and Trivariant by Huang et al. [12] that provide a good workspace volume. They are typically a combination of a parallel manipulator and a serial manipulator. As a hybrid alternative, a double parallel manipulator was proposed for enlarging the workspace by Lee et al. [13,14]. Tsai et al. [15] used a double parallel manipulator to enlarge the workspace of the manipulator and analyzed its kinematic properties. Furthermore, many researchers have noted that the redundancy of a double parallel manipulator can improve its abilities and performance, for example, enlarging the volume of the workspace. There are two main types of redundancy for the parallel manipulator: (a) kinematic redundancy and (b) actuation redundancy [16,17]. Zanganeh and Wang et al. [18][19][20] studied the kinematic redundancy of a parallel manipulator by analyzing its kinematics, and pointed out that the extra DOF not only allow for execution of the original output task, but also additional tasks such as increasing the workspace. It is notable that the Eclipse-II is a redundant actuation parallel manipulator. Redundant actuation has been proven to be a good method to enhance performance of parallel manipulators. Nokleby et al. [21] proved that redundant actuation can improve force capabilities. Kim et al. [22] investigated the redundantly actuated parallel manipulator and proved that redundant actuation of a parallel manipulator not only can improve force capabilities, but can also enhance the stiffness of a manipulator.
Li et al. [23] derived the conclusion that redundant actuation has little effect on the stiffness when the actuators are on prismatic joints and enhances stiffness value greatly when the actuators are on revolute joints.
Although the novel mechanisms are proposed one by one, they don't adapt for the flight simulator. Regarding existing parallel manipulators, the 3-DOF spherical parallel manipulator is a compact configuration with large rotational posture [24]. Therefore, given the above conclusions, a flight simulator of a fighter-aircraft based on a double and hybrid 3-RRS [25][26][27][28] spherical parallel manipulator was chosen for this study.
The remainder of this paper is organized as follows. In Section 2, the structure of the Twins flight simulator is described. In Section 3, the kinematic properties of the Twins are discussed, including the development of the direction-cosine matrixes in Section 3.1 with analytical spherical theory, the analysis of DOF in Section 3.2, and the development of the Jacobian matrix via screw theory in Section 3.3. Inverse kinematics and forward kinematics are conducted in Sections 4 and 5. Numerical examples are provided in Section 6. Lastly, the conclusions are discussed.

Architectural Description
As shown in Figure 1, the flight simulator is a hybrid manipulator developed on the basis of two identical 3-RRS spherical parallel mechanisms and one cabin with an independent DOF. Therefore, it has kinematic redundancy and a redundantly actuated manipulator as well. The basic makeup of the whole platform consists of two fixed platforms, one cabin, six RRS legs, two flange-bearing-gear decelerators, and eight motors. In the fixed coordinate frame, the x-axis is normal to plane C 1 ON, and extends outwards. The z-axis superposes on line MN, and extends upwards. The y-axis is defined by the right hand rule. In the moving coordinate frame, the u-axis is normal to plane A 1 ON, the w-axis superposes on line QV (the points Q and V are the geometrical center of triangles ∆A 1 ′A 2 ′A 3 ′ and ∆A 1 A 2 A 3 ), and extends outwards. The v-axis is defined by right hand rule.
Looking at Figure 2, based on the symmetrical rule [29], . α 1 is the angle between the B i -axis and C iaxis, and similarly α 1 ′ is the angle between the B i ′ axis and C i ′ axis. α 2 is the angle between the B i -axis and A iaxis, α 2 ′ is the angle between the B i ′ axis and A i ′ axis. β 1 is the angle between line MN and line OC i . β 2 is the angle between line QV and OA i .
Looking at Figure 3, to achieve 360° rotation of the cabin, two flange-bearing-gear reducers are used. The outer flange is fixed to the moving platform by bolts, and the inner flange is fixed to the cabin by bolts as well. Together, the outer flange and inner flange serve as a bearing. Three axis-fixed gears are fixed to the moving platform via shafts, and the center gear is linked to a motor.

Development of Direction-cosine Matrices
The direction cosines of a vector are the cosines of the angles between the vector and the three coordinate axes. Deriving direction cosines are required in this paper to perform inverse kinematics and obtain the Jacobian matrix. A 3-RRS spherical parallel manipulator and its mirror image are shown in Figure 4.
On the basis of analytical theory of spherical space [30], the direction cosines matrices of B i , B i ′ are built first. The fixed Cartesian coordinate frame of the two 3-RRS spherical parallel manipulators is shown in Figure 4, where the x-axis and x′-axis are inverses of each other, as are the z-axis and z′-axis and the y-axis and y′-axis. The direction cosines matrices of o B i and o B i ′ in the fixed Cartesian coordinate frame O-xyz are as follows: where θ i and θ i ′ are input angles of the C i -axis and C i ′-axis, and where φ 1 = 0°, φ 2 = γ 1 = 120°, and The rotation matrix R(z, φ i ) with reference to the z-axis is as follows: In the fixed Cartesian coordinate frame, the direction cosines matrices of °C i and °C i ′ can be given directly as follows: In the moving coordinate frame, the direction cosines matrices of A i and A i ′ are expressed as follows: (10) where T is the Euler angle rotation matrix with regard to the moving coordinate frame O ′ − uvw , and is depicted as follows: where λ, ε, υ are the yaw, pitch, and roll angles (see Figure 5).

Degrees of Freedom Analysis
When the cabin is fixed, the degrees of freedom of the mechanism can be obtained by using the modified formula of Grübler-Kutzbach: where M is the number of degrees of freedom, d is the number of common order, d = 6 -χ (χ is a common constraint), f i is the number of degrees of freedom of the ith joint, g is the number of joints, n is the number of bodies in the mechanism, and η is the number of over-constraints. The screw theory is used to analyze the number of over-constraints η, and the screws in limb A 1 B 1 C 1 are shown in Figure 6.
$ r 1 is defined as the reciprocal of screw $ 1 such that The reciprocal screw $ r 1 of limb A 1 B 1 C 1 is derived using Eqs. (20) and (21): The same method is used to obtain the reciprocal screws.

Jacobian Matrix
A Jacobian matrix is a mapping relationship between the angular velocity of active joints and the angular velocity of the moving platform. The establishment of a Jacobian matrix is fundamental to the analysis of forward kinematics in this paper. There are many methods for establishing a Jacobian matrix. In general, a Jacobian matrix is obtained by virtue of derivations of constraining equations. It is troublesome to obtain a Jacobian matrix by virtue of derivations if the constraining equations are too  many. For simplicity, in this paper, screw theory is used to build a Jacobian matrix [31][32][33] for which the moving platform angle velocity $ p is derived as follows: where and where ξ 1 , ξ 2 , ξ 3 , ξ 4 , ξ 5 are angular velocities associated with unit screws $ 11 ,$ 12 ,$ 13 ,$ 14 ,$ 15 .
The next step is to establish a relationship between ξ 1 and $ p while eliminating ξ 2 , ξ 3 , ξ 4 , ξ 5 . r $ 1 denotes a reciprocal screw and is derived by means of the following rules: (1) r $ 1 is through the center of the spherical joint of the A 1 -axis and B 1 -axis; (2) r $ 1 does not intersect with the C 1 -axis. where , J q is not a square matrix, so herein, the Jacobian matrix is defined as Rewriting Eq. (32) yields

Inverse Kinematics
In the equations of inverse kinematics, the posture parameters ( , ε, υ) of the moving platform are given. The objective is to calculate the input angle θ i and θ ′ 1 , i = 1, 2, 3. In the inverse kinematic analysis, the kinematic equations of the 3-RRS spherical parallel manipulator are the same as those of the 3-RRR spherical parallel manipulator, but the 3-RRR spherical parallel manipulator is an over-constrained mechanism that is difficult to practically assemble. That is why we chose a flight simulator model based on 3-RRS spherical parallel manipulators.
When considering the independently redundant DOF of the cabin, because there are three outputs and eight inputs, it is innumerable about the inverse solutions. If we ignore the independently redundant DOF of the cabin, and fix the independent DOF of cabin, the analysis of the inverse kinematics is as follows.
x i is defined as Eq. (42) Finally By repeating Eqs. (38)-(45), θ ′ i is derived by the same method: All forms of the combination according to Eq. (46) and Eq. (47) are shown in the Table 5 of Appendix 1.

Forward Kinematics
For the forward kinematics, without considering the independently redundant DOF of the cabin, the angles of the input θ i and θ ′ i , i = 1, 2, 3 are given, and the posture ( , ε, υ) of the moving platform is determined. However, considering the redundant actuation, the DOF of the moving platform of the flight simulator is three; thus, it can only choose three angles of the input at one time. Herein, we choose θ i (0 < θ i < π, i = 1, 2, 3. ) as the inputs, and the remaining three angles of the input θ ′ i (0 < θ ′ i < π, i = 1, 2, 3 ) are obtained via inverse kinematics. In general, the forward kinematics of a parallel manipulator is fairly complicated when using the analytical method, and many solutions are derived. Moreover, in practical application, the configuration of a manipulator is just one of the forward solutions. The motion of the flight simulator is continuous and the control system is based on a dynamical model of MIMO (multiple input multiple out) that needs to obtain the posture of the moving platform in real-time using the forward equations. Thus, an effective numerical method is employed to solve the forward equations. The numerical method was first used in Stewart's parallel manipulator [35].
The basic idea is to assume initial posture angles ω 0 = ( 0 , ε 0 , υ 0 ) T of the moving platform, and the (41) (42) (44) corresponding input angles p 0 = θ 0 1 , θ 0 2 , θ 0 3 T are obtained according to the initial posture ω 0 = ( 0 , ε 0 , υ 0 ) T by inverse kinematics. At the same time, giving three initial input angles chosen from the six input angles, assuming that ′ p 0 = ′ θ 1 , ′ θ 2 , ′ θ 3 T then defining Then The initial posture becomes The above process is repeated, and the iterative equation is The iteration is stopped by the constraint condition max �θ k i < δ, i = 1, 2, 3 . �θ k i is an element of �p (k) , and δ is the permissible error. J s must be nonsingular when the manipulator performs a continuous motion and the iteration must converge to a unique solution. Eq. (53) is similar to the Newton-Raphson iteration, but the physical and geometrical meanings are more significant, and the whole process of the algorithm is different. Appendix 2 shows a flow chart of the algorithm.

Numerical Examples
The parameters of the manipulator are set in Table 1. thus proving that the algorithm is correct. Herein, we must point out that as long as ω 0 is less than the range of the workspace, the result must converge to the correct solution, and the solution must be unique. However, the algorithm is like the Newton-Rapshon algorithm in that it is sensitive to the initial values. The results listed in Tables 2 and 3 show that the algorithm is correct. For θ ′ i , i = 1, 2, 3, obtained via the inverse equations, the posture of ω 0 = (0.1, 0.1, 0.1) T is shown in Figure 8 when the inverse solution is determined by the first row of Table 1 (identical to column 1 of Table 3).
As shown in Tables 2 and 3, the results are obtained numerically by the forward algorithm and analytically by the inverse equations. The time costs of the numerical forward algorithm are less than 3 ms using MATLAB software by Intel Core i7-3.2G CPU, thus proving that the algorithm of the iteration meets the typical real-time control requirement of less than 6 ms.
To consider the independent DOF of cabin, the forward solution is obtained in two steps. The first step is to obtain the posture ( , ε, υ) T by the numerical method.     Then the posture ( , ε, υ) T is changed by the cabin's spin with reference to the w-axis of the moving coordinate. Defining φ as the angle of the cabin's spin with reference to the w-axis, the rotation matrix with reference to the w-axis of the moving coordinate is Defining ( ′ , ε ′ , υ ′ ) T as the final posture by virtue of the spin of cabin, then the results of the numerical example are shown in Table 4.

Conclusions
(1) The flight simulator for a fighter-aircraft with a hybrid configuration of symmetrically double parallel manipulators named Twins is presented. (2) Twins is a multi-functional flight simulator. It can be used as a normal flight simulator like Stewart mechanism. And when the cabin spins that will be a flight simulator of fighter-aircraft. (3) Screw theory is used to establish the Jacobian matrix that simplifies the process of establishing the Jacobian by the means of closed equations. (4) A numerically forward kinematics is adopted by the inverse kinematics and Jacobian matrix and the method is more simple.