Modelling and Simulation of an Anti-Tank Guided Missilel

Modelling and Simulation of an Anti-Tank Guided Missilel

Modelling and Simulation of an Anti-Tank Guided Missile

M. Abdallah, A.N. Ouda,

Department of Guidance and Control Technical Research Center

Cairo, Egypt

Abstract In order to investigate and analyze the dynamic behavior of an Anti-tank guided missile system the aerodynamics should be studied and determined precisely as they are the main corner stone in designing a mathematical model of any even flying or ground vehicle. This paper presents a sufficient and a complete analysis of the total forces and moments affecting upon a missile body causing its spatial motion. Missile DATACOM is a software adopted to get a complete picture about the force and moment coefficients as a function of Mach number, angle of attack, side slip angle and the control fin deflections[1], then a complete nonlinear mathematical model is developed for our anti-tank guided missile. The simulated results yield from the mathematical model of the missile are validated with real shooting data, the comparison shows that the missile model is close to the real missile behavior. The simulated missile is then evaluated with different target scenarios (moving and stationary) and the missile is able to hit the target in each case. Lastly more investigation should be carried out on the control loop to enhance the missile behavior against other scenarios like uncertainty in aerodynamic calculations, unmodelled dynamics and degradation in missile main thrust.

Keywords Nonlinear, Modeling, Simulation, 6DoF, Anti-Tank Missile, Skid to Turn Missile..

1. INTRODUCTION

   The first problem faced by the designer of a missile guidance system is that of translating the missile tactical problem into specifications for the guidance system design. Simulation of the system by analog or digital computer is employed as an aid to the process of the missile design i.e. simulation is a continuing aid to designer throughout the duration of the missile program[2]. When the guidance system has been designed it can be evaluated by real flight recorded data. The evaluation process and redesign are carried out utilizing simulation on computers and in turn computers are considered to be the main component in missile design and implementation process.

   Simulation is a process of imitating the behavior of the actual missile system by the behavior of a set of physical equations governing the guided missile motion which solved on computers using the available numerical method yielding the behavior of the missile. This set constitutes a set of complicated differential equations involving nonlinearities of many kinds. These nonlinearities arise from aerodynamic equations of the missile or from such mechanical effects as limiting, dead- bands, backlash, and hysteresis effects. Therefore, the solution of these equations can be carried out either through [3]:

   • Reducing the complexity by considering some simplifications, taking into consideration that these simplifications must be similar to the full system.

   • Utilizing a great developments in computation means to solve this kind of complexity.

   Therefore, by the aid of appropriate computer the differential equations describing the guided missile system can be solved to give a vision to the behaviour of the missile depending on the approximation carried on the differential equations[4]. The solution of these set of the differential equations depending on the determination of force and moment coefficient which aerodynamically control the missile[5]. This paper is divided into three main sections as follows: the first section provide the total forces and moments affecting upon the missile also the mathematical representation for each force and moment. The coordinate reference frame is illustrated as well in this section. The second section shows the application of Missile DATCOM on our missile system yielding the coefficients of both forces and moments depending not only the physical dimension of the missile but also on some other parameters like Mach number, angle of attack, side slip angle and finally the deflection of the control fin located in the aft part of the missile. The validation of the simulated missile data with real shooting data is presented in the third section, the engagement of the missile with a target at different velocities is presented in this section as well. The last section includes the conclusion of this work and the investigations that have to be carried out to enhance the missile model.

