Numerical homogenization of concrete microstructures without explicit meshes
Life management of electric hydro or nuclear power plants requires to estimate long-term concrete properties on facilities, for obvious safety and serviceability reasons. Decades-old structures are foreseen to be operational for several more decades. As a large number of different concrete formulations are found in EDF facilities, empirical models based on many experiments cannot be an option for a large fleet of power plant buildings. To build predictive models, homogenization techniques offer an appealing alternative. To properly upscale creep, especially at long term, a rather precise description of the microstructure is required. However, the complexity of the morphology of concrete poses several challenges. In particular, concrete is formulated to maximize the packing density of the granular skeleton, leading to aggregates spanning several length scales with small inter particle spacings. Thus, explicit meshing of realistic concrete microstructures is either out of reach of current meshing algorithms or would produce a number of degrees of freedom far higher than the current generic FEM codes capabilities. This paper proposes a method to deal with complex matrix-inclusions microstructures such as the ones encountered at the mortar or concrete scales, without explicitly meshing them. The microstructure is superimposed to an independent mesh, which is a regular Cartesian grid. This inevitably yields so called 'gray elements', spanning across multiple phases. As the reliability of the estimate of the effective properties highly depends on the behavior affected to these gray elements, special attention is paid to them. As far as the question of the solvers is concerned, generic FEM codes are found to lack efficiency: they cannot reach high enough levels of discretization with classical free meshes, and they do not take advantage of the regular structure of the mesh. Thus, a specific finite differences/finite volumes solver has been developed. At first, generic off-the-shelf linear system solvers were used. To further improve the efficiency in terms of memory requirements, specific variants of the preconditioned conjugate gradient were implemented. This allowed to homogenize the conductivity of a concrete-like microstructure using more than 10{sup 9} degrees of freedom on a rather common hardware for 2010 (a PC embedding 48 GB of RAM). Taking benefit of the properties of the regular Cartesian grid we have also investigated a multi-level method to improve the CPU efficiency of the code.
- OSTI ID:
- 21576907
- Journal Information:
- Cement and Concrete Research, Vol. 41, Issue 12; Conference: International summit on cement hydration kinetics and modeling;CONMOD10: Symposium on concrete modelling, Quebec City, PB (Canada);Lausanne (Switzerland), 27-29 Jul 2009;22-25 Jun 2010; Other Information: DOI: 10.1016/j.cemconres.2011.03.023; PII: S0008-8846(11)00099-8; Copyright (c) 2011 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; ISSN 0008-8846
- Country of Publication:
- United States
- Language:
- English
Similar Records
Seismic velocity estimation from time migration
Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories
Related Subjects
ALGORITHMS
CONCRETES
CREEP
DEGREES OF FREEDOM
EFFICIENCY
MICROSTRUCTURE
MORPHOLOGY
MORTARS
NUCLEAR POWER PLANTS
PARTICLES
RELIABILITY
SAFETY
SKELETON
BODY
BUILDING MATERIALS
MATERIALS
MATHEMATICAL LOGIC
MECHANICAL PROPERTIES
NUCLEAR FACILITIES
ORGANS
POWER PLANTS
THERMAL POWER PLANTS