No description available.

Le contrôle de la sûreté des ouvrages hydrauliques est l’une des grandes priorité dans le domaine du génie civil et de l’ingénierie hydraulique. Durant sa vie, un ouvrage est soumis à des sollicitations variables hydromécaniques, physicochimiques et climatiques qui contribuent à son éventuelle détérioration. Parmi les phénomènes qui en résultent, l’érosion des sols sous toutes ses formes représente un enjeu majeur à comprendre, maîtriser et empêcher. L’objectif de ce travail est d’étudier l’érosion des sols par le Jet Erosion Test. La première partie est consacrée à la description des dispositifs expérimentaux, surtout le Jet Erosion Test (JET) développé à l’Ecole Centrale Paris permettant de mesurer directement quelques paramètres d’érosion. A partir des résultats du JET, à l’aide d’une loi d’érosion empirique, nous déduisons la contrainte de cisaillement critique, le coefficient d’érosion, la profondeur d’érosion d’équilibre. La deuxième partie du travail est consacrée à l’étude de l’influence des paramètres de compactage sur l’infiltration de l’eau et la résistance du sol, en utilisant le pénétromètre. Dans les troisième et quatrième parties, nous étudions l’influence des propriétés géotechniques du sol et celle des paramètres d’essai sur les paramètres d’érosion du sol. Les résultats obtenus montrent que les paramètres d’érosion sont influencés non seulement par les propriétés géotechniques du sol mais aussi par les paramètres d’essai. La dernière partie présente la synthèse entre les résultats des essais de pénétromètre et ceux des essai de JET, et tente de relier les paramètres d’érosion avec les propriétés géotechniques du sol. / Control of the safety of hydraulic structures is a major priority in the field of civil and hydraulic engineering. During its life, the hydraulic structure is submitted to variable hydromechanical, physicochemical and climatic loads that may contribute to its possible failure. Among the resulting phenomena, soil erosion under all its forms is a major challenge which it is important to understand, control and prevent. The objective of this work is to study soil erosion by the Jet Erosion Test. The first part is devoted to the description of the experimental devices, especially the Jet Erosion Test (JET) developed at the Ecole Centrale Paris to directly measure some erosion parameters. From the results of JET, using an empirical erosion law, we deduce the critical shear stress, the erosion coefficient, the equilibrium scour depth. The second part of this work is devoted to the study of the influence of compaction parameters on water infiltration and soil strength, using the penetrometer. In the third and fourth parts, we study the influence of the geotechnical properties of soil and of the test parameters on the erosion parameters of soil. The obtained results show that the erosion parameters are influenced not only by the geotechnical properties of soil but also by the test parameters. The final section presents a synthesis of the results of penetrometer tests and JET tests, and attempts to link the erosion parameters with the geotechnical properties of soil.

Dans le but d’améliorer les transistors bipolaires TBH SiGe, nous proposons d’étudier l’impact de la contrainte mécanique sur leurs performances. En effet, cette contrainte permet de libérer un degré de liberté supplémentaire pour améliorer les propriétés du transport grâce à un changement de la structure de bande d’énergie du semiconducteur. Ainsi, nous avons proposé de nouvelles architectures de composants basées sur l’ingénierie de la contrainte mécanique dans les semiconducteurs. Deux approches ont été utilisées dans cette étude pour générer la tension mécanique adéquate à l'intérieur du dispositif. La première approche consiste à appliquer une contrainte mécanique sur la base du transistor en utilisant une couche de SiGe extrinsèque. La seconde approche vise à appliquer une contrainte dans la région du collecteur en utilisant une couche contrainte. Les résultats obtenus montrent que cette méthode peut être une approche prometteuse pour améliorer les performances des TBH. / The Impact of strain engineering technology applied on Si BJT/SiGe HBT devices on the electrical properties and frequency response has been investigated. Strain technology can be used as an additional degree of freedom to enhance the carriers transport properties due to band structure changes and mobility enhancement. New concepts and novel device architectures that are based on strain engineering technology have been explored using TCAD modeling. Two approaches have been used in this study to generate the proper mechanical strain inside the device. The first approach was through introducing strain at the device’s base region using SiGe extrinsic stress layer. The second approach was through introducing strain at the device’s collector region using strain layers. The obtained results obviously show that strain engineering technology principle applied to BJT/HBT device can be a promising approach for further devices performance improvements.

