Modelling thermo-mechanical response of metal canister disposed in engineered barrier systems

Abstract

The main objective of this thesis is to describe and implement a rate-dependent (viscoplastic) constitutive equations for copper based canister using a single internal state variable and use it for modelling disposal schemes subjected to shear conditions. The model of Bodner and Partom (BP) is considered for this purpose (Bodner & Partom, 1975). To evaluate the copper's response to variations in strain rate throughout a temperature range, this model uses viscoplastic constitutive equations based on a single internal state variable that is a function of plastic work. The BP constitutive model has been implemented in CODE_BRIGHT computer code (DIT-UPC, 2021; Olivella et al., 1994; Olivella et al., 1996). The constitutive equation is verified against existing results and its validation is carried out by assessing the capability to reproduce experiments. For verification, the stress-strain response of a bar under uniaxial tensile conditions at a constant velocity was simulated with CODE_BRIGHT and the results compared with the solution obtained from Stealth Finite Difference Code and the semi-analytical solutions. In addition, it is analyzed the reaction of a canister installed inside an engineered barrier which undergoes hydration and swelling until it reaches full saturation prior to shearing. The implemented elasto-viscoplastic model actually represents the behavior of the canister in yielding conditions. Although the canister is very stiff and strong compared to the clay components where it will be emplaced in a clay volume shear deformations occurring on the engineered barrier system (EBS) will have an impact on the canister (Börgesson, 1986). Finally, a repository of nuclear spent fuel is analyzed under extreme conditions as it is a long-life system and has to resist all types of environmental conditions, for instance earthquakes and glaciations. After model implementation and verification in CODE BRIGHT this program has been used to model the canister-clay shear tests in 2D and 3D. In the Mock-up scale and Full-scale in 2D and 3D models, a comparison between the total stress method (mechanical analysis only) and the effective stress method (hydro-mechanical analysis) has been investigated. Experimental results from a canister-clay shear test (Börgesson, 1986) have been compared to the findings of the numerical model using CODE BRIGHT. Sensitivity analysis on viscosity for total stress calculations and permeability for effective stress calculations have been performed in the Mock-up scale and Full-scale. In the canister-clay shear test, a mesh sensitivity analysis and also a comparison between the updated mesh (Lagrangian method) and the fixed mesh method have been carried out. Finally, in two and three dimensions, the simulation of canister-clay shear tests based on impact of discontinuity (the discontinuity is in the rock and due to it, there is shearing in canister-rock system), as well as sensitivity analysis of parameters and mesh in mechanical and a coupled hydro-mechanical analysis, have been investigated.

El principal objetivo de esta tesis es describir e implementar una ecuación constitutiva dependiente de la velocidad (viscoplástica) para un contenedor a base de cobre usando una única variable de estado interna y usarlo para modelar esquemas de almacenamiento sujetos a condiciones de cizallamiento. Para ello se considera el modelo de Bodner y Partom (BP) (Bodner & Partom, 1975). Para evaluar la respuesta del cobre a las variaciones en la velocidad de deformación a lo largo de un rango de temperatura, este modelo utiliza ecuaciones constitutivas viscoplásticas basadas en una única variable de estado interno que es una función del trabajo plástico. El modelo constitutivo de BP se ha implementado en programa de ordenador CODE_BRIGHT (DIT-UPC, 2021; Olivella et al., 1994; Olivella et al., 1996). La ecuación constitutiva se verifica con los resultados existentes y su validación se lleva a cabo evaluando la capacidad para reproducir experimentos. Para la verificación, se simuló la respuesta tensión-deformación de una barra en condiciones de tracción uniaxial a velocidad constante con CODE_BRIGHT y se compararon los resultados con la solución obtenida de Stealth Finite Difference Code y las soluciones semi-analíticas. Además, se analiza la reacción de un contenedor instalado dentro de una barrera de ingeniería, que sufre hidratación e hinchazón hasta que alcanza la saturación completa antes del cizallamiento. El modelo elasto-viscoplástico implementado realmente representa el comportamiento del contenedor en condiciones de plastidad. Aunque el contenedor es muy rígido y resistente en comparación con los componentes de arcilla donde se colocará, las deformaciones por cizallamiento que ocurren en el sistema de barrera de ingeniería (EBS) tendrán un impacto en el contenedor (Börgesson, 1986). Por último, el contenedor de combustible nuclear gastado se podrá analizar en condiciones extremas, ya que es un sistema de larga duración y debe resistir todo tipo de condiciones ambientales, por ejemplo, terremotos y glaciaciones. Después de la implementación y verificación del modelo en CODE BRIGHT, este programa se ha utilizado para modelar las pruebas de cizallamiento en 2D y 3D. En la escala Mock-up y Full-scale en modelos 2D y 3D, se ha llevado a cabo una comparación entre el método de tensiones totales (solo análisis mecánico) y el método de tensiones efectivas (análisis hidromecánico). Los resultados experimentales de una prueba de cizallamiento canister en arcilla (Börgesson, 1986) se han comparado con los resultados del modelo numérico utilizando CODE BRIGHT. El análisis de sensibilidad sobre la viscosidad para los cálculos de tensión total y la permeabilidad para los cálculos de tensión efectiva se han realizado en la escala Mock-up y Full-scale. En la prueba de cizalla de contenedor-arcilla, se ha realizado un análisis de sensibilidad de la malla y también una comparación entre el método que actualiza la malla con los desplazamientos (método Lagrangiano) y el método de malla fija. Finalmente, en dos y tres dimensiones, se ha realizado la simulación de ensayos de cizallamiento contenedor-arcilla basados en el impacto de la presencia de una discontinuidad (la discontinuidad está en la roca y por ello hay cizallamiento en el sistema), así como análisis de sensibilidad de parámetros. y malla en análisis mecánico e hidromecánico acoplado