Numerical and Theoretical Analysis of the Inertia Effects and Interfacial Friction in SHPB Test Systems

The dynamic properties of materials should be analyzed for the material selection and safety design of robots used in the army and other protective structural applications. Split Hopkinson pressure bars (SHPB) is a widely used system for measuring the dynamic behavior of materials between 102 and 104 s−1 strain rates. In order to obtain accurate dynamic parameters of materials, the influences of friction and inertia should be considered in the SHPB tests. In this study, the effects of the friction conditions, specimen shape, and specimen configuration on the SHPB results are numerically investigated for rate-independent material, rate-dependent elastic-plastic material, and rate-dependent visco-elastic material. High-strength steel DP500 and polymethylmethacrylate are the representative materials for the latter two materials. The rate-independent material used the same elastic modulus and hardening modulus as the rate-dependent visco-elastic material but without strain rate effects for comparison. The impact velocities were 3 and 10 m/s. The results show that friction and inertia can produce a significant increase in the flow stress, and their effects are affected by impact velocities. Rate-dependent visco-elasticity material specimen is the most sensitive material to friction and inertia effects among these three materials (rate-independent material, rate-dependent elastic-plastic material, and rate-dependent visco-elastic material). A theoretical analysis based on the conservation of energy is conducted to quantitatively analyze the relationship between the stress measured in the specimen and friction as well as inertia effects. Furthermore, the methods to reduce the influence of friction and inertia effects on the experimental results are further analyzed.


