L'objectif principal de cette thèse concerne la mise en oeuvre d'une méthoded'éléments finis stabilisés pour la simulation des plasmas de fusion. Pour cela,nous avons d'abord dérivé les modèles magnétohydrodynamiques depuis lemodèle cinétique. Les modèles MHD sont généralement utilisés pour simuler lesinstabilités macroscopiques des plasmas. Nous nous sommes concentrés sur lemodèles de la MHD complète. Ensuite, l'approche numérique est décrite dans lecadre de la stabilisation Variationelle Multi-Échelles (VMS). Cette stabilisationvient ajouter un terme à la formulation faible pour mimer les effets des échellesnon-résolues sur celles résolues. Si les effets de ces sous-échelles ne sont paspris en compte lorsque l'on traite des écoulements dominés par convection,comme dans le cadre des plasmas de fusion, le schéma numérique conduit àdes résultats non-physiques. Une étude détaillée de l'instabilité de « Kinkinterne » a été faite ainsi qu'une étude préliminaire des plasmas avec point-Xayant pour but la validation du schéma numérique développé ici / The main objective of this thesis concerns the implementation of a robuststabilized finite element method for simulating fusion plasmas. For that, we firstderive the magnetohydrodynamic models from the kinetic model. MHD modelsare generally used for macroscopic simulations of plasma instabilities. Weconcentrate ou efforts on the full MHD model. Next, the numerical approach isdescribed in the context of the Variational Multi-Scale (VMS) stabilization. Thisstabilization comes to add a term to the weak formulation to mimics the effectsof the unresolved scales over the coarse scales. If the effects of these subscalesare not taken into account when dealing with fluxes dominated byconvection, as it is the cases for fusion plasmas, the numerical scheme canlead to unphysical results. A detailed study of the resistive internal kinkinstability has been done as well as an introductory study of the so called Xpointplasmas in order to validate the numerical scheme developed here

Simulations of particle systems at constant temperature may be used to estimate several of the system’s physical properties, and some require integration over very long time to be accurate. To achieve sufficient accuracy in finite time the choice of numerical scheme is important and we suggest to use the symplectic Euler method combined with a step in an Ornstein-Uhlenbeck process. This scheme is computationally very cheap and is often used in applications of molecular dynamics. This thesis strives to motivate the usage of the scheme due to the lack of theoretical results and comparisons to alternative methods. We conduct three numerical experiments to evaluate the scheme. The design of each experiment aims to expose weaknesses or strengths of the method. For both model problems and more realistic experiments are the results positive in favor of the method; the symplectic Euler method combined with an Ornstein- Uhlenbeck step does perform well over long times. / Simuleringar av partikelsystem vid konstant temperatur kan användas för att uppskatta flera av systemets fysiska egenskaper. Vissa klasser av egenskaper kräver integration över väldigt lång tid för att uppnå hög noggrannhet och för att uppnå detta i ändlig tid är valet av numerisk metod viktigt. Vi föreslår att använda den symplektiska Euler-metoden i kombination med ett implicit steg i en Ornstein-Uhlenbeck-process. Detta stegschema kräver låg beräkning jämfört med andra scheman och används redan i olika applikationer av molekyldynamik. Detta examensarbete eftersträvar att än mer motivera användandet av schemat, eftersom teoretiska resultat som stödjer metoder är få, och avsaknaden av tidigare liknande studier är betydlig. Vi genomför tre numeriska experiment för att pröva schemat. Under utformningen av experimenten har vi försökt att inkorporera olika fenomen som kan orsaka svårigheter för metoden för att exponera svagheter eller styrkor hos den. För båda modellproblem och för ett mer realistiskt experiment är resultaten positiva till schemats fördel; metoden att kombinera ett symplektisk Euler-steg med ett steg i Ornstein-Uhlenbeck-processen presterar bra över lång tid.

