Un modèle numérique basé sur des techniques de dynamique moléculaire a été développé pour étudier le comportement mécanique d'ensembles de fibres enchevêtrées. Pour limiter le nombre de degrés de liberté, les fibres sont discrétisées en une suite de segments. Le modèle permet de prendre en compte des paramètres géométriques (fraction volumique, facteur de forme, distribution d'orientation...) et des paramètres « matériaux » (flexibilité, frottement aux contacts...). Les comportements macroscopiques d'ensembles de fibres orientées aléatoirement ont été identifiés en fonction de ces paramètres, enrichissant ainsi les lois d'échelles simples de la littérature. Des informations microscopiques telles que l'évolution du nombre de contacts en fonction des sollicitations (compression isostatique, cisaillement) ont également pu être obtenues numériquement. Ceci a notamment permis de valider des techniques d'analyse d'image, ensuite appliquées à des images de tomographie aux rayons X de laines d'acier.