Multiscale computational homogenization methods with a gradient enhanced scheme based on the discontinuous Galerkin formulation

in Computer Methods in Applied Mechanics & Engineering (2013), 260

When considering problems of dimensions close to the characteristic length of the material, the size e ects can not be neglected and the classical (so–called first–order) multiscale computational ... [more ▼]

When considering problems of dimensions close to the characteristic length of the material, the size e ects can not be neglected and the classical (so–called first–order) multiscale computational homogenization scheme (FMCH) looses accuracy, motivating the use of a second–order multiscale computational homogenization (SMCH) scheme. This second–order scheme uses the classical continuum at the micro–scale while considering second–order continuum at the macro–scale. Although the theoretical background of the second–order continuum is increasing, the implementation into a finite element code is not straightforward because of the lack of high–order continuity of the shape functions. In this work, we propose a SMCH scheme relying on the discontinuous Galerkin (DG) method at the macro–scale, which simplifies the implementation of the method. Indeed, the DG method is a generalization of weak formulations allowing for inter-element discontinuities either at the C0 level or at the C1 level, and it can thus be used to constrain weakly the C1 continuity at the macro–scale. The C0 continuity can be either weakly constrained by using the DG method or strongly constrained by using usual C0 displacement–based finite elements. Therefore, two formulations can be used at the macro–scale: (i) the full–discontinuous Galerkin formulation (FDG) with weak C0 and C1 continuity enforcements, and (ii) the enriched discontinuous Galerkin formulation (EDG) with high–order term enrichment into the conventional C0 finite element framework. The micro–problem is formulated in terms of standard equilibrium and periodic boundary conditions. A parallel implementation in three dimensions for non–linear finite deformation problems is developed, showing that the proposed method can be integrated into conventional finite element codes in a straightforward and e cient way. [less ▲]

Imposing periodic boundary condition on arbitrary meshes by polynomial interpolation

in Computational Materials Science (2012), 55

In order to predict the effective properties of heterogeneous materials using the finite element approach, a boundary value problem (BVP) may be defined on a representative volume element (RVE) with ... [more ▼]

In order to predict the effective properties of heterogeneous materials using the finite element approach, a boundary value problem (BVP) may be defined on a representative volume element (RVE) with appropriate boundary conditions, among which periodic boundary condition is the most efficient in terms of convergence rate. The classical method to impose the periodic boundary condition requires the identical meshes on opposite RVE boundaries. This condition is not always easy to satisfy for arbitrary meshes. This work develops a new method based on polynomial interpolation that avoids the need of matching mesh condition on opposite RVE boundaries. [less ▲]

Projects in Fracture Simulations

Scientific conference (2011, September)