Special Issue

## Application and Policy Evaluation of New Energy and Materials

View this Special Issue

Research Article | Open Access

Volume 2021 |Article ID 5330068 | https://doi.org/10.1155/2021/5330068

Wei Wang, Ping Wang, Xuesong Liu, Zhibo Dong, Hongyuan Fang, "Mathematical Model for Charpy Impact Energy of V-Notch Specimens", Advances in Materials Science and Engineering, vol. 2021, Article ID 5330068, 7 pages, 2021. https://doi.org/10.1155/2021/5330068

# Mathematical Model for Charpy Impact Energy of V-Notch Specimens

Accepted15 Nov 2021
Published02 Dec 2021

#### Abstract

Firstly, by analyzing the response of Charpy V-notch specimen impacted by pendulum, the relationship between specimen geometry, material properties, and impact energy is established and simplified, and the mathematical model for evaluating impact energy of specimens with different sizes is established. Then, the effectiveness of the model through a series of impact tests is verified. Theoretical analysis and experimental results show that the relationship between ligament length and impact energy is quadratic, while the relationship between ligament thickness and impact energy is linear. In the derivation process, the intrinsic impact toughness is used to evaluate the toughness of materials. The mathematical model makes it possible to evaluate the impact energy of specimens with different sizes and provides a theoretical basis for evaluating the impact resistance of structures.

#### 1. Introduction

Impact energy characterizes the ability of a structure to consume energy through deformation and fracture under impact load. The impact energy Ak, a value measured by the Charpy impact test, is mainly used to control the quality of metallurgical and thermal processing products and determine the ductile-brittle transition temperature of materials [1]. However, the determination of the Ak value is required to use specimens with a certain standard size. The comparison between the impact energy of specimens of different sizes or nonstandard sizes is disapproved. Since it partly depends on the geometry of specimens, impact energy cannot be regarded as a mechanical property of the material itself. Meanwhile, the relationship between impact energy and specimen size is not clear, which cannot be applied to evaluate the impact bearing capacity of structures whose sizes are not the same as standard specimens.

Several methods have been proposed to isolate the impact of collision energy on geometry. Traditionally, impact toughness αk, the quotient of the impact energy, and the cross-sectional area of the notch are used as the parameters to evaluate the toughness of material. However, for material with plastic deformation, energy consumption does not occur on the notched section only, and energy consumption across the section is not uniform, indicating that the impact toughness is not relevant and accurate enough as a mechanical parameter of material. Mikhail et al. [2] studied the influence of different specimen sizes on the impact test results and presented a procedure for transforming data from subsize specimens to data from full-size specimens. Sokolove et al. [3] studied the correlation between the ductile-brittle transition temperatures measured with impact specimens of different sizes, developed a procedure for transforming data from subsize and the full-size specimens, and turned out a good agreement. Lucon et al. [4] found that there is a linear relation between the impact energy and the specimen size by studying the results of the impact test between the small-size specimen and the full-size specimen. Schill et al. [4] studied the correlation between the impact energy of a certain ferritic steel small-size specimen and a full-size standard impact specimen and evaluated the estimation effects of various conversion formulas. Konopik et al. [5] established the conversion relation between the impact test results of small-size impact specimens and full-size standard impact specimens based on the test results. Yang et al. [6] used GA-NN combined model to predict the Charpy impact energy, and R. Muscat et al. [7] used membership function modeling to predict the Charpy impact energy of heat-treated steel more accurately. Chaouadi and Fabry [8] did a preliminary study on the relation between the impact energy and the size of the simply supported rectangular cross-section structure under impact load, and the result of the study is that the impact energy is proportional to the section modulus in bending. Sreenivasan and Mannan [9] studied the impact energy of specimens with different cross-sectional shapes. He further verified that the impact energy is proportional to the section modulus in bending of rectangular cross-section structures and found that the regularity is not significant for nonrectangular cross-sections. Inspired by summarizing the existing research, we try to find the functional relationship between the impact sample size and the impact energy measurement results, so as to estimate the impact energy of any size sample.

In this study, based on the load-bearing form of the V-notch specimens subjected to the pendulum lateral impact, and the response process of structural deformation and failure, a mathematical model containing impact energy, mechanical properties, and geometrical parameters of the structure was established. In this way, the influence of specimen geometry to impact energy is extracted, and a new mechanical property to describe the impact load-bearing capability of a material is defined. Finally, the regression analyses based on impact test results of two kinds of steel were performed to verify the rationality of the model.

