Research on the Shock Response Characteristics of a Threaded Connection Using the Thin ‐ Layer Element Method

: Nonlinear factors such as the contact stiffness and friction damping at the threaded interface of a projectile–fuse system significantly affect the dynamic response characteristics. To obtain the dynamic response of the fuse body accurately during penetration, it is necessary to characterize these nonlinear factors reasonably. Because the existing structural dynamics software cannot effec-tively deal with nonlinear factors, the thin-layer element method was used to represent the nonlinear factors in this study. By combining the thread elastic model with thin-layer element principles, an effective method for determining the material parameters of the thin-layer element was established theoretically, which provided a different method of determining material parameters, not just relying on experiments. The accuracy of the material parameters was verified based on modal experiments with threaded tubes having different specifications. The errors were within 5%, indicating the reliability of the theoretical determination method for the material parameters. In addi-tion, projectile penetration into a semi-infinite concrete target was tested to verify the accuracy of the thin-layer element modeling. Compared with the ‘TIED’ constraint method, the resonant frequency obtained with the thin-layer element method was in better agreement with that of the experimental data. The maximum error decreased from 15.7 to 7.8%, indicating that the thin-layer element method could accurately represent the nonlinear factors. Thus, this study serves as a reference for accurately evaluating the dynamic response of the fuse body of a penetrator.


Introduction
As threaded connections are easy to install, disassemble, and maintain, the connection between a projectile and fuse is typically a threaded connection. During penetration of a projectile into a hard target, owing to the discontinuity and friction damping of the threaded connection interface [1], there is significant uncertainty in the transmission of overload from the projectile to the fuse, which can cause damage to the key electronic components inside the fuse. Finite element (FE) simulation is an effective method for studying the response of projectile-fuse systems under penetration overload. Because commonly used structural dynamics analysis software cannot deal with nonlinear factors, the threaded contact is usually simplified as a rigid connection for modeling. However, this causes the calculation result to be larger than the actual value. In the field of signal recognition for high-speed penetration of multi-layer hard targets, it is crucial to evaluate the vibration response characteristics of the projectile-fuse system accurately. To improve the accuracy of the simulation, the contact stiffness and friction damping of the threaded connection should be accurately characterized. Therefore, establishing an accurate modeling method for the threaded contact is of great significance for studying the vibration response of the projectile-fuse system. Much research has been conducted on the response mechanism of projectile-fuse systems under external loading. Zheng et al. [2] used a split Hopkinson pressure bar (SHPB) to investigate the loosening response of a threaded connection under a series of impact loadings. The correlation between the preload variation and the loosening rotation angle was established using a high-speed camera and digital image correlation (DIC). Zhang et al. [3] proposed a collision model for threaded connections and obtained calculation methods for the collision time and collision force. An FE model was established to study the collision process of the threaded interface between the projectile body and the fuse body. In addition, the wavelet decomposition of the measured deceleration signal indicated the existence of a thread collision frequency between the projectile and the fuse. Guo et al. [4] studied the vibration and energy dissipation efficiencies of flange connections under a series of impact loadings. High-frequency excitation amplified the nonlinear characteristics of the interface, resulting in shock-wave distortion and changes in the shock-wave frequency. As the excitation frequency decreased, the influence of the nonlinear characteristics of the interface also decreased. The above studies investigated the dynamic response of the projectile-fuse system experimentally, and demonstrated that the nonlinear characteristics of the connection interface had a significant impact on the dynamic response of the system. Some works in the literature [5][6][7] actively introduced nonlinear factors, such as contact stiffness and friction damping, which were in order to adjust the dynamic response characteristics of the combined system. During their works, they adopted a similar method which added coatings on the contact interface and changed the dynamic response characteristics of the combined system by adjusting the material properties of coatings.
In addition, much research has been conducted on the FE simulation of threaded connection interfaces. In particular, to characterize the nonlinear characteristics of the contact interface, the current treatment is to establish a fine solid-thread model in the explicit calculation [8,9]. The establishment of fine solid threads not only increases the complexity of the modeling, it also significantly increases the calculation time. Therefore, threaded connections are often ignored when the structural response is not of concern [10]. However, in the study of the structural response problem, the available structural dynamics software cannot deal with nonlinear factors, and ignoring the threaded connection will have a significant impact on the results. In the literature, the contact interface parameterization method has commonly been used to deal with this type of problem, e.g., through spring elements, general elements [11], the virtual material method [12], zero-thickness elements [13], and thin-layer elements (TLEs) [14]. Because of its easy of modeling and the convenience of its parameter determination, the TLE method has been widely used.
Desai [15] first proposed the concept of TLE modeling in 1984. The interface between the rock and soil was characterized by a TLE with normal contact stiffness and shear stiffness. Zhao et al. [16] used a TLE to represent the contact stiffness of a bolt joint. The variation in the contact stiffness at the bolt interface with the bolt preload was obtained. Lothar et al. [14] applied TLEs to represent the shrunken joints of a two-disk rotor in a generator. The damping and contact stiffness parameters obtained experimentally were coupled to the TLE model. An FE model was then used to predict the response characteristics of the generator. Alamdari et al. [17] introduced TLEs into the FE analysis of threaded pipes coupled to each other by a nut interface. The material parameters of the TLE were adjusted to minimize the residual between the numerical and experimental nonlinear frequency responses. This method significantly reduced the uncertainty of the connection modeling. Schmidt et al. [18] determined the material parameters of a TLE through general experiments. The parameters were coupled to the FE model to study the vibration and damping characteristics of steel composite structures. The numerical modal analysis results were in good agreement with the experimental results. In the above studies, the TLE method was mainly applied to bolt connections or other joint styles, and established at the contact interface to represent the contact stiffness and friction damping. To date, the TLE method has rarely been used to simulate threaded connections. Considering that the nonlinear factors of threaded connections also comprise the contact stiffness and frictional damping, the TLE method can be extended to threaded connections for further exploration.
In the application of the TLE method described above, the TLE material parameters were determined by general experiments [14,18] or identified using experimentally measured data [12,16,17,19]. Both of these two methods relied on experiments, which complicated the application of the TLE method.
In this study, the TLE method was used to characterize the threaded connection of a projectile-fuse system. The material parameters of the TLE were determined theoretically by combining the elastic model of the threaded connection with the basic principles of the TLE. The projectile penetration into a semi-infinite concrete target was tested. Frequency response analysis of the projectile-fuse system was performed using the TLE method. The accuracy and reliability of the TLE method were verified based on a concrete penetration test.