This study is devoted to the thermo-mechanical constitutive modeling for saturated stiff clays and the development of a corresponding efficient stress integration algorithm. The mechanical behavior of natural Boom Clay in isothermal conditions was first characterized. The Modified Cam Clay model (MCC) was then applied to simulate the natural Boom Clay behavior. It has been found that the MCC gives poor-quality predictions of the natural Boom Clay behavior. Thereby, an adapted Cam Clay model (ACC-2) was developed by introducing a new yield surface and a new plastic potential as well as a Two-surface plastic mechanism. This model allows satisfactory prediction of the main features of the mechanical behavior of natural Boom Clay. Moreover, the constitutive equations of this model can be formulated mathematically as in a classic elasto-plastic model. Thus, the classic stress integration algorithm can be applied. The thermal effects were considered by assessing the performance of some advanced thermo-mechanical models (Cui et al., 2000; Abuel-Naga et al., 2007; Laloui and François, 2008; 2009). It appears that all the three models can capture the main features of the thermo-mechanical behavior of saturated clays. However, each constitutive model has its own limitations or unclear points from the theoretical point of view. The stress integration algorithm of the thermo-mechanical model proposed by Cui et al. (2000) at the stress point level was also developed using a specifically designed adaptive time-stepping scheme. The computation time required to achieve a given accuracy is largely reduced with the adaptive sub-stepping considered for both mechanical and thermal loadings. A Two-surface thermo-mechanical model (TEAM model) was developed based on the Two-surface plastic mechanism. The proposed model extends the model of Cui et al. (2000) to a Two-surface formulation, considering the plastic strain coupling between the thermal and the mechanical loading paths. The simulation of drained tests shows that this model can capture the main thermo-mechanical features of natural Boom Clay along different loading paths. The TEAM model was finally extended to undrained conditions. After setting up an appropriate effective stress principle and defining a volumetric strain condition, the undrained heating process was analyzed. The validity of the thermo-hydro-mechanical constitutive equations was examined based on the data from typical tests

Les systèmes dynamiques qui évoluent sur un grand intervalle de temps (dynamique moléculaire, prédiction astronomique, turbulence...) occupent une place importante dans le domaine de la science de l'ingénieur. Leur résolution numérique constitue, jusqu'à l'heure actuelle, un défi. En effet, la simulation de la solution nécessite un solveur non seulement rapide mais aussi qui respecte les propriétés physiques du problème, pour garantir la stabilité. Dans cette thèse, on se propose d'étudier, vis-à-vis de cette problématique, un schéma d'intégration temporelle basée sur la décomposition de la solution en série temporelle, suivie de la technique de resommation de Borel des séries divergentes. On analyse alors la rapidité du schéma sur des problèmes modèles. Ensuite, on montre sa capacité à préserver la structure des équations (symplecticité, iso-spectralité, conservation de l'énergie...) à un ordre arbitrairement élevé. Par la suite, on applique le schéma à la résolution d'équations aux dérivées partielles issues de la mécanique, dont les équations de la chaleur, de Burgers et de Navier-Stokes bidimensionnelles. Pour cela, on associe le schéma à une méthode de discrétisation par éléments finis en espace. Enfin, dans le but de rendre l'algorithme plus robuste, on s'intéresse à la représentation de la somme de Borel par une série de factorielle généralisée. / Dynamical systems which evolve in a large time interval (molecular dynamic, astronomical prediction, turbulence…) take an important place in engineering science. Their numerical resolution has so far constituted a challenge. Indeed, the simulation of the solution requires a solver which is not only fast but also respects the physical properties of the problem, to ensure the stability. In this thesis, we propose to study, regarding this issue, a time integration scheme based on the decomposition of the solution into time series, followed by Borel's resummation technique of divergent series. We analyse the speed of scheme on model problems. Next, we show its capability to preserve the structure of the equation (symplecticity, iso-spectrality, conservation of energy…) up to an arbitrary high order. Thereafter, we use the scheme to resolve partial differential equations coming from mechanics, including the two-dimensional heat equation, Burger’s equation and the Navier-Stokes equation. To this aim, we choose a finite element method for space discretisation. Finally, and in order to make the algorithm more robust, we are interested in the representation of the Borel sum by a generalized factorials series.