Introduction
Designers, engineers, and manufacturers are always in search of new and better materials for performance improvements and cost reduction in their robot products to remain competitive in the market. This requirement has resulted in the limited selection of thousands of material choices for designers. Metal, composite material, and polymer are still the main choices for the structural design of robots, especially for their bearing components. In engineering applications, quadruped bionic robots [1], exoskeleton robots [2], and humanoid robots [3] are used in the army and other protective structural applications. Dynamic impact incidents-e.g., the falling of maintenance tools, unintentionally trampling, gravel impact, and hail impact-often occur, which is dangerous to the service safety of robots. Hence, it is essential to investigate the dynamic properties of different materials for a proper selection of materials and robot component safety design.
The split Hopkinson pressure bar (SHPB) is a widely used experimental system to investigate the dynamic compression properties of materials between strain rates of 10 2 and 10 4 s −1 [4,5]. There are many factors that can influence the accuracy of SHPB experimental results, such as the material response, test conditions and data processing [6][7][8]. Experimental results from SHPB experiments are the stress waves of materials under different strain rates [9,10]. A measurement error of less than 4% dynamic flow stress increment is needed if the strain rate effects are to be reliably characterized, and all sources of error must be examined closely [11,12]. There are two fundamental assumptions in SHPB compression tests, one of which is that the stress and deformation of the specimen is uniform and the other is that there is one-dimensional elastic stress wave propagation in the bars. The accuracy of the experimental results can be ensured by minimizing the influence of factors violating these two assumptions. Various factors have been determined to result in inaccurate SHPB experimental results, among which the inertia effect and interfacial friction are the most critical ones. The axial inertia results in a stress difference between the two ends of the specimen, which violates the first assumption of SHPB experiments. The radial inertia of the specimen during compression leads to multiaxial stress states in the specimen. The stress calculating method based on the second assumption will exceed the true stress of materials. In addition, the interfacial friction generates a radial-directional shear stress on the specimen end surface, which not only changes the stress state in the specimen from uni-axial to multi-axial but also affects the original shape of the specimen cross-section. When the stress waves cross variable cross-section areas, nonuniform deformation is promoted and the stress measurement is affected. Thus, evaluating the influence of inertia effects and interfacial friction in SHPB tests is vital for the accurate determination of the dynamic mechanical properties.
A theoretical model was firstly established by Davies and Hunter [13]. They found there was additional stress owing to the radial inertia effects in SHPB tests, and derived the actual stress calculation in a cylindrical specimen using the law of energy conservation during impact testing. A slenderness ratio was suggested for a specimen design to eliminate the difference. Based on Davies' analysis, Samanta [14] took the change rate of the material into consideration, as there will be an acceleration in the specimen during SHPB experiments even if the loading strain rate is constant. Gorham [15] released the axial freedom of specimen ends to the effect of axial and radial inertia effect in SHPB experiments. They also wrote the expression of specimen stress in terms of the mean stress at the two specimen-bar interfaces rather than just the output bar stress, which was more consistent with the experimental results. Sen et al. [16] modified Gorham's analysis to non-cylindrical specimens with a more complicated expression in strain acceleration terms. The results for non-cylindrical specimens are consistent with those of cylindrical specimens [24].
Various experimental analyses have been conducted by researchers. Jia [17] suggested that the size of the copper specimen should be reduced with the increasing strain rate, which requires a smaller and more precise SHPB experimental system. Casem [18] and Zhang [19] conducted SHPB experiments on polymeric materials with a low density and strength, and concluded that inertial effects have a significant effect on the stress level of polymeric materials in the initial deformation stage, resulting in the fluctuation and peak value of the output signal. Song [20] suggested an annular specimen design to reduce the effect of inertia and proposed an analytical solution for inertia calculation for such a specimen design. Bischoff [21] and Zhang [22] concluded that the inertial confinement is not significant in experiments with a strain rate up to 200 s −1 . However, it is generally (but not universally) accepted that inertial confinement plays an important role in the SHPB results when the strain rate is higher than 1000 s −1 ; According to experimental observation, Bertholf and Karnes [25] established a two-dimensional numerical model to analyze the effect of inertial and interfacial effects on the material dynamic properties. They concluded that a suitable specimen design and enough lubrication can minimize these factors' effects. Meng and Li [26] introduced two coefficients to measure the effects of friction and specimen size on the stress uniformity in the specimen when simulating the SHPB experiments by FEM to examine the effects of radial inertia and end friction. Zencker [27] conducted simulations which showed that the optimum specimen slenderness ratio proposed by Davies and Hunter [13] was not accurate as it did not consider the friction effects. Iwamoto and Yokoyama [28] conducted computational simulations to demonstrate the effects of inertia on SHPB measurements considering both rate-dependent and rate-independent material models to demonstrate the inertia effect in SHPB tests.
In the present study, the influence of friction effects and inertia on the dynamic mechanical properties is numerically analyzed, particularly for different materials and loading conditions. An analytical model is established to quantitatively analyze the interfacial friction and inertia effects in SHPB tests. Methods to minimize the influence of interfacial friction and inertia are proposed based on the simulation results and the analytical model.
The structure of the paper is as follows. In Section 2, the specimen deformation characteristics in SHPB experiments are theoretically analyzed. In Section 3, a series of finite element models of SHPB tests are established to study the friction effects and inertia effects. Friction coefficients between 0 and 0.50, the specimen shape (cubic and cylinder), and the specimen slenderness (0.5, 1, 15) are considered. In Section 4, an energy conservation analysis based on the theory by Davies and Hunter (1963) is carried out, and a quantitative equation for the relation of the stress measured in the SHPB to the interfacial friction and inertia effects including the basic material parameters and the experimental conditions is derived. In Section 5, the results of the numerical simulation and analytical model are discussed, based on which methods to minimize the influence of interfacial friction and inertia are proposed.

Principle of Split Hopkinson Pressure Bar
A conventional SHPB experimental device consists of a striker, an incident bar, and a transmit bar, as shown in Figure 1. The specimen is sandwiched between the incident bar and the transmit bar. At the beginning of experiments, the striker impacts the incident bar with an initial velocity V 0 . An incident stress pulse is generated and propagated into the incident bar. When the incident stress pulse reachs the interface between the specimen and bars, a reflected pulse σ r is generated in the incident bar and a transmitted pulse σ t is generated in the transmitted bar due to the stiffness mismatch between the bars and specimen. These three pulses will be recorded by the strain gauge on the incident bar and transmit bar. Based on the assumptions of SHPB, the stress , strain , and strain rate of the specimen can be calculated as follows: The elastic stress wave propagation speed C, incident pulse σ i , and the duration of the created incident stress pulse ∆t can be defined as follows: where E is Young's modulus of the bar, ρ is the mass density of the bar, V 0 is the velocity of the striker, and L is the length of the striker. Based on the assumptions of SHPB, the stress σ s , strain ε i , and strain rate . ε s of the specimen can be calculated as follows: (4)

