Eﬀect of Bending Deformation on the Lateral Force of Spinning Projectiles with Large Aspect Ratio

: The bending deformation can aﬀect the lateral force of spinning projectiles with large aspect ratios, thus interfering with their ﬂight stability. Based on the established spin–deformation coupling motion model, the unsteady Reynolds averaged Navier–Stokes (URANS) equations are solved to simulate the ﬂow over a large−aspect−ratio projectile undergoing spin and spin−defor-mation coupling motion by using the dual−time stepping method and dynamic mesh technique, obtaining the lateral force. Furtherly, the ﬂow mechanism is analyzed for the changed lateral force induced by the bending deformation. The results indicate that the variation of transient lateral force for the head of a projectile is consistent with that of the deformation−induced additional sideslip angle; aﬀected by the deformation−induced compression wave and expansion wave, the time−aver-aged lateral force for the middle of a projectile will be increased at small angles of a�ack, but changed li�le at large angles of a�ack; at small angles of a�ack, the change trend of transient lateral force for the tail of a projectile is similar to that of additional angle of a�ack caused by the deformation; at large angles of a�ack, the characteristic of phase lag is presented between the transient lateral force for the tail of a projectile and the additional sideslip angle.


Introduction
The shape of a large aspect ratio with a light structure is usually adopted by supersonic projectiles to increase their flight distance and payload, which consequently leads to the obvious structural deformation.The coupling of spinning, structural deformation and coning motion occurs in the flight process of spinning projectiles with a large aspect ratio, which brings uncertain effects on their maneuverability and stability.With the improvement of control precision and the decrease of stability margin, the influence of the bending deformation on aerodynamic characteristics has been widely considered [1][2][3].Therefore, it is necessary to investigate the aerodynamic characteristics of a projectile undergoing the spin−deformation coupling motion.
Due to the Magnus effect, lateral force is generated for a spinning projectile with angles of a ack [4,5].Although the lateral force is usually 1/100 to 1/10 of the normal force, the corresponding yaw moment always deviates the projectile from the plane of angle of a ack, which interferes with its flight stability [6,7].In 1955, Martin [8] began to theoretically explore the Magnus effect of rotating cylinders.Later, Iversen et al. [9,10] conducted a theoretical analysis of the Magnus effect on the body of revolution.Generally, theoretical research is usually suitable for spinning projectiles with simple shapes.However, the shape of projectiles is often complex practically, which makes it difficult to obtain the aerodynamic characteristics caused by the spinning motion using the theoretical method.With the development of computer technology and numerical methods, the numerical simulation based on computational fluid dynamics (CFD) has been widely used in investigations on the aerodynamic characteristics of spinning projectiles [11][12][13][14][15][16][17][18].For example, Simon et al. [11] used the Spalart-Allmaras (S-A) turbulence model to numerically simulate the flow field for the body of revolution.The results indicated that the Magnus effect can be effectively reduced with the appropriate location of the contra−spinning surface.Kla et al. [12] researched the flow field for the body of revolution with an aspect ratio of 6.37 based on the − turbulence model, and pointed out that the emergence of secondary vortices was a reason for the change in the locally lateral force.Yin et al. [13,14] considered the projectile AFF as the research object to analyze the influence of aerodynamic interference among canards, the body and fins on Magnus effect.Unfortunately, the present studies rarely involve the lateral force for the large−aspect−ratio projectile undergoing the spin−deformation coupling motion.
The complex aerodynamic effects will be induced by the spin−deformation coupling motion.Motivated to analyze the flight stability of the spinning projectile with deformation, the effective angle of a ack is usually employed to represent the deformation.As a result, the aerodynamic coefficient of the deformed projectile is obtained through multiplying the effective angle of a ack by the aerodynamic derivative [19][20][21][22].However, for this method, the effect of deformation on the flow field of the spinning projectile is ignored, while the nonlinear aerodynamic force caused by the deformation cannot be determined [23].Wu et al. [24] studied the aerodynamic characteristics of the body of revolution with a deformed head, and the results showed that the deformation of head affected the normal and lateral forces.Youn et al. [25] and Paul et al. [26] regarded the projectile as a rigid body with several hinges, and then analyzed the influence of the bending of the projectile at the hinged position on the aerodynamic characteristics.All these results demonstrate that it is necessary to research the effect of deformation on the aerodynamic characteristics of spinning projectiles.Unfortunately, the previous works are quite limited.On the one hand, the research object is mainly the body of revolution with a small aspect ratio.On the other hand, the deformation form of the projectile is too simplified to fully study the influence of the deformation on the aerodynamic characteristics.Therefore, it is urgent to systematically research the effect of deformation on the lateral force of spinning projectiles with large aspect ratios.
In order to investigate the influence of deformation on the lateral force of spinning projectiles with large aspect ratios, the spin−deformation coupling motion model is established in this paper.The Apache, a spinning projectile with a large aspect ratio, is chosen as the research object.Based on the unsteady numerical simulation method and dynamic mesh technology, the flow over the large−aspect−ratio projectile undergoing the spin−deformation coupling motion is simulated.Furtherly, the variation of lateral force is obtained, and the flow mechanism is analyzed.The conclusions provide effective guidance for the analysis of the aerodynamic characteristics and flight performance of the deformed spinning projectile with a large aspect ratio.

Computational Model and Grid
As a typical spinning projectile with a large aspect ratio, the Apache is used to explore the lateral force of the deformed projectile.Its model dimension is shown in Figure 1.The projectile diameter d is chosen as a standard caliber.The Apache comprises a 3d ogive nose, a 21.88d cylindrical body and four fins with an installation angle ( = 2 °).
The origin O of the inertial coordinate system, Oxyz, coincides with the initial position of head vertex.The axis Ox is consistent with the longitudinal axis of the projectile, and its positive direction is from the nose to the tail.The axis Oy is located in the vertical plane containing the axis Ox, which is perpendicular to the axis Ox.Its positive direction is upward.The axis Oz and the other two axes form the right−hand coordinate system.Figure 2 shows the spin direction, the fins and the roll angle.The spin is counterclockwise.The leeward fin at the initial time is defined as fin1.The others are defined as fin2, fin3 and fin4 in turn along the spin direction.The angle between the fin1 and the axis Oy is determined as the rolling angle .Figure 3 shows the mesh for the flow calculation.The outer domain and the inner domain constitute the whole computational domain, which are connected by the interface.Numerical interpolation is employed at the interface to ensure the flux conservation between the two domains.The dynamic mesh technique is used to adjust the mesh of the inner domain to adapt to the spin and deformation of the projectile, while the mesh of the outer domain remains unchanged.The boundary conditions in the forward and circumferential direction are freestream condition.The boundary in the projectile base direction is set to pressure−outlet.The projectile surface is set to no−slip and adiabatic wall condition.

Spin−Deformation Coupling Motion
The structural deformation of projectiles includes the bending and the twisting of the body and fins, which is extremely complicated.In this paper, only the deformation of the body is considered, and that of the fins is ignored.In Figure 4, the spin−−deformation coupling motion can be, respectively, decomposed into the spin and deformation.The where ( ) represents the deformation of the projectile and is described as where ( , ) denotes the deformation law of the projectile.By the coordinate transformation, it is obtained as Equations ( 1)-( 3) are the mathematical description for the spin−deformation coupling motion.Generally, the bending shape of the projectile is mainly manifested as its first order mode shape [21].Thus, the separation variable method can be exploited to obtain the deformation law of the projectile as where ( ) is the variation of deformation with time; and ( ) denotes the first order mode shape.To simplify the calculations, the ( ) is wri en as where A is the deformation amplitude.represents the initial phase, ∈ [0, 90 °].F denotes the frequency, which is equal to the spinning frequency of the projectile.Hence, the deformation law of the projectile is expressed as The first order mode shape of the Apache obtained by modal analysis is plo ed in Figure 5. Here, denotes the deformation of the head vertex.Once Equation ( 6) is determined, Equations (1)−(3) can be solved using the numerical method.x/L / head The spin−−deformation coupling motion trajectories are showed in Figure 6.In the figure, the and are the longitudinal and lateral components of the spin−deformation coupling motion, respectively.It is seen that the motion trajectory is an ellipse when Δ = 0 °, while the trajectory is the heart shape for Δ = 90 °.When Δ ∈ (0, 90 °), the trajectory gradually changes from an ellipse to a heart shape.For convenience, the motion trajectory is defined as (Δ ∈ [0, 90 °]).For example, represents the elliptical motion trajectory, while is the heart−shaped motion trajectory.Moreover, when the projectile does not deform (i.e., = 0), the trajectory is defined as .Remarkably, since is axisymmetric about the y axis, the time−averaged value of along the y axis is zero, and the induced time−averaged additional angle of a ack is also zero.Similarly, is axisymmetric about the z axis, so the induced time−averaged additional sideslip angle is also zero.

Numerical Method
For flow simulation, both large eddy simulation (LES) and URANS are available methods.The fundamental idea of LES is to directly simulate larger eddies that significantly contribute to the overall flow behavior, while modeling or filtering out smaller eddies that are dissipative and have less impact on the overall flow characteristics [27].LES can capture fine−scale flow features, thus enhancing its predictive precision for turbulence characteristics.Nonetheless, LES often demands fine grid resolutions, leading to a substantial increase in computational costs and the demand for computing resources.
In contrast, URANS significantly reduces the computational cost by computing time−averaged solutions of the Navier-Stokes equations to predict the mean behavior of fluid flow [27].Moreover, URANS typically necessitates specific grid density near the wall, while maintaining a relatively looser grid in other regions.Although it may not accurately capture fine−scale turbulent features, URANS can yield reasonable and reliable results in practical engineering applications.In fact, URANS is widely used in aerospace engineering to simulate the flow around aircraft, projectiles and other aerodynamic bodies due to its ability to provide a reasonable compromise between computational cost and accuracy for predicting lift, drag and other aerodynamic coefficients [28,29].Therefore, URANS is adopted in this study.

Governing Equations and Turbulence Model
In order to precisely simulate the unsteady flow generated by the spinning projectile with deformation and to assess alterations in the lateral force, the integral form of the three−dimensional URANS equations is adopted as the governing equation, as follows: where Ω is the control volume.∂Ω and n represent the boundary and the outer normal direction of the control volume.t is the physical time.H is the source terms, which maintains zero in this study.W, F and G, respectively, denote the conserved variables, the convective terms, and the viscous terms, which are expressed as where the stress tensor is expressed as where denotes the fluid velocity vector.denotes the velocity of mesh induced by the spinning motion and bending deformation, which is determined by the dynamic mesh technology.
is the heat flux., p, and E are, respectively, the density, the pressure, the dynamic viscosity and the total energy per unit mass.
The shear−stress transport (SST) − turbulence model proposed by Menter [30] is employed to calculate the flow of the deformed spinning projectiles.Besides the free flow independence of the − model in the far field, it effectively utilizes the robustness and accuracy of the − model in the near−wall field [3].Moreover, the turbulent viscosity in the SST − model is modified to involve the transport of turbulent shear stress.Consequently, the SST − model is competent to simulate the flow around a deformed spinning projectile.

Dynamic Mesh Method
In this paper, the spinning motion and bending deformation are involved for the simulation.Therefore, the dynamic mesh method is exploited to adjust the mesh in the fluid domain.The rigid−motion mesh method is employed to treat the mesh motion caused by the rigid motion, while the RBF mesh deformation method is used to solve the morphed boundary due to the bending deformation.

Rigid−Motion Mesh Method
For the rigid−motion mesh method, the mesh is updated based on the translation and rotation of the rigid body.It can ensure the topological structure and the quality of the mesh with great computational efficiency.x y z  r of mesh node P determined by the rigid motion at any time can be described as where ∆ is the displacement caused by the translation, which is consistent with that of the centroid; ∆ represents the displacement caused by the rotation, and its expression is where is the centroid position.T represents the transformation matrix, which can be wri en as where , and , respectively, denote the rolling angle, pitching angle and yawing angle.

RBF Mesh Deformation Method
By employing the RBF mesh deformation method proposed by de Boer et al. [31], the boundary deformation will be diffused to the interior mesh nodes in the fluid domain.Generally, the matrix equations whose dimensions are equal to the number of mesh nodes in the boundary are needed to solve for the RBF mesh deformation method.Hence, the computational cost is very high for a large number of mesh nodes in the boundary.In order to improve the computational efficiency, the RBF mesh deformation method based on dynamic control points in Ref. [32] is adopted to reduce computational cost with a strong robust.

Discretization Methods
The dual−time stepping method [33] not only guarantees precise computational outcomes but also significantly enhances the computational efficiency.Accordingly, the present study adopts the dual−time stepping method to address the unsteady flow.An unsteady preconditioning [34] is usable when applying the dual−time stepping method.This approach aims to boost the accuracy of the solution by refining the scaling of artificial dissipation.Additionally, it seeks to optimize computational efficiency by optimizing the number of sub−iterations needed during each time step.Here, a preconditioned pseudo−time−derivative term is introduced into Equation ( 7), as follows: where denotes the pseudo−time employed in the time−marching procedure.Γ and Q, respectively, denote the preconditioning matrix and the primitive variables [35].
The spatial discretization employs the second−order upwind scheme.Within this scheme, the face value is determined through the following expression: where and ∆ , respectively, represent the value at the cell center and its gradient in the upstream cell.⃗ signifies the displacement vector from the upstream cell centroid to the face centroid.
The time discretization uses the second−order time integration scheme.The specific formula for this scheme is expressed as: where represents a scalar quantity.Δ is the time step.The function L encompasses the spatial discretization.

Grid Independence
The grid independence research is carried out by utilizing the two types of grid with different amounts.The specific grid parameters are shown in Table 1.The topologies of the two type grid are the same, but their amounts are, respectively, 6.47 million and 12.82 million.The incoming Mach number is Ma = 3.0, the angle of a ack is = 12 ° and the Reynolds number is Re = 1.54 × 10 .The spin parameter ( = 2 ⁄ ) is 0.11, and the time step is chosen as ∆ = 1 × 10 s.Moreover, the spin−deformation coupling motion trajectory is , shown in Figure 6. and the pitching moment coefficient at different grids are similar, with the relative difference being less than 0.2%.Moreover, the relative difference between the time−averaged lateral force coefficient and the time−averaged yaw moment coefficient under different grids is within 8%.Obviously, the calculation results using the two types of different grids are almost the same, which satisfies the grid independence.Therefore, the coarse mesh is determined to simulate the flow for the large−aspect−ratio projectile undergoing spin−deformation coupling motion due to its high computational efficiency.

Time Step Independence
The coarse mesh is used for the validation of time step independence.Three time steps are set to ∆ = 5 × 10 s, ∆ = 1 × 10 s and ∆ = 5 × 10 s.The other calculation conditions are consistent with those in Section 5.1.Figure 7 indicates the aerodynamic coefficients obtained within a spin cycle at different time steps.It can be observed that the difference between the results calculated at ∆ and ∆ is small, while the results obtained at ∆ are significantly different from those calculated at ∆ and ∆ .Therefore, the time step for the numerical calculation is chosen as ∆ .

Comparison Validation
The experiment data in Ref. [36] is exploited to validate the numerical method.The number of grids is approximately 6.5 million.The computational conditions are Ma = 3.0, = 0.11 and Re = 1.54 × 10 .Figure 8 shows the comparison between the CFD−based numerical results and the experimental ones.It can be seen that the time−averaged lateral force coefficient and yaw moment coefficient , obtained by numerical calculation, are in good agreement with the experimental data.Therefore, the numerical method can be employed to simulate the flow for spinning projectiles with large aspect ratios.In order to further verify the accuracy, the numerical method proposed in this paper is used to compute the aerodynamic characteristics of a spinning projectile with elastic deformation as outlined in Ref. [3].This computation adheres to the original author's specified conditions, encompassing Ma = 3.0, = 15 °, Re = 5.8 × 10 and = 0.283.Figure 9 gives the comparison between the numerical result obtained in this paper and that obtained in Ref. [3].This comparison reveals a strong agreement between the two, providing additional confirmation of the accuracy of the numerical method presented in this paper.

Effect of Deformation on Lateral Force
The reference length is the length of the Apache.The relative position of the centroid and the head vertex before the engine working is 0.54 [36].However, after the engine work is completed, the centroid is moved forward by 10% relative to the initial position.Thus, is determined as 0.44 in this paper.The other computational conditions are determined as shown in Table 3, unless otherwise specified.On this basis, the flow is simulated when the projectile undergoes spin ( ) and the spin−deformation coupling motion ( , ∈ [0, 90 °]).The specific trajectory shape of is shown in Figure 6.

The Variation with Trajectory
Figure 10 exhibits the variation of the time−averaged lateral force coefficient for the projectile with the trajectory .As shown in Figure 10, the value of decreases continuously with the trajectory changing from to at the different angles of attack .That is, the value of is the largest when the trajectory is , while it is the smallest for the .11 shows the variation of the time−averaged lateral force coefficient for the projectile body with the trajectory .It is found that as the trajectory changes from to , the variations of and are similar at = 4 °, while their trends are opposite at = 12 °.However, it is worth noting that the extreme values of and are all obtained when the trajectory is or .Therefore, the following will focus on the analysis of the lateral force when the trajectory is or .

Figure 11.
Variation of time−averaged lateral force coefficient for the projectile body with .

The Variation with Angles of A ack
Figure 12 exhibits the variation of the time−averaged lateral force coefficient for the projectile with angles of a ack at different trajectories.As shown in Figure 12, the variation of is similar at different trajectories, which increases first and then decreases with .When the trajectory is , the value of is maximum.Although the induced time−averaged additional sideslip angle is zero for the , the is still changed by the deformation, especially at large .

The Variation with Deformation Amplitude
Figure 15 shows the variation of the time−averaged lateral force coefficient for the projectile with the deformation amplitude A. Obviously, the value of the increases linearly with the increase of the A for the .When the trajectory is , the changes nonlinearly with the A, and the trends are different at different .Specifically, the value of decreases with the increase of the A when = 4 ° , but the value increases at = 12 °. Figure 16 shows the variation of the time−averaged lateral force coefficient for the projectile body with A. It is observed that for the , the value of the increases linearly with the A at = 4 °, while it changes from negative to positive and increases continuously with the A when = 12 °.When the trajectory is , the variations of the are consistent with that of the in Figure 15b, which are nonlinearly correlated with the A.

Mechanism on the Lateral Force Induced by Deformation
As mentioned in Section 6, the lateral force of spinning projectiles with large aspect ratios is significantly changed due to the bending deformation.In this section, the flow is analyzed for the projectile experiencing the spin−deformation coupling motion, which is exploited to reveal the mechanism for the change of the lateral force.The trajectories and characteristic moment shown in Figure 17 are employed to explore the flow.In Figure 17a, for the , the roll angles at the moments P and Q are, respectively, and 180 °− .At the moment S, the additional angle of a ack is zero.In Figure 17b, when the trajectory is , the roll angles at the moments P and R are, respectively, and − .At the moment T, the additional sideslip angle is zero.

Compression and Expansion of the Fluid
When the projectile deforms, the compression and expansion waves are produced in the axial flow, which affects the aerodynamic characteristics.Figure 18 shows the compression wave and expansion wave induced by the deformation.It is seen that for the projectile with deformation, the high pressure gas after the shock wave continues to expand and accelerate in the head.When it reaches the middle or tail of the body, the compression wave is generated on the concave side of the body, which leads to the decrease of velocity and the increase of pressure for the local airflow.Whereas, the expansion wave is produced on the convex side of the body, which causes the increase of velocity and the decrease of pressure for the local airflow.Figure 19 shows the effects of the compression wave and expansion wave on the aerodynamic characteristics of the fins.In Figure 19a, affected by the compression wave and expansion wave, a pair of forces and with the same direction are generated for the fin2 and fin4 at the characteristic moment P (the roll angle is γ).As shown in Figure 19b, when the trajectory is , there exists the moment Q (the roll angle is 180 °− ) where the lateral force only appeared after the superposition of and , and .As shown in Figure 19c, when the trajectory is , there exists the moment R (the roll angle is − ) where the normal force only appeared by superimposing and , and .Moreover, the interference of the compression wave and expansion wave on the aerodynamic characteristics of the fins causes the lateral and normal pressure center to move backward, which induces the increase of the pitching moment and yaw moment of the projectile.

Change of Vortex and Pressure Distribution
In this section, the effect of the spin−deformation coupling motion on the vortex and pressure distribution is researched to explain the reason for the deformation−induced variation in the lateral force further.For convenience, the circumferential angle is defined as shown in Figure 20.The point P is the origin of (i.e., = °), and increases gradually along the counterclockwise direction.The value of is defined as

Effect of the Trajectory
Figure 21 shows distribution of the time−averaged lateral force coefficient ( / ) for the and .It is observed that for the , the distribution curves at different are presented as parabolas when / ≤ 0.85.Affected by the fins, the variation of ( / ) for the and is similar when / > 0.85, which increases first and then decreases.Notably, the ( / ) for the head of the projectile body is consistent at different when the trajectory is , which indicates it is only related to the deformation.However, the ( / ) for the middle or tail is inconsistent, which declares it is related to the deformation and the flow structures.Firstly, the mechanism for the lateral force of the head is analyzed.Figure 22 exhibits the variation of the lateral force coefficient .at the cross section x/L = 0.1 when the trajectory is .Meanwhile, the additional sideslip angle ∆ is also shown in Figure 22.Notably, the physical quantities involved are normalized in the figure, namely x/L where ( ) represents the original physical quantity, and ( ) is its normalized result.denotes the maximum value of ( ).As shown in Figure 22, the variation of . is completely consistent with that of ∆ at different , which indicates the lateral force for the head mainly depends on the additional slip angle caused by the bending deformation.Secondly, the mechanism for the lateral force of the middle is explored.Figure 23 shows the variation of the lateral force coefficient .at the cross section x/L = 0.5 when the trajectory is .Meanwhile, the lateral component of is also shown in Figure 23.As shown in Figure 23, the variation of .at = 4 ° is consistent with that of , while the variation of the two is different at = 12 °.This manifests that the lateral force of the middle is mainly related to the deformation at small when the projectile deforms, but the influence of deformation is weakened with the increase of .In fact, as shown in Figure 18, the compression wave and expansion wave induced by the deformation inthe pressure on the concave side of the middle, while the pressure on the convex side is decreased.This asymmetry of the pressure distribution results in the negative lateral force of the middle at small induced by the deformation.However, it can be seen from Figure 21 that the value of the lateral force of the middle changes li le at large , which indicates that the increase of will weaken the asymmetry of the pressure distribution.Finally, the mechanism for the lateral force of the tail is discussed.Figure 24 shows the variation of the lateral force coefficient .at the cross section x/L = 0.8 when the trajectory is .Meanwhile, the additional sideslip angle ∆ and angle of a ack ∆ are also shown in Figure 24.It is seen that the change trend of the .at = 4 ° is similar to that of ∆ , which indicates that even if the time−average value of ∆ is zero, it will also affect the generation of the lateral force at small .The variation of the .at = 12 ° is similar to that of ∆ , but there exists the characteristic of phase lag between them.The specific flow field is analyzed to further illustrate the mechanism for the lateral force of the tail induced by the deformation in the following.
Figure 25 and Figure 26, respectively, show the pressure contour, streamlines and the pressure difference between the left and right sides of the cross section x/L = 0.8 when is 4 ° or 12 °.As shown in Figure 25, the change of pressure difference within ∈ (0 °, 40 °) is very small, which indicates that the change of the separated vortex induced by the deformation has li le contribution to the lateral force at small .When the trajectory changes from to , the pressure difference within ∈ (40 °, 82 °) is positive due to the offset of the low−pressure region on the left and right sides.However, the negative pressure difference within > 82 ° is induced by the offset of the high−pressure region on the windward side, which plays a leading role in the generation of the lateral force.Therefore, the change of the high−pressure region on the windward side induced by the additional sideslip angle will lead to the negative lateral force of the tail.
Combined with the analysis of Figures 24 and 25, it can be seen that the coupling of the additional angle of a ack and spinning motion determines the variation of the transient lateral force of the tail at small angles of a ack, while the additional sideslip angle will increase its time−averaged value.
Comparing Figure 25 with Figure 26, the variation of the pressure and vortex system at the cross section x/L = 0.8 under = 12 ° is similar to that under = 4 ° when the projectile deforms.However, the separation vortex is stronger when = 12 °.Thus, the decreased distance between the left separation vortex and the projectile surface leads to a significant decrease in the pressure for the upper−left of the cross section, while the increased distance between the right separation vortex and the projectile surface results in a noticeable increase for the upper−right.Those lead to the positive pressure difference within ∈ (0 °, 75 °).Moreover, besides the offset of the low−pressure region on the left and right sides, the obvious decrease of pressure for the low−pressure region on the left side also causes the pressure difference within ∈ (75 °, 110 °) to be continuously positive due to the interference of the separation vortex.Although the offset of the high−pressure region on the windward side leads to the negative pressure difference within > 110 °, it does not take up a leading role in the generation of the lateral force.Hence, the change of the separation vortex is the main cause of the positive lateral force of the tail induced by the deformation at large angles of a ack.Moreover, according to the analysis of Figure 24, the change of the separation vortex induced by the deformation can also lead to the phase lag between the transient lateral force and the additional sideslip angle for the tail at large angles of a ack.x/L of the tail is more sensitive to the change of the angle of a ack.Hence, the negative ∆ within / ∈ (0.53, 1] has a greater influence on the lateral force of the projectile body.Form the variation of the lateral force with in Figure 12, the lateral force will decrease when = 4 ° due to the negative ∆ , while it will increase at = 12 °.

Effect of the Trajectory
Figure 29 shows the lateral force coefficient .for the cross section x/L = 0.8 under different when the trajectory is .The additional sideslip angle ∆ and angle of a ack ∆ are also shown in Figure 29.In Figure 29a, the variation of . is extremely similar to that of ∆ when = 4 °, which indicates that the additional angle of a ack is the main reason for the change of the transient lateral force of the tail at small angles of a ack.As shown in Figure 29b, the variation of the . is similar to that of the ∆ at = 12 °, and there exists the characteristic of phase lag between them.This is consistent with the relationship between the lateral force of the tail and the additional a ack angle or the additional sideslip angle shown in Figure 23.

Conclusions
The spin−deformation coupling motion model is established, and its specific trajectory is analyzed by the numerical method.On this basis, the URANS equations are solved to simulate the flow over a large−aspect−ratio projectile experiencing spin−deformation coupling motion using the unsteady numerical simulation method and the dynamic mesh   Normalization Data t/s technology, and the lateral force is obtained.Furtherly, the flow mechanism for the changed lateral force induced by the bending deformation is analyzed.The conclusions are drawn as follows: 1.The variation of the transient lateral force for the head of the projectile is consistent with that of the additional sideslip angle induced by the bending deformation, and it is almost independent of the angle of a ack. 2. The compression wave and expansion wave will increase the time−averaged lateral force for the middle of projectile at small angles of a ack, while they have li le effect on its value at large angles of a ack. 3.At small angles of a ack, the coupling of the additional angle of a ack and spinning motion is the main reason for the change of the transient lateral force for the tail of the projectile, while the additional sideslip angle will increase its time−averaged lateral force.4. At large angles of a ack, the separation vortex will lead to a phase lag between the transient lateral force for the tail of the projectile and the additional sideslip angle, which will increase the time−averaged lateral force for the tail.

Figure 1 .
Figure 1.Model dimension of the Apache.

