International audienceFor more than a decade, numerical methods for periodic elasticity, based on the fast Fourier transform, have been used successfully as alternatives to more conventional (fem, bem) numerical techniques for composites. These methods are based on the direct, point-wise, discretization of the Lippmann-Schwinger equation, and a subsequent truncation of underlying Fourier series required for the use of the fast Fourier transform. The basic FFT scheme is very attractive, because of its simplicity of implementation and use. However, it cannot handle pores or rigid inclusions, for which a specific (and significantly more involved) treatment is required. In the present paper, we propose a new FFT-based scheme which is as simple as the basic scheme, while remaining valid for infinite contrasts. Since we adopted an energy principle as an alternative to the Lippmann-Schwinger equation, our scheme is derived within a variational framework. As a by-product, it provides an energetically consistent rule for the homogenization of boundary voxels, a question which has been pending since the introduction of Fourier-based methods

To cite this version:Sébastien Brisard, Luc Dormieux. Combining Galerkin approximation techniques with the principle of Hashin and Shtrikman to derive a new FFT-based numerical method for the homogenization of composites. Computer Methods in Applied Mechanics and Engineering, Elsevier, 2012, 217-220, pp.197-212 AbstractWe report on the mathematical analysis of two different, FFT-based, numerical schemes for the homogenization of composite media within the framework of linear elasticity: the basic scheme of Suquet (1994, 1998), and the energy-based scheme of Brisard and Dormieux (2010). Casting these two schemes as Galerkin approximations of the same variational problem allows us to assert their well-posedness and convergence. More importantly, we extend in this work their domains of application, by relieving some stringent conditions on the reference material which were previously thought necessary. The origins of the flaws of each scheme are identified, and a third scheme is proposed, which seems to combine the strengths of the basic and energy-based schemes, while leaving out their weaknesses. Finally, a rule is proposed for handling heterogeneous pixels/voxels, a situation frequently met when images of real materials are used as input to these schemes.

No abstract

International audienceMorphological quantification of the complex structure of hierarchical geomaterials is of great relevance for Earth science and environmental engineering, among others. To date, methods that quantify the 3D morphology on length scales ranging from a few tens of nanometers to several hundred nanometers have had limited success. We demonstrate, for the first time, that it is possible to go beyond visualization and to extract quantitative morphological information from X-ray images in the aforementioned length scales. As examples, two different hierarchical geomaterials exhibiting complex porous structures ranging from nanometer to macroscopic scale are studied: a flocculated clay water suspension and two hydrated cement pastes. We show that from a single projection image it is possible to perform a direct computation of the ultra-small angle-scattering spectra. The predictions matched very well the experimental data obtained by the best ultra-small angle-scattering experimental setups as observed for the cement paste. In this context, we demonstrate that the structure of flocculated clay suspension exhibit two well-distinct regimes of aggregation, a dense mass fractal aggregation at short distance and a more open structure at large distance, which can be generated by a 3D reaction limited cluster-cluster aggregation process. For the first time, a high-resolution 3D image of fibrillar cement paste cluster was obtained from limited angle nanotomograph

No abstract