Simpliﬁed Model of a High Burnup Spent Nuclear Fuel Rod under Lateral Impact Considering a Stress-Based Failure Criterion

: The inventory of spent nuclear fuel (SNF) generated in nuclear power plants is continuously increasing, and it is very important to maintain the structural integrity of SNF for economical and efﬁcient management. The cladding surrounding nuclear fuel must be protected from physical and mechanical deterioration, which causes fuel rod breakage. In this study, the material properties of the simpliﬁed beam model of a SNF rod were calibrated for a drop accident evaluation by considering the pellet–clad interaction (PCI) of the high burnup fuel rod. In a horizontal drop, which is the most damaging during a drop accident of SNF, the stress in the cladding caused by the inertia action of the pellets has a great effect on the integrity of the fuel rod. The failure criterion for SNF was selected as the membrane plus bending stress through stress linearization in the cross-sections through the thickness of the cladding. Because the stress concentration in the cladding around the vicinity of the pellet–pellet interface cannot be simulated in a simpliﬁed beam model, a stress correction factor is derived through a comparison of the simpliﬁed model and detailed model. The applicability of the developed simpliﬁed model is checked through dynamic impact simulations. The developed model can be used in cask level analyses and is expected to be usefully utilized to evaluate the structural integrity of SNF under transport and in storage conditions.


Introduction
The spent nuclear fuel (SNF) discharged from a nuclear reactor is stored for several years to decades in wet storage and at an interim storage facility (ISF) for cooling before reprocessing or direct disposal [1]. The storage of SNF is consistently increasing worldwide, and dry storage facilities are operated in preparation for the saturation of wet storage [2]. It is essential that the SNF maintains a handleable condition without damage to the nuclear fuel rods before the final disposal or reprocessing. For this purpose, it is necessary to keep the structural integrity of the nuclear fuel rods from deteriorating, which could cause extensive damage to the cladding during the storage period [3][4][5]. Structural analysis using CAE (computer aided engineering) is important, but a single SNF cask contains dozens of fuel assemblies and thousands of fuel rods. Accordingly, the total number of elements of a transportation or storage cask is impractical if all the details of the fuel assemblies are to be modeled. The uncertainties lying in the physical-mechanical properties of SNF pose more difficulties to the structural evaluation of SNF. Choi et al. [6] verified the conditions under which delayed hydride cracking could occur for the integrity evaluation of the effect of hydride on the cladding during dry storage of SNF and compared the critical loads of delayed hydride cracking for the cladding with different hydrogen contents. In accordance with safety regulations, fuels must be protected from deterioration in physical-mechanical properties, resulting in damage to nuclear fuels [7,8]. However, even if the characteristics of SNF and the burnup are defined, the properties of SNF are uncertain when exposed to high temperatures and irradiation environments in a nuclear reactor. Jiang et al. [9,10] performed an extensive study on the impact of pellet-cladding interaction (PCI) on the fuel rod performance. Almomani et al. [11] studied the behavior of fuel rods about their buckling and bending resistance in drop accidents, using the properties of SNF and PCI as parameters, and recently analyzed the importance of the PCI on the structural integrity of SNF [12]. Lee et al. [13] studied the calibration of the simplified beam model properties of SNF by considering the PCI. However, these simplified model properties calibrated at a specific segment length as a reference have a large dependence on the length of the fuel rod being considered, showing that its applicability is limited. Thus, Lee and Kim [14] developed a new simplified model that was not limited by the length of the reference fuel rod segment.
In this study, a simplified model was developed for the horizontal drop condition, where the characteristics of the bending resistance of the fuel rods are dominant. It is known that this drop orientation is the most damaging to fuel rods. A new cladding failure criterion is proposed, which is based on linearized stress distribution through the thickness of the cladding. This criterion can account for the localized stress concentration and the through-thickness development of stress from which cladding failure can occur. The optimization for material property calibration of the simplified model was performed based on this stress-based criterion for two extreme interfacial conditions. The two extreme conditions are as follows: (1) a fully de-bonded condition is a gap between the pellets and the cladding when SNF was first stored, and (2) a fully bonded condition is where the pellets and the cladding are bonded through a physical-chemical reaction over time due to high temperatures and high pressures. To calibrate the properties of the simplified beam model, optimization of the properties of the simplified model is performed for the structural response of the detailed fuel rod model, using the equation of the same elasticplastic material as that of the Zircaloy-4 cladding. Development of the simplified model proceeds in four stages. The first is to select the interfacial bonding condition among the degradation characteristics of the SNF; the second is to analyze the fuel rod mechanical behaviors; the third is to develop a simplification model; and finally, an applicability test is conducted on the developed simplified model.

