A l'échelle de plusieurs dizaines de millions d'années, un système de subduction implique de grandes déformations de la plaque plongeante assimilée un solide viscoélastique, et du manteau supérieur assimilé à un fluide newtonien. L'objectif de ce travail est de développer une stratégie de couplage solide-fluide appliquée à l'étude de l'interaction lithosphère-asthénosphère. Cette stratégie est basée sur l'utilisation de maillages non-conformes aux interfaces et d'une méthode de domaines fictifs (MDF) pour la résolution du problème fluide. Pour l'efficience des modèles 3D, nous employons une formulation simplifiée de la méthode de domaines fictifs par multiplicateurs de Lagrange. La MDF développée est validée par des comparaisons avec des solutions analytiques qui montrent que la méthode est d'ordre 1. La stratégie de couplage est également validée par la comparaison avec d'autres méthodes de couplage solide-fluide. Une première étude est ensuite menée pour analyser l'influence de certains paramètres rhéologiques et cinématiques sur la dynamique d'une subduction contrôlée par les vitesses des plaques. Cette étude, en 2D, concerne plus spécifiquement le mécanisme de plissement périodique du slab lorsque celui-ci est ancré à 660 km de profondeur. Ce mécanisme induit des variations de pendage du slab générant des variations de l'état de contrainte de la plaque chevauchante. Un intérêt particulier est porté sur l'influence de la viscosité du manteau sur les plissements. Dans ce cadre, nous réalisons une application à la subduction andine. / Over the time scale of tens of millions of years, a subduction system involves large deformations of tectonics plates, as one plate sinks into the Earth's mantle. The aim of this work was to develop a soli-fluid coupling method applied to the lithosphere-asthenosphere interaction in the context of subduction zones. Plates were assumed to behave as viscoelastic bodies, while the upper mantle was assimilated to a newtonian fluid. The method developped here is based on the use of non-matching interface meshes and a fictitious domain method (FDM) for the fluid problem. To optimize the computational efficiency of 3D model, we used a simplified version of the Lagrange multipliers fictitious domain method. The developped FDM has been benchmarked with analytical solutions and we showed that this FDM is a first-order method. The coupling method has also been compared to other fluid-solid coupling methods using matching interfaces meshes. A first two-dimensional study was performed in order to evaluate the influence of some rheological and kinematic parameters on the dynamics of a subduction controlled by the velocity of the plates. This study aimed at investigating cyclic slab folding over a rigid 660 km depth transition zone. This folding mechanism induces variations in slab dip that generate variations in the stress state of the overriding plate. We focussed on the influence of the upper mantle viscosity on slab folding. We also applied this model to the Andean subduction zone. Several studies have determined a cyclic variation of the South-American tectonic regime (period of 30-40~Myrs) which may have been related to the slab dip evolution.