Les études présentées dans ce mémoire portent sur la simulation numérique et l'analyse physique des milieux granulaires immergés dans un fluide. Des développements numériques ont été réalisés pour coupler la méthode Lattice Boltzmann pour la dynamique du fluide avec la méthode Contact Dynamics en 2D et avec la méthode Molecular Dynamics en 3D pour la dynamique des grains. Ces outils numériques ont été utilisés pour étudier l'initiation des avalanches sur un plan incliné en fonction de la compacité initiale et de l'angle d'inclinaison en 3D. Les résultats sont en bon accord quantitatif avec les expériences et ont permis de mettre en évidence la stabilisation de la pente granulaire par une pression négative du fluide interstitielle induite par la dilatance, et l'évolution spatiotemporelle des grandeurs telles que la compacité et la déformation de cisaillement. Ces évolutions dans la phase de fluage qui précède la rupture de pente ont pu être mises à l'échelle par un modèle théorique incorporant la loi de Darcy et l'effet de la dilatance sur l'angle de frottement interne. L'analyse de la texture granulaire a révélé la distortion du réseau des contacts pendant le fluage et la saturation de l'anisotropie comme un critère de rupture. La propagation des avalanches granulaires a été étudiée dans une configuration 2D pour deux géométries différentes : 1) l'effondrement et l'étalement d'une colonne sous son propre poids, 2) l'étalement d'une pente sous l'effet d'une énergie cinétique injectée. Nous avons en particulier montré que la distance et la durée d'étalement obéissent à des lois de puissance en fonction du rapport d'aspect initial ou de l'énergie injectée. Le fluide exerce deux effets contradictoires : réduire les temps de relaxation et lubrifier les contacts. Ces effets ont été analysés dans le régime visqueux en fonction des conditions initiales et la viscosité du fluide.