Cette thèse traite des équations de Maxwell en domaine temporel. Le principal objectif est de proposer des méthodes de type éléments finis d'ordre élevé pour les équations de Maxwell et des schémas d'intégration en temps efficaces sur des maillages localement raffinés. Nous considérons des méthodes GDDT (Galerkine Discontinues en Domaine Temporel) s'appuyant sur une interpolation polynomiale d'ordre arbitrairement élevé des composantes du champ électromagnétique. Les méthodes GDDT pour les équations de Maxwell s'appuient le plus souvent sur des schémas d'intégration en temps explicites dont la condition de stabilité peut être très restrictive pour des maillages raffinés. Pour surmonter cette limitation, nous considérons des schémas en temps qui consistent à appliquer un schéma implicite localement, dans les régions raffinées, tout en préservant un schéma explicite sur le reste du maillage. Nous présentons une étude théorique complète et une comparaison de deux méthodes GDDT localement implicites. Des expériences numériques en 2D et 3D illustrent l'utilité des schémas proposés. Le traitement numérique de milieux de propagation complexes est également l'un des objectifs. Nous considérons l'interaction des ondes électromagnétiques avec les tissus biologiques qui est au cœur de nombreuses applications dans le domaine biomédical. La modélisation numérique nécessite alors de résoudre le système de Maxwell avec des modèles appropriés de dispersion. Nous formulons une méthode GDDT localement implicite pour le modèle de Debye et proposons une analyse théorique et numérique complète du schéma. / This work deals with the time-domain formulation of Maxwell's equations. The main objective is to propose high-order finite element type methods for the discretization of Maxwell's equations and efficient time integration methods on locally refined meshes. We consider Discontinuous Galerkin Time-Domain (DGTD) methods relying on an arbitrary high-order polynomial interpolation of the components of the electromagnetic field. Existing DGTD methods for Maxwell's equations often rely on explicit time integration schemes and are constrained by a stability condition that can be very restrictive on highly refined meshes. To overcome this limitation, we consider time integration schemes that consist in applying an implicit scheme locally i.e. in the refined regions of the mesh, while preserving an explicit scheme in the complementary part. We present a full theoretical study and a comparison of two locally implicit DGTD methods. Numerical experiments for 2D and 3D problems illustrate the usefulness of the proposed time integration schemes. The numerical treatment of complex propagation media is also one of the objectives. We consider the interaction of electromagnetic waves with biological tissues that is of interest to applications in biomedical domain. Numerical modeling then requires to solve the system of Maxwell's equations coupled to appropriate models of physical dispersion. We derive a locally implicit DGTD method for the Debye model and we achieve a full theoretical and numerical analysis of the resulting scheme.

Indiana University-Purdue University Indianapolis (IUPUI) / In this thesis, an explicit time integration code which integrates multibody dynamics and the discrete element method is used for modelling the excavation and moving operation of cohesive soft soil (such as mud and snow) by bulldozers. A soft cohesive soil material model (that includes normal and tangential inter-particle force models) is used that can account for soil compressibility, plasticity, fracture, friction, viscosity and gain in cohesive strength due to compression. In addition, a time relaxation sub-model for the soil plastic deformation and cohesive strength is added in order to account for loss in soil cohesive strength and reduced bulk density due to tension or removal of the compression. This is essential in earth moving applications since the soil that is dug typically becomes loose soil that has lower shear strength and lower bulk density (larger volume) than compacted soil. If the model does not account for loss of soil shear strength then the dug soil pile in front of the blade of a bulldozer will have an artificially high shear strength. A penalty technique is used to impose joint and normal contact constraints. An asperity-based friction model is used to model contact and joint friction. A Cartesian Eulerian grid contact search algorithm is used to allow fast contact detection between particles. A recursive bounding box contact search algorithm is used to allow fast contact detection between the particles and polygonal contact surfaces. A multibody dynamics bulldozer model is created which includes the chassis/body, C-frame, blade, wheels and hydraulic actuators. The components are modelled as rigid bodies and are connected using revolute and prismatic joints. Rotary actuators along with PD (Proportional-Derivative) controllers are used to drive the wheels. Linear actuators along with PD controllers are used to drive the hydraulic actuators. Polygonal contact surfaces are defined for the tires and blade to model the interaction between the soil and the bulldozer. Simulations of a bulldozer performing typical shallow digging operations in a cohesive soil are presented. The simulation of a rear wheel drive bulldozer shows that, it has a limited digging capacity compared to the 4-wheel drive bulldozer. The effect of the relaxation parameter can be easily observed from the variation in the Bulldozer's velocity. The higher the relaxation parameter, the higher is the bulldozer's velocity while it is crossing over the soil patch. For the low penetration depth run the bulldozer takes less time compared to high penetration depth. Also higher magnitudes of torques at front and rear wheels can be observed in case of high penetration depth. The model is used to predict the wheel torque, wheel speed, vehicle speed and actuator forces during shallow digging operations on three types of soils and at two blade penetration depths. The model presented can be used to predict the motion, loads and required actuators forces and to improve the design of the various bulldozer components such as the blade, tires, engine and hydraulic actuators.

