Les sols contaminés par les éléments traces potentiellement toxiques (PTTE) ont souvent des conséquences graves pour les écosystèmes terrestres. Plusieurs options de phytoremediaction ont été développées pour remédier les sols contaminés ; cependantl'efficacité et la capacité de ces techniques à réduire les concentrations excessives des éléments traces ou leur (phyto) disponibilité dans les sols contaminés doivent être évaluées Le présent travail s'est intéressé à étudier l'efficacité de deux options de de phytorémédiation, la phytostabilisation et la phytoextraction assistées par des amendements organiques et minéraux, à remédier les fortes concentrations de PTTE dans un sol naturel et dans un technosol contaminés. Les concentrations totales des éléments traces dissous ont été déterminées dans l'eau interstitielle du sol. L'intensité de l'exposition du sol a été évaluée par des capteurs DGT (gradient de diffusion dans les couches minces). Le phytodisponibilité des PTTE a été caractérisée par des tests de germination avec des haricots nains cultivés sur les sols contaminés pour lesquels les concentrations foliaires en éléments traces ont été déterminées. Ensuite un modèle de réseau de neurones artificiels a été appliqué pour comprendre les facteurs les plus pertinents sur la variabilité de la phytodisponibilité des PTTE. Les deux options ont étécapables de réduire les concentrations ou la phytodisponibilité des PTTE en présence des amendements. Les réseaux de neurones artificiels ont été très efficaces pour prédire les résultats manquants et pour déterminer les paramètres de contrôle de la variabilité de la phytodisponibilité des PTTE à partir des paramètres du sol.

Cette thèse porte sur le problème d'identification des systèmes non linéaires sur la base des modèles blocs. Deux types de modèles sont considérés ici, ceux de type Hammerstein et ceux de type Wiener. La plupart des solutions antérieures ont été élaborées sous des hypothèses assez contraignantes concernant le sous-système non linéaire du modèle. Celui-ci a souvent été supposé lisse (voire polynomial), inversible et sans mémoire. Les travaux présentés dans cette thèse tentent de repousser ces limites. Dans le cas des systèmes de type Hammerstein, nous élaborons un schéma d'identification qui ne fait aucune hypothèse sur l'élément non linéaire, à l'exception du fait qu'il est sans mémoire et l¥-stable ; ce schéma estime parfaitement (du moins dans la cas idéal) les paramètres du sous-système dynamique linéaire et un nuage de points de la caractéristique de l'élément non linéaire. Ce schéma est ensuite adapté au cas où l'on connaît la structure de l'élément non linéaire ; à ce propos, les non-linéarités statiques affines par morceaux et discontinues jouissent d'une attention particulière. Nous terminons la partie consacrée à l'identification des systèmes de Hammerstein, en abordant le problème des systèmes impliquant un élément non linéaire à mémoire. Deux familles d'éléments de cette nature sont considérées : celle comprenant des éléments hystérétiques non saturés et celle des éléments hystérisis-relais. Le problème est appréhendé à l'aide d'un schéma d'identification dont on établit la consistance en présence de perturbations assimilables à un bruit blanc appliqué en sortie du système. <br /> La dernière partie du mémoire est centrée sur l'identification des systèmes de Wiener, dont l'élément non linéaire n'est pas supposé inversible. A cet effet, nous présentons deux schémas d'identification de type fréquentiel et établissons leur consistance dans les mêmes conditions que précédemment concernant les perturbations. L'exigence d'excitation persistante occupe une place centrale dans cette thèse. Pour procurer cette propriété aux différents schémas d'identifications proposés, il a été fait appel à une famille de signaux d'excitation de type impulsionnelle. Dans ce cadre, un lemme technique est élaboré précisant, pour les systèmes linéaires, le lien entre cette famille de signaux et la propriété d'excitation persistante. L'adaptation de ce lemme au cas des systèmes non linéaires est illustrée dans les différents schémas d'identification.

