1. Introduction
Glass has been increasingly used as load bearing as well as non-load bearing members in engineering due to its transparent appearances and appealing functions. Among all of the glass types, the soda-lime-silica glass is the most prevalent, and its typical composition is 70% SiO
2, 10% CaO and 15% Na
2O [
1]. However, these glass components are brittle and vulnerable to impacts [
2]. Under an impact action, glass may fracture into pieces, resulting in catastrophic sequences. Moreover, the flying shards threaten the users greatly and may injury them seriously. Thus, the impact fracture and the post-failure behaviour of glass have attracted growing attention from both the academic and the engineering communities. Endeavours have been made to examine the impact failure of glass for a long time. Early experimental work can be traced to 1920s when Preston [
3] revealed that the breakage behaviour of annealed glass is directly related to the appearance of fracture surfaces. This opinion was further backed by Shand [
4]. Baird [
5] performed a test with a high-speed camera and investigated the stress waves travelling in a glass rod. Hu et al. [
6] provided test data for the spherical hard impact fracture of a rectangular glass plate with a polycarbonate backing. Kozłowski [
7] presented experimental results on glass balustrades subjected to soft double-tire impact in both the intact and the post-failure (with one glass ply damaged) states. The acceleration of the projectile, stress in glass, and the displacement of balustrades were revealed.
To record the impact fracture processes of glass, high-speed cameras and sensitive sensors are necessary. Further, experiments are expensive and not easy to perform. In view of this, numerical simulations have become popular. Among all of the computational approaches, the finite element method (FEM) and the discrete element method (DEM) are the most commonly used. Repetto et al. [
8] presented an FE model to evaluate the radial cracking in glass rods, allowing for the crack initiation and propagation in a tension-shear cohesive fashion. Sun et al. [
9] implemented a continuum damage mechanics (CDM) model into the commercial FEM package ABAQUS to solve the impact fracture behaviour of glass plies, and the resistance of glass subjects to stone impact was investigated. Grujicic et al. [
10] proposed an impact material model considering the high-strain-rate and high-pressure for soda-lime glass, and the model was embedded into the ABAQUS/Explicit. With the developed model, edge-on impact (EOI) on soda-lime glass with a spherical projectile was simulated. In Mohagheghian et al. [
11], the responses of chemically strengthened laminated glass window plates subjected to low velocity soft impacts were simulated with ABAQUS/Explicit, showing the energy absorption and peak force reduction. In Kozłowski et al. [
12,
13], 3D numerical studies on the behaviour of glass subjected to hard-body impact using ABAQUS were performed, reproducing the experiments and retrospectively determining the maximum principal stress in the glass during impact. In general, the FEM can predict the crack initiation and the propagation with reasonable accuracy. However, the crack branching and the post-failure behaviour are highly difficult to simulate with the FEM due to its intrinsic disadvantages.
In addition to the applications with the FEM, the DEM has also been employed in the glass impact fracture analysis. Spherical discrete elements were used by Zang et al. [
14], and fracture and fragmentation behaviours of monolithic and laminated glass under rigid ball impacts were investigated. Xu et al. [
15] proposed an adaptive combined DE/FE algorithm to simulate the fracture process of brittle material with plane stress assumption. The discrete elements were generated automatically in severely deformed regions, and the model was validated through glass impact fracture simulation. In Baraldi et al. [
16], the non-linear discrete element modelling approach was employed to predict the load-deflection relation of laminated glass from the initial cracking up to the collapse. A recent study [
17] examined the glass fracture and the interfacial debonding in laminated glass under impact using a cohesive DE/FE model.
Besides the FEM and the DEM, other numerical methods such as the Extended Finite Element Method (XFEM) [
18] and the Peridynamics [
19,
20] have also been employed in the glass fracture analyses. The combined finite-discrete element method (FDEM) [
21] is an extension of both the FEM and the DEM. It was proposed and developed by Munjiza in the 1990s [
22]. The structures are fully discretised into a number of discrete elements, in which an FE formulation is incorporated. Thus, accurate evaluation of the contact forces and the structural deformation can be obtained, leading to a rational manner of crack development. In this study, the FDEM program ‘Y’ using the Munjiza-NBS contact detection algorithm [
23] and a cohesive cracking model [
24] was employed to simulate the impact fracture of glass. Cracks are assumed to initiate and propagate along element edges [
25], and mesh bias can be relieved should fine mesh with irregular orientations be used [
26]. Further information on the FDEM can be referred to in the work of Munjiza and his co-workers [
27,
28]. The FDEM has previously been applied to the quasi-static and dynamic analysis of brittle solids [
29,
30,
31,
32,
33]. Early work [
34,
35] addressed the impact fracture simulations on the monolithic and the laminated glass with the FDEM. Hard body impact damage modes like flexural cracking, Hertzian cone failure, and punching were obtained by Chen et al. [
36] with the FDEM, and the breakage regimes were identified. Subsequent research addressed the capability of the FDEM in modelling the failure of laminated glass subjected to both hard and soft body impacts [
37,
38]. In Munjiza et al. [
39], impact fracture and the post-failure of 3D glass shells were presented using the FDEM elements. A recent paper [
40] introduced a generalised traction-separation model and simulated the impact fracture process of glass beams and plates.
Although there are some applications on the glass impact fracture using the FDEM, a comprehensive parametric study is still in demand and warranted. As a sequel to the previous work [
36], a parametric investigation into the impact fracture of monolithic glass (especially the soda-lime glass) using the FDEM was conducted in this study, and an insight into the glass impact fracture mechanism was provided. The emphasis of this research is given to the annealed glass (float glass). The layout for the rest of this paper is as follows. In
Section 2, the fundamentals of the FDEM are briefly introduced. The glass fracture model used in the FDEM simulation is addressed in
Section 3.
Section 4 presents numerical examples for verification purposes. In
Section 5, a convergence study on the impact fracture patterns, projectile velocity and the total kinetic energy of different meshes is conducted. Based on the selected base case, a comprehensive investigation was performed over a variety of parameters, including the impact velocity, the impact angle, the material properties of glass, etc. Finally, concluding remarks are reached in
Section 6.
3. Fracture Model of Glass
Basically, glass is ideally brittle and follows Mode I fracture [
41]. To account for the rupture of glass computationally, a cohesive fracture model based on Hillerborg’s [
42] is employed in this study. With this model, the brittleness of glass can be described, and the computational stability can be guaranteed. A crack is considered to develop once the bonding stress
σ reaches the tensile strength
ft, indicating the start of the damage. Meanwhile,
δ reaches its elastic limit, i.e.,
δ =
δp. With the increase of the deformation of the joint element, the bonding stress
σ decreases gradually. When
δ equals its ultimate value
δc,
σ = 0 and the cracking process is completed. This process is considered as ‘strain softening’, and the area between the softening curve and the coordinate axes equals the fracture energy G, as given by Equation (7).
Prior to reaching the tensile strength of the material, standard continuum formulation (FEM) is adopted. The cohesive fracture model is capable of capturing the damage through a weak constitutive law, enabling the cracks to be smeared out over the continuum. The breakage of glass fully depends on the stress distribution of the material since no pre-defined cracks or notches are necessary.
The complete relation between the bonding stress
σ and the joint element deformation
δ is given in Equation (8) according to Munjiza [
21].
The heuristic parameter
z has the form of
where
D is a fracture index within the interval [0, 1].
D = 0 suggests that the joint element is undeformed and intact, while
D = 1 implies the total damage of the joint element and the free movement of corresponding discrete element couples. The expression of index
D is given in Equation (10).
In Equation (9), parameters
a,
b and
c are constants. It is verified that, for glass,
a = 1.2,
b = −1.0 and
c = 1.0 [
36]. Based on the above definitions, the normalised strain softening curve of glass is given in
Figure 3.
6. Concluding Remarks
The FDEM has been used in this paper for parametric investigation of both the oblique and perpendicular impact fracture of glass. The basic theories of the FDEM and the cohesive fracture model of glass were briefly introduced. Numerical examples were shown for verification purposes.
A convergence study on the impact responses with different meshes was performed in
Section 5.1, demonstrating the reliability and robustness of the FDEM in simulating the impact fracture of glass. In
Section 5.2, four different types of damage, namely the local minor, bending, cone and punching failure, were captured when the impact velocity varied. The correlations between the impact velocity and the impact failure modes were revealed. The obtained data points agree with the established breakage regimes [
36] very well. The influence of the impact angle was examined in
Section 5.3, and it was concluded that the damage becomes more and more serious with the increase of the impact angle. Glass material properties like the tensile strength and fracture energy play parts in the impact failure responses, and their influences were studied in
Section 5.4. Different failure modes such as cone, half-cone, bending, spalling, etc. were obtained, and their relevance to glass material properties were summarised in
Table 7 and
Table 8. The threshold velocity for creating a dominant through-thickness bending crack against impact angles is investigated in
Section 5.5. The obtained threshold velocity curve is in a clear descending trend, and it fluctuates within a narrow band after
α = 35°.
Based on the parametric investigation in this study, guidance can be suggested for the design and manufacturing of glass, and comparative benchmark examples for both the oblique and the perpendicular impacts are also provided for future research.