La mécanoluminescence (ML) est un phénomène physique qui se traduit par l’émission de lumière suite à l’application d’une contrainte mécanique. Dans le cadre de ce travail de thèse, les matériaux mécanoluminescents vont être utilisés afin de visualiser les contraintes qui peuvent être générées dans un verre lors d’une sollicitation mécanique. Nous avons élaboré un composite constitué de particules mécanoluminescentes SrAl2O4:Eu, Dy incorporées dans une matrice vitreuse. Nous nous sommes ensuite intéressés à l’origine et au mécanisme de ce phénomène en réalisant une étude mécanique sur le composite. Ces essais nous ont permis de développer un modèle empirique qui relie la contrainte à l’intensité émise. La structure cristalline des particules jouant un rôle clef dans la ML, nous avons également mené une étude théorique de cette structure, ce qui nous a permis de corréler les observations macroscopiques faites sur le composite aux phénomènes ayant lieu à une échelle atomique. / Mechanoluminescence (ML) is the phenomenon of light emission generated by a mechanical loading. Mechanoluminescent materials were used to visualize the stress in glasses during mechanical loadings. A mechanoluminescent particulate glass composite was synthesized by incorporating SrAl2O4:Eu, Dy particles in a glass matrix for the first time. We performed a mechanical study on the composite to determine the origin and the mechanism of this phenomenon. We developed a theoretical model which accurately predicts the ML intensity changes induced by a complex mechanical loading. The crystal structure of the particles has a key role in the ML and we carried out a theoretical study of this structure to correlate the macroscopic observations made on the composite to the phenomena which occur at the atomic scale.

Les verres métalliques massifs à base de zirconium, développés depuis la fin des années 80, présentent des propriétés mécaniques très intéressantes. Ils peuvent être envisagés pour de multiples applications y compris à des températures élevées et il est donc intéressant d'étudier leur résistance à l'oxydation dans de telles conditions. L'objectif fondamental de cette étude consiste à mieux comprendre le rôle de divers paramètre thermodynamiques et chimiques sur le comportement à l'oxydation des verres métalliques à base de zirconium à des températures intermédiaires sous air sec, à déterminer les contraintes résiduelles au sein de la couche d'oxyde formée, en comparaison avec des alliages amorphe cristallisés après un traitement de recuit. La cinétique d'oxydation de ces verres et la structure cristalline de la couche d'oxyde ZrO2 dépend de la température et de la durée d'oxydation : pour des durées courtes d'oxydation et pour une température légèrement inférieure à Tg, la cinétique d'évolution est parabolique alors que si l'échantillon est oxydé à T > Tg, la loi cinétique peut être décomposée en deux parties. Les mêmes alliages cristallisés après un traitement de recuit, s'oxydent selon une loi parabolique quelque soit la température. Pour des durées d'oxydation longues à une température proche de Tg, les lois cinétiques deviennent plus complexes et la cristallisation du verre métallique pouvant avoir lieu au cours des essais d'oxydation. De même la structure cristalline des couches d'oxyde dépend de la température d'oxydation. Ainsi, pour T < Tg, la couche d'oxyde des alliages amorphes est uniquement constituée de la zircone tétragonale alors que pour ces mêmes alliages portés à une température T > Tg et pour les alliages cristallisés, on observe de la zircone monoclinique. Les mécanismes d'oxydation dépendent de la température d'oxydation. Pour les températures où les alliages restent vitreux, la formation de l'oxyde est due à la fois à la diffusion des ions d'oxygène vers l'intérieur mais également à la diffusion des ions de zirconium vers l'extérieur de la couche. En revanche, pour les alliages en cours de cristallisation, la diffusion des ions de zirconium est progressivement freinée jusqu'à devenue nulle lorsque la cristallisation est totale, l'oxydation étant alors uniquement contrôlée par la diffusion interne des ions d'oxygène.Les contraintes résiduelles correspondantes au sein de couches d'oxyde sont de compression, aussi bien pour les alliages amorphes que pour les alliages cristallisés, et qu'elles augmentent linéairement avec l'épaisseur de la couche. Les contraintes de croissance lors de l'oxydation et résiduelle après l'oxydation sont fortement influencées par le changement de phase de la zircone.