Target Fuel Rod
The design of the fuel rod assembly considered in this study refers to the CE 16 by 16 assembly model in Figure 1, and the design specifications are shown in Table 1. A nuclear fuel rod consists of fuel pellets and Zircaloy-4 cladding surrounding it. In a fresh fuel rod, there exists a small gap between the pellets and the cladding, but during reactor operation with high burnup, the gap is filled with a fused material or compressed through residual stress caused by the mechanical interaction of the two interfaces. Because the pellet-cladding interface (PCI) has a strong influence on the dynamic response of SNF, two extreme interfacial conditions were considered in this study. The target fuel rod considered in this study was a high burnup fuel (HBF) with a burnup rate of over 45 GWd/MTU.

Material Model of the Degraded HBF Cladding and Pellets
The properties of the HBF are functions of parameters such as the burnup and storage conditions, and important factors that influence the HBF properties include the hydrogen concentration, irradiation level, and temperature. Based on much experimental data, a detailed model was developed to predict the mechanical properties of the pellets and cladding as functions of the temperature, neutron fluence, hydrogen content, and burnup [16,17]. In this study, the parameters of the fuel were as follows: a discharge burnup of 60 GWd/MTU, a temperature of 573 K corresponding to a 5 year dry storage condition, a neutron fluence of 11.4 × 10 25 n/m 2 , and a total hydrogen concentration of 352 ppm. A detailed explanation about the choice of the above parameters can be found in [11]. The cladding material is described as an elastoplastic material that can be mathematically formulated using Hook's law and power law, as in following equations: where σ y is the yield stress (MPa); K is the strength coefficient; E is Young's modulus (MPa); n is the strain hardening exponent; m is the strain rate exponent; . ε is the strain rate(/s); σ e+p is the true stress; and UE is the uniform elongation. The parameters in the model for the reference fuel rods considered in this study are summarized in Table 2. At 573 K, with only circumferential hydride, the fracture strain of the target HBF Zircaoly-4 cladding is known to be 1~4% [18], and the failure stress corresponding to a 1% plastic strain is 749 MPa. In this study, a stress of 749 MPa was selected as the failure criterion of the fuel cladding [19,20].

Determination of the Failure of the Fuel Cladding
The failure strain of irradiated Zircaloy-4 cladding was obtained by experimental measurement using the uniaxial tension test or pressurized tube test with defueled cladding specimens. In both tests, an almost uniform tensile stress and strain are induced through the thickness of the cladding. This stress state and deformation pattern are quite different from those observed in actual fuel rods with fuel pallets. When a fuel rod is subject to bending, a very complicated 3D stress state is generated within the vicinity of the pelletpellet interface shown in Figure 2. A very localized stress concentration is observed in the area where the edges of the pellets contact the cladding. If the failure strain criteria obtained from the uniaxial tension test or biaxial tube test is applied, very early failure is predicted at the beginning stage of loading with premature through-thickness stress or strain development. The localized contact stress or strain is self-limiting and is not a good candidate for the cladding failure criteria.
where is the membrane stress (MPa); t is the thickness of the section; along the path; is the coordinate along the path; and is the uniform int between endpoints. In Equation (5), and are the bending stress at th the section in Figure 3. and are the total thickness stress at the en B.  Therefore, a new cladding failure criterion is proposed in this work, which can account for the localized stress concentration and the through-thickness stress development. The failure of the cladding is determined by the membrane plus bending stress generated through the thickness of the cladding, which can be calculated by a process called stress linearization along the stress classification line shown in Figure 3. The fuel rod failure criterion was based on the aforementioned stress of 749 MPa of the membrane plus bending contribution of the von Mises stress of the cladding.