Experimental Setup
Several modal experiments of threaded tubes were conducted to investigate the influence of the nonlinear factors of the thread interface on the vibration response characteristics. The specimens were composed of two hollow cylindrical pipes of 45 steel connected by threads. Figure 1a shows the geometries of the threaded tubes. Five sets of threaded tubes (40_3.0, 40_2.0, 40_1.5, 30_2.0, and 50_2.0) were designed, as listed in Table 1. All of the threaded tubes (6H/6g) had the same machining precision and assembly accuracy; only the thread length and pitch varied. The label 40_3.0 refers to the thread length (40 mm) and pitch (3.0 mm). In addition, threadless specimens (L30, L40, and L50) were designed and tested in the same group, as shown in Figure 1b. The label L30 represents a threadless specimen with an overlap length of 30 mm.    Figure 2 shows the details of the test setup. The threaded tubes were suspended using elastic ropes to simulate free boundary conditions [18]. An electromagnetic shaker (HEV-20) was used to excite the specimens. The excitation point was selected at location Z1, and two three-directional accelerometers (ENDEVCO 2250AM1-10) were attached at positions A1 and A2 of the threaded tubes. The vibration test system (LMS Inc., Leuven, Belgium) consisted of a data acquisition instrument and vibration analysis software (LMS Test.Lab). The 'MIMO Sweep & Stepped Sine Testing' module of the LMS Test.Lab software was applied for auxiliary measurement and analysis, and a 16-channel data acquisition instrument was used to receive and process the signal. The sampling rate was set to 20.48 kHz. In Figure 2, the swept signal generated by the vibration test system is transferred to the electromagnetic shaker after power amplification. The shaker then receives the swept signal to generate an excitation force. The excitation force is collected by a force transducer (PCB 208C01) installed at the end of the shaker and sent back to the vibration test system. Thus, a feedback loop is formed to maintain a stable excitation force.