Les travaux de thèse traitent du problème de la planification de manoeuvres pour le maintien à poste de satellites géostationnaires équipés de tuyères électriques (à poussée faible). Nous évaluons l'opportunité de substituer une telle planification à celle traditionnellement utilisée pour les satellites géostationnaires équipés de tuyères chimiques (à poussée forte). Dès son apparition, la technologie des systèmes de propulsion à poussée faible a rencontré un vif intérêt auprès des agences et des sociétés spatiales. Grâce à sa haute impulsion spécifique (qui implique une basse consommation de carburant), cette technologie est devenue très compétitive par rapport à la technologie traditionnelle des propulseurs chimiques à poussée forte, surtout dans les phases de transfert et rendez-vous des missions spatiales. Pendant la définition des missions à poussée faible, les analyses de faisabilité des phases de transfert et rendez-vous (via la solution de problèmes d'optimisation de trajectoire) ont été réalisées avec des solutions d'optimisation alternatives. En effet, pendant ces phases, il est nécessaire d'activer les systèmes de propulsion à faible poussée sur des longues portions du temps de transfert. Par conséquent, les problèmes d'optimisation de trajectoire à poussée forte (typiquement formulés en temps discret) ont été remplacés par des problèmes d'optimisation de trajectoire à poussée faible formulés en temps continu et résolus par des techniques de contrôle en temps continu. Le premier objectif de cette thèse est de comprendre quel est l'impact de la technologie à faible poussée lors de l'analyse de faisabilité de la phase de maintien à poste de satellites géostationnaires. Nous étudions en particulier l'impact de l'utilisation des systèmes de propulsion à faible poussée sur la planification de manoeuvres et sur la boucle entière de maintien à poste géostationnaire. L'étude consiste à déduire si la planification de manoeuvres à poussée faible est compétitive au regard des stratégies classiques de planification couramment employées pour des manoeuvres à poussée forte. Généralement, les stratégies classiques à long terme pour le maintien à poste sont déduites de modèles de propagation d'orbite simplifiés (en fonctions des paramètres orbitaux moyennés) par la conjonction des trois facteurs suivants : la forte poussée des propulseurs, la dimension de la fenêtre de maintien à poste pas très contraignante ainsi que la possibilité d'exécuter des manoeuvres à basse fréquence. Dans le cadre de cette thèse, compte tenu du faible niveau des poussées et des contraintes strictes en position (fenêtres de maintien à poste petites), nous considérons comme plus appropriés l'hypothèse d'une plus haute fréquence de manoeuvres et l'utilisation d'un modèle de propagation d'orbite en fonction de paramètres osculateurs. Pour la planification de manoeuvres, nous proposons une solution par approche directe : le problème de maintien à poste en tant que problème de contrôle optimal est discrétisé et traduit en un problème d'optimisation paramétrique. Deux techniques différentes d'optimisation sont proposées : l'optimisation sous contraintes à horizon fixe et celle à horizon glissant. Cette deuxième technique est appliquée aux équations linéarisées du mouvement préalablement transformées via un changement de variable à la Lyapunov sur l'état des déviations des paramètres équinoxiaux osculateurs. Cette transformation de Lyapunov définit des nouveaux paramètres orbitaux. Elle rend le processus de planification plus compréhensible du point de vue du contrôle et plus facile à implémenter d'un point de vue numérique, grâce aux concepts de platitude et inclusion différentielles. Les résultats de la planification de manoeuvres à poussée faible sont obtenus dans un premier temps en fonction des changements de vitesse, dans un deuxième temps en fonction des forces engendrées par les tuyères des systèmes de propulsion classiques. Le but est de déterminer la solution la plus efficace en conditions nominales et en cas de panne d'un des propulseurs. Le problème du positionnement simultané de plusieurs satellites dans une même grande fenêtre de maintien à poste n'est pas adressé explicitement. Il est implicitement résolu en proposant une technique fine de contrôle pour maintenir chaque satellite à poste dans une fenêtre de dimension très petite.

