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.