Application of Pontryagin’s Minimum Principle in Optimum Time of Missile Manoeuvring CAUCHY - JOURNAL OF PURE AND APPLIED MATHEMATICS Volume 4(3) (2016), Pages 107-111 Submitted: June, 28 2016 Reviewed: October, 21 2016 Accepted: November, 29 2016 DOI: http://dx.doi.org/10.18860/ca.v4i3.3534 Application of Pontryagin’s Minimum Principle in Optimum Time of Missile Manoeuvring Sari Cahyaningtias1 ,Subchan2 1Adi Buana University 2Kalimantan Institute of Technology Email: scahyaningtias@gmail.com, s.subchan@gmail.com ABSTRACT Missile is a guided weapon and designed to protect outermost island from a thread of other country. It, commonly, is used as self defense. This research presented surface-to-surface missile in final dive manoeuvre for fixed target. Furthermore, it was proposed manoeuvring based on unmanned aerial vehicle (UAV), autopilot system, which needs accuration and minimum both time and thrust of missile while attacking object. This paper introduced pontryagin’s Minimum Principle, which is useable to solve the problem. The numerical solution showed that trajectory of the missile is split it up in three sub-intervals; flight, climbing, and diving. The numerical simulation showed that the missile must climb in order to satisfy the final dive condition and the optimum time of a missile depend on initial condition of the altitude and the terminal velocity. Keywords: missile, optimal control, pontryagin’s minimum principle, thrust, angle of attack. INTRODUCTION Missile is a military rocket which has automatic control system to target an object and suit the direction. Missile is one of the unmanned aeral vehicle (UAV) which uses both in military and civilian needed such as controlling forest from illegal logging and sea from illegal fishing. It clearly is known that UAV can be handled using autopilot system to follow the references trajectory which could minimized both the risk and cost [1] Trajectory tracking of surface-to-surface missile is divided in three sub-interval, namely: minimum altitude flight, climbing, and diving [2] [3]. Restrictiveness of fuel supply in maoeuvring can be handled by minimizing time manoeuvring to attack the target. Mathematics modelling of mass centre missile follow several assumption: the earth is spherical and non-rotation and ignoring lateral movement of the missile. Some boundary problems also were considered in this research such as massa of missile remains constant π‘‘π‘š 𝑑𝑑 β‰… 0; there was no any obstacle while manoeuvring; The thrust and angle of attack are the controller of this problem which should be attacked fixed target. Trajectory problem of missile has been studied by many researches both using classical control theory and modified optimal control. In 2013, Wang and Dong using hybrid algorithm to intercept trajectory optimization for multi-stage air defense missile which the missile always flies in the plane formed by the earth centre, launch point, and the predicted intercept point [4]. Maopeng at. al. In 2014, integrated missile guidance and control law using adaptive fuzzy sliding mode control. They also conducted stability of close loop system [5]. http://dx.doi.org/10.18860/ca.v4i3.3534 mailto:scahyaningtias@gmail.com mailto:s.subchan@gmail.com Application of Pontryagin’s Minimum Principle in Optimum Time of Missile Manoeuvring Sari Cahyaningtias 108 This paper presents computational result of the optimal trajectory of missile which minimized the time manoeuvring along the optimal trajectory. Furthermore, Vertical dive manoeuvre was chosen as final maneuvers by considering effect of its attacking. It would have shot to pieces the target permanently. Optimal control is applied by forcing the thrust and angle of attack to get optimum time using Pontryagin’s Minimum Principle. METHODS Purpose of this research is to attack motionless target precisely and to minimize time manoeuvring. Here was given the definition of missile axes and angles. Figure 1 Definition of missile axes and angles Based on Figure 1, The dynamic equations of a point mass missile, which may be written as [6]: 𝑑𝑉 𝑑𝑑 = 1 π‘š [𝑇 cos 𝛼 βˆ’ 𝐷] βˆ’ 𝑔 sin 𝛾 (1) [7] 𝑑𝛾 𝑑𝑑 = 1 π‘šπ‘‰ [𝐿 + 𝑇 sin 𝛼] βˆ’ 𝑔 𝑉 cos 𝛾 (2) 𝑑π‘₯ 𝑑𝑑 = 𝑉 cos 𝛾 (3) π‘‘β„Ž 𝑑𝑑 = 𝑉 sin 𝛾 (4) There are four state variables, namely: speed 𝑉, 𝛾 flight path angle, horizontal position π‘₯, and altitude β„Ž. The control variables are Thrust 𝑇 and angle of attack 𝛼. Axial aerodynamics forces and normal aerodynamics forces are function of altitude, speed, and angle of attack, which are defined as: 𝐷(β„Ž, 𝑉, 𝛼) = 1 2 πΆπ·πœŒπ‘‰ 2𝑆 𝐢𝐷 = 𝐴1𝛼 2 + 𝐴2𝛼 + 𝐴3 𝐿(β„Ž, 𝑉, 𝛼) = 1 2 𝐢𝐿 πœŒπ‘‰ 2𝑆 𝐢𝐿 = 𝐡1𝛼 + 𝐡2 where 𝜌, air density, is given by 𝜌 = 𝐢1β„Ž 2 + 𝐢2β„Ž + 𝐢3 whereas S reference area of missile, π‘š Massa and 𝑔 gravity. Table 1. Physical Modelling Parameter Parameter Nilai Unit M 1005 Kg G 9.81 m/sΒ² S 0.3376 mΒ² A1 -1.9431 A2 -0.1499 A3 0.2359 B1 21.9 B2 0 C1 3.312 10-9 Kg m-5 C2 -1.142 10-4 Kg m-4 C3 1.224 Kg m-3 Table 2. Boundary Condition and Constraints Parameter Nilai Unit V(0) 210 m/s2 𝛾(0) 0 Rad π‘₯(0) 0 m β„Ž(0) 30, 100, 200 m 𝑉(𝑑𝑓 ) 290, 310 m/s 2 Application of Pontryagin’s Minimum Principle in Optimum Time of Missile Manoeuvring Sari Cahyaningtias 2 𝛾(𝑑𝑓 ) βˆ’ πœ‹ 2 rad π‘₯(𝑑𝑓 ) 10000 meter π‘‡π‘šπ‘Žπ‘˜π‘  6000 N π‘‡π‘šπ‘–π‘› 1000 N Boundary of control π‘‡π‘šπ‘–π‘› ≀ 𝑇 ≀ π‘‡π‘šπ‘Žπ‘˜π‘  The problem of this research is how to get the minimum time of missile to attack the fixed target in vertical dive manoeuvre. Performance index of this problem is function of time which is given by: J = ∫ dt tf t0 (5) where 𝑑 time of manoeuvring, 𝑑0 ≀ 𝑑 ≀ 𝑑𝑓 in which 𝑑0 initial time and 𝑑𝑓 final time. RESULTS AND DISCUSSION Mathematical analysis showed that some parts of those solution couldn’t be gained using analytic method. That is why, computational analysis was applied in that case. This paper used optimal control Toolbox, namely Miser3 which is developed by Jennings [8] The trajectory of surface-to-surface missile with vertical dive manoeuvre is split of three sub- interval: 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 500 1000 1500 2000 2500 3000 WAKTU (Normalisasi) K E T IN G G IA N ( M e te r) V(tf)=310 V(tf)=290 Figure 4 Performance of angle of attack missile toward time Figure 5 Performance of thrust of missile toward time Figure 2 Performance of missile’s altitude toward normalization time 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000 6500 WAKTU (Normalisasi) G A Y A D O R O N G ( N e w to n ) V(tf)=310 V(tf)=290 Figure 3 Performance of missile’s velocity toward nirmalization time 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 -0.08 -0.06 -0.04 -0.02 0 0.02 0.04 0.06 0.08 WAKTU (Normalisasi) S U D U T S E R A N G ( ra d ia n ) V(tf)=310 V(tf)=290 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 200 220 240 260 280 300 320 WAKTU (Normalisasi) K E C E P A T A N ( m /d e ti k ) V(tf)=310 V(tf)=290 109 Application of Pontryagin’s Minimum Principle in Optimum Time of Missile Manoeuvring Sari Cahyaningtias 110 1. Minimum altitude flight The missile launches at current altitude and therefore the altitude constraints are active directly at the start of the manoeuvre. In this case the speed of the missile remains constant on the minimum value hmin (see Figures 2) until the missile must start climbing while the thrust is on the maximum value. The flight time of first arc depends to the final-speed. Figure 5 shows that the thurst of missile remains at maximum, 6000 N, while the speed slightly higher until the next sub-interval in Figure 3. 2. Climbing The Speed of missile after flight in minimum altitude decreased gradually until reached the optimum altitude to start manoeuvre which is vertical dive manoeuvre. The thrust of missile remains maximum until time-0.6 (normalization time), furthermore there suddenly switched into lowest boundary, 1000 N. It showed that missile was turning down and starting dive manoeuvre. 3. Diving The missile must gain the power to reach the target therefore the speed increase rapidly since the initial diving speed is lower than the final speed. The normal acceleration is saturated on the minimum value for this arc. Whereas the thrust remains constant around 1000 N. The angle of attack convergenly is moving at –pi/40 radian. Here presents computational solution of minimizing time problem depends on initial altitude and final speed. Table 3. Optimum time by computation Initial Altitude (meter) Final Speed (m/s) Objective Function (s) 30 290 48.0781385 310 44.8815557 100 290 47.9329181 310 44.8199641 200 290 47.7656314 310 44.6765635 CONCLUSION In this paper, analysis of surface-to-surface missile which considered objective and dynamical system showed that trajectory tracking of this missile was split up in three sub- intervals, namely flying, climbing, and diving in which had different performances of speed. Those paths were result by controlling thrust and angle of attack. First path, the missile launched at h=30m and remained stable in this altitude. Second part, velocity of missile decreased dramatically while climbing to reach optimum altitude before manoeuvring. Furthermore the last path, dive manoeuvring, the speed of missile got faster rapidly to final speed. It would be happened while missile was at 30 meters above the target with angle of attack’s -90 deg. Travel time of missile was influenced by starting position while it was shoot and final speed of missile itself. The higher launch site and the bigger final speed of missile affected the shorter its time travel. Another variable that should be considered in missile manoeuvring is energy. It can be minimize in objective function. This problem can be conducted in future reserach. Application of Pontryagin’s Minimum Principle in Optimum Time of Missile Manoeuvring Sari Cahyaningtias 111 REFERENCES [1] S. Subchan, "Pythagorean Hodograph Path Planning for Tracking Airborne Contaminant using Sensor Swarn," in International Instrumentation and Measurement Technology Conference Victoria, Canada, 2008. [2] S. Subchan, "Trajectory happing of Surface-to-Surface Missile with Terminal Impact Angle Constraint," Makara Teknologi Indonesia University, pp. 65-70, 2007. [3] S. Subchan and R. Zbikowski, Computational Optimal Control: Tools and Practice, UK: John Wiley and Sons Ltd, 2009. [4] F. B. Wang and C. H. Dong, "Fast Intercept Trajectory Optimization for Multi-Stage Air Defense Missile using Hybrid Algorthm," Elsivier Procedia engineering, pp. 447-456, 2013. [5] e. a. Maoping, "Backstepping Design of Missile Guidance and Control Based on Adaptive Fuzzy Sliding Mode Control," Chinese Journal of Aeronautics, pp. 634-642, 2014. [6] G. Siouris, Missile Guidance and Control Systems, USA: Springer, 2003. [7] D. S. Naidu, Optimal Control Systems, USA: CRC Press LLC, 2002. [8] L. S. Jennings, M. E. Fisher, K. L. Teo and C. J. Goh, Miser3 Optimal Control Software, Australia: The University Of Western Australia, 2002.