Dans le contexte industriel actuel, où la simulation numérique joue un rôle majeur, de nombreux outils sont développés afin de rendre les calculs les plus performants et exacts possibles en utilisant les ressources numériques de façon optimale. Parmi ces outils, ceux non-intrusifs, c’est-à-dire ne modifiant pas les codes commerciaux disponibles mais permettant d’utiliser des méthodes de résolution avancées telles que l’analyse isogéométrique ou les couplages multi-échelles, apparaissent parmi les plus attirants pour les industriels. L’objectif de cette thèse est ainsi de coupler l’Analyse IsoGéométrique (AIG) et la Méthode des Éléments Finis (MEF) standard pour l’analyse de détails structuraux par une approche non-intrusive et certifiée. Dans un premier temps, on développe un lien global approché entre les fonctions de Lagrange, classiquement utilisées en éléments finis et les fonctions NURBS bases de l’AIG, ce qui permet d’implémenter des analyses isogéométriques dans un code industriel EF vu comme une boîte noire. Au travers d’exemples linéaires et non-linéaires implémentés dans le code industriel Code_Aster de EDF, nous démontrons l’efficacité de ce pont AIG\MEF et les possibilités d’applications industrielles. Il est aussi démontré que ce lien permet de simplifier l’implémentation du couplage non-intrusif entre un problème global isogéométrique et un problème local éléments finis. Ensuite, le concept de couplage non-intrusif entre les méthodes étant ainsi possible, une stratégie d’adaptation est mise en place afin de certifier ce couplage vis-à-vis d’une quantité d’intérêt. Cette stratégie d’adaptation est basée sur des méthodes d’estimation d’erreur a posteriori. Un estimateur global et des indicateurs d’erreur d’itération, de modèle et de discrétisation permettent de piloter la définition du problème couplé. La méthode des résidus est utilisée pour évaluer ces erreurs dans des cas linéaires, et une extension aux problèmes non-linéaires via le concept d’Erreur en Relation de Comportement (ERC) est proposée. / In the current industrial context where the numerical simulation plays a major role, a large amount of tools are developed in order to perform accurate and effective simulations using as less numerical resources as possible. Among all these tools, the non-intrusive ones which do not modify the existing structure of commercial softwares but allowing the use of advanced solving methods, such as isogeometric analysis or multi-scale coupling, are the more attractive to the industry. The goal of these thesis works is thus the coupling of the Isogeometric Analysis (IGA) with the Finite Element Method (FEM) to analyse structural details with a non-intrusive and certified approach. First, we develop an approximate global link between the Lagrange functions, commonly used in the FEM, and the NURBS functions on which the IGA is based. It’s allowed the implementation of isogeometric analysis in an existing finite element industrial software considering as a black-box. Through linear and nonlinear examples implemented in the industrial software Code_Aster of EDF, we show the efficiency of the IGA\FEM bridge and all the industrial applications that can be made. This link is also a key to simplify the non-intrusive coupling between a global isogeometric problem and a local finite element problem. Then, as the non-intrusive coupling between both methods is possible, an adaptive process is introduced in order to certify this coupling regarding a quantity of interest. This adaptive strategy is based on a posteriori error estimation. A global estimator and indicators of iteration, model and discretization error sources are computed to control the definition of the coupled problem. Residual base methods are performed to estimated errors for linear cases, an extension to the concept of constitutive relation errors is also initiated for non-linear problems.

Des sites de visionnement de contenu audio-vidéo en temps-réel comme YouTube sont devenus très populaires. Le téléchargement des fichiers audio/vidéo consomme une quantité importante de bande passante des réseaux Internet. L’utilisation de codecs à bas débit permet de compresser la taille des fichiers transmis afin de consommer moins de bande passante. La conséquence est une diminution de la qualité de ce qui est transmis. Une diminution de qualité mène à l’apparition de défauts perceptibles dans les fichiers. Ces défauts sont appelés des artifices de compression. L’utilisation d’un algorithme de post-traitement sur les fichiers sonores pourrait augmenter la qualité perçue de la musique transmise en corrigeant certains artifices à la réception, sans toutefois consommer davantage de bande passante. Pour rehausser la qualité subjective des fichiers sonores, il est d’abord nécessaire de déterminer quelles caractéristiques dégradent la qualité perceptuelle. Le présent projet a donc pour objectif le développement d’un algorithme capable de localiser et de corriger de façon non intrusive, un artifice provoqué par des discontinuités et des incohérences au niveau des harmoniques qui dégrade la qualité objective dans les signaux sonores compressés à bas débits (8 – 12 kilobits par seconde).

