Dynamic behavior of a crack in a functionally graded piezoelectric strip bonded to two dissimilar half piezoelectric material planes

Summary. The dynamic behavior of a crack in a functionally graded piezoelectric material (FGPM) strip bonded to two half dissimilar piezoelectric material planes subjected to combined harmonic anti-plane shear wave and in-plane electrical loading was studied under the limited permeable and permeable electric boundary conditions. It was assumed that the elastic stiffness, piezoelectric constant and dielectric permittivity of the functionally graded piezoelectric layer vary continuously along the thickness of the strip. By using the Fourier transform, the problem can be solved with a set of dual integral equations in which the unknown variables are the jumps of the displacements and the electric potentials across the crack surfaces. In solving the dual integral equations, the jumps of the displacements and the electric potentials across the crack surfaces were expanded in a series of Jacobi polynomials. Numerical results illustrate the effects of the gradient parameter of FGPM, electric loading, wave number, thickness of FGPM strip and electric boundary conditions on the dynamic stress intensity factors (SIFs).     

在反平面变形下功能梯度压电板条的应力强度因子分析

利用Fourier变换技术将混合边值问题化为对偶积分方程,求解对偶积分方程得到应力强度因子以及电位移强度因子的表达式。最后通过数值计算讨论了材料参数、载荷条件以及裂纹的几何参数等对功能梯度压电材料中裂纹尖端应力强度因子的影响。     

Dynamic stress intensity factors of two 3D rectangular permeable cracks in a transversely isotropic piezoelectric material under a time‐harmonic elastic P‐wave

The dynamic behavior of two 3D rectangular permeable cracks in a transversely isotropic piezoelectric material is investigated under an incident harmonic stress wave by using the generalized Almansi's theorem and the Schmidt method. The problem is formulated through double Fourier transform into three pairs of dual integral equations with the displacement jumps across the crack surfaces as the unknown variables. To solve the dual integral equations, the displacement jumps across the crack surfaces are directly expanded as a series of Jacobi polynomials. Finally, the relations among the dynamic stress field and the dynamic electric displacement filed near the crack edges are obtained, and the effects of the shape of the rectangular crack, the characteristics of the harmonic wave, and the distance between two rectangular cracks on the stress and the electric intensity factors in a piezoelectric composite material are analyzed. Copyright © 2015 John Wiley & Sons, Ltd.     

Crack propagating in a functionally graded strip under the plane loading

In the present paper, a finite crack with constant length (Yoffe type crack) propagating in the functionally graded strip under the plane loading is investigated by means of the Schmidt method. By using the Fourier transform and defining the jumps of displacement components across the crack surface as the unknown functions, two pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement components across the crack surface are expanded in a series of Jacobi polynomial. Numerical examples are provided to show the effects of the material properties, the thickness of the functionally graded strip, and speed of the crack propagating upon the dynamic fracture behavior.     

On the Moving Griffith Crack in a Nonhomogeneous Orthotropic Strip

A finite crack with constant length (Yoffe type crack) propagating in the functionally graded orthotropic strip under the plane loading is investigated by means of the Schmidt method. By using the Fourier transform and defining the jumps of displacement components across the crack surface as the unknown functions, two pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement components across the crack surface are expanded in a series of Jacobi polynomial. Numerical examples are provided to show the effects of material properties, the thickness of the functionally graded orthotropic strip and the speed of the crack propagating upon the dynamic fracture behavior.     

Non-local theory solution of a mode-I crack in a piezoelectric/piezomagnetic composite material plane

The non-local theory solution of a mode-I permeable crack in a piezoelectric/piezomagnetic composite material plane was given by using the generalized Almansi’s theorem and the Schmidt method in this paper. The problem was formulated through Fourier transform into two pairs of dual integral equations, in which the unknown variables are the displacement jumps across the crack surfaces. To solve the dual integral equations, the displacement jumps across the crack surfaces were directly expanded as a series of Jacobi polynomials. Numerical examples were provided to show the effects of the crack length and the lattice parameter on the stress field, the electric displacement field and the magnetic flux field near the crack tips. Unlike the classical elasticity solutions, it is found that no stress, electric displacement and magnetic flux singularities are present at the crack tips in piezoelectric/piezomagnetic composite materials. The non-local elastic solution yields a finite hoop stress at the crack tip, thus allowing us to use the maximum stress as a fracture criterion.     

