Effects of Grain Orientation on Stress State near Grain Boundary of Austenitic Stainless Steel Bicrystals

/e 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. /e 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./eMises 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. /e strain inconsistency in bicrystals is affected by the mismatch of two neighboring grains. /e larger the elastic modulus differences between two neighboring grains caused by misorientation, the larger the strain inconsistency in the bicrystal.


Introduction
Stress corrosion cracking (SCC) of structural components is one of the major factors that affect the reliability and integrity of the nuclear power plants (NPPs).Indeed, intergranular stress corrosion cracking (IGSCC) has occurred in many internal components made up of stainless steel, nickel alloy, and their weld metals [1][2][3][4].e slip/dissolution-oxidation model has been widely accepted as a reasonable description of SCC in an oxygenated aqueous system [5][6][7].In this model, SCC is regarded as the result of localized oxidation enhanced by stress/strain, and the crack tip strain rate is expressed as a unique parameter to cause the physical degradation of the interfacial film.Many attempts have been made to acquire the stress state of the crack tip due to its importance [6,7]; however, the stress and strain are obtained macroscopically based on the assumption that the materials are isotropic and homogeneous in terms of elastic deformation when the materials have random crystallographic and morphologic texture.In fact, the crack tip stress state is significantly affected by the microstructure of materials due to the misorientation and anisotropic of each component crystal, and the crack tip would have a nonuniform stress at the microstructural level even under a uniform remote stress condition [8][9][10][11].e inhomogeneity stress distribution caused by anisotropic of grains would play an important role in the initiation and propagation of SCC [12].Researches have shown that the initiation of IGSCC in nickel alloys in high-temperature water greatly depends on inclination of the grain boundary (GB) to the tensile axis [13]; the cracks tend to be initiated easily at the GB whose inclination is perpendicular to the tensile axis and suggested that the local stress perpendicular to the GB governs the cracking behavior; this result seems to be consistent with the theory.However, other experiments and observations show that cracks could travel in different directions, such as at an angle perpendicular or even parallel to the applied load [14][15][16].
In order to precisely predict cracking characteristics, researchers have sought to acquire the local stress state near an individual GB with random shape and orientation of the crystal grains, and to find the driving force for cracking [17].e test of bicrystals with a single symmetrical 111-tilt GB using the slow strain rate technique (SSRT) shows that the susceptibility to IGSCC depends on the misorientation rather than on the GB energy.
e local stress concentration at large-angle GBs could be attributed to a high susceptibility to the IGSCC [18].Single GBs in real alloys were also generated and monitored by using cantilevers manufactured by focused ion beam (FIB) and tested by nanoindentation, which was expected to provide a more accurate method of measuring the crack growth and investigate the dependence of SCC resistance of the single GB on GB character and composition [19][20][21][22].However, it is difficult to manufacture a specific single GB with respect to the random orientation of neighboring grains and also to accurately measure the stress and strain distribution near the GB.us, a three-dimensional anisotropic finite element method (FEM) is performed in this paper to determine how the crystal orientation affects the stress state near the GB of the bicrystal.

Finite Element Model
2.1.Geometry Model.Two types of GBs are considered in bicrystals with respect to the relationship between the GB and tensile axis.e type-I GB is perpendicular to the tensile axis as shown in Figure 1, and the type-II GB is parallel to the tensile axis as shown in Figure 2. e average grain size of 316L stainless steel varies from 17 μm to 200 μm with different aging time, temperatures, and other factors in literature [23,24]; the dimension of each component grain in the bicrystal with type-I GB is assumed to be 20 μm × 20 μm × 50 μm, and each grain in the bicrystal with type-II GB is of the dimension 10 μm × 20 μm × 50 μm in this study.

Material Propriety.
Stainless steel 316L is a polycrystalline material composed of austenitic crystal which has the face-centred cubic (FCC) configuration.e single crystal of 316L is orthotropic, and the elastic anisotropy is modeled using the elastic constants d 11 � 204.6 GPa, d 12 � 137.7 GPa, and d 44 � 126.2 GPa [25].When the crystallographic axial system coincides with the global coordinate, the elastic stress and strain relationship of the single crystal is described by the generalized Hooke's law: where where σ and ε denote the stress tensor and strain tensor, respectively, and D is the stiffness matrix.In a polycrystalline material, the crystal orientation is random and always differs from the global coordinate.e stress and strain should be recalculated with the rotation stiffness matrix D xyz as where A is the rotation matrix, which is given as  2 Advances in Materials Science and Engineering where l i , m i , and n i (i � 1, 2, and 3) are the direction cosines of X, Y, and Z axes between the global coordinate and crystal coordinate.e elastic modulus E, Poisson's ratio μ, and shear modulus G are equal, respectively, in the three main axes, and they can be calculated by elastic constants: By substituting d 11 , d 12 , and d 44 into ( 7)-( 9), it can be obtained that E � 94.06 GPa, G � 126.2 GPa, and μ � 0.402.
e elastic modulus of the single crystal in <110> and <111> crystal orientations can be described with E, G, and μ in the three main axes as According to (10) and (11), the elastic modulus along <110> and <111> orientations is E 110 � 193.81 GPa and E 111 � 299.77GPa, respectively.e highest elastic stiffness is along the <111> orientation, and the lowest stiffness is along the <100> orientation, as shown in Figure 3.For different crystallographic directions, Considering the relationship between [100], [110], and [111] orientations and the tensile axis, three types of grains    In order to deeply understand the GB characterization, three types of incompatibility bicrystal models were built by combining grains with different orientations mentioned above, as listed in Table 2. e orientation of each grain is defined in the finite element model by rotating the material coordinate.