Cette thèse porte sur les relations entre la limite supérieure de la zone sismogène et la construction et la structuration du prisme de Nankai, au Sud-Est du Japon. La détermination des contraintes dans le prisme à partir des données d’ovalisation de puits montrent une permutation des contraintes au dessus d’un hors séquence, appelé splay fault, situé à la limite de la zone sismogène. Le régime de contrainte est extensif dans cette zone et décrochant dans le reste du prisme. L’étude de la compaction des sédiments à partir des données de porosité corrigée de l’effet des argiles permet de quantifier l’érosion récente des sédiments de pentes, et montre l’existence d’une érosion plus ancienne au dessus des chevauchements majeurs. L’étude des styles de déformations des sédiments montre que l’érosion peut rendre compte de la répartition des structures de déformation dans le prisme. Enfin, la modélisation mécanique de la marge montre que l’activité de la splay fault, très faible depuis le dernier million d’années, et la présence d’une zone d’extension en arrière de cette faille s’expliquent par l’arrêt de l’accrétion et un décollement faible dans le prisme frontal. / This study is focused on the upper limit of the seismogenic zone in the Kumano transect of the Nankai margin in relation with the construction of the wedge. The first step was to determine the stress orientations and amplitude within the wedge from borehole breakouts analysis. We show the striking presence of a trench normal extension zone behind the splay fault in a wedge dominated by a strike-slip stress regime. Locally an extensional stress regime is also observed in the slope sediment of the outer part of the wedge. In a second part of this work, we used porosity corrected from water bound to clay minerals to study the compaction state of the sediments and to quantify the amount of erosion in the slope sediments. Results also highlight the existence of past erosion related to the activity of the splay fault. In a third part, the observation of deformation style within sediments showed that erosion could explain the distribution of dilatant and compactive structures within the wedge. Finally, we show, based on mechanical modeling, that the splay fault slow slip rate over the last million year, and the presence of an extension zone landward of the splay fault, could both be explained by a pause in the accretion and a weak décollement beneath the outer wedge.

Le stockage des céréales se fait en milieu rural au nord du Bénin dans des greniers (silo) en terre qui ont une contenance limitée. Les nombreuses tentatives d’introduction de magasins modernes ont essuyé le rejet des populations. Le présent travail a pour but d’étudier la possibilité d’augmenter le volume de stockage du grenier de type Yom sans changer ni sa forme architecturale, ni le composite argile-paille-décoction de néré. La caractérisation expérimentale du comportement mécanique du banco, matériau de construction de ces greniers a permis de définir le meilleur dosage pour des caractéristiques maximales. En faisant varier la proportion de paille et de décoction de néré, l’influence sur la terre crue de la paille et de la décoction est mise en évidence. La simulation de l’action des grains, à la vidange et au remplissage sur la paroi du grenier et l’analyse des déformations induites montrent qu’il y a possibilité de construire des greniers de plus grande taille. Une optimisation paramétrique des dimensions a été faite. L'effet des échauffements successifs sur la paroi du grenier a été pris en compte à travers la modélisation des transferts thermiques. La modélisation des profils de température dans diverses coupes orientées du silo montre que le flux est bien amorti et les conditions d'ambiance dans le grenier garantissent une bonne conservation. Il est donc possible de construire des greniers de grande capacité de stockage qui serviraient de greniers communautaires ou banques de céréales à disposition des coopératives villageoises. Cette solution règlerait le problème de l’insécurité alimentaire au niveau du monde rural / The storage of cereals is done in rural areas in northern Benin in earth silo that have a limited capacity. The numerous attempts to introduce modern shops have been rejected by the population. This work aims to study the possibility of increasing the storage volume of Yom type silo without changing its architectural form, nor the clay-straw-nere composite. The experimental characterization of the mechanical behaviour of the banco, the construction material of these silo made it possible to define the best dosage for maximum characteristics. By varying the proportion of straw and decoction of nere, the influence on the raw earth of the straw and the decoction is highlighted. The simulation of grain action (filling and discharging) on the silo wall and the analysis of the induced deformations shows that there is a possibility of building larger earth silo. A parametric optimization of the dimensions was made. The effect of successive heating on the silo wall has been considered by modelling heat transfer. The modelling of the temperature profiles in various silo-oriented sections shows that the flow is well damped and the ambient conditions in the attic guarantee good conservation. It is therefore possible to build earth silo with large storage capacity that could be used as community silo or cereal banks for village cooperatives. This solution would solve the problem of food insecurity in rural areas

La programmation par contraintes intègre des algorithmes de résolution génériques dans des langages de modélisation déclaratifs basés sur les contraintes : ces langages permettent de décrire des problèmes combinatoires sous la forme d’un ensemble de variables devant prendre leurs valeurs dans des domaines en satisfaisant des contraintes. Nous introduisons dans cette thèse un algorithme de résolution générique paramétré par : — une stratégie d’exploration de l’espace de recherche, à choisir parmi, chronological backtracking, conflict directed backjumping, conflict directed backjumping with reordering, dynamic backtracking, decision repair, et backtracking with tree decomposition ; — une heuristique de choix de variables, à choisir parmi, min-domain/ddeg et min-domain/wdeg ; — une technique de propagation de contraintes, à choisir parmi, forward checking et maintaining arc consistency. Ainsi, cet algorithme générique s’instancie en vingt-quatre configurations différentes ; certaines correspondant à des algorithmes connus, d’autres étant nouvelles. Ces vingt- quatre configurations ont été comparées expérimentalement sur un benchmark de plus de mille instances, chaque configuration étant exécutée plusieurs fois sur chaque instance pour tenir compte du non déterminisme des exécutions. Des tests statistiques sont utilisés pour comparer les performances. Cette évaluation expérimentale a permis de mieux comprendre la complémentarité des différents mécanismes de résolution, avec une attention particulière portée sur la capacité à tirer parti de la structure des instances pour accélérer la résolution. Nous identifions notamment treize configurations complémentaires telles que chaque instance de notre benchmark est bien résolue par au moins une des treize configurations. Une deuxième contribution de la thèse est d’introduire un sélecteur capable de choisir automatiquement la meilleure configuration de notre algorithme générique pour chaque nouvelle instance à résoudre : nous décrivons chaque instance par un ensemble de descripteurs et nous utilisons des techniques d’apprentissage automatique pour construire un modèle de choix de configuration à partir de ces descripteurs. Sachant que l’apprentissage est généralement plus difficile quand il y a beaucoup de configurations, nous exprimons le problème du choix du sous-ensemble de configurations pouvant être sélectionnées comme un problème de couverture d’ensemble et nous comparons deux critères de choix : le premier vise à maximiser le nombre d’instances résolues par au moins une configuration et le second vise à maximiser le nombre d’instances pour lesquelles il y a une bonne configuration disponible. Nous montrons expérimentalement que la seconde stratégie obtient généralement de meilleurs résultats, et que le sélecteur obtient de meilleures performances que chacune de nos vingt-quatre configurations initiales. / Constraint programming integrates generic solving algorithms within declarative languages based on constraints : these languages allow us to describe combinatorial problems as a set of variables which have to take their values in domains while satisfying constraints. Numerous real-life problems can be modelled in such a way, as for instance, planification problems, scheduling problems, . . . These problems are NP-complete in the general case of finite domains. We introduce in this work a generic solving algorithm parameterized by : — a strategy for exploring the search space, to be chosen from the following six, chronological backtracking, conflict directed backjumping, conflict directed backjumping with reordering, dynamic backtracking, decision repair, and backtracking with tree decomposition ; — a variable ordering heuristic, to be chosen from the following two, min-domain/ddeg and min-domain/wdeg ; — a constraint propagation technique, to be chosen from the following two, forward checking and maintaining arc consistency. Thus, this algorithm leads to 24 different configurations ; some corresponding to already known algorithms, others being new. These 24 configurations have been com- pared experimentally on a benchmark of more than a thousand instances, each configuration being executed several times to take into account the non-determinism of the executions, and a statistical test has been used to compare performances. This experimental evaluation allowed us to better understand the complementarity of the different solving mechanisms, with a focus on the ability to exploit the structure of the instances to speed up the solving process. We identify 13 complementary configurations such that every instance of our benchmark is well solved by at least one of the 13 configurations. A second contribution of this work is to introduce a selector able to choose automatically the best configuration of our generic solver for each new instance to be solved : we describe each instance by a set of features and we use machine learning techniques to build a model to choose a configuration based on these features. Knowing that the learning process is generally harder when there are many configurations to choose from, we state the problem of choosing a subset of configurations that can be picked as a set covering problem and we compare two criterion : the first one aims to maximize the number of instances solved by at least one configuration and the second one aims to maximize the number of instances for which there is a good configuration available. We show experimentally that the second strategy obtains generally better results and that the selector obtains better performances than each of the 24 initial configurations.

