Enriched contact finite elements for stable peeling computations

During peeling of a soft elastic strip from a substrate, strong adhesional forces act locally inside the peeling zone. It is shown here that when a standard contact finite element (FE) formulation is used to compute the peeling process, a large mesh refinement is required since the numerical solution procedure becomes unstable otherwise. To improve this situation, several different efficient enrichment strategies are presented that provide stable solution algorithms for comparably coarse meshes. The enrichment is based on the introduction of additional unknowns inside the contact elements discretizing the slave surface. These are chosen in order to improve the approximation of the peeling forces, while keeping the overall number of degrees of freedom low. If needed, these additional unknowns can be condensed out locally. The enrichment formulation is developed for both 2D and 3D nonlinear FE formulations. The new enrichment technique is applied to the peeling computation of a gecko spatula. The proposed enriched contact element formulations are also investigated in sliding computations. Copyright © 2011 John Wiley & Sons, Ltd.

[1]  N. Rizzo,et al.  Characterization of the structure and composition of gecko adhesive setae , 2006, Journal of The Royal Society Interface.

[2]  Ted Belytschko,et al.  Arbitrary discontinuities in finite elements , 2001 .

[3]  Ralph Spolenak,et al.  Evidence for capillarity contributions to gecko adhesion from single spatula nanomechanical measurements. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[4]  P. Wriggers.  Computational contact mechanics , 2012 .

[5]  Roger A. Sauer,et al.  An atomic interaction-based continuum model for adhesive contact mechanics , 2007 .

[6]  Giorgio Zavarise,et al.  Modeling of mixed-mode debonding in the peel test applied to superficial reinforcements , 2008 .

[7]  Tod A. Laursen,et al.  Two dimensional mortar contact methods for large deformation frictional sliding , 2005 .

[8]  Giulio Alfano,et al.  Adaptive hierarchical enrichment for delamination fracture using a decohesive zone model , 2002 .

[9]  Peter Wriggers,et al.  Mortar based frictional contact formulation for higher order interpolations using the moving friction cone , 2006 .

[10]  Yu Tian,et al.  Adhesion and friction in gecko toe attachment and detachment , 2006, Proceedings of the National Academy of Sciences.

[11]  Dominique Chamoret,et al.  New smoothing procedures in contact mechanics , 2004 .

[12]  Shaker A. Meguid,et al.  On the modelling of smooth contact surfaces using cubic splines , 2001 .

[13]  Long Chen.  FINITE ELEMENT METHOD , 2013 .

[14]  Roger A. Sauer,et al.  The Peeling Behavior of Thin Films with Finite Bending Stiffness and the Implications on Gecko Adhesion , 2011 .

[15]  Roger A. Sauer,et al.  Computational Contact Formulations for Soft Body Adhesion , 2012 .

[16]  Roger A. Sauer,et al.  A contact mechanics model for quasi‐continua , 2007 .

[17]  Xiaopeng Xu,et al.  Numerical simulations of fast crack growth in brittle solids , 1994 .

[18]  Tod A. Laursen,et al.  A segment-to-segment mortar contact method for quadratic elements and large deformations , 2008 .

[19]  M. Puso,et al.  A 3D contact smoothing method using Gregory patches , 2002 .

[20]  Roger A. Sauer,et al.  A detailed 3D finite element analysis of the peeling behaviour of a gecko spatula , 2013, Computer methods in biomechanics and biomedical engineering.

[21]  P. Wriggers,et al.  Smooth C1‐interpolations for two‐dimensional frictional contact problems , 2001 .

[22]  Peter Wriggers,et al.  Formulation and analysis of a three-dimensional finite element implementation for adhesive contact at the nanoscale , 2009 .

[23]  T. Laursen.  Computational Contact and Impact Mechanics: Fundamentals of Modeling Interfacial Phenomena in Nonlinear Finite Element Analysis , 2002 .

[24]  Roger A Sauer,et al.  Multiscale modelling and simulation of the deformation and adhesion of a single gecko seta , 2009, Computer methods in biomechanics and biomedical engineering.

[25]  J. C. Simo,et al.  Geometrically non‐linear enhanced strain mixed methods and the method of incompatible modes , 1992 .

[26]  P. Wriggers,et al.  A C1-continuous formulation for 3D finite deformation frictional contact , 2002 .

[27]  M Mohammad Samimi,et al.  An enriched cohesive zone model for delamination in brittle interfaces , 2009 .

[28]  D. Braess.  BOOK REVIEW: Finite Elements: Theory, fast solvers and applications in solid mechanics, 2nd edn , 2002 .

[29]  A. Klarbring,et al.  Rigid contact modelled by CAD surface , 1990 .

[30]  G.A.O. Davies,et al.  Decohesion finite element with enriched basis functions for delamination , 2009 .

[31]  Alexander Konyukhov,et al.  Incorporation of contact for high-order finite elements in covariant form , 2009 .

[32]  Pavel Neuzil,et al.  The nature of the gecko lizard adhesive force. , 2005, Biophysical journal.

[33]  Ted Diehl,et al.  On using a penalty-based cohesive-zone finite element approach, Part I: Elastic solution benchmarks , 2008 .

[34]  Ernst Rank,et al.  A comparison of the h-, p-, hp-, and rp-version of the FEM for the solution of the 2D Hertzian contact problem , 2010 .

[35]  John W. Hutchinson,et al.  Interface strength, work of adhesion and plasticity in the peel test , 1998 .

[36]  Roger A. Sauer,et al.  A computational contact model for nanoscale rubber adhesion , 2009 .

[37]  T. Laursen,et al.  A framework for development of surface smoothing procedures in large deformation frictional contact analysis , 2001 .

[38]  E. Rank,et al.  Erratum to: A comparison of the h-, p-, hp-, and rp-version of the FEM for the solution of the 2D Hertzian contact problem , 2010 .

[39]  Gerhard A. Holzapfel,et al.  Cn continuous modelling of smooth contact surfaces using NURBS and application to 2D problems , 2003 .