Boundary Condition and Loading Process.
e 8-node linear brick C3D8R element was used in the finite element model, and a fine mesh was carried out at GBs to acquire more accurate stress and strain.e surface perpendicular to direction 1 of grain A is fixed initially with U 1 � U 2 � U 3 � 0, followed by a constant displacement U 1 equal to 0.5 along direction 1 applied on the surface of grain B, and this causes a total stain of 0.5% along direction 1 for the whole bicrystal.Advances in Materials Science and Engineering models are shown in Figure 4, which illustrates that the heterogeneity of two adjacent grains will lead to a stress concentration at the GB and the largest stress exists in the grain with higher elastic modulus along the tensile axis.e Mises stress distributions near the GB are different in the three bicrystals, which depend on the crystallographic orientation of the two adjacent grains and GBs.

Results and Discussion
e maximum Mises stress appears on the GB of model III bicrystal.
More intuitive stress distribution was shown in Figure 5, which illustrates the stress distribution on different cross sections of model I bicrystal.On the cross section ±1 μm apart from the GB, the Mises stress shows a plateau in the middle of the cross section and a high gradient close to the surfaces, and the Mises stress gradient is higher along axis 3 than axis 2. With the cross section away from the GB, the high-stress area in the cross section reduces, accompanied by the reduction of the Mises stress gradient.e stress inconsistency at the GB induced by misorientation disappears when the cross section is far from the GB.
Let ε xA and ε xB represent the average strain ε xx of component crystals A and B along axis 1 and let ε A and ε B represent the average strain of all elements in contact with the GB in component crystals A and B, respectively.e characteristic of deformation near the GB of three incompatibility bicrystal models was evaluated, as listed in Table 3. ε xA /ε xB is bigger than 1, which shows that the average strain of grain A is higher than that of grain B, and the average strain on the GB of grain A is also higher than that of grain B as represented by ε A /ε B . is indicates that the strain in the softer component (grain A) is higher than that in the harder component (grain B). e values of ε A /ε xA denote that the strain near the GB is lower than the average strain of this grain in the softer component, while the strain near the GB is higher than the average strain of this grain in the harder component as denoted by ε B /ε xB .
e strain inconsistency between two neighboring grains and between the GB and the entire grain is affected by the mismatch of two neighboring grains.In model II bicrystal, the elastic modulus difference of two grains is the largest; thus, the strain inconsistency in this bicrystal is the most serious.
To study the stress state caused by misorientation between grain A and grain B, a stress concentration factor λ is introduced to normalize the stress as where the actual stress and average absolute stress are observed along two paths, as shown in Figure 6.Path 1 penetrates the GB from grain A to grain B and path 2 locates on the GB. e average absolute stress is calculated along path 1 and path 2, respectively.Figure 7 illustrates the stress concentration factor for path 1, and the Mises stress concentration factor is approximately flat in the middle and increases or decreases close to the edges in the softer and harder components, respectively.e Mises stress is higher in the harder component than in the softer component in each model.By defining the stress difference coefficient κ as where σ A and σ B are the stress in grain A and grain B, respectively, and the Mises stress difference for each model was demonstrated in Figure 8. e maximum stress   Advances in Materials Science and Engineering difference appears at the middle of the GB and reduces gradually away from the middle in model II, which has the largest elastic modulus difference along the tensile axis, while the stress difference coefficient behaves constant in the middle and reduces slightly close to the edges in the other models.
e minimum stress difference near the GB of model III has the smallest elastic modulus difference between the two adjacent grains along the tensile axis.
As shown in Figures 9 and 10, the Mises stress concentration factor and Mises stress difference coefficient for incompatibility model I and model II have the same tendency along path 2, which clearly demonstrates that a stress concentration exists near the GB and a uniform stress distribution in the regions is far from the GB in both grains.In the two harder components, the stress concentration factors are approximately 1.6 and 1.4, respectively, and both reduce to 1 far from the GB.In the softer grain near the GB, the stress concentration factor is below 1 and increases to 1 far from the GB. e stress difference coefficient for both the models is higher at the GB and reduces to 1 far from the GB.
Figure 11 shows a different Mises stress distribution in both grains along path 2 in model III.Stress concentration exists at the GBs in both grains, and it is more significant in the harder grain.e elastic modulus difference in model III is the smallest; thus, the stress difference coefficient at the GB is the lowest and also reduces to 1 far from the GB. is indicates that the stress far away from the GB equals each other and is not affected by the grain orientation and GB.us, the Mises stress in the middle of grains may be equal if the adjacent grain size is larger enough.

