L'objectif principal de cette thèse est la caractérisation des dynamiques hautes fréquance en milieu côtier et en particulier leurs extrêmes, par l'intermédiaire de l'étude de série temporelle biogéochimiques à long terme enregistrées par des systèmes automatisés. Les bases de données hautes fréquences utilisées dans cette étude proviennent majoritairement du programme MAREL, qui a été mise en oeuvre par l'IFREMER. Des séries temporelles basses fréquences provenant des programmes de surveillance du littoral SOMLIT (CNRS, INSU) et SRN (Ifremer) sont mises à contribution pour appuyer l'importance des systèmes automatisés. La méthode EMD (Empirical Mode Decomposition) nous a servi de base dans de nombreuses analyses pour étudier ces séries temporelles. Nous avons aussi utilisé des méthodes plus classiques empruntées aux domaines de l'analyse numérique et de la turbulence. Cette études se décompose en 3 parties, et plusieurs annexes. Les matériels et méthodes sont présentés dans la première partie. Dans la seconde partie, la méthode EMD nous a permis de mettre en avant les fortes fluctuations contenues dans les blooms, ainsi que de mener des analyses spectrales grâce à un couplage avec la transformée de Hilbert. L'analyse en composante principale (ACP) a mis en avant les principaux forçages exercés sur la production primaire et les profils de température SOMLIT laissent supposer un impact de la stratification sur l'intensité des blooms. Dans la troisième partie, nous avons mené une étude comparative entre les données basses fréquences et hautes fréquences. Et deux méthodes de cross-corrélation (TDIC et co-spectre) nous ont permis de définir une échelle caractéristique de transition entre les températures de la Manche occidentale et orientale. En annexe, nous avons testé la robustesse de différentes méthodes d'analyses spectrales quant au manque de données dans les séries temporelles, qui est un problème inhérent aux bases de données enregistrées par des systèmes automatisés, et nous avons reproduit un article qui est en cours de soumission. / The main objective of this thesis is the characterization of high frequency dynamics in coastal areas and in particular their extremes, through the study of long-term biodeochemical time series registered by automated systems. The majority of high-frequency data sets used in this study came from MAREL program. The low-frequency time series from coastal monitoring programs SOMLIT (CNRS, INSU) and SRN (Ifremer) are employed to support the importance of automated systems. The EMD (Empirical Mode decomposition) method has provided a basis for us to study several of these time series. We also have used some methods more classical borrowed from numerical analysis field and turbulence. This study is organized in three chapters, and several appendices. The first chapter is devoted to the material and method. In the second chapter, using the EMD method we have highlighted the strong fluctuations contained in the blooms, and we have performed spectral analyzes. The principal component analysis (PCA) highlighted the main forcing exerted on primary production and SOMLIT temperature profiles suggest an impact of stratification on the intensity of blooms. In the third chapter, we conducted a comparative study between low-frequency and high-frequency data. Two cross-correlation methods (TDIC and co-spectra) allowed us to define a characteristic transition scale between the temperatures of the western and eastern English Cahnnel. In appendices we tested the robustness of different spectral analysis methods about the missing data in the time series, which is an underlying problem in the database registered by automated systems, and we reproduce a paper, which is under submission.

