Effect of Wood Medium on Dispersion Parameters of Prefabricated Spherical Fragments in Forest

: In order to study the inﬂuence of the forest environment on the killing power of pre-fabricated fragments, the relationship between the wood medium in the forest and the residual velocity of prefabricated fragments and a simple calculation model were obtained. The impact test of 300 mm thick pine target plate was carried out with 6 mm and 11 mm diameter tungsten alloy spherical fragments. The LS-DYNA ﬁnite element modeling of wood was carried out and compared with the experimental parameters, which veriﬁed the reliability of the numerical simulation method on the residual velocity of fragments. Based on this, a simple mathematical calculation method for the forest environment is constructed to calculate the penetration track length of fragments in the forest environment. Combined with the model and the formula for calculating the residual velocity of wood, the impact of forests on fragments within a certain radius can be evaluated. According to the characteristics of the forest environment, the boundary effects of the multi-layer spacing targets and trees were further studied. The research shows that the wood reduces the fragment power mainly by affecting the penetration length. The inﬂuence of forest density, tree diameter, and other parameters on the fragment velocity attenuation in the forest environment is analyzed. Using this method, the inﬂuence of forest environmental parameters on fragment dispersion parameters can be simply evaluated.


Introduction
As a natural environment widely distributed on the land, the latitude and geographical environment make the species of trees very rich, and their density changes greatly, generally between 0.4~0.8 g/cm 3 . In a battlefield environment, the forest is generally used for military purposes such as concealing troop movements and military equipment. When the kill-blasting warhead acts on the forest environment, fragments may penetrate trees in the process of dispersion. Trees affect the flying motion of fragments, resulting in the reduction of flying distance and residual speed of fragments. As a factor to be considered for all types of weapons in the battlefield environment, environmental media can generally be divided into natural environment and artificial environment. The natural environment includes soil, trees, air, water, etc., and the artificial environment includes concrete materials, metal materials, etc. For the penetration study of metal target and concrete targets, Xu et al. [1] studied the failure behavior and damage mechanism of tungsten alloy fragments penetrating low carbon steel targets at high speed. The impact of a shaped charge jet on the penetration, performance, and penetration mechanism of different metals was studied, and the penetration depths were compared [2]. It can be seen that because tanks and armored vehicles and other military equipment use metal protection, the penetration of metal and composite targets was studied in depth. Similarly, there have been many penetration theories and calculation formulas about concrete in the attack on various types of fortifications. Fu, H. [3] proposed a formula based on a penetration theory of thick concrete, an acceleration theory, and a model of projectile penetration of concrete. The calculated values were found to be consistent with the experimental values. Dong, YX [4] used finite element software to analyze the ballistic characteristics of flat-nosed projectiles on concrete, their penetration depth, residual velocity, and other indicators.
At present, in the research of projectile penetration, the target plates are mostly reinforced concrete and metal targets; there is less research on wood-medium penetration and high-speed impact, and the research is mostly about the mechanical properties of wood materials, such as fatigue and low-speed impact. The impact resistance of three kinds of wood materials was studied with experiments, and the ultimate speed and energy absorption of different materials were obtained [5]. By changing the types of wood in the sandwich panel, the damage difference of the sandwich panel under different impact energies was studied, and the influence of the sandwich panel on resisting damage was analyzed [6]. The mechanical behavior of palm wood particle-reinforced high-density polyethylene under high strain-rate compression and ballistic impact load was investigated, and it was found that the mechanical properties can be enhanced [7]. The damage mode of the wood sandwich panel under low velocity impact was studied, and it was found that a light wood core panel could have a greater impact load [8]. Due to the anisotropy and low density of wood target plate, and the fact that its structure is fibrous material, which is different from most man-made materials, its structural characteristics and strength are special, and its damage state is different when impacted [9]. The material density, tree diameter, skin thickness, inner core thickness, and other parameters of different trees are also different. In addition, the impact of wood moisture content in the natural environment may have different effects on the penetration process of fragments. The long-term bearing capacity of wood is far lower than its temporary bearing capacity. This is because under long-term load-bearing conditions, the wood will experience fiber creep, which will cause large deformation after accumulation, thus reducing the wood-bearing capacity [10]. Therefore, it is necessary to analyze the anti-destruction ability of wood materials under instantaneous impact.
In order to study the penetration mechanism of high-speed fragments into wood and their influence on the residual velocity of fragments, the laws of forest density, tree diameter and other parameters on the residual velocity and killing radius of fragments were obtained, and the penetration mechanism of fragments impacting wood materials and the formula calculation method of penetration resistance were analyzed and calculated. The method of comparison between experiment and simulation is used to compare the failure mode and residual velocity of wood materials, and the method of comparison between modeling calculation and theoretical calculation is used to verify the reliability of theoretical calculation value of average penetration length. Through the comparison between the actual value and several residual velocity formulas, the calculation model of the residual velocity of light wood materials was obtained.

