Abstract

Coating often plays a role in monitoring and protecting substrates in engineering applications. Interface cracks between the coating and the substrate can lead to crack growth under the action of external loading and will cause device failure. In this paper, the behavior of a fine-grained piezoelectric coating/substrate with a Griffith interface crack under steady-state thermal loading is studied. The temperature field, displacement field, and electric field of the coupling of thermal and electromechanical problems are constructed via integral transformation and the principle of superposition. Thus, problems are transformed into a system of singular integral equations, and the expressions of thermal intensity factor, thermal stress intensity factor, and electric displacement intensity factor are obtained. We used a numerical calculation and a system of singular equations to obtain the relationship of strength factor with material parameters, coating thickness, and crack size.

1. Introduction

Piezoelectric materials have played an important role in the production of intelligent structures. Due to strength demand, we often bond another piezoelectric material on the surface of the device to protect it or monitor the device. However, piezoelectric composites are prone to cracks at the interface during the process of fabricating defects and load conditions. This can even lead to structural failure. Therefore, problems of interface fracture of the piezoelectric composite under all kinds of loads have attracted wide concern [19]. Qin and Mai [10] analyzed the interface crack problem of a piezoelectric bimaterial subjected to combined thermal, mechanical, and electrical loads via Stroh’s formalism and the singular integral equation method. However, the influence of material size on the thermal intensity factor has not been considered in the numerical analysis. Ueda [11, 12] investigated the problem of a parallel crack in a piezoelectric strip under thermoelectric loading and the FGM material strip containing an embedded crack or an edge crack perpendicular to its boundaries. Wang and Noda [13] studied the problems of the piezoelectric material strip with a Griffith interface crack under thermal loading and discussed the effects of polarization direction, crack size, and location on the thermal strength factor. Kang et al. [14] considered the influences of piezoelectric laminated composite shells with interlaminar stresses under electrical, thermal, and mechanical loads. They concluded that the interlaminar shear stress of the laminated shell can be reduced by selecting the appropriate electric field value. Cook and Vel [15] examined a comprehensive multiscale analysis of laminated plates with integrated piezoelectric fiber composite actuators. They used the coupling method of microscale field variables and macroscale field variables to solve the three-dimensional macroscopic equilibrium equation of the piezoelectric laminated plate under arbitrary boundary conditions solved. Li et al. [16] presented a semianalytical technique to solve fracture behaviors of piezoelectric composites under thermal loading. This method could analytically represent the resulting stress and electric displacement distribution along the radial direction obtained. Kuuttia and Virkkunen [17] justified the surface crack behavior under periodic thermal loads. The finite element method and weight function method could simulate the crack behavior over the entire load cycle. Mojahedin et al. [18] considered the mechanical behavior of a solid circular plate which is made of saturated and unsaturated porous material with piezoelectric actuators under thermal loading. They discussed the effect of porosity on the thermal and mechanical stability of the material. Herrmann and Loboda [19] performed fracture analysis on an electrically impermeable interface crack with contact zones in thermopiezoelectrical bimaterials. Using the admissible directions of the heat and the electrical fluxes, the dependencies of the electrical intensity factors on the intensities of the thermal and electrical fluxes were discussed.

Many researchers have studied the interface fracture problems of piezoelectric materials under loading [2024]. However, all of polycrystalline materials studied above are composed of multidomain large grains, and such sizes cannot meet the requirements of the current intelligent devices, but fine-grained piezoelectric materials are gaining increasing interest. However, there are few studies on the mechanical properties of fine-grained piezoelectric materials.

As an important fracture parameter that can measure structural safety, the intensity factor plays a crucial role in checking the safety of structures. In this article, a Griffith interface crack in fine-grained piezoelectric coating/substrate under steady-state thermal loading is established. An integral transformation method and superposition principle of the solution of the equation are used to transform the thermal load problem into a singular integral equation to obtain the intensity factor conveniently. Furthermore, the interaction between the intensity factor and material parameters is investigated. The results indicate that the larger elastic modulus and thinner coating thickness improve the safety of the coating/substrate structure.

2. Problem Formulation

Figure 1 shows that a fine-grained ceramic powder coating and the substrate are bond together by plasma spraying, while the crack with a length of is along the interface between the coating and the substrate. A fine-grained piezoelectric coating/substrate structure was obtained by the polarization treatment of the coating. The coating and substrate are polarized along the -axis, and both are transversely isotropic. The thicknesses of the coating and substrate are and , respectively. We assumed that the crack faces remain thermally and electrically insulated, and and are the environment temperatures.