Cette thèse a pour objet la simulation d'une zone de mélange turbulente issue de l'instabilité de Richtmyer-Meshkov à l'aide d'un modèle à fonction de densité de probabilité (PDF). Nous analysons plus particulièrement la prise en charge par le modèle PDF du transport de l'énergie cinétique turbulente dans la zone de mélange.Dans cette optique, nous commençons par mettre en avant le lien existant entre les statistiques en un point de l'écoulement et ses conditions initiales aux grandes échelles. Ce lien s'exprime à travers le principe de permanence des grandes échelles, et permet d'établir des prédictions pour certaines grandeurs de la zone de mélange, telles que son taux de croissance ou son anisotropie.Nous dérivons ensuite un modèle PDF de Langevin capable de restituer cette dépendance aux conditions initiales. Ce modèle est ensuite validé en le comparant à des résultats issus de simulations aux grandes échelles (LES).Enfin, une analyse asymptotique du modèle proposé permet d'éclairer notre compréhension du transport turbulent. Un régime de diffusion est mis en évidence, et l'expression du coefficient de diffusion associé à ce régime atteste l'influence de la permanence des grandes échelles sur le transport turbulent.Tout au long de cette thèse, nous nous sommes appuyés sur des résultats issus de simulations de Monte Carlo du modèle de Langevin. A cet effet, nous avons développé une méthode spécifique eulérienne et à l'avons comparé à des alternatives lagrangiennes. / The aim of the thesis is to simulate a turbulent mixing zone resulting from the Richtmyer-Meshkov instability using a probability density function (PDF) model. An emphasis is put on the analysis of the turbulent kinetic energy transport.To this end, we first highlight the link existing between the one-point statistics of the flow and its initial conditions at large scales. This link is expressed through the principle of permanence of large eddies, and allows to establish predictions for quantities of the mixing zone, such as its growth rate or its anisotropy.We then derive a Langevin PDF model which is able to reproduce this dependency of the statistics on the initial conditions. This model is then validated by comparing it against large eddy simulations (LES).Finally, an asymptotic analysis of the derived model helps to improve our understanding of the turbulent transport. A diffusion regime is identified, and the expression of the diffusion coefficient associated with this regime confirms the influence of the permanence of large eddies on the turbulent transport.Throughout this thesis, our numerical results were based on Monte Carlo simulations for the Langevin model. In this regard, we proceeded to the development of a specific Eulerian method and its comparison with Lagrangian counterparts.

La thèse porte sur le développement des techniques non linéaires robustes de stabilisation et commande des systèmes avec perturbations de model. D’abord, on introduit les concepts de base de stabilité et stabilisabilité robuste dans le contexte des systèmes non linéaires. Ensuite, on présente une méthodologie de stabilisation par retour d’état en présence d’incertitudes qui ne sont pas dans l’image de la commande («unmatched»). L’approche récursive du «backstepping» permet de compenser les perturbations «unmatched» et de construire une fonction de Lyapunov contrôlée robuste, utilisable pour le calcul ultérieur d’un compensateur des incertitudes dans l’image de la commande («matched»). Le contrôleur obtenu est appelé «recursive Lyapunov redesign». Ensuite, on introduit la technique de stabilisation par «Immersion & Invariance» comme outil pour rendre un donné contrôleur non linéaire, robuste par rapport à dynamiques non modelées. La première technique de contrôle non linéaire robuste proposée est appliquée au projet d’un autopilote pour un missile air-air et au développement d’une loi de commande d’attitude pour un satellite avec appendices flexibles. L’efficacité du «recursive Lyapunov redesign» est mis en évidence dans le deux cas d’étude considérés. En parallèle, on propose une méthode systématique de calcul des termes incertains basée sur un modèle déterministe d’incertitude. La partie finale du travail de thèse est relative à la stabilisation des systèmes sous échantillonnage. En particulier, on reformule, dans le contexte digital, la technique d’Immersion et Invariance. En premier lieu, on propose des solutions constructives en temps continu dans le cas d’une classe spéciale des systèmes en forme triangulaire «feedback form», au moyen de «backstepping» et d’arguments de domination non linéaire. L’implantation numérique est basée sur une loi multi-échelles, dont l’existence est garantie pour la classe des systèmes considérée. Le contrôleur digital assure la propriété d’attractivité et des trajectoires bornées. La loi de commande, calculée par approximation finie d’un développement asymptotique, est validée en simulation de deux exemples académiques et deux systèmes physiques, le pendule inversé sur un chariot et le satellite rigide. / The dissertation deals with the problems of stabilization and control of nonlinear systems with deterministic model uncertainties. First, in the context of uncertain systems analysis, we introduce and explain the basic concepts of robust stability and stabilizability. Then, we propose a method of stabilization via state-feedback in presence of unmatched uncertainties in the dynamics. The recursive backstepping approach allows to compensate the uncertain terms acting outside the control span and to construct a robust control Lyapunov function, which is exploited in the subsequent design of a compensator for the matched uncertainties. The obtained controller is called recursive Lyapunov redesign. Next, we introduce the stabilization technique through Immersion \& Invariance (I\&I) as a tool to improve the robustness of a given nonlinear controller with respect to unmodeled dynamics. The recursive Lyapunov redesign is then applied to the attitude stabilization of a spacecraft with flexible appendages and to the autopilot design of an asymmetric air-to-air missile. Contextually, we develop a systematic method to rapidly evaluate the aerodynamic perturbation terms exploiting the deterministic model of the uncertainty. The effectiveness of the proposed controller is highlighted through several simulations in the second case-study considered. In the final part of the work, the technique of I\& I is reformulated in the digital setting in the case of a special class of systems in feedback form, for which constructive continuous-time solutions exist, by means of backstepping and nonlinear domination arguments. The sampled-data implementation is based on a multi-rate control solution, whose existence is guaranteed for the class of systems considered. The digital controller guarantees, under sampling, the properties of manifold attractivity and trajectory boundedness. The control law, computed by finite approximation of a series expansion, is finally validated through numerical simulations in two academic examples and in two case-studies, namely the cart-pendulum system and the rigid spacecraft.

