MULTISCALE FINITE ELEMENT METHOD FOR HETEROGENEOUS MEDIA WITH MICROSTRUCTURES: CRACK PROPAGATION IN A POROUS MEDIUM

A novel multiscale FE scheme is proposed for simulation of crack propagation in heterogeneous media. A fine scale mesh is constructed to model the crack tip and the microstructures around the near-tip field, while a coarse scale mesh is introduced to model the far field, wherein the effect of the microstructures is averaged through the homogenization theory to yield the average macroscopic constitutive relationship. The so-called variable node elements are subsequently employed to connect the fine scale zone in a seamless way to the coarse scale zone. The moving crack tip is modeled in a straightforward manner with the conventional quarter-point singular element. Several numerical examples are presented to demonstrate the scheme as an effective tool for multiscale simulation of crack propagation in consideration of the microstructures in a heterogeneous body.

[1]  Jae Hyuk Lim,et al.  A new computational approach to contact mechanics using variable‐node finite elements , 2008 .

[2]  Laurent Champaney,et al.  A multiscale extended finite element method for crack propagation , 2008 .

[3]  Jae Hyuk Lim,et al.  Variable‐node elements for non‐matching meshes by means of MLS (moving least‐square) scheme , 2007 .

[4]  Ted Belytschko,et al.  A multiscale projection method for macro/microcrack simulations , 2007 .

[5]  Jae Hyuk Lim,et al.  MLS (moving least square)-based finite elements for three-dimensional nonmatching meshes and adaptive mesh refinement , 2007 .

[6]  Somnath Ghosh,et al.  Concurrent multi-level model for damage evolution in microstructurally debonding composites , 2007 .

[7]  S. Im,et al.  (4+n)-noded Moving Least Square(MLS)-based finite elements for mesh gradation , 2007 .

[8]  M. H. Aliabadi,et al.  Multi-scale boundary element modelling of material degradation and fracture , 2007 .

[9]  Wing Kam Liu,et al.  Multiresolution analysis for material design , 2006 .

[10]  S. Im,et al.  MLS‐based variable‐node elements compatible with quadratic interpolation. Part II: application for finite crack element , 2006 .

[11]  S. Im,et al.  On the computation of the near-tip stress intensities for three-dimensional wedges via two-state M-integral , 2003 .

[12]  T. Belytschko,et al.  Extended finite element method for cohesive crack growth , 2002 .

[13]  Insu Jeon,et al.  The role of higher order eigenfields in elastic–plastic cracks , 2001 .

[14]  Hyun Gyu Kim,et al.  Mode decomposition of three-dimensional mixed-mode cracks via two-state integrals , 2001 .

[15]  S. Im,et al.  An application of two-state M-integral for computing the intensity of the singular near-tip field for a generic wedge , 2000 .

[16]  Ted Belytschko,et al.  Elastic crack growth in finite elements with minimal remeshing , 1999 .