Finite element model for linear-elastic mixed mode loading using adaptive mesh strategy

An adaptive mesh finite element model has been developed to predict the crack propagation direction as well as to calculate the stress intensity factors (SIFs), under linear-elastic assumption for mixed mode loading application. The finite element mesh is generated using the advancing front method. In order to suit the requirements of the fracture analysis, the generation of the background mesh and the construction of singular elements have been added to the developed program. The adaptive remeshing process is carried out based on the posteriori stress error norm scheme to obtain an optimal mesh. Previous works of the authors have proposed techniques for adaptive mesh generation of 2D cracked models. Facilitated by the singular elements, the displacement extrapolation technique is employed to calculate the SIF. The fracture is modeled by the splitting node approach and the trajectory follows the successive linear extensions of each crack increment. The SIFs values for two different case studies were estimated and validated by direct comparisons with other researchers work.

[1]  Michael H. Santare,et al.  Numerical Calculation of Stress Intensity Factors in Functionally Graded Materials , 2000 .

[2]  Paul A. Wawrzynek,et al.  Quasi-automatic simulation of crack propagation for 2D LEFM problems , 1996 .

[3]  A. B. de Morais,et al.  Calculation of stress intensity factors by the force method , 2007 .

[4]  Jaime Planas,et al.  KI evaluation by the displacement extrapolation technique , 2000 .

[5]  Ahmad Kamal Ariffin,et al.  An adaptive finite element procedure for crack propagation analysis , 2007 .

[6]  Ahmad Kamal Ariffin,et al.  Finite element simulation of stress intensity factors in elastic-plastic crack growth , 2006 .

[7]  D. Munz,et al.  Fracture tests for ceramics under mode-I, mode-II and mixed-mode loading , 1995 .

[8]  Y. Mutoh,et al.  A general mixed-mode brittle fracture criterion for cracked materials , 2006 .

[9]  Rainald Löhner,et al.  Automatic unstructured grid generators , 1997 .

[10]  Luiz Fernando Martha,et al.  NUMERICAL ESTIMATION OF FRACTURE PARAMETERS IN ELASTIC AND ELASTIC-PLASTIC ANALYSIS , 2000 .

[11]  P.M.S.T. de Castro,et al.  Dual boundary element analysis of cracked plates: post-processing implementation of the singularity subtraction technique , 2004 .

[12]  Stéphane Bordas,et al.  Enriched finite elements and level sets for damage tolerance assessment of complex structures , 2006 .

[13]  Louis Napoleon George Filon,et al.  On an Approximate Solution for the Bending of a Beam of Rectangular Cross-Section under any System of Load, with Special Reference to Points of Concentrated or Discontinuous Loading , 1903 .

[14]  Sutthisak Phongthanapanich,et al.  Adaptive Delaunay triangulation with object-oriented programming for crack propagation analysis , 2004 .

[15]  André P. Vieira,et al.  Subcritical crack growth: the microscopic origin of Paris' law. , 2008, Physical review letters.

[16]  R. Barsoum On the use of isoparametric finite elements in linear fracture mechanics , 1976 .

[17]  Ted Belytschko,et al.  Level set crack propagation modelling in the element free galerkin method. , 2001 .

[18]  D. M. Tracey,et al.  The natural isoparametric triangle versus collapsed quadrilateral for elastic crack analysis , 1976, International Journal of Fracture.

[19]  Giorgio Monti,et al.  IV European Conference on Computational Mechanics , 2010 .

[20]  Kimmo Berg,et al.  European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS) , 2004 .

[21]  O. C. Zienkiewicz,et al.  The Finite Element Method: Its Basis and Fundamentals , 2005 .

[22]  M. Boniecki,et al.  Subcritical crack growth in zirconia-toughened alumina (ZTA) ceramics , 2006 .