#### 2. Mathematical Model to Evaluate Impact Energy

##### 2.1. Dynamic Response Process of V-Notched Specimen Subjected to Lateral Impact

The deformation mode of V-notch specimen subjected to lateral impact, as shown in Figure 1, can be regarded as a dynamic 3-point bending. After impacted by the pendulum, the dynamic response process of the specimen mainly includes three parts before completely broken: elastic deformation, plastic deformation, and crack propagation. Based on the plane assumption and the analysis of the stress-strain state of the bending deformation, during the elastic deformation, the kinetic energy of the pendulum is transformed into elastic potential energy, which is distributed in the entire specimen. In the plastic deformation stage, the plastic deformation locally occurs near the section which contains the impact position. The point at the tip of the V-notch has the maximum tensile stress and strain. When the bending deformation develops to a certain extent, cracks should occur at this point. After crack initiation, the crack shall propagate in the direction most conducive to its propagation until the sample is completely broken.

Through the analysis of the dynamic response process of the specimen after the impact, it can be known that the resistance capability of the specimen structure to impact load is related to the bending resistance in the elastic and plastic deformation stage and the energy absorption during the crack propagation.

##### 2.2. Construction of Mathematical Model of V-Notch Impact Energy

Based on the above analysis, it can be concluded that the impact energy consists of three parts: the elastic deformation energy Ee, the plastic deformation energy of the structure Ep, and the energy consumed by crack propagation in fracture section Ec. The impact energy of the specimen can be expressed as the sum of these energies, that is,

In the bending process, energy consumed in the deformation process is the product of bending moment and deflection angle increment. Meanwhile, the stress of the material will change due to strain strengthening effect and strain rate effect, which will lead to the real-time change of the bending moment which resists the deformation. Therefore, bending moment can be regarded as a function of deflection angle M(θ). The energy consumed by crack propagation is the sum of the energy consumed by each area unit on the fracture section. The energy consumed on the area unit can be expressed as a function of position γ(x, y). To sum up, each item at the right side of equation (1) can be expressed as

Among them, is the real-time bending moment of the specimen during bending deformation; is the energy consumed on the area units at different locations of the cross-section when the crack grows; is the rotation angle of the specimen at the end of the elastic deformation; is the angle of the specimen when crack initiated at the backside of the specimen; is the area of the crack propagation section.

Thus, equation (1) can be expressed as

In the elastic deformation stage, for the simply supported rectangular cross-section structure, the bending moment of the cross-section is

Among them, is the maximum stress on the section, in the elastic stage; its size is related to the bending deformation; is the width of the cross-section of the notch center; is the thickness of the cross-section of the notch center.

Substituting equation (6) to equation (2) gives

In the elastic stage, the maximum stress and strain on the notch cross-section increase linearly with the deflection angle until the elastic limit is reached, so it can be obtained from equation (7) that

Among them, is the yield strength of the material.

According to mechanics of materials, the deflection angle at the end of elastic deformation stage is

Among them, is the span between the support ends of the simply supported structure. is young’s modulus of the material.

Substituting equation (9) into equation (8), the energy dissipation in the elastic stage iswhere is the strain when the material reaches the elastic limit.

Regarding αe as an elastic index reflects the ability of material to absorb energy in elastic stage,

Substituting equation (11) into equation (10) gives

When the maximum stress on the cross-section reaches the yield strength of the material, the elastic deformation generally turns into plastic deformation. After entering the plastic deformation stage, considering strain strengthening effects and strain rate effect, mechanical properties of the material show a high degree of nonlinear characteristics during the dynamic deformation process [10]. Therefore, the real bending moment at the bending section changes in real time, behaving as a function of deflection angle θ. The real-time bending moment of cross-section in plastic deformation stage can be written as a product of correction coefficient kr(θ) and the ultimate plastic moment of cross-section Mp, that is,

The ultimate plastic bending moment Mp is calculated by the static yield strength of the material and the geometric size of the notched section, reflecting the ultimate bending moment when the material on the section theoretically yields completely under static bending. Mp does not change in the deformation process. The actual change of bending moment during deformation is considered by multiplying the correction coefficient function kr(θ).

Substituting equation (13) to equation (3), we can get

