Stress corrosion cracking (SCC) is an important destruction form of materials such as stainless steel, nickel-based alloy and their welded components in nuclear reactor pressure vessels and pipes. The existing popular quantitative prediction models of SCC crack growth rate are mainly influenced by fracture toughness values K_Jc or J_c. In particular, the composite constraint, containing the in-plane constraints and out-of-plane constraints around the crack front, has a significant influence on the fracture toughness of structures in nuclear power plants. Since the plastic strain gradient is a characterization parameter of the quantitative prediction model for crack growth rate, it may be a characterization parameter of composite constraint. On the basis of the experimental data at a low temperature of alloy steel 22NiMoCr3-7 used in nuclear pressure vessels, the gradient of equivalent plastic strain D_PEEQ around the crack fronts at different constraint levels was calculated using the finite element method, which introduces a new non-dimensional constraint parameter D_p, to uniformly characterize the in-plane and out-of-plane constraint effects. Compared with constraint parameters A_PEEQ or A_p, the process of obtaining parameters D_PEEQ or D_p is much simpler and easier. In a wide range, a single correlation curve was drawn between parameter D_p and normalized fracture toughness values K_Jc/K_ref or J_c/J_ref of specimens at a low or high constraint level. Therefore, regardless of whether the constraint levels of the structures or standard specimens are low or high, constraint parameter D_p can be used to measure their fracture toughness. To build an evaluation method that has structural integrity and safety while containing the composite constraint effects, in addition to accurate theoretical interpretation, further verification experiments, numerical simulations and detailed discussions are still needed.

Hydrogen embrittlement, as one of the major concerns for austenitic stainless steel, is closely linked to the diffusion of hydrogen through the grain boundary of materials. The phenomenon is still not well understood yet, especially the full interaction between hydrogen diffusion and the misorientation of the grains. This work aimed at the development of a robust numerical strategy to model the full coupling of the hydrogen diffusion and the anisotropic behavior of crystals in 316 stainless steel. A constitutive model, which allows easy incorporation of crystal orientation, various loading conditions, and arbitrary model geometries, was established by using the finite element package ABAQUS. The study focuses on three different bicrystal models composed of misoriented crystals, and the results indicate that the redistribution of hydrogen is significant closely to the grain boundary, and the redistribution is driven by the hydrostatic pressure caused by the misorientation of two neighboring grains. A higher elastic modulus ratio along the tensile direction will lead to a higher hydrogen concentration difference in the two grains equidistant from the grain boundary. The hydrogen concentration shows a high value in the crystal along the direction with stiff elastic modulus. Moreover, there exists a large hydrogen concentration gradient in a narrow region very close to the grain boundary to balance the concentration difference of the neighboring grains.

The material constants calculation models for hyperbolic-sine creep model were proposed. The material constants used in hyperbolic-sine creep model for 316 stainless steel were calculated due to the models proposed and experimental data in the temperature range from 873K to 1023K. The relationships between material constants of 316 stainless steel creep model and temperature were obtained by curve fitting. The creep rate predict model of 316 stainless steel with only stress and temperature was also developed, the creep rates predicted were in good agreement with experimental data.

The stress state at the crack tip of structural components in nuclear power plants used in the SCC quantitative prediction models is based on the assumption that the polycrystalline material is isotropic and homogeneous at present. However, the crystals in polycrystalline materials are anisotropic with different orientations, which would induce a nonuniform stress to cause the initiation and propagation of SCC. By using a finite element method, the elastic responses of anisotropic behaviors of austenitic stainless steel bicrystals are studied. The results indicate that the stress distribution near GBs depends strongly on the crystal orientation. A larger Mises stress concentration exists on the GB with larger stiffness along the load direction. The Mises stress difference is higher in the bicrystal with bigger elastic modulus difference of two neighboring grains along the tensile axis. In the bicrystal with GB perpendicular to the tensile axis, the grain orientation has little effects on the Mises stress far from the GB in both grains. The strain inconsistency in bicrystals is affected by the mismatch of two neighboring grains. The larger the elastic modulus differences between two neighboring grains caused by misorientation, the larger the strain inconsistency in the bicrystal.