Cette thèse est consacrée à la modélisation par transition d'échelles d'une large classe de composites particulaires fortement chargés tels que les propergols solides. L'approche (AM) repose en amont sur une schématisation géométrique et cinématique inspirée des travaux de Christoffersen (1983). L'objectif des présents travaux est double : prouver l'applicabilité de l'AM à la viscohyperélasticité (comportement de la matrice des élastomères chargés) et évaluer quantitativement ses performances. Pour traiter le 1er point, l'AM est appliquée à un composite aléatoire à matrice viscohyperélastique, généré numériquement. On montre le caractère direct de la résolution du problème de localisation-homogénéisation grâce à un algorithme opérant dans l'espace-temps réel. Les résultats obtenus sont qualitativement corrects. Concernant le 2ème point, les effets des hypothèses cinématiques propres à l'AM sont testés au travers de comparaisons entre résultats (locaux et globaux) AM et éléments finis (EF) sur des microstructures périodiques (simple et complexe) satisfaisant la schématisation géométrique et pour des comportements de phase hyperélastiques et viscohyperélastiques. Un certain nombre d'atouts et de points d'amélioration de l'AM sont ainsi dégagés. Enfin, un programme transversal de confrontation des estimations à des résultats expérimentaux et à des calculs EF sur un propergol réel est élaboré et mis en œuvre sur un composite énergétique : la butalite 400. Chaque étape –essais sur composite et constituants, caractérisation morphologique par tomographie, maillage EF automatique de la microstructure réelle, détermination des VER (plusieurs centaines de grains) relatifs aux deux méthodes (AM et EF)– est détaillée. De multiples confrontations entre les résultats expérimentaux et les premiers résultats numériques (AM et EF) permettent de préciser les perspectives à entreprendre pour terminer la validation.

Cette thèse traite de la modélisation multi-échelle de composites particulaires fortement chargés. La méthode d’estimation, qualifiée d’“Approche Morphologique” (A.M.), repose sur une double schématisation géométrique et cinématique du composite permettant de fournir la réponse aux deux échelles. Afin d’évaluer les capacités prédictives de l’A.M. en élasticité linéaire avec évolution de l’endommagement, l’A.M. est testée vis-à-vis de ses aptitudes à rendre compte des effets de taille et d’interaction de particules sur la chronologie de décohésion. Pour cela, différentes microstructures périodiques simples, aléatoires monomodales et bimodale générées numériquement sont considérées. Les résultats obtenus sont cohérents avec les données de la littérature : la décohésion des grosses particules précède celle des plus petites et est d’autant plus précoce que le taux de charges est important. Puis, l’objectif est de coupler deux non-linéarités traitées séparément dans deux versions antérieures de l’A.M : l’endommagement par décohésion charges/matrice et les grandes déformations. La formulation du problème de localisation-homogénéisation est reprise à la source de manière analytique. Le critère de nucléation de défauts est étendu en transformations finies. Le problème obtenu, fortement non-linéaire, est résolu numériquement via un algorithme de Newton-Raphson. Les étapes sous-jacentes à la résolution (calcul de la matrice tangente, codage en langage Python®) sont explicitées. Des évaluations progressives (matériaux sain et endommagé)permettent de valider la mise en oeuvre numérique. Les effets de taille et d’interaction sont alors restitués en transformations finies. / This study is devoted to multi-scale modeling of highly-filled particulate composites.This method, the “Morphological Approach” (M.A.), is based on a geometrical and kinematicalschematization which allows the access to both local fields and homogenized response. In order toevaluate the predictive capacities of the M.A. considering a linear elastic behavior for the constituentsand evolution of damage, analysis is performed regarding the ability of the M.A. to accountfor particle size and interaction effects on debonding chronology. For that purpose, simple periodic,random monomodal and bimodal microstructures are considered. The results are consistent withliterature data : debonding of large particles occurs before the one of smaller particles and thehigher the particle volume fraction, the sooner the debonding. Finally, the objective is to operatethe coupling of two non linearities which were separately studied in previous versions of the M.A. :debonding between particles and matrix, and finite strains. The whole analytical background of theapproach is reconsidered in order to define the localization-homogenization problem. The nucleationcriterion is extended to the finite strains context. The final problem, strongly non linear, is numericallysolved through a Newton-Raphson algorithm. The different solving steps (jacobian matrix,coding with Python®) are developed. Progressive evaluations (sound and damage materials) allowthe validation of numerical implementation. Then, size and interaction effects are reproduced infinite strains.