La prédiction du comportement à long terme des matériaux cimentaires est un enjeu majeur pour contribuer à l'étude de la durabilité des structures précontraintes. Le présent travail porte sur l'utilisation de la méthode de l'inclusion équivalente, approche d'homogénéisation multi-échelle simplifiée, pour la prédiction du fluage dans ces matériaux. Le fluage est modélisé par la viscoélasticité linéaire sans vieillissement. La méthode de l'inclusion équivalente permet de contourner certaines difficultés et limitations que présentent les approches classiques. Pour les matériaux cimentaires, fortement hétérogènes, les approches multiéchelles classiques sont ou bien numériquement lourdes et très complexes à mettre en œuvre, ou bien pas suffisamment détaillées pour prendre en compte les spécificités d'une microstructure. La méthode de l'inclusion équivalente présente un juste-milieu et permet de calculer des microstructures simplifiées de type matrice-inclusions et de fournir des estimations ou des bornes sur le comportement homogénéisé. Sous sa forme variationnelle, la méthode de l'inclusion équivalente n'a jusqu'alors été mise en œuvre que pour des inclusions de forme sphérique. Le présent travail propose d'étendre cette méthode à des inclusions de forme ellipsoïdale dont la variation de l'élancement permet de modéliser de nouveaux éléments asphériques tels que les fissures, les fibres et les cristaux de portlandite. Cette complexification de la géométrie a un impact sur le temps de calcul, qui est amplifié dans le cadre du fluage. Le second volet du travail porte alors sur l'extension de la méthode de l'inclusion équivalente à la viscoélasticité linéaire sans vieillissement par l'intermédiaire de la transformée de Laplace-Carson. Une méthodologie efficace (tant du point de vue de la précision que de celui du temps de calcul) est finalement proposée pour effectuer l'inversion numérique de cette transformée / The prediction of long-term behaviour of cementitious materials is a major concern which contributs to the study of the durability of prestressed structures. This work focuses on the use of the equivalent inclusion method, simplified multi-scale homogenization approach, for the prediction of creep in these materials. Creep is modelled by the non-ageing linear viscoelasticity. The equivalent inclusion method overcomes certain difficulties and limitations posed by conventional approaches. For cementitious materials (highly heterogeneous), conventional multi-scale approaches are, either digitally heavy and complex to implement, or not sufficiently detailed to take into account the specificities of a microstructure. The equivalent inclusion method presents a middle way and allows the calculation of simplified matrix-inclusion type microstructures and to provide estimates or bounds on the homogenized behaviour.Under its variational form, the equivalent inclusion method has, up to now, been implemented only for spherical inclusions. This work proposes to extend this method to ellipsoidal inclusions whose variation of slenderness allows the modelling of new aspheric elements such as cracks, fibers and portlandite crystals. Such enrichment of the geometry has an impact on the computation time, that is amplified in the context of creep. The second aspect of the work then applies to the extension of the equivalent inclusion method to the non-ageing linear viscoelasticity by means of the Laplace-Carson transform. An effective methodology (both from the viewpoint of precision and calculation time) is finally proposed to perform the numerical inversion of this transform