The constitutive equations for the elastic field are [12]

In Equations (1), (2), (3), (4) and (5), the and are the displacement components; is the electric potential; is the temperature change; , , and are stress components; and are electric displacement components;, , and are elastic modulus; , , and , are piezoelectric and dielectric constants, respectively. and are the stress-temperature coefficients; is the pyroelectric constant. The superscript () stands for the fine-grained piezoelectric coating and piezoelectric substrate, respectively.

Assume that the temperature satisfies the Fourier heat conduction equation, as follows:where , , and are coefficients of thermal conductivity. .

The equations of equilibrium are [12]

Without loss of generality, the boundary conditions can be written as

The thermal boundary conditions are specified aswhere and are environment temperatures, is the thermal flow, and and are coefficients of thermal conductivity.

Since the temperature field does not depend on the stress and displacement field, it can therefore be calculated separately, and then, the thermal stress and potential displacement field can be solved.

3. Solution to the Problem

3.1. Temperature Field

By using the Fourier integral transform, the solution of temperature field can be expressed aswhere , , and are the unknown functions .

To solve the problem easily, the density function is defined as follows:

According to Equation (11a), the dislocation density function should make the following single-valued conditions valid:

Substituting formula (12) into Equations (11a), (11b), (11c) and (11d), we obtain

Equations (16), (17) and (18) can be written in the form of a matrix as follows:where

Solving the linear Equation (19), we get

The expressions of and , , are given in Appendix A. Obviously, the unknown functions , , , and depend on . Once is determined, the temperature field could be obtained.

By substituting Equation (12) into Equation (11c), we obtainwhere

Furthermore, Equation (23) can be written aswhere

Through changing the order of the integration on the cracks’ surface, we obtainwhere

Equation (26) is a singular integral equation with Cauchy kernel. The density function in Equation (26) has the square root-type singularity. Thus, we can express as follows:where is a continuous function defined in the interval .

The thermal flow intensity factors are defined by [13].

According to Equation (26), the singular part of is

Substituting Equation (28) into Equations (31) and (32), we obtain

3.2. Thermal Stress and Electric Displacement Fields

Similar to the solution of linear equation, the solutions to Equations (7), (8) and (9) of equilibrium for displacements and electric potential consist of a homogeneous solution and a particular solution.

Equations (7), (8) and (9) can be written into homogeneous forms as follows:

By using the Fourier integral transform, the solutions of Equations (34), (35) and (36) can be expressed as [13]where , , , and , , are unknown functions to be determined.

The special solutions of Equations (7), (8) and (9) can be expressed aswhere , , , , , and satisfy the expression as follows:where

Therefore, the solution of Equations (7), (8) and (9) can be written as

Substituting Equations (45), (46) and (47) into Equations (2), (3), and (5), the thermal stress and electric displacement fields are obtained:where

In order to solve the thermal stress and electric displacement fields, the dislocation density functions are introduced:

Substituting Equations (52), (53) and (54) into Equations (10a), (10b), (10c), (10d), (10e), (10f), (10g), (10h), (10i), (10j), (10k) and (10l), we obtainwhere

The expressions of , , and , , , are given in Appendix B. Similar to the process of solving the temperature field, once , , and are determined, the thermal stress and electric displacement fields could be obtained.

Substituting Equations (55), (56), (57), (58), (59) and (60) into Equations (48), (49) and (50), we get

Equations (64), (65) and (66) are divided into two parts, as follows:where

Obviously, it is easy to solve the solution of Equations (71), (72) and (73) according to Equations (39), (40), (41) and (42). Next, we need to transform Equations (68), (69) and (70) into singular integral equations of the first kind with Cauchy kernel.

In another form, integral Equations (68), (69) and (70) will be

The expressions of and , , are given in Appendix C.

Substituting Equations (61), (62) and (63) into Equations (74), (75) and (76), we getwhere

On the cracks’ surface,where

The expressions of and , , are given in Appendix C. Ultimately, Equations (68), (69) and (70) are reduced to coupling the first kind Cauchy singular integral equations.

In Equations (61), (62) and (63), the unknown function , with the following conditions:

