Research Article | Open Access

F. A. Abd El-Salam, S. E. Abd El-Bar, "Computation of the Different Errors in the Ballistic Missiles Range", *International Scholarly Research Notices*, vol. 2011, Article ID 349737, 16 pages, 2011. https://doi.org/10.5402/2011/349737

# Computation of the Different Errors in the Ballistic Missiles Range

**Academic Editor:**Z. Huang

#### Abstract

The ranges of the ballistic missile trajectories are very sensitive to any kind of errors. Most of the missile trajectory is a part of an elliptical orbit. In this work, the missile problem is stated. The variations in the orbital elements are derived using Lagrange planetary equations. Explicit expressions for the errors in the missile range due to the in-orbit plane changes are derived. Explicit expressions for the errors in the missile range due to the out-of-orbit plane changes are derived when the burnout point is assumed on the equator.

#### 1. Introduction

The fundamental problem of astrodynamics is the orbit determination and orbit correction. For a spacecraft moving under the influence of gravitational field of Earth in free space (no air drag), the trajectory is an ellipse with the center of Earth lying at one of the foci of the ellipse. This constitutes a standard two-body-central-force problem, which has been treated, in detail, in many standard textbooks [1, 2]. The trick is to first reduce the problem to two dimensions by showing that the trajectory always lies in a plane perpendicular to the angular momentum vector. Then, the problem is set up in plane-polar coordinates. Angular momentum is conserved, and the problem, effectively, reduces to one-dimensional problem [3].

A ballistic missile is a missile that follows a suborbital ballistic flight path with the ballistic missile objective of delivering one or more warheads (often nuclear) to a predetermined target. The missile is only guided during the relatively brief initial powered phase of flight, and its course is subsequently governed by the laws of orbital mechanics and ballistics. To date, ballistic missiles have been propelled during powered flight by chemical rocket engines of various types. Therefore, the ballistic missile trajectory consists of three parts; see Figure 1(a). (1) *The powered flight portion*, sometimes called boost phase, takes usually from 3 to 5 minutes (shorter for a solid rocket than for a liquid-propellant rocket); the altitude of the missile at the end of this phase is typically 150 to 400βkm depending on the trajectory chosen, and the typical burnout speed is 7βkm/s. (2)* The free-flight portion*, or the midcourse phase which constitutes most of the flight time, takes approximately 25 minutes. It is a part of an elliptic orbit with a vertical major axis; the apogee is at an altitude of approximately 1,200βkm; the semimajor axis is between 3,186βkm and 6,372βkm; the projection of the orbit on the Earth's surface is close to a great circle, slightly displaced due to Earth rotation during the time of flight; the missile may release several independent warheads, and penetration aids such as metallic-coated balloons, aluminum chaff, and full-scale warhead decoys. (3) *The reentry phase* during which energy is dissipated as a result of friction with the atmosphere, starts at an ill-defined point at an altitude of 100βkm. It takes about 2 minutes to impact at a speed of up to 4βkm/s (for early ICBMs less than 1βkm/s).

**(a)**

**(b)**

Ballistic missiles can be launched from fixed sites or mobile launchers, including vehicles (transporter erector launchers, TELs), aircraft, ships, and submarines. The powered flight portion can last from a few tens of seconds to several minutes and can consist of multiple rocket stages. When in space and no more thrust is provided, the missile enters free flight. In order to cover large distances, ballistic missiles are usually launched into a high suborbital spaceflight; for intercontinental missiles, the highest altitude (apogee) reached during free flight is about 1200βkm. The reentry stage begins at an altitude where atmospheric drag plays a significant role in missile trajectory and lasts until missile impact.

#### 2. Types of Ballistic Missiles