Lors de l’infection par un pathogène intracellulaire, l’organisme déclenche une réponse immunitaire spécifique dont les acteurs principaux sont les lymphocytes T-CD8. Ces cellules sont responsables de l’éradication de ce type d’infections et de la constitution du répertoire immunitaire de l’individu. Les processus qui composent la réponse immunitaire se répartissent sur plusieurs échelles physiques inter-connectées (échelle intracellulaire, échelle d’une cellule, échelle de la population de cellules). La réponse immunitaire est donc un processus complexe, pour lequel il est difficile d’observer ou de mesurer les liens entre les différents phénomènes mis en jeu. Nous proposons trois modèles mathématiques multi-échelles de la réponse immunitaire, construits avec des formalismes différents mais liés par une même idée : faire dépendre le comportement des cellules TCD8 de leur contenu intracellulaire. Pour chaque modèle, nous présentons, si possible, sa construction à partir des hypothèses biologiques sélectionnées, son étude mathématique et la capacité du modèle à reproduire la réponse immunitaire au travers de simulations numériques. Les modèles que nous proposons reproduisent qualitativement et quantitativement la réponse immunitaire T-CD8 et constituent ainsi de bons outils préliminaires pour la compréhension de ce phénomène biologique. / Upon infection by an intracellular pathogen, the organism triggers a specific immune response,mainly driven by the CD8 T cells. These cells are responsible for the eradication of this type of infections and the constitution of the immune repertoire of the individual. The immune response is constituted by many processes which act over several interconnected physical scales (intracellular scale, single cell scale, cell population scale). This biological phenomenon is therefore a complex process, for which it is difficult to observe or measure the links between the different processes involved. We propose three multiscale mathematical models of the CD8 immune response, built with different formalisms but related by the same idea : to make the behavior of the CD8 T cells depend on their intracellular content. For each model, we present, if possible, its construction process based on selected biological hypothesis, its mathematical study and its ability to reproduce the immune response using numerical simulations. The models we propose succesfully reproduce qualitatively and quantitatively the CD8 immune response and thus constitute useful tools to further investigate this biological phenomenon.

