Research Article  Open Access
Wei Wang, Ping Wang, Xuesong Liu, Zhibo Dong, Hongyuan Fang, "Mathematical Model for Charpy Impact Energy of VNotch 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 VNotch Specimens
Abstract
Firstly, by analyzing the response of Charpy Vnotch 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 A_{k}, a value measured by the Charpy impact test, is mainly used to control the quality of metallurgical and thermal processing products and determine the ductilebrittle transition temperature of materials [1]. However, the determination of the A_{k} 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 crosssectional 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 fullsize specimens. Sokolove et al. [3] studied the correlation between the ductilebrittle transition temperatures measured with impact specimens of different sizes, developed a procedure for transforming data from subsize and the fullsize 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 smallsize specimen and the fullsize specimen. Schill et al. [4] studied the correlation between the impact energy of a certain ferritic steel smallsize specimen and a fullsize 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 smallsize impact specimens and fullsize standard impact specimens based on the test results. Yang et al. [6] used GANN 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 heattreated 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 crosssection 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 crosssectional shapes. He further verified that the impact energy is proportional to the section modulus in bending of rectangular crosssection structures and found that the regularity is not significant for nonrectangular crosssections. 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 loadbearing form of the Vnotch 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 loadbearing 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 VNotched Specimen Subjected to Lateral Impact
The deformation mode of Vnotch specimen subjected to lateral impact, as shown in Figure 1, can be regarded as a dynamic 3point 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 stressstrain 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 Vnotch 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 VNotch Impact Energy
Based on the above analysis, it can be concluded that the impact energy consists of three parts: the elastic deformation energy E_{e}, the plastic deformation energy of the structure E_{p,} and the energy consumed by crack propagation in fracture section E_{c}. 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 realtime 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 realtime bending moment of the specimen during bending deformation; is the energy consumed on the area units at different locations of the crosssection 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 crosssection structure, the bending moment of the crosssection 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 crosssection of the notch center; is the thickness of the crosssection of the notch center.
Substituting equation (6) to equation (2) gives
In the elastic stage, the maximum stress and strain on the notch crosssection 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 crosssection 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 realtime bending moment of crosssection in plastic deformation stage can be written as a product of correction coefficient k_{r}(θ) and the ultimate plastic moment of crosssection M_{p}, that is,
The ultimate plastic bending moment M_{p} 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. M_{p} does not change in the deformation process. The actual change of bending moment during deformation is considered by multiplying the correction coefficient function k_{r}(θ).
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 crosssection 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 crosssection 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 l_{h} shall change during the deformation. However, for mathematical simplicity, some fixed values of l are recommended. In the case of rectangular crosssection beam, l_{h} 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 crosssectional 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 crosssection 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 crosssection specimens with different sizes and regression analyses of the impact test results were carried out.
The impact test was conducted at room temperature, and JB300B 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 Vnotch, 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.
Ultrahigh strength steel 921A is chosen for impact test in this study; its basic mechanical properties are measured as shown in Table 1.

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.

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).

Nonlinear double independent variable regression analyses were conducted with software Mathematica, and the analysis results are shown in Table 4. The Rsquare 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.

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.
(a)
(b)
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 A_{k} and α_{k} are influenced by the geometry of specimen. Thus, the comparison of A_{k} 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 K_{Ic} and K_{Id}.
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 Vshaped 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 crosssection. The rules related to the secondtime 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 ductilebrittle transition temperature and highspeed 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
 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
 Y. Nakamura, M. Yamaguchi, and M. Okubo, “Instrumented charpy impact test of epoxy resin filled with irregularshaped silica particles,” Polymer Engineering & ence, vol. 33, no. 5, pp. 279–284, 2010. View at: Google Scholar
 H. Kurishita, H. Kayano, M. Narui, M. Yamazaki, Y. Kano, and I. Shibahara, “Effects of VNotch 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
 Y. Takashima, M. Ohata, and F. Minami, “Evaluation method for fracture mechanicsbased material toughness from Charpy impact test,” Materials Science Forum, vol. 512, pp. 61–66, 2006. View at: Publisher Site  Google Scholar
 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
 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
 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. 13, pp. 49–57, 2002. View at: Publisher Site  Google Scholar
 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
 P. R. Sreenivasan and S. L. Mannan, “Plastic nfactor for threepoint 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
 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
 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. 13, pp. 190–193, 2007. View at: Publisher Site  Google Scholar
 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
 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
 T. Iwadate and H. Takemata, “Prediction method of fracture toughness KIC transition curve from Charpy VNotch 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
 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
 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
 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
 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
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.