L'etude porte sur la modélisation du comportement mécanique de matériaux granulaires cohésifs sous faibles niveaux de chargement. Le sujet est abordée par des approches numérique (en éléments discrets), expérimentale et théorique. Un code de calcul tridimensionnel (nommé tapio-K) a été mis au point avec comme spécificité d'intégrer plusieurs lois de cohésion. La confrontation d'essais expérimentaux et numériques de cisaillement direct d'échantillons de billes de verre a permis la validation du code de calcul. En outre, cette confrontation a mis en évidence la nécessité de réaliser des essais faiblement contraints afin de juger de l'influence de la cohésion capillaire. Un nouveau dispositif a été réalisé à cette fin. Les mesures obtenues à l'aide de ce dispositif, les simulations numériques, ainsi qu'une analyse théorique ont permis de mettre en évidence l'importance de la distribution du liquide au sein d'un matériau granulaire humide (phase liquide discontinue) sur la résistance au cisaillement. Enfin, quelques cas tests d'application (compression simple, vidange d'un matériau granulaire humide, traficabilité d'un sol de surface) sont présentés pour montrer les possibilités du code de calcul.

Ce travail de thèse est consacré à la modélisation et à l'étude numérique de la texture et du comportement mécanique des matériaux granulaires à grande polydispersité. La polydispersité est considérée en termes d'étalement et de courbure de la distribution de tailles. Afin d'aboutir à des distributions de tailles représentatives, une méthode de génération ainsi qu'un modèle de courbe granulométrique sont présentés. La mise en place de grands échantillons denses par une méthode géométrique a permis de réaliser une étude systématique de l'influence de la polydispersité sur la texture. Le comportement mécanique est étudié par simulations numériques discrètes sur la base de la méthode de Dynamique des Contacts. Les échantillons sont successivement soumis à une relaxation uniaxiale et à un cisaillement simple avec des conditions aux limites périodiques afin d'atteindre respectivement l'état d'équilibre statique et l'état stationnaire après un cisaillement long. Nos résultats montrent que la distribution uniforme par fraction volumique conduit à la compacité la plus élevée. Dans l'état stationnaire, la résistance au cisaillement est indépendante de la polydispersité. Une analyse systématique de la texture permet de montrer que cette propriété résulte, d'une part, de la compensation mutuelle entre les anisotropies des directions de contact et des longueurs des vecteurs branches et, d'autre part, du fait que les forces fortes sont essentiellement transmise à travers les plus grosses particules. Enfin, en présence de forces d'adhésion, la cohésion de Coulomb augmente avec l'étalement granulométrique.