Cette thèse constitue une analyse de la dynamique du bois mort en rivière par une approche de bilan ligneux multi-scalaire dans une rivière semi-alluviale de région froide: la rivière Saint-Jean, Gaspé. L’étude de la variabilité spatiale et temporelle de la dynamique du bois mort a été réalisée grâce à une approche méthodologique combinant quatre années de terrain et par l’analyse de documents historiques. Les rivières de la péninsule gaspésienne produisent et transportent annuellement de grandes quantités de bois mort. Cette production provient à la fois de la puissance spécifique élevée des cours d’eau et de berges composées de sédiments peu cohésifs et comportant une ripisylve arborescente généralement dense. Jusqu’en 2015, le delta de la rivière Saint-Jean comportait plusieurs embâcles de bois de très grandes tailles. Ces embâcles se mettent en place depuis les années 1960, ils représentent une opportunité unique de quantifier et d’appliquer un budget ligneux pour cerner les variables clefs liées à la dynamique du bois mort à multiples échelles spatiales et temporelles.Nos résultats montrent que la quasi-totalité des bois morts en rivière est produite par la migration latérale et par l’influence de la morphologie. Au niveau des accumulations dans le corridor fluvial, deux zones accumulent la majorité des bois et ces premiers résultats montrent une mobilité importante qui peut varier fortement d'une année à l'autre. Au niveau de la mobilité, l’analyse d’imagerie vidéo de trois évènements différents ont montré que l’intensité (nombre de bois par minute) du transport de bois pouvait être jusqu’à dix fois supérieures durant un évènement avec débâcle mécanique glacielle, comparativement à un évènement hydroclimatique en eau libre.Finalement, nous avons réalisé une analyse par bilans ligneux et une analyse de la trajectoire éco-hydrophormologique sur une période de plus de 50 ans. Nous avons quantifié chaque composante (intrant, accumulée et sortie) d'un budget ligneux à multiples échelles spatio-temporelles. Autant à l’échelle interannuelle, qu’à l’échelle décennale, la dynamique du bois mort comporte des périodes où l’entrée de bois, le stockage et la mobilité diffèrent selon les contextes éco-hydromorphologiques. La trajectoire éco-hydromorphologiques suggère une augmentation de la dynamique fluviale en raison d'un changement significatif dans l'hydrologie conduisant à une plus grande production et mobilité des bois morts et par l'augmentation des volumes accumulés dans le corridor de la RSJ depuis 2004. / This thesis is an analysis of the dynamics of large wood in river by a multi-scale large wood budget approach in a semi-alluvial river of cold regions: the Saint-Jean River, Gaspé. The study of the spatial and temporal variability of the dynamics of large wood was carried out through a methodological approach combining four years of field and by analyzing historical documents. The rivers of the Gaspé Peninsula produce annually and carry large amounts of large wood. This production comes from the high specific power of rivers and by banks composed of noncohesive sediment and having a generally dense riparian forest tree. Until 2015, the Delta of the Saint-Jean River had several very large jams. These jams are put in place since the 1960s, they represent a unique opportunity to quantify and apply a wood budget and to identify key variables related to the dynamics of large wood at multiple spatial and temporal scales.Our results show that almost all large wood in river is produced by lateral migration and by the influence of the morphology. For the accumulations in the river corridor, two areas accumulate the majority of wood and these first results show a significant mobility, which can fluctuate substantially from year to year. In terms of mobility, video analysis of three different events showed that the intensity of the transport (number of large wood per minute) can be higher to ten times during an event with mechanical ice-breakup, compared to an open water.hydroclimatic event.Finally, we conducted an analysis by large wood budget and analysis of eco-hydromorphological trajectory over more than 50 years. We quantified each component (input, output and accumulated) of a large wood budget at multiple spatial and temporal scales. At the interannual scale or decadal, scale, the dynamics of large wood have periods when the input, storage and mobility of large wood differ according to eco-hydromorphological contexts. The eco-hydromorphological trajectory suggests an increase in river dynamics due to a significant change in the hydrology, resulting in higher production and mobility of large wood and increased in volumes accumulated in the corridor of the Saint-Jean River since the last decade.

Notre thèse porte sur la Modélisation d'Entreprise et plus particulièrement sur une vision Ingénierie Système de celle-ci afin de faciliter son application « à minima » en entreprise. Alors que la modélisation d'entreprise arrive à maturité scientifique et technologique, force est de constater que sa pratique est encore trop peu courante. En outre, il n'existe pas à proprement parler de tâche de modélisation, de maintenance et de gestion de modèles dans l'entreprise. La question s'est alors posée de la relation d'interopération entre les systèmes, d'un Système-Entreprise évolutif, le plus souvent composé de façon ad-hoc de sous-systèmes (COTS), qui ont pour particularité d'encapsuler leur propre Savoir-Faire et leurs Ingénieries respectives, en incluant le Système-Produit (son besoin, son projet, ...) à la source de chaque recomposition « à la demande » du Système-Entreprise. Un des objectifs de nos travaux est donc de faciliter l'interopération entre les ingénieries de ces systèmes interopérants (autour du produit A Faire) et l'Ingénierie du Système-Projet d'Entreprise (Pour Faire). Nous nous sommes alors inspirés de la proposition de Kuras pour formaliser notre Système-Entreprise (en fait son Ingénierie) dans une sorte de « Système de Systèmes d'Information » (SdSI) en tenant compte de ses aspects récursifs et multi-échelles dans un contexte organisationnel évolutif. De plus, en nous appuyant sur la théorie des patterns, nous avons identifié puis dérivé le « composite-pattern» pour mieux décrire la nature « Tout et Partie » de l'Holon Système-Entreprise utilisé par Kuras pour extraire le patron de conceptualisation le plus adéquat. Parmi tous les points de vue d'un système, cette dérivation nous permet de représenter ceux du Système A Faire et ceux du Système Pour Faire avec son contenu d'Ingénierie associée guidée par le cadre de modélisation de Zachman. Nous avons outillé notre proposition de cadre d'Ingénierie Système Basée sur des Modèles (ISBM) en étendant le méta-modèle de l'outil « MEGA Modelling Suite » ainsi qu'en développant les interfaces utilisateurs nécessaires et nous l'avons appliqué à l'Ingénierie d'un Système de Traçabilité d'un Produit dans le contexte d'un scénario PLM plausible entre l'AIPL et DIMEG.