Interfacial Effect in SHPB Tests
In a SHPB experiments, for the sufficiently lubricated contact condition the friction coefficient on the interface is zero. The specimen is uniformly shortening in the axial direction (Z direction) and freely expanding in the radial direction (γ direction), as is shown in Figure 2   Based on the assumptions of SHPB, the stress , strain , and strain rate of the specimen can be calculated as follows:

Interfacial Effect in SHPB Tests
In a SHPB experiments, for the sufficiently lubricated contact condition the friction coefficient on the interface is zero. The specimen is uniformly shortening in the axial direction (Z direction) and freely expanding in the radial direction ( direction), as is shown in Figure 2 [30].
For the frictional contact condition, the radial expansion of the specimen ends is a constraint due to the friction forces caused by contact pressure. Consequently, the cross-section area of the specimen varies in the axial direction. For an infinitesimal element along the axial of a deformed specimen, a stress pulse propagates from one side to another side, in which the area varies. There will be a stress difference between the two ends due to the stress equivalence requirement in the dynamic tests. The difference is determined by the interfacial coefficient and specimen diameter [28,30,31].

Finite Element Model of SHPB
A series of nonlinear finite element models using the commercial software ABAQUS/Explicit are established to simulate the effects of interfacial friction and inertia on the dynamic response for different materials and loading conditions. For the frictional contact condition, the radial expansion of the specimen ends is a constraint due to the friction forces caused by contact pressure. Consequently, the cross-section area of the specimen varies in the axial direction. For an infinitesimal element along the axial of a deformed specimen, a stress pulse σ i propagates from one side to another side, in which the area varies. There will be a stress difference between the two ends due to the stress equivalence requirement in the dynamic tests. The difference is determined by the interfacial coefficient and specimen diameter [28,30,31].

Finite Element Model of SHPB
A series of nonlinear finite element models using the commercial software ABAQUS/Explicit are established to simulate the effects of interfacial friction and inertia on the dynamic response for different materials and loading conditions.
Since the geometries of bars and specimens are axisymmetric, a quarter of the model was used throughout the simulation work. Both the incident bar and transmit bar are 2000 mm long and modelled with 494,400 C3D8R elements. The striker is 200 mm long with 30,800 C3D8R elements. The dynamic load is applied by imposing the initial velocity on the striker. Two initial impact velocities are considered: 3 and 10 m/s. The stress waves are recorded by the strain gauges as the position in the experiments. Surface-to-surface contact is set between the specimen and bars to prevent interpenetration. Two specimen shapes are considered: a cylinder specimen and a cubic specimen (Figure 3a,b). The friction coefficients are set between 0 and 0.5 with an increment of 0.1. Here, 0 represents ideal lubrication and 0.5 is a large friction coefficient. The range 0-0.5 can cover the possible range of friction coefficients in SHPB tests [30]. To investigate the effect of specimen slenderness on the dynamic response, three specimen slenderness are considered: 0.5, 1, and 1.5 (Figure 3c,e).

