Abstract
In this study, the effect of the axial gap on the mechanical response of a cartridge-loaded propellant grain under vibration loads is investigated. The wide strain rate range of uniaxial compression tests () on the composite modified double base (CMDB) propellant was carried out by using a universal testing machine, a hydraulic testing machine, and a split Hopkinson pressure bar system, respectively. A linear viscoelastic constitutive model of the CMDB propellant was developed by using the experimental measurements. The results show the studied CMDB propellant has a strong strain rate dependence, exhibiting an initial linear elasticity followed by a strain hardening region. The dynamic process of collision between the propellant grain and the motor case in the axial direction induced by vibration loads was simulated with the developed constitutive model by using the finite element method. The effects of the gap size between the propellant grain and the case and the vibration frequency on the mechanical response of the grain were studied. This shows that with a constant vibration frequency, the stress of the grain increases first and then decreases with increasing gap size. Moreover, the stress increases with increasing vibration loads.
1. Introduction
In solid rocket motors (SRMs), one of the methods of holding the grain in the case is using cartridge-loaded propellant grains [1–6]. In cartridge-loaded solid rocket motors, the propellant grains are manufactured separately from the case and then loaded into or assembled into the case. The propellant grain is normally fixed by the forward and backwards support inside motor cases [7–11]. The cartridge-loaded SRM has the advantage of low cost and ease of inspection and is commonly used in small missiles and a few medium-sized motors [12–14]. However, due to temperature variation, aging, and other issues, the size of the propellant grain could decrease gradually, leading to an axial gap between the propellant grain and the support [15–17]. This gap causes the propellant grain and the support gasket to collide with each other under vibration conditions, possibly resulting in the destruction of the structural integrity of the propellant grain. The vibration load during the transport subjected to the propellant has a certain periodicity, and the stress and strain of grain may exceed the capacity range of the solid propellant causing crack propagation in the propellant [18, 19]. These cracks can lead to the ignition failures of grain and even the explosion of the solid rocket motor. Therefore, it is important to evaluate the effect of the axial gap on the mechanical response of cartridge-loaded solid propellant grains under vibration loads.
Composite-modified double-base propellants (CMDBs) have been widely used in cartridge-loaded SRMs, owing to their excellent properties of low characteristic signals, good combustion performance, and exceptional processability [7, 20, 21]. CMDB propellant is a polymer composite material filled with a certain number of solid granules. The mechanical properties of the CMDB propellant have a great dependency on the strain rate [16, 22–25]. For a comprehensive understanding of the effect of the axial gap on the mechanical response of the CMDB propellant grain under vibration conditions, the deformation of the solid propellant material and the vibration conditions to which the propellant grain is subjected should be modelled accurately. A resume of literature reveals that there have been various studies to investigate the reliability of solid propellants under vibration loads [26–30]. For instance, Yilmaz et al. presented a methodology considering damage accumulation for the assessment of the service life of SRMs under random storage and transportation loads [18]. Kohsetsu proposed a simplified method to generate the structural vibration model of an SRM based on the mixed finite element method (M-FEM) for the analysis of carrier rocket vibration [31]. Kunz built a type of linear cumulative damage (LCD) model of a solid propellant under a variety of loading histories [32]. Zhang et al. used pavement inequality excitation to calculate the vibration response of an SRM by numerical simulation [33]. Cao et al. studied the structural integrity in the process of SRM launching and transporting by the vibration data monitored on the practical transport process of SRMs [34]. Huang and Zhang analyzed the effects of temperature and low-frequency vibration loads on the stress and strain of SRMs based on monitored temperature and vibration data [35].
The aforementioned studies of vibration load subjected to solid propellants are mainly about the reliability, vibration test, and new damage constitutive model development of solid propellants [17, 32, 36, 37]. However, mechanical response analysis for cartridge-loaded CMDB propellant grains with axial gaps under vibration loads has rarely been studied. To better understand the collision process caused by the gap between the case and the propellant grain under vibration loads, not only is the viscoelastic constitutive model of the propellant capable of accurately predicting the mechanical behavior of the propellant over a wide range of strain rates needed, but the influence of the vibration frequency and gap size on the collision process also needs to be studied.
Therefore, the objective of the present work is to investigate the effect of the axial gap on the mechanical response of a cartridge-loaded CMDB propellant grain under different vibration conditions. First, a linear viscoelastic constitutive model of the CMDB propellant was developed based on axial compression measurements conducted over a large range of strain rates. Then, the dynamic collision process of a cartridge-loaded propellant grain under vibration loading was simulated by using the developed constitutive model with the finite element method. The dynamic stress–strain distribution of the grain was obtained. The effects of the axial gap between the case and the propellant grain and the vibration frequency on the stress response of the propellant grain were discussed in detail.
2. Methodology
2.1. Material Properties
The components of the CMDB propellant studied in the present work are shown in Table 1. The nitrocellulose (NC) and nitroglycerin (NG) of the propellant are used as a double base binder, cyclotetramethylene tetranitramine (HMX) is used as an oxidizer, and the others are used as additives, including lead/copper-based compound salt as a combustion catalyst and dimethyl diphenylurea as a stabilizer. The mean size of the HMX particles used is approximately 90 μm. The degree of esterification of the NC is approximately 11.9%. Figure 1 shows the microstructure of the CMDB propellant, indicating that the HMX filler is homogeneously distributed in the binder matrix. Before the test, all specimens were stored in a vacuum oven at 20°C for 24 h to eliminate the residual stress and humidity of the propellant.