It is difficult to obtain the analytical form of the strain strengthening effect and strain rate effect of the material in the process of plastic deformation, but the integral result can be replaced by the average value multiplied by the total deformation, that is,

Among them, is the mean value in the process of change .

Substituting equation (15) to equation (14), it gives

For a rectangular cross-section structure, the ultimate moment of the section is

Different from elastic deformation, the plastic deformation of the specimen after impact is limited to a local position. The plastic area of a rectangular cross-section specimen is assumed to be two symmetrically distributed triangles, as shown in Figure 2.

When plastic deformation occurs, the relation between deflection angle increment and curvature iswhere is the largest length of plastic area, as shown in Figure 2.

It should be noted that the local length of plastic area lh shall change during the deformation. However, for mathematical simplicity, some fixed values of l are recommended. In the case of rectangular cross-section beam, lh roughly equals to the thickness of bending section [11,12].

The relation between the maximum strain in the bending section and the curvature iswhere is the distance from the notch to the neutral surface.

If εmax increases to fracture strain εf, crack initiation occurs, and fracture strain εf is a mechanical property of material itself. Let εmax = εf; from equation (18) and (19), the deflection angle increment when crack initiation occurs can be derived

From equations (16)∼(20), the energy consumed by plastic deformation can be derived

Let

So, equation (21) can be written as

Among them, is the intrinsic impact toughness, which is defined by equation (22) to reflect the toughness of material itself in the dynamic process.

Similarly, the average value is used to simplify the expression of energy consumed during crack propagation, which can be written as the product of average energy per unit area and fracture cross-sectional area. Therefore, equation (4) can be simplified as follows:

Among them, is the average energy consumption per unit area during crack propagation, reflecting the ability of the material to resist crack propagation.

Substituting equations (12), (23), and (24) to equation (5), the simplified mathematical model of the impact work for rectangular cross-section structure can be obtained

Taking the local length of plastic deformation [11,12], we substitute it to equation (25) and derive

#### 3. Experimental Verification of Impact Work Mathematical Model

The mathematical model described by equation (26) was verified by the impact experiments of 921A and Q235B steel.

##### 3.1. Impact Test Equipment and Test Materials

In order to verify the validity of equation (26), a series of impact tests on rectangular cross-section specimens with different sizes and regression analyses of the impact test results were carried out.

The impact test was conducted at room temperature, and JB-300B pendulum impact tester was selected to carry out impact loading on the specimen. The geometry of the specimens is shown as Figure 3. The span between supports L was 40 mm; the length of all specimens was 55 mm. To control the crack position, a V-notch, whose angle is 45° and depth is 2 mm, was made at the center of each specimen. The radius of curvature at the bottom of the notch is 0.25 mm. Ligament length h is the width of specimen which really takes part in the baring. All specimens were processed by wire cutting and polished to make the roughness meet ASTM A370 standard requirements. Oil pollution was removed with acetone.

Ultra-high strength steel 921A is chosen for impact test in this study; its basic mechanical properties are measured as shown in Table 1.

 Material (A) σs (MPa) σb (MPa) e (%) 921 665 729 25.5

Three impact tests were conducted on each impact specimen of the same size and material, and the average of the three tests was taken as the final result.

##### 3.2. Analysis of Impact Test Results

The impact test results of specimens made of 921A steel are shown in Table 2.

 Specimen number (mm) h (mm) A (J) A1 10 3 30 A2 10 4 46 A3 10 6 102 A4 10 6 108 A5 10 6.5 110 A6 10 7 130 A7 10 8 202 A8 10 8 210 A9 10 9 248 A10 6 8 88 A11 6.5 8 96 A12 8 8 120 A13 9.5 8 166 A14 10 8 202 A15 12 8 228

According to the physical meaning of each parameter in equation (26), boundary conditions of parameters were given in Table 3 for the regression analysis. The intrinsic impact toughness αt and average energy consumption per unit area γave must be positive. The value of αe can be calculated by substituting the yield strength and yield strain of each material in equation (11).

 αt (mJ/mm3) αe (mJ/mm3) γave (mJ/mm2) >0 0.0554 >0