Interfacial Effect in SHPB Tests
In the SHPB experiments, when the incident stress pulse crosses a variable cross-sectional specimen, the confining effect due to friction will bring a non-uniform stress distribution in the specimen axial direction. The value of the transmitted stress pulse is influenced by the radius difference in the cross section area and the friction coefficient. As the assumption of uniform stress distribution is challenged, the results obtained by the traditional SHPB processing method become unreliable. The analysis model of the inertia effect of the specimen is shown in Figure 4. and are the diameter and the thickness of a specimen, respectively. S is the reference position in the axial coordinate; and are the velocities at specimen ends. is the deformation velocity of a specimen; P1 and P2 are the pressure forces on the interaction surface between the specimen and bars. and are the friction forces caused by P1 and P2. The striker and incident/transmitted bars are modelled by linear elastic materials with a Young's modulus of 150 GPa and Possion's ratio of 0.3. Specimens are modeled with rate-independent and rate-dependent elastic-plastic material and rate-dependent viscoelasticity material, respectively.
For rate-independent material, the constituting relationship is described by elastic-plastic materials.
where A and n are the material parameters shown in Table 1. For rate-dependent elastic-plastic material, the flow stress is described by the Johnson-Cook constitutive model using the experimental results in Refs. [27,28,[30][31][32], expressed as follows: where T is the temperature parameter; T room is the room temperature and also the reference temperature; T melt is the melt temperature; A, B, C, n, and m are the material parameters and are shown in Table 2. For rate-dependent viscoelasticity material, the flow stress is described in the Zhu-Wang-Tang (ZWT) constitute model whose parameters are obtained from experimental testing [33] and expressed as follows: where E 0 , E 1 , E 2 , α, β, θ 1 , and θ 2 are the material constants listed in Table 3. This model is implemented with the explicit dynamic finite element software ABAQUS with the user subroutine to define material behavior (VUMAT). Table 3. Parameters for the rate-dependent viscoelasticity material.

Interfacial Effect in SHPB Tests
In the SHPB experiments, when the incident stress pulse crosses a variable cross-sectional specimen, the confining effect due to friction will bring a non-uniform stress distribution in the specimen axial direction. The value of the transmitted stress pulse is influenced by the radius difference in the cross section area and the friction coefficient. As the assumption of uniform stress distribution is challenged, the results obtained by the traditional SHPB processing method become unreliable. The analysis model of the inertia effect of the specimen is shown in Figure 4.

Interfacial Effect in SHPB Tests
In the SHPB experiments, when the incident stress pulse crosses a variable cross-sectional specimen, the confining effect due to friction will bring a non-uniform stress distribution in the specimen axial direction. The value of the transmitted stress pulse is influenced by the radius difference in the cross section area and the friction coefficient. As the assumption of uniform stress distribution is challenged, the results obtained by the traditional SHPB processing method become unreliable. The analysis model of the inertia effect of the specimen is shown in Figure 4. and are the diameter and the thickness of a specimen, respectively. S is the reference position in the axial coordinate; and are the velocities at specimen ends. is the deformation velocity of a specimen; P1 and P2 are the pressure forces on the interaction surface between the specimen and bars. and are the friction forces caused by P1 and P2.
The deformation velocity of sample can be expressed as: Assuming that the strain rate in the cross-section area of the specimen is uniform, the velocity field of the specimen in the axial direction can be expressed as: The deformation velocity v d of sample can be expressed as: Assuming that the strain rate in the cross-section area of the specimen is uniform, the velocity field of the specimen in the axial direction can be expressed as: Assuming that the volume of the specimen is constant during compression, the velocity field of the specimen in radial direction can be expressed as: Assuming that the interfacial friction has little influence on the velocity field, temporal change rate of the kinetic energy . E k of the specimen can be expressed as: The temporal change rate of the deformation energy . E p can be expressed as: In which E s is the Young's modulus of materials. The temporal rate of change of the external work can be expressed as: Assuming that the friction is uniformly distributed on the end surface, the temporal rate of change of the external work can be expressed as: where η is the friction coefficient. The equation of motion in the axial direction is: where v s is the rigid body motion velocity of the specimen. Thus, Assuming that the deformation process of the specimen is adiabatic, the conservation of energy can be expressed as: Substituting Equations (9)-(15) into (16), it can be concluded in Equation (17) if the average stress (0.5 * (P 1 + P 2 )) is used to measure the specimen stress in SHPB experiments.
It can be concluded that the strain acceleration as well as the geometrical properties of the specimen are the key factors causing the inertia effects in the SHPB test, so the resulting dynamic strength enhancement can be simply expressed using only the dominant terms, as shown in: The interfacial friction can be expressed as: The inertia effect can be expressed as: .

Influence of Interfacial Friction on SHPB Dynamic Results
The dynamic response of the cylinder specimen with a slenderness of 0.66 under SHPB compression for different materials is numerically analyzed by considering the friction coefficient set between 0 and 0.5 with an increment 0.1. The results of the stress pulse of the incident bar and transmit bar are presented in Figures 5-7).
For rate-independent materials under a 3 m/s impact velocity, the duration of the waves is essentially not affected by the friction coefficient and the transmit stress levels increase by about 13% when the friction coefficient increases from 0 to 0.5. As the impact velocity is 10 m/s, the duration of the waves is slightly affected and the transmit stress pulses increase by about 9%. For rate-dependent elastic plastic materials under a 3 m/s impact velocity, the duration of the waves and the transmit stress levels are affected slightly by the friction coefficient. As the impact velocity increases to 10 m/s, the duration of the waves is still not affected but the transmit stress levels increase by about 10%, with the friction coefficient increasing from 0 to 0.5. For rate-dependent viscoelasticity materials under a 3 m/s impact velocity, the duration of the waves slightly decreases with the friction coefficient increases and the transmit stress levels increase by about 19% when the friction coefficient increases from 0 to 0.5. The growth value of the transmit stress level is about 16% when the impact velocity is 10 m/s. For both the 3 and 10 m/s cases, there is a cross when the viscosity deformation process, after which point the reflected stress level for the case with a smaller friction coefficient was found to be with a higher transmit stress level than the case with a larger friction coefficient.