Ballistic missiles are categorized according to their range, the maximum distance measured along the surface of the Earth's ellipsoid from the point of launch of a ballistic missile to the point of impact of the last element of its payload. Various schemes are used by different countries to categorize the ranges of ballistic missiles as follows: tactical ballistic missile: range between about 150βkm and 300βkm, battlefield range ballistic missile (BRBM), range less than 200βkm, theatre ballistic missile (TBM): range between 300βkm and 3500βkm, short-range ballistic missile (SRBM): range 1000βkm or less, medium-range ballistic missile (MRBM): range between 1000βkm and 3500βkm, intermediate-range ballistic missile (IRBM) or long-range ballistic missile (LRBM): range between 3500βkm and 5500βkm, intercontinental ballistic missile (ICBM): range greater than 5500βkm, and submarine-launched ballistic missile (SLBM): launched from ballistic missile submarines (SSBNs), and all current designs have intercontinental range.

Short and medium-range missiles are often collectively referred to as theater or tactical ballistic missiles (TBMs). Long- and medium-range ballistic missiles are generally designed to deliver nuclear weapons, because their payload is too limited for conventional explosives to be efficient (though the U.S. may be evaluating the idea of a conventionally armed ICBM for near-instant global air strike capability despite the high costs). The flight phases are like those for ICBMs except with no exoatmospheric phase for missiles with ranges less than about 350βkm.

Sometimes, the designers of the ballistic missiles need to perform maneuvers in flight or make unexpected changes in direction and range; this type is known as a quasi-ballistic missile or a semi ballistic missile. At a lower trajectory than a ballistic missile, a quasi-ballistic missile can maintain higher speed, thus allowing its target less time to react to the attack, at the cost of reduced range.

#### 3. Literature Survey

McFarland [4] treated ballistic missile problem by modeling spherical Earth, Earth rotation, and addition of atmospheric drag using the state transition matrix. Forden [5] described the integration of the three degrees of freedom equations of motion, and approximations are made to the aerodynamic for simulating ballistic missiles. Bao and Murray [6] improved the ground range calculation of a ballistic missile trajectory on a nonrotating oblate Earth. Isaacson and Vaughan [7] described a method of estimating and predicting ballistic missile trajectories using a Kalman filter over a spherical, nonrotating Earth. They determined uncertainties in the missile launch point and missile position during flight. Harlin and Cicci [8] developed a method for the determination of the trajectory of a ballistic missile over a rotating, spherical Earth given only the launch position and impact point. The iterative solution presented uses a state transition matrix to correct the initial conditions of the ballistic missile state vector based upon deviations from a desired set of final conditions. AkgΓΌl and Karasoy [9] developed a trajectory prediction program to predict the full trajectory of a tactical ballistic missile. Vinh et al. [10] obtained a minimum-fuel interception of a satellite, or a ballistic missile, in elliptic trajectory in a Newtonian central force field, via Lawdenβs theory of primer vector. Kamal [11] developed an algorithm includes detection of cross-range error using Lambert scheme in free space in the absence of atmospheric drags. Bhowmik and Sadhukhan [12] investigated the advantages and performance of extended Kalman filter for the estimation of nonlinear system, where linearization takes place about a trajectory that was continually updated with the state estimates resulting from the measurement. They took tactile ballistic missile reentry problem as a nonlinear system model and extended Kalman filter technique is used to estimate the positions and velocities at the and direction at different values of ballistic coefficients. Kamal [13] presented an innovative adaptive scheme which was called βthe multistage Lambert schemeβ. Liu and Chen [14] presented a novel tracking algorithm by integrating input estimation and modified probabilistic data association filter to identify warhead among objects separation from the reentry vehicle in a clear environment.

#### 4. Statement of the Problem

Our ICBM problem concerns with the determination of the free-flight range angle taking into account the perturbation in the orbital elements. Let us define the dimensionless parameter as , where is the magnitude of the position vector of the missile relative to the Earth, is the missile speed at any point in its orbit, is the corresponding missile circular speed at this point, and , where is the gravitational constant, is the mass of the Earth, and is the missile mass. From the orbital mechanics of two body and the symmetry of the free-flight portion shown in Figures 1(a) and 1(b), we have [15] where is the flight path angle, the true anomaly, and the dimensionless parameter at the burnout point.

