# Experimental Investigation and Discrete Element Modelling of Composite Hollow Spheres Subjected to Dynamic Fracture

Authors: Arthur Coré; Jean-Benoît Kopp; Philippe Viot; Jean-Luc Charles; Fréderic Dau
Journal: International Journal of Polymer Science
 (2017)

---

## Abstract

This paper deals with the characterization and the numerical modelling of the collapse of composite hollow spherical structures developed to absorb energy during high velocity impacts. The structure is composed of hollow spheres (ϕ=2–30 mm) made of epoxy resin and mineral powder. First of all, quasi-static and dynamic (v=5 mm·min−1 to v=2 m·s−1) compression tests are conducted at room temperature on a single sphere to study energy dissipation mechanisms. Fracture of the material appears to be predominant. A numerical model based on the discrete element method is investigated to simulate the single sphere crushing. The stress-strain-time relationship of the material based on the Ree-Eyring law is numerically implemented. The DEM modelling takes naturally into account the dynamic fracture and the crack path computed is close to the one observed experimentally in uniaxial compression. Eventually, high velocity impacts (v>100 m·s−1) of a hollow sphere on a rigid surface are conducted with an air cannon. The numerical results are in good agreement with the experimental data and demonstrate the ability of the present model to correctly describe the mechanical behavior of brittle materials at high strain rate.


---

## Body

## 1. Introduction

Hollow sphere structure (HSS) belongs to cellular solids that have been studied recently for their multiple properties [1].Compared to foam structure, HSS has both closed and open porosity. Its applications can be of various types including energy absorber [2], acoustic damping [3], and thermal insulation [4]. In our case, HSS aims to absorb soft impacts energy on an airliner cockpit. Most of impacts are due to bird-strikes during take-off or landing at high velocity (vmax=175 m·s−1). Protective structures have to dissipate the kinetic energy of the projectile by all possible means: friction [5], plastic deformation [6], or fracture [7]. Many materials and structures have been tested to absorb energy of external impacts such as aluminium honeycombs [8], polymer [9], and metallic [10] foams and composite materials [7]. CELPACT [11] and MANSART [12] projects have recently focused on cellular structures designed for energy absorption. HSS is now investigated through the SAMBA (Shock Absorber Material for Bird-shield Application) project because of its promises in terms of specific energy dissipated (J·kg−1) during impact [13]. Hollow spheres can also easily fill a sandwich structure to be placed in front of the plane. The commonly used bird-shield structure measures approximately one square meter and ten centimeters thick which represents 10–12 kg.Hollow spheres used in this study are made of commercial epoxy resin, a thermosetting polymer subjected to crack propagation [14]. In fact, rapid crack propagation (RCP) is observed to be predominant at macroscopic scale to dissipate the available energy stored in the structure [15]. As it is commonly known, most thermosetting polymers have a strain rate and temperature dependence on their mechanical properties (Young’s modulus, yield strain) [16, 17]. An influence of the compression velocity can be predicted.In the interest of studying the fracture on polymer hollow spheres, a numerical model is investigated. It is based on the discrete element method, often used for rocks or soils modelling [18]. It reveals to be an interesting way to model the mechanical behavior of brittle materials. It allows large displacements and strain and the fracture of the material. It offers an alternative to the FEM by allowing a natural propagation of the crack [19]. Indeed the DEM does not require an adaptive meshing refinement at the crack tip contrary to some finite element methods (FEM) [20, 21]. The use of “classical” continuum models such as the FEM was well adopted to simulate the crushing of elastoplastic hollow spheres [22–25]. Despite several methods (elements erosion, cohesive zones, extended finite element method, etc.) developed to represent the fracture of a material [26–28], a discrete model is chosen. Continuous material behavior can now be simulated with the DEM by using cohesive beams linking each discrete element [29]. The crack closure can be taken into account by activating and computing contacts between discrete elements.Before optimizing an assembly of HSS (which corresponds to the macroscopic scale) as shown in Figure1(a) to dissipate more energy when crushed, the present paper focuses on the investigation of one sphere crushing (mesoscopic scale). The sphere behavior under impact, which remains important to study, is often quasi-invisible at macro-scale. It leads nevertheless to a better understanding of energy dissipation mechanisms. That is why, an experimental and numerical approach is proposed to take into account the fragmentation on a single composite hollow sphere. Experimental part will present the structure response as F(δ) (force displacement curve) as well as the characterization of the constitutive material for different strain rates and temperatures. A fracture criterion based on the principal stress and changing according to the Ree-Eyring formalism [30] is implemented in the DEM. Eventually, high velocity impacts are conducted with an air cannon. Hollow spheres are fired on a rigid surface, and multicracks are observed and compared to the numerical results.Figure 1
Hollow sphere structure (a), hollow sphere’s parameters (b), and an optical microscope capture of the constitutive material (c).


## 2. Materials and Methods

### 2.1. Hollow Spheres

Hollow spheres are manufactured and patented [32] by ATECA SAS. The constitutive material of the hollow spheres is composed of an epoxy resin with mineral powder to mainly prevent the spheres from sticking together during the manufacture process. It can also increase the structural rigidity and the material yield stress. Powder aggregates are in high proportion in the matrix. The grain size varies from 10 μm to 100 μm as it can be observed with the help of microscope observations; see Figure 1(c). The presence of this reinforcement increases the density of the material from 1200 kg·m−3 (resin only) to 1880 kg·m−3. Hollow spheres are manufactured by coating the resin on expanded polystyrene foam. They are then heated to polymerize the composite layer and to melt the polystyrene. Therefore, the spheres can be considered as hollow, the polystyrene remaining being negligible compared to the sphere’s shell. The hollow sphere diameter depends on the size of the polystyrene ball. It varies within 1≤Dext≤32 mm. The thickness ratio is defined by Rt=e/rm, where e is the shell thickness and rm the average radius (rm=rin+rext/2), in Figure 1(b).Tested hollow spheres present a diameter ofDext=29.7±0.1 mm, a radius ratio Rt=0.08, and a mass m=6.30±0.5 g. For clarity, only results obtained on one geometry are presented although experiments with different radius ratios (Rt=0.04) have been conducted.


### 2.2. Experimental Procedure