Influence of Specimen Shape on SHPB Dynamic Results
The dynamic response of specimens with different geometry shapes under SHPB compression for rate-independent material, rate-dependent elastic-plastic material, and rate-dependent viscoelasticity material are numerically analyzed by considering the cubic specimen and cylinder specimen (both specimens with 0.66 slenderness). Two initial impact velocities were considered, 3 and 10 m/s. The contact condition is set frictionless. The results of the stress pulses of the incident and transmitted bars are presented in Figures 8-11.
For rate-independent materials, the duration of the waves is not affected significantly by the specimen geometry for both the 3 and 10 m/s impact velocities. When the impact velocity is 3 m/s, the transmit stress level of the cubic specimen is 30% higher than that of the cylinder specimen, and the value decreases to 16% when the impact velocity is 10 m/s. For rate-dependent elastic plastic materials, the duration of the waves is still not affected by the specimen shape for the 3 and 10 m/s impact velocities. When the impact velocity is 3 m/s, the transmit stress level of the cubic specimen is 22.5% higher than that of the cylinder specimen, and the value decreases to 11.7% when the impact velocity is 10 m/s. For rate-dependent viscoelasticity materials, similarly the duration of the waves is not affected by the specimen shape for both impact velocities. When the impact velocity is 3 m/s, the transmit stress level of the cubic specimen is 27% higher than that of the cylinder specimen, and the value decreases to 16% when the impact velocity is 10 m/s. Besides this, the transmit wave shape of the cylinder specimen is a triangle, but it transforms into square wave when the specimen shape is cubic. The interfacial friction brings extra stress and strain rate variation in the SHPB specimen, whose value is related to the testing material properties and impact velocities. Rate-dependent viscoelasticity material is the most sensitive to friction effects among the three materials. However, the increase rate of the transmitted stress value is not influenced by the impact velocity; rate-dependent elastic-plastic material is sensitive to the interfacial friction. The increase rate is significantly influenced by the impact velocity; rate-independent material is influenced by interfacial friction, and its effect is independent of the impact velocity. The conventional SHPB data computed method without consideration of the effect of interfacial friction leads to a falsely higher flow stress for the testing material regardless of the analysis method (three-wave method or two-wave method).

Influence of Specimen Shape on SHPB Dynamic Results
The dynamic response of specimens with different geometry shapes under SHPB compression for rate-independent material, rate-dependent elastic-plastic material, and rate-dependent viscoelasticity material are numerically analyzed by considering the cubic specimen and cylinder specimen (both specimens with 0.66 slenderness). Two initial impact velocities were considered, 3 and 10 m/s. The contact condition is set frictionless. The results of the stress pulses of the incident and transmitted bars are presented in Figures 8-11. for testing materials. Rate-dependent viscoelasticity material is the most sensitive to specimen shape among the three materials. However, the increase rate of the transmitted stress value decreases with the increase in the impact velocity. Rate-dependent elastic-plastic material is sensitive to the specimen shape, whose transmitted stress level increases significantly with the impact velocity increasing; rateindependent material is influenced by interfacial friction, and its effect decreases with the impact velocity.   among the three materials. However, the increase rate of the transmitted stress value decreases with the increase in the impact velocity. Rate-dependent elastic-plastic material is sensitive to the specimen shape, whose transmitted stress level increases significantly with the impact velocity increasing; rateindependent material is influenced by interfacial friction, and its effect decreases with the impact velocity.