Equation (4.1) gives the free-flight range by which the free-flight range angle can be computed for any given combination of burnout conditions , , and . The problem can now be specified as βgiven a particular launch point and target, it is required to calculate , and knowing , , can be calculatedβ.

Equation (4.2) gives the flight-path angle, which provides two trajectories to the missile. The trajectory corresponding to the larger value of is called the high trajectory and to the smaller value is the low trajectory. The nature of the trajectory, high or low, depends primarily on the value of . This is obvious from (4.2), where if, is always less than 180Β°; otherwise, the right side of (4.2) exceeds 1, and both high and low trajectories are possible,, one trajectory is circular, and for , both high and low trajectories are possible ( for low), while high trajectory only is possible for , and low trajectory skims Earth ( for high),, (4.2) yields one positive and one negative value for regardless of range. The low trajectory, corresponding to the negative value, is not practical, since it would penetrate the Earth.

When the right hand side of (4.2) equals unity, we obtain a single trajectory called the maximum range trajectory, , and the flight path angle is .

#### 5. In-Orbit-Plane Changes

The variations in the parameters due to changes in the orbital elements take place in the plane of the orbit. Therefore, in what follows, we will compute the errors in due to changes in the mentioned orbital elements. To do this, we need first the following partial derivatives:

##### 5.1. Error in due to the Change in the Semimajor Axis

We can write the change in the free flight range angle due to the change in the semimajor axis as follows: The required derivatives are given by The involved products are Using the Lagrange planetary equations, we computed taking into account the oblate model of the Earth (retaining the zonal harmonics up to ). The integration is performed between and . Substitution of the obtained expression into (5.2) yields where nonvanishing coefficients are given by

##### 5.2. Error in due to the Change in the Eccentricity

We can write the change in the free flight range angle due to the change in the eccentricity as follows: The required derivatives are given by The bracket vanishes when , which is impossible due to the fact that the eccentricity is a real value.

The involved products are Using the Lagrange planetary equations, we computed (retaining the zonal harmonics up to ). The integration is performed between and . Substitution of the obtained expression into (5.7) yields where where nonvanishing coefficients are given by

#### 6. Out-of-Orbit Plane Changes

All out-of-orbit plane changes, for example, , will cause a cross-range errors .

##### 6.1. Error in due to the Change in the Ascending Node

For the sake of the simplicity, let us take the burnout point on the equator. For some reason, it was displaced by an amount, . This displacement could be interpreted as a change in the longitude of the ascending node if the rest of the orbital elements were kept fixed. Due to this change, a cross-range error,ββ, at impact occurs, the value of which is obtained, in a similar manner as the previous subsection, by applying the law of cosines. Hence, Since , are small angles, then we have Using the Lagrange planetary equations, we computed (retaining the zonal harmonics up to ). The integration is performed between and . Substitution of the obtained expression into (6.2) yields where where nonvanishing coefficients are given by

##### 6.2. Error in due to the Change in the Inclination

Again, assume that the burnout point is on the equator and that the actual launch azimuth differs from the intended value by an amount . This amount could be interpreted as a change in the orbital inclination if all other orbital elements were kept fixed. Due to this change, a cross-range error,ββ, at impact occurs, the value of which is obtained, in a similar manner as the previous subsection, by applying the law of cosines. Hence, Since are small angles, then we have Using the Lagrange planetary equations, we computed (retaining the zonal harmonics up to ). The integration is performed between and . Substitution of the obtained expression into (6.7) yields where where nonvanishing coefficients are given by

#### 7. Conclusions and Future Work

Due to the high sensitivity of the ballistic missile range to the different kinds of errors, we computed the explicit expressions for the errors in the missile range due to the in-orbit plane changes. We derived explicit expressions for the errors in the missile range due to the out-of-orbit plane changes when the burnout point is assumed on the equator. In a forthcoming work, we aim to generalize this situation. Also, we aim to do the corresponding algorithms and give numerical examples.