Modal Test of the Threadless Structure
The swept-sine signal can ensure the repeatability of the experiment and is easily used to measure the natural frequency and frequency response function [20]. Therefore, the swept-sine signal was selected as the excitation signal. The swept range of 100 to 5000 Hz could cover the first three order frequencies of the structure. The swept speed and excitation level were 5 Hz/s and 0.1 N, respectively. Figure 3 illustrates the frequency response function (FRF) measured at location A2 of the threadless specimen (L40). The first-, second-, and third-order modal frequencies are 938.26, 2669.86, and 4582.51 Hz, respectively. The obtained FRF has sharp peaks, and the peaks are symmetrical with respect to the resonance, indicating that frictional damping is negligible and the FRF has prominent linear characteristics. For comparison with the experimental results, we obtained the first three order frequencies of the threadless specimen (L40) using the MD/Nastran software, as indicated in the second row of Table 2. The errors between the test and simulation are 1.1, 1.3, and −0.2%, respectively. The difference can thus be considered negligible, which confirms that the threadless specimen has no nonlinear factors. Because of the negligible difference between the test and FEM in the threadless structure, we used the simulation results instead of experimental results to investigate the threadless structures (L30 and L50), as indicated in Table 2   There are three states associated with a mechanical joint subjected to increased loading: sticking, micro-slip, and macro-slip. The transition from one state to another alters the interface stiffness and frictional dissipation, which corresponds to the hysteresis mechanism [21]. The thread contact surface appears to stick under a low excitation force. However, with the gradual increase in the excitation, the contact status transforms from sticking to micro-slip or even macro-slip. The effect of nonlinear factors on the vibration behavior will gradually increase. The contact characteristics of threaded connections are generally affected by several factors, such as the roughness, lubrication, and thread preload. In this experiment, the thread preload was determined as a factor in the investigation of the two connection states (tight and loose). Because the thread preload cannot be accurately quantified, the loose thread state was obtained by turning the tight thread back one turn. The excitation levels were 0.8 and 0.1 N. To improve the measurement precision, the segmented swept method was applied to the threaded specimen (40_2.0). The excitation frequencies were divided into three sections: 100-1500, 2000-3500, and 3500-5000 Hz. The swept speed was 5 Hz/s. Table 3 presents the measurement results. In Table 3, each natural frequency of the threaded tube decreases by varying degrees compared with the threadless specimen (L40), and the natural frequency of the loosely connected thread decreases more significantly. In addition, when the excitation level increases from 0.1 to 0.8 N, the natural frequencies of the tightly connected threads change slightly, whereas the natural frequencies of the loosely connected threads change significantly. This indicates that the increase in the excitation level does not change the contact state of the tight thread; however, the contact state of the loose thread changes from sticking to micro-slip or macro-slip. Figure 4 shows the FRF curves of the threadless specimen (L40) and threaded specimen (40_2.0) obtained in the experiments and FEM. It can be seen that the FRF of the threadless specimen obtained with the FEM appears similar to that obtained in the experiment. The first natural frequency is 948.5 Hz, with an error of 1.1% compared with the experiment. In addition, compared with the threadless structure, the measured natural frequency of the threaded specimen decreases from 938.26 to 815.84 and 753.92 Hz. The first natural frequency of the loosely connected thread decreases more obviously, indicating that the loose connection state is more significantly affected by nonlinear factors than the tight connection state.
The FRF of the loosely connected thread is compared with that of the tightly connected thread. It can be seen that variation in the excitation level will result in the transition of the thread contact state. With increasing excitation force, the first natural frequency of the loosely connected thread decreases from 753.92 to 748.84 Hz, and the amplitude of the FRF peak decreases from 1.32 to 1.18. This results from the transition from sticking to the micro-slip or macro-slip state of the thread contact with the increase in the excitation level. Accordingly, the frictional damping increases, causing the amplitude and natural frequency to decrease simultaneously. For the tight connection, the first natural frequency decreases slightly from 815.84 to 814.18 Hz, and the amplitude of the FRF peak increases from 1.55 to 2.02 as the excitation force increases. Thus, the tight connection exhibits an entirely different form from that of the loose connection. It is concluded that the tightly connected threaded contact is initially in the sticking state, and the increase in the excitation level does not alter the thread contact state. Therefore, to ensure the simulation precision, it is not advisable to treat the threaded contact as a 'TIED' constraint. The thread contact state should instead be accurately characterized according to the actual situation of the thread contact.

