Mode I fracture in adhesively-bonded joints: A mesh-size independent modelling approach using cohesive elements

•A Cohesive Zone Model has been used to simulate mode I fracture in adhesive joints.•Mesh-size independent numerical solutions are obtained with quadratic FEA models.•Quadratic cohesive elements accurately predict the energy dissipated in the FPZ.•The proposed quadratic element formulation results i...

Full description

Saved in:
Bibliographic Details
Published in:Engineering fracture mechanics 2014-01, Vol.115, p.73-95
Main Authors: Álvarez, D., Blackman, B.R.K., Guild, F.J., Kinloch, A.J.
Format: Article
Language:English
Subjects:
Adhesive joints
Adhesives
Cohesion
Cohesive Zone Modelling
Computer simulation
Crack propagation
Finite element method
Finite-element analysis
Fracture
Fracture mechanics
Industrial applications
Mathematical models
Mesh-dependency/mesh-size
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:•A Cohesive Zone Model has been used to simulate mode I fracture in adhesive joints.•Mesh-size independent numerical solutions are obtained with quadratic FEA models.•Quadratic cohesive elements accurately predict the energy dissipated in the FPZ.•The proposed quadratic element formulation results in shorter simulation times. In recent years cohesive elements, coupled with a finite-element analysis (FEA) approach, have become increasingly popular for simulating both delamination in composite materials and fracture in adhesively-bonded joints. However, the industrial application of Cohesive Zone Models to model large and complex structures has been hindered by the requirement of extremely fine meshes along the crack propagation path. In the present work two-dimensional linear and quadratic (i.e. second-order) cohesive elements to model crack initiation and growth have been implemented in Abaqus using a user subroutine. These elements, which have a modified topology that allows a user-defined number of integration points, have been employed to model the fracture response of various mode I test specimens consisting of metallic substrates bonded with a structural film-adhesive. The effects of the mesh-density, element order and number of integration points on the numerical solution have been investigated. Whilst the linear models have shown the typical mesh-size dependent behaviour, the results obtained with their quadratic counterparts have been found to be independent of the element size. Furthermore, it is shown that increasing the number of integration points improves the stability, convergence and smoothness of the solutions. The mesh-size independent response obtained with the quadratic models arises from more accurate simulation of the deformed profile of the substrates and a more accurate calculation of the energy dissipated in the process zone due to damage. Overall, it is demonstrated that the quadratic cohesive-element formulation enables the use of much coarser meshes, resulting in shorter simulation times, and will therefore allow an increase in the industrial application of Cohesive Zone Models.
ISSN:0013-7944
1873-7315
DOI:10.1016/j.engfracmech.2013.10.005