La prédiction du comportement à long terme des matériaux cimentaires est un enjeu majeur pour contribuer à l'étude de la durabilité des structures précontraintes. Le présent travail porte sur l'utilisation de la méthode de l'inclusion équivalente, approche d'homogénéisation multi-échelle simplifiée, pour la prédiction du fluage dans ces matériaux. Le fluage est modélisé par la viscoélasticité linéaire sans vieillissement. La méthode de l'inclusion équivalente permet de contourner certaines difficultés et limitations que présentent les approches classiques. Pour les matériaux cimentaires, fortement hétérogènes, les approches multiéchelles classiques sont ou bien numériquement lourdes et très complexes à mettre en œuvre, ou bien pas suffisamment détaillées pour prendre en compte les spécificités d'une microstructure. La méthode de l'inclusion équivalente présente un juste-milieu et permet de calculer des microstructures simplifiées de type matrice-inclusions et de fournir des estimations ou des bornes sur le comportement homogénéisé. Sous sa forme variationnelle, la méthode de l'inclusion équivalente n'a jusqu'alors été mise en œuvre que pour des inclusions de forme sphérique. Le présent travail propose d'étendre cette méthode à des inclusions de forme ellipsoïdale dont la variation de l'élancement permet de modéliser de nouveaux éléments asphériques tels que les fissures, les fibres et les cristaux de portlandite. Cette complexification de la géométrie a un impact sur le temps de calcul, qui est amplifié dans le cadre du fluage. Le second volet du travail porte alors sur l'extension de la méthode de l'inclusion équivalente à la viscoélasticité linéaire sans vieillissement par l'intermédiaire de la transformée de Laplace-Carson. Une méthodologie efficace (tant du point de vue de la précision que de celui du temps de calcul) est finalement proposée pour effectuer l'inversion numérique de cette transformée / The prediction of long-term behaviour of cementitious materials is a major concern which contributs to the study of the durability of prestressed structures. This work focuses on the use of the equivalent inclusion method, simplified multi-scale homogenization approach, for the prediction of creep in these materials. Creep is modelled by the non-ageing linear viscoelasticity. The equivalent inclusion method overcomes certain difficulties and limitations posed by conventional approaches. For cementitious materials (highly heterogeneous), conventional multi-scale approaches are, either digitally heavy and complex to implement, or not sufficiently detailed to take into account the specificities of a microstructure. The equivalent inclusion method presents a middle way and allows the calculation of simplified matrix-inclusion type microstructures and to provide estimates or bounds on the homogenized behaviour.Under its variational form, the equivalent inclusion method has, up to now, been implemented only for spherical inclusions. This work proposes to extend this method to ellipsoidal inclusions whose variation of slenderness allows the modelling of new aspheric elements such as cracks, fibers and portlandite crystals. Such enrichment of the geometry has an impact on the computation time, that is amplified in the context of creep. The second aspect of the work then applies to the extension of the equivalent inclusion method to the non-ageing linear viscoelasticity by means of the Laplace-Carson transform. An effective methodology (both from the viewpoint of precision and calculation time) is finally proposed to perform the numerical inversion of this transform

Les travaux effectués au cours de cette thèse portent principalement sur la construction de nouvelles bornes du comportement effectif des matériaux biphasés de type matrice-inclusions à comportement linéaire élastique. Dans un premier temps, afin d’encadrer le comportement effectif, nous présentons une nouvelle approche numérique, inspirée des travaux de Huet (J. Mech. Phys. Solids 1990 ; 38:813-41), qui repose sur le calcul des comportements apparents associés à des volumes élémentaires (VE) non-carrés construits à partir d'assemblages de cellules de Voronoï, chaque cellule contenant une inclusion entourée de matrice. De tels VE non-carrés permettent d'éviter l'application directe des CL sur les inclusions à l’origine d’une surestimation artificielle des comportements apparents. En utilisant les théorèmes énergétiques de l'élasticité linéaire et des procédures de moyennisation appropriées portant sur les comportements apparents, un nouvel encadrement du comportement effectif est obtenu. Son application au cas d'un composite biphasé, constitué d'une matrice isotrope et de fibres cylindriques parallèles et identiques distribuées aléatoirement dans le plan transverse, conduit à des bornes plus resserrées que celles obtenues par Huet. En nous appuyant sur cette nouvelle procédure numérique, nous avons ensuite réalisé une étude statistique des comportements apparents à l'aide de simulations de type Monté Carlo. Puis, à partir des tendances issues de cette étude statistique, nous avons proposé et mis en œuvre de nouveaux critères de tailles de VER. / This work is devoted to the derivation of improved bounds for the effective behavior of random linear elastic matrix-inclusions composites. In order to bounds their effective behavior, we present a new numerical approach, inspired by the works of Huet (J. Mech. Phys. Solids 1990 ; 38:813-41), which relies on the computation of the apparent behaviors associated to non square (or non cubic) volume elements (VEs) comprised of Voronoï cells assemblages, each cell being composed of a single inclusion surrounded by the matrix. Such non-square VEs forbid any direct application of boundary conditions to particles which is responsible for the artificial overestimation of the apparent behaviors observed for square VEs. By making used of the classical bounding theorems for linear elasticity and appropriate averaging procedures, new bounds are derived from ensemble averages of the apparent behavior associated with non square VEs. Their application to a two-phase composite composed of an isotropic matrix and aligned identical fibers randomly and isotropically distributed in the transverse plane leads to sharper bounds than those obtained by Huet. Then, by making use of this new numerical approach, a statistical study of the apparent behavior is carried out by means of Monte Carlo simulations. Subsequently, relying on the trends derived from this study, some proposals to define RVE criteria are presented.