Influence of Specimen Configuration on SHPB Dynamic Results
The dynamic response of the cylinder specimen with different slenderness under SHPB compression for rate-independent material, rate-dependent elastic-plastic material, and ratedependent viscoelasticity material are numerically analyzed by considering 0.5, 1, 1.5. The contact condition is set frictionless. The results of the stress pulses of the incident and transmitted bars are presented in Figures 12-14.
For rate-independent materials, the duration of the waves is not affected by the specimen slenderness regardless of the impact velocities. When the impact velocity is 3 m/s, the transmit stress level of the specimen first increases then decreases with the specimen slenderness increasing. The maximum stress value is 31% higher than the minimum stress value. When the impact velocity is 10 m/s, the transmit stress level shows a similar trend and the ratio of increment is 44%. For ratedependent elastic plastic materials, the duration of the waves and the transmit stress level are not affected by the specimen slenderness for different impact velocities. For rate-dependent viscoelasticity materials, similar observations as for rate-independent materials are noticed. The duration of the waves is not affected by the specimen slenderness. The transmit stress level of the specimen first increases then decreases with the specimen slenderness increasing. The ratio of increment is 18% for the 3 m/s impact velocity and 17% for the 10 m/s impact velocity.
The specimen geometry has a significant effect on the reflected and transmitted wave stress levels for rate-independent materials and rate-dependent viscoelasticity materials, but not for ratedependent elastic plastic materials. The transmit stress value first increases then decreases with the increasing specimen slenderness. For rate-independent materials, the duration of the waves is not affected significantly by the specimen geometry for both the 3 and 10 m/s impact velocities. When the impact velocity is 3 m/s, the transmit stress level of the cubic specimen is 30% higher than that of the cylinder specimen, and the value decreases to 16% when the impact velocity is 10 m/s. For rate-dependent elastic plastic materials, the duration of the waves is still not affected by the specimen shape for the 3 and 10 m/s impact velocities. When the impact velocity is 3 m/s, the transmit stress level of the cubic specimen is 22.5% higher than that of the cylinder specimen, and the value decreases to 11.7% when the impact velocity is 10 m/s. For rate-dependent viscoelasticity materials, similarly the duration of the waves is not affected by the specimen shape for both impact velocities. When the impact velocity is 3 m/s, the transmit stress level of the cubic specimen is 27% higher than that of the cylinder specimen, and the value decreases to 16% when the impact velocity is 10 m/s. Besides this, the transmit wave shape of the cylinder specimen is a triangle, but it transforms into square wave when the specimen shape is cubic. Figure 11 shows the stress distribution in the specimen with different materials and impact velocities. For all cases, the stress concentration was located at the corner of the specimen. The rate-dependent viscoelasticity materials show the most severe stress concentration when the impact velocity is 10 m/s.
The specimen geometry has a significant effect on the reflected and transmitted wave stress levels, but not the stress wave duration. Cubic specimens present a higher stress than cylinder specimens regardless of the testing material properties. This is due to the fact that the stress waves would converge at the endpoints of the cuboid specimen and result in a mistaken higher flow stress for testing materials. Rate-dependent viscoelasticity material is the most sensitive to specimen shape among the three materials. However, the increase rate of the transmitted stress value decreases with the increase in the impact velocity. Rate-dependent elastic-plastic material is sensitive to the specimen shape, whose transmitted stress level increases significantly with the impact velocity increasing; rate-independent material is influenced by interfacial friction, and its effect decreases with the impact velocity.

Influence of Specimen Configuration on SHPB Dynamic Results
The dynamic response of the cylinder specimen with different slenderness under SHPB compression for rate-independent material, rate-dependent elastic-plastic material, and rate-dependent viscoelasticity material are numerically analyzed by considering 0.5, 1, 1.5. The contact condition is set frictionless. The results of the stress pulses of the incident and transmitted bars are presented in Figures 12-14.