Les verres métalliques massifs à base de zirconium, développés depuis la fin des années 80, présentent des propriétés mécaniques très intéressantes. Ils peuvent être envisagés pour de multiples applications y compris à des températures élevées et il est donc intéressant d’étudier leur résistance à l’oxydation dans de telles conditions. L’objectif fondamental de cette étude consiste à mieux comprendre le rôle de divers paramètre thermodynamiques et chimiques sur le comportement à l’oxydation des verres métalliques à base de zirconium à des températures intermédiaires sous air sec, à déterminer les contraintes résiduelles au sein de la couche d’oxyde formée, en comparaison avec des alliages amorphe cristallisés après un traitement de recuit. La cinétique d’oxydation de ces verres et la structure cristalline de la couche d’oxyde ZrO2 dépend de la température et de la durée d’oxydation : pour des durées courtes d’oxydation et pour une température légèrement inférieure à Tg, la cinétique d’évolution est parabolique alors que si l’échantillon est oxydé à T > Tg, la loi cinétique peut être décomposée en deux parties. Les mêmes alliages cristallisés après un traitement de recuit, s’oxydent selon une loi parabolique quelque soit la température. Pour des durées d’oxydation longues à une température proche de Tg, les lois cinétiques deviennent plus complexes et la cristallisation du verre métallique pouvant avoir lieu au cours des essais d’oxydation. De même la structure cristalline des couches d’oxyde dépend de la température d’oxydation. Ainsi, pour T < Tg, la couche d’oxyde des alliages amorphes est uniquement constituée de la zircone tétragonale alors que pour ces mêmes alliages portés à une température T > Tg et pour les alliages cristallisés, on observe de la zircone monoclinique. Les mécanismes d’oxydation dépendent de la température d’oxydation. Pour les températures où les alliages restent vitreux, la formation de l’oxyde est due à la fois à la diffusion des ions d’oxygène vers l’intérieur mais également à la diffusion des ions de zirconium vers l’extérieur de la couche. En revanche, pour les alliages en cours de cristallisation, la diffusion des ions de zirconium est progressivement freinée jusqu’à devenue nulle lorsque la cristallisation est totale, l’oxydation étant alors uniquement contrôlée par la diffusion interne des ions d’oxygène.Les contraintes résiduelles correspondantes au sein de couches d’oxyde sont de compression, aussi bien pour les alliages amorphes que pour les alliages cristallisés, et qu’elles augmentent linéairement avec l’épaisseur de la couche. Les contraintes de croissance lors de l’oxydation et résiduelle après l’oxydation sont fortement influencées par le changement de phase de la zircone. / The Zr-based bulk metallic glasses, developed since the late 1980s, have very interesting mechanical properties, which can be considered for many applications including working under oxidizing atmosphere conditions at high temperatures. It is therefore interesting to study their oxidation resistance and to characterize the oxide scale formed on alloys’ surface. The fundamental objective of this thesis is to enhance the understanding of the role of various thermodynamic and chemistry parameters on the oxidation behaviour of the Zr-based bulk metallic glasses at high temperature under dry air, to determine the residual stresses in the oxide layer, in comparison with their crystalline alloys with the same chemical composition after an annealing treatment. The oxidation kinetics of these glasses and the crystalline structure of oxide scale ZrO2 depend on the temperature and the oxidation duration: for short periods of oxidation or at a temperature below Tg, the kinetics follows a parabolic law, whereas, if the sample is oxidized at T > Tg, the kinetics can be divided into two parts. The crystalline counterparts are oxidized by a parabolic rule whatever the temperature; for long oxidation duration at a temperature close to Tg, the kinetics becomes more complex because of the crystallisation of the glasses during the oxidation tests. Also the crystalline structure of the oxide layers depends on the oxidation temperature: the oxide layer consists only in tetragonal Zirconia at T < Tg, while monoclinic Zirconia was formed at higher temperature. The mechanism of the formation of the oxide scale is due to both the interior diffusion of Oxygen ions and the external diffusion of Zirconium ions. However the diffusion of Zirconium ions slows gradually during the crystallisation process of the glass matrix. When the crystallisation is completed, the formation of Zirconia is controlled by only the internal diffusion of oxygen ions. The corresponding residual stresses existing in the oxide layer are compressive, both for the metallic glass and counterpart, and they increase linearly with the thickness of the oxide layer. The growth and residual stresses are strongly influenced by the phase transition of Zirconia during oxidation process.