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.