La technique de renforcement des sols par inclusions linéaires rigides, c'est-à-dire capables de reprendre des efforts de flexion et cisaillement, connaît depuis quelques décennies un essor important bien que la simulation et le calcul de tels ouvrages reste un problème difficile, en raison de la très forte hétérogénéité du sol renforcé combinée au grand nombre d'inclusions mis en jeu dans de tels ouvrages. Une modélisation multi phasique a récemment été développée permettant d'aborder le dimensionnement de tels ouvrages de géotechnique avec des temps de calcul considérablement réduits par rapport à ceux d'une modélisation numérique classique, où les inclusions de renforcement doivent être discrétisées individuellement. On s'intéresse dans le cadre de ce travail à la mise en oeuvre de ce modèle dans le cadre du comportement élastoplastique des différents constituants, lorsque les effets de flexion et de cisaillement des inclusions sont pris en compte. Après avoir développé un certain nombre de solutions analytiques de référence, un outil numérique basé sur la méthode des éléments finis est mis au point, puis appliqué à la simulation du comportement de différentes structures renforcées, telles que les radiers de fondation sur groupes de pieux, ou la stabilisation de pentes instables par l'incorporation de telles inclusions.