An interface crack in a functionally graded piezoelectric bi-layer under anti-plane shear impact

The transient response of an interface crack between two dissimilar functionally graded piezoelectric material (FGPM) layers under anti-plane shear impact loading is analyzed using the integral transform method. The properties of the FGPM layers vary continuously along the thickness, and the two layers are connected weak-discontinuously. Laplace transform and Fourier transform are used to reduce the problem to two sets of dual integral equations, which are then expressed to the Fredholm integral equations of the second kind. Numerical values on the dynamic energy release rate are presented for the FGPM to show the effects on the electric loading, variation and gradient of material properties, and thickness of layers. Following things are helpful to increase the resistance of transient fracture of interface crack in FGPMs: (a) increase of the material properties from the interface to the upper or lower free surface; (b) decrease of weak discontinuity at the interface; (c) increase of the gradient of material properties; (d) certain direction and magnitude of the electric loading; and (e) increase of the thickness of the FGPM layer.     

A penny-shaped crack in a functionally graded piezoelectric strip under thermal loading

The mixed-mode thermoelectromechanical fracture problem for a functionally graded piezoelectric material (FGPM) strip with a penny-shaped crack is considered. It is assumed that the thermoelectroelastic properties of the strip vary continuously along the thickness of the strip, and that the strip is under thermal loading. The crack faces are supposed to be insulated thermally and electrically. The thermal and electromechanical problems are reduced to singular integral equations and solved numerically. The stress and electric displacement intensity factors are presented for different crack size, crack position and material nonhomogeneity.     

Investigation of anti-plane shear behavior of a Griffith crack in a piezoelectric material by using the non-local theory

In this paper, the behavior of a Griffith crack in a piezoelectric material under anti-plane shear loading is investigated by using the non-local theory for impermeable crack surface conditions. By using the Fourier transform, the problem can be solved with two pairs of dual integral equations. These equations are solved using Schmidt method. Numerical examples are provided. Contrary to the previous results, it is found that no stress and electric displacement singularity is present at the crack tip.     

Non-Local Theory Solution for an Anti-Plane Shear Permeable Crack in Functionally Graded Piezoelectric Materials

In this paper, the non-local theory solution of a Griffith crack in functionally graded piezoelectric materials under the anti-plane shear loading is obtained for the permeable electric boundary conditions, in which the material properties vary exponentially with coordinate parallel to the crack. The present problem can be solved by using the Fourier transform and the technique of dual integral equation, in which the unknown variable is the jump of displacement across the crack surfaces, not the dislocation density function. To solve the dual integral equations, the jump of the displacement across the crack surfaces is directly expanded in a series of Jacobi polynomials. From the solution of the present paper, it is found that no stress and electric displacement singularities are present near the crack tips. The stress fields are finite near the crack tips, thus allows us to use the maximum stress as a fracture criterion. The finite stresses and the electric displacements at the crack tips depend on the crack length, the functionally graded parameter and the lattice parameter of the materials, respectively. On the other hand, the angular variations of the strain energy density function are examined to associate their stationary value with locations of possible fracture initiation.     

Geometrically nonlinear static and free vibration analysis of functionally graded piezoelectric plates