Comparison of Penetration Experiment and Simulation
The 12.7 mm ballistic gun was used for the test of fragment penetration into a wood medium. Three pine target plates, with parameters of 300 mm × 300 mm × 100 mm, were placed in the target box, and two velocity measuring targets were set in front of and behind the target box and connected with a six-channel velocity meter to measure the initial velocity and residual velocity of the fragments. Figures 1 and 2 show the schematic diagram and physical diagram of the fragment penetration test. Appl. Sci. 2022, 12,   The ballistic gun test was carried out on 11 mm and 6 mm tungsten alloy spherical fragments using the above device. In Figure 2, the left side is the fragment, cartridge holder and cartridge, and the right side ① is the front velocity target; ② This is a rearspeed measuring target; ③ This is a box for placing three wood target plates; ④ In order to connect the wires of the tachometer, the on-off test method was adopted for the test speed measurement. Different initial velocities were obtained by controlling the amount of charge in the cartridge.
The penetration test of tungsten alloy fragments into the wood target was conducted 16 times in the speed range of 589.55 m/s~1363.24 m/s. Table 1 shows the initial and residual velocities of the fragments, where 0 represents no penetration.  The ballistic gun test was carried out on 11 mm and 6 mm tungsten alloy spherical fragments using the above device. In Figure 2, the left side is the fragment, cartridge holder and cartridge, and the right side ① is the front velocity target; ② This is a rearspeed measuring target; ③ This is a box for placing three wood target plates; ④ In order to connect the wires of the tachometer, the on-off test method was adopted for the test speed measurement. Different initial velocities were obtained by controlling the amount of charge in the cartridge.
The penetration test of tungsten alloy fragments into the wood target was conducted 16 times in the speed range of 589.55 m/s~1363.24 m/s. Table 1 shows the initial and residual velocities of the fragments, where 0 represents no penetration.  The ballistic gun test was carried out on 11 mm and 6 mm tungsten alloy spherical fragments using the above device. In Figure 2, the left side is the fragment, cartridge holder and cartridge, and the right side 1 is the front velocity target; 2 This is a rear-speed measuring target; 3 This is a box for placing three wood target plates; 4 In order to connect the wires of the tachometer, the on-off test method was adopted for the test speed measurement. Different initial velocities were obtained by controlling the amount of charge in the cartridge.
The penetration test of tungsten alloy fragments into the wood target was conducted 16 times in the speed range of 589.55 m/s~1363.24 m/s. Table 1 shows the initial and residual velocities of the fragments, where 0 represents no penetration. It can be seen from the damaged state of the wood in Figure 3 that the front-end face of the wood leaves a hole with a diameter similar to the fragment after penetration, and there is no obvious damage expansion effect around the hole; The rear-end face is basically the same as the front-end face. Some wood chips are brought out at the mouth, and some wood is peeled off along the axial surface of the wood. It can be seen that when the wood was penetrated by small caliber projectiles, the removal of the penetration trajectory did not damage the additional part of the wood, and the trajectory had no obvious deviation. Therefore, the damage of tungsten alloy spherical fragments to the wood was only a tunnel with a diameter similar to that of the fragments. Appl It can be seen from the damaged state of the wood in Figure 3 that the front-end face of the wood leaves a hole with a diameter similar to the fragment after penetration, and there is no obvious damage expansion effect around the hole; The rear-end face is basically the same as the front-end face. Some wood chips are brought out at the mouth, and some wood is peeled off along the axial surface of the wood. It can be seen that when the wood was penetrated by small caliber projectiles, the removal of the penetration trajectory did not damage the additional part of the wood, and the trajectory had no obvious deviation. Therefore, the damage of tungsten alloy spherical fragments to the wood was only a tunnel with a diameter similar to that of the fragments. Using mat in LS-DYNA MAT_WOOD material parameters, which include parameters of various strength indexes of wood, can also add erosion contact, and can be modified according to different strength indexes in radial and axial directions caused by different wood anisotropy. The material parameters of wood mainly include density, elastic Using mat in LS-DYNA MAT_WOOD material parameters, which include parameters of various strength indexes of wood, can also add erosion contact, and can be modified according to different strength indexes in radial and axial directions caused by different wood anisotropy. The material parameters of wood mainly include density, elastic modulus, shear modulus, tensile strength, compressive strength and shear strength, and other parameters are shown in Table 2. In Table 2, ρ is the density of wood, E // is the parallel modulus of elasticity, E ⊥ is the vertical modulus of elasticity, σ st is the parallel tensile yield stress, σ sc// is the parallel compression yield stress, σ st⊥ is the vertical tensile yield stress, σ sc⊥ is the vertical compression yield stress, and µ is Poisson's ratio.
From the photos of the wood target plate, it can be seen that the cutting plane is perpendicular to the radial direction, which is consistent with the state of most fragments in the actual situation. Based on this, a model with the same parameters as the experiment was established in the simulation, and the residual velocity was taken as the main parameter for comparison, as shown in Table 3. From the comparison of simulation and experimental parameters in Table 3, it can be seen that the material model of wood has a high degree of explanation for residual velocity and failure mode, wherein the maximum error of residual velocity is 22.18% and the minimum error is 0.14%. The failure mode in the simulation is the same as that in the experiment. A little wood flies out on the back of the target plate, as shown in Figure 4.
From the comparison of simulation and experimental parameters in Table 3, it can be seen that the material model of wood has a high degree of explanation for residual velocity and failure mode, wherein the maximum error of residual velocity is 22.18% and the minimum error is 0.14%. The failure mode in the simulation is the same as that in the experiment. A little wood flies out on the back of the target plate, as shown in Figure 4.

