Energetic versus maximally-dissipative local solutions of a quasi-static rate-independent mixed-mode delamination model

A quasi-static rate-independent model of delamination of linearly elastic bodies at small strains, sensitive to mode of delamination, using interfacial damage and interfacial plasticity as two internal parameters, is further developed with the aim to extract representations typically employed in engineering interface-models, i.e. fracture envelope and fracture energy dependence on the mode mixity, which are suitable for the model fitting to experimental data. Moreover, two concepts of solutions are implemented: globally stable energy-conserving solutions or stress-driven maximally-dissipative local solutions, arising by the fully implicit or by a semi-implicit time-stepping procedures, respectively, both yielding numerically stable and convergent time-discretizations. Spatial discretization is performed by the symmetric Galerkin boundary-element method (SGBEM). Alternating quadratic programming is implemented to cope with, respectively, global or local, energy-minimizations in the computation of the time-discretized solutions. Sample 2D numerical examples document applicability of the model as well as efficiency of the SGBEM numerical implementation and facilitate comparison of the two mentioned solution concepts.