In this paper, nonlinear static and free vibration analysis of functionally graded piezoelectric plates has been carried out using finite element method under different sets of mechanical and electrical loadings. The plate with functionally graded piezoelectric material (FGPM) is assumed to be graded through the thickness by a simple power law distribution in terms of the volume fractions of the constituents. Only the geometrical nonlinearity has been taken into account and electric potential is assumed to be quadratic across the FGPM plate thickness. The governing equations are obtained using potential energy and Hamilton’s principle that includes elastic and piezoelectric effects. The finite element model is derived based on constitutive equation of piezoelectric material accounting for coupling between elasticity and electric effect using higher order plate elements. The present finite element is modeled with displacement components and electric potential as nodal degrees of freedom. Results are presented for two constituent FGPM plate under different mechanical boundary conditions. Numerical results for PZT-4/PZT-5H plate are given in dimensionless graphical forms. Effects of material composition and boundary conditions on nonlinear response are also studied. The numerical results obtained by the present model are in good agreement with the available solutions reported in the literature.     

On the dynamic propagation of an anti-plane shear crack in a functionally graded piezoelectric strip

Summary. In this paper, a finite crack propagating at constant speed in a functionally graded piezoelectric material (FGPM) is studied. It is assumed that the electroelastic material properties of the FGPM vary continuously according to exponential gradients along the thickness of the strip, and that the strip is under anti-plane shear mechanical and in-plane electrical loads. The analysis is conducted on the electrically unified (natural) crack boundary condition, which is related to the ellipsoidal crack parameters. By using the Fourier transform, the problem is reduced to the solutions of Fredholm integral equations of the second kind. Numerical results for the stress intensity factor and crack sliding displacement are presented to show the influences of the elliptic crack parameters, crack propagation speed, electric field, FGPM gradation, crack length, and electromechanical coupling coefficient. It reveals that there are considerable differences between traditional electric crack models and the present unified crack model.     

Transient response of a functionally graded piezoelectric strip with a penny-shaped crack under electric time-dependent loading

Summary The dynamic fracture problem for a functionally graded piezoelectric material (FGPM) strip containing a penny-shaped crack parallel to the free boundaries is considered in this study. It is assumed that the electroelastic properties of the strip vary continuously along the thickness direction of the strip, and that the strip is under time-dependent electric load. Integral transform techniques and dislocation density functions are employed to reduce the problem to the solutions of a system of singular integral equations. The stress and electric displacement intensity factors versus time are presented for various values of dimensionless parameters representing the crack size, the crack location and the material nonhomogeneity.     

Fracture parameters of a penny-shaped crack at the interface of a piezoelectric bi-material system

A penny-shaped crack at the interface of a piezoelectric bi-material system is considered. Analytical general solutions based on Hankel integral transforms are used to formulate the mixed-boundary value problem corresponding to an interfacial crack and the problem is reduced to a system of singular integral equations. The integral equations are further reduced to two systems of algebraic equations with the aid of Jacobi polynomials and Chebyshev polynomials. Thereafter, the exact expressions for the stress intensity factors and the electric displacement intensity factor at the tip of a crack are obtained. Selected numerical results are presented for various bi-material systems to portray the significant features of crack tip fracture parameters and their dependence on material properties, poling orientation and electric loading.     

Interface crack problem of functionally graded piezoelectric materials: effects of the position of electromechanical impact and gradient

This paper analyzes the transient response of the dynamic stress intensity factor for an interfacial crack of a functionally graded piezoelectric material (FGPM) coated on the surface of a homogeneous piezoelectric substrate. Different from previous analyses, this study mainly considers a realistic situation when electromechanical loadings are suddenly applied at the material surface. Obtained results are compared with those when the crack surfaces are directly loaded by the same impacts. By using the integral transform method, the problem is reduced to solving two singular integral equations. It is found that dynamic stress intensity factors are significantly amplified and reduced depending on the negative and positive gradient for electromechanical impacts at the material surface. Positive or negative electric impact also decreases or increases the overshoot of the dynamic stress intensity factor. It is suggested that designing an FGPM with a positive gradient index is safer than a negative gradient index.     

Basic solutions of two parallel mode-I cracks or four parallel mode-I cracks in the piezoelectric materials