Fragment Penetration Calculation Based on Forest Environment
When a warhead is detonated in the forest, various parameters of the warhead and the forest will cause different initial variables of fragment scattering, such as the height H of warhead initiation, the size and direction v of fragment velocity vector, the main wood species caused by forest environment, tree diameter D, and forest density [11] ρ, etc. Among them, due to the difference of tree species, the strength and density of bark and inner core of trees are different, and the thickness of bark is also different [12], so the difference between bark and inner core is ignored in this study. The initial velocity vector of the fragment is decomposed into a vector v xy parallel to the ground and a vector v z perpendicular to the ground, where v xy mainly controls the flight direction of the fragment in the horizontal direction; v z mainly controls the flight time of fragments in the vertical direction. It can be seen that v z and H determine the flight time of fragments in the air, and then affect the flight distance of fragments. See Figure 5 for the fragment penetrating trees and the relationship between the variable parameters of fragments.

Fragment Penetration Calculation Based on Forest Environment
When a warhead is detonated in the forest, various parameters of the warhead and the forest will cause different initial variables of fragment scattering, such as the height H of warhead initiation, the size and direction v of fragment velocity vector, the main wood species caused by forest environment, tree diameter D, and forest density [11] ρ, etc. Among them, due to the difference of tree species, the strength and density of bark and inner core of trees are different, and the thickness of bark is also different [12], so the difference between bark and inner core is ignored in this study. The initial velocity vector of the fragment is decomposed into a vector vxy parallel to the ground and a vector vz perpendicular to the ground, where vxy mainly controls the flight direction of the fragment in the horizontal direction; vz mainly controls the flight time of fragments in the vertical direction. It can be seen that vz and H determine the flight time of fragments in the air, and then affect the flight distance of fragments. See Figure 5 for the fragment penetrating trees and the relationship between the variable parameters of fragments. In the x-y plane, the fragments may penetrate the wood in the forest in all directions, based on the forest density ρ and tree diameter D. assuming that the trees are evenly distributed in the forest, the number of trees per unit area can be determined according to the stand density (n/m 2 ), so as to determine the land area of each tree. If each tree is placed in the center of its occupied area, the following forest average distribution diagram can be obtained.
Taking the 1/4 model as an example, the lower left corner is the explosion center, the stand density is 1/m 2 , and the tree diameter is 300 mm. From the top view, the trajectory In the x-y plane, the fragments may penetrate the wood in the forest in all directions, based on the forest density ρ and tree diameter D. assuming that the trees are evenly distributed in the forest, the number of trees per unit area can be determined according to the stand density (n/m 2 ), so as to determine the land area of each tree. If each tree is placed in the center of its occupied area, the following forest average distribution diagram can be obtained.
Taking the 1/4 model as an example, the lower left corner is the explosion center, the stand density is 1/m 2 , and the tree diameter is 300 mm. From the top view, the trajectory of the fragments will fly around without considering the vertical movement time until v xy is 0. By numbering the fragment trajectories at 3 • intervals, it can be seen that the lengths of different fragments penetrating the wood in the same radius area are different, some through four trees, and some through only one tree. Although this difference can reflect the maximum damage radius and the minimum damage radius of fragments in the forest environment, individual differences cannot reflect the comprehensive damage degree of all fragments. The trajectory of fragments penetrating trees is shown in Figure 6.  In this paper, a method of average distribution of wood area is adopted. The penetration trajectory area of fragments in the x-y plane is equivalent to the product of fragment diameter d and trajectory length L. The wood area S penetrated in the penetration area can be calculated through the trajectory area it passes through. Ignoring the irregular part of the tree penetration surface, the path length L of the wood medium can be obtained by dividing the penetration area S by fragment diameter d, as shown in formulas (1) to (3).
The calculation formula of the average area of trees is: Among ρ is the density of trees per unit area (n/m 2 ), and D is the average diameter of trees.
The penetration area of the fragment with diameter d on the horizontal plane is: R is the radius length of fragment flight. The average penetration path of fragments is: In the above formula, the following table p represents penetration and w represents wood.
In order to verify the accuracy of the calculation method, the above method is used to verify different parameters in the following three cases: In this paper, a method of average distribution of wood area is adopted. The penetration trajectory area of fragments in the x-y plane is equivalent to the product of fragment diameter d and trajectory length L. The wood area S penetrated in the penetration area can be calculated through the trajectory area it passes through. Ignoring the irregular part of the tree penetration surface, the path length L of the wood medium can be obtained by dividing the penetration area S by fragment diameter d, as shown in formulas (1) to (3).
The calculation formula of the average area of trees is: Among ρ is the density of trees per unit area (n/m 2 ), and D is the average diameter of trees.
The penetration area of the fragment with diameter d on the horizontal plane is: R is the radius length of fragment flight. The average penetration path of fragments is: In the above formula, the following table p represents penetration and w represents wood. In order to verify the accuracy of the calculation method, the above method is used to verify different parameters in the following three cases:

1.
Set the wood diameter and stand density as fixed values, change the horizontal dispersion angle between two adjacent fragments, and compare the average penetration trajectory length of fragments with formula (3).
Set different wood diameters without changing the stand density and scattering angle, and compare the average penetration trajectory length of fragments with formula (3).
The average penetration trajectory length of fragments is compared with formula (3) by setting different stand densities without changing the tree diameter and scattering angle. (1/m 2 , 0.5/m 2 , 0.25/m 2 ) The comparison results are shown in Table 4. It can be seen from the error between the actual value and the calculated value corresponding to the three variables of different scatter angles, tree diameters, and stand densities that the calculated value was generally larger than the actual value, and the error rate increased with the increase in the scatter angle. It decreased with the increase in tree diameter, and with the increase in stand density. It can be seen that the accuracy of the simple calculation formula was closely related to the three variables. Therefore, on the one hand, the appropriate range of use of each variable could be obtained to improve the accuracy. On the other hand, the correction coefficient w for the three variables could be added to Equation (3). The dimensionless expression of w is: Using the penetration area of the fragment, the length of the fragment penetrating the wood in the horizontal direction can be obtained, that is, the thickness of the target plate of the fragment penetrating the wood. The overall penetration path L xyz can be obtained by using the penetration length L xy of the horizontal plane and the angle β between the fragment velocity vector and the horizontal plane, as shown in formula (5).
According to the experimental results of 11 mm tungsten alloy fragments and 6 mm fragments, the relationship between the residual velocity and the final velocity is an exponential function when the thickness H of the target plate is constant, which is similar to the velocity attenuation law of the projectile in the air.
Due to the low density and low quality of wood, the velocity attenuation formula of fragments in air is considered first [13]. When fragments fly in the air, they are affected by gravity and air resistance, and gravity makes the flight trajectory of fragments bend. Air resistance causes the attenuation of fragment velocity. Since the distance from the fragment to the target is not too long and the time is very short, the influence of gravity can be ignored, and the fragment trajectory can be approximated as a straight line.
The differential equation of fragment motion is: where, C D is the aerodynamic resistance coefficient; ρ a is the air density; S is the windward display area of fragments; m v and t are fragment mass, fragment velocity and time respectively. The above formula is transformed into the differential form of velocity v and distance x: To integrate the above formula, there are: It can be seen from the above formula that the relationship between v x and v 0 should be linear when the thickness x of the target plate and the mass m of the projectile are constant. However, it can be seen from the simulation results that the growth of v x is nonlinear with the increase in v 0 , as shown in Figure 7. Therefore, the velocity attenuation formula in air cannot be applied to the case of penetrating wood.
Appl. Sci. 2022, 12, x FOR PEER REVIEW 9 of 14 To integrate the above formula, there are: So we get: It can be seen from the above formula that the relationship between vx and v0 should be linear when the thickness x of the target plate and the mass m of the projectile are constant. However, it can be seen from the simulation results that the growth of vx is nonlinear with the increase in v0, as shown in Figure 7. Therefore, the velocity attenuation formula in air cannot be applied to the case of penetrating wood. According to the hypothesis of Poncelet [14], the equation of penetration depth of a stable moving projectile varying with velocity can be obtained by integration: The result is: where s represents the distance along the straight trajectory. The speed after passing through the obstacle with thickness h is: where v is the residual velocity; v0 is the initial speed; c1 and c3 are constants; A is the crosssectional area of the projectile; h is the thickness of target plate; and m is the fragment mass. Since the prefabricated fragments are spherical, the above formula can be simplified According to the hypothesis of Poncelet [14], the equation of penetration depth of a stable moving projectile varying with velocity can be obtained by integration: The result is: where s represents the distance along the straight trajectory. The speed after passing through the obstacle with thickness h is: where v is the residual velocity; v 0 is the initial speed; c 1 and c 3 are constants; A is the crosssectional area of the projectile; h is the thickness of target plate; and m is the fragment mass.
Since the prefabricated fragments are spherical, the above formula can be simplified as: r is the radius of the fragment. According to the Poncelet formula, the difference between the two parameters c 1 and c 3 obtained from the initial and residual velocity data of 11 mm and 6 mm projectiles is too large, and the difference between c 1 and c 3 is more obvious. This may be caused by experimental errors, and there is a problem of an insufficient number of samples. Therefore, the simulation values of the two fragments are fitted to obtain the following relationship between v 2 and v 0 2 , as shown in Figure 8. r is the radius of the fragment. According to the Poncelet formula, the difference between the two parameters c1 and c3 obtained from the initial and residual velocity data of 11 mm and 6 mm projectiles is too large, and the difference between c1 and c3 is more obvious. This may be caused by experimental errors, and there is a problem of an insufficient number of samples. Therefore, the simulation values of the two fragments are fitted to obtain the following relationship between v 2 and v0 2 , as shown in Figure 8. From the results of linear fitting, it can be seen that there is an obvious linear relationship between the two, so the c3 value of the two fragments can be calculated. Through comparison, it is found that the difference ratio between the two is not large, so the c3 parameter value of the dried pine wood is 325; However, the c1 value obtained by the Poncelet formula is too different and cannot be regarded as consistent. This should be the result of fragment parameters. Therefore, the fragment radius r in the latter half of the formula is modified and r 2.2 is added to the denominator term and substituted into the latter half of the formula. At this time, the c1 value obtained is −2.6. Thus, a modified formula about the residual velocity of wood based on the Poncelet formula is obtained: After the residual velocity calculation formula is obtained, the influence of forest on fragment velocity can be obtained according to the scattering radius, stand density, and tree diameter.