Uniaxial compressive tests on hollow spheres were conducted on two machines: one classical compression machine (Zwick Roell Z250) for quasi-static tests and an original fly wheel (Figure2) for dynamic tests [31]. The rotational movement of the wheel is transformed into a translational movement. The velocity of the moving plates can be considered as constant due to the high moment of inertia of the wheel. In real conditions of a bird strike, velocity up to 175 m·s−1 can be reported. The maximum macroscopic strain rate of the sandwich structure is then close to approximately 2000 s−1. The macroscopic strain rate of ≈50 s−1 achievable with the fly wheel is supposed to be a good indication to show the influence of a dynamic solicitation. Conventional compression Hopkinson’s bars have not been able to be directly used because of the high heterogeneity of the HSS. High velocity impacts are conducted with a 35 mm in diameter and 1 m length air cannon. An air booster is used to pressurize an air tank at 60 bars. Hollow spheres are placed at the entry of the cannon and are fired by releasing the pressure with a solenoid valve. The impact surface is made of steel and is considered perfectly rigid. Impact velocities up to 130 m·s−1 were measured by image analysis. Force sensors of 10 kN were used in quasi-static and dynamic compression tests. Force as a function of time for the high velocity impacts was measured with a 150 kN force sensor (Kistler 9377c). A high speed camera (Photron SA-5) is used to capture the hollow sphere crushing and the propagation of cracks. 75000 frames per second with a resolution of 320∗264 pixels is set for uniaxial compression tests. Due to the very short time of the high velocity impacts, the camera is set to 150000 frames per second with a resolution of 192∗184 pixels.Figure 2
Fly wheel device [31].To characterize the constitutive material of the hollow sphere, compression tests were conducted on a compression testing machine (Zwick Roell Z250) equipped with a climatic chamber. High temperature is reached by thermal resistance and low temperature is reached by adding liquid nitrogen. A control loop keeps the temperature constant during the tests. A 250 kN force sensor is used. Hopkinson bars (25 mm in diameter, aluminium) allow conducting high strain rate characterization [33, 34] at room temperature.


### 2.3. The Discrete Element Method

The DEM was originally developed for granular assemblies [35] but its applications have been expanded these last years to other areas such as tribology [36], fluid mechanics [37], and continuum mechanics [38]. The GranOO workbench [39, 40], developed in the laboratory, is naturally chosen. GranOO is based on an explicit time integration scheme (velocity Verlet scheme) [41] which permits a reasonable computational cost for high velocity simulations. Discrete element spatial position and velocity are calculated using the acceleration determined by the second Newton’s law of motion. The mass of each element is found so that it represents the density of the material modelled. Linear positions and orientations are, respectively, associated with resultant forces and torques. Quaternions are used for elements rotations instead of Euler angles for computing efficiency.Discrete elements can interact only by contact or can be connected by springs or 3D beams. As established in [29], DEM models using cohesive beams to link discrete elements are appropriate to model continuous materials. The idea is to build a hollow sphere made of thousands of discrete elements; see Figure 3. Appendix A gives the methodology to find the microscopic parameters (i.e., physical parameters of cohesive beams) of the model.Figure 3
Hollow sphere modelled with discrete elements and cohesive beams.The aptitude of the DEM to simulate the fracture of materials is possible by breaking the beams according to a failure criterion. The simplest fracture model is based on the Rankine stress criterion applied for each beam:(1)σ=12σmax+σmax2+4τmax2,where σ is the maximal equivalent Rankine stress, σmax the maximum normal stress, and τmax the maximal shear stress. But, as it was shown before in [42], this criterion failed to simulate accurately the failure of brittle material. A novel failure criterion was recently developed [42], no longer on each beam but on each discrete element. It consists in an equivalent Cauchy stress tensor (similar to what can be obtained with the finite element method) expressed at the center of each discrete element. The stress tensor is calculated taking into account the forces generated by all the beams connected to the element under consideration. It is computed as follows:(2)σ-i=12Ωi12∑jr→ij⊗f→ij+f→ij⊗r→ij,where(i)
⊗ is the tensor product between two vectors;(ii)
σ-i is the equivalent Cauchy stress tensor of the discrete element i;(iii)
Ωi is the volume of the discrete element i;(iv)
f→ij is the interaction force between i and j adjacent DE;(v)
r→ij is the relative position vector of the center of the two adjacent discrete elements i and j in a global basis BG.The equivalent Cauchy stress tensor can be diagonalized in order to find the principal stresses in the principal basisBG:(3)σ=σ11σ12σ13σ21σ22σ23σ31σ32σ33BG=σI000σII000σIIIBG.A failure criterion based on the Cauchy stress tensor can be defined, the maximum principal stress criterion, computed here:(4)σprincipal  stress=max⁡σI,σII,σIII.In order to take into account the time-temperature superposition principle, the strain rate has to be evaluated on each element. The strain calculation in DEM is not as direct as it can be for continuous modelling (e.g., finite element). The Lagrangian formulation for small deformationsεij=1/2∂ui/∂uj+∂uj/∂ui is indeed not well adapted in discrete domains. Specific models have been implemented to overpass this problem [43]. Most of these techniques rely on an equivalent discrete-continuum approach. Delaunay triangulation [44, 45] or Dirichlet tessellation [46] has made it possible to compute a local strain in discrete arrangements. The construction of a domain around the elements is nevertheless necessary which increases the computational time. A simpler method is proposed in this paper; it relies on the Cauchy stress tensor. Even if the stress tensor and the strain tensor in discrete domains show a strong duality, it is not complete since expression for both tensors does not belong to the same domain. The strain tensor takes into account the space between elements whereas the stress tensor is calculated for one element. Since only an equivalent strain rate is sought, the true strain is not essential. Similar reasoning was done for the DEM modelling of concrete [47] where an interaction strain between elements is calculated to implement a strain rate dependent model. In that case, it is possible to compute the equivalent strain in each principal direction:(5)εI=1EσI-νσII+σIII,εII=1EσII-νσI+σIII,εIII=1EσIII-νσI+σII.The strain rate is then calculated knowing the strain on each element at each iteration:(6)ε˙it=εit-εit-ΔtΔtwith Δt the time step.The failure criterion becomes(7)σprincipal  stress=max⁡σIε˙I,T,σIIε˙I,T,σIIIε˙I,T.The temperature is considered constant in the simulation; no heat conduction is performed.Once the element considered reaches the stress limit, the element breaks. All its beams are disabled and they do not interact anymore in the simulation. The broken element remains in the domain giving the right microscopic brittle behavior so that the energy conservation is verified. A brittle calibration shall be carried out in the same manner as the elastic calibration to find the limit stresses values.


## 2.1. Hollow Spheres