A l'Ouest d'Ardebil, le Sabalan (4820 ml. massif volcanique central représente un ensemble subactif complexe. Le socle tectonisé de ce massif est constitué par des séries volcaniques d'âge Crétacé supérieur au Miocène, qui prolongent la ceinture volcanique d'Iran central. Après une phase tectonique ayant fortement plissé les assises anté-Crétacé supérieur (phase laramienne), des mouvements d'extension (andins? ), localisés principalement dans le Nord de cette région, provoquent la montée d'un magma alcalin sous-saturé (blairmorites) dont les épanchements sont associés à des calcaires pélagiques sénoniens. Les mouvements orogéniques du Danien permettent une émersion assez générale à laquelle fait suite une phase d'expansion importante. Cette phase, accompagnée de basculements traduit l'ouverture d'un "rift" intracontinental qui se poursuivra durant tout l'Eocène. Par des fissures profondes, l'ascension d'un magma alcalin donne naissance à des formations volcaniques d'épaisseur et d'étendue considérables. La différenciation complexe et poussée de ce magma se traduit par des laves très diversifiées (de latites basaltiques à des dacites).appartenant à l'association shoshonitique. Une transgression, du Nord vers le Sud, associée à l'extension du "rift", atteint son développement maximum au Lutétien. La fin de l'Eocène est caractérisée par l'arrêt de l'expansion, c'est-à-dire par la cicatrisation du "rift", accompagnée par une nouvelle émersion générale. Au début de l'Oligocène, une anatexie de la croûte se traduit par une phase de granitisation profonde, consécutive à la cicatrisation du "rift" Eocène. Le "rift", ainsi que l'anatexie,peuvent être l'expression de la présence d'un "dôme thermique" sous l'Azerbaïdjan oriental durant toute la période éocène. La granitisation oligocène se traduit par la mise en place d'énormes volumes de roches intrusives (monzonitesl, auxquelles sont liés quelques épanchements acides réduits. La fin de l'Oligocène est caractérisée par une phase tectonique cassante, durant laquelle s'individualisent de grands horsts d'orientation grossièrement Est-Ouest. Au Miocène et au Plio-Quaternaire, des mouvements tectoniques (contemporains de la surrection de la chaîne du Caucase) provoquent une intense fracturation des horsts. Des appareils volcaniques ponctuels se mettent en place le long de ces fractures miocènes, de direction méridienne. Sur l'un, des horsts se développe le Sabalan, affecté, au cours de son édification,par l'effondrement d'une caldera. La construction au strato-volcan anté-caldera s'est faite en cinq épisodes successifs, précédés chacun de phases explosives réduites. Les caractères volcanologiques des séries sont constants: dynamisme de type strombolien accompagné de coulées très étendues, avec un paroxisme effusif pour les séries intermédiaires.

Le concept de l'avion plus électrique a conduit les industriels à focaliser une partie de leurs efforts sur les risques liés à l'utilisation de la haute tension dans un environnement sévère (Basse pression, large plage de température et d'humidité …). Les risques associés sont l'existence et le développement de décharges partielles conduisant à terme à la défaillance du système dans lequel, elles se produisent. Considérant cette problématique, le groupe Liebherr Elektronik GmbH, en collaboration avec le laboratoire Laplace, a lancé l'étude d'une méthode permettant de détecter les décharges partielles dans des convertisseurs destinés à des applications aéronautiques. Ce sont les résultats de ce travail qui font l'objet de cette thèse. La première partie fait état des connaissances actuelles sur les décharges partielles, de leur nature physique aux méthodes de détection utilisées. Dans la seconde partie, trois phases de mesures validant l'efficacité de la méthode étudiée sont présentées et leurs résultats sont discutés. Enfin, dans une troisième partie, les conclusions de nos travaux et leurs perspectives sont présentées. / The more electrical airplane concept led industrial companies to focus a part of their efforts on risks linked to the use of high voltage in a severe environment (Low pressure, wide range of temperature and humidity …). Associated risks are the existence and the growing of partial discharges ultimately leading to the breakdown of the system in which they occur. Considering this problematic, the Liebherr Elektronik GmbH group, in collaboration with the Laplace laboratory, launched the study of a method allowing partial discharges detection in converters intended to be used in aeronautical applications. The results of this work are the subject of this thesis. The first part brings the status of the current knowledge about partial discharges from their physical nature to the detection methods which are used. In the second part, three measurement phases assessing the efficiency of the studied method are introduced and their results discussed. Finally, in the third part, the conclusionsof our works and their perspectives are presented.

Cette thèse s'inscrit dans le cadre des méthodes de couplage non-intrusives qui permettent de concilier deux types de logiciels : logiciels commerciaux généralistes et logiciels dédiés. L'objectif est d'appliquer un tel couplage pour le traitement des assemblages de plaques stratifiées.En effet, le comportement aux environs des zones critiques de ces matériaux composites nécessite une prise en compte fine (3D) des phénomènes en jeu alors que la pratique industrielle consiste à représenter les structures par des modèles simplifiés de plaques ou de coques.La méthode non-intrusive développée ici permet ainsi de « patcher » virtuellement le modèle plaque par des modèles 3D locaux, de coupler fortement le modèle global avec les patchs, sans modifier les maillages ni les modèles.Diverses stratégies de raccords sont envisagées, ainsi que des améliorations à la technique de couplage. / This thesis is part of the nonintrusive coupling methods framework which allow to merge to kind of software : general commercial software and dedicated software.The aim is to apply such coupling for the analysis of assembly of laminated plates.In fact, the behavior in composite material around critical zones requires to take in account complex phenomena (in 3D), whereas the industrial way of work is to simplify the structure with plate or shell model in order to handle very large structures.The nonintrusive coupling method developed here allow to virtually apply a 3D patch to the plate model in order to enriched the global description without any modification of the global meshes, models, or solver.Various methods of coupling and improvement to the classical nonintrusive coupling are proposed.

Les assemblages sont des éléments critiques pour les structures industrielles. De fortes non-linéarités de type contact frottant, ainsi que des précharges mal maîtrisées rendent complexe tout dimensionnement précis. Présents en très grand nombre sur les structures industrielles (quelques millions pour un A380), cela implique de rafiner les modèles localement et donc de gérer des problèmes numé-riques de très grandes tailles. Les nombreuses interfaces de contact frottant sont des sources de difficultés de convergence pour les simulations numériques. Il est donc nécessaire de faire appel à des méthodes robustes. Il s’agit d’utiliser des méthodes itératives de décomposition de domaine, permettant de gérer des modèles numériques extrêmement grands, couplées à des techniques adaptées afin de prendre en compte les non-linéarités de contact aux interfaces entre sous-domaines. Ces méthodes de décomposition de domaine restent encore très peu utilisées dans un cadre industriel. Des développements internes aux codes éléments finis sont souvent nécessaires et freinent ce transfert du monde académique au monde industriel.Nous proposons, dans ces travaux de thèse, une mise-en-oeuvre non intrusive de ces méthodes de décomposition de domaine : c’est-à-dire sans développement au sein du code source. En particulier, nous nous intéressons à la méthode Latin dont la philosophie est particulièrement adaptée aux problèmes non linéaires. La structure est décomposée en sous-domaines reliés entre eux au travers d’interfaces. Avec la méthode Latin, les non-linéarités sont résolues séparément des aspects linéaires. La résolution est basée sur un schéma itératif à deux directions de recherche qui font dialoguer les problèmes linéaires globaux etles problèmes locaux non linéaires.Au cours de ces années de thèse, nous avons développé un outil totalement non intrusif sous Code_Aster permettant de résoudre par une technique de décomposition de domaine mixte des problèmes d’assemblage. Les difficultés posées par le caractère mixte de la méthode Latin sont résolues par l’introduction d’une direction de recherche non locale. Des conditions de Robin sur les interfaces des sous-domaines sont alors prises en compte simplement sans modifier les sources de Code_Aster. Nous avons proposé une réécriture algébrique de l’approche multi-échelle assurant l’extensibilité de la méthode. Nous nous sommes aussi intéressés à coupler la méthode Latin en décomposition de domaine à un algorithme de Krylov. Appliqué uniquement à un problème sous-structuré avec interfaces parfaites, ce couplage permet d’accélérer la convergence. Des structures préchargées avec de nombreuses interfaces de contact frottant ont été traitées. Des simulations qui n’auraient pu être menées par un calcul direct sous Code_Aster ont été réalisées via cette stratégie de décomposition de domaine non intrusive. / Abstract : Assemblies are critical elements for industrial structures. Strong non-linearities such as frictional contact, as well as poorly controlled preloads make complex all accurate sizing. Present in large numbers on industrial structures (a few million for an A380), this involves managing numerical problems of very large size. The numerous interfaces of frictional contact are sources of difficulties of convergence for the numerical simulations. It is therefore necessary to use robust but also reliable methods. The use of iterative methods based on domain decomposition allows to manage extremely large numerical models. This needs to be coupled with adaptedtechniques in order to take into account the nonlinearities of contact at the interfaces between subdomains. These methods of domain decomposition are still scarcely used in industries. Internal developments in finite element codes are often necessary, and thus restrain this transfer from the academic world to the industrial world.In this thesis, we propose a non-intrusive implementation of these methods of domain decomposition : that is, without development within the source code. In particular, we are interested in the Latin method whose philosophy is particularly adapted to nonlinear problems. It consists in decomposing the structure into sub-domains that are connected through interfaces. With the Latin method the non-linearities are solved separately from the linear differential aspects. Then the resolution is based on an iterative scheme with two search directions that make the global linear problems and the nonlinear local problems dialogue.During this thesis, a totally non-intrusive tool was developed in Code_Aster to solve assembly problems by a mixed domain decomposition technique. The difficulties posed by the mixed aspect of the Latin method are solved by the introduction of a non-local search direction. Robin conditions on the subdomain interfaces are taken into account simply without modifying the sources of Code_Aster. We proposed an algebraic rewriting of the multi-scale approach ensuring the extensibility of the method. We were also interested in coupling the Latin method in domain decomposition to a Krylov algorithm. Applied only to a substructured problem with perfect interfaces, this coupling accelerates the convergence. Preloaded structures with numerous contact interfaces have been processed. Simulations that could not be carried out by a direct computationwith Code_Aster were performed via this non-intrusive domain decomposition strategy.

Cette thèse vise à proposer des outils innovants pour le calcul de structures aéronautiques évoluant à haute température. En effet, les régimes de fonctionnement des moteurs actuels conduisent à des évolutions élasto-viscoplastiques généralisées dans les pièces métalliques et l’utilisation de modèles simplifiés (élastiques) n’est plus totalement satisfaisante en terme de précision, même en phase de préconception. De même, la géométrie complexe permettant le refroidissement continu des pièces (micro-perforations) doit être prise en compte de manière exacte. Les techniques de calcul standard pour ce genre de problème conduiraient à des simulations lentes et peu flexibles (la moindre modification entraînant une remise en œuvre complète de la chaîne de calcul). Plus précisément, cette thèse étend les méthodes de type global/local non-intrusives au cas de la viscoplasticité généralisée en utilisant deux échelles de temps et d'espace, chacune adaptée aux phénomènes locaux et globaux à capturer. La méthode est ensuite étendue au calcul de nombreux cycles complexes de chargement, par des techniques de saut de cycles. Le schéma de couplage en temps permet alors une adaptation locale du pas de temps par sous-domaine. Des techniques d’accélération de convergence sont proposées, à l’échelle d’un incrément puis à celle de la succession de cycles (sauts de cycles). Ces développements permettent d’obtenir rapidement et précisément une estimation du cycle limite qui alimente un modèle de durée de vie. Le couplage non-intrusif est réalisé dans un script de programmation pilotant un code commercial (dans notre cas le langage Python et Abaqus/Standard). La méthode a été appliquée sur des plateformes de calculs industrielles, en réutilisant directement des maillages et les mises en données issues de modèles intervenant plus tôt dans la chaîne de calcul. Un cas métier, issu d’un bureau d’études de Safran Aircraft Engines, a pu être traité. / This thesis consists in developing innovating tools destined to the simulation of aeronautical structures evolving at high temperature. Indeed, working rates of current engines lead to an elasto-viscoplastic evolution generalized in metallic parts and the use of simplified models (linear elastic) are no longer totally satisfying in term of accuracy, even in initial design process. Likewise, the complex geometry allowing the continuous cool down process of parts (micro-perforations) has to be exactly taken into account. The standard computation techniques dedicated to this kind of models would lead to slow simulations with a lack of flexibility (the slightest modifications leading to restart the whole design process of the computation chain).More precisely, this thesis extends the noninvasive global/local methods to the framework of viscoplasticity generalized to the whole structure, using two scales in time and space, each one adapted to global and local phenomena to capture. The method is then extended to the computation of high number of complex load cycles, by skipped cycles techniques. The time coupling scheme lets then a local adaptation of time steps per subdomain. Convergence acceleration techniques are also set up, first for one time step and then through several load cycles (skipped cycles). These developments conduct to obtain quickly an evaluation of the limit cycle providing data to a lifetime expectancy model.The noninvasive coupling is realized in a programming language script managing the commercial software (respectively in our case Python and Abaqus/Standard). The method has been applied on industrial computational platforms, by reusing directly meshes and data from previous engineering tasks appearing earlier in the computational chain. A genuine test case from a Safran Aircraft Engines design office, was performed successfully.

Les séries de l'anticlinal des Oulad Khalouk sont des séries détritiques continentales de l'Ordovicien supérieur. déposées en bordure de la plate-forme saharienne dans la zone de pro-delta circalittorale. L'orogenèse hercynienne, gouvernée par une tectonique de blocs de socle, débute par une phase de plissement synschisteux de direction axiale NE-SW, post-Namurien et ante-Westphalien moyen. Une deuxième phase cassante reprend les structures par des décrochements et des virgations : elle doit être rattachée à la mise en place du granite des Zaër (298 ma). Le métamorphisme régional culmine au droit du granite, à la limite épizone-mésozone; la mise en place du diapir granitique s'accompagne d'un métamorphisme de contact ainsi que d'une schistosité de strain-slip S2 reprenant la schistosité de flux SO-1. Les minéralisations en pechblende sont datées de l'Hercynien et correspondrait à un gisement primaire de type " contact ", hydrothermal de basse-température. La radioactivité de surface est due au piégeage de l'uranium sur des phases d' oxy*hydroxydes de fer.

Le projet Européen " GENIUS " ambitionne de développer les méthodologies génériques pour le diagnostic de systèmes piles à combustible à haute température de type oxyde solide (SOFC). Le travail de cette thèse s'intègre dans ce projet ; il a pour objectif la mise en oeuvre d'un outil de diagnostic en utilisant le stack comme capteur spécial pour détecter et identifierles défaillances dans les sous-systèmes du stack SOFC.Trois algorithmes de diagnostic ont été développés, se basant respectivement sur la méthode de classification k-means, la technique de décomposition du signal en ondelettes ainsi que la modélisation par réseau Bayésien. Le premier algorithme sert au diagnostic ex-situ et est appliqué pour traiter les donnés issues des essais de polarisation. Il permet de déterminer les variables de réponse significatives qui indiquent l'état de santé du stack. L'indice Silhouette a été calculé comme mesure de qualité de classification afin de trouver le nombre optimal de classes dans la base de données.La détection de défaut en temps réel peut se réaliser par le deuxième algorithme. Puisque le stack est employé en tant que capteur, son état de santé doit être vérifié préalablement. La transformée des ondelettes a été utilisée pour décomposer les signaux de tension de la pile SOFC dans le but de chercher les variables caractéristiques permettant d'indiquer l'état desanté de la pile et également assez discriminatives pour différentier les conditions d'opération normales et anormales.Afin d'identifier le défaut du système lorsqu'une condition d'opération anormale s'est détectée, les paramètres opérationnelles réelles du stack doivent être estimés. Un réseau Bayésien a donc été développé pour accomplir ce travail.Enfin, tous les algorithmes ont été validés avec les bases de données expérimentales provenant de systèmes SOFC variés, afin de tester leur généricité.

Le projet Européen « GENIUS » ambitionne de développer les méthodologies génériques pour le diagnostic de systèmes piles à combustible à haute température de type oxyde solide (SOFC). Le travail de cette thèse s’intègre dans ce projet ; il a pour objectif la mise en oeuvre d’un outil de diagnostic en utilisant le stack comme capteur spécial pour détecter et identifierles défaillances dans les sous-systèmes du stack SOFC.Trois algorithmes de diagnostic ont été développés, se basant respectivement sur la méthode de classification k-means, la technique de décomposition du signal en ondelettes ainsi que la modélisation par réseau Bayésien. Le premier algorithme sert au diagnostic ex-situ et est appliqué pour traiter les donnés issues des essais de polarisation. Il permet de déterminer les variables de réponse significatives qui indiquent l’état de santé du stack. L’indice Silhouette a été calculé comme mesure de qualité de classification afin de trouver le nombre optimal de classes dans la base de données.La détection de défaut en temps réel peut se réaliser par le deuxième algorithme. Puisque le stack est employé en tant que capteur, son état de santé doit être vérifié préalablement. La transformée des ondelettes a été utilisée pour décomposer les signaux de tension de la pile SOFC dans le but de chercher les variables caractéristiques permettant d’indiquer l’état desanté de la pile et également assez discriminatives pour différentier les conditions d’opération normales et anormales.Afin d’identifier le défaut du système lorsqu’une condition d’opération anormale s’est détectée, les paramètres opérationnelles réelles du stack doivent être estimés. Un réseau Bayésien a donc été développé pour accomplir ce travail.Enfin, tous les algorithmes ont été validés avec les bases de données expérimentales provenant de systèmes SOFC variés, afin de tester leur généricité. / The EU-project “GENIUS” is targeted at the investigation of generic diagnosis methodologies for different Solid Oxide Fuel Cell (SOFC) systems. The Ph.D study presented in this thesis was integrated into this project; it aims to develop a diagnostic tool for SOFC system fault detection and identification based on validated diagnostic algorithms, through applying theSOFC stack as a sensor.In this context, three algorithms, based on the k-means clustering technique, the wavelet transform and the Bayesian method, respectively, have been developed. The first algorithm serves for ex-situ diagnosis. It works on the classification of the polarization measurements of the stack, aiming to figure out the significant response variables that are able to indicate the state of health of the stack. The parameter “Silhouette” has been used to evaluate the classification solutions in order to determine the optimal number of classes/patterns to retain from the studied database.The second algorithm allows the on-line fault detection. The wavelet transform has been used to decompose the SOFC’s voltage signals for the purpose of finding out the effective feature variables that are discriminative for distinguishing the normal and abnormal operating conditions of the system. Considering the SOFC as a sensor, its reliability must be verifiedbeforehand. Thus, the feature variables are also required to be indicative to the state of health of the stack.When the stack is found being operated improperly, the actual operating parameters should be estimated so as to identify the system fault. To achieve this goal, a Bayesian network has been proposed serving as a meta-model of the stack to accomplish the estimation. At the end, the databases originated from different SOFC systems have been used to validate these three algorithms and assess their generalizability.