#### Acknowledgments

The authors are deeply indebted to the Professor Dr. M. K. Ahmed, the professor of space dynamics at Cairo University, Faculty of Science, Department of Astronomy, for his valuable discussions and critical comments and ideas that help us to finalize this work. This research work was supported by a Grant no. (627) from the deanship of the scientific research at Taibah university, Al-Madinah Al-Munawwarah, Saudi Arabia.

#### References

- H. Goldstein,
*Classical Mechanics*, Addison-Wesley, Reading, Mass, USA, 2nd edition, 1981. View at: Zentralblatt MATH - J. B. Marion,
*Classical Dynamics of Particles and Systems*, Academic Press, New York, NY, USA, 2nd edition, 1970. - J. T. Wu, βOrbit determination by solving for gravity parameters with multiple arc data,β
*Journal of Guidance, Control, and Dynamics*, vol. 15, no. 2, pp. 304β313, 1992. View at: Publisher Site | Google Scholar | Zentralblatt MATH - J. S. McFarland, βModeling the ballistic missile problem with the state transition matrix: an analysis of trajectories including a rotating earth and atmospheric drag,β
*Semester Report*, Virginia Tech, Blacksburg, Va, USA, 2004. View at: Google Scholar - G. Forden, βGUI missile flyout: a general program for simulating ballistic missiles,β
*Science and Global Security*, vol. 15, no. 2, pp. 133β146, 2007. View at: Publisher Site | Google Scholar - U. N. Bao and E. D. Murray, βComputation of effective ground range using an oblate earth model,β
*The Journal of the Astronautically Sciences*, vol. 51, no. 3, pp. 291β305, 2003. View at: Google Scholar - J. A. Isaacson and D. R. Vaughan,
*Estimation and Prediction of Ballistic Missile Trajectories*, RAND Corporation, Santa Monica, Calif, USA, 1996. - W. J. Harlin and D. A. Cicci, βBallistic missile trajectory prediction using a state transition matrix,β
*Applied Mathematics and Computation*, vol. 188, pp. 1832β1847, 2007. View at: Publisher Site | Google Scholar - A. AkgΓΌl and S. Karasoy, βDevelopment of a tactical ballistic missile trajectory prediction tool,β
*Istanbul University Journal of Electrical & Electronic Engineering*, vol. 5, no. 2, pp. 1463β1467, 2005. View at: Google Scholar - X. Vinh, T. Kabamba, and T. Takehira, βOptimal interception of a maneuvering long-range missile,β
*Acta Astronautica*, vol. 48, no. 1, pp. 1β19, 2001. View at: Publisher Site | Google Scholar - S. A. Kamal, βCross range error in the lambert scheme,β in
*Proceedings of the 10th National Aeronautical Conference*, S. R. Sheikh, Ed., pp. 255β263, College of Aeronautical Engineering, PAF Academy, Risalpur, NWFP, Pakistan, April 2006. View at: Google Scholar - S. Bhowmik and C. Sadhukhan, βApplication of extended kalman filter to tactical ballistic missile re-entry problem,β 2007. View at: Google Scholar
- S. A. Kamal, βThe multi-stage-lambert scheme for steering a satellite-launch vehicle (SLV),β in
*Proceedings of the 12th IEEE International Multitopic Conference*, M. K. Anisx, M. K. Khan, and S. J. H. Zaidi, Eds., pp. 294β300, Bahria University, Karachi, Pakistan, December 2008. View at: Google Scholar - C. Y. Liu and C. T. Chen, βTracking the warhead among objects separation from the reentry vehicle in a clear environment,β
*Defence Science Journal*, vol. 59, no. 2, pp. 113β125, 2009. View at: Google Scholar - R. R. Bate, D. D. Mueller, and J. E. White,
*Fundamentals of Astrodynamics*, Dover, New York, NY, USA, 1971.

#### Copyright

Copyright Β© 2011 F. A. Abd El-Salam and S. E. Abd El-Bar. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.