AD1089736
AD1089736
AD1089736
Disclaimers
The findings in this report are not to be construed as an official Department of the
Army position unless so designated by other authorized documents.
14. ABSTRACT
Spin balancing is used to determine the mass imbalance of rotating machinery parts. Recently, the US Army Combat
Capabilities Development Command Army Research Laboratory commissioned the development and construction of a
spin-balance machine for measuring the mass imbalance of projectiles. Mass imbalance is one source of error that contributes
to increased target impact dispersion and reduced hit probability. Such a machine would be useful in assessing and quantifying
the mass imbalance of as-manufactured projectiles or recovered postlaunch projectiles. The current report presents the theory
behind the operation of these machines and how they are used to measure the static and dynamic imbalance of projectiles. The
report also provides analysis that links the measured static and dynamic imbalance to the trajectory deflection produced by the
mass imbalance. Results of measured imbalance of a set of calibration pins are presented to benchmark the measurement
capability. The methodology is also used to predict the trajectory deflection of a 20-mm projectile with a mass asymmetry,
and the predictions are validated with previously published results.
15. SUBJECT TERMS
projectile mass imbalance, dispersion, trajectory deflection, spin balancing, lateral throw-off
17. LIMITATION 18. NUMBER 19a. NAME OF RESPONSIBLE PERSON
16. SECURITY CLASSIFICATION OF: OF OF
Paul Weinacht
ABSTRACT PAGES
a. REPORT b. ABSTRACT c. THIS PAGE 19b. TELEPHONE NUMBER (Include area code)
UU 44
Unclassified Unclassified Unclassified (410) 306-0800
Standard Form 298 (Rev. 8/98)
Prescribed by ANSI Std. Z39.18
ii
Contents
List of Figures iv
List of Tables iv
Acknowledgment v
1. Introduction 1
7. Conclusion 29
8. References 30
Distribution List 37
iii
List of Figures
Fig. 1 ARL spin-balance machine ................................................................... 2
Fig. 2 Close-up of ARL spin-balance machine ............................................... 3
Fig. 3 Planform view of spin-balance machine setup (facing front of spin-
balancing machine) ............................................................................... 4
Fig. 4 Cross-sectional view of spin-balance machine setup viewed from left
side of fixture ........................................................................................ 5
Fig. 5 Pickup forces on left and right pickups................................................. 9
Fig. 6 Ballistic- and body-fixed coordinate systems ..................................... 17
Fig. 7 Circumferential angles in body-fixed and spin-balance coordinate
systems (view from nose) ................................................................... 18
Fig. 8 20-mm projectile with mass offset ...................................................... 22
Fig. 9 Schematic of calibration Pin C (dimensions in inches) ...................... 25
Fig. 10 Calibration pins ................................................................................... 26
List of Tables
Table 1 Sample parameters for spin balance measurement ............................... 9
Table 2 Hypothetical (theoretical) imbalance and spin-balance correction ..... 23
Table 3 Static and couple imbalances for 20-mm projectile ............................ 24
Table 4 Lateral throw-off and aerodynamics jump for 20-mm projectile ....... 24
Table 5 Physical properties of mass asymmetry for the four calibration pins . 26
Table 6 Measured and theoretical mass corrections at left and right correction
planes for the four pins ....................................................................... 27
Table 7 Measured and theoretical mass corrections at left and right correction
planes .................................................................................................. 28
Table 8 Measured and theoretical static and couple imbalances for the four
calibration pins .................................................................................... 28
Table 9 Comparison of measured and theoretical values of static and couple
imbalance ............................................................................................ 29
iv
Acknowledgment
The original inspiration for the investigation and eventual procurement of the
spin-balance machines by US Army Combat Capabilities Development Command
Army Research Laboratory (ARL) was made by Jim Garner, Flight Sciences
Branch, ARL. Subsequent application of these machines to measure the spin
imbalance of projectiles led to the need to more fully understand the operation of
these machines and how the measured imbalance quantities could be more directly
tied to projectile trajectory deflection. This is the subject of this report.
v
1. Introduction
1
measurements of the spin imbalance of a set of calibration pins are presented in
Section 6 to further benchmark the measurement capability.
2
Fig. 2 Close-up of ARL spin-balance machine
3
Fig. 3 Planform view of spin-balance machine setup (facing front of spin-balancing
machine)
The part being measured rides on two V-blocks that are adjusted so that the part
rides smoothly. The precise location of the V-blocks is not measured in the
spin-balance procedure and does not directly enter into the measurement
calculations. There does not appear to be a requirement that the cylindrical portion
of the projectile ride on the V-blocks. In other words, the projectile could ride on
the forward ogive and/or the aft boattail. A belt is used to rotate the part in the
V-blocks. The belt provides a constant downward force 𝑃𝑃 at a location 𝐿𝐿𝑃𝑃 on the
measured part. The weight of the part 𝑊𝑊 acts at the center of gravity of the part
𝐿𝐿𝐶𝐶𝐶𝐶 . The belt force and part weight and their location appear in the analysis
discussed below, but do not directly play into the measurement of the imbalance to
the extent that these are quantities would need to be explicitly specified for the
spin-balance machine measurements.
The spin-balance machine measures the forces 𝑅𝑅𝑙𝑙 and 𝑅𝑅𝑟𝑟 at the two pickup
locations 𝐿𝐿𝑙𝑙 and 𝐿𝐿𝑟𝑟 . Based on the measured forces from the pickups, the
spin-balance machine computes the imbalance. In many applications, the
spin-balance machine would be used to correct the imbalance on the part. This can
be done by adding or removing mass from the part at two locations (correction
planes) along the length of the part shown as 𝐿𝐿𝑙𝑙𝑙𝑙 and 𝐿𝐿𝑟𝑟𝑟𝑟 in Fig. 3. Since the
correction is typically made at the surface of the part, the diameter of the part at the
correction plane 𝐷𝐷𝑙𝑙𝑙𝑙 and 𝐷𝐷𝑟𝑟𝑟𝑟 is also specified. The corrections at each correction
plane are also made at specified circumferential locations ∅𝑙𝑙𝑙𝑙 and ∅𝑟𝑟𝑟𝑟 , as shown in
Fig. 4. The correction planes are suitable locations selected for the imbalance to be
4
corrected without interfering with the operation of the part. For example, for
projectile applications the cylindrical portion of the projectile may be a poor
location for the correction to be made because this portion of the projectile contacts
the gun barrel during the firing process, and any correction may be eliminated by
the engraving process during firing. If the purpose of the imbalance measurement
was simply to define the degree of imbalance, the selection of the location of the
correction planes as well as the diameter at the correction plane are somewhat
arbitrary.
Fig. 4 Cross-sectional view of spin-balance machine setup viewed from left side of fixture
Figure 3 shows the part being measured with the nose to the left when viewed from
the front of the machine. The fixture is capable of measuring parts and correcting
imbalances with the nose oriented to either the left or right. However, if the part is
measured with the nose right, care must be taken in analyzing the resulting
measurements, particularly if they are used to quantify the ballistic performance,
since the analysis presented here assumes that the nose is oriented to the left. At a
minimum, the measured angular orientations ∅𝑙𝑙𝑙𝑙 and ∅𝑟𝑟𝑟𝑟 will likely have to be
reversed in sign.
5
3. Theoretical Background of the Spin-Balance Machine
Operation
A given part may have an arbitrary number of mass asymmetries located at various
locations throughout the part. Although it may be difficult to determine the exact
magnitude and location of each asymmetry, the analysis here provides the first step
in representing the net effect of an arbitrary number of asymmetries within the part.
A summation of the forces in the vertical direction and the moments acting normal
to the vertical plane at a point 𝐿𝐿1 along the centerline of the part is shown in Eqs. 1
and 2.
𝑖𝑖=𝑁𝑁
6
𝑖𝑖=𝑁𝑁
�𝑀𝑀𝑍𝑍 |𝐿𝐿𝑙𝑙 = −𝑃𝑃(𝐿𝐿𝑃𝑃 − 𝐿𝐿𝑙𝑙 ) − 𝑊𝑊(𝐿𝐿𝐶𝐶𝐶𝐶 − 𝐿𝐿𝑙𝑙 ) + 𝑅𝑅𝑟𝑟 (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) + � 𝐹𝐹𝑖𝑖 (𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙 )cos(∅𝑖𝑖 + ∅) = 0 (2)
𝑖𝑖=1
𝐹𝐹𝑖𝑖 is the centrifugal force due the individual mass asymmetry, 𝜔𝜔 is the rotation
speed of the part, 𝑚𝑚𝑖𝑖 is the mass of each individual asymmetry and 𝑟𝑟𝑖𝑖 is the radial
location of the asymmetry from the axis of rotation, ∅𝑖𝑖 is the circumferential
location of the asymmetry from the vertical axis, and ∅ is the circumferential
location of the part from the vertical axis as it rotates about its longitudinal axis.
The mass asymmetry is considered a positive mass, but voids in a part could be
considered as negative mass. 𝑈𝑈𝑖𝑖 is the imbalance due to the individual mass
asymmetry.
Using Eqs. 1 and 2, the pickup forces 𝑅𝑅𝑙𝑙 and 𝑅𝑅𝑟𝑟 can be determined.
𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑃𝑃 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝐶𝐶𝐶𝐶 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙 )
𝑅𝑅𝑟𝑟 = 𝑃𝑃 + 𝑊𝑊 − � 𝐹𝐹𝑖𝑖 cos(∅𝑖𝑖 + ∅) (4)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
𝑖𝑖=1
𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝐶𝐶𝐶𝐶 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙 )
𝑅𝑅𝑙𝑙 = 𝑃𝑃 �1 − � + 𝑊𝑊 �1 − � − � 𝐹𝐹𝑖𝑖 �1 − � cos(∅𝑖𝑖 + ∅) (5)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿1 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
𝑖𝑖=1
7
𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙 )
� 𝑟𝑟 � = − � 𝐹𝐹𝑖𝑖 sin(∅𝑖𝑖 )
𝑅𝑅�𝑟𝑟 sin�∅ (9)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
𝑖𝑖=1
2 2
𝑖𝑖=𝑁𝑁 𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙 )
𝑅𝑅�𝑟𝑟 = ��� 𝐹𝐹𝑖𝑖 cos(∅𝑖𝑖 ) � + �� 𝐹𝐹𝑖𝑖 sin(∅𝑖𝑖 ) � (10)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
𝑖𝑖=1 𝑖𝑖=1
𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙 )
� 𝑙𝑙 � = − � 𝐹𝐹𝑖𝑖 cos(∅𝑖𝑖 ) �1 −
𝑅𝑅�𝑙𝑙 cos�∅ � (11)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
𝑖𝑖=1
𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙 )
� 𝑙𝑙 � = − � 𝐹𝐹𝑖𝑖 sin(∅𝑖𝑖 ) �1 −
𝑅𝑅�𝑙𝑙 sin�∅ � (12)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
𝑖𝑖=1
2 2
𝑖𝑖=𝑁𝑁 𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙 )
𝑅𝑅�𝑙𝑙 = ��� 𝐹𝐹𝑖𝑖 cos(∅𝑖𝑖 ) �1 − �� + �� 𝐹𝐹𝑖𝑖 sin(∅𝑖𝑖 ) �1 − �� (13)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
𝑖𝑖=1 𝑖𝑖=1
�𝑟𝑟 and �
The phase angles of the imbalance ∅ ∅𝑙𝑙 are determined from the appropriate
application of the arctangent function using Eqs. 8 and 9 or Eqs. 11 and 12,
respectively.
The pickup forces for a hypothetical mass offset case that contains four individual
mass asymmetries (parameter listed in Table 1) are shown in Fig. 5. As shown in
the previous analysis, there is a mean contribution from the belt tension and part
weight and a fluctuating component from the mass asymmetries that is sinusoidal
in nature. The amplitude of the fluctuating component of the pickup forces varies
with the square of the rotation rate of the spin-balance machine. As a result,
sufficient spin rate must be applied for a measurable fluctuating component to be
accurately measured.
8
Table 1 Sample parameters for spin-balance measurement
9
3.2 Equations of Motion of Mass Offsets at Two Correction
Planes
�𝑀𝑀𝑍𝑍 |𝐿𝐿𝑙𝑙 = −𝑃𝑃(𝐿𝐿𝑃𝑃 − 𝐿𝐿𝑙𝑙 ) − 𝑊𝑊(𝐿𝐿𝐶𝐶𝐶𝐶 − 𝐿𝐿𝑙𝑙 ) + 𝑅𝑅𝑟𝑟 (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) + 𝐹𝐹𝑙𝑙𝑙𝑙 (𝐿𝐿𝑙𝑙𝑙𝑙 − 𝐿𝐿𝑙𝑙 ) cos(∅𝑙𝑙𝑙𝑙 + ∅)
(15)
+𝐹𝐹𝑟𝑟𝑟𝑟 (𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) cos(∅𝑟𝑟𝑟𝑟 + ∅) = 0
Here, 𝐹𝐹𝑙𝑙𝑙𝑙 and 𝐹𝐹𝑟𝑟𝑟𝑟 are the forces due to the mass imbalance at the left and right
correction planes; 𝑚𝑚𝑙𝑙𝑙𝑙 and 𝑚𝑚𝑟𝑟𝑟𝑟 are the equivalent masses of the imbalance at the
left and right correction planes; 𝐿𝐿𝑙𝑙𝑙𝑙 and 𝐿𝐿𝑟𝑟𝑟𝑟 are the axial locations of the left and
right correction planes; and 𝑟𝑟𝑙𝑙𝑙𝑙 and 𝑟𝑟𝑟𝑟𝑟𝑟 and ∅𝑙𝑙𝑙𝑙 and ∅𝑟𝑟𝑟𝑟 are the radial and
circumferential locations, respectively, of the imbalance at the left and right
correction planes. The axial location and radial distance (actually input as a
diameter) at both correction planes are required to be input into the spin-balance
machine software. If the spin-balance machine is to be used to measure the
imbalance and no correction is to be made to the part, the selection of the axial
location and radial distance of the correction at both correction planes is somewhat
arbitrary:
𝐹𝐹𝑙𝑙𝑙𝑙 = 𝜔𝜔2 𝑚𝑚𝑙𝑙𝑙𝑙 𝑟𝑟𝑙𝑙𝑙𝑙 = 𝜔𝜔2 𝑈𝑈𝑙𝑙𝑙𝑙 (16)
The pickup forces due to mass imbalance at the two corrections planes is shown in
the following equations:
(𝐿𝐿𝑃𝑃 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝐶𝐶𝐶𝐶 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑙𝑙𝑙𝑙 − 𝐿𝐿𝑙𝑙 )
𝑅𝑅𝑟𝑟 = 𝑃𝑃 + 𝑊𝑊 − 𝐹𝐹𝑙𝑙𝑙𝑙 cos(∅𝑙𝑙𝑙𝑙 + ∅)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
(18)
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
− 𝐹𝐹𝑟𝑟𝑟𝑟 cos(∅𝑟𝑟𝑟𝑟 + ∅)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
10
(𝐿𝐿𝑃𝑃 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝐶𝐶𝐶𝐶 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑙𝑙𝑙𝑙 − 𝐿𝐿𝑙𝑙 )
𝑅𝑅𝑙𝑙 = 𝑃𝑃 �1 − � + 𝑊𝑊 �1 − � − 𝐹𝐹𝑙𝑙𝑙𝑙 �1 − � cos(∅𝑙𝑙𝑙𝑙 + ∅)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
(19)
− 𝐹𝐹𝑟𝑟𝑟𝑟 �1 − � cos(∅𝑟𝑟𝑟𝑟 + ∅)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
The equivalent mass imbalance at the two correction planes can be related to the
arbitrary mass imbalance, as shown in Eqs. 20–25:
𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑟𝑟𝑟𝑟 )
𝐹𝐹𝑙𝑙𝑙𝑙 cos(∅𝑙𝑙𝑙𝑙 ) = − � 𝐹𝐹𝑖𝑖 cos(∅𝑖𝑖 ) (20)
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
𝑖𝑖=1
𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑟𝑟𝑟𝑟 )
𝐹𝐹𝑙𝑙𝑙𝑙 sin(∅𝑙𝑙𝑙𝑙 ) = − � 𝐹𝐹𝑖𝑖 sin(∅𝑖𝑖 ) (21)
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
𝑖𝑖=1
𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
𝐹𝐹𝑟𝑟𝑟𝑟 cos(∅𝑟𝑟𝑟𝑟 ) = � 𝐹𝐹𝑖𝑖 cos(∅𝑖𝑖 ) (22)
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
𝑖𝑖=1
𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
𝐹𝐹𝑟𝑟𝑟𝑟 sin(∅𝑟𝑟𝑟𝑟 ) = � 𝐹𝐹𝑖𝑖 sin(∅𝑖𝑖 ) (23)
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
𝑖𝑖=1
2 2
𝑖𝑖=𝑁𝑁 𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑟𝑟𝑟𝑟 ) (𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑟𝑟𝑟𝑟 )
𝐹𝐹𝑙𝑙𝑙𝑙 = ��� 𝐹𝐹𝑖𝑖 cos(∅𝑖𝑖 ) � + �� 𝐹𝐹𝑖𝑖 sin(∅𝑖𝑖 ) � (24)
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 ) (𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
𝑖𝑖=1 𝑖𝑖=1
2 2
𝑖𝑖=𝑁𝑁 𝑖𝑖=𝑁𝑁
(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙𝑙𝑙 ) (𝐿𝐿𝑖𝑖 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
𝐹𝐹𝑟𝑟𝑟𝑟 = ��� 𝐹𝐹𝑖𝑖 cos(∅𝑖𝑖 ) � + �� 𝐹𝐹𝑖𝑖 sin(∅𝑖𝑖 ) � (25)
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 ) (𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
𝑖𝑖=1 𝑖𝑖=1
The phase angles of the imbalance ∅𝑙𝑙𝑙𝑙 and ∅𝑟𝑟𝑟𝑟 are determined from the appropriate
application of the arctangent function using Eqs. 20 and 21 or Eqs. 22 and 23,
respectively.
As shown in Eqs. 6 and 7, the pickup forces contain a constant force due to the
force from the belt P and the weight of the part W. There is also a fluctuating
component of the pickup forces (Eqs. 6 and 7). As shown in Eqs. 26–29, the
� 𝑟𝑟 and ∅
amplitude of the pickup forces 𝑅𝑅�𝑙𝑙 and 𝑅𝑅�𝑟𝑟 and associated phase angles ∅ � 𝐿𝐿
can be used to determine the mass imbalance 𝐹𝐹𝑙𝑙𝑙𝑙 and 𝐹𝐹𝑟𝑟𝑟𝑟 and associated phase
angles, which serve as one of the primary measures of the imbalance of the part and
are the basis of establishing the mass correction required to balance the part.
Presumably, these calculations are done internally within the spin-balance machine.
If the two correction planes are located at the same location as the two pickup
11
locations, as expected the pickup forces 𝑅𝑅�𝑙𝑙 and 𝑅𝑅�𝑟𝑟 are opposite in sign to the forces
due to the mass imbalance 𝐹𝐹𝑙𝑙𝑙𝑙 and 𝐹𝐹𝑟𝑟𝑟𝑟 .
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑟𝑟𝑟𝑟 ) (𝐿𝐿 − 𝐿𝐿𝑟𝑟𝑟𝑟 )
𝐹𝐹𝑙𝑙𝑙𝑙 cos(∅𝑙𝑙𝑙𝑙 ) = 𝑅𝑅�𝑟𝑟 � 𝑟𝑟 � +𝑅𝑅�𝑙𝑙 𝑙𝑙
cos�∅ � 𝑙𝑙 �
cos�∅ (26)
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 ) (𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑟𝑟𝑟𝑟 ) (𝐿𝐿 − 𝐿𝐿𝑟𝑟𝑟𝑟 )
𝐹𝐹𝑙𝑙𝑙𝑙 sin(∅𝑙𝑙𝑙𝑙 ) = 𝑅𝑅�𝑟𝑟 � 𝑟𝑟 � +𝑅𝑅�𝑙𝑙 𝑙𝑙
sin�∅ � 𝑙𝑙 �
sin�∅ (27)
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 ) (𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 ) (𝐿𝐿 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
𝐹𝐹𝑟𝑟𝑟𝑟 cos(∅𝑟𝑟𝑟𝑟 ) = −𝑅𝑅�𝑟𝑟 � 𝑟𝑟 � −𝑅𝑅�𝑙𝑙 𝑙𝑙
cos�∅ � 𝑙𝑙 �
cos�∅ (28)
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 ) (𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 ) (𝐿𝐿 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
𝐹𝐹𝑟𝑟𝑟𝑟 sin(∅𝑟𝑟𝑟𝑟 ) = −𝑅𝑅�𝑟𝑟 � 𝑟𝑟 � −𝑅𝑅�𝑙𝑙 𝑙𝑙
sin�∅ � 𝑙𝑙 �
sin�∅ (29)
(𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 ) (𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝑙𝑙𝑙𝑙 )
As before, the magnitude of the imbalance forces 𝐹𝐹𝑙𝑙𝑙𝑙 and 𝐹𝐹𝑟𝑟𝑟𝑟 is determined from
the root sum squares of Eqs. 26 and 27 and Eqs. 28 and 29, respectively, while the
phase angles of the imbalance ∅𝑙𝑙𝑙𝑙 and ∅𝑟𝑟𝑟𝑟 are determined from the appropriate
application of the arctangent function using Eqs. 26 and 27 or Eqs. 28 and 29,
respectively.
The analysis in this section presents the equivalent representation of the mass
imbalance of a part due to an arbitrary number of mass imbalances using individual
mass imbalances located two correction planes. Note that the imbalances presented
here represent the actual mass imbalance of the part, not the correction required to
balance the part. In the analysis here, the mass imbalances are considered positive
for mass addition and negative for mass removal (or mass void). The imbalance of
the part is characterized in terms of the imbalance masses, 𝑚𝑚𝑙𝑙𝑙𝑙 and 𝑚𝑚𝑟𝑟𝑟𝑟 at the left
and right correction planes, the imbalance phase angles ∅𝑙𝑙𝑙𝑙 and ∅𝑟𝑟𝑟𝑟 , and the axial
and radial location of the corrections, 𝐿𝐿𝑙𝑙𝑙𝑙 and 𝐿𝐿𝑟𝑟𝑟𝑟 and 𝑟𝑟𝑙𝑙𝑙𝑙 and 𝑟𝑟𝑟𝑟𝑟𝑟 . The imbalance
𝑈𝑈𝑙𝑙𝑙𝑙 and 𝑈𝑈𝑟𝑟𝑟𝑟 and imbalance mass 𝑚𝑚𝑙𝑙𝑙𝑙 and 𝑚𝑚𝑟𝑟𝑟𝑟 are determined from Eqs. 16 and 17
using the imbalance forces 𝐹𝐹𝑙𝑙𝑙𝑙 and 𝐹𝐹𝑟𝑟𝑟𝑟 determined from the pickup forces.
The spin-balance machine typically outputs the mass correction rather than the
actual mass imbalance. There are two modes for the correction for the spin-balance
machine, either mass addition or mass removal. To characterize the mass imbalance
properly, if spin-balance machine is used in mass addition mode, the output
spin-balance corrections correspond to a negative mass imbalance at the same
locations or a positive mass imbalance located at 180° to the imbalance phase
angles ∅𝑙𝑙𝑙𝑙 and ∅𝑟𝑟𝑟𝑟 .
12
3.3 Equations of Motion of Dynamic Imbalance due to a Static
and Couple Imbalance
�𝑀𝑀𝑍𝑍 |𝐿𝐿𝑙𝑙 = −𝑃𝑃(𝐿𝐿𝑃𝑃 − 𝐿𝐿𝑙𝑙 ) − 𝑊𝑊(𝐿𝐿𝐶𝐶𝐶𝐶 − 𝐿𝐿𝑙𝑙 ) + 𝑅𝑅𝑟𝑟 (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) + 𝐹𝐹𝑠𝑠 (𝐿𝐿𝐶𝐶𝐶𝐶 − 𝐿𝐿𝑙𝑙 ) cos(∅𝑆𝑆 + ∅)
(31)
+𝑀𝑀𝐶𝐶 sin(∅𝑐𝑐 + ∅) = 0
Here, 𝐹𝐹𝑆𝑆 is the centrifugal force produced by the static imbalance 𝑈𝑈𝑆𝑆 . The static
imbalance 𝑈𝑈𝑆𝑆 can be represented by a mass imbalance 𝑚𝑚𝑆𝑆 located at a radial
distance 𝑟𝑟𝑆𝑆 from the geometric axis. The static imbalance is located at an angle ∅𝑆𝑆
with respect to the y-axis.
𝐹𝐹𝑠𝑠 = 𝜔𝜔2 𝑈𝑈𝑠𝑠 = 𝜔𝜔2 𝑚𝑚𝑠𝑠 𝑟𝑟𝑠𝑠 (32)
𝑀𝑀𝐶𝐶 is the couple produced by the couple imbalance 𝑈𝑈𝐶𝐶 . The vector direction of the
couple imbalance is oriented at an angle phase angles ∅𝐶𝐶 with respect to the rotation
angle of the part ∅. The couple imbalance 𝑈𝑈𝐶𝐶 can be represented by two masses 𝑚𝑚𝐶𝐶
separated by an axial distance 𝐿𝐿𝐶𝐶 located at opposing circumferential locations at a
radial distance 𝑟𝑟𝐶𝐶 across the longitudinal axis (Eq. 33). Note that the vector
direction of the couple is normal to the plane formed by the longitudinal axis of the
part and the two masses. Other representations of the couple are possible as long as
the same couple imbalance is generated. The axial location of the couple along the
geometric axis is arbitrary.
𝑀𝑀𝐶𝐶 = 𝜔𝜔2 𝑈𝑈𝐶𝐶 = 𝜔𝜔2 𝑚𝑚𝐶𝐶 𝑟𝑟𝐶𝐶 𝐿𝐿𝐶𝐶 (33)
The pickup forces due to mass imbalance at the static and couple imbalance is
shown in Eqs. 34 and 35.
13
(𝐿𝐿𝑃𝑃 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝐶𝐶𝐶𝐶 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝐶𝐶𝐶𝐶 − 𝐿𝐿𝑙𝑙 )
𝑅𝑅𝑙𝑙 = 𝑃𝑃 �1 − � + 𝑊𝑊 �1 − � − 𝐹𝐹𝑆𝑆 �1 − � cos(∅𝑆𝑆 + ∅)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 ) (𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
1
(34)
+ 𝑀𝑀𝐶𝐶 sin(∅𝐶𝐶 + ∅)
(𝐿𝐿𝑟𝑟 − 𝐿𝐿𝑙𝑙 )
Comparing Eq. 1 with Eq. 30, the force due to the static imbalance can be related
to the forces due to an arbitrary number of mass asymmetries within the part:
𝑖𝑖=𝑁𝑁
𝑖𝑖=𝑁𝑁
Using Eqs. 36 and 37, the relation between the static imbalance and the imbalance
from an arbitrary number of mass asymmetries can be obtained:
𝑖𝑖=𝑁𝑁
𝑖𝑖=𝑁𝑁
Similarly, by comparing Eqs. 2 and 31, the relation between the couple imbalance
and the moment about the center of gravity due to an arbitrary number of mass
imbalances within the part can be found:
𝑖𝑖=𝑁𝑁
𝑖𝑖=𝑁𝑁
14
𝑖𝑖=𝑁𝑁
The magnitude of 𝐹𝐹𝑆𝑆 , 𝑈𝑈𝑆𝑆 , 𝑀𝑀𝐶𝐶 , and 𝑈𝑈𝑐𝑐 can be determined from the root sum squares
of Eqs. 36–43, and the phase angles ∅𝑆𝑆 and ∅𝐶𝐶 are determined from the appropriate
application of the arctangent function using Eqs. 36–43.
or equivalently,
𝑈𝑈𝑆𝑆 cos(∅𝑆𝑆 ) = 𝑈𝑈𝑙𝑙𝑙𝑙 cos(∅𝑙𝑙𝑙𝑙 ) + 𝑈𝑈𝑟𝑟𝑟𝑟 cos(∅𝑟𝑟𝑟𝑟 ) (46)
The magnitude of the static imbalance force 𝐹𝐹𝑠𝑠 and the static imbalance 𝑈𝑈𝑠𝑠 is
determined from the root sum squares of Eqs. 44 and 45 or Eqs. 46 and 47,
respectively, and the phase angle of the static imbalance ∅𝑠𝑠 is determined from the
appropriate application of the arctangent function using Eqs. 44 and 45 or 46 and
47. The phase angle of the static imbalance ∅𝑠𝑠 essentially provides the orientation
of the center of gravity of the unbalanced part relative to the coordinate systems.
Similarly, comparing Eqs. 15 and 31, the couple due to the couple imbalance can
be related to the forces due to the imbalance at the two correction planes:
𝑀𝑀𝐶𝐶 cos(∅𝐶𝐶 ) = −𝐹𝐹𝑙𝑙𝑙𝑙 (𝐿𝐿𝑙𝑙𝑙𝑙 − 𝐿𝐿𝐶𝐶𝐶𝐶 )sin(∅𝑙𝑙𝑙𝑙 ) − 𝐹𝐹𝑟𝑟𝑟𝑟 (𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝐶𝐶𝐶𝐶 )sin(∅𝑟𝑟𝑟𝑟 ) (48)
𝑀𝑀𝐶𝐶 sin(∅𝐶𝐶 ) = 𝐹𝐹𝑙𝑙𝑙𝑙 (𝐿𝐿𝑙𝑙𝑙𝑙 − 𝐿𝐿𝐶𝐶𝐶𝐶 )cos(∅𝑙𝑙𝑙𝑙 ) + 𝐹𝐹𝑟𝑟𝑟𝑟 (𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝐶𝐶𝐶𝐶 )cos(∅𝑟𝑟𝑟𝑟 ) (49)
or equivalently,
𝑈𝑈𝐶𝐶 cos(∅𝐶𝐶 ) = −𝑈𝑈𝑙𝑙𝑙𝑙 (𝐿𝐿𝑙𝑙𝑙𝑙 − 𝐿𝐿𝐶𝐶𝐶𝐶 )sin(∅𝑙𝑙𝑙𝑙 ) − 𝑈𝑈𝑟𝑟𝑟𝑟 (𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝐶𝐶𝐶𝐶 )sin(∅𝑟𝑟𝑟𝑟 ) (50)
𝑈𝑈𝐶𝐶 sin(∅𝐶𝐶 ) = 𝑈𝑈𝑙𝑙𝑙𝑙 (𝐿𝐿𝑙𝑙𝑙𝑙 − 𝐿𝐿𝐶𝐶𝐶𝐶 )cos(∅𝑙𝑙𝑙𝑙 ) + 𝑈𝑈𝑟𝑟𝑟𝑟 (𝐿𝐿𝑟𝑟𝑟𝑟 − 𝐿𝐿𝐶𝐶𝐶𝐶 )cos(∅𝑟𝑟𝑟𝑟 ) (51)
The magnitude of 𝑀𝑀𝐶𝐶 and 𝑈𝑈𝑐𝑐 can be determined from the root sum squares of
Eqs. 48–51, and the phase angles ∅𝐶𝐶 is determined from the appropriate application
of the arctangent function using Eqs. 48–51.
15
As stated in Section 3.2, the spin-balance machine typically outputs the mass
correction rather than the actual mass imbalance. If the spin-balance machine is
used in mass addition mode, the actual imbalance angles at the two correction
planes used to compute the static and couple imbalance angles as shown is out of
phase is 180° (𝜋𝜋 radians) with the measured correction angles from the spin-balance
machine.
Section 3.3 provides a means of determining the static and couple imbalance from
the spin-balance machine measurements of the imbalances at the two correction
planes. This formulation assumes that the static imbalance is located at the center
of gravity and the remaining imbalance is attributable to the couple imbalance. The
resulting couple imbalance can be applied at any axial location along the part.
ARL’s spin-balance machine also provides a static and couple balance computation
in one of the output tabulations that essentially breaks down the imbalance at the
two correction planes into a static imbalance component and a couple imbalance
component. Effectively, this locates the static imbalance at the midpoint between
the two correction planes (which is not necessarily located at the center of gravity
of the part). The remaining imbalance of the part is represented as a couple
imbalance, which is a function of the location of the static imbalance. Thus, the
couple imbalance from the spin-balance machine’s tabulated output will not be
equivalent to the couple imbalance computed using Eqs. 50 and 51 unless the center
of gravity is located at the midpoint between the two correction planes. Since
ballistic applications assume the static imbalance is located at the same axial
location as the center of gravity, it is recommended that for ballistics applications,
the mass imbalance at the two correction planes be used to directly compute the
static and couple imbalance (Eqs. 46, 47, 50, and 51) rather than the static couple
imbalance output from the spin-balance machine.
16
4. Effect of Mass Imbalance on Projectile Motion
McCoy3 presents derivations for both the lateral throw-off and aerodynamic jump
due to mass asymmetry. The net deviation (or jump) due to both effects is additive,
although each component can differ in direction. McCoy’s ballistic coordinate
system is used here. This system orients the x-axis in the downrange direction, the
y-axis oriented up, and the z-axis oriented to the right when looking downrange, as
shown in Fig. 6. For simplicity, the nominal line of fire is oriented along the x-axis.
17
There is also a second body-fixed coordinate system in which the physical
properties (for example, moments of inertia and mass asymmetries) of the projectile
are characterized. Note the x- and z-axes of this coordinate system are oriented in
different directions than the coordinate system used in spin-balance measurements
shown in Fig. 3. The changes in the orientations have been accounted for in the
analysis that follows.
In the following analysis, the most significant difference in the body-fixed and
spin-balance coordinate systems is with respect to the circumferential angle
measurements in both coordinate systems, as shown in Fig. 7. Effectively, the
circumferential angle measurement involves a sign change when going from one
coordinate system to another. It is assumed that the circumferential angles from
spin-balance measurements are used as the actual measured values, and the required
sign change is incorporated in the equations in Section 4.1.
18
projectile leaves the gun tube. For right-hand gun-tube twist, if the offset in the
center of gravity produced by the static imbalance is oriented upward in the vertical
plane, and the deflection due to throw-off will be to the right when looking
downrange. The deflection can be computed as shown in Eq. 52. Here, n is the twist
of the rifling (calibers/turn: the number of body diameters traveled in one complete
revolution of the projectile [positive for right-hand twist]), ϵ is the magnitude of the
lateral displacement of the center of gravity, D is the projectile diameter, and 𝛷𝛷𝑆𝑆0
is the circumferential orientation of the center of gravity at muzzle exit relative to
the ballistic coordinate frame. The throw-off produces a vertical velocity 𝑉𝑉𝑌𝑌0
(positive upward) and/or a lateral velocity 𝑉𝑉𝑍𝑍0 (positive to the gunner’s right) at the
muzzle that results in an angular deflection off the nominal trajectory when scaled
by the forward velocity 𝑉𝑉𝑋𝑋0 of the projectile.
19
projectile is a random event and the throw-off can be computed from the first term
alone:
𝐿𝐿𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵
Φ𝑆𝑆0 = (−∅𝑆𝑆 ) + (∅𝑐𝑐ℎ𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 + 2𝜋𝜋𝜋𝜋 𝐷𝐷
) (54)
Here 𝐶𝐶𝐿𝐿𝛼𝛼 and 𝐶𝐶𝑀𝑀𝛼𝛼 are the aerodynamic lift coefficient and pitching moment
coefficient slopes respectively, 𝜉𝜉0 and 𝜉𝜉0′ are the complex yaw and yawing rate at
the muzzle and 𝑘𝑘𝑡𝑡 and 𝑃𝑃 are defined by Eqs. 56 and 57:
𝐼𝐼𝑡𝑡
𝑘𝑘𝑡𝑡2 = (56)
𝑚𝑚𝐷𝐷2
𝐼𝐼𝑥𝑥 𝑝𝑝𝑝𝑝
𝑃𝑃 = � � (57)
𝐼𝐼𝑡𝑡 𝑉𝑉
𝑝𝑝𝑝𝑝
𝐼𝐼𝑥𝑥 and 𝐼𝐼𝑡𝑡 are the axial and transverse moments of inertia, � 𝑉𝑉 � is the nondimensional
spin rate of the projectile, 𝐷𝐷 is the projectile reference diameter, 𝑚𝑚 is the projectile
mass, 𝑝𝑝 is the projectile spin rate, and 𝑉𝑉 is the magnitude of the projectile velocity.
The aerodynamic jump shown in Eq. 55 includes the classical aerodynamic jump
for a symmetric projectile as well as the additional term (shown as the last term)
20
that accounts for the jump due to the mass asymmetry. Here, 𝑛𝑛 is the twist of the
rifling (calibers/turn: the number of body diameters traveled for one complete
revolution of the projectile [positive for right-hand twist]), 𝐾𝐾𝑇𝑇 is amplitude of the
yawing motion produced by the couple imbalance, and 𝜙𝜙0 is the orientation of the
yaw arm produced by the asymmetry at muzzle exit relative to the vertical
coordinate, with positive angles determined from the right-hand rule in the ballistic
coordinate frame.
Due to the imbalance, the projectile will rotate about its longitudinal principal axis,
producing an additional component to the yawing motion. The amplitude of this
additional component of yawing motion can be determined from the principal axis
tilt produced by the couple imbalance, as shown in Eq. 58. (Additional details can
be found in the Appendix.) Here, 𝑈𝑈𝐶𝐶 is the magnitude of the couple imbalance.
Representation of the principal axis tilt in terms of the couple imbalance represents
a generalization of the derivation by McCoy3 for the dynamic imbalance produced
by a single small mass removal:
𝑈𝑈𝐶𝐶
𝐾𝐾𝑇𝑇 = (58)
𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥
The aerodynamic jump due to the couple imbalance 𝑈𝑈𝐶𝐶 is shown in Eq. 59:
Similar to the lateral throw-off, the aerodynamic jump is a function of the angular
orientation of the couple imbalance when the projectile exits the muzzle, as
presented in Eq. 60. The first term in Eq. 60 (−∅𝐶𝐶 ) accounts for the correct vector
orientation of the couple imbalance relative the ballistic coordinate frame using the
orientation of the couple imbalance ∅𝐶𝐶 as measured in the spin-balance fixture
coordinate frame. The second term accounts for the initial orientation of the
projectile when inserted in the chamber relative to its position in the spin-balance
fixture ∅𝑐𝑐ℎ𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 as well as the rotation of the projectile as it travels down the barrel
length 𝐿𝐿𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵 . If the contribution of dynamic imbalance to the dispersion of a group
of bullets is required to be determined, the initial orientation of the projectile is a
random event and the aerodynamic jump can be computed from the first term alone.
𝐿𝐿𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵𝐵
𝜙𝜙0 = (−∅𝐶𝐶 + 𝜋𝜋/2) + (∅𝑐𝑐ℎ𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 + 2𝜋𝜋𝜋𝜋 𝐷𝐷
) (60)
While the orientation of the projectile in the weapon chamber can often be
considered a random variable, mass imbalance generally produces both a static and
21
dynamic imbalance, and the amplitude and phases angles associated with both types
of mass imbalance are directly related and should not be treated as independent
random variables relative to each other.
22
are shown in Table 2. Conforming to the spin-balance coordinate system with the
nose oriented to the left, the left correction planes were located at the cone cylinder
junction, and the right correction plane was located at the center of gravity of the
mass imbalance. The computed imbalance of the part is shown in the upper half of
Table 2 as obtained from appropriate application of Eqs. 16, 17, and 20–25.
Mass
Mass imbalance Mass imbalance Mass imbalance
Circumferential imbalance
at left correction phase angle at at right
location of mass phase angle at
plane 𝑼𝑼𝒍𝒍𝒍𝒍 left correction correction plane
offset, φi (deg) right correction
(mg–cm) plane ∅𝒍𝒍𝒍𝒍 (deg) 𝑼𝑼𝒓𝒓𝒓𝒓 (mg-cm)
plane ∅𝒓𝒓𝒓𝒓 (deg)
180° (down) 0 0 1583 0
270° (gunner’s right) 0 0 1583 90
0° (up) 0 0 1583 180
90° (gunner’s left) 0 0 1583 –90
Circumferential Mass Mass correction Mass correction Mass correction
location of mass correction phase angle at left (addition) at right phase angle at
offset, φi (deg) (addition) at left correction plane correction plane right correction
correction plane (deg) (mg) plane (deg)
(mg)
180ׄ° (down) 0 0 1586 180
270° (gunner’s right) 0 0 1586 270
0° (up) 0 0 1586 0
90° (gunner’s left) 0 0 1586 90
23
correction plane because the center of gravity of the imbalance is centered over the
right correction plane.
Once the correction plane imbalance magnitudes and phase angles are obtained
from the spin-balance mass correction and phase angles, the static and couple
imbalance can be computed using Eqs. 46, 47, 50, and 51. The results are shown in
Table 3.
Static
Static Couple Couple
Circumferential imbalance
imbalance imbalance imbalance
location of mass offset, phase
𝑼𝑼𝒔𝒔 𝑼𝑼𝒄𝒄 phase angle
φi (deg) angle ∅𝑺𝑺
(mg-cm) (mg– m2) ∅𝑪𝑪 (deg)
(deg)
180° (down) 1583 0 4331 90
270° (gunner’s right) 1583 90 4331 180
0° (up) 1583 180 4331 –90
90° (gunner’s left) 1583 270 4331 0
Using the static and couple imbalance, the lateral throw-off and aerodynamic jump
were computed, as shown in Table 4. The results are virtually identical to McCoy’s
results. The results also show the relative differences between the more accurate
computation of aerodynamic jump (Eq. 59) and the simplified computation
(Eq. 61). The difference in the two approach is about 10% for this example.
Lateral
Aero jump Aerodynamic
Circumferential Lateral throw-off
Aerodynamic direction jump (mil) –
location of mass throw- direction
jump relative to simplified
offset, φi off relative to
(mil) ballistic computation,
(deg) (mil) ballistic
coordinates Eq. 61
coordinates
180° (down) 1.12 Right 1.25 Right 1.14
270° (gunner’s right) 1.12 Up 1.25 Up 1.14
0° (up) 1.12 Left 1.25 Left 1.14
90° (gunner’s left) 1.12 Down 1.25 Down 1.14
24
6. Results from Imbalance Measurements of Calibration Pins
Calibration pins were constructed to verify the calibration and accuracy of the
spin-balance machine. Because of the simple geometry of the part, theoretical
calculations of the imbalance were also obtained and compared with the
spin-balance measurements. Each of the pins were cylindrical in shape with mass
imbalances drilled into the pins at fixed locations. Figure 9 shows a drawing of one
of the pins (Pin C). The same hole dimensions were used of each of the pins. With
the exception of Pin B, which had a hole drilled at the center, the remaining pins
had holes drilled at the same outboard locations.
Figure 10 shows a drawing of each of the four calibration pins used in the study.
Pins A and C had asymmetries drilled at the outboard end of the pins. These pins
produced both static and dynamic imbalances. (Pin A had two additional holes
drilled that opposed each other and did not contribute to the imbalance.) Pin B had
a hole drilled in the middle of the pin, yielding only a static imbalance. Pin D had
two holes drilled at opposing ends of the pin and opposite circumference angles,
producing a pure couple imbalance with no static imbalance.
25
Fig. 10 Calibration pins
Each of the pins were weighed before and after the holes were drilled. This allowed
the actual weight of the asymmetry to be measured. The weight measurement of
the cylindrical pins before being drilled also allowed the density of the material to
be computed. The drilled holes were not perfectly cylindrical because of the curved
top produced by the cylindrical shape of the pin itself. The actual volume of the
drilled hole was computed numerically and found to be about 5% less than a right
cylinder 0.1 inch long. Theoretical predictions of the mass of the hole were based
on the measured density and the nominal dimensions from the drawing. The
measured mass of the asymmetry was on average 4% higher than the theoretical
value. This was assumed to be due to machining tolerances. To account for the mass
difference, the hole depth in the theoretical calculations was increased, as shown in
Table 5. The increased depth of the hole decreased the theoretical location of the
center of gravity of the mass asymmetry by less than 2% of the nominal location.
This allowed the theoretical mass of the asymmetry to be equal to the measured
mass.
Table 5 Physical properties of mass asymmetry for the four calibration pins
Asymmetry CG
Measured weight of Hole depth
Pin from pin
asymmetry (mg) (cm/inches)
centerline (cm)
Nominal 161.5 (theoretical) 0.254 / 0.1 0.341
A 170 0.2667 / 0.105 0.335
B 168 0.2637 / 0.1038 0.3365
C 167 0.2621 / 0.1032 0.3373
D 166.5 0.2616 / 0.103 0.3376
26
Table 6 presents the measured and theoretical mass corrections at the left and right
correction planes. The correction was assumed to be made at the outer diameter of
the pins. Note that the location of the correction planes was placed at same location
as the pickup locations, so the mass correction would be proportional to the
amplitude of the fluctuating part of the pickup forces. The theoretical imbalance at
the left and right correction plane was out of phase with each other by 180° except
for Pin B, where the imbalances at the left and right correction plane were
in-phase. The measured and theoretical phase angle difference between the left and
right correction plane differed by less than a few degrees. The measured values of
the imbalance for Pin B at both correction planes shows a deviation from the
expected symmetry that should exist for this configuration.
Table 6 Measured and theoretical mass corrections at left and right correction planes for
the four pins
Measured Theoretical
Measured left Theoretical left
right mass right mass
Pin mass correction mass correction
correction correction
(mg) (mg)
(mg) (mg)
A 6.5 30.0 149.0 149.9
A 6.6 30.0 148.0 149.9
B 72.9 59.5 59.2 59.5
B 73.0 59.5 59.4 59.5
C 153.2 148.2 8.5 29.6
C 155.7 148.2 8.4 29.6
D 155.8 177.5 156.9 177.5
D 156.7 177.5 156.1 177.5
Average before
0.8 0 1.0 0
drilling holes
Pins A and C essentially had the same mass asymmetries, but the asymmetries were
located at the opposite ends of the pins. Both pins should show the same measured
imbalance when swapped end-for-end or when compared with each other in the
same orientation. Table 7 presents the measured and theoretical mass imbalances
at the left and right correction planes for Pins A and C with the asymmetries in their
nominal location (right and left, respectively, for Pin A and Pin C) and with the
asymmetries reversed (left and right, respectively for Pin A and Pin C). Overall, the
measured mass imbalance yielded similar results for Pin A and C irrespective of
the orientation of the pin when it was measured. (Note that Pins A and C are slightly
different because of the slight mass differences in the mass of material removed
when drilled. These differences are accounted for in the theoretical predictions.)
The differences between the measured and theoretical mass imbalances are most
apparent when the imbalance is small.
27
Table 7 Measured and theoretical mass corrections at left and right correction planes
Table 8 Measured and theoretical static and couple imbalances for the four calibration pins
Measured Theoretical
Measured static Theoretical static
couple couple
Pin imbalance imbalance
imbalance imbalance
(mg–cm) (mg–cm)
(mg–cm2) (mg–cm2)
A 67.68 56.95 75.05 86.79
A 67.17 56.95 74.60 86.79
B 62.71 56.53 6.64 0.00
B 62.88 56.53 6.60 0.00
C 68.77 56.33 77.99 85.85
C 70.01 56.33 79.15 85.85
D 1.40 0.00 150.90 171.33
D 1.33 0.00 150.95 171.33
Average before
drilling holes 0.74 0 0.4 0
Similar results were found for Pins A and C when the pins were swapped end-for-
end, as shown in Table 9. This set of test results were obtained on a different
occasion than the results in Table 8. The measured results between pins occasion
differences in the results are seen when the Table 9 results are compared with the
results in Table 8.
28
Table 9 Comparison of measured and theoretical values of static and couple imbalance
Measured Theoretical
Measured static Theoretical static
couple couple
Pin imbalance imbalance
imbalance imbalance
(mg–cm) (mg–cm)
(mg–cm2) (mg–cm2)
A (Nominal) 62.13 56.95 76.36 86.79
A (Nominal) 62.68 56.95 76.48 86.79
A (Reversed) 64.49 56.95 81.18 86.79
A (Reversed) 64.87 56.95 81.18 86.79
C (Nominal) 63.56 56.33 81.17 85.85
C (Nominal) 64.13 56.33 80.68 85.85
C (Reversed) 63.34 56.33 74.53 85.85
C (Reversed) 63.34 56.33 74.44 85.85
Overall, the calibration pins provided some measure of validation of the results
obtained with the spin-balance machine’s ability to measure the static and couple
imbalance, although there is some indication that improvements in the calibration
of the pickup sensors may be desired.
7. Conclusion
Spin-balance machines typically measure the applied loads at two locations that
support the rotating part. The fluctuating component of the applied loads can be
directly related to the net effect of the centrifugal forces produced by the individual
mass asymmetries within the part. The analysis presented here shows that a
projectile with an arbitrary number of mass imbalances can be represented by two
strategically located mass imbalances or by a static imbalance and a couple
imbalance. The relationship between each of the three types of imbalance is
presented. Representation of the imbalance as mass imbalances at two locations
allows for a straightforward means of correcting the imbalance.
The imbalance can be directly related to the deflection of the trajectory of the
projectile due to the imbalance. While the spin-balance machine used by ARL
typically represents the imbalance by mass imbalance at two correction planes, the
trajectory deflection is more easily determined from the static and couple
imbalances although the static and couple imbalances are easily determined from
the two mass imbalances using the analysis shown here.
29
8. References
3. McCoy RL. Modern exterior ballistics. Atglen (PA): Schiffer Military History;
1999.
30
Appendix. Derivation of Aerodynamic Jump from Couple
Imbalance
31
A pure dynamic imbalance will produce a contribution to the yawing motion of the
projectile and a deflection to the trajectory due to aerodynamic jump from this
induced yawing motion. The dynamic imbalance causes the projectile to spin about
its longitudinal principal axis rather than its geometric longitudinal axis, producing
an angle of attack and subsequent yawing motion. The combination of yawing
motion due to the dynamic imbalance along with the normal epicyclic yawing
motion of a symmetric projectile produces what is referred to as a tricyclic yawing
motion.
To determine the aerodynamic jump due to the dynamic imbalance, the principle
axis tilt relative to the geometric axis must be determined. The moments of inertia
in the body-fixed (geometric) axis frame can be represented by the inertia matrix in
Eq. A-1. For a symmetric projectile, the diagonal terms consist of the axial moment
of inertia ( 𝐼𝐼𝑥𝑥𝑥𝑥 = 𝐼𝐼𝑥𝑥 ) and the transverse moment of inertia (𝐼𝐼𝑦𝑦𝑦𝑦 = 𝐼𝐼𝑧𝑧𝑧𝑧 = 𝐼𝐼𝑡𝑡 ) and the
off-diagonal terms are zero.
32
𝑖𝑖=𝑛𝑛
𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆
𝐼𝐼𝑥𝑥𝑥𝑥 = 𝐼𝐼𝑥𝑥 + � 𝑚𝑚𝑖𝑖 (𝑦𝑦𝑖𝑖2 + 𝑧𝑧𝑖𝑖2 )
𝑖𝑖=1
𝑖𝑖=𝑛𝑛
𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆
𝐼𝐼𝑦𝑦𝑦𝑦 = 𝐼𝐼𝑡𝑡 + � 𝑚𝑚𝑖𝑖 (𝑥𝑥𝑖𝑖2 + 𝑧𝑧𝑖𝑖2 )
𝑖𝑖=1
𝑖𝑖=𝑛𝑛
𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆
𝐼𝐼𝑧𝑧𝑧𝑧 = 𝐼𝐼𝑡𝑡 + � 𝑚𝑚𝑖𝑖 (𝑥𝑥𝑖𝑖2 + 𝑦𝑦𝑖𝑖2 )
𝑖𝑖=1
(A-2)
𝑖𝑖=𝑛𝑛
The principal moments of inertia and the principal directions can be obtained from
the eigenvalues and eigenvectors associated with the inertia matrix. The primary
interest here is the principal direction along the spin axis 𝑃𝑃�⃗1, as this will allow the
principal axis tilt to be determined. The eigenvalue associated with the spin axis is
approximately equal to 𝐼𝐼𝑥𝑥𝑥𝑥 for small mass asymmetries. Computing the eigenvector
associated with this eigenvalue allows the direction of the principal axis to be
determined:
1 𝐼𝐼𝑥𝑥𝑥𝑥 𝐼𝐼𝑥𝑥𝑥𝑥
𝑃𝑃�⃗1 = �𝚤𝚤̃ + � � 𝚥𝚥̃ + � � 𝑘𝑘��
𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥 𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥
𝐼𝐼 2 2 (A-3)
�1 + � 𝑥𝑥𝑥𝑥 � + � 𝐼𝐼𝑥𝑥𝑥𝑥 �
𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥 𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥
As shown in Eq. A-4, the couple imbalance is related to the combined effect of the
individual imbalances as measured in the spin-balance fixture coordinates, as
shown in Eqs. 42 and 43 from Section 3.3 of the main report.
𝑖𝑖=𝑛𝑛
33
The off-diagonal moment of inertia terms 𝐼𝐼𝑥𝑥𝑥𝑥 and 𝐼𝐼𝑥𝑥𝑥𝑥 can be related to the couple
imbalance measured on the spin-balance fixture. The location of the asymmetries
(𝑥𝑥𝑖𝑖 , 𝑦𝑦𝑖𝑖 , 𝑧𝑧𝑖𝑖 ) in body-fixed ballistic coordinate frame can be transformed to their
locations in the spin-balance coordinate frame ( 𝐿𝐿𝑖𝑖 , 𝑟𝑟𝑖𝑖 , ∅𝑖𝑖 ), as seen in Figs. 6 and 7
of the main report.
𝑥𝑥𝑖𝑖 = −(𝐿𝐿𝑖𝑖 − 𝐿𝐿𝐶𝐶𝐶𝐶 )
𝑖𝑖=𝑛𝑛 𝑖𝑖=𝑛𝑛
𝐼𝐼𝑥𝑥𝑥𝑥 = � 𝑚𝑚𝑖𝑖 𝑥𝑥𝑖𝑖 𝑧𝑧𝑖𝑖 = + � 𝑚𝑚𝑖𝑖 (𝐿𝐿𝑖𝑖 −𝐿𝐿𝐶𝐶𝐶𝐶 )𝑟𝑟𝑖𝑖 sin( ∅𝑖𝑖 ) = −𝑈𝑈𝐶𝐶 cos(∅𝐶𝐶 ) = −𝑈𝑈𝐶𝐶 cos(−∅𝐶𝐶 )
𝑖𝑖=1 𝑖𝑖=1
Considering just the component of the yawing motion due to the dynamic
imbalance in the absence of the epicyclic yawing motion, the freestream velocity
vector would be aligned with the principal axis associated with the spin axis 𝑃𝑃�⃗1 .
�⃗
𝑉𝑉 𝑉𝑉𝑥𝑥 𝑉𝑉𝑦𝑦 𝑉𝑉𝑧𝑧
= 𝚤𝚤̃ + 𝚥𝚥̃ + 𝑘𝑘� = −𝑃𝑃������⃗1
�⃗ � �𝑉𝑉
�𝑉𝑉 �⃗ � �⃗ �
�𝑉𝑉 �⃗ �
�𝑉𝑉
−1 𝐼𝐼𝑥𝑥𝑥𝑥
= �𝚤𝚤̃ + � � 𝚥𝚥̃
2 2
𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥
𝐼𝐼
�1 + � 𝑥𝑥𝑥𝑥 � + � 𝐼𝐼𝑥𝑥𝑥𝑥 �
𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥 𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥
𝐼𝐼𝑥𝑥𝑥𝑥
+� � 𝑘𝑘�� (A-6)
𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥
−1 𝑈𝑈𝐶𝐶 sin(−∅𝐶𝐶 )
= �𝚤𝚤̃ + � � 𝚥𝚥̃
2 2
𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥
𝐼𝐼
�1 + � 𝑥𝑥𝑥𝑥 � + � 𝐼𝐼𝑥𝑥𝑥𝑥 �
𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥 𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥
−𝑈𝑈𝐶𝐶 cos(−∅𝐶𝐶 )
+� � 𝑘𝑘��
𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥
For small asymmetries, the resulting angle of attack components are small and the
complex yaw can be written as shown in Eq. A-7.
𝑉𝑉𝑌𝑌 + 𝑖𝑖𝑉𝑉𝑧𝑧
𝜉𝜉𝑇𝑇 = = 𝛼𝛼𝑇𝑇 + 𝑖𝑖𝛽𝛽𝑇𝑇 = 𝐾𝐾𝑇𝑇 𝑒𝑒 𝑖𝑖𝜙𝜙𝑇𝑇
�⃗
�𝑉𝑉 �
(A-7)
−𝑈𝑈𝐶𝐶
=� � (sin(−∅𝐶𝐶 ) − 𝑖𝑖 cos(−∅𝐶𝐶 ))
𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥
34
From Eq. A-7, the magnitude of the yaw 𝐾𝐾𝑇𝑇 due to the couple imbalance is obtained
and shown in Eq. A-8, and the relation between the orientation of the yaw arm in
the ballistic frame to the orientation of the vector orientation of the couple in the
spin-balance frame resulting from the dynamic imbalance is shown in Eq. A-9.
𝑈𝑈𝐶𝐶
𝐾𝐾𝑇𝑇 = � � (A-8)
𝐼𝐼𝑡𝑡 − 𝐼𝐼𝑥𝑥
∅𝐶𝐶 = 𝜙𝜙0 − 𝜋𝜋/2 (A-9)
McCoy 1 provided a derivation of the aerodynamic jump of a projectile with a mass
asymmetry, as shown in Eq. A-10. Similar to the lateral throw-off, the aerodynamic
jump produces a vertical velocity 𝑉𝑉𝑌𝑌0 (positive upward) and/or a lateral velocity 𝑉𝑉𝑍𝑍0
(positive to the gunner’s right) at the muzzle that results in an angular deflection
off the nominal trajectory when scaled by the forward velocity 𝑉𝑉𝑋𝑋0 of the projectile.
The aerodynamic jump includes the classical aerodynamic jump for a symmetric
projectile as well as the additional term (shown as the last term of Eq. A-10), which
accounts for the jump due to the mass asymmetry:
𝐼𝐼𝑡𝑡
𝑘𝑘𝑡𝑡2 = (A-11)
𝑚𝑚𝐷𝐷2
𝐼𝐼𝑥𝑥 𝑝𝑝𝑝𝑝
𝑃𝑃 = � � (A-12)
𝐼𝐼𝑡𝑡 𝑉𝑉
𝑝𝑝𝑝𝑝
𝐼𝐼𝑥𝑥 and 𝐼𝐼𝑡𝑡 are the axial and transverse moments of inertia, � 𝑉𝑉 � is the nondimensional
spin rate of the projectile, 𝐷𝐷 is the projectile reference diameter, 𝑚𝑚 is the projectile
mass, 𝑝𝑝 is the projectile spin rate, and 𝑉𝑉 is the magnitude of the projectile velocity.
The last term in Eq. A-10 represents the aerodynamic jump contribution due to the
mass asymmetry. Here, n is the twist of the rifling (calibers/turn: the number of
body diameters traveled for one complete revolution of the projectile [positive for
right-hand twist]), 𝐾𝐾𝑇𝑇 is amplitude of the yawing motion produced by the couple
imbalance, and 𝜙𝜙0 is the orientation of the yaw arm produced by the asymmetry at
1
McCoy RL. Modern exterior ballistics. Atglen (PA): Schiffer Military History; 1999.
35
muzzle exit relative to the vertical coordinate with positive angles determined from
the right-hand rule in the ballistic coordinate frame.
Substituting Eqs. A-8 and A-9 into the last term of Eq. A-10, the aerodynamic jump
due to the mass asymmetry is obtained, as shown in Eq. A-13. Here the
aerodynamics jump is related to the magnitude and angular orientation of the couple
imbalance 𝑈𝑈𝐶𝐶 and ∅𝐶𝐶 , respectively, as measured from the spin-balancing procedure.
McCoy 1 presented a derivation of the aerodynamic jump for the dynamic imbalance
produced by a single small mass removal. Eq. A-13 represents the aerodynamic
jump due to an arbitrary distribution of multiple mass asymmetries in terms of the
couple produced by the asymmetries and represents a generalization of the result
presented by McCoy.
1
McCoy RL. Modern exterior ballistics. Atglen (PA): Schiffer Military History; 1999.
36
1 DEFENSE TECHNICAL
(PDF) INFORMATION CTR
DTIC OCA
1 CCDC ARL
(PDF) FCDD RLD CL
TECH LIB
4 CCDC ARL
(PDF) FCDD RLW LE
P WEINACHT
T PUCKETT
G OBERLIN
FCDD RLW LH
J GARNER
37