Simulation Comparison of Different Penetration Modes
Based on the material parameters of wood, this chapter expands and analyzes the influence of the edge effect of wood medium and the number of interval target layers on the residual velocity.

Edge Effect of Wood Penetration
From the fragment scattering trajectory shown in Figure 6, it can be seen that there are a large number of fragment trajectories passing through the edge of the tree, that is, from the cross-section of the top view, the chord tangent length of the trajectory is very short, and the wood located outside the trajectory is very low in strength, so it is more easily damaged. At this time, it is obviously unscientific to compare the length of the penetration trajectory to the thickness of the target plate [15,16]. Therefore, based on the wood parameter in LS-DYNA, the position of velocity penetration was changed, so as to change the thickness of trees outside the track. The track length in each case was measured and From the results of linear fitting, it can be seen that there is an obvious linear relationship between the two, so the c 3 value of the two fragments can be calculated. Through comparison, it is found that the difference ratio between the two is not large, so the c 3 parameter value of the dried pine wood is 325; However, the c 1 value obtained by the Poncelet formula is too different and cannot be regarded as consistent. This should be the result of fragment parameters. Therefore, the fragment radius r in the latter half of the formula is modified and r 2.2 is added to the denominator term and substituted into the latter half of the formula. At this time, the c 1 value obtained is −2.6. Thus, a modified formula about the residual velocity of wood based on the Poncelet formula is obtained: After the residual velocity calculation formula is obtained, the influence of forest on fragment velocity can be obtained according to the scattering radius, stand density, and tree diameter.

Simulation Comparison of Different Penetration Modes
Based on the material parameters of wood, this chapter expands and analyzes the influence of the edge effect of wood medium and the number of interval target layers on the residual velocity.

Edge Effect of Wood Penetration
From the fragment scattering trajectory shown in Figure 6, it can be seen that there are a large number of fragment trajectories passing through the edge of the tree, that is, from the cross-section of the top view, the chord tangent length of the trajectory is very short, and the wood located outside the trajectory is very low in strength, so it is more easily damaged. At this time, it is obviously unscientific to compare the length of the penetration trajectory to the thickness of the target plate [15,16]. Therefore, based on the wood parameter in LS-DYNA, the position of velocity penetration was changed, so as to change the thickness of trees outside the track. The track length in each case was measured and compared with the remaining velocity under the target plate with the same thickness, and the relationship between the outside thickness and the phase difference ratio was observed. See Table 5 for the comparison of the two cases.  Figure 9 is a schematic diagram of the relationship between the difference ratio values of the two cases corresponding to different outer thickness. With the increase in the outer thickness of the trees, the residual velocity of penetrating the trees is always greater than the residual velocity of the normal target plate with the same thickness, and the difference value generally shows a trend that fluctuates upwards. However, it should also be seen that the maximum difference is only 19.2 m/s, so it can be considered that the edge effect has little influence on the high-speed impact fragments.
Appl. Sci. 2022, 12, x FOR PEER REVIEW 11 of 14 the relationship between the outside thickness and the phase difference ratio was observed. See Table 5 for the comparison of the two cases.  Figure 9 is a schematic diagram of the relationship between the difference ratio values of the two cases corresponding to different outer thickness. With the increase in the outer thickness of the trees, the residual velocity of penetrating the trees is always greater than the residual velocity of the normal target plate with the same thickness, and the difference value generally shows a trend that fluctuates upwards. However, it should also be seen that the maximum difference is only 19.2 m/s, so it can be considered that the edge effect has little influence on the high-speed impact fragments. In conclusion, due to the edge effect of wood penetration, the impact can not be ignored when the outer thickness d is greater than a certain value. Therefore, the penetration trajectory of fragments to trees is still an important basis for the residual velocity of fragments. In conclusion, due to the edge effect of wood penetration, the impact can not be ignored when the outer thickness d is greater than a certain value. Therefore, the pene-tration trajectory of fragments to trees is still an important basis for the residual velocity of fragments.

Influence of the Number of Target Layers on the Residual Velocity
Due to the particularity of the forest environment, the penetration factors in the target area are affected by the tree parameters. The most significant of these is the large number of trees on the fragment penetration trajectory caused by the large forest density, which makes the fragments penetrate multi-layer spaced targets. It can be seen from the above inference that the penetration type in the forest environment is generally interval target, and the relationship between the thickness of interval target and the overall target needs to be further calculated [17][18][19][20]. In order to analyze the influence of the number of spaced target layers on the residual velocity under the same target plate thickness, the initial velocity of spherical fragments is set to be 1000 m/s, the fragment diameter is 10 mm, and the total thickness of wood target plates is 500 mm. The target plate layers are divided into six types: 1, 3, 5, 7, 9, and 11. The residual velocity of fragments after penetration is counted separately.
It can be seen from the analysis results in Table 6 that with the increase in the number of target layers, the residual velocity also increases, that is, the number of wood target plates weakens the impact of the thickness of the target plate, and the smaller the thickness of the target plate, the smaller the impact resistance, so that the residual velocity of fragments increases continuously [21,22]. When the number of target plates is 11, the difference ratio reaches 23.89%.

Conclusions
Through the penetration test of tungsten alloy spherical fragments into wood targets and the construction of wood parameters in finite element analysis, we achieved a preliminary understanding of the calculation of residual velocity of fragments penetrating wood. By analyzing the relationship between the initial parameters of fragments, residual velocity, stand density and average diameter of trees in the forest, the Poncelet formula was used to obtain the calculation formula of residual velocity, and a calculation model for fragment track length was proposed. The main conclusions of this paper are as follows: 1.
Based on the Poncelet formula, it is feasible to calculate the residual velocity of the spherical fragment penetrating the wood target, which is verified by the comparison with the experiment.

2.
Based on the forest parameters such as stand density, tree diameter, and fragment flying angle, a simple calculation model of penetration path length is obtained, which has good applicability when the flying angle is small and the stand density is large.

3.
The residual velocity increases with the increase in the number of spacer targets on the fragment penetration trajectory, and it changes most obviously when the number of layers is seven. The wood target is less affected by the boundary effect. When the thickness d outside the penetration trajectory decreases continuously, the influence of the wood boundary effect on the residual velocity is relatively stable, and the maximum difference is 2.67%.