(a) ×1000

(b) ×5000
2.2. Cartridge-Loaded Solid Motor
In this paper, the solid rocket motor studied is a cartridge-loaded type with an end-burning propellant grain, as shown in Figure 2. The solid rocket motor consists of a motor case, an end propellant grain, the insulation, and the motor case.

The motor case of the solid motor is made of steel material with a thickness of 4 mm, which is considered an elastic material. To study the effect of gap size on the collision process, the axial gap between the case and the propellant grain is prospectively set to 0 mm, 0.3 mm, mm, 1 mm, 2 mm, 3 mm, π mm, 4 mm, 5 mm, mm and 2π mm. Table 2 lists the mechanical properties of the case and grain used in the simulation. The geometrical details of the motor case and the propellant grain are shown in Table 3.
2.3. Numerical Setup
2.3.1. Constitutive Model
CMDB propellants are a type of particle-filled polymer composite material. The mechanical properties of the CMDB propellant are affected by the strain rate. As the collision process between the propellant grain and the motor cases is a typical dynamic impact issue, an accurate constitutive model of propellant grain under dynamic conditions must be first established. Assume that the propellant is homogeneous and linear viscoelastic. The Poisson ratio of propellant is constant. In the present work, the Prony series constitutive model applied a wide strain rate was established by using the data of uniaxial compression tests at a wide strain rate ().
The results of uniaxial compression tests show that the CMDB propellant exhibits linear mechanical behavior before the yield point. Therefore, a single integral linear viscoelastic constitutive model which is based on the Boltzmann superposition principle and generalized Maxwell model is used to describe the stress–strain behavior before the yield point of the CMDB propellant. The formula of the constitutive model used in the study is
is the true stress, and is the elastic modulus. In the developed constitutive model, the elastic modulus is usually expressed in the form of a Prony series.
is the equilibrium modulus; and are the modulus and relaxation time of each Prony series component.
The constitutive model of the CMDB propellant before the yield point is
It is considered that the propellant specimen deforms at a constant strain rate in the uniaxial compression test.
2.3.2. Basic Assumptions and Meshing
First, the following approximations have been made. In the present work, the energy loss caused by the collision and the effect of the temperature change caused by the collision are not considered in the simulation. The contact surfaces between the case and the propellant grain are set as frictionless. To improve the computational efficiency, in our simulation, the geometry of the studied solid rocket motor was simplified to an asymmetric two-dimensional model consisting only of the motor case and the propellant grain, as illustrated in Figure 3.

The uniform structure quadrilateral mesh is used into the grain and case. Precalculations with mesh sizes of 0.5 mm, 1 mm, and 1.5 mm are initially carried out to make sure the mesh size is small enough for an accurate calculation. Figure 4 shows the displacement of the propellant grain’s top center at various mesh sizes. The axial displacement curves obtained using three different mesh sizes are very similar, illustrating that the mesh size between 0.5 and 1.5 mm is tiny enough to obtain accurate results. The propellant grain mesh size is set to 1 mm for all the conditions described in the next sections. The mesh of the case is divided more finely with a size of 0.5 mm. The number of propellant grain’s total elements is 4480, and the number of motor case’s total element is 5375 potentially. It should be noted that the insulation is not considered in our simulations. If the insulation is considered, the calculated von Mises stress will decrease due to the cushion action of insulation. However, the existence of insulation will not change the tendency of the effects of the gap size and vibration frequency. Therefore, the insulation was not considered in the simulation model.

2.3.3. Boundary and Load Conditions
In our simulation, the vibration load is achieved by adding a periodic movement of the motor casing impacting the fuel grain. The periodic displacement applied to the case follows the equation below:
and are the amplitude and frequency of the vibration, respectively. The displacement curve of the case is shown in Figure 5. In the whole vibration process, the motor case drives the propellant grain to move. By simulating the motion of the grain and the motor case under vibration load, the stress of the contact surfaces between the grain and the case in the collision process is calculated.

2.4. Uniaxial Compression Experiments
To establish the constitutive model of the CMDB propellant suitable for simulating the mechanical response of the propellant grain under vibration loads. A split Hopkinson pressure bar, an Instron VHS 160/100-20 high-speed hydraulic servo testing machine, and an Instron 4505 universal testing machine were used to carry out uniaxial compression tests of the CMDB propellant at low, medium, and high strain rates. A wide range of strain rate uniaxial compression tests of the CMDB propellant were performed. The measurement matrix of the strain rate is outlined in Table 4. Each strain rate test was performed at least five times to confirm its correctness and reliability, and the average of those results served as the final experimental measurement values.
The specimens were produced as 16 mm in diameter and 20 mm in length in the low strain rate compression tests. In order to remove any variability in the results brought by the selection of dimensions, the specimens in the intermediate strain rate tests with dimensions identical to those used for low strain rate compression testing were tested. Figure 6(a) depicts the intermediate strain rate compression testing device and specimens. The specimens had a diameter of 10 mm and a length of 5 mm, and they were cylindrical for the high strain rate compression test. The length-to-diameter ratio was designed as 0.5. Figure 6(b) shows a schematic representation of the intermediate strain rate compression testing device and specimens. By such a design, it can minimize wave attenuation in the strain signals recorded [16, 38] and the effects of radial and longitudinal inertia in the specimen [39]. The CMDB propellant which was used in all tests was in its original state. In our earlier study by Yang et al., the hydraulic testing device that was used here is described in great detail. [25].

(a) Intermediate strain rate compression

(b) High strain rate compression
3. Results and Discussion
3.1. Compression Test
The obtained stress–strain curves of the CMDB propellant under different strain rates are shown in Figure 7. The strain rate range is from to .

In Figure 7, similar characteristics can be found in the stress-stain curve of the CMDB propellant at a wide range of strain rates: initial linear elasticity, then strain hardening, and the effect of strain hardening increases with the strain rate. These characteristics show the substantial ductile behavior of the CMDB propellant at a wide range of strain rates. Furthermore, the mechanical properties of the CMDB propellant are strain-rate dependent. Moreover, the yield stress of the CMDB propellant increases with increasing strain rate. These characteristics should be considered in the developed constitutive model.
The result of uniaxial compression tests of the CMDB propellant shown in Figure 7 was used to fit by formula (3). The nonlinear least square regression model was used during the fitting process to get the minimal value of the objective function. The objective function is
is the number of strain rates; is the number of data points in the single strain rate; is an unknown parameter in the Prony series; is the fitting stress; and is the test stress corresponding to the present strain rate and strain. Table 5 shows the Prony series fitted from CMDB propellant uniaxial compression test data. Figure 8 shows the comparison of the strain–stress curves predicted by the developed constitutive model and those obtained by tests. The results calculated by the developed constitutive model and the test results show good agreement.

3.2. Mechanical Response of Propellant Grain
By using the developed constitutive model, the dynamic collision process between the case and the propellant grain under a vibration load of 20 Hz and an amplitude of 1 mm is simulated. In this section, the radial gap size is set to 5 mm, and the axial gap size is set to 1 mm.
As illustrated in Figure 9, when the motor moves with a cosine function, the motor case and the propellant grain will collide periodically. It should be noted that the stress on the bottom or the top surface of the propellant grain changes over time and is not strictly sinusoidal periodic. Figure 10 shows the maximum stress appears at the first third or fourth collision between the grain and the case; then, the stress decreases gradually and tends to a stable periodic value after 3 collisions.


The maximum radial displacement distribution of the grain is illustrated in Figure 11. The maximum radial displacement of the top surface is . The maximum radial displacement of the bottom surface is . As the radial gap of the case and grain is 1 mm, the radial deformation of the grain is far less than the width of the radial gap, so the grain was not in contact with the case during the collision. This is crucial to avoid friction ignition induced by the collision between the grain and the motor case.

(a)