Nonlinear double independent variable regression analyses were conducted with software Mathematica, and the analysis results are shown in Table 4. The R-square value is close to 1, indicating the effectiveness of fitting. However, the value of is close to 0, and the value is close to 1, showing that the value of is extremely insignificant. Compared with the energy dissipating during elastic deformation and the plastic deformation energy, the surface energy of materials used in this research is much smaller than the first two with a couple of orders of magnitude differences [13]. Therefore, as a result of regression analysis, the insignificance of γave is reasonable, suggesting that γave may be ignored under low accuracy requirement.

 Fitting result value R-squared αt 579.765 0.0002 0.989375 αe 0.0554 — γave 0.00034 9.99482

Substituting sizes of the specimens and the fitting results of αe and αt to equations (12) and (23), respectively, we found out that for materials with appreciable plasticity, the energy absorbed by plastic deformation is much larger than that of elastic deformation, which was commonly accepted by existing researches [14].

Substituting the fitting results of 921A steel into equation (26), we obtained the relationship between a and H (Figure 4(a)) and the relationship between a and (Figure 4(b)). Meanwhile, impact test results were also marked in the corresponding coordinates for comparation, as shown in Figure 4.

It can be seen from Figure 4 that the results of impact test were all distributed near the curve base on equation (26). The average difference between the results evaluated with equation (26) and those obtained by impact tests is 10.54%. Considering the considerable discreteness of impact test itself, the accuracy and validity of the mathematical model described by equation (26) can be acceptable. Furthermore, the relation between ligament length h and impact energy A was quadratic approximately and the relationship between thickness W and impact energy a is approximately linear, which is in good agreement with the mathematical model [15].

#### 4. Discussion

The mathematical model of impact energy illuminated the relation between mechanical properties of the material itself, the geometric dimensions of the structure, and the impact work. Thus, it is possible to estimate the impact energy of different size structures with known mechanical properties parameters of the materials. If the impact energy of a certain structure is taken as an index to evaluated impact bearing capability of the structure, this model can be used as a design basis or a tool of reliability evaluation [16].

Furthermore, the commonly used parameters Ak and αk are influenced by the geometry of specimen. Thus, the comparison of Ak or αk between specimens of different or nonstandard sizes cannot reflect their difference in toughness [17]. The intrinsic impact toughness αt in this paper has no relation with the geometry of specimen, which can be regarded as a mechanical property of material itself. Furthermore, αt can be obtained by regression analysis based on the mathematical model constructed in this paper [18]. Considering the convenience of impact test, αt has the potential to be applied to evaluate the toughness of materials such as KIc and KId.

It is worth noting that the mathematical model of impact energy in this paper is established for materials with obvious plasticity. For materials with little plasticity, the mathematical model proposed in this study is not applicable [18]. The mathematical model proposed in this study did not consider the inertia effect of the structure under higher rate impact which needs to be further studied.

#### 5. Conclusion

(1)By analyzing the dynamic response process of the V-shaped gap reduced by the transverse impact of the placed hammer, a mathematical model can be established to estimate the impact energy of different size fragments. This study reveals that the impact energy is linearly related to the width of the notch cross-section. The rules related to the second-time nonlinearity of the thickness have been verified by a series of impact defects.(2)The inherent impact toughness is independent of the size and geometry of the material, which can be called the inherent impact toughness of the material itself. When comparing the volumes of different materials, it is no longer necessary to use specimens of the same standard size for impact test.(3)When the structure is used below the ductile-brittle transition temperature and high-speed impact load, the mathematical model established in this study cannot be used in theory, which needs to be considered in the structural inertia.

#### Data Availability

The labeled dataset used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare no conflicts of interest.

#### Acknowledgments

This work was supported by Youth Program of National Natural Science Foundation of China (51605116).

#### References