In this paper, the stress and the electric intensity factors of two parallel mode-I cracks or four parallel mode-I cracks in the piezoelectric materials were examined by means of the Schmidt method for the permeable electric boundary conditions. The present problem can be solved by using the Fourier transform and the technique of dual integral equation, in which the unknown variables are the jumps of displacements across the crack surfaces, not the dislocation density functions. To solve the dual integral equations, the displacement jumps are directly expanded in a series of Jacobi polynomials. Finally, the effects of the distance between two parallel cracks and the distance between two collinear cracks on the stress and the electric intensity factors in the piezoelectric materials are analyzed. These results can be used for the strength evaluation of the piezoelectric materials with multi-cracks.     

Dugdale plastic zone of a penny-shaped crack in a piezoelectric material under axisymmetric loading

The Dugdale plastic zone ahead of a penny-shaped crack in a piezoelectric material, subjected to electric and axisymmetric mechanical loadings, is evaluated analytically. Hankel transform is employed to reduce the mixed boundary-value problem of the penny-shaped crack to dual integral equations, which are solved exactly under the assumption of electrically permeable crack face conditions. A closed-form solution to the mixed boundary-value problem is obtained to predict the relationship between the length of the plastic zone and the applied loading. The stress distribution in and outside of the yield zone has been derived analytically, and the crack opening displacement has been investigated. The electric displacement has a constant value in the strip yield zone. The current Dugdale crack model leads to non-singular stress and electric fields near the crack front, and it is observed that the material properties affect the crack opening displacement.     

Transient thermoelectroelastic response of a functionally graded piezoelectric strip with a penny-shaped crack

Transient response of a penny-shaped crack in a plate of a functionally graded piezoelectric material (FGPM) is studied under thermal shock loading conditions. It is assumed that the thermoelectroelastic properties of the strip vary continuously along the thickness of the strip, and that the crack faces are completely insulated. By using both the Laplace and Hankel transforms, the thermal and electromechanical problems are reduced to a singular integral equation and a system of singular integral equations which are solved numerically. The intensity factors vs. time for various crack size, crack position and material nonhomogeneity are obtained.     

Mechanical and electrical yielding for an electrically insulated crack in an interlayer between piezoelectric materials

A plane problem for a crack in a thin ductile layer between two piezoelectric materials under remote electromechanical loading is considered. The piezoelectric substrates are of the same material. They are modeled by half-spaces. For the interlayer thickness tending to zero, electrical and mechanical yielding zones of different unknown lengths are introduced at the crack continuations. These zones are modeled by jumps of electric potential and mechanical displacement jumps, respectively. A vector Hilbert problem is formulated and solved exactly. The lengths of the yield zones are found from the conditions of finiteness of the normal stress and of the electrical displacement. The crack opening and the electric potential jumps at the initial crack tip are found in closed form for the cases of electrical yield zone longer and shorter than the mechanical one. Simple formulae for energy release rate are presented for Griffith crack and the developed electromechanical yield models. Also, an approximate model for eliminating the mechanical and electrical singularities is suggested. Numerical results are presented for different values of yield parameters and of remote electromechanical loading. A good agreement between the energy release rate values obtained for a Griffith crack, exact and approximate yield zone models are established.     

Electroelastic analysis for a Griffith crack interacting with a coated inclusion in piezoelectric solid

A Mode III Griffith crack interacting with a coated inclusion in piezoelectric media is investigated. The crack, the coated inclusion are embedded in an infinitely extended piezoelectric matrix media, with the crack being along the radial direction of the inclusion. In the study, three different piezoelectric material phases are involved: the inclusion, the coating layer, and the matrix. A far-field loading condition is considered. During the solution procedure, the crack is simulated as a continuous distribution of screw dislocations. By using the solution of a screw dislocation near a coated inclusion in piezoelectric media as the Green function, the problem is formulated into a set of singular integral equations, which are solved by numerical method. The stress and electric displacement intensity factors are derived in terms of the asymptotic values of the dislocation density functions evaluated from the integral equations. Numerical examples are given for various material constants combinations and geometric parameters.