Nous abordons le développement d'une nouvelle génération de méthodes numériques pour la résolution des EDP évolutives qui modélisent des phénomènes multi-échelles en temps et en espace issus de divers domaines applicatifs. La raideur associée à ce type de problème, que ce soit via le terme source chimique qui présente un large spectre d'échelles de temps caractéristiques ou encore via la présence de fort gradients très localisés associés aux fronts de réaction, implique en général de sévères difficultés numériques. En conséquence, il s'agit de développer des méthodes qui garantissent la précision des résultats en présence de forte raideur en s'appuyant sur des outils théoriques solides, tout en permettant une implémentation aussi efficace. Même si nous étendons ces idées à des systèmes plus généraux par la suite, ce travail se focalise sur les systèmes de réaction-diffusion raides. La base de la stratégie numérique s'appuie sur une décomposition d'opérateur spécifique, dont le pas de temps est choisi de manière à respecter un niveau de précision donné par la physique du problème, et pour laquelle chaque sous-pas utilise un intégrateur temporel d'ordre élevé dédié. Ce schéma numérique est ensuite couplé à une approche de multirésolution spatiale adaptative permettant une représentation de la solution sur un maillage dynamique adapté. L'ensemble de cette stratégie a conduit au développement du code de simulation générique 1D/2D/3D académique MBARETE de manière à évaluer les développements théoriques et numériques dans le contexte de configurations pratiques raides issue de plusieurs domaines d'application. L'efficacité algorithmique de la méthode est démontrée par la simulation d'ondes de réaction raides dans le domaine de la dynamique chimique non-linéaire et dans celui de l'ingénierie biomédicale pour la simulation des accidents vasculaires cérébraux caractérisée par un terme source "chimique complexe''. Pour étendre l'approche à des applications plus complexes et plus fortement instationnaires, nous introduisons pour la première fois une technique de séparation d'opérateur avec pas de temps adaptatif qui permet d'atteindre une précision donnée garantie malgré la raideur des EDP. La méthode de résolution adaptative en temps et en espace qui en résulte, étendue au cas convectif, permet une description consistante de problèmes impliquant une très large palette d'échelles de temps et d'espace et des scénarios physiques très différents, que ce soit la propagation des décharges répétitives pulsées nanoseconde dans le domaine des plasmas ou bien l'allumage et la propagation de flammes dans celui de la combustion. L'objectif de la thèse est l'obtention d'un solveur numérique qui permet la résolution des EDP raides avec contrôle de la précision du calcul en se basant sur des outils d'analyse numérique rigoureux, et en utilisant des moyens de calculs standard. Quelques études complémentaires sont aussi présentées comme la parallélisation temporelle, des techniques de parallélisation à mémoire partagée et des outils de caractérisation mathématique des schémas de type séparation d'opérateur.