Der Vortrag beschreibt wie sich mittels der unterschiedlichen Berechnungsverfahren zur Lösung dynamischer Strukturpobleme der Einschlag eines idealisierten Bruchstücks in eine Schutzwand berechnen lässt. Dies wird mittels zweier kommerzieller FEM-Programme beschrieben: a.) Creo Simulate nutzt zur Lösung die Methode der modalen Superposition, d.h., es können nur lineare dynamische Systeme mit rein modaler Dämpfung berechnet werden. Kontakt zwischen zwei Bauteilen lässt sich damit nicht erfassen. Die unbekannte Kraft-Zeit-Funktion des Einschlagvorganges muss also geeignet abgeschätzt und als äußere Last auf die Schutzwand aufgebracht werden. Je dynamischer der Einschlagvorgang, desto eher wird der Gültigkeitsbereich des zugrunde liegenden linearen Modells verlassen. b.) Abaqus/Explicit nutzt ein direktes Zeitintegrationsverfahren zur schrittweisen Lösung der zugrunde liegenden Differentialgleichung, die keine tangentiale Steifigkeitsmatrix benötigt. Damit können sowohl Materialnichtlinearitäten als auch Kontakt geeignet erfasst und damit die Kraft-Zeit-Funktion des Einschlages ermittelt werden. Auch bei extrem hochdynamischen Vorgängen liefert diese Methode ein gutes Ergebnis. Es müssen dafür jedoch weit mehr Werkstoffdaten bekannt sein, um das nichtlineare elasto-plastische Materialverhalten mit Schädigungseffekten korrekt zu beschreiben. Die Schwierigkeiten der Werkstoffdatenbestimmung werden in den Grundlagen erläutert. / The presentation describes how to analyze the impact of an idealized fragment into a stell protective panel with different dynamic analysis methods. Two different commercial Finite Element codes are used for this: a.) Creo Simulate: This code uses the method of modal superposition for analyzing the dynamic response of linear dynamic systems. Therefore, only modal damping and no contact can be used. The unknown force-vs.-time curve of the impact event cannot be computed, but must be assumed and applied as external force to the steel protective panel. As more dynamic the impact, as sooner the range of validity of the underlying linear model is left. b.) Abaqus/Explicit: This code uses a direct integration method for an incremental (step by step) solution of the underlying differential equation, which does not need a tangential stiffness matrix. In this way, matieral nonlinearities as well as contact can be obtained as one result of the FEM analysis. Even for extremely high-dynamic impacts, good results can be obtained. But, the nonlinear elasto-plastic material behavior with damage initiation and damage evolution must be characterized with a lot of effort. The principal difficulties of the material characterization are described.

