Study on the Fracture Process Zone near the Mode I Dynamic Crack Tip

Evaluation of the shape and size of the fracture process zone near the mode I dynamic crack tip is still a problem unsolved completely at present. The research on the relationship between the fracture process zone and crack velocity near the mode I dynamic crack tip is quite limited, and some researchers have also developed experimental methods or numerical methods. In this research, based on the theory of elastodynamics and the complex stress function method, an approximate method for solving the mode I dynamic crack problem was proposed. The fracture process zone near the mode I dynamic crack tip was analyzed. The results showed that the areas of the fracture process zone determined based on the approximate method are nearly the same as the results obtained based on the well-known stress fields. The approximate method could provide a good reference for determining the fracture process zone near the mode I dynamic crack tip since no analytic methods had been found for evaluating the fracture process zone near the dynamic crack tip to the authors’ knowledge.


Introduction
e study on the fracture process zone has been mainly focused on the fracture process zone near the static crack tip, while the analysis of the fracture process zone near the dynamic crack tip is very complicated due to the inertia and the dynamic effect of materials. For the fracture process zone near the static crack tip, some analytical models, such as the cohesion model [1], the virtual crack model [2], and twoparameter model [3], treat the fracture process zone as linear in form. e cell model proposed by Broberg [4][5][6][7] well explained the existence of the fracture process zone at the crack tip and provided a very visual explanation for future researchers. At present, the research studies on the dynamic crack tip fracture process zone are very limited, and some researchers have also developed experimental methods or numerical methods. For the determination of the fracture process zone near the dynamic crack tip, although Broberg [8] and Fineberg [9] concluded that the fracture process zone at the dynamic crack tip increases with the increase of the crack velocity, the authors have searched relevant literature and found no relevant analytical methods yet. e determination of the profile of the fracture process zone at the dynamic crack tip is still an unsolved problem. e integral transformations with dual integral equation method and Wiener-Hopf method are difficult to allow for obtaining analytical solutions due to the complexity of calculations. e complex function method is much simpler than the integral transformation methods because it does not need to carry out the complex inversion operations.
In this study, based on the theory of elastodynamics and the complex stress function method, an approximate method for solving the dynamic crack problem was proposed. e fracture process zones at the mode I dynamic crack tip were analyzed. e contours and areas of the fracture process zones were determined by using the maximum tensile stress criterion for different crack velocities. Finally, the results based on the approximate method proposed in this paper were compared with the results calculated based on the well-known stress field at the crack tip to verify the feasibility and correctness of the approximate method.

Method
e crack tip is assumed to move along the direction of xaxis at a constant speed V. e far-field tension σ is applied, see Figure 1. e fixed coordinate systems (x 1 , y, t) and the moving coordinate system (x, y, t) have the following relationships: e coordinate transformation gives For plane elastic problem, the equation of motion takes the following form: where ρ is mass density and u and v are displacements along x and y axes. e relationships between stresses and strains give [10] where μ and λ are Lame's constants and θ is the first invariant of strain tensor, that is, where Since the analytical solution of equation (16) cannot be obtained accurately, and the expression of the stress function of the dynamic crack cannot be determined, the stress field near the mode I crack tip cannot be determined directly by solving equation (16). e controlling equation of the static problem is one biharmonic equation, and the stress function method for solving static fracture process zone near the crack tip is already very complicated. Equation (16) is very similar in form to the biharmonic equation. Hence, we consider using the method for solving static crack problems to solve the dynamic problem of the stress field near the crack tip. Let So, equation (16) can be rearranged as follows: So, equation (19) can be rearranged as follows: It can be seen that equation (21) is a fourth-order partial differential equation about x and y. It is very difficult to solve equation (21) without some efforts. However, the second term on the left side of the equation contains parameter β, and the values of β at different crack velocity are shown in Table 1.
It can be seen from Table 1 that β increases with the crack velocity.
e maximum value of β is reached when the Rayleigh wave velocity is approached. However, the generally observed maximum crack propagation velocity in experiments is 0.34-0.45 times of the Rayleigh wave velocity, and the value of β is in the range of −0.00025 to −0.00139. e value of β is very small, and the second term of equation (21) can be neglected. Hence, the following equation can be used to explore the approximate solution of equation (16): At this time, equation (15) is the biharmonic equation in the coordinate system (xy) and has the same form as that for the static crack. en, we can find the approximate solution for the dynamic crack by using the stress function method.
For mode I crack, Westergaard [11] assumed the stress function taking the following form: e complex stress function is [12,13] Mathematical Problems in Engineering en, the stress fields near the mode I crack tip can be calculated: e size of the FPZ can be calculated theoretically based on the appropriate criterion and the crack tip stress field. e maximum tensile stress criterion for brittle materials can be expressed as follows: where σ u is the tensile strength and is assumed to be independent of the crack velocity. When more than one term in the expansion of the stress function is taken, the expressions for the stress field based on the stress function method at the crack tip are shown in Appendix.