(b)
During the first collision of the grain and case, the von Mises stress distribution variation versus time is shown in Figure 12. The first collision between the grain and case starts at 45.6 ms and ends at 46.2 ms. The total interaction time of the collision time is approximately 0.6 ms, indicating that the instantaneous collision time is very short. During the collision process, the stress wave propagates in the grain quickly. The maximum stress points are located at the center of the top and bottom surfaces of the grain. According to our calculation, the maximum stress at the center of the top surface is 3.0 MPa, and the maximum stress at the center of the bottom surface is 3.6 MPa. Therefore, the center of the bottom and top surfaces of the grain are the danger points.

During the first collision of the grain and case, the strain distribution of the propellant grain is shown in Figure 13. It can be found that the maximum strain points are located at the center of the contact surfaces between the grain and case. The maximum strain at the center of the top surface is , and the maximum strain at the bottom surface is . It should be noted that in our simulation, the effect of friction is not considered. If the friction effect between the grain and the case is considered, the calculated von Mises stress will possibly decrease due to the limiting action of friction on the grain along the radial direction.

3.3. Effect of the Gap Size
To study the effect of the gap size on the mechanical response of the propellant grain, simulations are performed by changing the gap size to 0 mm, 0.3 mm, mm, 1 mm, 2 mm, 3 mm, π mm, 4 mm, 5 mm, mm and 2π mm with the same vibration load of frequency of 20 Hz and amplitude of 1 mm.
Figure 14 illustrates the stress variation of the top and bottom of the propellant grain with different radial gap sizes. This shows that the stress value has a significant dependency on the gap size. When there is no gap between grain and case, i.e., there is no relative motion between them. The stress value varies synchronously with the vibration load, and the stress variation is mainly due to the stress wave propagation inside the grain. The result is shown in Figure 15(a). The maximum stress value is small, with a value of approximately 0.006 MPa. When the gap size increases, the stress value increases dramatically due to the collision between the propellant and the case. As the collision strength is mainly dependent on the impact velocity, when the gap size increases, the relative velocity increases, yielding higher stress values, which are shown in Figures 15(b) and 15(c). Meanwhile, as the motion velocity of the motor case follows a sinuous function, when the gap size reaches 2π, the case and the grain move in the same direction, leading to a small relative velocity at the time of the collision. Therefore, the maximum stress value becomes tiny, as shown in Figure 15(d).

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)
The maximum stresses of the top and bottom of the grain at the different gas sizes are shown in Figure 16. The maximum stress of the grain increases gradually with increasing gap size. When the gap is 6 mm, the maximum stress value is 3.55 MPa. When the gap is 2π mm, at the first collision, the relative velocity is 0 between the grain and the case. Therefore, the maximum stress value decreases sharply with a maximum value of 0.30 MPa, which is close to the case in which there is no gap between the propellant and the motor case.

3.4. Effect of Vibration Frequency
To study the effect of the vibration frequency on the mechanical response of the propellant grain, simulations are performed by changing the vibration frequency from 1 Hz to 50 Hz with the same gap size of 0.3 mm and amplitude of 1 mm.
Figure 14 shows the stress variation of the center points of the top and bottom surfaces at different vibration frequencies. With increasing frequency, the number of collisions of the grain and the case increases, and the collision period is shortened. Figure 17 shows the frequency-stress curves. This figure shows that with increasing frequency, the maximum stress of the grain also increases.

4. Conclusion
Based on the low, intermediate, and high strain rate uniaxial compression measurements of the CMDB propellant, a developed constitutive model which is linear viscoelastic was established. The mechanical response of a cartridge-loaded CMDB propellant grain with an axial gap under a vibration load was simulated and analyzed. The effects of the gap size and the vibration frequency on the response of the propellant grain under vibration loading were discussed. The main conclusions are as follows: (i)The mechanical properties of the CMDB propellant have a strong strain rate dependency. It exhibits an initial linear elasticity, then strain hardening, and the effect of strain hardening increases with a higher strain rate(ii)Due to the existence of the axial gap between the case and the grain, the grain collided with the case during the periodic movement of the case. Due to the existence of the gap, the stress value of the propellant increased significantly compared with the situation without radial gaps(iii)By analyzing the stress and strain distribution of grains in the process of vibration, it is found that the top and bottom of the grain are dangerous positions. Moreover, as the gap size increases, the stress of the propellant grain first increases and then decreases. There exists a gap size that makes the grain obtain the maximum stress value. As the vibration frequency increases, the stress of the grain also increases
Data Availability
Data are available within the manuscript.
Conflicts of Interest
The authors declare that they have no conflicts of interest.