Modal Tests of Threaded Structures with Different Specifications
Because the loosely connected thread is a non-working condition, the tight connection state was adopted for subsequent calculations of the thread connection axial stiffness.
To study the influence of the thread pitch and length on the vibration response characteristics, several modal experiments were performed for tightly connected thread structures with different thread pitches and lengths. The excitation level had little effect on the results of the tightly connected thread. Therefore, the excitation level was selected as 0.1 N. Moreover, other experimental conditions were the same as those described in Section 2.3.1. Table 4 lists the first three order frequencies of the threaded specimens in the modal tests. It can be seen that the change in pitch has a negligible effect on the frequency of each order. Each natural frequency varies slightly as the pitch decreases. However, each natural frequency increases significantly as the thread length increases. This results from the increase in the structural stiffness of the threaded connection with the increase in the thread length.

TLE Parameter Determination
Based on the results of modal experiments, nonlinear factors such as the contact stiffness and friction damping significantly affect the dynamic response characteristics of the threaded structure. To ensure simulation precision, these nonlinear factors should be accurately characterized during the modeling process. The TLE was adopted to represent the contact stiffness and frictional damping of the threaded connection. Based on the thread elastic model, the material parameters of the TLE were determined theoretically. Figure 5 illustrates the threaded connection, including the nut body and bolt body. The thread is a critical component in connecting the bolt body and nut body. The interaction between the bolt body and nut body is transferred by the thread body. Among the analytical methods for the thread load distribution, the Sopwith method [22] and Yamatoto method [23] have received extensive recognition. These two methods are similar. They both regard the thread as a cantilever beam, and reflect the thread load distribution by the thread axial load distribution. On this basis, Zhang et al. [24] considered the load distribution of threaded connections and proposed an accurate calculation method for the axial stiffness of a threaded connection. The axial stiffness of the threaded connection can be expressed as follows:

Axial Stiffness of the Threaded Connection
where b E and n E are the Young's modulus of the bolt body and nut body, respectively; b S and n S are the cross-sectional area of the bolt body and nut body, respectively; and l represents the thread length.  can be expressed as follows: where ba k and na k represent the axial stiffness of the unit length for the threads on the bolt and nut, respectively. These can be expressed as follows: where ba  and na  represent the elastic deflection caused by the unit force per unit width along the axial direction, and  is the lead angle. Zhang et al. [24] presented a calculation method for ba  and na  without considering the contact friction of the threaded teeth. Lu et al. [25] proposed a calculation method for ba  and na  considering the friction force of the threaded teeth. As the consideration of friction complied with the actual situation, the method considering the friction force [25] was adopted for the present calculation.  Figure 6 shows the axial deflection of the TLE under axial force. For the threaded connection under axial force F, as shown in Figure 6, the shear stress  generated by the axial force is:

Determination of the TLE Material Parameters
where  is the shear angle, A is the contact area, u is the axial displacement, and t is the thickness of the TLE. From Equation (4), the axial force F can be expressed as: Substituting Equation (1) into Equation (5), the shear modulus s G of the TLE can be expressed as: where  is calculated using Equation (2).
where the subscripts n, x, and y refer to the contact normal and local x-and y-tangential directions, respectively. The TLE can be simulated as an isotropic material [19,27]. In this situation, parameters 11 c and 55 66 c c  are the Young's modulus ( t E ) and shear modulus ( t G ), respectively, which define the normal and tangential stiffness of the joint.