2. PROPOSED MODEL

   In this section, the physical configuration of the missile and the mathematical representation of both aerodynamic forces and moments acting on missile fuselage and the equations of motions are described as follows:

   1. Body Model

      The missile body has a plane-symmetrical configuration with a set of four orthogonal rear fins in X configuration responsible for the control action on the missile body. These configurations are also useful in case of deploying guidance and control techniques to maintain stability control and enhance hitting accuracy [6]. Fig. 1(a, b) shows the airframe shape of the missile fuselage and the X-form fin configuration of the bodys tail.

      Fig.1 (a)

      1. Acting Forces and Moments

         The aerodynamic forces which does not pass through the c.g. appears on the control surfaces when they are deflected and originates a moment resulting in missile rotation around c.g. this rotation of the missile frame with respect to the trajectory

         changes the missile angle of attack () and the side slip angle

         () and in turn the missile is subjected to lifting and lateral

         forces[10]. So a complete investigation on the forces and moments acting and affecting the missile during its flight is represented in the upcoming subsections.

         1. Thrust Force (T)

            This force is acting on the missile due to the thrust of the rocket motor even booster or sustainer and the only main component is acting along the missile longitudinal axis

            (1 )[11]

            i.e.

            T=1 (4)

         2. Missile Weight (G)

            The missile weight is determined by the following relationship:[11]

            Fig.1 (b)

            Fig. 1. (a) The missile body airframe shape, (b) X-form fin configuration

            Where:

            = (5)

   2. Reference Frames

      To formulate a mathematical model for a missile, the equations of motion should be defined these are associated to a set of reference frames or coordinate systems. On account of this, the following three orthogonal coordinate systems were used to formulate the mathematical model[7]:

      1. OX1Y1Z1 body-fixed reference frame with its origin at the missiles centre of mass.

      2. OXYZ velocity reference frame.

      3. OXgYgZg Earth-referenced system with its origin at the missile centre of gravity at the instant of lunching.

      The body coordinate system is transformed to the fixed

      (earth) coordinate system throughout the Eulers angles (, ,)

      .instantons missile mass

      the vector of gravity acceleration

      1. Aerodynamic forces

         It is usually distributed into velocity coordinate system, which related to the directions of the missile. The components of aerodynamic forces are given by[12]:

         =

         =

         = (6)

         Where:

         1. drag force.

         2. lift force.

         3. Side force.

            such that the transformation matrix (

            ) equals8]:

            S characteristic area

            T =

            g

            CC SS SCC SC + CSS C CC CS

            CS CS + SSC CC SSS

            q Dynamic pressure

            The aerodynamic force coefficients can be written as

            = (, , , , )

            Furthermore the Velocity and Body coordinate systems are

            transformed using the Side Slip angle () and the Angle of

            Attack (), such that the transformation matrix from Body to

            Where:

            = (, , , , )

            = (, , , , ) (7)

            b

            T =

            CC CS S

            S C 0

            M Mach number =

            . missile velocity

            SC SS C

            Velocity coordinate systems are as follow [9].

            Also the Velocity and ground reference frames can be related

            to each other by Eulers angles (, , ) and consequently the

            sound velocity at the missile position

            rudder fin deflection

            . elevation fin deflection

      2. Acting Moments on the Missile Body

   transformation matrix from Velocity to Ground is

   +

   [3]

   The main affecting moment on the missile body (X form missile type) is the aerodynamic moment. The

   T =

   g

   Where,

   +

   C, S are cos, sin etc

   aerodynamic moment originates from the aerodynamic forces that dont pass through the missile c.g. This moment has three components; pitch, yaw and roll moments. Usually the aerodynamic moment is given by

   its components along the axes of board coordinate system as follows[12]:

   1 = 1

   1 = 1

   1 = 1 (8)

   Where:

   , , .. Characteristic linear dimensions of the

   missile

   1, 1, 1dimensionless aerodynamic moment

   coefficient

   The aerodynamic moments coefficients are dependent upon the body aerodynamic shape; for the underlying missile, the aerodynamic moments coefficients are given as follows:

   Where

   m is the instantaneous missile mass

   v is the missile velocity

   H is the angular momentum

   The following assumptions are considered during derivation of the equations of missile motion:

   1. The mass of the missile is remains constant during any particular dynamic analysis.

   2. The missile is assumed to be a rigid body which means that any two points on or within the airframe remains fixed w.r.t each other.

   3. The earth is an inertial reference frames Then

   ()

   + 1

   1

   =

   = ( ) + ( ) (11)

   1 = 1

   1

   So, the linearized equations of missile motion are obtained as

   + 1 1 +

   follows

   1 = 1

   1

   1

   F m(V V

   V )

   + 1 1 + (9)

   x1

   x1 z1 y1

   y1 z1

   (12)

   1 = 1

   Where:

   1

   1

   Fy1 m(Vy1 Vx1z1 Vz1x1 )

   mx , y , z

   is the aerodynamic derivative due to fin

   z1

   z1 y1 x1

   x1 y1

   r ,a

   1 1 1

   deflections

   Aslo

   F m(V V

   V )

   F mg

   T X

   wx1 , wy1 , wz1 is the angular velocity components.

   wx1 ,wy1 ,wz1

   x x x

   Fy mg y Ty Y

   (13)

   1 1 1

   mx , y , z

   is the aerodynamic derivative due to angular

   F mg

   T Z

   velocity

   z z z

   From Equations (12) and (13), we can get that:

   Lastly the total force and moments acting on the missile body

   can be illustrated on Fig. (2).

   Tx1

   • Cx1Sq mgx1

     m(Vx1

     Vz1y1

     Vy1z1 )

     Ty1 Cy1Sq mg y1 m(Vy1 Vx1z1 Vz1x1 )

     z1 x1 x1 Vz1 y1 x1 x1 y1

     T C Sq mg m( V V ) (14)

     but

     Tx1 T ,Ty1 Tz1 0

     gx1 g sin

     gy1 g cos cos

     gz1 g sin cos

     From the previous equation it can be considered that:

     Fig. 2. Acting forces and moments on the missile fuselage

     Vx1

     Tx Cx Sq g m m

     sin

     Vy1z1

     Vz1y1

     (15)

     • C Sq

       V y g

       V

       V

       D. Equation of the missile motion

       The equations describing the missile motion in the space can

       y1 m

       • C Sq

       cos cos

       z1 x1

       x1 z1

       be derived easily from Newtons second law of motion which

       z

       Vz1

       m

   • g cos sin Vx1y1 Vy1x1

   states that the summation of all external forces acting on a body

   must be equal to the time rate of change of its momentum and the summation of all external moments acting on a body must be equal to the time rate of change of its momentum[13]:

   = ()

   The linear displacement in GCS (kinematic translation of missile c.g. motion is as follows:

   = 1 () ( ) 1

   + ( + ) 1

   = 1 () () 1 () 1

   = 1 () + ( + ) 1

   =

   (10)

   ( ) 1

   (16)

   The aerodynamic force and moments coefficients are mainly

   function of angle of attack and side slip angle . Thus it is

   i j k

   necessary to show the dependence of theses angles upon the

   H x1

   y1

   z1

   velocity components as follows:

   H H H

   x y z

   (23)

   V

   tan1 ( y1 )

   Vx1

   sin1 ( Vz1 )

   V

   (17)

   Ixy Ixz Iyz 0,

   so

   The guided missile is assumed to be a massless body that is firmly attached to the board coordinate system. The position of the BCS w.r.t the ground coordinate system is determined by means of (,,), while the missile rotation around its c.g. corresponds to the rotation of the BCS around its origin, which is positioned at the missile c.g. Thus the angular velocity vector

   H [(Izz Iyy )z1y1]i[(Ixx Izz )z1 x1] j[(Iyy Ixx )x1y1]k

   Thus

   M ax1 I xx x1 (Izz I yy )z1y1

   M ay1 I yy y1 (I xx Izz )z1x1

   (24)

   ( of missile rotation w.r.t GCS is given by its components in

   the directions of BCS axes as = [1 11] which can be

   related to Eulers angles through the direction cosine matrix as

   follows

   M az1 Izz z1 (I yy I xx )x1y1

   . From the above discussion it is found that the missile equations of motion are function of not only the missile

   x1

   sin

   0 1

   (18)

   physical properties but also the aerodynamic forces and

   y1 cos cos

   cos sin

   sin

   cos

   0

   0

   moment coefficients. Those coefficients can be calculated using

   The Missile DATCOM software. The results of the

   z1

   First element yields and by multiplying 2nd element by sin and 3rd element by cos and add to yield , then multiply 2nd element by cos and 3rd by sin and subtract to yield the above matrix equation can lead to the kinematic equations of missile rotation around its c.g. as fllows

   (y1 cos z1 sin ) cos ec

   aerodynamic coefficients and the missile trajectory with different target engagements will be shown in the next section.

3. SIMULATION AND RESULTS

   The Missile DATCOM accepts an input text file including the missile physical parameters and geometric characteristics (the shape and type of the wing/tail aerofoil, the position of the centre of gravity, the location of wing or tail, etc.), the flight conditions, the mass properties, and so on[1]. The main

   y1 sin z1 cos

   x1 (y1 cos z1 sin ) tan

   (19)

   parameters of the physical configuration of the missile and the

   flight conditions used in this study are included in the input file of Missile DATCOM and are given in Table II. The missiles body was divided into stations as to supply detailed information

   From Newtons second law of rotary motion which states

   that the summation of all external moments acting on a body must be equal to the time rate of change of its moment of momentum i.e.

   of the physical geometry of the missile.

   TABLE II. MISSILE DATCOM INPUT CONFIGURATION FILE

   PARAMETERS

   M d( H)

   dt

   Parameter	Value [unit]
   Reference area (Sref)	0.0.01736[m2]
   Longitudinal reference length (c)	1.51 [m]
   Lateral reference length (b)	0.152 [m]
   Airspeed (VT)	0.1, 0.2, 0.3, 0.4, 0.5,0.6,0.7,0.8 [M]
   Angle of attack ()	-4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8 ,9,10,11,12,13,14,15[°]
   Angle of sideslip ()	0 [°]
   Elevator Fin deflection ()	2,-6 [°]
   Rudder Fin deflection ()	6,-6 [°]

   (20)

   The moments acting upon the missile during its flight are given from the following vector equation

   M dH H H H dt

   (21)

   For two plane of symmetry, like the underlying missile

   Hx Ixx x1

   And

   Hy I yy y1

   Hz Izz z1

   (22)

   The output file of Missile DATCOM gives the total aerodynamic forces and moments coefficients as multidimensional pages. Each page corresponds to a selected value of the Mach number, the sideslip angle, the altitude and the rear fin deflection.

   1. Lift Force Aerodynamic Coefficient (1)

      For efficient flight at subsonic speeds, it is necessary to use

      wings that have as large a span as possible to reduce the drag. Additionally, it is advantageous to use aero foil section that are slightly curved. The missile is symmetrical in pitch and yaw planes. The coefficient of lifting force is a function of the angle of attack (), rear fin deflection () according to the following relations:

      y1 y1o y1 y1

      C C C *sin C *

      (25)

      Fig. 5. Side Force Coefficient

      Where 10 is the coefficient of lifting force when ==0.

      For the underlying system it is equal to zero because the missile

      is symmetrical about its longitudinal axis. From Missile DAT COM output file we can get the lift force coefficient at every rear fin deflection as shown Fig.3.

   2. Drag Force Coefficient (1 )

      Trajectory calculations and autopilot design require the

      missile configuration and a method by which the drag coefficient can be efficiently calculated by responsible accuracy. Essentially, there are two separate contributions to the drag force on the missile through its flight, one due to the

      pressure and viscous forces on the missile (0) while the other is associated with the production of lift ( ) such that:

      1 = 0 + (26)

      So, the total drag force coefficient can be also yielded from

      DAT COM as follows in Fig.4.

   3. Side Force Coefficient ( )

      The side force originates from the deflection of the rudder

      fins as shown in Fig.5.

      Fig. 3. Lift Force Aerodynamic Coefficient (a,b) due to elevator deflection, (c,d) due to rudder deflection

      Fig. 4. Drag Force Coefficient

   4. Pitching Moment Coefficient (1)

      The moment arises from non-zero angle of attack () and

      when the aerodynamic force is related to the missile c.g. it causes the missile to rotate in the pitch plane and consequently it is called pitch moment and can be expressed as follows:

      1 = 1 (27)

      The DATCOM result for pitching moment coefficient due to

      elevator fin deflections is shown in Fig.6.

   5. Yawing Moment Coefficient (1)

   Like pitching moment; yawing moment is also arisen from

   the non-zero side slip angle () and consequently the side force

   causes the missile to rotate around its yaw axis (1) and can be expressed as follows:

   1 = 1 (28)

   The DATCOM result for yawing moment coefficient due to

   rudder fin deflections is shown in Fig.7.

   Fig. 6. Pitching Moment Coefficient

   Fig. 7. Yawing moment Coefficient

   These coefficients are then implemented in the full mathematical model of the missile as well as the equations of motion to get the missile trajectory.

   1. Flow Chart of the missile Simulation model

      The set of differential equations governing the missile spatial motion is solved using Rung Kutta-4 method. First the initial target and missile parameters are manually entered to the program and depending on the pitch and yaw loops illustrated in Fig.8 the error is filtered, compensated and the steering commands are generating causing the missile rear fins to deflect

      with ( ) in pitch plane and ( ) in yaw plane.

      According to the missile fin deflections (), Mach number (M), angle of attack () and side slip angle () the corresponding force and moment coefficients are obtained using lookup table yielding the new missile position and new angular position in yaw plane () and pitch plane (), these calculations are carried out till the missile interception yielding the miss distance in pitch and yaw planes or missile ground hit occurs.

   2. Validation of the Missile Model

   The missile model is validated by comparing its output flight path trajectory with real recorded shootings results. Fig.9 and Fig.10 show the radial velocities and the trajectories between real and simulated missile data. Fig.11 (a,b,c,d) show the engagement of the missile with a target in different scenarios.

   Fig. 8. Flow Chart of the missile simulation model

   Fig. 9. Comparison between simulated and real missile data

   Fig .10. Comparison between simulated and real missile data

   Fig .11 (a) Yaw plane trajectory with target velocity 13m/sec

   Fig .11 (b) Yaw plane trajectory with target velocity -5m/sec

   Fig .11 (c) Yaw plane trajectory with stationary target

   Fig .11 (d) Pitch plane trajectory

   From the previous figures it is shown that the missile is able to hit the target at different velocities and due to Bang-to- Bang movements of the rear control fin the missile goes up and down in pitch plane until hitting the target.

4. CONCLUSIONS AND FUTURE WORK

This paper has presented the total forces and moments that affecting on the missile body and also the coordinate reference frames are then defined. Taking into consideration the forces and moments affecting on the missile fuselage, the equations of motion that governing the missile spatial motion are derived. A special software technique (Missile DATCOM) is adopted to calculate the coefficients of forces and moments as a function of Mach number, angle of attack, side slip angle and contrl fin deflections. The complete 6 DOF is created based on the

DATCOM results such that the simulated results show a satisfactory behaviour when validated with real shooting results. However the simulated missile success to attack the missile at different scenarios but it fails to hit it with thrust degradation so more investigation on the controller loop will be carried out in the future work.

REFERENCES

[1] USAF, "Missile DATCOM " USAF, ed., 1998.

[2] Islam H. Elnady, Ahmed M. Kamel, Yehia Z. Elhalwag, Modeling and Simulation of an Aerial Gliding Body in Free-Fall, the International Journal of Engineering Research and Technology vol. 7, no. 8, 2018.

[3] G. El-Sheikh, Guidance System and Theory, Cairo, 2010.

[4] L. S. a. Callaborators, Principles of Guided Missile Design, 1955.

[5] A.N.Ouda, Adaptive Autopilot Design for Flying Vehicle, the International journal of dynamic and control (IJDAC), 2017.

[6] A. M. Kamal, Modeling, Analysis and Identification of Airplane Flight Dynamics, Military Technical College, Cairo, 2014.

[7] J. D. Anderson, FUNDAMENTALS OF AERODYNAMICS, 1991.

[8] A. H. A. Elsayed Ahmed, A. N. Ouda, H. Eldin Hussein Ahmed, H. Mohamed Abd-Elkader, Modeling of a Small Unmanned Aerial Vehicle, the International Journal of advanced robotic and automation vol. 4, no. 126, 2015.

[9] A. E. E.-H. Shady A. Kader, A. N. Ouda, Quadcopter System Modeling and Autopilot Synthesis, the International Journal of Engineering Research and Technology, vol. 3, no. 11, 2014.

[10] P. Zarcan, Tactical and Strategic Missile Guidance, Progress in Astronautics and Aeronautic. AIAA, 1994.

[11] R.M.Goreckt, A baseline 6 degree of freedom mathematical model of generic missile, DSTO Systems Sciences Laboratory Australia 2003.

[12] N. Jack N, Missile Aerodynamics, Mountain View, USA, 1988.

[13] G. M. Siouris, Missile Guidance and Control Systems, Ohio,USA, 2003.