Paper received: 2011-12-20, paper accepted: 2012-04-03 © 2012 Journal of Mechanical Engineering. All rights reserved. Experimental and Simulation Testing of Flight Spin Stability for Small Caliber Cannon Projectile Momcilo Milinovic1* - Damir Jerkovic2 - Olivera Jeremic1 - Mitar Kovac2 1 University of Belgrade, Faculty of Mechanical Engineering, Serbia 2 University of Defense, Military Academy, Serbia The basic aim of this paper is to consider correlations of stability flight criteria, derived as the relations of aerodynamic coefficients and derivatives, on the model of a small caliber cannon spin stabilized projectile. Model of stability criteria calculations are performed by experimentally testing of aerodynamic data in the wind tunnel, and composed with the semi-empirical data, both applied on the flight trajectory stability simulation test. Authors' wind tunnel tests and calculated values of aerodynamic coefficients, as the function of Mach numbers of projectile model, are presented in the simulation flight trajectories stability criteria. The comparative analysis of experimental and calculated aerodynamic coefficients of projectile model is done, and refers to the stability flight criteria. Calculation of projectile aerodynamic Magnus moment derivatives, with other aerodynamic representatives, is used as the critical stability factors testing data vs. flight Mach numbers. Influences of this derivative absence and presence on the model sequence of the flight trajectory are presented for the estimation of angles of attack damping and stability factors. Simulation tests are presented for the supersonic and subsonic integral flight velocities and spin damping data. Research is realized due to the considerations of further projectiles correction possibilities on trajectory and other new applications, vs. existing of unreliable lateral moments. Keywords: aerodynamic coefficients spin stabilized small caliber cannon projectile, gyroscopic stability factor, dynamic stability factor, damping stability coefficients 0 NOMENCLATURE CM pitching moment coefficient, Mv/QSd, [-], CMa derivative of pitching moment coefficient, dC^/da, [-], CMa derivative of pitch damping moment coefficient due to q*, dCM/dq*, [-], Cm derivative of pitch damping moment coefficient due to a , dCM/doc, [-], CL lift (normal) force coefficient, -FJQS, [-], CLa derivative of lift force coefficient, dCL/da, [-], CN yawing moment coefficient, MJQSd, [-], CMpa derivative of Magnus moment coefficient, d2CN/dp*da, [-], CD drag (axial) force coefficient, FJQS, [-], CD0 zero angle drag coefficient, (CD)a=0 , [-], d reference diameter (caliber), Ix axial moment of inertia, [kgm2], Iy transverse moment of inertia, [kgm2], rx relative axial radius of gyration, reversed to the caliber, (IJmd2)1'2, [-], ry relative transverse radius of gyration reversed to the caliber, (Iy/md2)12, [-], Fx, Fy, Fz forces along x, y, z axes, [N], Mx, My, Mz moments about x, y, z axes, [Nm], m mass of projectile, [kg], p spin rate, [s-1], p* reduced spin rate, pd/V, [-], Q dynamic pressure, pV2/2, [Pa], q pitch rate, [s-1], q* reduced pitch rate, qd/V, [-], S reference area, nd2/4, [m2], a angle of attack (pitch), [°], P angle of sideslip (yaw), [°], at total yaw angle, approx. (a2+P2)1/2, [°], £ complex angle of attack, a+ifi, [°], p atmospheric density, [kg/m3], V velocity vector of projectile, [kg/s], u, v, w components of velocity along x, y, z axes, [m/s], H, P, T, M, G Murphy's coefficients, [-], E = (pSd)/(2m) reduction mass expression, [-], Sg gyroscopic stability factor, [-], Sd dynamic stability factor, [-], A1,2 damping stability coefficients, [-]. denotes reduction of coefficients and derivatives C* = CijE. 1 INTRODUCTION In the last twenty years the modern ammunition design, extended precision technology applications of guidance and control on the lower calibers of tactical ammunition. Analyses to be considered for guidance redesigning are especially challenging for the anti-aircraft (AA) cannon ammunition of smaller dimensions because of specific properties of projectiles flight, which have constraining possibilities of guidance technologies applications, [1] to [5]. Unguided AA projectile is, primarily, used for air targets, shooting with ballistic trajectories, but new considerations may also suppose its applications for the ground targets engaged from air and ground platforms in the so called close support operations. Ballistic trajectory of AA gun projectile have significant changes of supersonic transonic and subsonic velocities during flight, from the very high initial to the much lower terminal values in the target impact point. Stable free ballistic flight of projectiles and stability criteria, determines attitudes of projectile axes towards the trajectory, [6] to [8]. It depends on aerodynamic shape, sensitivities on the drag and lift forces, aerodynamic moments their derivatives and spin rate stabilization efficiency vs. flight Mach numbers, but also of less discovered lateral Magnus moment and other artificial designed moments coupled with them, which changes on the flight trajectory. The importance of experimental aerodynamically accurate estimations coupled with numerical stability criteria simulations is required because the spin stabilized symmetric AA cannon type projectile, exposes sensitive effects on the projectile stability which could influence conclusions about its novel applications and the redesign in the modernization considering processes. A comparative analysis of calculated and experimental parameters is well known, proceeding for a preliminary estimation of all kinds of flight bodies, including flight vehicles with typical gyro - aerodynamics correlation loads, as in the [9] and [10] as rotary wings flight objects. Spin gyro - aerodynamics correlation loads, which strongly influenced on stability of flight projectiles, are referencing on the second class of flight bodies, known as the projectiles. This paper deals with the aerodynamics of the spin projectile which is tested in the wind tunnel. The tested aerodynamics is compared with numerically predicted aerodynamic coefficients and derivatives. Both are used for the stability criteria analysis vs. Mach numbers.The simulated data are used from the AA gun model trajectory designed in [11] to [13]. 2 SPIN STABILIZED CRITERIA AND DATA REQUIREMENTS The key property of spin stabilized projectiles flight could be determined by the quality of its space spherical motion around the projectile body gravity center, during the flight on the ballistic trajectory. A very high spin of longitudinal axis and the low spin of the lateral disturbance provides gyro moment to stabilize projectile pitch. Composed spherical oscillations of body axes and direction of velocity vector, is question of stabilization [6] and [7]. As a result these lateral oscillations decrease perturbed amplitudes of total angle of attack if spin stabilization is successful, or increase if the projectile has not realized enough initial rpm by spin to form gyro-moment for damping. Complex variable of the total angle of attack and its lateral angular motion due to the projectile body can be described by linearized differential equation derived by dimensionless distance instead of time as the main argument [6] in the form: ^ + (H - iP)'-(M + iPT) = -iPG, (1) where H, P, T, M and G represents the so called Murphy's coefficient, [6] developed from the aerodynamic and dynamic solutions of the revolution body with gyroscopic low and high spin coupled motion. Since Eq. (1) is complex but linear, its solution is given as: M^"10 ) + S*. (2) Each value of the angular spin angles as the frequencies ® j, and damping coefficients Xj in homogeneous solution, will vary with the relative magnitudes of the H, P, T and M which could be expressed as the basic function derived from the projectile designed form its aerodynamic, flight dynamic, inertial and all over flight performances. With an aim to separate conditions of damping abilities of a complex angle of attack expressed by j and gyroscopic frequencies effect ® j, the conditions of projectile's dynamic behavior of inertial axes, are enough precise approximately described in [6], [7] and [14], by assuming that 2' >> \X1. This provides using a real solution to be applied and considered in the form: j 2 P ±vP2 - 4M j = 1,2, (3) and damping coefficient separated from joint solution as: ^ =-1 j 2 " P(2T-£) H + i = VP2 - 4M j = 1,2. (4) Real values of both Aj and ® j Eqs. (3) and (4) determine gyroscopic stability criterion in the form: S 4M ~Pr' where, Sg > 1 as the obvious condition, from the Eq. (5), for the real roots in Eq. (3). Using aerodynamic coefficient [6] and [7] these expressions are determined by: t2 2 Ix P where, g 2pIySdV 2Cm P = L.. Pd V m = C'M j;2. (6) (7) (8) The value of Sg is known as the gyroscopic stability factor. A coupled requirement for the stability is that velocity vector angle of attack has to be damped. This is satisfied by the lim£ — 0, which S— requires a necessity that Xj and X2 are both negative Xj < 0, j = 1, 2. This second condition evolves to the so called dynamic criterion of stability Sd as referencing value in Eq. (4), in the form: with values, S _ 2T d H T ~ CLa + CMp a ' H = Cl- CD -(cMq+CM a ) —-2 ~y ' (9) (10) (11) where using aerodynamic properties [1], expression is: 2 (CLa + rx CMpa) (CMq + CM CL a CD -r (12) The condition that damping coefficients Aj have to be negative, relates to the following inequalities, respectively: and H > 0, p 2 (Sd -i)2 P1 - 4M < 1. (13) (14) These Eqs. (13) and (14) coupled with Eq. (5) gyroscopic conditions Sg > 1, finally give stability expression inequality as: ■f < Sd (2-Sd). (15) Eq. (15) describe general criteria of dynamic stability for any axis-symmetric projectile spin or fin equipped. According to [6] to [8], the stability criteria Eqs. (6) and (12) are developed in relation to aerodynamic coefficients and derivatives coupled by inertial and dynamic properties of the body motion. Eqs. (7), (8), (10) and (11) appears as the conditionally for the aerodynamic behavior of projectile's forces and moments as the loadings expressed by the Murphy's coefficients, in the Eq. (1) vs. flight Mach numbers, angle of attack, and the so called derivatives in the changes during flight, [6]. Spin stabilized projectiles during the real flight, change damping coefficients Eq. (4) of characteristic Eq. (2), Xj and X2, and also factors of stability, Sg, given by Eqs. (5) and (6), and Sd, given in Eqs. (9) and (12), which is further derived. The calculated data are used in changes of stability criteria estimations at the expected muzzle distances after firing from the cannon barrel, [13]. For AA cannon projectile of small caliber this spin stabilized behavior vs. flight Mach number will be considered further as a result of experimental testing of measured aerodynamic coefficients and appropriate semi empirical and theoretical data in the stability criteria Eqs. (6) and (12). The aim is to discover influence of the longitudinal position of lateral aerodynamic Magnus force over derivative of Magnus moment coefficient CMpa, on the spin stabilizing factors, which changes its values vs. flight Mach number from high supersonic to the high subsonic values. Considerations have been made using the condition in Eqs. (6) and (12), by the complex simulation in redesigned basic software showed in [8], [11] and [13]. 3 BASE EXPERIMENTAL EQUIPMENT AND TESTING MODEL SET UP The base experimental equipment used for aerodynamic simulation was the three-sonic wind tunnel test facility T-38, [15]. The tunnel is a blow down pressurized type with a 1.5*1.5 m square test section, aimed for subsonic and supersonic tests, [15]. The tunnel was fully equipped by appropriate equipment to simulate flight flow conditions. Mach number can be set and regulated to within 0.5% of the required value. Total pressure in the test is within 1.1*105 to 15*105 Pa regulated to 0.3% errors of real flight conditions. Run times are in the range 6 to 60 s, depending on Mach number and total pressure. The facility supports, step-by-step model movement y Sd = and continuous movement of model ("sweep mode") during the measurements. Simulation and experimentally tested model of projectile in the tunnel is shown on the Fig. 1, [13]. Characteristic values of mass and dimensions, approximately, corresponds to the AA cannon 40 mm HE unguided projectile, Table 1. The model is supported in the test section by a tail sting mounted on a pitch simulation mechanism by which desired aerodynamic angle can be achieved. 5.0* Fig. 1. Model of spin stabilized projectile Table 1. Model of projectile [13] Parameter Symbol, unit Value Reference diameter (cal) d [mm] 40 Total length l [-] ~ 5* Nose length li [-] ~ 2.5* Ogive radius Ro [-] ~ 20* Boat tail length la [-] ~ 0.5* Center of mass XCG [-] ~ 3.3* Mass of projectile m [kg] ~1.0 Axial inertia moment Ix [kg m2] -2.1-10-4 Transversal inertia moment Iy [kg m2] -2.3-10-3 * relative values as number of reference diameter (caliber) Tests on the model were performed in the Mach number range from 0.2 to 3.0 (14 different values of Mach number). The simulated total angles of attack at, redefined to the angle of attack in vertical plane, a, (pitch angle, as a good approximation to the total yaw angle, pp. 33, [6]), were in the interval, -10 to +10°, [12] and [13]. Test conditions were given in Table 2. Aerodynamic forces and moments of the model were measured by ABLE 1.00 MKXXIIIA internal six-component strain gauge balance, [13] and [15]. The nominal load range of the balance was 2800 N for normal, 620 N for side forces, 134 N for axial force, 145 Nm for pitching, 26 Nm for yawing moment and 17 Nm, for static spin damping moment; the accuracy was approximately 0.25% F.S. for each component. Instrumentation and data recording were performed after each run using the standard T38-APS software [15] in several stages, i.e.: Data acquisition system interfacing and signal normalization; Determination of flow parameters in the test section of the wind tunnel; Determination of the model position (orientation) is relative to test section and airflow. Determination of non-dimensional aerodynamic coefficients of forces and moments for each stage has been performed by appropriate tunnel tests, different software modules, [17]. Table 2. Test conditions [13] Parameter Symbol, unit Value static pressure Ps [Pa] 0.2x105 to 2.2x105 stagnation pressure P0 [Pa] 2.3x105 to 6x105* atm. temperature Tatm [K] « 280 Mach number Ma [-] 0.2 to 3 Re number Re [-] 0.5x106 to 2x106 angle of attack (pitch) a [°] -10 to 10 * for all Mach numbers was 2.3-105 Pa, since for Ma = 2.5, p0 =4 105 Pa and Ma = 3, p0 = 6105 Pa. 4 EXPERIMENTAL RESULTS AND SPIN STABILIZED MODELING 4.1 Aerodynamic Coefficients and Derivatives The characteristic functions of aerodynamic coefficients in relation to the flight Mach numbers and angles of attack a, measured in the author's testing, [13] in the wind tunnel is presented in Figs. 2 to 5, and Fig. 6 represents additional data supposed from [14]. The influenced aerodynamic coefficients tested on the projectile model are aerodynamic coefficient of drag force, Fig. 2, derivative of lift force coefficient presented in Fig. 3, derivative of pitching moment coefficient given in Fig. 4, and all vs. flight Mach number. Flight Mach numbers were corresponding to the projectile flight velocities on the modeling trajectories. Aerodynamic prediction or calculations of coefficients and derivatives are determined with two semi-empirical methods: ADK0 for zero angle drag coefficient and ADK1 for others, [13]. Research was developed according to [11] and [16] to [18]. Data of these predictions for appropriate aerodynamic coefficients are also presented in the same figures and compared with the above mentioned experimental values. Aerodynamic coefficients vs. side slip component of total yaw angle p, was not tested experimentally and further considerations took these effects in integral yaw angle, by semi-empirical predictions denoted as ADK1 in the previously mentioned references. Measurements of the so called dynamic derivatives of pitch damp coefficient in wind tunnel facilities require complex and expensive testing equipment, and improvement of the test model design, which were not used in these experiments. Experiments of pitch derivatives, vs. flight Mach numbers are approximately determined in the tunnel test simulations using the test model pitch motion, with threshold ability of angular rate in sweep-mode as was 2 degrees per second, [12] and [13]. Values of derivatives of aerodynamic coefficient CMq+ CMa , realized in these experiments are presented in Fig. 5, by dot-points curve realized in the singular flight test runs, vs. 14 values of simulated flight Mach numbers. A comparison of these values with the data calculated by ADK1 prediction is also presented in the same figure. o 3 Fig. 2. Predicted values (ADK0) and experimental values (exp.) of drag coefficients 5 4 E 3 8 2 •■exp. —ADK1 ------ * * ..........o ... y ■■ -o-O"O0 0.5 1.5 Ma [- 2.5 Fig. 3. Predicted values (ADK1) and experimental values (exp.) of derivative of lift force coefficient The derivative of Magnus moment coefficient, CMpa, which could not be simulated by the wind tunnel was estimated by calculations from the approximated model ADK1, [13] according to the developed estimations in [6], [7] and [14], Fig. 6. Variations are accepted based on data in [6], [7], [11] and [14]. It was the base challenge in estimation because Magnus force effects influenced the stability similarly as any other lateral force and is composed of flight Mach number effect and spin peripheral velocity designed in the stream flow. This aerodynamic loading coefficient based on the relative small values of real Magnus forces could make undetermined problems to the stability of projectiles if the force lateral position along longitudinal axes is not well known, or vary vs. flight Mach numbers changes. Flight Mach numbers and other data are changed on the simulated trajectory by the six degrees of freedom software (6DOF) [6], [8], [11] and [13] modeling. 8 6 - 4 U 2 *exp. -ADK1 Q...... ,9 i,rO.....« .... ......... -o. ,. . .c 0 0.5 1 1.5 Ma [- 2 2.5 3 Fig. 4. Predicted values (ADK1) and experimental values (exp.) of derivative of pitching moment coefficient •♦•exp.sim. —ADK1 - PO> # HH] 0 0.5 1 1.5 2 2.5 3 Ma [-] Fig. 5. Predicted values (ADK1) and simulated experimental values (exp.sim.) of dynamic derivative of pitch damping moment 1.5 0.5 l n ii ii ii 11 —calc.det. --calc.undet. -*-calc.det.2 ■A-calc.undet.2 \ i i i i ( H' i.': i.'1. 1.5 2.5 Fig. 6. Predicted values (ADK1) of derivative of Magnus moment coefficient [14] The flight conditions are given in Table 3. Aerodynamic data have been estimated by ADK0, ADK1 and presented from experimental tests.The estimated main derivatives were assumed as the challenge for sensitivity tests of the so called stability vs. flight Mach number, represented velocities on any type of trajectory. 1 2 3 5 2 3 Table 3. Flight simulation condition [13] Parameter Symbol, unit Value pressure Pa [Pa] 1.013 105 temperature Ta [K] 288.13 muzzle velocity V [m/s] 985 x component of muzzle angular velocity - spin P [s-1] 6300 y component of muzzle angular velocity q [s-1] 1.0 ** z component of muzzle angular velocity r [s-1] 7.7 ** Table range angle (gun elevation) fo [°] 5 to 20 ** according to the research of the initial fire disturbing conditions, [11] and [13] 4.2 Spin Stability Parameters Modeling Qualitative evaluation of projectile stability is determined through an analysis of the simulated data using software 6DOF, [11] and [13], expressed in the following: • absolute values of total angle of attack |at| vs. time, Fig. 7, which corresponds to the solution of in Eq. (2), with zero initial £g, • damping behavior A1,2, of high and law spin of complex |£| module, vs. stability factors relation Figs. 8a and b, all vs. flight Mach number, • gyroscopic stability factor Sg and criteria relation, Fig. 9, • spin angular velocity damping p, Fig. 10, • dynamic stability factor Sd and criteria, Figs. 11 and 12, • total stability criteria Fig. 13. —calc. .....exp. i S i: ! i H'y n i •¡Pi < ii ï ; iîfï 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 t [s] Fig. 7. Absolute simulated angle of attack on the initial part of trajectory Two groups of aerodynamic data in flight projectile stability modeling are used, as well as calculated data, (ADK1) denoted in Figs. 8 to 13, as the (calc.) and experimental aerodynamic data used from Figs. 2 to 6, composed in the matrix form adapted for the software, and denoted as the (exp.) roots simulations. Experimental rooted data, of a, corresponded to the tunnel measured values vs. angles of attack a. Complex yaw simulations taken from the calculated data and with a total yaw angle at, composed from vertical (pitch) angle of attack a are simulated to prove the damping effect of velocity direction to the projectile body, which was preliminarily designed by simulation without derivative of Magnus moment coefficient, which means Magnus force point is in the center of mass. Representative damping quality should be exposed at the very beginning of flight within less than half a second, Fig. 7. These influences caused by Murphy's coefficient H were dominant aerodynamic derivatives calculated estimations and experiments shown in Fig. 5, as the summary dynamic derivatives of coefficient CMq+ CMa, in both, supersonic and subsonic flight Mach numbers, expressed deviations of more than 30 percent for calculated and experimental data. This causes a strong influence on Eq. (11), in the simulation done using static tunnel test data, which was performed by inappropriate derivative measurements. Differences of experimental and calculated data from the initial flight Mach numbers of about 3 to 1.5, are caused in the derivative of pitch damp moment coefficient CMa, Fig. 4, which strongly influenced Sg estimations by calculated and experimental data, Fig. 9. This was significantly exposed for the high initial flight Mach numbers. Approval is shown in Fig. 9 where both gyro stability curves factors and criteria of the calculated and experimental roots make a crossing in the Mach number of 1.5 values. This is in correlation to the values of Fig. 4. Other related influences on the differences of the angle like the Magnus force and their lateral position along projectile longitudinal axis, affected the side-slip component ft in the total yaw angle at were not considered but were performed through influencing estimations in the stability factors and damping performances in the gyro stabilization. The angle damping behavior vs. flight Mach numbers, Fig. 7 correlated with the measurement in flight tunnel tests and calculated data with the derivative of Magnus moment coefficient CMpa, Fig. 6, influences are presented in Fig. 8b, expresses negative and 8a, one negative and one positive values of damping coefficients Eq. (4). These data have shown disastrous influencing of derivative of Magnus moment coefficient CMpa, which was tested, using Fig. 6 values by 100% variations, keeping the same other conditions in the stability tests. This means more an unknown position of Magnus force along projectile axis than their intensity values. The 2 0.5f \ ! -0.5> -1> ; -1.5- s f ^ I i -2.5-- —(calc) (calc) -a-^1 (exp) •a-^2 (exp) a) Ma [-] (cac (CMpa /2) (cac) (CMpa /2) (exp) (CMpa /2) 2 (exp) (CMpa/2) 'ïa; ...... b) Ma [-] Fig. 8. Damping coefficients vs. Mach number; a) for full value of derivative of Magnus moment coefficient, b) for half value of derivative of Magnus moment coefficient 7 is 6 5 ° 2 —S (calc) g ■ e-S (exp) ...A • 'o a) 0.5 „0.6 .g | 0.5 £ o £ 0.4 -S 0 ^ 0.3 0 1 ^ 0.2 1.5 2 Ma [-] 2.5 "0.1 J». — 1/S (calc) g - ■ •■ 1/Sg (exp) b) 1.5 2 Ma [-] 2.5 Fig. 9. Gyroscopic stability relation vs. Mach number; a) gyroscopic stability factor, b) gyroscopic stability criteria differences between experimental and the calculated data vanished by decreasing this derivative vs. Mach number distribution, Fig. 8b. This effect is caused by coupling with estimations of H, Murphy's coefficient, Eq. (10), which remains positive values, Eq. (13), but vary, making influences on the X12, more on the high 6500 6000 5500 5000 ^4500 ^ 4000 3500 3000 --(1) exp. 00=5° -----(2) exp. 00=10° - .......(3)exp. 00=15° - —(4) calc. 00=5° -—(5) calc. 00=10° .......(6) calc. 00=15° (4k A(f) (5)x „-(2) ^ (3) Ma [-] Fig. 10. Angular spin velocity vs. Mach number 'E 1 -is 'E & Ü 0 -S/2-Sd ) (calc) (CMpa ) ■"■Sd(2-Sd) (exp) (CMpa) ..» *........ a) 1.5 2 Ma [-] 'C 1 'E & it 0 -Sß-Sd ) (calc) (CMpa /2) ■°Sß-Sd) (exp) (CMpa 2 1.5 2.5 b) Ma [-] Fig. 11. Dynamic stability relation criteria vs. Mach number; a) for full value of derivative of Magnus moment coefficient, b) for half value of derivative of Magnus moment coefficient frequency damping coefficient A1 then on the X2, low frequency damping coefficient, both expressed in the in Eq. (4). Gyroscopic stability factor determined by the changes of flight Mach number and the data of the projectile model are presented in Fig. 9 demonstrates flight gyro-stability factor as a relation of Eq. (6) vs. Mach numbers not violated by the derivatives. The spin velocity, projectile performances are taken in both data types of simulations representative trajectory with equal values vs. Mach number, represented in Fig. 10, and corresponding with projectile model data in Tables 1 and 3. Changes of spin were not tested experimentally. The most significant influences of the derivative of Magnus moment coefficient CMpa, Fig. 6, were 2500 5 1.5 2 2.5 3 1.5 2 2.5 3 x 10 -0.5 2 -1.5 5 1.5 2 2.5 3 5 2.5 3 2 4 3 5 2 3 1 3 1 3 exposed as disastrous on the estimated calculations of dynamic stability factors. In the presented Figs. 11 and 12, the condition of dynamic stability Eq. (15), and stability factor given by Eq. (12) with Magnus moment derivative influences were tested. Influencing of the derivative of Magnus moment coefficient CMpa, which varied 100% in Fig. 6, making the stability of the dynamic factor bigger than 2 thus, making Eq. (15) negative, Fig. 12, and violating stability. Double less values of CMpa vs. flight Mach numbers rearranged the stability of dynamic factors, but making neutral stability approximately. Effects are related to the values of H also by modestly reliable data, taken as the sum of derivatives CMq+ CMa . The sum is shown as influence on the dynamic stability in Fig. 12, with respectable differences between the simulated stability criteria in Eq. (12), with experimental and calculated data vs. Mach number. Values are exposed as the coupled with CLa influences and the above mentioned derivatives. Fig. 12. Influences of derivative of Magnus moment coefficient on dynamic stability factor t -l- ...,!.. 1 /(5)/(6) ............... ,»-*............ "•'■'.'."...........1..................... 1-,—J: • -» y« - . - --— m / (4)' / —(1) 1/Sg (calc) -»-(2) 1/S (exp) .....(3) Sd(2-SJ (calc) -+-(4) S/2-SJ (exp) ^(5) Sd(2-SJ (calc)(CMpa/2) -■(6) SJ(2V (exP)(CMpa/2) (3)/ (2) (l) Ma [-] Fig. 13. Total stability factors equation influences for full and half derivative of Magnus moment coefficient vs. Mach number denoted expected total stability of projectile in the flight. The coupled effect in Fig. 13 proves that the dominant role of the estimated positive values CMpa of about 0.8 to 0.4, in supersonic and subsonic flight has a negative effect on the projectile flight stability. Table 4. Approximate estimations of aerodynamic parameters' influences on stability parameters [13] AD Manner of AD Stability parameter parameter parameter Sensitivity of stability parameter increase to 20% increasing 1/Sg increase to 20% Sd(2 - Sd) increase to 5% p decrease 3 to 8% ¿2 decrease to 20% decreasing 1/Sg decrease to 20% Sd(2 - Sd) decrease to 5% p increase 3 to 8% increasing ¿2 increase to 20% C La Sd(2 - Sd) increase to 5% decreasing Ä2 decrease to 20% Sd(2 - Sd) decrease to 6% CMa increasing 1/Sg increase to 15% decreasing 1/Sg decrease to 15% increasing Ai decrease to 30% CMq+ Sd(2 - Sd) decrease 2.5 to 7% C Ma decreasing Ai increase to 30% Sd(2 - Sd) increase to 8% increasing Ai, A2 increase to 30% C ^Mpa Sd(2 - Sd) increase to 25% decreasing Ai, A2 increase to 30% Sd(2 - Sd) increase to 25% These data correspond with [13] but do not fully correspond with [14], and are directed to general conclusions in [6], which are that Magnus moment effect could be expected to be suspicious if it is taken from the data which are not proved by real measurements which are missed in the published data. The paper confirms estimations of instability boundaries and area of main derivative coupling influenced on the flight. Lateral force in the center of gravity proves determination of the stable flight. An Increase of lateral moment derivatives caused by a variation of the force center along the projectile axis is not suggested. Sensitivity of stability factors tested in [13] and in this paper is presented in Table 4. Dynamic stability factor Sd and gyroscopic stability factor Sg related by Eq. (15) are presented in Fig. 13 as the final value of stability testing influences vs. Mach number tests and simulations. The main differences of the experimental and calculated data 5 CONCLUSION The stability of AA gun projectile model is affected by aerodynamic coefficients and derivatives through their steady state testing data. The test performed by a C 0 0.5 1.5 2 2.5 3 an approximately real designed model of projectile was used, and the obtained measurements are valid but not enough for the full reliable estimations of flight stability. The trait of the change of the based calculated values of aerodynamic coefficients coincides with the experimental results of aerodynamic coefficients but does not satisfy the required conclusions about the stability without derivative measurements. A comparative calculation model of coefficient estimations is tested to prove possible indeterminations which could appear as significant, but were not present in the measurement on the tunnel tests at approximately steady state conditions. Further research of ammunition correction functions could be possible for preliminary consideration by simulated stability methodology of lateral forces required in correction by the method of coefficient derivatives shown in this paper. Damping efficiency of projectile initial magnitude of angle velocity vector at (total yaw angle) is about 0.4 seconds after launching, which corresponds to the expected if lateral force is in the center of mass. Also, sensitivity tests give satisfied frame values of gyroscopic stability and dynamic stability influences for the approximated values of Magnus moments and other similar disturbances representatively arranged in the coefficient derivatives values. Further research should comprise of simulations and a test to develop the best way for ammunition guidance considerations based on the behavior of the main projectile axes and velocity vector during flight. Further dynamic testing of aerodynamic derivatives requires appropriate tunnel facility including equipment for a good simulation of the angular motion for ammunition models tests. 6 AKNOWLEDGMENT The results presented in this paper are made within project No III47029 financed by Ministry of Science and Technology Development of the Republic of Serbia. 7 REFERENCES [1] Regan, F.J. (1975). Aeroballistics of terminally corrected spinning projectile (TCSP). Journal of spacecraft, vol. 12, no. 12, p. 733-738. [2] Morrison, P.H., Amberntson, D.S. (1977). Guidance and control of a cannon-launched guided projectile. Journal of spacecraft, vol. 14, no. 6, p. 328-334. [3] Zarchan, P. (1990). Tactical and strategic missile guidance. Progress in Astronautics and Aeronautics, AIAA, vol. 124. [4] Vogt, R., Glebocki, R. (2000). Impulse control of anti-tank mortar missile. RTO AVT Symposium on Active Control Technology for Enhanced Performance Operational Capabilities of Military Aircraft, Land Vehicles and Sea Vehicles, Braunschweig. [5] Gupta, S.K., Saxena, S., Singhal, A., Ghosh, A.K. (2008). Trajectory correction flight control system using pulse jet on artillery rocket. Defense Science Journal, vol. 58, no. 1, p. 15-33. [6] McCoy, R.L. (1999). Modern exterior ballistics. Schiffer Military History, Atglen. [7] Carlucci, D.E., Jacobson, S.S. (2008). Ballistics -Theory and design of guns and ammunition. CRC Press, Boca Raton. [8] Regodic, D. (2006). Exterior ballistics. Military academy, Belgrade. (in Serbian) [9] Mitrovic, Č., Bengin, A., Cvetkovic, D., Bekric, D. (2010). An optimal main helicopter rotor projection model obtained by viscous effect and unsteady lift stimulation. Strojniški vestnik - Journal of Mechanical Engineering, vol. 56, no. 6, p. 357-367. [10] Petrovic, Z., Stupar, S., Kostic, I., Simonovic, A. (2010). Determination of a light helicopter flight performance at the preliminary design stage. Strojniški vestnik - Journal of Mechanical Engineering, vol. 56, no. 9, p. 535-543. [11] Regodic, D. (2003). Exercises in exterior ballistics. Military academy, Belgrade. (in Serbian) [12] Jerkovic, D., Samardzic, M. (2008). The aerodynamic characteristics determination of classic symmetric projectile. The 5th international symposium about design in mechanical engineering, KOD, p. 275-282. [13] Jerkovic, D. (2009). The influence of aerodynamic coefficients on the motion of axis-symmetrical body. Faculty of technical sciences University of Novi Sad. (in Serbian) [14] Regan, F.J., Schermerhorn, V.L. (1971). Aeroballistic evaluation and computer stability analysis of the US Navy 20-mm general purpose projectile. NOLTR, p. 71-95. [15] Elfstrom, G.M., Medved, B. (1986). The Yugoslav 1.5 m trisonic blowdown wind tunnel. 14th Aerodynamic Testing Conference. West Palm Beach, AIAA Paper 8624726, p. 89-95. [16] Krasnov, N.F. (1985). Aerodynamics 1 - Fundamentals of Theory. Mir Publishers, Moscow. [17] Krasnov, N.F. (1985). Aerodynamics 2 - Methods of Aerodynamic Calculations. Mir Publishers, Moscow. [18] Jankovic, S. (1979). Aerodynamics of projectiles, Faculty of mechanical engineering, Belgrade (in Serbian) [19] McCoy, R.L. (1985). Aerodynamic and flight dynamic characteristics of the new family of 5.56 mm NATO ammunition. Memorandum report BRL-MR-3476, US Army Ballistic Research Laboratory, Aberdeen Proving Ground, Maryland.