Der Vortrag beschreibt wie sich mittels der unterschiedlichen Berechnungsverfahren zur Lösung dynamischer Strukturpobleme der Einschlag eines idealisierten Bruchstücks in eine Schutzwand berechnen lässt. Dies wird mittels zweier kommerzieller FEM-Programme beschrieben: a.) Creo Simulate nutzt zur Lösung die Methode der modalen Superposition, d.h., es können nur lineare dynamische Systeme mit rein modaler Dämpfung berechnet werden. Kontakt zwischen zwei Bauteilen lässt sich damit nicht erfassen. Die unbekannte Kraft-Zeit-Funktion des Einschlagvorganges muss also geeignet abgeschätzt und als äußere Last auf die Schutzwand aufgebracht werden. Je dynamischer der Einschlagvorgang, desto eher wird der Gültigkeitsbereich des zugrunde liegenden linearen Modells verlassen. b.) Abaqus/Explicit nutzt ein direktes Zeitintegrationsverfahren zur schrittweisen Lösung der zugrunde liegenden Differentialgleichung, die keine tangentiale Steifigkeitsmatrix benötigt. Damit können sowohl Materialnichtlinearitäten als auch Kontakt geeignet erfasst und damit die Kraft-Zeit-Funktion des Einschlages ermittelt werden. Auch bei extrem hochdynamischen Vorgängen liefert diese Methode ein gutes Ergebnis. Es müssen dafür jedoch weit mehr Werkstoffdaten bekannt sein, um das nichtlineare elasto-plastische Materialverhalten mit Schädigungseffekten korrekt zu beschreiben. Die Schwierigkeiten der Werkstoffdatenbestimmung werden in den Grundlagen erläutert. / The presentation describes how to analyze the impact of an idealized fragment into a stell protective panel with different dynamic analysis methods. Two different commercial Finite Element codes are used for this: a.) Creo Simulate: This code uses the method of modal superposition for analyzing the dynamic response of linear dynamic systems. Therefore, only modal damping and no contact can be used. The unknown force-vs.-time curve of the impact event cannot be computed, but must be assumed and applied as external force to the steel protective panel. As more dynamic the impact, as sooner the range of validity of the underlying linear model is left. b.) Abaqus/Explicit: This code uses a direct integration method for an incremental (step by step) solution of the underlying differential equation, which does not need a tangential stiffness matrix. In this way, matieral nonlinearities as well as contact can be obtained as one result of the FEM analysis. Even for extremely high-dynamic impacts, good results can be obtained. But, the nonlinear elasto-plastic material behavior with damage initiation and damage evolution must be characterized with a lot of effort. The principal difficulties of the material characterization are described.

Fires following earthquake (FFE) have historically produced enormous post-earthquake damage and losses in terms of lives, buildings and economic costs, like the San Francisco earthquake (1906), the Kobe earthquake (1995), the Turkey earthquake (2011), the Tohoku earthquake (2011) and the Christchurch earthquakes (2011). The structural fire performance can worsen significantly because the fire acts on a structure damaged by the seismic event. On these premises, the purpose of this work is the investigation of the experimental and numerical response of structural and non-structural components of steel structures subjected to fire following earthquake (FFE) to increase the knowledge and provide a robust framework for hybrid fire testing and hybrid fire following earthquake testing. A partitioned algorithm to test a real case study with substructuring techniques was developed. The framework is developed in MATLAB and it is also based on the implementation of nonlinear finite elements to model the effects of earthquake forces and post-earthquake effects such as fire and thermal loads on structures. These elements should be able to capture geometrical and mechanical non-linearities to deal with large displacements. Two numerical validation procedures of the partitioned algorithm simulating two virtual hybrid fire testing and one virtual hybrid seismic testing were carried out. Two sets of experimental tests in two different laboratories were performed to provide valuable data for the calibration and comparison of numerical finite element case studies reproducing the conditions used in the tests. Another goal of this thesis is to develop a fire following earthquake numerical framework based on a modified version of the OpenSees software and several scripts developed in MATLAB to perform probabilistic analyses of structures subjected to FFE. A new material class, namely SteelFFEThermal, was implemented to simulate the steel behaviour subjected to FFE events.