Hollow spheres are manufactured and patented [32] by ATECA SAS. The constitutive material of the hollow spheres is composed of an epoxy resin with mineral powder to mainly prevent the spheres from sticking together during the manufacture process. It can also increase the structural rigidity and the material yield stress. Powder aggregates are in high proportion in the matrix. The grain size varies from 10 μm to 100 μm as it can be observed with the help of microscope observations; see Figure 1(c). The presence of this reinforcement increases the density of the material from 1200 kg·m−3 (resin only) to 1880 kg·m−3. Hollow spheres are manufactured by coating the resin on expanded polystyrene foam. They are then heated to polymerize the composite layer and to melt the polystyrene. Therefore, the spheres can be considered as hollow, the polystyrene remaining being negligible compared to the sphere’s shell. The hollow sphere diameter depends on the size of the polystyrene ball. It varies within 1≤Dext≤32 mm. The thickness ratio is defined by Rt=e/rm, where e is the shell thickness and rm the average radius (rm=rin+rext/2), in Figure 1(b).Tested hollow spheres present a diameter ofDext=29.7±0.1 mm, a radius ratio Rt=0.08, and a mass m=6.30±0.5 g. For clarity, only results obtained on one geometry are presented although experiments with different radius ratios (Rt=0.04) have been conducted.


## 2.2. Experimental Procedure

Uniaxial compressive tests on hollow spheres were conducted on two machines: one classical compression machine (Zwick Roell Z250) for quasi-static tests and an original fly wheel (Figure2) for dynamic tests [31]. The rotational movement of the wheel is transformed into a translational movement. The velocity of the moving plates can be considered as constant due to the high moment of inertia of the wheel. In real conditions of a bird strike, velocity up to 175 m·s−1 can be reported. The maximum macroscopic strain rate of the sandwich structure is then close to approximately 2000 s−1. The macroscopic strain rate of ≈50 s−1 achievable with the fly wheel is supposed to be a good indication to show the influence of a dynamic solicitation. Conventional compression Hopkinson’s bars have not been able to be directly used because of the high heterogeneity of the HSS. High velocity impacts are conducted with a 35 mm in diameter and 1 m length air cannon. An air booster is used to pressurize an air tank at 60 bars. Hollow spheres are placed at the entry of the cannon and are fired by releasing the pressure with a solenoid valve. The impact surface is made of steel and is considered perfectly rigid. Impact velocities up to 130 m·s−1 were measured by image analysis. Force sensors of 10 kN were used in quasi-static and dynamic compression tests. Force as a function of time for the high velocity impacts was measured with a 150 kN force sensor (Kistler 9377c). A high speed camera (Photron SA-5) is used to capture the hollow sphere crushing and the propagation of cracks. 75000 frames per second with a resolution of 320∗264 pixels is set for uniaxial compression tests. Due to the very short time of the high velocity impacts, the camera is set to 150000 frames per second with a resolution of 192∗184 pixels.Figure 2
Fly wheel device [31].To characterize the constitutive material of the hollow sphere, compression tests were conducted on a compression testing machine (Zwick Roell Z250) equipped with a climatic chamber. High temperature is reached by thermal resistance and low temperature is reached by adding liquid nitrogen. A control loop keeps the temperature constant during the tests. A 250 kN force sensor is used. Hopkinson bars (25 mm in diameter, aluminium) allow conducting high strain rate characterization [33, 34] at room temperature.


## 2.3. The Discrete Element Method

The DEM was originally developed for granular assemblies [35] but its applications have been expanded these last years to other areas such as tribology [36], fluid mechanics [37], and continuum mechanics [38]. The GranOO workbench [39, 40], developed in the laboratory, is naturally chosen. GranOO is based on an explicit time integration scheme (velocity Verlet scheme) [41] which permits a reasonable computational cost for high velocity simulations. Discrete element spatial position and velocity are calculated using the acceleration determined by the second Newton’s law of motion. The mass of each element is found so that it represents the density of the material modelled. Linear positions and orientations are, respectively, associated with resultant forces and torques. Quaternions are used for elements rotations instead of Euler angles for computing efficiency.Discrete elements can interact only by contact or can be connected by springs or 3D beams. As established in [29], DEM models using cohesive beams to link discrete elements are appropriate to model continuous materials. The idea is to build a hollow sphere made of thousands of discrete elements; see Figure 3. Appendix A gives the methodology to find the microscopic parameters (i.e., physical parameters of cohesive beams) of the model.Figure 3
Hollow sphere modelled with discrete elements and cohesive beams.The aptitude of the DEM to simulate the fracture of materials is possible by breaking the beams according to a failure criterion. The simplest fracture model is based on the Rankine stress criterion applied for each beam:(1)σ=12σmax+σmax2+4τmax2,where σ is the maximal equivalent Rankine stress, σmax the maximum normal stress, and τmax the maximal shear stress. But, as it was shown before in [42], this criterion failed to simulate accurately the failure of brittle material. A novel failure criterion was recently developed [42], no longer on each beam but on each discrete element. It consists in an equivalent Cauchy stress tensor (similar to what can be obtained with the finite element method) expressed at the center of each discrete element. The stress tensor is calculated taking into account the forces generated by all the beams connected to the element under consideration. It is computed as follows:(2)σ-i=12Ωi12∑jr→ij⊗f→ij+f→ij⊗r→ij,where(i)
⊗ is the tensor product between two vectors;(ii)
σ-i is the equivalent Cauchy stress tensor of the discrete element i;(iii)
Ωi is the volume of the discrete element i;(iv)
f→ij is the interaction force between i and j adjacent DE;(v)
r→ij is the relative position vector of the center of the two adjacent discrete elements i and j in a global basis BG.The equivalent Cauchy stress tensor can be diagonalized in order to find the principal stresses in the principal basisBG:(3)σ=σ11σ12σ13σ21σ22σ23σ31σ32σ33BG=σI000σII000σIIIBG.A failure criterion based on the Cauchy stress tensor can be defined, the maximum principal stress criterion, computed here:(4)σprincipal  stress=max⁡σI,σII,σIII.In order to take into account the time-temperature superposition principle, the strain rate has to be evaluated on each element. The strain calculation in DEM is not as direct as it can be for continuous modelling (e.g., finite element). The Lagrangian formulation for small deformationsεij=1/2∂ui/∂uj+∂uj/∂ui is indeed not well adapted in discrete domains. Specific models have been implemented to overpass this problem [43]. Most of these techniques rely on an equivalent discrete-continuum approach. Delaunay triangulation [44, 45] or Dirichlet tessellation [46] has made it possible to compute a local strain in discrete arrangements. The construction of a domain around the elements is nevertheless necessary which increases the computational time. A simpler method is proposed in this paper; it relies on the Cauchy stress tensor. Even if the stress tensor and the strain tensor in discrete domains show a strong duality, it is not complete since expression for both tensors does not belong to the same domain. The strain tensor takes into account the space between elements whereas the stress tensor is calculated for one element. Since only an equivalent strain rate is sought, the true strain is not essential. Similar reasoning was done for the DEM modelling of concrete [47] where an interaction strain between elements is calculated to implement a strain rate dependent model. In that case, it is possible to compute the equivalent strain in each principal direction:(5)εI=1EσI-νσII+σIII,εII=1EσII-νσI+σIII,εIII=1EσIII-νσI+σII.The strain rate is then calculated knowing the strain on each element at each iteration:(6)ε˙it=εit-εit-ΔtΔtwith Δt the time step.The failure criterion becomes(7)σprincipal  stress=max⁡σIε˙I,T,σIIε˙I,T,σIIIε˙I,T.The temperature is considered constant in the simulation; no heat conduction is performed.Once the element considered reaches the stress limit, the element breaks. All its beams are disabled and they do not interact anymore in the simulation. The broken element remains in the domain giving the right microscopic brittle behavior so that the energy conservation is verified. A brittle calibration shall be carried out in the same manner as the elastic calibration to find the limit stresses values.


## 3. Numerical Calibration

Elastic and fracture parameters are used in the DEM model to reproduce the elastobrittle behavior of the material. These parameters are numerically calibrated according to the properties of the material: Young’s modulus, Poisson ratio, and the yield stress.


### 3.1. Experimental Compression Tests for Numerical Calibration

In order to calibrate the numerical model, preliminary quasi-static and dynamic compressive tests on the material’s hollow sphere are conducted. Compression cylindrical specimens were provided to determine directly the mechanical properties such as Young’s modulus and the failure force of this type of structure. These specimens were moulded and they underwent the same baking cycle as hollow spheres.To calculate Young’s modulus strain gauges were placed on cylindrical specimens of 30 mm in diameter and 20 mm in height to measure the true strain. The geometry was chosen to ensure a good measurement of the strain according to the strain gauges used. Young’s modulus is only evaluated in quasi-static test (the crosshead velocity was set at 5 mm·min−1) and at room temperature (T=20°C). A value of 10.3±0.3 GPa was found which is three times more than the pure epoxy resin. As expected, the reinforcements stiffen significantly the material while reducing the deformation at failure. The composite material showed a brittle behavior when compressed as expected for thermosetting polymer matrices at ambient temperature. After the elastic phase, samples present fractures (surfaces oriented at 90° to 45° as a function of the compression direction is observed, in Figure 6(a)). The elastic curve goes up to the failure of the material and no plasticity seems to occur. The fracture surfaces do not exhibit permanent deformation.Compressive tests in temperature were conducted on cylindrical specimen of 12.7 mm in diameter and 25.4 mm in height (ASTM D695). Temperatures range from −40°C to +80°C by step of 20°C. The crosshead velocity was set to 5 mm·min−1 (ε˙=0.003 s−1) and 500 mm·min−1 (ε˙=0.3 s−1).High strain rate tests (ε˙≈2000 s−1) were conducted on Hopkinson’s bars at room temperature. Cylindrical specimens of 10 mm in diameter and 5 mm in height were used.Results in terms of failure stress as a function of temperature and for different strain rate are presented in Figure4. The material shows a high sensibility to temperature: failure stress ranges from 50 MPa to more than 200 MPa between −40°C and 80°C for a strain rate of ε˙=0.003 s−1. The strain rate effect is also important; the failure stress in dynamic (ε˙≈2000 s−1) is more than twice as high compared to quasi-static compression.Figure 4
Failure stress of the material in compression for different strain rates and temperatures.In order to predict the yield stress of polymers, physical or phenomenological laws have been studied, especially for thermosetting resin [16, 48, 49]. The original Eyring law considers only the process of pure flow and neglects nonviscous effect (isotherm adiabatic heating or high elasticity). Validated for moderate strain rates [50], this theory is not suitable for very high strain rates. The generalized theory of Ree-Eyring [30] is more adapted and was successfully applied for epoxy resins [51]:(8)σyT=43kV0αν0αln⁡32ν0αJ0αε˙+QαRT+43kV0βν0βsinh⁡34ν0βJ0βε˙exp⁡QβRT-1=Aαln⁡2Cαε˙+QαRT+Aβsinh⁡Cβε˙exp⁡QβRT-1with(i)
σy the yield stress,(ii)
T the absolute temperature,(iii)
ε˙ the strain rate,(iv)
k the Boltzmann constant,(v)
R the universal gas constant,(vi)
V0α and V0β the activation volumes for each process of α and β,(vii)
ν0α and ν0β the elementary shear,(viii)
J0α and J0β the rate constants,(ix)
Aα,β, Cα,β, and Qα,β the activation energies for each process of α and β.Ree-Eyring constants are identified by numerical optimization. The Downhill Simplex Method [52] is used in Python to find the constants that will fit to the experimental data. Constants initial values of a similar material [51] are set before optimization. New values found are close to the initial ones and are presented in Table 1. The failure stress as a function of strain rate and temperature according to the Ree-Eyring law with the constants identified is plotted in Figure 5.Table 1
Ree-Eyring constants identified.
Parameters
Values
Units
A
α
0.019
MPa⋅K−1
C
α
1.13 × 10−31
s−1
Q
α
230
kJ⋅mol−1
A
β
0.010
MPa⋅K−1
C
β
7.06 × 10−5
s−1
Q
β
25.99
kJ⋅mol−1Figure 5
Failure stress as a function of strain rate and temperature according to the Ree-Eyring law. Experimental data is plotted as dark dots.Figure 6
Epoxy cylinder after failure in compression for experimental (a) and DEM simulation with principal stress criterion (b).
(a)
(b)


### 3.2. Brittle Calibration

Prior to the brittle calibration, the elastic calibration is conducted according to the methodology presented in AppendixB. This process remains the same for any type of elastic material and was successfully achieved before. Microscopic values found, corresponding to the material studied here, are shown in Appendix B.Brittle calibration is then carried out. The brittle parameter (stress limit) is scanned over a wide range of values in order to find the right microscopic value leading to the correct macroscopic value. The geometry tested was chosen to be the same as in the experiment (a cylinder withH=0.02 m, D=0.03 m is chosen). Firstly, the strain rate effect is not taken into account as only the correspondence between macroscopic and microscopic yield stress values is sought. Therefore, force displacement curves can be directly compared between experimental and numerical results. The sample is compressed between two plates. A very low loading velocity is applied to be in quasi-static condition (kinetic energy is negligible). The top and bottom faces of cylinders have been numerically smoothed so that numerical roughness (which is much larger than the cylindrical specimen at macroscopic scale) does not interfere with the simulation. To permit a crack closure and the sliding between cracked faces, the contact is activated between each element which are linked by a broken beam. The stiffness contact, computed by the penalty method, is equal to the rigidity of the beam in compression:(9)Econtact=EμSbeam2×rDE=Eμπrbeam22×rDE=EμπrDERμ22with Sbeam the beam section, Eμ Young’s modulus of the beam, rbeam the radius of the beam, rDE the mean radius of the two elements connected, and Rμ the radius ratio.Friction coefficients from 0.25 to 0.5 are encountered in literature for polymers; see [53] for pure epoxy in contact with steel. The failure force value varies from 1% with a fluctuation of the friction coefficient value between 0.3 and 0.5. The authors decide therefore to set a friction coefficient of 0.3 for both contacts between each discrete element (cracked faces) and between discrete elements and the plates. Static and dynamic friction coefficients are supposed to be the same as their values are close. The result with the principal stress criteria is plotted in Figure 7. The macroscopic yield stress is found to be 0.67×σμ. For example, at T=20°C and ε˙=0.003 s−1, σμ=195 MPa and σmacro=131 MPa. Qualitative results can be seen in Figure 6 where elements colour shows their horizontal displacement after failure (ΔL=1.25 mm). Failure crack paths are well represented by the DEM model: crack surfaces are oriented at 90° to 45°. The force decreases faster after the main failure which is due to the broken elements that are ejected from the crack surfaces. Indeed, the third body created by the fracture is thicker (minimum the size of the discrete elements) and more mobile in the simulation. This induces a less important stress take-up.Figure 7
Force displacement curve for brittle calibration.Numerical compression tests have been conducted at different velocities to validate the strain rate dependent model (based on the general theory of Ree-Eyring). Cylindrical geometries (12.7 mm in diameter and 25.4 mm in height) with 2000 discrete elements where used. This small number of elements does not give accurate results but allows decreasing the computational time to run low velocities compression tests. Results for both the strain rate dependent and strain rate independent models are presented in Figure8. Below 10 m·s−1, the failure force is quasi-constant for the strain rate independent model but increases accordingly to the Ree-Eyring law for the strain rate dependent model.Figure 8
Failure force for different compression velocities on a cylindrical specimen: comparison of the strain rate dependent and strain rate independent model.Above 10 m·s−1, the specimen is not in an equilibrium state, inertia effects are not negligible anymore and the failure stress increased for both models. The equilibrium time can be calculated with this equation:(10)tequi=π2ρL2Ewith E Young’s modulus, L the length of the specimen, and ρ the density of the material. A value of t≈34μs is found in this case. Considering that the failure is reached at a strain of 1.5%, the velocity threshold occurs at v=11 m·s−1 which is consistent with the results found for both models.


## 3.1. Experimental Compression Tests for Numerical Calibration

In order to calibrate the numerical model, preliminary quasi-static and dynamic compressive tests on the material’s hollow sphere are conducted. Compression cylindrical specimens were provided to determine directly the mechanical properties such as Young’s modulus and the failure force of this type of structure. These specimens were moulded and they underwent the same baking cycle as hollow spheres.To calculate Young’s modulus strain gauges were placed on cylindrical specimens of 30 mm in diameter and 20 mm in height to measure the true strain. The geometry was chosen to ensure a good measurement of the strain according to the strain gauges used. Young’s modulus is only evaluated in quasi-static test (the crosshead velocity was set at 5 mm·min−1) and at room temperature (T=20°C). A value of 10.3±0.3 GPa was found which is three times more than the pure epoxy resin. As expected, the reinforcements stiffen significantly the material while reducing the deformation at failure. The composite material showed a brittle behavior when compressed as expected for thermosetting polymer matrices at ambient temperature. After the elastic phase, samples present fractures (surfaces oriented at 90° to 45° as a function of the compression direction is observed, in Figure 6(a)). The elastic curve goes up to the failure of the material and no plasticity seems to occur. The fracture surfaces do not exhibit permanent deformation.Compressive tests in temperature were conducted on cylindrical specimen of 12.7 mm in diameter and 25.4 mm in height (ASTM D695). Temperatures range from −40°C to +80°C by step of 20°C. The crosshead velocity was set to 5 mm·min−1 (ε˙=0.003 s−1) and 500 mm·min−1 (ε˙=0.3 s−1).High strain rate tests (ε˙≈2000 s−1) were conducted on Hopkinson’s bars at room temperature. Cylindrical specimens of 10 mm in diameter and 5 mm in height were used.Results in terms of failure stress as a function of temperature and for different strain rate are presented in Figure4. The material shows a high sensibility to temperature: failure stress ranges from 50 MPa to more than 200 MPa between −40°C and 80°C for a strain rate of ε˙=0.003 s−1. The strain rate effect is also important; the failure stress in dynamic (ε˙≈2000 s−1) is more than twice as high compared to quasi-static compression.Figure 4
Failure stress of the material in compression for different strain rates and temperatures.In order to predict the yield stress of polymers, physical or phenomenological laws have been studied, especially for thermosetting resin [16, 48, 49]. The original Eyring law considers only the process of pure flow and neglects nonviscous effect (isotherm adiabatic heating or high elasticity). Validated for moderate strain rates [50], this theory is not suitable for very high strain rates. The generalized theory of Ree-Eyring [30] is more adapted and was successfully applied for epoxy resins [51]:(8)σyT=43kV0αν0αln⁡32ν0αJ0αε˙+QαRT+43kV0βν0βsinh⁡34ν0βJ0βε˙exp⁡QβRT-1=Aαln⁡2Cαε˙+QαRT+Aβsinh⁡Cβε˙exp⁡QβRT-1with(i)
σy the yield stress,(ii)
T the absolute temperature,(iii)
ε˙ the strain rate,(iv)
k the Boltzmann constant,(v)
R the universal gas constant,(vi)
V0α and V0β the activation volumes for each process of α and β,(vii)
ν0α and ν0β the elementary shear,(viii)
J0α and J0β the rate constants,(ix)
Aα,β, Cα,β, and Qα,β the activation energies for each process of α and β.Ree-Eyring constants are identified by numerical optimization. The Downhill Simplex Method [52] is used in Python to find the constants that will fit to the experimental data. Constants initial values of a similar material [51] are set before optimization. New values found are close to the initial ones and are presented in Table 1. The failure stress as a function of strain rate and temperature according to the Ree-Eyring law with the constants identified is plotted in Figure 5.Table 1
Ree-Eyring constants identified.
Parameters
Values
Units
A
α
0.019
MPa⋅K−1
C
α
1.13 × 10−31
s−1
Q
α
230
kJ⋅mol−1
A
β
0.010
MPa⋅K−1
C
β
7.06 × 10−5
s−1
Q
β
25.99
kJ⋅mol−1Figure 5
Failure stress as a function of strain rate and temperature according to the Ree-Eyring law. Experimental data is plotted as dark dots.Figure 6
Epoxy cylinder after failure in compression for experimental (a) and DEM simulation with principal stress criterion (b).
(a)
(b)


## 3.2. Brittle Calibration

Prior to the brittle calibration, the elastic calibration is conducted according to the methodology presented in AppendixB. This process remains the same for any type of elastic material and was successfully achieved before. Microscopic values found, corresponding to the material studied here, are shown in Appendix B.Brittle calibration is then carried out. The brittle parameter (stress limit) is scanned over a wide range of values in order to find the right microscopic value leading to the correct macroscopic value. The geometry tested was chosen to be the same as in the experiment (a cylinder withH=0.02 m, D=0.03 m is chosen). Firstly, the strain rate effect is not taken into account as only the correspondence between macroscopic and microscopic yield stress values is sought. Therefore, force displacement curves can be directly compared between experimental and numerical results. The sample is compressed between two plates. A very low loading velocity is applied to be in quasi-static condition (kinetic energy is negligible). The top and bottom faces of cylinders have been numerically smoothed so that numerical roughness (which is much larger than the cylindrical specimen at macroscopic scale) does not interfere with the simulation. To permit a crack closure and the sliding between cracked faces, the contact is activated between each element which are linked by a broken beam. The stiffness contact, computed by the penalty method, is equal to the rigidity of the beam in compression:(9)Econtact=EμSbeam2×rDE=Eμπrbeam22×rDE=EμπrDERμ22with Sbeam the beam section, Eμ Young’s modulus of the beam, rbeam the radius of the beam, rDE the mean radius of the two elements connected, and Rμ the radius ratio.Friction coefficients from 0.25 to 0.5 are encountered in literature for polymers; see [53] for pure epoxy in contact with steel. The failure force value varies from 1% with a fluctuation of the friction coefficient value between 0.3 and 0.5. The authors decide therefore to set a friction coefficient of 0.3 for both contacts between each discrete element (cracked faces) and between discrete elements and the plates. Static and dynamic friction coefficients are supposed to be the same as their values are close. The result with the principal stress criteria is plotted in Figure 7. The macroscopic yield stress is found to be 0.67×σμ. For example, at T=20°C and ε˙=0.003 s−1, σμ=195 MPa and σmacro=131 MPa. Qualitative results can be seen in Figure 6 where elements colour shows their horizontal displacement after failure (ΔL=1.25 mm). Failure crack paths are well represented by the DEM model: crack surfaces are oriented at 90° to 45°. The force decreases faster after the main failure which is due to the broken elements that are ejected from the crack surfaces. Indeed, the third body created by the fracture is thicker (minimum the size of the discrete elements) and more mobile in the simulation. This induces a less important stress take-up.Figure 7
Force displacement curve for brittle calibration.Numerical compression tests have been conducted at different velocities to validate the strain rate dependent model (based on the general theory of Ree-Eyring). Cylindrical geometries (12.7 mm in diameter and 25.4 mm in height) with 2000 discrete elements where used. This small number of elements does not give accurate results but allows decreasing the computational time to run low velocities compression tests. Results for both the strain rate dependent and strain rate independent models are presented in Figure8. Below 10 m·s−1, the failure force is quasi-constant for the strain rate independent model but increases accordingly to the Ree-Eyring law for the strain rate dependent model.Figure 8
Failure force for different compression velocities on a cylindrical specimen: comparison of the strain rate dependent and strain rate independent model.Above 10 m·s−1, the specimen is not in an equilibrium state, inertia effects are not negligible anymore and the failure stress increased for both models. The equilibrium time can be calculated with this equation:(10)tequi=π2ρL2Ewith E Young’s modulus, L the length of the specimen, and ρ the density of the material. A value of t≈34μs is found in this case. Considering that the failure is reached at a strain of 1.5%, the velocity threshold occurs at v=11 m·s−1 which is consistent with the results found for both models.


## 4. Results

### 4.1. Compressive Tests on Hollow Spheres

Quasi-static (Figure9) and dynamic compression tests on hollow spheres have been conducted and compared in Figure 10. The dispersion is plotted as the area between the maximum and minimum force curve during compression. As expected, the structure of a hollow sphere presents strain rate dependence. The increase in failure force is more than 50% higher in dynamic compression (v≈1.3 m·s−1). A significant dispersion is observed both in quasi-static and dynamic compression. The standard deviation (SD) for the failure force is 110 N in quasi-static and 160 N in dynamic compression. And for the displacement at failure, it is, respectively, 0.24 mm and 0.20 mm. This dispersion may have two origins: dispersion associated with the accuracy of force and displacement sensors and the dispersion due to geometrical defaults of hollow spheres induced by the manufacturing process. The two first causes can be quantified. In quasi-static compression, the compressive speed is the same and less than 1% accurate. In dynamic compression, the measurement is estimated to be 5% accurate. Geometrical defaults of hollow spheres are important. Indeed, a mass dispersion of about ±10% can be measured and hollow spheres sphericity is not perfect. The failure force is therefore highly dependent of the orientation of the sphere between plates.Figure 9
Quasi-static compression of hollow spheres (not all tests are presented).Figure 10
Quasi-static and dynamic compression of hollow spheres.The failure process of the sphere remains the same in quasi-static or dynamic compression. It can be described in three stages:(1)
Elastic phase: at the beginning of the compression (ΔL/D=0–0.04 of relative longitudinal displacement, where ΔL is the displacement of the moving plate and D the diameter of the sphere), a linear elastic behavior can be considered as no drop of force occurs.(2)
Failure: a great drop of force is observed forΔL/D=0.04–0.06 which corresponds to rapid crack propagation. The circumferential stress in the hollow sphere’s shell is mainly positive (tensile stress) far from the point of load application when the structure is under uniaxial compression [54]. The opening mode I of the crack is then expected and is confirmed here. Before failure, the structure reaches its maximum force which is on average 1037 ± 160 N in dynamic and 752 ± 110 N in quasi-static compression. The crack propagates vertically (between the two plates) on a half perimeter, in Figures 9 and 11. The F(δ) curves plotted in the paper only represent the first two phases. The displacement of the plates was stopped after the main failure.(3)
Multiple cracks: once the hollow sphere is vertically cracked, the resultant force is close to zero. Then walls of the sphere are bending and buckling causing the proliferation of cracks oriented in every direction.Figure 11
Dynamic fracture of a hollow sphere: experiment (a) and simulation (b).
(a)
(b)


### 4.2. Numerical Modelling of Hollow Spheres

Once the hollow sphere was created (see AppendixC for more details), micro-parameters found in the elastic and brittle calibration were used (Table 2). The velocity of the moving plate was set to 1.3 m·s−1, close to the average impact velocity for experimental tests. The principal stress criterion with the strain rate dependence was employed and compared to the experimental results. The force displacement curve of the dynamic experimental data and the numerical result is plotted in Figure 12. The force at failure is 1010 N, close to the experimental results. The displacement at failure is 1.29 mm, 6% more, and the energy at failure is 0.71 J, 9% less than in the experiments.Table 2
Macroscopic and microscopic elastic parameters.
Young’s modulus
Poisson ratio
Radius ratio
Macroscopic
E
M
=
10.3 GPa
ν
M
=
0.3
None
Microscopic
E
μ
=
356 GPa
ν
μ
=
0.3
R
μ
=
0.3Figure 12
Comparison between experimental and numerical (principal stress criterion) results in dynamic compression.In addition to quantitative results, the numerical modelling permits to obtain excellent qualitative results. The crack propagates in the same way: the crack is frank, rectilinear, and directed vertically between the two plates (Figure13). The mean crack tip velocity computed is approximately 200 m·s−1, close to the experimental results (a˙≈215 m·s−1). The position of the crack tip as a function of time was measured directly on experimental and numerical images.Figure 13
Zoom in the microcracks before the main crack on experimental data (a) and on numerical results (b) where broken elements are in yellow.
(a)
(b)Microcracks are also observed at the top and bottom of the hollow sphere; see Figure13. It shows that, before the main crack, hollow spheres are subjected to permanent deformation and damage. The numerical model takes into account microcracks and allows quantifying them, one can notice the number of broken elements associated with a small drop of force which is not observable with experimental data, in Figure 14. It is now possible to analyse the position of the microcracks from inside the hollow sphere. Cracks are concentrated at the bottom and top center from inside the hollow sphere but are more spread at the outside face, as one can see the experimental data. Indeed, hollow sphere’s shells are bending and microcracks are observed where tensile stress occurs.Figure 14
F
(
δ
) curve with the broken discrete element counting during compression.


### 4.3. Validation for High Strain Rate

Eventually, high velocity tests are conducted with an air cannon. Qualitative results of both the strain rate dependent and strain rate independent model compared to the experimental results are presented in Figure15. In the experiment, 33 μs after impact, several cracks (8 to 10) are observed near the point of impact. Cracks propagate as far as 2/3 of the hollow sphere (66 μs). These scraps of spheres then bend and break (99 μs). The remaining dome crashes in turn. The strain rate independent model fails to reproduce the crushing of the hollow sphere at high strain rate: the hollow sphere breaks instantly when it impacts. On the contrary, the strain rate dependent model takes into account the different phases of crushing.Figure 15
Hollow sphere impacting a rigid surface at 124 m⋅s−1. From top to bottom: experimental results, strain rate independent, and strain rate dependent models.The impact force has been experimentally measured on a reduced number of tests. Indeed, an impact location outside the center of the force sensor caused nonaccurate results that cannot be exploited. Figure16 presents the impact force (in red) of a hollow sphere impacting a rigid surface at 70 m·s−1.Figure 16
Hollow sphere impacting a rigid surface at 70 m⋅s−1. Comparison of the force of impact as a function of time for the experimental and numerical results. Raw and smooth (by moving average) numerical results are presented.The DEM model is then used to reproduce the impact at the same velocity. The noisy computed force is due to the progressive impact of all the discrete elements. For clarity of results, a moving average technique is conducted for the numerical data. Raw results (in transparency) and smooth results are presented in Figure16. The maximal force is close for both the strain rate dependent model and the experiment (respectively, 3350 N and 3696 N for the smooth results). The decrease of force after the peak is similar until 150 μs after impact. The strain rate independent model has a lower maximal force (2444 N, -34% compared to the experimental results).The numerical results are in good agreement with the experimental data and justify the use of a strain rate dependent law.


## 4.1. Compressive Tests on Hollow Spheres

Quasi-static (Figure9) and dynamic compression tests on hollow spheres have been conducted and compared in Figure 10. The dispersion is plotted as the area between the maximum and minimum force curve during compression. As expected, the structure of a hollow sphere presents strain rate dependence. The increase in failure force is more than 50% higher in dynamic compression (v≈1.3 m·s−1). A significant dispersion is observed both in quasi-static and dynamic compression. The standard deviation (SD) for the failure force is 110 N in quasi-static and 160 N in dynamic compression. And for the displacement at failure, it is, respectively, 0.24 mm and 0.20 mm. This dispersion may have two origins: dispersion associated with the accuracy of force and displacement sensors and the dispersion due to geometrical defaults of hollow spheres induced by the manufacturing process. The two first causes can be quantified. In quasi-static compression, the compressive speed is the same and less than 1% accurate. In dynamic compression, the measurement is estimated to be 5% accurate. Geometrical defaults of hollow spheres are important. Indeed, a mass dispersion of about ±10% can be measured and hollow spheres sphericity is not perfect. The failure force is therefore highly dependent of the orientation of the sphere between plates.Figure 9
Quasi-static compression of hollow spheres (not all tests are presented).Figure 10
Quasi-static and dynamic compression of hollow spheres.The failure process of the sphere remains the same in quasi-static or dynamic compression. It can be described in three stages:(1)
Elastic phase: at the beginning of the compression (ΔL/D=0–0.04 of relative longitudinal displacement, where ΔL is the displacement of the moving plate and D the diameter of the sphere), a linear elastic behavior can be considered as no drop of force occurs.(2)
Failure: a great drop of force is observed forΔL/D=0.04–0.06 which corresponds to rapid crack propagation. The circumferential stress in the hollow sphere’s shell is mainly positive (tensile stress) far from the point of load application when the structure is under uniaxial compression [54]. The opening mode I of the crack is then expected and is confirmed here. Before failure, the structure reaches its maximum force which is on average 1037 ± 160 N in dynamic and 752 ± 110 N in quasi-static compression. The crack propagates vertically (between the two plates) on a half perimeter, in Figures 9 and 11. The F(δ) curves plotted in the paper only represent the first two phases. The displacement of the plates was stopped after the main failure.(3)
Multiple cracks: once the hollow sphere is vertically cracked, the resultant force is close to zero. Then walls of the sphere are bending and buckling causing the proliferation of cracks oriented in every direction.Figure 11
Dynamic fracture of a hollow sphere: experiment (a) and simulation (b).
(a)
(b)