Les activités de recherche qui ont été développées s'inscrivent dans une logique d'amélioration de la modélisation/simulation des fluides denses multiconstituants non réactifs, aussi bien d'un point de vue statique que dynamique, en utilisant une démarche de type échelle microscopique vers échelle macroscopique (bottom-up). Un accent tout particulier à été mis sur l'étude des fluides d'intérêt pétrolier dans des conditions de gisement et ce par le biais de simulations numériques de la dynamique moléculaire classique (DM).<br />Les premiers travaux réalisés ont porté essentiellement sur la simulation numérique de la DM appliquée à l'étude de la thermodiffusion, ou effet Soret, en fluide libre et en milieu poreux. De ce travail est notamment ressortit que le confinement géométrique n'affecte pas le facteur de thermodiffusion et que l'effet Soret dans les mélanges d'alcanes est principalement du à un effet de masse.<br />Dans l'optique mentionnée dans l'introduction, trois thèmes ont été ensuite abordés. Deux thèmes ont trait à la dynamique moléculaire classique, l'un concernant l'influence de la description des interactions sur les propriétés thermophysiques dans les fluides et l'autre la simulation/modélisation de propriétés de transport du fluide de Lennard-Jones et son application aux corps réels. En particulier il a été montré qu'un état correspondant existait entre les potentiels de type Mie et Exponentiel au niveau des propriétés de transport et que Lennard-Jones (LJ) était un bon compromis pour les corps simples. De même ont été développées des corrélations basées sur le fluide de LJ permettant d'estimer les propriétés de transport. Cette approche s'est avérée particulièrement efficace quant à l'estimation de la viscosité des gaz acides (H2S, CO2). Le troisième thème concerne la modélisation macroscopique de la dynamique de fluides multiconstituants en milieux poreux pour lequel une approche de type équation Darcy généralisée par espèce a été proposée. <br />Les activités en cours concernent la prédiction/modélisation de fluides plus complexes avec ou sans interfaces mais également des aspects liés aux simulations multi-échelles/multi-physiques au travers d'un couplage direct entre mécanique des fluides numérique et DM. En particulier à été proposée une approche permettant l'étude de la thermogravitation dans les gisements par simulation moléculaire.

Cette thèse propose différentes méthodes théoriques et numériques pour simuler à coût réduit le comportement des plasmas ou des faisceaux de particules chargées sous l'action d'un champ magnétique fort. Outre le champ magnétique externe, chaque particule est soumise à un champ électromagnétique créé par les particules elles-mêmes. Dans les modèles cinétiques, les particules sont représentées par une fonction de distribution f(x,v,t) qui vérifie l'équation de Vlasov. Afin de déterminer le champ électromagnétique, cette équation est couplée aux équations de Maxwell ou de Poisson. L'aspect champ magnétique fort est alors pris en compte par un adimensionnement adéquat qui fait apparaître un paramètre de perturbation singulière 1/ε. Le premier chapitre de cette thèse est une introduction à la fusion contrôlée par confinement magnétique dans les Tokamaks. Le second chapitre est consacré à la théorie gyrocinétique géométrique. Cette théorie repose sur la géométrie différentielle et la dynamique des systèmes hamiltoniens. L'objectif est de faire une succession de changements de coordonnées afin de se ramener à un système proche du centre-guide historique dans lequel les expressions de la matrice de Poisson et du Hamiltonien permettent une réduction de la dimension des trajectoires. Le chapitre 3 met en pratique les mêmes techniques sur un autre problème, la modélisation paraxiale d'un faisceau de particules chargées. Le dernier chapitre est dédié à un schéma numérique basé sur un intégrateur exponentiel en vitesse. Ce schéma a pour objectif d'approcher numériquement des solutions fortement oscillantes avec une méthode Particle-In-Cell en utilisant un pas de temps beaucoup plus grand que la période d'oscillation rapide. Il est testé sur une équation de Vlasov linéaire ainsi que sur le système de Vlasov-Poisson.

Cette thèse pose les principes d'une utilisation de la réduction de modèle dans un contexte de décomposition de domaine en dynamique transitoire des structures. L'intégration temporelle est explicite, dédiée à la prise en compte de phénomènes rapides tels que des impacts. L'approche multi-domaines mise en oeuvre permet de découpler les échelles de temps entre les sous-domaines, par l'intermédiaire d'une gestion des incompatibilités de maillage aux frontières et d'une stratégie à plusieurs pas de temps. La réduction est introduite dans le formalisme de façon totalement générique via la projection de l'équilibre sur un sous-domaine donné sur une base modale locale et sans modification de l'algorithme générale, les opérateurs éléments finis initiaux étant remplacés par leur projection sur la base et les relations de continuité aux frontières entre les domaines étant reconstruites dans l'espace physique. La cas particulier de la non-linéarité géométrique issue de sous-domaines en vibration libres de blocage et subissant des rotations d'ensemble d'amplitude finie fait l'objet d'une analyse spécifique, aboutissant à une réduction à deux échelles. La connexion entre un sous-domaine ainsi projeté et les autres sous-domaines présentent des particularités au niveau du coût numérique du calcul des efforts d'interface, ce qui fait l'objet de l'élaboration d'une technique de résolution totalement nouvelle. La pertinence et les performances, en résolution séquentielle, du formalisme proposé sont illustrées sur des exemples issus de l'industrie en support de la thèse.