Methods to Minimize the Effect of Inertia Effects and End Friction
Thus, ways to minimize the effect of inertia effects and end friction can be considered as follows: In order to eliminate the effect of term, the coefficient should be approximately zero. Thus, the optimum specimen slenderness ratio should be 1/3 during experiments, which is also the criterion derived by Davies and Hunter in 1963. However, the specimen slenderness ratio increases with the specimen deformation increasing. The criterion derived by Davies and Hunter is only suitable at the experiment's beginning when the deformation is small. Since it is impossible to ensure that the (a)

Methods to Minimize the Effect of Inertia Effects and End Friction
Thus, ways to minimize the effect of inertia effects and end friction can be considered as follows: In order to eliminate the effect of term, the coefficient should be approximately zero. Thus, the optimum specimen slenderness ratio should be 1/3 during experiments, which is also the criterion derived by Davies and Hunter in 1963. However, the specimen slenderness ratio increases with the specimen deformation increasing. The criterion derived by Davies and Hunter is only suitable at the experiment's beginning when the deformation is small. Since it is impossible to ensure that the (a)

Methods to Minimize the Effect of Inertia Effects and End Friction
Thus, ways to minimize the effect of inertia effects and end friction can be considered as follows: In order to eliminate the effect of term, the coefficient should be approximately zero. Thus, the optimum specimen slenderness ratio should be 1/3 during experiments, which is also the criterion derived by Davies and Hunter in 1963. However, the specimen slenderness ratio increases with the specimen deformation increasing. The criterion derived by Davies and Hunter is only suitable at the experiment's beginning when the deformation is small. Since it is impossible to ensure that the (a) For rate-independent materials, the duration of the waves is not affected by the specimen slenderness regardless of the impact velocities. When the impact velocity is 3 m/s, the transmit stress level of the specimen first increases then decreases with the specimen slenderness increasing. The maximum stress value is 31% higher than the minimum stress value. When the impact velocity is 10 m/s, the transmit stress level shows a similar trend and the ratio of increment is 44%. For rate-dependent elastic plastic materials, the duration of the waves and the transmit stress level are not affected by the specimen slenderness for different impact velocities. For rate-dependent viscoelasticity materials, similar observations as for rate-independent materials are noticed. The duration of the waves is not affected by the specimen slenderness. The transmit stress level of the specimen first increases then decreases with the specimen slenderness increasing. The ratio of increment is 18% for the 3 m/s impact velocity and 17% for the 10 m/s impact velocity.
The specimen geometry has a significant effect on the reflected and transmitted wave stress levels for rate-independent materials and rate-dependent viscoelasticity materials, but not for rate-dependent elastic plastic materials. The transmit stress value first increases then decreases with the increasing specimen slenderness.

Methods to Minimize the Effect of Inertia Effects and End Friction
Thus, ways to minimize the effect of inertia effects and end friction can be considered as follows: In order to eliminate the effect of term, the coefficient should be approximately zero. Thus, the optimum specimen slenderness ratio should be 1/3 during experiments, which is also the criterion derived by Davies and Hunter in 1963. However, the specimen slenderness ratio increases with the specimen deformation increasing. The criterion derived by Davies and Hunter is only suitable at the experiment's beginning when the deformation is small. Since it is impossible to ensure that the slenderness ratio stays constant during the experiments, a new criterion which minimizes the absolute sum of the coefficient can be applied to obtain the optimum specimen slenderness ratio. According to the constant strain rate hypothesis and the plastic incompressibility hypothesis, the specimen thickness and diameter can be expressed by its final strain as: Thus, the optimum specimen slenderness ratio to minimize the inertia effects can be expressed as: The optimum specimen slenderness ratio is related to its final strain, whose relationship can be fitted by Equation (25), and expressed in Figure 15: The effect of interfacial friction can be ignored if its influence is less than 5%. Thus, the friction coefficient can be calculated as: The above analysis is an efficient method to assess the influence of interfacial friction. A more general method to eliminate the interfacial friction effect is to correct the experimental results according to the theoretical analysis.
The optimum specimen slenderness ratio is related to its final strain, whose relationship can be fitted by Equation (25), and expressed in Figure 15: The effect of interfacial friction can be ignored if its influence is less than 5%. Thus, the friction coefficient can be calculated as: The above analysis is an efficient method to assess the influence of interfacial friction. A more general method to eliminate the interfacial friction effect is to correct the experimental results according to the theoretical analysis.