1. K. K. Mathur, A. Needleman, and V. Tvergaard, “3D analysis of failure modes in the Charpy impact test,” Modelling and Simulation in Materials Science and Engineering, vol. 2, no. 3A, pp. 617–635, 1994. View at: Publisher Site | Google Scholar
2. Y. Nakamura, M. Yamaguchi, and M. Okubo, “Instrumented charpy impact test of epoxy resin filled with irregular-shaped silica particles,” Polymer Engineering & ence, vol. 33, no. 5, pp. 279–284, 2010. View at: Google Scholar
3. H. Kurishita, H. Kayano, M. Narui, M. Yamazaki, Y. Kano, and I. Shibahara, “Effects of V-Notch dimensions on Charpy impact test results for differently sized miniature specimens of ferritic steel,” Materials Transactions, JIM, vol. 34, no. 11, pp. 1042–1052, 1993. View at: Publisher Site | Google Scholar
4. Y. Takashima, M. Ohata, and F. Minami, “Evaluation method for fracture mechanics-based material toughness from Charpy impact test,” Materials Science Forum, vol. 512, pp. 61–66, 2006. View at: Publisher Site | Google Scholar
5. C. H. Lee, H. S. Shin, K. T. Park, and S. H. Yang, “Evaluation of the applicability of structural steels to cold regions by the Charpy impact test,” Journal of Korean Society of Steel Construction, vol. 23, no. 4, pp. 483–491, 2011. View at: Google Scholar
6. K. Toshiro, Y. Isamu, and N. Mitsuo, “Evaluation of dynamic fracture toughness parameters by instrumented Charpy impact test,” Engineering Fracture Mechanics, vol. 24, no. 5, pp. 773–782, 1986. View at: Publisher Site | Google Scholar
7. S. H. Kim, Y. W. Park, S. S. Kang, and H. D. Chung, “Estimation of fracture toughness transition curves of RPV steels from Charpy impact test data,” Nuclear Engineering and Design, vol. 212, no. 1-3, pp. 49–57, 2002. View at: Publisher Site | Google Scholar
8. R. Chaouadi and A. Fabry, “On the utilization of the instrumented Charpy impact test for characterizing the flow and fracture behavior of reactor pressure vessel steels,” From Charpy To Present Impact Testing, vol. 30, pp. 103–117, 2002. View at: Publisher Site | Google Scholar
9. P. R. Sreenivasan and S. L. Mannan, “Plastic n-factor for three-point bend specimens: analysis of instrumented Charpy impact test results for AISI 308 weld and AISI 316 stainless steels,” International Journal of Fracture, vol. 101, no. 3, pp. 215–228, 2000. View at: Publisher Site | Google Scholar
10. N. S. Cannon and D. S. Gelles, “Charpy impact test results from low activation alloys irradiated to 10 dpa at 365°C,” Journal of Nuclear Materials, vol. 186, no. 1, pp. 68–76, 1991. View at: Publisher Site | Google Scholar
11. H. Kikuchi, M. Harada, K. Ara, Y. Kamada, S. Kobayashi, and S. Takahashi, “Development of apparatus for magnetic measurements of Charpy impact test pieces,” Journal of Materials Processing Technology, vol. 181, no. 1-3, pp. 190–193, 2007. View at: Publisher Site | Google Scholar
12. M. M. Shokrieh, M. A. Torabizadeh, and A. Fereidoon, “Dynamic failure behavior of glass/epoxy composites under low temperature using Charpy impact test method,” Indian Journal of Engineering & Materials ences, vol. 18, no. 3, pp. 211–220, 2011. View at: Google Scholar
13. O. S. Lee and S. K. Hong, “Dynamic fracture characteristics of highly brittle materials by using instrumented Charpy impact test,” KSME International Journal, vol. 11, no. 5, pp. 513–520, 1997. View at: Publisher Site | Google Scholar
14. T. Iwadate and H. Takemata, “Prediction method of fracture toughness KIC transition curve from Charpy V-Notch impact test results,” Journal of the Society of Materials Science, Japan, vol. 41, no. 467, pp. 1241–1247, 1992. View at: Publisher Site | Google Scholar
15. G. V. Kozlov, R. A. Shetov, and A. K. Mikitayev, “Determination of the limit of forced elasticity in the Charpy impact test,” Polymer Science U.S.S.R., vol. 29, no. 9, pp. 2205–2207, 1987. View at: Publisher Site | Google Scholar
16. K. Seo, F. Nogata, and J.-i. Masaki, “The absorbed energy in Charpy impact test for specimens with mechanical heterogeneity,” Journal of the Society of Materials Science, Japan, vol. 35, no. 399, pp. 1405–1410, 1986. View at: Publisher Site | Google Scholar
17. J. A. Wang and N. Rao, “New methodologies for developing radiation embrittlement models and trend curves of the Charpy impact test data,” Journal of ASTM International, vol. 1, no. 9, p. 19, 2004. View at: Publisher Site | Google Scholar
18. T. Otsuka, K. Hashizume, and M. Sugisaki, “Charpy impact test of oxidized and hydrogenated zircaloy using a thin strip specimen,” Journal of Nuclear Science and Technology, vol. 41, no. 3, pp. 247–251, 2004. View at: Publisher Site | Google Scholar

Copyright © 2021 Wei Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.