Considering that the solution of the singular integral equation will use Chebychev polynomials which are defined in the interval [−1, 1], we introduce the following normalized quantities:

Using Equations (84) and (85), integral Equations (79), (80) and (81) may be expressed as

The crack-tip behavior can be characterized by square-root singularity weighed down by the oscillation. However, because the zone of oscillation is very small [25], the unknown function , , can be expressed as follows:where , , are continuous functions defined in the interval .

Using the Lobatto-Chebyshev collocation method, one can transform Equations (79), (80), (81) and (83) into a system of algebraic equations:where

is the node number of quadrature formula; and are the zero points of the first and second kinds of Chebyshev polynomials. According to Theocaris and Ioakimidis [26] as well as Lu [27], the convergence rate of the Lobatto-Chebyshev collocation method is first order (i.e., , where stands for step size).

The intensity factors are defined by [13]

Substituting Equation (87) into Equations (90), (91), (92), (93), (94) and (95), we obtain

According to the numerical process adopted in [28], the thermal stress and electric displacement fields can be calculated finally.

4. Numerical Process and Discussion

In what follows, we assume that the piezoelectric substrate is cadmium selenide piezoelectric ceramic. Its parameters are given as [24]where superscript 2 represents the substrate material. Assume that the coefficients of thermal conductivity are , .

4.1. Verification

We chose a special case to verify the foregoing derivation. We first assumed that the fine-grained piezoelectric coating and substrate are the same material. The problem in Figure 1 then becomes that of a Griffith crack in a single material plate with thermal loading.

Therefore, in Equation (26), becomes

Thus, the thermal flow intensity factors (29) and (30) can be written as follows:

Because the crack is symmetrical about the -axis, we get

In Figure 2, the effects of on thermal flow intensity factor for different strip height are displayed (where ). It is shown that with an increase in the thickness of the material plate, the size of the thermal strength factor decreases gradually and tends to stabilize. The results show that the thickness of the sheet has a great influence on the thermal strength factor, which is consistent with the conclusion of Wang and Noda [13].

4.2. Thermal Intensity Factors of Fine-Grained Piezoelectric Coating/Substrate Structure
4.2.1. Studies of the Effect of Elastic Modulus on Thermal Intensity Factor

Here, we select material 1, which is a fine-grained coating made of cadmium selenide piezoelectric ceramic; material 2 is normal cadmium selenide piezoelectric ceramic (where 1 and 2 represent the fine-grained coating and substrate material, respectively). This means that material 1 and material 2 have approximately the same piezoelectric and dielectric constants under suitable conditions, and the elastic modulus of material 1 is larger than that of material 2.

Figures 35 show the effects of the elastic modulus on the thermal intensity factor for different (where , , , , and ). These data indicate that with an increase in the crack length, the thermal stress intensity factor and the electric displacement intensity factor increase, and with the increase of the ratio of the elastic modulus between the fine-grained piezoelectric coating and substrate. The peaks of the thermal stress intensity factor and the electric displacement intensity factor increase correspondingly with the increase in the ratio of the elastic modulus between the fine-grained piezoelectric coating and substrate. Concurrently, when increases from 2, the growth rates of and are greater than that of , and the interface fracture will eventually occur with a continuous increase in .

4.2.2. Studies of the Effect of the Thickness of Coating on Thermal Intensity Factor

Figures 68 illustrate the effects of the thickness of the coating on thermal intensity factor for different (where , , , , and ). This shows that with an increase in crack length, the thermal stress intensity factor and the electric displacement intensity factor increase monotonically. The thermal stress intensity factor and the electric displacement intensity factor have a higher peak value when the coating thickness is small. However, when , the influence of on is not significant enough. The influence of on is obviously significant with an increase of .

5. Conclusions

The Griffith crack problem of fine-grained piezoelectric coating/substrates under steady-state thermal loading is studied. Based on the principle of integral transformation and superposition of equation solutions, the thermomechanical problem is transformed into a solution of thermal stress and electric displacement intensity factors. We then analyzed the effects of crack geometry, coating thickness, and material modulus of elasticity on the fracture behavior of the crack under thermal load. The results show that a larger elastic modulus and thinner coating thickness improve the safety of the coating/substrate structure. This reflects the advantages of fine-grained structure materials.

Appendix

A.

B.

C.

Data Availability

The program codes that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no competing interests.

Acknowledgments

This study was funded by the National Natural Science Foundation of China (11972019).