Strategy of the Model Simplification
The purpose of developing the simplified model is to develop a fuel r can be used in transport or storage cask analyses to evaluate the structu nuclear fuel rods and to calculate the fuel damage ratio under accident cond the cask drop impact, the fuel rods receive the inertial load from the impa and the pinch loads that occur due to collision with other fuel rods and the cask. In the work of Lee and Kim [14], a simplified model was developed t Stress linearization is a process of decomposing the total stress generated in an arbitrary cross-section of a structure into three stress components, as shown in Figure 4, and evaluated by linearizing the stress distribution in the thickness direction of the structure [21]. The membrane values of the stress components are computed using Equation (3), and it does not transfer any moments. Assuming that the whole thickness between point A and point B in Figure 5 are divided uniformly into n intervals, the thickness membrane stress can be calculated as the average thickness stress by solving Equation (4). The linear bending values of the stress components at the endpoints of the sections are computed using Equation (5), and thickness surface stress does not have a peak stress contribution. Hence, the thickness bending stress is that of the sum of the thickness membrane and bending stresses at the total stress at each point, as shown in Equation (6), and interpolates the bending stress linearly [22,23].
where σ m is the membrane stress (MPa); t is the thickness of the section; σ is the stress along the path; x is the coordinate along the path; and n is the uniform intervals divided between endpoints. In Equation (5), σ b A and σ b B are the bending stress at the endpoints of the section in Figure 3. σ tA and σ tB are the total thickness stress at the endpoints A and B.

Strategy of the Model Simplification
The purpose of developing the simplified model is to develop a fuel rod model that can be used in transport or storage cask analyses to evaluate the structural integrity of nuclear fuel rods and to calculate the fuel damage ratio under accident conditions. During the cask drop impact, the fuel rods receive the inertial load from the impact acceleration and the pinch loads that occur due to collision with other fuel rods and the baskets in the cask. In the work of Lee and Kim [14], a simplified model was developed that considered the deflection of the fuel rods under a pure bending moment, and it was demonstrated that the simplified beam model is not dependent on the length of the reference fuel rod segment for calibration. A 1% plastic strain calculated from the curvature of the fuel rod was used as the failure criterion of the cladding. However, this criterion cannot account for the effect of the localized stress concentration in the vicinity of the pellet-pellet interface, which obviously promotes the failure of the cladding.
In this work, the framework of the model calibration proposed by Lee and Kim [14] is modified using a new cladding failure criterion. The strategy of the model calibration is as follows: 1) the fuel rod is simplified into a hollow beam with an isotropic material; 2) the flexural rigidity of the simplified beam should be matched with that of the detailed model so that the deflection and pinch load can be estimated accurately; 3) because an isotropic beam cannot express the localized contact stress, a factor named the stress correction factor is proposed so that the same failure criterion can be applied to the simplified model. To implement a model calibration procedure that complies to the above-mentioned strategy, a reference fuel rod segment of 56.7 mm in length was chosen, and the critical bending moment was calculated using static analyses with the detailed model. The critical bending moment is the one that produces 749 MPa of membrane plus bending stress through a thickness of the cladding. The parameters of the simplified beam were calibrated in such a way that the deflection of the simplified model is matched with that of the detailed model under this critical bending moment. This procedure is applied to the fully bonded fuel rod and fully de-bonded fuel rod separately to account for the effect of the PCI.

Strategy of the Model Simplification
The purpose of developing the simplified model is to develop a fuel rod model that can be used in transport or storage cask analyses to evaluate the structural integrity of nuclear fuel rods and to calculate the fuel damage ratio under accident conditions. During the cask drop impact, the fuel rods receive the inertial load from the impact acceleration and the pinch loads that occur due to collision with other fuel rods and the baskets in the cask. In the work of Lee and Kim [14], a simplified model was developed that considered the deflection of the fuel rods under a pure bending moment, and it was demonstrated that the simplified beam model is not dependent on the length of the reference fuel rod segment for calibration. A 1% plastic strain calculated from the curvature of the fuel rod was used as the failure criterion of the cladding. However, this criterion cannot account for the effect of the localized stress concentration in the vicinity of the pellet-pellet interface, which obviously promotes the failure of the cladding.
In this work, the framework of the model calibration proposed by Lee and Kim [14] is modified using a new cladding failure criterion. The strategy of the model calibration is as follows: (1) the fuel rod is simplified into a hollow beam with an isotropic material; (2) the flexural rigidity of the simplified beam should be matched with that of the detailed model so that the deflection and pinch load can be estimated accurately; (3) because an isotropic beam cannot express the localized contact stress, a factor named the stress correction factor is proposed so that the same failure criterion can be applied to the simplified model. To implement a model calibration procedure that complies to the above-mentioned strategy, a reference fuel rod segment of 56.7 mm in length was chosen, and the critical bending moment was calculated using static analyses with the detailed model. The critical bending moment is the one that produces 749 MPa of membrane plus bending stress through a thickness of the cladding. The parameters of the simplified beam were calibrated in such a way that the deflection of the simplified model is matched with that of the detailed model under this critical bending moment. This procedure is applied to the fully bonded fuel rod and fully de-bonded fuel rod separately to account for the effect of the PCI.