Figure 2 .
Figure 2. Definitions of spin direction, fins and roll angle.
local coordinate system Ox b y b z b is a moving coordinate system whose origin O coincides with the center of the cross section.The direction of the axis Ox b is consistent with that of the axis Ox.And, the axis Oy b is located in the longitudinal symmetry plane of projectiles, which is perpendicular to the axis Ox b .Then, the axis Oz b and the other two axes form the right−hand coordinate system.The position of point P( , , ) at the time t is, respectively, defined as P ( , , ) and P ( , , ) in the Oxyz and Ox b y b z b .The relationship between P ( , , ) and P( , , ) can be expressed as

Figure 4 .
Figure 4. Diagram of the spin−deformation coupling motion.

Figure 5 .
Figure 5. First order mode shape of the Apache.
the initial position of mesh node P, and the position  

Figure 7 .
Figure 7. Aerodynamic coefficients within a spin cycle at different time steps.(a) and ; (b) and .

Figure 10 .
Figure 10.Variation of the time−averaged lateral force coefficient for the projectile with the trajectory .

Figure 12 .
Figure 12.Variation of time−averaged lateral force coefficient for the projectile with .

Figure 13
Figure13shows the variation of the time−averaged lateral force coefficient for the projectile body with at different trajectories.It is found that the difference between and is small at the same when the trajectory is or .Moreover, although the is larger than the at same the when the trajectory is , the variation of the . Thus, the body determines the variation of the time−averaged lateral force coefficient of the projectile with at different trajectories.

Figure 13 .
Figure 13.Variation of time−averaged lateral force coefficient for the projectile body with .

Figure 14 depictsFigure 14 .
Figure 14  depicts the variation of the time−averaged lateral force coefficient for the fins with at different trajectories.It is seen that the variation of the for all fins with is consistent when the trajectory is , while the values of fin2 and fin4 are the largest.When the trajectory is , the of all tails change li le around the zero at different .

Figure 15 .
Figure 15.Variation of time−averaged lateral force coefficient for the projectile with deformation amplitude A. (a) ; (b) .

Figure 16 .
Figure 16.Variation of time−averaged lateral force coefficient for the projectile body with deformation amplitude A. (a) ; (b) .

Figure 17 .
Figure 17.Characteristic moment for the trajectory (a) and (b) .

Figure 18 .
Figure 18.Compression wave and expansion wave induced by the deformation.V and P, respectively, represent velocity and pressure.

Figure 19 .
Figure 19.Effect of compression wave and expansion wave on the aerodynamic force of the fins.(a) Moment P; (b) moment Q; (c) moment R.

Figure 20 .
Figure 20.Diagram of the circumferential angle in the cross section.

Figure 22 .
Figure 22.Lateral force coefficient .and additional slip angle ∆ at the cross section x/L = 0.1 for the .

Figure 23 .
Figure 23.Lateral force coefficient .and lateral component at the cross section x/L = 0.5 for the .

Figure 24 .Figure 25 .Figure 26 .
Figure 24.Lateral force coefficient., additional sideslip angle ∆ and angle of a ack ∆ at the cross section x/L = 0.8 for the .

Figure 27
Figure27shows the distribution of the time−averaged lateral force coefficient ( / ) for the and .It is seen that the variations of the at different are similar when / ≤ 0.85 for the and , which increase with / .Due to the larger value of the ( / ) of the tail at same , the additional angle of a ack greatly affects the lateral force of the tail.That causes the obvious change of the ( / ) of the tail for the .More specifically, its value decreases under = 4 °, but it increases when = 12 °.

Figure 28
Figure28shows the distribution of the additional angle of a ack ∆ induced by the .In Figure28, the distribution curve of the ∆ is approximately central symmetry about the point (0.53, 0).The ∆ decreases from ∆ to zero when / ∈ [0, 0.53], while it decreases from zero to −∆ for / ∈ (0.53, 1].As mentioned above, the lateral force

Figure 28 .Figure 29 .
Figure 28.Distribution of additional angle of a ack induced by the trajectory .

Table 2
displays the relative difference of time−averaged aerodynamic coefficients obtained by the two type grid.The values of the time−averaged normal force coefficient

Table 2 .
Relative difference of aerodynamic coefficients.