Ce travail de thèse concerne la prédiction des structures en béton soumises à des chargements extrêmes. Il s'intéresse plus particulièrement au comportement du béton sous fort confinement où la contrainte peut atteindre des niveaux de l'ordre du giga Pascal. La modélisation de ce comportement doit être capable de reproduire la compaction irréversible. Pour ce faire, deux lois de comportement élasto-plastique - endommageable ont été développées et implantées dans un code aux éléments discrets. Les paramètres utilisés dans ces lois sont calibrés par les simulations des essais de traction/compression uniaxial, des essais hydrostatiques et triaxiaux. Une fois les paramètres calibrés, la loi montrant le meilleur agrément avec l'expérience a été choisie pour la prédiction de la réponse du béton sous différents niveaux de confinement. Les résultats du modèle sont analysés non seulement à l'échelle macroscopique mais également à l'échelle de l'élément discret. La nécessité de prendre en compte une loi d'interaction de type élasto-plastique-endommageable est aussi montrée. La deuxième partie du travail de thèse développe une méthode de couplage entre le modèle éléments discrets et un modèle d'écoulement compressible en tenant compte des mécanismes physiques fondamentaux d'interaction entre l'écoulement interne et les particules solides d'un matériau poreux. Le problème d'écoulement est résolu par une méthode en volumes finis, où le volume est discrétisé en tétraèdres issus d'une triangulation régulière de Delaunay. Notre modèle est une adaptation aux fluides compressibles d'un modèle développé initialement pour les écoulements incompressibles. Ce couplage a été utilisé pour simuler le comportement triaxial des bétons humides et saturés sous différents niveaux de confinement. Les résultats nous montrent une bonne reproduction du comportement non-drainé du béton saturé sous faible confinement. Pour fort confinement, les simulations ne se rapprochent des résultats expérimentaux qu'au prix d'une compressibilité du fluide plus faible que celle de l'eau. Par ailleurs, la contrainte effective était une variable pertinente pour décrire le comportement du béton humide par un état limite intrinsèque indépendant du degré de saturation.

La modélisation des problèmes d'interactions sol-structure et fluide-structure couvre plusieurs domaines de recherche très actifs qui traitent une multitude d'aspects tels que la géométrie non bornée du sol et dans certains cas du fluide stocké, les effets dissipatifs visqueux et radiatifs, l'application du chargement sismique, le choix des variables de base, les propriétés algébriques des systèmes d'équations résultant du couplage,... etc. Dans le présent travail, différents modèles numériques de couplage sol-structure et fluide-structure ont été examinés. Les limites de troncature géométrique du sol et du fluide on été traitées avec des éléments infinis dont les performances ont été comparées à celles des conditions de radiations. Le problème de vibrations libres couplées des systèmes fluide-structure a été résolu en introduisant de nouvelles techniques de symétrisation efficaces. De plus, une nouvelle formulation symétrique en éléments de frontière a été proposée. Cette formulation permet de produire une matrice symétrique définie positive et aboutit ainsi à un système algébrique similaire à celui qui découle de la discrétisation en éléments finis. La matrice bâtie dite "raideur équivalente" peut facilement être assemblée ou couplée avec les matrices de la formulation en éléments finis. Toutes les applications qui ont servi soit à comparer des modèles soit à valider les programmes développés, ont été effectué es dans le cas des barrages poids en béton. Ce cas constitue un problème de couplage fluide-sol-structure typique

Ce travail présente une méthode mixte expérimentale-numérique traitant des mesures mécaniques réalisées sur structures pour l'identification de propriétés matériau viscoélastiques linéaires. Les caractéristiques de ces matériaux sont complexes et dépendent de la fréquence et nous cherchons à les identifier par une approche directe et non paramétrique. Pour cela, deux méthodes d'identification adaptées au traitement des mesures de champs de réponses en fréquences sont présentées. Cellesci utilisent les données fréquentielles, sans analyse modale expérimentale. Les champs expérimentaux sont obtenus par vibromètrie laser et les champs numérique par la méthode des éléments finis. Nous utilisons un modèle volumique adapté à l'étude des structures minces et épaisses pour l'analyse numérique. Les applications se basent sur des plaques simples ou multicouches. Les résultats sont corrélés avec des mesures DMA et des identifications à partir de données modales identifiées.