Detailed Finite Element Model and Results of the Static Analyses
The detailed fuel rod model used in this study is shown in Figure 5. A half model of a fuel rod segment was constructed using the symmetry condition. A total of 24,768 elements of an eight-node linear brick and incompatible modes were used in this model. To see the stress distribution between the inner surface of the cladding contacting the corner of the pellets in detail, a dense mesh was used in the corresponding cladding area shown in Figure 6. In the fully bonded condition, the outer surface of the pellets and the inner surface of the cladding were given the tie constraints, and the fully de-bonded condition was simulated using the frictionless contact conditions. The detailed modeling conditions for the two interfacial conditions are shown in Table 3. The bending moment applied to a fuel rod is given in Figure 7. A more detailed explanation of interfacial bonding and its effect can be found in [9,10]. The hoop stress distribution caused by the rod internal pressure was not considered in this work. It can be obtained using code such as FRAPCON [24] and will be investigated in our future work. elements of an eight-node linear brick and incompatible modes were used To see the stress distribution between the inner surface of the cladding c corner of the pellets in detail, a dense mesh was used in the corresponding shown in Figure 6. In the fully bonded condition, the outer surface of the p inner surface of the cladding were given the tie constraints, and the ful condition was simulated using the frictionless contact conditions. The deta conditions for the two interfacial conditions are shown in Table 3. The ben applied to a fuel rod is given in Figure 7. A more detailed explanation bonding and its effect can be found in [9,10]. The hoop stress distribution rod internal pressure was not considered in this work. It can be obtained us as FRAPCON [24] and will be investigated in our future work.    Static analyses for the two interfacial bonding conditions were perfor critical moment that generates a membrane plus bending stress of 749 MP thickness of cladding was derived using linear interpolation. The results o at which the failure occurs are 31,413 N·mm for the fully bonded conditio elements of an eight-node linear brick and incompatible modes were used To see the stress distribution between the inner surface of the cladding c corner of the pellets in detail, a dense mesh was used in the corresponding shown in Figure 6. In the fully bonded condition, the outer surface of the p inner surface of the cladding were given the tie constraints, and the ful condition was simulated using the frictionless contact conditions. The deta conditions for the two interfacial conditions are shown in Table 3. The ben applied to a fuel rod is given in Figure 7. A more detailed explanation bonding and its effect can be found in [9,10]. The hoop stress distribution rod internal pressure was not considered in this work. It can be obtained us as FRAPCON [24] and will be investigated in our future work.    Static analyses for the two interfacial bonding conditions were perfor critical moment that generates a membrane plus bending stress of 749 MP thickness of cladding was derived using linear interpolation. The results o at which the failure occurs are 31,413 N·mm for the fully bonded conditio  corner of the pellets in detail, a dense mesh was used in the corresponding shown in Figure 6. In the fully bonded condition, the outer surface of the inner surface of the cladding were given the tie constraints, and the fu condition was simulated using the frictionless contact conditions. The deta conditions for the two interfacial conditions are shown in Table 3. The ben applied to a fuel rod is given in Figure 7. A more detailed explanation bonding and its effect can be found in [9,10]. The hoop stress distribution rod internal pressure was not considered in this work. It can be obtained u as FRAPCON [24] and will be investigated in our future work.     Static analyses for the two interfacial bonding conditions were performed, and the critical moment that generates a membrane plus bending stress of 749 MPa through the thickness of cladding was derived using linear interpolation. The results of the moment at which the failure occurs are 31,413 N·mm for the fully bonded condition and 22,410 N·mm for the fully de-bonded condition. The deformed shapes and stress linearization are shown in Figure 8.
the curvature of the fuel rods, and they were 31,919 N·mm for the fully bonded case and 27,975 N·mm for the fully de-bonded case. It was observed that the failure criterion used in this work predicts almost the same critical moment for the fully bonded fuel rods, while smaller critical moments are predicted for the fully de-bonded case. This result implies that the linearized stress-based criterion of this work takes more account of the localized stress compared to the curvature-based criteria of [14].

