1. Introduction
Low impedance materials such as elastomers are finding new applications in industries such as automotive, biomedical, aerospace, etc. A typical dynamic application of elastomers is shock absorbers for automotive and aerospace vehicles. Since hyperelastic materials like elastomers show dependency on the strain rate at which the loading is applied, researchers and engineers need to rely on properties determined under static and quasi-static conditions for such applications. Therefore, accurate models and techniques are needed to characterize elastomers in the dynamic range. In order to generate strain rate dependent constitutive models for these materials, testing at different strain rates is required.
The most common method to characterize materials at a strain rates greater than 100/s–10,000/s is the Kolsky bar technique. This apparatus has been widely used to dynamically characterize materials such as metals [
1], ceramics [
2], foam [
3], composites [
4], and smart materials [
5]. The Kolsky bar is an ideal choice to characterize the stress strain response of material as a function of strain rate. A family of curves can be generated for dynamic models to accurately predict performance in impact applications. This can be done by fitting the response curves to a strain rate–dependent constitutive equation in order to derive the material parameters.
In compression, this apparatus is widely known as the Split Hopkinson pressure bar (SHPB), named after John Hopkinson (1849–1898) and his son, Bertram Hopkinson (1849–1898), known for demonstrating the transmission of stress waves in iron wires [
6]. They conducted plate impact experiments to measure the pressure-time curves by a detonation or a bullet impact. Davies (1948) [
7] studied this technique critically and discussed the dispersion of stress waves through long rods. In 1949, Kolsky modified the SHPB by adding an elastic bar on both sides (incident or input bar and reflected or output bar) and measured stress strain responses of materials under impact-loading conditions [
8]. Kolsky also analyzed the importance of radial inertia in the specimen during deformation, which is susceptible when the strain rate is changed very rapidly [
9]. The transmission of stress waves was recorded using a condenser microphone. Later, a strain-gauging technique was introduced for this purpose [
10]. Krafft et al (1954) modified the Kolsky bar design to generate the loading through an accelerated projectile launch [
10]. Since then, this technique has been modified and extended to incorporate tension, torsion and bend-testing of materials at high strain rates. The technique has been reviewed by Follansbee (1995) [
11], Nemat-Nasser (1991) [
12], Gray et al (2000) [
13], and Field et al (2004) [
14].
The finite element modelling of Kolsky bar has been extensively reported in previous studies, e.g., for the determination of incident pulse shaping design [
15,
16], modelling of pulse shaper [
17], for numerical verification of the apparatus design [
18], misalignment effect in the loading bars [
19], and investigation on wave dispersion in bars [
20]. This study aims to provide comprehensive methodology to model elastomers in a Kolsky bar experiment with a focus on specimen geometry optimization, experimental verification, and pulse shaper geometry effect.
A conventional Kolsky bar apparatus needs to be modified to test soft materials. Since soft materials are characterized by low stiffness and low elastic wave velocities, considerations such as impedance mismatch and uniform response of the specimens have to be taken into account; otherwise, it could result in a distorted and weak signal. To overcome these effects, essential modifications are essential, such as using low impedance bar materials, hollow rods [
21], and pulse shapers [
22].
In this paper, we examine several techniques as well as other important design parameters and requirements. A design for dynamic characterization of elastomers is presented based on the modifications needed for its use with soft materials. Investigations are made on the use of hollow transmission bar, pulse shaping technique and homogenous deformation of the specimen. Also, in order to achieve a constant strain-rate the loading pulse needs to be modified using a pulse shaper. Numerical simulations are performed using commercial explicit dynamic package (Abaqus Explicit). Various design parameters such as specimen geometry and pulse shaper have also been optimized using numerical simulations and experimental results. Finally, the optimized pulse shaper and geometry design are incorporated into the Kolsky bar for the experimental verification. A family of curves for the dynamic stress–strain responses of nitrile butadiene rubber (NBR) specimen has been generated through the modified Kolsky bar technique based on the optimized design of specimen and pulse shaper.
2. Kolsky Bar Design and Setup
2.1. Conventional Kolsky Bar
A Kolsky bar experiment is the most commonly high strain rate testing equipment used to determine the dynamic response of various engineering materials. Although there is no standard design of this apparatus, most researchers use an apparatus that includes the following components (as shown in
Figure 1):
two long cylindrical elastic bars (incident, transmission);
a striker launching mechanism (usually achieved by a gas gun);
force/strain measuring sensors;
bearing and base setup, for the axial alignment of bars; and
data acquisition system.
In a Kolsky bar, a sample is sandwiched between the incident and transmission bar. When the striker bar impacts the incident bar, an elastic compressive stress wave is generated. This stress wave travels through the incident bar. Upon reaching the incident bar and specimen interface, part of the wave is reflected and some part of it is transmitted through the specimen into the transmission bar (
Figure 2). The pulse profile for the incident, transmitted, and reflected waves are recorded as a function of time. This pulse profile is further analyzed to determine the dynamic stress strain response of the specimen [
23]. A Kolsky bar setup is customized in different ways to allow for the testing of a wide variety of materials.
The traditional selection of the bar material is usually a high-strength structural material such as steel, aluminum, or nickel alloy. These high strength materials are used because the yield strength of the bar determines the maximum attainable stress. Some researchers have also used a selection of low strength materials for bars such as magnesium (45 GPa) or polymeric bars (<20 GPa) [
24,
25,
26]. However, the disadvantage of using polymeric bars is their low heat conductivity and greater degree of dispersion and attenuation as compared to metallic bars.
In an ideal scenario, the striker bar is chosen to be the same diameter and material as of the incident and transmission bar [
8]. The impacting velocity and length of the striker dictates the loading time of the pulse. This loading time is critical when testing low impedance materials. It has to be greater than the time taken for the specimen to reach stress equilibrium [
8]. Under equilibrium, the amplitude of the reflected wave is directly proportional to maximum achievable strain rate. Similarly, the amplitude of the transmitted pulse is also directly proportional to the stress developed within the specimen [
8,
27].
2.2. Theory of the Kolsky Bar
The working principle of the Kolsky bar is based on the one-dimensional propagation of the compressive wave through the bars. This one-dimensional compressive elastic wave is related to the velocities that the particles experience at either interface of the bar and specimen by the following equations:
where
ɛi,
ɛr,
ɛt are the measured incident, reflected, and transmitted strain pulses recorded from the strain gages.
Cb is the elastic wave speed of the bar material.
Using the above equations, the relation for averaged engineering strain integrated over the pulse duration
t in specimen having the length
Ls is
and the stresses at both of the interfaces of the specimen can be determined as follow:
where
Ab and
As are the bar and specimen cross section areas and
Eb is the Young’s Modulus of both bars.
In order to determine the stress-strain response, we assume that the equilibrated state of stress is in the specimen (
σ1 =
σ2). Making this assumption we arrive at the following relation:
Thus Equations (3) and (5) can be solved to determine the stress, strain rate, and strain profile in the specimen.
2.3. Kolsky Bar Design Considerations
The performance specification of the Kolsky bar apparatus is centered on the choice of specimen size and the max desired strain rate. The strain in the incident pulse can be determined by the impact velocity of the striker by the following relation:
where
V denotes the striker impact velocity.
A particular striking velocity is chosen to achieve the desired strain rate. The maximum amount of stress level and strain rate that can be achieved in the specimen, with a striking velocity of
V is given by the following relations:
Higher levels of stress can be achieved by increasing the impact velocity or by increasing the area mismatch Ab/As. Similarly, reducing specimen length and increasing the impact velocity would increase the maximum attainable strain rate.
When selecting the bar and specimen dimensions, certain design requirements have to be met for a valid Kolsky bar.
Lb/
Db ratio of at least 20 is required to allow the generation of one-dimensional wave propagation. Typically, this ratio is selected to be in order of 100 [
28]. The diameter of the bar should be around 2–4 times the diameter of the specimen. Selecting a large ratio for
Db/
Ds is desirable for achieving high stress but can result in punch problems. This ratio should be kept such that the bar surface does not act like a plane for the specimen.
Ls/
Ds is usually selected in between 0.6 and 1 [
28] (the
Ls/
Ds ratio also depends on the type of material being tested). Soft materials have shown dependency on the length of the specimen, which could result in axial and radial inertia. The striker bar is made up of the same material and diameter as of the incident and transmitter bars to avoid wave dispersion.
2.3.1. Uniform Deformation and Inertial Effects
When designing a Kolsky bar for soft materials, as discussed above, the measurements are based on the assumption of uniform deformation in the specimen. Due to the low strength, stiffness and acoustic impedance of soft materials, their mechanical response depends upon the compressive stress wave profile. Due to the low wave speed in specimen as compared to that of the bar, it takes more time for the specimen to achieve stress equilibration. The specimen can be considered to be in a state of stress equilibrium by allowing several reverberations of the stress wave inside the specimen. It requires an elastic wave to ring up almost three times inside the specimen in order to achieve near ideal stress distribution [
29]. The time for a one-way trip of stress wave through the specimen can be calculated as follows:
where
Ls is the length of the specimen and
Cs is the characteristic elastic wave velocity.
It can be deduced from Equation (13) that the length of the specimen is the most important variable in achieving a uniform state of stress in the specimen. Reducing the length of the specimen would help to achieve this state quickly, but on the contrary, it has been observed that a small value of Ls could result in the increased frictional effect at the specimen bar interfaces, creating distortion in the recorded signal. The non-uniform deformation also results in the axial and radial inertia which is evident when testing very soft materials like elastomeric foams and gels. This induced inertia is dependent on both the loading conditions and on the specimen geometry and density. These effects cannot be ignored when optimizing the specimen dimensions.
In addition, high strain rates can also produce inertia induced stress, especially when the specimen experiences large deformation. This inertia induced stress can greatly affect the strain profiles recorded by the strain gauges [
30].
2.3.2. Low Transmitted Signals
Another challenge associated with testing softer materials using the Kolsky bar experiment is that the amplitude of the transmitted pulse is very low [
31]. Because of the impedance mismatch between the specimen and bar materials, almost the entire the stress wave is reflected back, and only a very small portion of it is transmitted. This weak signal becomes difficult to measure by regular strain sensor and more responsive and sensitive strain gages are required. One way to address this problem is to use low impedance bars, e.g., polymer bars. Using polymer bars allows for greater wave dispersion, and spectral analysis techniques are therefore required to filter out the dispersion effect. Another approach is to use a hollow transmitter bar to help increase the amplitude of the transmitted signal up to several times [
32].
2.3.3. Loading Profile
Elastomers show dependency on the applied strain rate, generating different stress–strain responses for different loading rates. As a result, special care must be taken to ensure that the stress pulse profile is achieved at a constant strain rate with low rate of initial loading. A low initial rate of loading is required since acceleration induced inertia is inevitable in the initial period needed to reach a particular strain rate. A common way of avoiding this acceleration induced inertia is to use a constant rate of loading pulse with extended rise time.
2.3.4. Specimen Size Effect
Specimen size needs to be properly optimized to mitigate the effects of associated radial and axial inertial responses, achieving uniform deformation and attaining a desired strain level. The specimen diameter has to be smaller than the bar diameter. Upon the assumption that the maximum diameter during deformation does not exceed the bar diameter (for the choice of desirable strain), the required diameter of the specimen can be calculated as follows:
For example, a bar with a diameter of 20 mm and a specimen diameter of 8 mm can achieve a maximum engineering strain close to 0.8. In addition to the selection of the specimen diameter, appropriate length of the specimen needs to be determined. This also contributes toward visible inertial effects and the stress equilibration. A thinner specimen helps in achieving a state of equilibrated stress quicker [
33] but also produces visible frictional effects.
A specimen can be considered to be in equilibrated stress by comparing the stress profile on both of the faces of the specimen. Previous researchers have determined that equilibrium is reached if the ratio of the difference between the stress from either end “Δ
σ” and the mean stress in the specimen “
σm” is less than 0.05 [
34].
Song et al. [
33] demonstrated that a thin specimen maybe one of the requirements in Kolsky bar experiment but it is not sufficient enough to achieve equilibrium. Specimen design can also help in avoiding inertia in a radial direction. It was also demonstrated that this inertia reaches maximum at the center and diminishes at the edges. Therefore, using a hollow specimen could also help in addressing this challenge [
35].
2.3.5. Alignment of the Bars
A very important constraint in successfully conducting the Kolsky bar experiment is the alignment of the bars. It is known that misalignment could cause distortion in the recorded signal which in some cases cannot be differentiated from the actual strain signals. Because the goal is to record micro level strains, caution must be taken when setting up the apparatus. Kareem et al [
19] quantified the effect of misalignment on the recorded stress strain response of the specimen. It was experimentally verified that in order to minimize this distortion, the straightness of the bars has to be at least 0.08 mm per 205 m, centerless grinded to a target diameter with a tolerance of ±0.025 mm, and finally have a perfectly perpendicular lubricated end face with a tolerance of ±0.03°.
2.4. Experimental Setup
For the experimental design of the Kolsky bar apparatus, the specimen geometry was optimized to be Ds = 8 mm and Ls = 1.6 mm. The length and the diameter were selected to be, Lb = 2000 mm and Db = 20 mm respectively. These selected values satisfy the design parameter ratios that were discussed in the previous section. We get Lb/Db = 100, Db/Ds = 2.5, and Ls/Ds = 0.2. It can be noted that the Ls/Ds ratio is different when compared to the discussed range in the previous section (between 0.6 and 1). The specimen geometry was optimized carefully by taking into account the effect of Ls on the uniform deformation, and the equilibration of stress. Based on Equation (14) and the selection of Db and Ds, we can achieve a strain level up to 0.8. The striker diameter Dst is kept equal to the incident and the transmitter rods to avoid wave dispersions. The length of striker was selected to be Lst = 500 mm. Relatively longer Lst accounts for stress equilibration, as the loading pulse duration needs to be three times more than it takes the wave to travel back and forth inside the specimen. By using the relation , it can be seen that the choice of Lst would generate a loading pulse of duration 200 µs. Therefore, it is enough for the specimen to attain uniform stress distribution and respond according to it accordingly.
From the available options of bar material, Aluminum 6061 Anodized was selected. It is evident from Equation (7) that either Young’s modulus,
Eb or
Ab/
As must be reduced in order to increase the transmitted strain. This is the main reason for selecting Aluminum as bar material in lieu of other higher strength materials that have been reportedly used in Kolsky setups. It is recommended that the striker, transmitter, and incident bar are made from the same material. Because of the impedance mismatch and very weak transmitted wave, a hollow transmitter bar was chosen to increase the amplitude of the transmitted wave. When using hollow bars, modifications need to be made to the Equation (9). Accounting, for the change in the cross-sectional area of the bars with the use of a hollow transmission bar, the following relation can be derived [
32]:
To ensure uniform deformation, a pulse shaper was attached between the incident and striker bar impacting surface. Control and smoothness of the loading pulse can be achieved by using a pulse shaper as it absorbs the high frequency components of the pulse. Different choices of pulse shapers are used by researchers, such as polymer disks [
32], copper disks [
36] and tissue paper. In our study, pulse shaper geometry was optimized experimentally. Different dimensions of copper shaper were analyzed experimentally to determine their effect on the loading pulse profile and on homogenous deformation of elastomers.
As already discussed above, the maximum achievable strain rate depends upon the impact velocity of the striker. The most common way of accelerating the striker is by using a two staged compressed gas gun. The main supply of high-pressure gas or a compressor is used to fill the reservoir. In our testing configuration, we used a one-liter sampling cylinder with a pressure rating of up to 3000 psi. Pressure is released using a solenoid valve, with a pressure rating of 750 psi. As the velocity of the striker depends upon the initial gas pressure, choosing a high-pressure reservoir would generate a wide range of velocities in the striker. Having control over a wide range of velocities is also needed because the striking velocity relates to the strain rate achieved in the specimen. In our testing configuration, the striker could be accelerated to value up to 80 m/s upon exiting the barrel.
5. Conclusions
For the characterization of elastomers at strain rates of 100 s−1–10,000 s−1, a modified Kolsky bar was developed. In the design of such an apparatus, important considerations such as homogenous deformation of the test specimen, hollow transmission rod, and pulse shaping techniques were taken into account. A complete methodology to numerically model Kolsky bar in FEA is also described, which was used to optimize the specimen geometry and to study the effect of homogenous deformation considering various dimensions for the specimen. Experimental and numerical analysis was made on the pulse shaping technique to analyze the effect of thickness and diameter of the pulse shaper on the incident pulse. Moreover, resultant incident pulses generated from commonly used materials for pulse shapers were also taken into account and a design for a hybrid elastomer/copper pulse shaper is also discussed. Finally, employing the optimized geometry of the specimen and the pulse shaper, response of NBR rubber is investigated at various different strain rates.