Results
By substituting equation (25) into equation (26), the implicit equation of the FPZ can be easily obtained. As an example, we take the following material parameters: the longitudinal wave velocity c 1 � 5370 m/s, the shear wave velocity c 2 � 3180 m/s, and the Rayleigh wave velocity c � 2900 m/s. e contours of the FPZ constructed by taking different terms of the stress function of equation (24) for different crack velocities are shown in Figure 2.
From Figure 2, it can be seen that the fracture process zone increases with the increasing crack velocity. e higher the velocity is, the more significant the change of the fracture process zone is. e area of the fracture process zone at 0.38c is 1.026 times that at 0.19c, and the area of the fracture process zone at 0.76c is 1.34 times that at 0.95c, or 1.657 times that at 0.19c. Taking the velocity of 0.38c as an example, the calculated areas of the fracture process zone are shown in Table 2. e area of the fracture process zone calculated by the first two terms of the stress function is the largest. e more the stress function terms are taken, the smaller the change rate of the area of the fracture process zone at the crack tip is, which also indicates that the calculation results become more and more stable.

Discussion
e approximate method is proposed in the second part by ignoring the influence of the less-significant parameters. e well-known stress field was used to calculate the fracture process zone at the crack tip, and the results by using the approximate method were compared to the results by using the well-known stress fields to verify the correctness of the approximate method. e expressions for the mode I stress fields are given as follows [14]: where If equation (27) is substituted into equation (26), we can obtain the equation for the fracture process zone at the crack tip r � f (θ). e higher-order terms of the stress field have little effect on the stress field near the crack tip of mode I, but the second term has a significant effect on the fracture process zone at the crack tip. erefore, the first term and the first two terms in the stress field expansion expression are considered, respectively, to calculate the fracture process zone at the crack tip of mode I. K 2 0 is a parameter similar to the stress intensity factor, which is difficult to determine. Because the dependence of T-stress on the velocity for dynamic crack is not clear, the ratio of T-stress to the far-field stress is 0.526 referred to the work by Seed and Nowell [15].
When the first term and the first two terms of the stress function expansion are taken respectively, the comparisons of the areas of the fracture process zone at the mode I crack tip calculated based on the approximate method and the well-known stress field expression are shown in Table 3. It can be seen that the results of the mode I crack tip fracture process zone calculated by the approximate method are highly consistent with those of the crack tip fracture process zone calculated based on the known stress field, and both increase with the increase of the crack velocity. When the velocity is relatively high, the area of the fracture process zone at the mode I crack tip calculated by the approximate method is slightly smaller than that calculated based on the well-known stress field, which is believed to be caused by neglecting the influence of small parameters β. e comparison of the calculation results of the two methods in Table 3 verifies the correctness of the approximate method in this paper for evaluating the fracture process zone at the tip of mode I dynamic crack, which can provide a reference for the study of the mode I dynamic crack problem.

Conclusions
An approximate method was proposed by transforming the governing equation and disregarding the less-significant parameter. e areas and sizes of the fracture process zone of the mode I dynamic crack were evaluated by the approximate method. e correctness of the method was verified by comparing the results calculated by the approximate method and the well-known stress fields. e conclusions can be drawn as follows: (1) e fracture process zone near the mode I dynamic crack tip distributed symmetrically about the crack surface. e higher the crack propagation velocity, the larger the size of the FPZ for the mode I crack.
When the velocity approached the Rayleigh wave velocity, the fracture process zone changed more significantly.
(2) e fracture process zone can be calculated by using this approximate method with different terms of the stress function. e sizes and areas of the fracture process zone calculated by the proposed method were very close to those calculated based on the wellknown stress field. erefore, the proposed method is suitable for the mode I dynamic crack. (3) Although the well-known stress field is widely accepted, it is just an asymptotic solution for the mode I dynamic crack tip, i.e., n � 1. When n > 1, the parameter K 2 0 cannot be determined except by experiment. is means that it is difficult to use the well-known stress field expression to evaluate the FPZ for the dynamic crack tip for n > 1. However, the approximate method can be used to obtain the results with different accuracies by taking different terms from the expansion of the stress function. It is also noteworthy that the method is very easy to implement for determining the shape and size of the FPZ of the mode I dynamic crack.

Appendix
In this appendix, the authors will show the expression of the stress field used in the calculations.
So, the stress fields for taking the first term of the expansion of the stress function are shown as follows: When the first two terms of the expansions of the stress functions are taken, the stress fields near mode are e expressions of the stress fields for the mode crack when taking first three terms are shown as follows: