An improved interfacial bonding model for material interface modeling
Introduction
The bulk mechanical behavior of a material is largely determined by its microstructures, e.g. grain and grain boundaries in polycrystalline materials, bi-material interface in composite materials. For instance, in natural staggered composites such as bone and nacre, the brittle tablets are bonded by softer interfaces capable of dissipating a significant amount of energy [1], which makes the material remarkably strong and tough. Therefore, how to characterize and model these interfacial zones has been the focus of intense research. A significant amount of research efforts are being dedicated to develop interfacial zone models, mainly in the simulation of materials failure. The coupled atomistic/continuum interface zone models were developed for the analysis of dynamically propagating crack of interfaces [2], [3], [4]. Gao and Klein [5] proposed a virtual internal bond (VIB) model with randomized cohesive interactions between material particles. This VIB model incorporates an atomic cohesive force law into the constitutive model of materials for modeling deformation and failure in the interfacial region. Additionally, different cohesive zone models were developed to describe interfacial behaviors for material failure analysis [6], [7], [8], [9], [10].
Among different cohesive zone models, the exponential cohesive zone model [7] is one of the most popular interfacial zone models. The constitutive relationship of cohesive fracture is described by a potential in the model, which characterizes the physical debonding behavior. Although the exponential cohesive zone model has gained much popularity in material failure simulations, it has several limitations. It is often based on the assumption that the normal fracture energy equals the shear fracture energy [11], [12], [13]. This assumption is often not consistent with the experiment proof. In fact, multiple experimental studies indicated that the fracture energies in modes I and mode II are different, e.g. Araki et al. [14], Benzeggagh and Kenane [15], Dollhofer et al. [16], Pang [17], Warrior et al. [18] and Yang et al. [19]. Furthermore, when the interface is under a large normal compression condition, the maximum shear traction will become negative and this does not appear to be realistic [20]. In addition, with large tangential separation, the maximum normal repulsive traction will decrease to zero, which might result in surface penetration of two contact surfaces under large compressive displacement.
The original exponential cohesive zone model [7] has been extended and altered by many researchers for different applications. An irreversible exponential cohesive zone model that uses an effective opening displacement was developed by Ortiz and Pandolfi [21] to consider different ratio of tractions along sliding and normal direction under mixed-mode failure. Later, Zhang and Paulino [22] extended the original exponential cohesive zone model to functionally graded materials (FGMs) modeling, which considers the influence of material gradation on crack initiation in mixed-mode fracture problem. Van den Bosch et al. [20] adopted the original exponential cohesive zone model as a mixed-mode exponential cohesive zone model with different normal and shear fracture energy. Recently, Zeng and Li [23] developed a multiscale cohesive zone model, in which the cohesive laws can be obtained from atomic lattice structures.
In this paper, an improved interfacial bonding model was proposed to address the aforementioned problems and to study the interfacial interactions in biological materials, especially to model the organic interface of extrafibrillar matrix in bone. The mechanical responses of bone not only depend on its microstructure, but also depend on different loading conditions [24]. Current experimental studies on bone fracture are mainly in tension or bending tests because it is easy to conduct those tests. Only limited knowledge is available on the mechanical response of bone under compressive loading. In fact, bones in life are usually loaded in compression although they can fail at any loading direction [25], [26]. A few experiments have been conducted to understand the mechanism of bone failure in compression [24], [27], [28]. However, due to the complex character of bone failure under compression [24], [25], it is difficult to pinpoint the key characteristic of bone failure under compression, e.g. shear damage or slippage interaction between collagen and mineral phase causing the irreversible deformation. Therefore, it is necessary to develop numerical models to study bone failure mechanisms under compressive loading. This model was developed from a potential function and it characterizes different potentials/fracture energies, different interfacial strengths and describes attractive and repulsive behaviors of interfacial interactions. This improved interfacial bonding model not only preserves all essential features of an improved exponential cohesive zone model [20], it is also physically realistic with interface under both tension and compression condition. Furthermore, the proposed interfacial bonding model was verified through a simple compression test in a standard hexagonal structure. Ultimately, the proposed interfacial bonding model was employed to study the mechanical behavior of the extrafibrillar matrix in bone.
The paper is organized in seven Sections: in Section 2, the traditional exponential cohesive zone model was reviewed; in Section 3, the improved interfacial bonding model was developed; in Section 4, path dependence of work-of-separation of the improved interfacial bonding model was studied; Section 5 verified the proposed model by analytical solutions; in Section 6, a fracture simulation of extrafibrillar matrix in bone was presented; and Section 7 concluded the present work.
Section snippets
Exponential cohesive zone model
Based on a fit to atomistic calculations, the specific fracture energy φ on inter-surfaces between bulks in the exponential cohesive zone model is given by [7]:where δn and δt represent normal and shear characteristic lengths related to the debonding process, respectively and r is defined as:where is the value of Δn after complete shear separation with Tn = 0.
To yield the realistic results of mixed-mode condition, it is necessary
The development of an improved interfacial bonding model
A potential function usually has a minimum at the equilibrium position in that the interaction force (derivative of the potential) must be attractive when distance is larger than equilibrium position and repulsive when distance is smaller than equilibrium position. To overcome the limitations of the original exponential cohesive zone model [7], the following interface potential function is proposed with four additional variables :
Path dependence of work-of-separation
The amount of energy dissipated strongly depends on separation paths when the normal and shear fracture energy in the interfacial zone are different. To evaluate the energy variation with respect to separation path, it is necessary to study the work-of-separation [20]. Many experimental data illustrates that the work-of-separation in normal direction is not equal to the work in tangential direction, e.g. Araki et al. [14], Benzeggagh and Kenane [15], Dollhofer et al. [16], Pang [17], Warrior et
Analytical analysis
The proposed interfacial bonding model was verified through analytical solution when the interface is in linear elastic deformation. The analytical solution of stress-strain relation is derived based on a simple compression test in a standard hexagonal structure, as shown in Fig. 8. The hexagonal plates will slide each other along the interface during compressive loading process. The red springs between bulks in Fig. 8(a) represent the interface traction.
From the deformed configuration in
Numerical simulation of extrafibrillar matrix in bone
Studying the fracture mechanism in bone is attracting the attention of engineering researchers, due to its highly hierarchical structure and the exceptional mechanical and load-bearing properties, e.g. Rho et al. [36], Weiner and Wagner [37] and Weiner et al. [38]. Bone poses various levels of hierarchical structural organization from macroscale to nanoscale [36], [39]. The mechanism of bone failure has been extensively studied, however, the underlying mechanism of plastic deformation in bone
Conclusions
In this work, we have reported an improved interfacial bonding model. The proposed model was examined in mixed-mode loading and the physically realistic coupling behavior was obtained. The work-of-separation analysis implies that the proposed interfacial bonding model can capture the work-of-separation in mixed mode, which is a realistic process for work transformation from normal to shear or from shear to normal. The proposed improved interfacial bonding model has the following properties:
- (1)
By
Acknowledgements
This work is partially supported by a grant from National Institutes of Health (Grant No. R21AR066925), a grant from National Science Foundation (Grant No. CMMI-1538448), and a grant from the University of Texas at San Antonio, Office of the Vice President for Research. We are grateful to Dr. Hai-Chao Han for valuable discussions in the derivation of the analytical solution.
References (61)
- et al.
Toughness amplification in natural composites
J Mech Phys Solids
(2011) - et al.
A bridging domain method for coupling continua with molecular dynamics
Comput Methods Appl Mech Eng
(2004) - et al.
Numerical simulation of crack growth in an isotropic solid with randomized internal cohesive bonds
J Mech Phys Solids
(1998) - et al.
The relation between crack growth resistance and fracture process parameters in elastic-plastic solids
J Mech Phys Solids
(1992) - et al.
Numerical simulations of fast crack growth in brittle solids
J Mech Phys Solids
(1994) - et al.
Computational modelling of impact damage in brittle materials
Int J Solids Struct
(1996) - et al.
Impact-induced delamination of composites: a 2D simulation
Compos Part B: Eng
(1998) - et al.
A unified potential-based cohesive model of mixed-mode fracture
J Mech Phys Solids
(2009) - et al.
Cohesive element modeling of viscoelastic fracture: application to peel testing of polymers
Int J Solids Struct
(2000) - et al.
Computational analysis of thin coating layer failure using a cohesive model and gradient plasticity
Eng Fract Mech
(2003)
Fracture toughness for mixed mode I/II of epoxy resin
Acta Mater
Measurement of mixed-mode delamination fracture toughness of unidirectional glass/epoxy composites with mixed-mode bending apparatus
Compos Sci Technol
Fracture mechanical characterisation of mixed-mode toughness of thermoplast/glass interfaces
Comput Mater Sci
Mixed mode fracture analysis and toughness of adhesive joints
Eng Fract Mech
Numerical simulations of adhesively-bonded beams failing with extensive plastic deformation
J Mech Phys Solids
An improved description of the exponential Xu and Needleman cohesive zone law for mixed-mode decohesion
Eng Fract Mech
Cohesive zone modeling of dynamic failure in homogeneous and functionally graded materials
Int J Plast
A multiscale cohesive zone model and simulations of fractures
Comput Methods Appl Mech Eng
Variability and anisotropy of mechanical behavior of cortical bone in tension and compression
J Mech Behav Biomed Mater
Principles of determination and verification of muscle forces in the human musculoskeletal system: muscle forces to minimise bending stress
J Biomech
Mechanisms governing the inelastic deformation of cortical bone and application to trabecular bone
Acta Biomater
Differences in the mechanical behavior of cortical bone between compression and tension when subjected to progressive loading
J Mech Behav Biomed Mater
An irreversible cohesive zone model for interface fatigue crack growth simulation
Eng Fract Mech
Cohesive zone representations of failure between elastic or rigid solids and ductile solids
Eng Fract Mech
Implicit-explicit finite elements in nonlinear transient analysis
Comput Methods Appl Mech Eng
Failure of mineralized collagen fibrils: modeling the role of collagen cross-linking
J Biomech
Effect of mineral–collagen interfacial behavior on the microdamage progression in bone using a probabilistic cohesive finite element model
J Mech Behav Biomed Mater
Mechanical properties and the hierarchical structure of bone
Med Eng Phys
Lamellar bone: structure–function relations
J Struct Biol
Mechanical properties of nanostructure of biological materials
J Mech Phys Solids
Cited by (13)
Research on trapezoidal shape optimization of laser remanufacturing heterogeneous interface under fatigue load
2022, International Journal of FatigueCitation Excerpt :The interface fatigue crack propagation also showed a strong orientation preference [24]. In the initial stage of the homogeneous material propagation, the mode I direction was generally selected [25]. However, the interface crack propagation was always under the action of two to three stress intensity factors.
Adherend surface roughness effect on the mechanical response of silicone-based adhesive joints
2020, Engineering Fracture MechanicsCitation Excerpt :In this area, there are several techniques to attach subparts to each other, including adhesives, especially for bonding dissimilar materials. Adhesives are widely used materials in both engineering and biology industries [1–4] and the enhancement of interfacial properties is one of the key research topics [5–7]. The measurement of the interfacial mechanical properties between adhesive and adherend, which is essential to assess an optimized superior interface resistance, can be performed through peeling tests.
A nonlocal method for modeling interfaces: Numerical simulation of decohesion and sliding at grain boundaries
2020, Computer Methods in Applied Mechanics and EngineeringA heterogeneous orientation criterion for crack modelling in cortical bone using a phantom-node approach
2018, Finite Elements in Analysis and DesignCitation Excerpt :Their results claim that low-stiffness osteons (newly formed) may toughen cortical bone tissue as microcracks tend to propagate towards them, limiting their growth [35]. Some authors have developed cohesive zone models to simulate material interface behaviour at different length scales in bone [9,36–39]. For example, Lin et al. [36] recently defined a cohesive zone model to define the mechanical behaviour at the nanoscale of the extrafibrillar matrix in bone and the interfacial interactions in a simplified hexagonal-based model under compression loading [36].
Calculation of the critical energy release rate G<inf>c</inf> of the cement line in cortical bone combining experimental tests and finite element models
2017, Engineering Fracture MechanicsCitation Excerpt :Other authors have simulated material interface behavior at different scales in bone by means of cohesive zone models [4,15–20]. Lin et al. [18] recently defined a cohesive zone model to define the mechanical behavior at the nano scale of the extrafibrillar matrix in bone. Pereira et al. [19] studied the influence of hydration on the fracture properties of bovine cortical bone at macroscopic scale by carrying out DCB tests.
Computational modeling and simulation of epithelial wound closure
2023, Scientific Reports