Abstract

In this section, the dynamic propagation behavior of a penny-shaped interface crack in piezoelectric bimaterials is analyzed. The objective of this paper is to use the boundary conditions of the penny-shaped interface crack to study the dynamic propagation of the crack under the action of load, so as to provide some valuable implications for the fracture mechanics of the piezoelectric bimaterials and simulate the interface crack between piezoelectric bimaterials, it is necessary to establish a suitable model and give appropriate boundary conditions according to the actual situation. The elastic displacement and potential equations are constructed according to the structural characteristics of the circular crack. In the case of a given displacement or stress, the Laplace transform and Hankel transform are used to simplify the problem into an integral equation with unknown functions. According to the boundary conditions, the corresponding unknowns are obtained, and the closed solution is derived. The results show that the fracture toughness of a penny-shaped interface crack in piezoelectric bimaterials is related to the thickness of the material, the impact time, the material characteristics, and the electric field. At the same time, it can be found that different materials have different roles in the crack propagation, so it is very important to study the crack opening displacement (COD) intensity factor of the crack for safety design.

1. Introduction

1.1. Background

Piezoelectric materials are widely used to make all kinds of transducers and sensors because of their piezoelectric effects. They are commonly used in electronic, excited light, ultrasonic, hydroacoustic, microacoustic, red external, navigation, biological, and other technical fields. However, the brittleness of piezoelectric ceramic material in mechanical properties will lead to the generation and expansion of cracks from the stress concentrations, leading to the failure of components. It is necessary to analyze the crack propagation process of piezoelectric materials to improve the operational performance of high-voltage electric ceramic elements and predict their active life.

In order to meet the practical needs, experts have used experimental, theoretical, and simulation methods to study the influence factors of piezoelectric materials in the different situations [1, 2]. The forced vibration of a piezoelectric plate with initial stress and the factors affecting the dynamic stability of a prestressed piezoelectric plate are studied [3, 4].

Based on the above situation, research on the fracture mechanics of piezoelectric ceramics has received attention from the experts [513]. Influence of incomplete bonding on the dynamic response of prestressed sandwich plate–strip with elastic layers and a piezoelectric core was studied by Daşdemir [14]. The influence of poling direction and imperfection defects of the prestressed system with a piezoelectric core bonded to elastic faces was analyzed [15]. In most cases, the crack is simulated by a penny shape [16, 17]. A 3D problem of a half-space with penny-shaped cracks has been studied [18]. Zikun [19] solved a penny-shaped crack problem and obtained analytical expressions of the stress field and electric displacement field near the crack tip. Using the Dugdale hypothesis and Hankel transform theory, Danyluk et al. [20] analyzed the penny-shaped crack in a thick transversely isotropic elastic layer. Kogan et al. [21] obtained the stress intensity factors of a penny-shaped crack. The plasticity of a penny-shaped Dugdale crack tip in an infinite elastic medium has also been estimated [22].

A penny-shaped crack has been considered in a piezoelectric medium to analyze the coupling behavior [23]. Penny-shaped cracks in 3D piezoelectric media are analyzed by the extended displacement discontinuity method [24]. The Dugdale plastic zone of a penny-shaped crack under axisymmetric loading has been analyzed [25]. Wu at al. [26] studied a penny-shaped crack in a piezoelectric layer sandwiched between two elastic layers with the electrical saturation and mechanical yielding zones. The crack propagation behavior of an infinite 1D hexagonal piezoelectric quasicrystal plate with a penny-shaped dielectric crack has been analyzed [27].

The crack propagation will change into a dynamic behavior with time under an explosion or impact load. Wang et al. [28] investigated the penny-shaped interface crack configuration in orthotropic multilayers under dynamic torsional loading by utilizing Laplace transform and the Hankel transforms technique. The dynamic behavior of a penny-shaped crack in a magnetoelectroelastic materials has been analyzed in detail [29]. A point force method was proposed for obtaining the transient response of dynamic penny-shaped cracks in multilayer sandwich composites [30]. The dynamic behavior of a magnetoelectroelastic material with a penny-shaped dielectric crack under impact loading has been determined [31]. A mathematical expression for the dynamic behavior (open function circular) of a permeable penny-shaped crack in an infinitely porous elastic solid has been presented [32].

In order to overcome the weak characteristics of the piezoelectric effect of a single piezoelectric material and better expand the piezoelectric effect, the piezoelectric double material has emerged, which is better applied in intelligent structures such as transducers, sensors, and drivers, and it has also produced great application and economic value. At the same time, because of its superiority, it has also been studied and analyzed by many experts and scholars [3339]. Under loading conditions, the interface of piezoelectric bimaterials is prone to fracture due to surface offset. In addition, crack propagation becomes a dynamic behavior when piezoelectric bimaterials are impacted in practice [40]. There are few studies on the dynamic change of a penny-shaped interface crack in piezoelectric bimaterials. Therefore, the dynamic analysis of a penny-shaped interface crack in piezoelectric bimaterials is very important in designing practical engineering applications.

It is well-known that defects such as interface cracks, holes, and dislocations seriously affect the mechanical behavior and change their strength. In addition, the interface crack boundary conditions are of great significance for selecting fracture criteria and predicting crack propagation. In addition, the penny-shaped interface crack, which is the research target of this paper, is a common type of crack, and the dynamic propagation to be studied is also the most practical and research value. However, the dynamic propagation of a penny-shaped interface crack is difficult in experiments and simulations. Therefore, in order to meet the needs of development and safety, the objective of this paper is to use the boundary conditions of the penny-shaped interface crack to study the dynamic propagation of the crack under the action of load, so as to provide some valuable implications for the fracture mechanics of the piezoelectric bimaterials.

1.2. Outline

After the introduction, the model and constitutive equations are described. Section 3 discusses the calculation methods under the Laplace transform and Hankel transform. In Section 4, the dynamic field intensity factors of crack-tip propagation are derived. In Section 5, numerical results are obtained based on the model and calculation. The final section presents some conclusions drawn from this study.

2. Problem Statement and Formulation

In Figure 1, a penny-shaped interface crack made in piezoelectric materials one and two is considered. The thickness of piezoelectric Materials 1 and 2 are and , respectively. The radius of the penny-shaped interface crack is . Let us analyze the case of the poling direction along the -axis in polar the coordinates . The center of the crack is the origin of coordinates, and its region is . The names of the parameters can be seen in nomenclature.

According to Tiersten [41], thess balance equation (divergence equation) of the dynamic behavior of a penny-shaped interface crack in piezoelectric bimaterials can be expressed as follows:

The constitutive equations (equations of state) [41] can be written as follows:

3. Solution to the Problem

Inserting Equation (2) into Equation (1), we have

In order to further solve the problem, the piezoelectric bimaterials are in a mechanical and electrical static state when . The initial conditions for piezoelectric bimaterials are zero at , which can be expressed by the formulas as follows:

For piezoelectric media, the boundary conditions of the crack surface are well known. In this paper, we consider the boundary conditions in the presence of a circular coin-type interface crack. For piezoelectric bimaterials, we assume that there are impact stresses and electric fields above and below, and it is sufficient to solve the corresponding problems presented in the following boundary conditions according to symmetry. The boundary conditions under impacts of elastic stress and electric field at are as follows:

In this paper, the defect problem of piezoelectric bimaterials with finite thickness is studied. According to symmetry, the upper and lower materials are continuous in the plane (crack plane). On the other hand, applying symmetry to the plane at , we get

To solve, the following is the Laplace transform concerning the time where denotes the Bromwich path of integration.

Using the Laplace transform, Equations (7)–(15) can be transformed into

From Equations (17)–(19), the elastic displacement and electric potential by unknown functions and ( and ) in the transformed domain can be obtained as follows:where , , , and are zeroth order and first order Bessel functions of the first kind, respectively, known constants and are determined by the Hankel and Laplace transform.

By substituting Equations (26)–(28) to Equations (3)–(5), we find satisfy the followingwhereand meet this condition

From Equation (2) and Equations (26)–(28), the solutions can be obtained as follows:where

According to Li and Lee [42], the relationship between and is as follows:

According to the boundary conditions, we have

Inserting boundary conditions Equations (23) and (24) into Equations (27) and (28), respectively, they are

Introducing new functions and , we have

From Equations (42)–(44), one can see that the unknown functions can be obtained as follows:

From the study of Abramowitz and Stegun [43], we have

By applying Equation (46) and substituting Equations (27) and (28) into Equations (23) and (24), respectively, the values of elastic displacement and electric potential on the can be obtained as follows:

It is moreover, inserting Equations (32) and (34) into Equations (21) and (22), respectively, we arrive at

Then we have

Using what we already know [42]

We finally obtain

By substituting Equations (43) and (44) into Equations (54) and (55), respectively, one can obtain the following results

Using the known relation [42]where is the Dirac delta function.

By applying Equation (58), Equations (56) and (57) can be further read as follows:where

Let us consider the following limiting form

Equations (59) and (60) can be further written as follows:where

For the convenience of calculation, the following transformations are introduced as follows:

Hence, Equations (65) and (66) can be expressed as follows:

4. Dynamic Field Intensity Factors

From the perspective of fracture mechanics, the dynamic field intensity factor is an essential factor in characterize the penny-shaped crack. From Equation (32), we have

After some calculation, the expressions for and can be obtained as follows:

We note that

The dynamic field intensity factor in the Laplace transform domain can be defined as follows:

Similarly, the field intensity factor is related to the propagation displacement of the crack tip and the electric potential, so it can be defined as follows:

According to Equations (69)–(71), one can obtain

5. Numerical Examples

