Continuum Damage Dynamic Model Combined with Transient Elastic Equation and Heat Conduction Equation to Solve RPV Stress

Abstract

The development of the world cannot be separated from energy: the energy crisis has become a major challenge in this era, and nuclear energy has been applied to many fields. This paper mainly studies the stress change of reaction pressure vessels (RPV). We established several different physical models to solve the same mechanical problem. Numerical methods range from 1D to 3D; the 1D model is mainly based on the mechanical equilibrium equations established by the internal pressure of RPV, the hoop stress, and the axial stress. We found that the hoop stress is twice the axial stress; this model is a rough estimate. For 2D RPV mechanical simulation, we proposed a new method, which combined the continuum damage dynamic model with the transient cross-section finite element method (CDDM-TCFEM). The advantage is that the temperature and shear strain can be linked by the damage factor effect on the elastic model and Poission ratio. The results show that with the increase of temperature (damage factor μ^,d^), the Young’s modulus decreases point by point, and the Poisson’s ratio increases with the increase of temperature (damage factor μ^,Et). The advantage of the CDDM-TCFEM is that the calculation efficiency is high. However, it is unable to obtain the overall mechanical cloud map. In order to solve this problem, we established the axisymmetric finite element model, and the results show that the stress value at both ends of RPV is significantly greater than that in the middle of the container. Meanwhile, the shape changes of 2D and 3D RPV are calculated and visualized. Finally, a 3D thermal–mechanical coupling model is established, and the cloud map of strain and displacement are also visualized. We found that the stress of the vessel wall near the nozzle decreases gradually from the inside surface to the outside, and the hoop stress is slightly larger than the axial stress. The main contribution of this paper is to establish a CDDM-TCFEM model considering the influence of temperature on elastic modulus and Poission ratio. It can dynamically describe the stress change of RPV; we have given the fitting formula of the internal temperature and pressure of RPV changing with time. We also establish a 3D coupling model and use the adaptive mesh to discretize the pipe. The numerical discrete theory of FDM-FEM is given, and the numerical results are visualized well. In addition, we have given error estimation for h-type and p-type adaptive meshes. So, our research can provide mechanical theoretical support for nuclear energy safety applications and RPV design.