Dans de nombreuses applications industrielles, les structures en vibration à dimensionner sont en contact avec un fluide (fluide autour des coques des bateaux, réservoirs, échangeurs de chaleur dans les centrales, l'industrie automobile, etc). Cependant, le comportement dynamique de la structure peut être modifié de façon importante par la présence du fluide. Le dimensionnement doit donc prendre en compte les effets de l'interaction fluide-structure.Ces applications nécessitent un couplage efficace. En outre, l'analyse dynamique des systèmes industriels est souvent coûteuse du point de vue numérique. Pour les modèles éléments finis des problèmes couplés fluide-structure, l'importance de la réduction de la taille devient évidente car les degrés de liberté du fluide seront ajoutés à ceux de la structure. Des méthodes de réduction du modèle seront utilisées pour réduire la taille des matrices obtenues.Traditionnellement, l'étude de ces systèmes couplés est fondée sur une démarche déterministe dans laquelle l'ensemble des paramètres utilisés dans le modèle prennent une valeur fixe.Par contre, il suffit d'avoir procédé à quelques expérimentations pour se rendre compte des limites d'une telle modélisation, d'où la nécessité de la prise en compte des incertitudes sur les paramètres du système couplé.Ce travail de thèse s'articule autour de trois études principales. La première consiste à mener une étude déterministe numérique et analytique des problèmes vibro-acoustiques sans réduction de modèles. Cette dernière est basée sur une formulation non symétrique déplacement/pression et une formulation symétrique déplacement/pression et potentiel des vitesses. Dans la deuxième étude, on propose deux méthodes de réduction du modèle : analyse et synthèse modales pour la résolution des problèmes vibro-acoustiques des grandes tailles des systèmes couplés modélisés par la méthode des éléments finis. La méthode de synthèse modale développée couple une méthode de sous-structuration dynamique de type Craig et Bampton et une méthode de sous domaines acoustiques.Enfin, pour tenir compte des incertitudes sur les paramètres du système couplé, on a développé dans la troisième étude une méthode numérique stochastique de synthèse modale étendue à une étude de fiabilité basée sur les approches FORM et SORM pour la résolution de ces problèmes. Ces démarches vont nous permettre de résoudre les problèmes vibro-acoustiques, sans utiliser les méthodes classiques, qui consistent à faire un calcul modal direct allié à la simulation de Monte Carlo demandant un coup de temps très élevé.Plusieurs exemples académiques et industriels ont été traités pour valider les approches proposées.L'étude numérique est conduite en utilisant un code élaboré sous MATLAB couplé au code commercial ANSYS afin d'évaluer la fiabilité du système couplé. La confrontation des résultats numériques, analytiques et expérimentaux nous permet de valider conjointement le processus de calcul et les méthodes proposées dans le domaine de l'analyse fréquentielle et l'étude fiabiliste des structures immergées. D'un point de vue industriel, ces méthodes visent à promouvoir l'introduction de la culture de l'incertain dans les métiers de la conception et encouragent la construction d'un modèle fiable et robuste pour les problèmes d'interaction fluide-structure.

Ce travail tente de développer des outils numériques efficaces pour une approche plus rationnelle et moins empirique du dimensionnement à la ruine des ouvrages de génie civil. Contrairement aux approches traditionnelles reposant sur une combinaison de calculs élastiques, l'adoption de coefficients de sécurité et une vérification locale des sections critiques, la théorie du calcul à la rupture nous semble être un outil prometteur pour une évaluation plus rigoureuse de la sécurité des ouvrages. Dans cette thèse, nous proposons de mettre en œuvre numériquement les approches statique par l'intérieur et cinématique par l'extérieur du calcul à la rupture à l'aide d'éléments finis dédiés pour des structures de plaque en flexion et de coque en interaction membrane-flexion. Le problème d'optimisation correspondant est ensuite résolu à l'aide du développement, relativement récents, de solveurs de programmation conique particulièrement efficaces. Les outils développés sont également étendus au contexte de l'homogénéisation périodique en calcul à la rupture, qui constitue un moyen performant de traiter le cas des structures présentant une forte hétérogénéité de matériaux. Des procédures numériques sont spécifiquement développées afin de déterminer puis d'utiliser dans un calcul de structure des critères de résistance homogènes équivalents. Enfin, les potentialités de l'approche par le calcul à la rupture sont illustrées sur deux exemples complexes d'ingénierie : l'étude de la stabilité au feu de panneaux en béton armé de grande hauteur ainsi que le calcul de la marquise de la gare d'Austerlitz / This work aims at developping efficient numerical tools for a more rational and less empirical assessment of civil engineering structures yield design. As opposed to traditionnal methodologies relying on combinations of elastic computations, safety coefficients and local checking of critical members, the yield design theory seems to be a very promising tool for a more rigourous evaluation of structural safety. Lower bound static and upper bound kinematic approaches of the yield design theory are performed numerically using dedicated finite elements for plates in bending and shells in membrane-bending interaction. Corresponding optimization problems are then solved using very efficient conic programming solvers. The proposed tools are also extended to the framework of periodic homogenization in yield design, which enables to tackle the case of strong material heterogeneities. Numerical procedures are specifically tailored to compute equivalent homogeneous strength criteria and to use them, in a second step, in a computation at the structural level. Finally, the potentialities of the yield design approach are illustrated on two complex engineering problems : the stability assessment of high-rise reinforced concrete panels in fire conditions and the computation of the Paris-Austerlitz railway station canopy

Les structures en béton et béton armé de grandes dimensions, en particulier les enceintes de confinement, peuvent être sujettes à de la fissuration localisée suite à leur vieillissement ou dans le cas d’une forte sollicitation (APRP, par exemple). Afin d’optimiser les actions de maintenance, il est nécessaire de disposer d’un modèle prédictif de l’endommagement du béton. Ce phénomène se produit à une échelle matériau relativement petite et un modèle prédictif nécessite un maillage fin et une loi de comportement non linéaire. Hors ce type de modélisation ne peut être directement appliquée sur une structure de génie civil de grande échelle, le calcul étant trop lourd pour les machines actuelles.Une méthode de calcul est proposée, qui concentre l’effort de calcul sur les zones d’intérêt (parties endommagées) de la structure en éliminant les zones non endommagées. L’objectif est ainsi d’utiliser la puissance de calcul disponible pour la caractérisation des propriétés des fissures notamment. Cette approche utilise la méthode de condensation statique de Guyan pour ramener les zones élastiques à un ensemble de conditions aux limites appliquées aux bornes des zones d’intérêt. Lorsque le système évolue, un système de critères permet de promouvoir à la volée des zones élastiques en zones d’intérêt si de l’endommagement y apparaît. Cette méthode de condensation adaptative permet de réduire la dimension du problème non linéaire sans altérer la qualité des résultats par rapport à un calcul complet de référence. Cependant, une modélisation classique ne permet pas de prendre en compte les divers aléas impactant le comportement de la structure : propriétés mécaniques, géométrie, chargement… Afin de mieux caractériser ce comportement en tenant compte des incertitudes, la méthode de condensation adaptative proposée est couplée avec une approche de collocation stochastique. Chaque calcul déterministe nécessaire pour caractériser les incertitudes sur les grandeurs d’intérêt de la structure est ainsi réduit et les étapes de prétraitement nécessaires à la condensation sont elles-mêmes mutualisées via une deuxième collocation. L’approche proposée permet ainsi de produire pour un coût de calcul limité des densités de probabilités des grandeurs d’intérêt d’une grande structure. Les stratégies de résolution proposées rendent accessibles à l’échelle locale une modélisation plus fine que celle qui pourrait s’appliquer sur l’ensemble de la structure. Afin de bénéficier d’une meilleure représentativité à cette échelle, il est nécessaire de représenter les effets tridimensionnels des hétérogénéités. Dans le domaine du génie civil et nucléaire, cela concerne au premier chef les câbles de précontrainte, traditionnellement représentés en unidimensionnel. Une approche est donc proposée, qui s’appuie sur un maillage et une modélisation 1D pour reconstruire un volume équivalent au câble et retransmettre les efforts et rigidités dans le volume de béton. Elle combine la représentativité d’un modèle 3D complet et conforme des câbles lorsque le maillage s’affine et la facilité d’utilisation et paramétrage d’un modèle 1D. L’applicabilité des méthodes proposées à une structure de génie civil de grande échelle est évaluée sur un modèle numérique d’une maquette à l’échelle 1/3 de l’enceinte de confinement interne d’un réacteur de type REP 1300 MWe à double paroi. / Large-scale concrete and reinforced concrete structures, and in particular containment buildings, may undergo localized cracking when they age or endure strong loadings (LOCA for instance). In order to optimize the maintenance actions, a predictive model of concrete damage is required. This phenomenon takes place at a rather small material scale and a predictive model requires a refined mesh and a nonlinear constitutive law. This type of modelling cannot be applied directly on a large-scale civil engineering structure, as the computational load would be too heavy for the existing machines.A simulation method is proposed to focus the computational effort on the areas of interest (damaged parts) of the structure while eliminating the undamaged areas. It aims at using the available computing power for the characterization of crack properties in particular. This approach uses Guyan’s static condensation technique to reduce the elastic areas to a set of boundary conditions applied to the areas of interest. When the system evolves, a set of criteria allows to promote on the fly the elastic areas to areas of interest if damage appears. This adaptive condensation technique allows to reduce the dimension of a nonlinear problem without degrading the quality of the results when compared to a full reference simulation.However, a classical modelling does not allow to take into account the various unknowns which will impact the structural behaviour: mechanical properties, geometry, loading… In order to better characterize this behaviour while taking into account the various uncertainties, the proposed adaptive condensation method is coupled with a stochastic collocation approach. Each deterministic simulation required for the characterization of the uncertainties on the structural quantities of interest is therefore reduced and the pre-processing steps necessary to the condensation technique are also reduced using a second collocation. The proposed approach allows to produce for a reduced computational cost the probability density functions of the quantities of interest of a large structure.The proposed calculation strategies give access at the local scale to a modelling finer than what would be applicable to the full structure. In order to improve the representativeness at this scale, the tridimensional effects of the heterogeneities must be taken into account. In the civil and nuclear engineering field, one of the main issues is the modelling of prestressing tendons, usually modelled in one dimension. A new approach is proposed, which uses a 1D mesh and model to build a volume equivalent to the tendon and redistribute the forces and stiffnesses in the concrete. It combines the representativeness of a full conform 3D modelling of the tendon when the mesh is refined and the ease of use of the 1D approaches.The applicability of the proposed methodologies to a large-scale civil engineering structure is evaluated using a numerical model of a 1/3 mock-up of a double-wall containment building of a PWR 1300 MWe nuclear reactor.

Cette thèse s’intéresse à la déformabilité des renforts à bases de fibres de lin lors de leur mise en forme, au cours du procédé RTM (Resin Tranfer Molding), notamment lorsque les formes considérées sont complexes. La compréhension des défauts et mécanismes de déformations de ces renforts est alors nécessaire. S’appuyant sur un banc d’essai d’emboutissage spécifique, une caractérisation globale (défauts de mise en forme) et locale (déformations des mèches par méthode optique) des préformes a été réalisée, pour différents paramètres procédés (pressions serre-flan et orientation du renfort) et matériaux (architecture). De cette étude, un défaut spécifique lors de la mise en forme a été mis en évidence : le bouclage, sous la forme de flambement hors plan des mèches. Afin de contribuer à sa compréhension, peu décrite dans la littérature, un banc de caractérisation, indépendamment du procédé, a été conçu pour reproduire le phénomène dans le plan et identifier les différents mécanismes à l’origine de ce défaut. Une méthode interférométrique pour quantifier ce phénomène est proposée et une étude paramétrique est réalisée (tensions appliquées, angle de courbure, orientation, architecture du renfort). Au sein du renfort, le comportement des mèches présente de fortes non-linéarités en tension. Un modèle macroscopique par éléments finis est proposé afin de mettre en avant l’importance de ces non-linéarités sur les résultats numérique issus de l’opération de formage. Une étude de sensibilité sous différentes conditions procédés de ces non-linéarités est également apportée. / This study concerns the manufacturing process of composite material from woven flax reinforcement and specifically the preforming stage of the RTM (Resin Transfer Molding) process, with complex geometries. During the process several deformation modes take place and several defects may appear and it is therefore important to understand the mechanisms controlling their appearance. A specific sheet forming bench has been used to characterise the formed shapes at the macroscopic scale (global defects) as well as at the mesoscopic scale (tow strains). Various process parameters (blank holder pressures and reinforcement orientation) and the local tensile strains were measured using an optical method. From the observed defects, the tow buckling, out of plane bending of tows, was particularly investigated as this one is not very much described in the litterature. An experimental setup independent of the forming process was designed to specifically study this defect and to quantify the key parameters controlling its appearance. An interferometric method was proposed to measure the size of the buckles with accuracy and a parametric study (applied tensions, bending angles, orientation and architecture of the reinforcement) has been completed. Within the reinforcement, the tensile behaviour of tows presents strong non-linearities. A finite elements model at the macroscopic scale is proposed to highlight the impact of these non-linearities on the forming modeling results. A sensitivity study under several process conditions has been realized.

Les impacts sur structures composites peuvent fortement diminuer leur résistance résiduelle sans laisser de marque visible sur la surface extérieure. Dans le domaine aéronautique, un seuil minimum de détection de l’endommagement d’impact est défini, basé sur l’indentation permanente laissée par l’impact. En deçà de ce seuil, la structure doit résister à un chargement défini : c’est la notion de tolérance aux dommages d’impact. Dimensionner numériquement une structure composite en tenant compte des aspects détectabilité et tolérance aux dommages nécessite donc de savoir modéliser à la fois l’impact, l’indentation permanente et la résistance résiduelle sous compression.Ces travaux se focalisent sur la modélisation numérique des composites stratifiés formés de plis unidirectionnels. L’objectif est d’établir un modèle prédictif de la tenue résiduelle après impact. Une étude expérimentale a été réalisée afin d’étudier le scénario d’endommagement à l’impact et sous compression après impact (CAI), et de fournir des résultats expérimentaux pour valider les modèles numériques.Une modélisation par éléments finis avec une approche de type Discrete Ply Model (DPM) est effectuée, basée sur des travaux précédents. Le modèle d'impact est amélioré et validé sur différentes séquences d'empilement pour assurer la robustesse du modèle. Des essais de flexion trois points spécifiques sont réalisés pour apporter une meilleure compréhension de la formation de l'indentation permanente. Un nouveau modèle d'indentation permanente est alors proposé et appliqué dans le modèle d'impact. Enfin, un modèle de CAI est construit pour prédire la résistance résiduelle. Les trois étapes : impact, indentation et CAI sont combinées au sein d’un unique modèle. / Impacts on composite structures can greatly reduce their residual strength without leaving a visible mark on the outer surface. In aeronautics, a minimum detection threshold of the impact damage is defined, based on the permanent indentation left after impact. Below this threshold, the structure must withstand a defined load: it is the notion of impact damage tolerance. The numerical design of a composite structure taking into account aspects such as detectability and damage tolerance thus requires to know how to represent impact, permanent indentation and residual strength under compression.This work focuses on the numerical modeling of composite laminates made of unidirectional plies. The objective is to develop a predictive model of post-impact residual strength. An experimental study was conducted to investigate the damage scenario during impact and compression after impact (CAI), and provide experimental data to validate the simulations.A finite element modeling with a Discrete Ply Model (DPM) approach is performed based on previous work done at the laboratory. The impact model is improved and validated on different stacking sequences to ensure the robustness of the model. Specific three-point bending tests are performed to have a better understanding of the formation of permanent indentation. A new model of permanent indentation is then proposed and applied in the impact model. Finally, a model is built to predict CAI residual strength. The three steps: impact, indentation and CAI are combined into a single model.