Parameter Calibration for the Simplified Model
The material of the simplified model was modeled as an elastic-plastic material described with the same equation used for Zircaloy-4 in Section 2. Young's modulus of the simplified model was calculated using Equation (7) with the analyses results of the detailed model as follows: where is Young's modulus of the simplified beam model (GPa); is the bending moment (N·mm) of a small, arbitrary value under which the fuel rod shows a linear In the work of Lee and Kim [14], the critical bending moments were calculated from the curvature of the fuel rods, and they were 31,919 N·mm for the fully bonded case and 27,975 N·mm for the fully de-bonded case. It was observed that the failure criterion used in this work predicts almost the same critical moment for the fully bonded fuel rods, while smaller critical moments are predicted for the fully de-bonded case. This result implies that the linearized stress-based criterion of this work takes more account of the localized stress compared to the curvature-based criteria of [14].

Parameter Calibration for the Simplified Model
The material of the simplified model was modeled as an elastic-plastic material described with the same equation used for Zircaloy-4 in Section 2. Young's modulus of the simplified model was calculated using Equation (7) with the analyses results of the detailed model as follows: where E s is Young's modulus of the simplified beam model (GPa); M is the bending moment (N·mm) of a small, arbitrary value under which the fuel rod shows a linear behavior; L is the length of the fuel rod segment (mm); I is the area moment of the inertia for the simplified beam model (mm 4 ); and δ d is the vertical deflection of the detailed fuel rod under M. The area moment of the inertia I is determined from the calibration parameters, and as a result, the E s also varies with the choice of calibration parameters. For the calibration of the properties of the simplified beam model, ISGITH [25], the optimization program, was used to optimize the calibration parameters to minimize the discrepancy in the behavior of the simplified model and the detailed model. The strength coefficient (MPa), the strain hardening exponent n, and the thickness of the hollow shaft beam T were considered as the parameters of the calibration. The calibration procedure was formulated into an optimization problem as follows: Fine K, n and T such that where U E is the displacement of the fuel rod in the elastic region, and U C is the displacement of the fuel rod under the critical moment. The subscripts d and s denote the detailed model and the simplified model, respectively. The optimization procedure is depicted in Figure 9.
behavior; L is the length of the fuel rod segment (mm); is the area moment of the i for the simplified beam model (mm 4 ); and is the vertical deflection of the detaile rod under . The area moment of the inertia is determined from the calib parameters, and as a result, the also varies with the choice of calibration parame For the calibration of the properties of the simplified beam model, ISGITH [25 optimization program, was used to optimize the calibration parameters to minimiz discrepancy in the behavior of the simplified model and the detailed model. The str coefficient (MPa), the strain hardening exponent , and the thickness of the hollow beam were considered as the parameters of the calibration. The calibration proc was formulated into an optimization problem as follows: where is the displacement of the fuel rod in the elastic region, and displacement of the fuel rod under the critical moment. The subscripts and d the detailed model and the simplified model, respectively. The optimization proced depicted in Figure 9. and . The material properties model is then inputted into the ABAQUS material and the displacement of the simplified model required for calculating the optimiz objective functions is calculated. The displacement output is the input to the ISIGTH the design variables are updated until the objective function is minimized throug adaptive simulated annealing (ASA) method [26]. A total of 2,784 iterations for the bonded condition and 242 iterations for the fully de-bonded condition were perfo before optimal solutions were obtained. The optimization results corresponding to PCI condition are shown in Table 4. The material model derived through the optim design variables was applied to the simplified model, and the deformed shape o detailed model and the simplified model are compared in Figure 10. Figure 11 show moment-displacement curve. It is shown that the deflection of the simplified m agrees well with those of the detailed models at the evaluation points ( , ), alth a slight discrepancy is observed in the plastic region. In the comparison of the mem plus bending stress of the von Mises stress in the simplified model, which is the fu failure criterion of the detailed model proposed in this study, it has a deviation of 4 the fully bonded condition and 67% in the fully de-bonded condition. An idealized beam cannot express the localized stress concentration observed in actual fuel rods The structural response of the simplified beam model was calculated using the ABAQUS/Standard (V2016, Dassault Systèmes Simulia Corp, Warwick, RI, USA). Using MATLAB (R2016, Mathworks, Natick, MA, USA), the stress-strain relationship model is created by inputting the Zircaloy-4 property equation and calibration parameters K, n, and T. The material properties model is then inputted into the ABAQUS material table, and the displacement of the simplified model required for calculating the optimization objective functions is calculated. The displacement output is the input to the ISIGTH, and the design variables are updated until the objective function is minimized through the adaptive simulated annealing (ASA) method [26]. A total of 2,784 iterations for the fully bonded condition and 242 iterations for the fully de-bonded condition were performed before optimal solutions were obtained. The optimization results corresponding to each PCI condition are shown in Table 4. The material model derived through the optimized design variables was applied to the simplified model, and the deformed shape of the detailed model and the simplified model are compared in Figure 10. Figure 11 shows the moment-displacement curve. It is shown that the deflection of the simplified models agrees well with those of the detailed models at the evaluation points (U E , U C ), although a slight discrepancy is observed in the plastic region. In the comparison of the membrane plus bending stress of the von Mises stress in the simplified model, which is the fuel rod failure criterion of the detailed model proposed in this study, it has a deviation of 44% in the fully bonded condition and 67% in the fully de-bonded condition. An idealized Euler beam cannot express the localized stress concentration observed in actual fuel rods, and this discrepancy in the stress was predictable. In this study, the error ratio was used to calculate the stress correction factor to compensate for the stress results in the simplified models. Before that, it is necessary to check whether the von Mises stress result of the simplified model has a consistent error ratio for the various segment lengths in the two interface bonding conditions for the failure criterion of 749 MPa. this discrepancy in the stress was predictable. In this study, the error ratio was used to calculate the stress correction factor to compensate for the stress results in the simplified models. Before that, it is necessary to check whether the von Mises stress result of the simplified model has a consistent error ratio for the various segment lengths in the two interface bonding conditions for the failure criterion of 749 MPa.   this discrepancy in the stress was predictable. In this study, the error ratio was used to calculate the stress correction factor to compensate for the stress results in the simplified models. Before that, it is necessary to check whether the von Mises stress result of the simplified model has a consistent error ratio for the various segment lengths in the two interface bonding conditions for the failure criterion of 749 MPa.

Discussions on the Failure Limit Criteria of the Cladding with Various Segment Lengths
A simplified model considering the bending load was developed for a fuel rod segment with a length of 56.71 mm containing five fuel pellets. The detailed and simplified models of five fuel rod segments with different lengths were constructed and analyzed to investigate whether the calibrated material parameters and error ratio are applicable to fuel rods with different lengths. The beam parameters found in Section 3.3 were applied to all the simplified models. The 749 MPa membrane plus bending stress selected as the failure criterion at the critical moment and displacement of a simplified model with isotropic materials was compared for the two interfacial bonding conditions and summarized in Tables 5 and 6. Table 5. Application of the simplified model parameters for the bonded condition to different lengths of fuel rod segments.

Cases
Number   Table 6. Application of the simplified model parameters for the de-bonded condition to different lengths of fuel rod segments.  The structural response of the segments with various lengths shows almost identical deviation levels in both interfacial bonding conditions. At the critical moment, the maximum displacement deviation of the detailed model and the simplified model was 1.8% in the fully bonded condition and 3.1% in the fully de-bonded condition. When comparing the membrane plus bending stress, the fully bonded condition showed a consistent 44% deviation, and the fully de-bonded condition showed a consistent 68% deviation. Thus, it is shown that the correlation of the behavior of models as well as the deviation percentages in the linearized stress are insensitivity to the choice of the length of the reference fuel rod segment for calibration. From this result, it was confirmed that the error ratio of 44% and 68% can be used to find the stress correction factor to compensate the stress calculation results using the simplified beam model, regardless of the lengths of the fuel rods being considered. The stress correction factor of the simplified model for the fully bonded condition is 1.78 and 3.06 for the fully de-bonded condition. It should be noted that these factors are applicable only to the fuel rods with the same specification with the one in the current study, which is subject to pure bending. The applicability of these factors and the factors those have effects on them will be investigated in our future study.

Models for Impact Simulation and Impact Conditions
Dynamic impact simulations were performed using the detailed models and simplified models of fuel rods to check the applicability of the calibrated beam parameters. The procedure in this section is similar to the procedure introduced in [14]. In this section, the effect of different failure criteria and the usefulness of the stress correction factor were mainly investigated. The detailed model and simplified model for the dynamic impact analysis for the two interface bonding conditions are shown in Figure 12. A fuel rod with a length of 113.41 mm that included 10 pellets was used as a reference model to examine the dynamic response during the drop impact. The fuel rod segment was constrained to Zircaloy-4 blocks at both ends with a coupling constraint. The fuel rod and Zircaloy-4 block impact a rigid surface in the horizontal orientation, as shown in Figure 12.

Models for Impact Simulation and Impact Conditions
Dynamic impact simulations were performed using the detailed models and simplified models of fuel rods to check the applicability of the calibrated beam parameters. The procedure in this section is similar to the procedure introduced in [14]. In this section, the effect of different failure criteria and the usefulness of the stress correction factor were mainly investigated. The detailed model and simplified model for the dynamic impact analysis for the two interface bonding conditions are shown in Figure  12. A fuel rod with a length of 113.41 mm that included 10 pellets was used as a reference model to examine the dynamic response during the drop impact. The fuel rod segment was constrained to Zircaloy-4 blocks at both ends with a coupling constraint. The fuel rod and Zircaloy-4 block impact a rigid surface in the horizontal orientation, as shown in Figure 12. ABAQUS/Explicit was used for the dynamic impact simulations, and the drop heights in the range of 0.3 m to 1.5 m were transformed into the initial velocity using Equation (9) as follows: where v is the initial velocity (mm/s); g is the gravitational acceleration mm/s 2 , and h is the drop heights (mm).

Results and Discussion
Dynamic impact simulations were performed considering four drop heights for both of the interfacial conditions. The results of the maximum deflection of the beam and membrane plus bending stress are summarized in Tables 7 and 8. It was observed that the deviation of the maximum deflection of the detailed model and the simplified model is 9.9% for the fully bonded condition and 19.0% for the fully de-bonded condition. During impact, a shear force is exerted to the fuel rods and a bending moment resulting from the inertia of mainly the fuel pellets, and the fully de-bonded fuel rods have less shear resistance than that of the fully bonded fuel rods. It is believed to be the cause of the bigger deviation of the deflection results in the fully de-bonded case, and its results are consistent with those reported in [14]. At the same drop height of 1 m, the maximum displacement deviation of the detailed model of the two interfacial bonding conditions is 25%, confirming that the PCI has a significant effect on the fuel rod behavior in a drop accident.  The maximum membrane plus bending stresses during the impacts are summarized in Tables 7 and 8. The numbers in parenthesis are those corrected by multiplying the stress correction factor derived in Section 3 to the results of the simplified models. When the stress compensation factors were applied, the deviation of the stress predictions in the detailed model and the simplified model ranged up to 31% in the fully bonded case and up to 18% in the fully de-bonded case. However, those errors diminish to almost 0 as the drop height approaches the critical drop height, where the membrane plus bending stress becomes 749 MPa. For the fully bonded case, the critical drop height is predicted around 1.5 m, and for the fully de-bonded case, it is expected to be bigger than 1.5 m. It is very notable that the critical drop height is bigger in the fully de-bonded case than in the fully bonded case, which is different from the conclusions reported in other works [11,13,14].
One reason for this phenomenon is that the von Mises stress is used in the failure criteria for the stress linearization. The von Mises stress and the Tresca stress are governed by the shear component in a stress state and are suitable for the prediction of ductile failure. In the fully bonded condition, a much bigger shear stress is developed in the cladding within the vicinity of the pellet-clad interface than in the fully de-bonded case. Another reason is that the stress is more localized in the fully bonded case. Stress is concentrated in the narrow region around the pellet-pellet interface, which is relatively flexible compared to the other area of the cladding, which is bonded to the pellets. This is clearly seen in Figure 8. It can be seen that the stress is more uniformly distributed to the length of the cladding in the fully de-bonded case. The results of this section might be controversial, and it is believed that the failure criteria introduced in this study can be reasonably and conservatively applied to the prediction of cladding failure as long as its failure mechanism is a ductile failure. The validity of this statement can be checked if experimental data on fuel rod failure become available.
The deflection time history is summarized in Table 9. The solid lines with legend "solid" are the responses of the detailed model. It is seen that the simplified models are more flexible than the detailed model under the dynamic load, and this trend is more obvious in the fully de-bonded case. As discussed above, it is due to the difference in the shear resistance of the two interfacial conditions.

Discussion and Conclusions
In this study, a new failure criterion for SNF cladding was introduced. The membrane plus bending stress of the von Mises stress was used for the failure criterion,

Discussion and Conclusions
In this study, a new failure criterion for SNF cladding was introduced. The membrane plus bending stress of the von Mises stress was used for the failure criterion,

Discussion and Conclusions
In this study, a new failure criterion for SNF cladding was introduced. The membrane plus bending stress of the von Mises stress was used for the failure criterion, and the value is 749 MPa for the reference fuel rod considered in this work. The simplified

Discussion and Conclusions
In this study, a new failure criterion for SNF cladding was introduced. The membrane plus bending stress of the von Mises stress was used for the failure criterion, and the value is 749 MPa for the reference fuel rod considered in this work. The simplified

Discussion and Conclusions
In this study, a new failure criterion for SNF cladding was introduced. The membrane plus bending stress of the von Mises stress was used for the failure criterion, and the value is 749 MPa for the reference fuel rod considered in this work. The simplified

Discussion and Conclusions
In this study, a new failure criterion for SNF cladding was introduced. The membrane plus bending stress of the von Mises stress was used for the failure criterion, and the value is 749 MPa for the reference fuel rod considered in this work. The simplified

Discussion and Conclusions
In this study, a new failure criterion for SNF cladding was introduced. The membrane plus bending stress of the von Mises stress was used for the failure criterion, and the value is 749 MPa for the reference fuel rod considered in this work. The simplified

Discussion and Conclusions
In this study, a new failure criterion for SNF cladding was introduced. The membrane plus bending stress of the von Mises stress was used for the failure criterion, and the value is 749 MPa for the reference fuel rod considered in this work. The simplified

Discussion and Conclusions
In this study, a new failure criterion for SNF cladding was introduced. The membrane plus bending stress of the von Mises stress was used for the failure criterion, and the value is 749 MPa for the reference fuel rod considered in this work. The simplified beam models were developed based on parameter calibration using design optimization. The purpose of these models is to predict the failure of fuel rods in cask level analyses. To account for the effect of the PCMI, two separate models were developed for the fully bonded pellet-clad interface condition and the fully de-bonded interface condition. Because the beam models cannot simulate the concentrated stress state, stress correction factors were derived, which can be used in the prediction of the maximum membrane plus bending stress in the simplified models. The applicability of the simplified models and the stress correction factors were verified in the dynamic impact analyses, considering various drop heights. The following conclusions were drawn. - In a comparison with the previous work of the authors, it is shown that the stressbased failure criterion is more conservative than the strain-based criteria based on the curvature of the deflected fuel rod. - The stress correction factors are very effective in the failure prediction of the simplified beam models. However, the specific values of these factors can vary due to a number of factors, such as the shape of the fuel pellets, the material properties of the cladding and pellets, the gaps between the pellets and cladding, etc. - In the dynamic impact simulations, the stress prediction using the simplified model and the stress correction factor show good accuracy as the drop height approaches the critical drop height. This is natural, because the simplified models were developed focusing on the failure point of the fuel rods. -Although the fully bonded fuel rod deflects less than the fully de-bonded case, earlier failure is predicted with the stress-based failure criterion in the dynamic simulations. It is mainly due to the development of the shear stress in the interface of the fuel pellets and cladding in the fully bonded case. The validity of the failure criterion of the fuel cladding under a complicated 3D stress state should be confirmed with experimental data. - The discrepancy of simplified model behavior can be reduced by adopting material models with more parameters to be calibrated together with more powerful fitting methods. It is expected that the correlation in the small loading situation can be improved significantly with such efforts.
Author Contributions: S.L. was the principal investigator of this research and contributed to the research planning, verification of the results and the manuscript writing. S.K. contributed to the survey of the relevant literature, scientific computing using ABAQUS and preparation of the data, figures and tables in the manuscript and the manuscript writing. All authors have read and agreed to the published version of the manuscript.