Reliability Verification of the TLE Modeling
To investigate the simulation accuracy of the TLE method, we performed a modal analysis of threaded structures with different pitch and thread lengths using the TLE method. From Equation (6), the tangential moduli of the TLE with different pitch and thread lengths were obtained. The Young's modulus could be obtained according to 16,28], where G is the tangential modulus and  ( = 0.269) is Poisson's ratio, as indicated in Table 5. Shokrollahi and Adel [29] proposed a new concept for the contact influence area. It is assumed that the joint effects are extended to a region on both sides of the contact interface, which is referred to as the 'joint affected region'. This study extends the thread contact influence area t/2 in the normal direction according to the actual thread tooth height, and a TLE with a thickness of t was used to simulate the thread influence area. The thread of the threaded tubes was an ISO metric triangular thread. For the convenience of modeling, the TLE thickness was determined to be 2 mm according to the height of the thread teeth, and the TLE length was the thread length. An FE model of the threaded structure was created in MD/NASTRAN, consisting of three parts: an internal thread tube, external thread tube, and TLE, as shown in Figure 7. The internal and external thread tubes and TLE were created using an eight-node solid element (CHEXA), and the contacts between the TLE and threaded tubes were defined as 'tied' constraints. Table 6 lists the material parameters of the threaded tubes [24]. The threaded tubes and TLE were described using a linear elastic material model. The TLE material parameters for different specifications of the threaded structures are listed in Table 5.
Internal thread tube TLE External thread tube   Table 7 lists the first three order frequencies of the threaded structures obtained in the experiments and with the TLE method. Compared with the experiments, the maximum error of the TLE modeling method is 3.6%, and the minimum is -0.08%. All of the errors are within 5%, which satisfies the precision requirements. In summary, the theoretical determination method for the TLE material parameters is reasonable and reliable for threads with different pitches and thread lengths.

Penetration Test
To further verify the validity and reliability of the TLE method proposed in this study, a test of projectile penetration into a semi-infinite concrete target was performed using a 152 mm caliber one-stage light-gas gun. The projectile was machined from 30CrMnSiNi2A high-strength steel, and the necessary heat treatment was performed. Figure 8 shows the projectile geometry. The projectile had an ogive nose with caliber-radiushead (CRH) 3 and a nominal mass of 4.2 kg. A deceleration measurement device was attached to the projectile tail using a threaded connection for acquisition of the deceleration. An accelerometer and electronic components were sealed inside the test device. The structural material of the device housing should be the same as that of the projectile to avoid the influence of different contact surface stiffnesses on the deceleration signal. The specification of the thread body connecting the projectile and measuring device was M42 × 2-6H/6g-33. The thread pitch and length were 2 and 33 mm, respectively. The sampling rate and trigger threshold of the accelerometer were set at 120 kHz and 8000 g, respectively, and a sleep time of 50 min was set to prevent false triggering of the measurement device. A C40 concrete target plate with a diameter of 120 cm and a thickness of 60 cm was used in the experiment. The measured unconfined compressive strength was 50 MPa. The projectile was accelerated along the gun barrel by the instant emission of highpressure gas. After being separated from the sabots, the projectile impacted the concrete target plate vertically at a high speed of 320 m/s. The projectile target encounter conditions, including the striking velocity and incident attitude angle, were measured using a high-speed camera (Photron/SA5). The deceleration signal was acquired by a deceleration measurement device. The air pressure of the light-gas gun was adjusted to 4.4 MPa according to the projectile mass and striking velocity. After the experiment, the deceleration data in the acceleration data recorder were extracted. Figure 9 presents the deceleration-time curve in the experiment. In Figure 9, the raw experimental signal was low-pass filtered using cut-off frequencies of 40, 20, and 10 kHz. The results indicate that the signal frequency was concentrated within 40 kHz. Therefore, we selected the measured signal spectrum within 40 kHz for further research. Figure 10 shows the fast Fourier-transform (FFT) spectrum of the experimental deceleration signal. The first peak of the curve corresponds to the rigid body deceleration with a frequency of 354 Hz, indicating the response frequency of the projectile-fuse system to the penetration loading source. According to the classic theory of vibration, the dominant frequency of the penetration excitation signal is approximately 354 Hz. The rigid body deceleration was obtained by low-pass filtering of the experimental deceleration with a cut-off frequency of 1500 Hz (Figure 9). The rigid body deceleration represents the real-time resistance signal of the projectile nose. Figure 10 presents the FFT spectrum of rigid body deceleration (red solid line). It can be seen that only one peak exists in the FFT spectrum, with a peak frequency of 354 Hz. This confirms the hypothesis that the frequency of 354 Hz is the dominant frequency of the penetration excitation signal. Moreover, the frequencies at other peaks are roughly integer multiples of the base frequency (7919 Hz), indicating the existence of free vibration with a base frequency of 7919 Hz in the deceleration signal.   Figure 11 shows the high-pass filtering curve of the deceleration signal with a 1500 Hz cut-off frequency. The change in amplitude of the signal is fitted to an exponential function, indicating that the high-frequency region of the signal is free-damping vibration. The decay rate is generally regarded as an essential indicator for evaluating the amplitude attenuation rate, and it can be expressed as follows: where  is the decay rate,  is the logarithmic decay rate, i A and 1 i A  are the amplitudes of two adjacent periods,  is the damping ratio, n  is the natural frequency of the vibration system, and d T is the damping period.
The decay curve in Figure 11 contains five decay cycles, and the damping ratio  is approximately 0.06 according to Equation (8). The damping ratio  of the penetration vibration response has a significant influence on the signal aliasing. When the damping ratio is relatively small, the decay time of the vibration increases. During penetration of the multi-layer hard target, the vibration signal has not yet been completely attenuated. The residual signal adheres to the penetration signal of the next target layer, which affects the identification of the overload signal. The 'NASTRAN/NORMAL MODES' was used to carry out modal analysis of the projectile-fuse system. The mode shape near a frequency of 7919 Hz is shown in Figure 11. The vibration form is first-order axial vibration, and the frequency is 8264.2 Hz. Because the FE software cannot deal with the nonlinear factors of the interface contact reasonably, the calculation result will be overestimated. The experimental deceleration signal was band-pass filtered with a cut-off frequency of 6000-9000 Hz, as shown in Figure 11 (blue solid line). The appearance of a complete sine wave indicates the existence of a base frequency. As shown in Figure 11, the first-order axial vibration of the sine wave ceases before the end of penetration, and the end time is approximately 0.8 ms. This indicates that the base frequency is not the dominant frequency of the signal, which can also be seen in the FFT spectrum ( Figure 10). The intersection of the damping curve and centerline (0 g) is defined as the endpoint of the vibration response region. It can be seen that the vibration response terminates at approximately 1.4 ms. In Figure 9, the vibration amplitude superimposed on the rigid-body deceleration curve is significantly reduced after 1.4 ms, indicating that the vibration response has ceased. First order axial vibration Figure 11. High-pass filtering curve of the deceleration signal with a 1500 Hz cut-off frequency (source elaborated by authors).

Determination of the TLE Material Parameters for the Projectile-Fuse Thread
The threaded body between the projectile and the fuse comprises an ISO metric triangle thread, and the specification is M42 × 2-6H/6g-33. The materials of the bolt body and nut body are the same as those of the projectile, and the specific parameters are listed in Table 8. Figure 12 shows the influence of the thread friction coefficient on the axial stiffness and TLE shear modulus. As shown in Figure 12, the thread axial stiffness and shear modulus gradually increase as the friction coefficient increases. The thread contact friction coefficient [30] is set to 0.08. The resulting thread connection stiffness and material parameters of the TLE are listed in Table 9. Only the axial modes are excited during the penetration event. Thus, the axial stiffness associated with the shear modulus should be considered, and the radial stiffness associated with the Young's modulus can be ignored. The Young's modulus of the TLE can be determined directly based on the relationship between the shear modulus and Young's modulus.

Nodal Force Distribution of the Projectile Nose
To simulate the real penetration excitation environment, the penetration excitation force should be accurately applied to the projectile nose. The projectile has a symmetrical structure. Therefore, the radial forces cancel each other out for vertical penetration, and only the excitation force in the axial direction should be considered. In the frequency response analysis of the projectile-fuse system, the axial penetration resistance was usually evenly distributed to the nodes [31]. However, this cannot accurately represent the excitation of the projectile nose and affects the dominant frequency distribution. Thus, we propose to assign the axial penetration resistance to the nodes of the projectile nose according to the node position. Figure 13 shows the force analysis on the oval nose of the projectile. Forrestal and Luk [32] proposed the relationship between the radial stress on the cavity surface and the cavity expansion velocity based on the cavity expansion theory. The relationship between the radial stress on the projectile surface and the rigid body velocity is expressed as: where n  is the radial stress at the projectile surface, ' c f is the unconfined compressive strength of the concrete target, z V is the rigid body velocity, and  is the target density.
A and B are target material coefficients, which can be obtained from triaxial tests of the target material. The value of B is generally within a small range around 1.0 [32,33]. For concrete materials, B is approximately 1.0, and A (A = S) is related to the compressive strength of the concrete material [34]. Based on many experiments, Li and Chen [35] obtained the relationship between S and ' c f as  For an area of ds at the projectile nose surface, the axial resistance caused by the radial stress on the projectile surface can be expressed as: From Equation (10), it can be seen that the axial resistance at any point is a thirdorder function of sin  . To investigate the relationship between The projectile was meshed using the mapped method. The numbers of circumferential and axial nodes are s and n, respectively ( Figure 13).
where F1 is the axial resistance of the node at the tip of the projectile, Fn is the axial resistance of the node at the arc-column interface of the projectile, F is the rigid body resistance in the axial direction, p m is the projectile mass, and a is the rigid body deceleration.
From Equation (10), the axial resistance of any node at the projectile nose can be expressed as follows:

TLE Modeling
The NASTRAN/FREQUENCY RESPONSE software was used to simulate the penetration response of the projectile-fuse system. The projectile body and fuse body are connected by threads. The vibration response of the fuse cannot be accurately predicted because of the uncertain factors at the threaded contact. The primary uncertainty in threaded connections arises from the contact stiffness and friction damping [30]. Because the commonly used dynamic analysis software cannot correctly deal with nonlinear factors, the TLE was used to characterize the contact stiffness between threaded contacts. In addition, the thread friction damping was described by setting a reasonable damping coefficient for the TLE. According to the height of the thread tooth, the thickness of the TLE was determined to be 2 mm. The length of the TLE was 33 mm, i.e., the threaded length. Figure 15 shows the FE model of the projectile-fuse system. The FE model was composed of three parts: the projectile body, fuse body, and TLE. All of the entities were modeled using eight-node hexahedral elements (CHEXA), and each node had six degrees of freedom. For the convenience of constraint treatment at the contact interface, the element sizes of the projectile body, fuse body, and TLE should be consistent at the thread contact interface. The projectile and fuse were characterized using a linear elastic model. The material parameters are listed in Table 8. The TLE was also characterized by a linear elastic material model, and Table 9 lists the material parameters. The contacts between the two pairs (TLE and projectile, TLE and fuse) were defined as 'tied' constraints. At the projectile nose, the numbers of circumferential and axial nodes were 120 and 44, respectively. To satisfy the quality requirements of the hexahedral mesh, the number of circumferential nodes at the tip of the projectile nose was generally reduced during mesh generation. Thus, we allocated the penetration excitation forces by combining the node position with the number of circumferential nodes. In the present simulation, F in Equation (12) refers to the experimental rigid-body deceleration. The magnitudes of each nodal force were obtained using Equation (12). The projectile had a free state without any boundary constraints.
In addition, the damping was incorporated by the complex stiffness matrix * K domain with dissipation multipliers αi and βi for the material and joint damping, respectively, and the associated element stiffness matrices.
where K is a real-value stiffness matrix, and n and m are the number of different materials and joints with different properties used in the model, respectively. Material damping is applied to the elements belonging to the components of the projectile-fuse system, and joint damping is realized by applying damping to the TLEs [13]. In the present calculation, the material damping α is 5 × 10 −4 , and the joint damping β is 0.03, determined through a trial-and-error method. The Lanczos algorithm for the modal superposition method was applied for analysis and solution. The number of desired roots was set to 100. To meet the 40 kHz output frequency range, the length and number of output frequency steps were 100 Hz and 400 respectively. One node at the middle partition (accelerometer position) of the fuse was selected as the frequency response point (Figure 15). Figure 16 shows a comparison of the nodal force distribution methods, including an equal distribution and the distribution according to Equation (12). The frequency response resonance peaks of the nodal force distribution using Equation (12) are concentrated in the frequency band of 30-40 kHz, while the resonance peaks of the equal distribution method are concentrated around the first-order frequency. The concentration region of the resonance peaks for the distribution method according to Equation (12) is consistent with the experimental FFT spectrum. In addition, the amplitude variation of the resonance peaks exhibits a similar trend as that obtained in the experiment. Thus, the nodal force distribution method using Equation (12) is more reasonable. Figure 17 shows the FRF curve of the TLE and TIED modeling methods. Compared to the experimental FFT spectrum, it can be seen that the two modeling methods have a similar frequency distribution as the experiment, and only at the frequency of 20,284 Hz are the resonance peaks obtained by the two simulation methods not noticeable. This is because the vibration at the frequency of 20,284 Hz is a non-axial vibration mode. Only the nodal force along the axial direction was applied in the FE model. Therefore, this vibration mode was not completely excited. During penetration, it is impossible to guarantee the complete vertical penetration of the projectile. Thus, the resonance peak (20,284 Hz) appears in the experimental spectrum. Compared with the 'TIED' method, the amplitudes of the FRF curve at each frequency obtained with the TLE method are closer to those of the experimental spectrum. Table 10 lists the resonance frequencies of the TLE method, TIED method, and experiment. Each vibration frequency obtained with the TLE method appears closer to the experimental value. The maximum error is 7.8%, indicating that the TLE method can accurately simulate the contact interface of the threaded connection by defining a reasonable damping coefficient of the TLE.

Conclusions
The TLE method was proposed to represent the contact stiffness and friction damping of the threaded connection of a projectile-fuse system. Combining the thread elastic model and TLE principle, a method for determining the TLE material parameters was established. The accuracy and reliability of the material parameters were verified based on modal experiments of threaded tubes. The TLE modeling method was applied to evaluate the vibration response characteristics of the projectile-fuse system. The simulation results were in good agreement with the experimental results for a projectile penetrating a semi-infinite concrete target. The results can be summarized as follows: (1) Several modal tests of threaded tubes were performed to investigate the nonlinear response characteristics of threaded connections. The results indicate that the frictional damping and contact stiffness had a significant impact on the dynamic characteristics, causing the resonance frequency and resonance peak to decrease.
(2) An effective method for determining the TLE material parameters was established based on the thread elastic model and TLE principle. Threaded tubes with different specifications were used to verify the accuracy and reliability of the material parameters. The errors were within 5%, indicating the accuracy of the material parameters. (3) The experimental deceleration signal mainly consisted of the shock response and vibration response. The shock response frequency depended on the frequency of the penetration excitation signal. The vibration response was a damped vibration with the first-order axial vibration frequency as the base frequency. (4) A reasonable force distribution method for a projectile nose node was proposed.
Compared with the equal distribution method, the proposed method could better simulate the penetration excitation force. This study thus provides a reference for the reasonable distribution of the penetration excitation force in dynamic analyses of projectile-fuse systems. (5) Frequency response analysis of the projectile-fuse system was performed using the TLE. The maximum error in the resonant frequencies between the TLE method and experiment was 7.8%, whereas the maximum error between the existing 'TIED' constraint method and the experiment was 15.7%. In addition, the resonant peak distribution obtained by the TLE method was in better agreement with the frequency distribution of the experimental deceleration signal.