Le contexte d'étude est celui de la Compatibilité ÉlectroMagnétique (CEM). L'objectif de la CEM est, comme son nom l'indique, d'assurer la compatibilité entre une source de perturbation électromagnétique et un système électronique victime. Or, la prédiction de ces niveaux de perturbation ne peut pas s'effectuer à l'aide d'un simple calcul analytique, en raison de la géométrie qui est généralement complexe pour le système que l'on étudie, tel que le champ à l'intérieur d'un cockpit d'avion par exemple. En conséquence, nous sommes contraints d'employer des méthodes numériques, dans le but de prédire ce niveau de couplage entre les sources et les victimes. Parmi les nombreuses méthodes numériques existantes à ce jour, les méthodes Multi-Domaines (MD) sont très prisées. En effet, elles offrent la liberté aux utilisateurs de choisir la méthode numérique la plus adaptée, en fonction de la zone géométrique à calculer. Au sein de ces méthodes MD, la " Domain Decomposition Method " (DDM) présente l'avantage supplémentaire de découpler chacun de ces domaines. En conséquence, la DDM est particulièrement intéressante, vis-à-vis des méthodes concurrentes, en particulier sur l'aspect du coût numérique. Pour preuve, l'ONERA continue de développer cette méthode qui ne cesse de montrer son efficacité depuis plusieurs années, notamment pour le domaine des Surfaces Équivalentes Radar (SER) et des antennes. L'objectif de l'étude est de tirer profit des avantages de cette méthode pour des problématiques de CEM. Jusqu'à maintenant, de nombreuses applications de CEM, traitées par le code DDM, fournissaient des résultats fortement bruités. Même pour des problématiques électromagnétiques très simples, des problèmes subsistaient, sans explication convaincante. Ceci justifie cette étude. Le but de cette thèse est de pouvoir appliquer ce formalisme DDM à des problématiques de CEM. Dans cette optique, nous avons été amenés à redéfinir un certain nombre de conventions, qui interviennent au sein de la DDM. Par ailleurs, nous avons développé un modèle spécifique pour les ouvertures, qui sont des voies de couplage privilégiées par les ondes, à l'intérieur des cavités que représentent les blindages. Comme les ouvertures sont, en pratique, de petites dimensions devant la longueur d'onde, on s'est intéressé à un modèle quasi-statique. Nous proposons alors un modèle, qui a été implémenté, puis validé. Suite à ce modèle, nous avons développé une méthode originale, basée sur un calcul en deux étapes, permettant de ne plus discrétiser le support des ouvertures dans les calculs 3D.