Numerical simulation of reinforced concrete nuclear containment under extreme loads

A finite element model for the non-linear dynamic analysis of a reinforced concrete (RC) containment shell of a nuclear power plant subjected to extreme loads such as impact and earthquake is presented in this work. The impact is modeled by using an uncoupled approach in which a load function is applied at the impact zone. The earthquake load is modeled by prescribing ground accelerations at the base of the structure. The nuclear containment is discretized spatially by using 20-node brick finite elements. The concrete in compression is modeled by using a modified $Dr{\ddot{u}}cker$-Prager elasto-plastic constitutive law where strain rate effects are considered. Cracking of concrete is modeled by using a smeared cracking approach where the tension-stiffening effect is included via a strain-softening rule. A model based on fracture mechanics, using the concept of constant fracture energy release, is used to relate the strain softening effect to the element size in order to guaranty mesh independency in the numerical prediction. The reinforcing bars are represented by incorporated membrane elements with a von Mises elasto-plastic law. Two benchmarks are used to verify the numerical implementation of the present model. Results are presented graphically in terms of displacement histories and cracking patterns. Finally, the influence of the shear transfer model used for cracked concrete as well as the effect due to a base slab incorporation in the numerical modeling are analyzed.