A corotational hybrid-Trefftz stress formulation for modelling cohesive cracks

This paper presents a corotational hybrid-Trefftz formulation for modelling propagating cohesive cracks and captures moderate rotations but small strains. The formulation is characterised by the fact that stresses are approximated within the domain of the element and the stiffness can be expressed via a boundary rather than a domain integral. Thus, compared to their FEM counterpart, hybrid-Trefftz stress elements exhibit faster convergence of the stress fields. Furthermore, the displacements are approximated on element boundaries and the displacement basis is defined independently on each element interface. Thus, the overall bandwidth of the stiffness matrix is very small and computationally efficient to solve. A corotational formulation for hybrid-Trefftz elements is also introduced in order to capture the effect of geometric nonlinearities in the form of moderate rotations. The model’s performance is demonstrated on three examples, illustrating crack propagation and the influence of geometric nonlinearities.

[1]  C. Rankin,et al.  Finite rotation analysis and consistent linearization using projectors , 1991 .

[2]  Kim-Chuan Toh,et al.  Fast iterative solution of large undrained soil‐structure interaction problems , 2003 .

[3]  Carlos A. Felippa,et al.  A unified formulation of small-strain corotational finite elements: I. Theory , 2005 .

[4]  J. A. Freitas,et al.  Formulation of elastostatic hybrid-Trefftz stress elements , 1998 .

[5]  A. P. Zieliński,et al.  On trial functions applied in the generalized Trefftz method , 1995 .

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

[7]  C. G. Broyden,et al.  The convergence of an algorithm for solving sparse nonlinear systems , 1971 .

[8]  Eros Roberto Grau.  TEIXEIRA DE FREITAS , 2002 .

[9]  A. Hillerborg,et al.  Analysis of crack formation and crack growth in concrete by means of fracture mechanics and finite elements , 1976 .

[10]  James Demmel,et al.  An Asynchronous Parallel Supernodal Algorithm for Sparse Gaussian Elimination , 1997, SIAM J. Matrix Anal. Appl..

[11]  de R René Borst,et al.  On the numerical integration of interface elements , 1993 .

[12]  J. Oliver.  MODELLING STRONG DISCONTINUITIES IN SOLID MECHANICS VIA STRAIN SOFTENING CONSTITUTIVE EQUATIONS. PART 2: NUMERICAL SIMULATION , 1996 .

[13]  M. V. Vliet.  Size effect in Tensile Fracture of Concrete and Rock , 2000 .

[14]  Garth N. Wells,et al.  Discontinuous modelling of strain localisation and failure , 2001 .

[15]  J. Oliver.  MODELLING STRONG DISCONTINUITIES IN SOLID MECHANICS VIA STRAIN SOFTENING CONSTITUTIVE EQUATIONS. PART 1: FUNDAMENTALS , 1996 .

[16]  Martinus Gertrudis Auntonius Tijssens.  On the cohesive surface methodology for fracture of brittle heterogeneous solids : computational and material modeling : proefschrift , 2000 .