## 4.2. Numerical Modelling of Hollow Spheres

Once the hollow sphere was created (see AppendixC for more details), micro-parameters found in the elastic and brittle calibration were used (Table 2). The velocity of the moving plate was set to 1.3 m·s−1, close to the average impact velocity for experimental tests. The principal stress criterion with the strain rate dependence was employed and compared to the experimental results. The force displacement curve of the dynamic experimental data and the numerical result is plotted in Figure 12. The force at failure is 1010 N, close to the experimental results. The displacement at failure is 1.29 mm, 6% more, and the energy at failure is 0.71 J, 9% less than in the experiments.Table 2
Macroscopic and microscopic elastic parameters.
Young’s modulus
Poisson ratio
Radius ratio
Macroscopic
E
M
=
10.3 GPa
ν
M
=
0.3
None
Microscopic
E
μ
=
356 GPa
ν
μ
=
0.3
R
μ
=
0.3Figure 12
Comparison between experimental and numerical (principal stress criterion) results in dynamic compression.In addition to quantitative results, the numerical modelling permits to obtain excellent qualitative results. The crack propagates in the same way: the crack is frank, rectilinear, and directed vertically between the two plates (Figure13). The mean crack tip velocity computed is approximately 200 m·s−1, close to the experimental results (a˙≈215 m·s−1). The position of the crack tip as a function of time was measured directly on experimental and numerical images.Figure 13
Zoom in the microcracks before the main crack on experimental data (a) and on numerical results (b) where broken elements are in yellow.
(a)
(b)Microcracks are also observed at the top and bottom of the hollow sphere; see Figure13. It shows that, before the main crack, hollow spheres are subjected to permanent deformation and damage. The numerical model takes into account microcracks and allows quantifying them, one can notice the number of broken elements associated with a small drop of force which is not observable with experimental data, in Figure 14. It is now possible to analyse the position of the microcracks from inside the hollow sphere. Cracks are concentrated at the bottom and top center from inside the hollow sphere but are more spread at the outside face, as one can see the experimental data. Indeed, hollow sphere’s shells are bending and microcracks are observed where tensile stress occurs.Figure 14
F
(
δ
) curve with the broken discrete element counting during compression.


## 4.3. Validation for High Strain Rate

Eventually, high velocity tests are conducted with an air cannon. Qualitative results of both the strain rate dependent and strain rate independent model compared to the experimental results are presented in Figure15. In the experiment, 33 μs after impact, several cracks (8 to 10) are observed near the point of impact. Cracks propagate as far as 2/3 of the hollow sphere (66 μs). These scraps of spheres then bend and break (99 μs). The remaining dome crashes in turn. The strain rate independent model fails to reproduce the crushing of the hollow sphere at high strain rate: the hollow sphere breaks instantly when it impacts. On the contrary, the strain rate dependent model takes into account the different phases of crushing.Figure 15
Hollow sphere impacting a rigid surface at 124 m⋅s−1. From top to bottom: experimental results, strain rate independent, and strain rate dependent models.The impact force has been experimentally measured on a reduced number of tests. Indeed, an impact location outside the center of the force sensor caused nonaccurate results that cannot be exploited. Figure16 presents the impact force (in red) of a hollow sphere impacting a rigid surface at 70 m·s−1.Figure 16
Hollow sphere impacting a rigid surface at 70 m⋅s−1. Comparison of the force of impact as a function of time for the experimental and numerical results. Raw and smooth (by moving average) numerical results are presented.The DEM model is then used to reproduce the impact at the same velocity. The noisy computed force is due to the progressive impact of all the discrete elements. For clarity of results, a moving average technique is conducted for the numerical data. Raw results (in transparency) and smooth results are presented in Figure16. The maximal force is close for both the strain rate dependent model and the experiment (respectively, 3350 N and 3696 N for the smooth results). The decrease of force after the peak is similar until 150 μs after impact. The strain rate independent model has a lower maximal force (2444 N, -34% compared to the experimental results).The numerical results are in good agreement with the experimental data and justify the use of a strain rate dependent law.


## 5. Discussion and Conclusions

The fracture behavior of composite hollow sphere under quasi-static and dynamic compression was investigated using an experimental and a numerical approach. Experiments have been conducted at room temperature and at different compressive velocities (5 mm·min−1 to 130 m·s−1). A numerical model based on the discrete element method which takes into account the strain rate dependence is proposed.Compression tests on composite hollow spheres have shown that the compression velocity influences the response of the material/structure. As expected for polymer materials [51, 55], the more the compressive velocity, the more the force at failure (+50%) and the energy at failure (+30%). The material is subjected to dynamic fracture to dissipate the free available energy stored in the structure: a principal crack propagates at high speed (a˙≈215 m·s−1) and initiates near the sphere summits.The numerical model using DEM has been developed to simulate the dynamic fracture and to consider the time-temperature superposition principle of the polymer. The constants of the Ree-Eyring law are identified according to compression tests at different strain rates (ε˙=0.003 s−1 to ε˙=0.3 s−1) and temperatures (−40°C to +80°C). The strain rate is successfully computed in DEM through the stress tensor. By taking into account the strain rate dependence on the yield stress of the material, the model allows simulating accurately the fracture of the material for dynamic compression, both on a cylindrical specimen and on one hollow sphere. The maximal force numerically measured before the principal crack propagates in the hollow is close to the experimental data (difference of 3%).The model highlights the phenomena of microcracks coalescence prior to the main failure. It makes it possible to quantify the energy dissipated by fracture before and after the failure of the hollow sphere. Moreover, the DEM shows all its potential for high velocity impacts (v=50 m·s−1 to 130 m·s−1) where multicracks, large displacements, and contacts contribute majorly to the results. It is shown that the use of a strain rate dependent law is essential to properly simulate the ruin of the structure. This study highlights a DEM approach to model a dynamic fracture mechanism with taking into account the time dependency of polymer materials.Eventually, the objective is to develop a macroscopic model of the hollow spheres crushing. Thousands of hollow spheres will interact between each other. Each sphere will be modelled by one discrete element and the interaction between each element will contain the constitutive law identified with the model presented in this paper. A promising multiscale model will eventually take advantage of the DEM to simulate fracture on the microscopic scale (one to five hollow spheres) and to describe the mechanical behavior of the whole structure on another higher scale. Optimization could be performed at the macroscopic scale by finding the right microscopic parameters to dissipate more energy during an impact.



---

*Source: 7638482-2017-06-06.xml*