Bicrystal with Type-II Grain Boundary. As shown in
Figure 12, a big difference of the Mises stress distribution on the GB of grain B for type-II bicrystal is observed, which intuitively demonstrates the misorientation effects of two adjacent grains.e Mises stress along path 1 increases with a large gradient initially and then reaches a plateau with a small gradient.e Mises stress on path 2 seems to be equal in the central region and has a small gradient close to the edges.e stress difference coefficient along path 2 in Figure 13 shows that a stress concentration exists at the GB and the most serious stress concentration appears in model II bicrystal, which has the largest elastic modulus difference along the load axis.e stress concentration is the smallest in model III bicrystal with the smallest elastic modulus difference along the load axis.e stress components σ 22 , σ 21 , and σ 23 on the GB along path 2 are further studied.e tensile stress component σ 22 is perpendicular to the GB plane, which could lead to the fracture of bicrystal at GB. e shear stress components σ 21    Advances in Materials Science and Engineering larger than the shear stress in the central region, and then, they all reduce gradually from the central region to edges and result in a larger shear stress at edges.e stress components are much higher in grain B than in grain A, which demonstrates the stress concentration in grain B.
e angle between shear stress and load direction is in a narrow range of [−1 °, 6 °], which indicates that the stress component σ 21 dominates the shear stress τ and the shear stress is almost parallel to the [110] direction.e symmetrical characteristic of stress components in grain A and grain B could be induced by the physical property symmetry of two grains with respect to tensile direction.
As shown in Figure 15(a), the tensile stress and shear stress along path 2 in grain A of model II are similar to those of model I, due to the same grain orientation of grain A in model II and model I. e tensile stress is flat with small values, and the shear stress is a little bigger than tensile stress in the central region and increases rapidly close to the edges.Due to the symmetry of grain A and the asymmetry of grain B, the stress components are approximately symmetry in grain A. As shown in Figure 15(b), the shear stress is about 2 times the tensile stress along path 2; thus, the shear failure appears.e angle between shear stress and tensile direction is in the range of [−150 °, −135 °], which is close to the angle between [012] direction and tensile direction on planes 1-3 or (-12-1).
A bigger difference between tensile stress and shear stress in grain A and grain B is observed in Figures 14 and 15, which indicates that a high stress concentration exists on the 8 Advances in Materials Science and Engineering GB of grain B, which has greater stiffness along the load direction in model I and model II.However, the situation is different in model III as shown in Figure 16, even the Mises stress is higher in grain B with greater stiffness, and the average tensile stress and shear stress in grain A are about 3.8 and 1.1 times than those in grain B, respectively.

Conclusions
A constitutive model has been presented to assess the elastic response of anisotropic behaviors of the austenitic stainless steel crystal.e FE method has been proposed to determine how the crystal orientation affects the stress state near the GB. e results indicate that the stress distribution near the GB depends strongly on the crystal orientation.
For bicrystals with the GB perpendicular or parallel to the tensile axis, elastic anisotropy causes additional stress in grain boundary regions and leads to the stress concentration on the GB of grains with larger stiffness along the load direction.e larger the elastic modulus difference of two neighboring grains along the tensile axis, the bigger the Mises stress difference appears at the GB.
With a GB perpendicular to tensile axes, a uniform Mises stress distribution in the region far from the GB in both grains is observed, which indicates that the grain orientation has little effects on the Mises stress distribution in this region.e Mises stress in the middle of adjacent grains may be equal if the grain size is larger enough.
e strain in the softer component is higher than that in the harder component; the strain near the GB is lower than Advances in Materials Science and Engineering the average strain of this grain in the softer component, while the strain near the GB is higher than the average strain of this grain in the harder component.e strain inconsistency between two neighboring grains and between the GB and the entire grain is affected by the mismatch of two neighboring grains.e larger the elastic modulus differences between two neighboring grains caused by misorientation, the larger the strain inconsistency in the bicrystals.

Figure 1 :
Figure 1: Type-I grain boundary: GB perpendicular to the tensile axis.

Figure 2 :
Figure 2: Type-II grain boundary: GB parallel to the tensile axis.

Figure 6 :
Figure 6: Locations of the two paths to determine the stress concentration factor.

Figure 7 :
Figure 7: Mises stress concentration factor on the GB along path 1.

Figure 11 :
Figure 11: Mises stress concentration factor and Mises stress difference coefficient of model III along path 2.

Figure 12 :
Figure 12: Mises stress distribution on the GB of grain B for type-II bicrystal.

Figure 16 :Figure 15 :
Figure 16: Stress components on model III GB: (a) grain A and (b) grain B.

Table 1 :
Crystal orientations considered for grains.

Table 2 :
Incompatibility bicrystal models with three types of GBs.

Table 3 :
e comparison of strain in each grain of the incompatibility models.