When the dynamic field intensity factors exceed the corresponding critical value, the crack expands. Some numerical calculations are shown to describe the dynamic behavior of a penny-shaped crack. The material coefficients of the numerical examples are shown in Table 1.

Let us take the inverse Laplace transform as Stehfest [44] proposed. In our example, the stress impact load is selected to be 4.2 MPa. The parameters relating the electrical and mechanical loadings are expressed by and . Let us define the dimensional function as , where is the static COD intensity factor ( is in Equation (81)).

When the double material degenerates into a single material, we can calculate the corresponding value (Appendix B). Figure 2 is drawn to describe the trend of relative to (). The two piezoelectric single materials are PZT-6B and BaTiO3, respectively. The of both materials decreases monotonically and tends to change smoothly with increasing . On the other hand the changes in the two curves are shown in Figure 2. We can see that the change in the value of the material BaTiO3 is greater than that of material PZT-6B, which in a particular sense implies the correctness of our results.

For the plot with change in Figure 3, a thinner piezoelectric layer can be found to cause an increase in , which means that the contribution of the electric field to the fracture toughness is evident for the thin piezoelectric layer in Figure 3. When is also insensitive.

Figure 4(a)4(c) discuss the influence of the applied electric fields and on . Under the impact problem, this time normalization is denoted by , where is the shear wave velocity. In Figure 4, for a given pressure increases with increasing of electric fields and . At a fixed time, it can be found that an increase in will lead to an increase in , which is consistent with the result shown in Figure 3. For a fixed electric field , an increase in leads to an increase in . However, we can see that the increase in flattens out in the later period. From Figures 24 show that different materials play different roles in the crack propagation.

Figure 5 shows that when the substrate is constant, the thicker the material 1, the less influence it has. In other words, when the thickness of Material 2 is constant, the thickness of Material 1 is small, and the COD intensity factor changes greatly, indicating that the thinner of Material 1 when appropriate, the more conducive it is to safe design.

It can be seen from Figures 3 and 4 that both positive and negative electric fields will affect the propagation of cracks. At the same time, it can be found that the COD intensity factor value is large under a positive electric field, while the COD intensity factor value is small under a negative electric field, which indicates that the positive electric field will cause the expansion of the piezoelectric material, while the negative electric field will cause the contraction of the piezoelectric material. This conclusion is consistent with that of the study by Li and Lee [42]. In other words, crack propagation can be controlled by the regulation of electric field.

For the plot with in Figure 6, a thinner piezoelectric layer can be found to cause an increase in . In addition, observations of Figures 46 suggest that the electromechanical coupling coefficient plays an important role in crack propagation: the stronger the electromechanical coupling effect, the smaller is.

All the c figures in Figures 46 studied the crack growth of different piezoelectric bimaterials. By comparing the images, we can find that when Materials 1 and 2 are different, the variation trend of COD intensity of crack growth is basically the same as that when Materials 1 and 2 are the same. That is, different material parameters will affect the change of COD intensity factor, so in practical applications, the appropriate material should be selected according to the actual situation.

It can be seen from Figures 5 and 6 that when the thickness of Material 2 is constant, the thickness of Material 1 is small, and the COD intensity factor changes greatly, indicating that the thinner the material when appropriate, the more conducive it is to safe design.

6. Conclusion

Based on the piezoelectric theory, the dynamic penny-shaped interface crack propagation of piezoelectric bimaterials is analyzed. The boundary conditions are transformed into a nonlinear Fredholm integral equation by using the Hankel transformation technique. Numerical solutions are given. According to the given surface displacement and stress of the layer, the corresponding models are constructed, and the displacement functions that meet the conditions are established. Solving the model, the influences of the electric field, impact time and layer thickness on the dynamic COD are analyzed accordingly. The results show that the COD increases with decreasing . At the same time, when the speed and crack size are fixed, the COD increases with increasing impact time, and tends to be flat after reaching the peak value. The stronger the electromechanical coupling effect is, the smaller is. At the same time, different materials have different roles in crack propagation, so it is very important to study the COD factor of cracks for safety design.

Appendix

A. Elements of Matrix [C]

B. The Element Value of Matrix [C] in Single Material

The value of corresponds to the value of for the matrix in Appendix A.

Nomenclature

:Elastic constants
:Dielectric constants
:Radial and axial components of displacements vector
:Electric displacement
:Shear wave velocity
:Materials 1 and 2
:Electric loading
:Piezoelectric constants
:Density
:Stress tensor
:Step function
:Impact time
:Mechanical loading
:Dielectric load normalization.

Data Availability

All the numerical calculated data used to support the findings of this study can be obtained by calculating the equations in the paper, and piezoelectric material parameters are taken from the study of Li [40] and Li and Lee [42].

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (11972019 and 52275567), Shanxi Provincial Key Research and Development Project (202102090301027), and Shanxi Postgraduate Innovation Project (2022Y668).