Conclusions
Distinguished with the available literature, friction effects and inertia effects were analyzed theoretically and numerically for different materials in this study. Typical materials were represented by rate-independent material, rate-dependent elastic-plastic material, and rate-dependent viscoelastic material. Inertia effects were investigated by two factors: specimen shape and specimen slenderness. Two impact velocities of 3 and 10m/s were considered. An analytical model was established to quantitatively analyze the interfacial friction and inertia effects in the SHPB tests. Methods to minimize the influence of interfacial friction and inertia were proposed based on the simulation results and analytical model. This study provided material selection suggestions for army Figure 15. Relation between the optimum specimen slenderness ratio and the final strain of the specimen.

Conclusions
Distinguished with the available literature, friction effects and inertia effects were analyzed theoretically and numerically for different materials in this study. Typical materials were represented by rate-independent material, rate-dependent elastic-plastic material, and rate-dependent visco-elastic material. Inertia effects were investigated by two factors: specimen shape and specimen slenderness. Two impact velocities of 3 and 10m/s were considered. An analytical model was established to quantitatively analyze the interfacial friction and inertia effects in the SHPB tests. Methods to minimize the influence of interfacial friction and inertia were proposed based on the simulation results and analytical model. This study provided material selection suggestions for army and other protective robot designers when they analyze the dynamic properties of materials for a safety design. The main conclusions are as follows: The interfacial friction brings extra stress for the SHPB specimen, whose values is related to the testing material properties and impact velocities. Rate-dependent viscoelasticity material is the most sensitive to friction effects among the three materials. However, the increase rate of the transmitted stress value is not influenced by the impact velocity; there is a cross in the viscosity deformation process, after which point the reflected stress level for the case with a smaller friction coefficient was found to be with a higher transmit stress level than the case with a larger friction coefficient. The rate-dependent elastic-plastic material is sensitive to the interfacial friction. The increase rate is significantly influenced by the impact velocity; rate-independent material is influenced by interfacial friction, and its effect is independent of the impact velocity.
The specimen geometry has a significant effect on the reflected and transmitted wave stress levels, but not the stress wave duration. Cubic specimens present a higher stress than cylinder specimens regardless of the testing material properties. Rate-dependent viscoelasticity material is the most sensitive to specimen shape among the three materials. However, the increase rate of the transmitted stress value decreases with the impact velocity increasing; rate-dependent elastic-plastic material is sensitive to the specimen shape, whose transmitted stress level increases significantly with the impact velocity increasing. Rate-independent material is influenced by interfacial friction, and its effect decreases with the impact velocity.
The result of impact response of different materials with different interfacial and geometry can be used in the structural design of these materials, especially under impact loads. Under the same loading condition, the difference in the mechanical response of different materials can be attributed to the difference in microstructure. This result can provide some theoretical basis for researchers to improve the material properties.
The specimen geometry has a significant effect on the reflected and transmitted wave stress levels for rate-independent materials and rate-dependent viscoelasticity materials, but not for rate-dependent elastic plastic materials. The transmit stress value first increases then decreases with the increasing specimen slenderness.
The influence of the friction and inertia effect on the SHPB test results is related to the specimen size, deformation strain, strain rate, density of the specimen, and friction coefficient of the interface. Methods to reduce the influence of the inertia effect on the experimental results can be achieved by reducing the influence of these factors, such as using constant strain rate loading experiment technology and a lubricated contact interface. The optimum specimen slenderness ratio is related to its final strain, whose relationship can be described by a linear function. The interface friction is a very complex dynamic process; when the object is deformed in the friction experiments, the process will be even more complex. However, this needs to be researched, since it always occurs in daily life.