Transient analysis of dynamic crack propagation in piezoelectric materials

Abstract

In this paper, the transient analysis of semi‐infinite propagating cracks in piezoelectric materials subjected to dynamic anti‐plane concentrated body force is investigated. The crack surface is assumed to be covered with an infinitesimally thin, perfectly conducting electrode that is grounded. In analyzing this problem, it has characteristic lengths and a direct attempt towards solving this problem by transform and Wiener‐Hopf techniques (Noble, 1958) is not applicable. In order to solve this problem, a new fundamental solution for propagating cracks in piezoelectric materials is first established and the transient response of the propagating crack is obtained by superposition of the fundamental solution in the Laplace transform domain. The fundamental solution to be used is the responses of applying exponentially distributed traction in the Laplace transform domain on the propagating crack surface. Taking into account the quasi‐static approximation, exact analytical transient solutions for the dynamic stress intensity factor and the dynamic electric displacement intensity factor are obtained by using the Cagniard‐de Hoop method (Cagnard, 1939; de Hoop, 1960) of Laplace inversion and are expressed in explicit forms. Numerical calculations of dynamic intensity factors are evaluated and the results are discussed in detail. The transient solutions for stationary cracks have been shown to approach the corresponding static values after the shear wave of the piezoelectric material has passed the crack tip.     

Effect of electromechanical coupling on the dynamic interaction of cracks in piezoelectric materials

Summary This article provides a comprehensive treatment of the dynamic interaction between two arbitrarily located and oriented cracks in a piezoelectric medium under steady-state inplane electrical and antiplane mechanical loads. Using an impermeable condition along the crack surfaces, a fundamental dynamic solution was developed for the single crack problem. In this fundamental solution, the single crack problem was treated using Fourier transform and the appropriate singular integral equations. The fundamental solution was then implemented into a pseudo-incident wave method to account for the interaction between the cracks. Numerical examples are provided to show the effect of the geometry of the cracks, the material constants, the frequency of the incident wave and the applied electrical field upon the dynamic stress intensity factors. The results show the significant effect of electromechanical coupling upon the stress intensity factor at the crack tip.     

Finite element analysis of piezoelectric corner configurations and cracks accounting for different electrical permeabilities

Based on eigenfunctions of asymptotic singular electro-elastic fields obtained from a kind of ad hoc finite element method [Chen MC, Zhu JJ, Sze KY. Finite element analysis of piezoelectric elasticity with singular inplane electroelastic fields. Engng Fract Mech 2006;73(7):855-68], a super corner-tip element model is established from the generalized Hellinger-Reissner variational functional and then incorporated into the regular hybrid-stress finite element to determine the coefficients of asymptotic singular electro-elastic fields near a corner-tip. The focus of this paper is not to discuss the well-known behavior of electrically impermeable and permeable (usually it means fully permeable, hereinafter the same) cracks but analyze the limited permeable crack-like corner configurations embedded in the piezoelectric materials, i.e., study the influence of a dielectric medium inside the corner on the singular electro-elastic fields near the corner-tip. The boundary conditions of the impermeable or permeable corner can be considered as simple approximations representing upper and lower bounds for the electrical energy penetrating the corner. Benchmark examples on the piezoelectric crack problems show that present method yields satisfactory results with fewer elements than existing finite element methods do. As application, a piezoelectric corner configuration accounting for the limited permeable boundary condition is investigated, and it is found that the limited permeable assumption is necessary for corners with very small notch angles.     

Finite Element Techniques for Dynamic Crack Analysis in Piezoelectrics

This paper deals with numerical methods, developed to analyze plane stationary cracks in piezoelectric structures under dynamic electromechanical loading conditions. In the first part an explicit finite element scheme is presented, which has been developed to solve the transient coupled electromechanical boundary value problem. A special technique is implemented in the algorithm, accounting for the limited electrical permeability of the crack. In contrast to well known algorithms for static calculations it does not require any iteration. In order to calculate dynamic stress and electric displacement intensity factors for arbitrary crack configurations, the interaction integral is generalized for electromechanical problems. The efficient applicability and the high accuracy of the implementations are demonstrated by numerical examples, giving insight into several effects occuring with dynamically loaded cracks in piezoelectrics     

Electro-elastic stress analysis for a wedge-shaped crack interacting with a screw dislocation in piezoelectric solid

The electro-elastic stress investigation on the interaction problem of a screw dislocation near the tip of a semi-infinite wedge-shaped crack in piezoelectric material has been carried out. Explicit closed-form analytical solutions are obtained for the stress intensity factor (SIF) and the electric displacement intensity factor (EDIF) of the crack, as well as the force on dislocation. The derivation is based on the conformal mapping method and the perturbation technique. The dislocation has Burger's vector normal to the isotropic basal plane, with a line force and a line charge being applied at the core of the dislocation. The influence of the location and the wedge angle of the crack on the image force of the dislocation has been discussed in detail. At the same time the effect of the dislocation on the crack behavior has been also examined under different configurations. Two types of PZT materials are used to numerically illustrate the influences of the wedge angle and the location of the dislocation on the image force and the crack intensities. Results obtained in the current study can be fully reduced to various special cases available in the literatures.     

Electro-elastic behavior induced by an external circular crack in a piezoelectric material

The electro-elastic problem of a transversely isotropic piezoelectric material with a flat crack occupying the outside of a circle perpendicular to the poling axis is considered in this paper. By using the Hankel transform technique, a mixed boundary value problem associated with the considered problem is solved analytically. The results are presented in closed form both for impermeable crack and for permeable crack. A full field solution is given, i.e., explicit expressions for electro-elastic field at any point in the entire piezoelectric space, as well as field intensity factors near the crack front, are determined. A numerical example for a cracked PZT-5H ceramic is given, and the effects of applied electric fields on elastic and electric behaviors are presented graphically.     

Electromechanical model of periodic cracks in piezoelectric materials

The electro-elastic problem for a periodic array of cracks in a piezoelectric medium subjected to coupled electro-mechanical loads is investigated. The mixed boundary value problem, which is formulated directly in terms of the crack surface displacements and electrical potentials, results in a system of hyper-singular integral equations in which the unknown functions are the crack surface displacement and electric potential. Numerical results include the crack surface displacement and the stress and electric intensity factors for the entire range of possible periodic crack spacing and medium size. The central contribution of this paper is the development of an analytical model that predicts crack-spacing effect. The resulting model is validated by a 2D finite element analysis.     

BEM for crack-hole problems in thermopiezoelectric materials

The boundary element formulation for analysing interaction between a hole and multiple cracks in piezoelectric materials is presented. Using Green's function for hole problems and variational principle, a boundary element model (BEM) for a 2-D thermopiezoelectric solid with cracks and holes has been developed and used to calculate stress intensity factors of the crack-hole problem. In BEM, the boundary condition on the hole circumference is satisfied a priori by Green's function, and is not involved in the boundary element equations. The method is applicable to multiple-crack problems in both finite and infinite solids. Numerical results for stress and electric displacement intensity factors at a particular crack tip in a crack-hole system of piezoelectric materials are presented to illustrate the application of the proposed formulation.     

Coupled behaviour of interacting dielectric cracks in piezoelectric materials

Existing studies indicate that the commonly used electrically impermeable and permeable crack models may be inadequate in evaluating the fracture behaviour of piezoelectric materials in some cases. In this paper, a dielectric crack model based on the real electric boundary condition is used to study the electromechanical behaviour of interacting cracks arbitrarily oriented in an infinite piezoelectric medium. The electric boundary condition along the crack surfaces is governed by the opening displacement of the cracks. The formulation of this nonlinear problem is based on modelling the cracks using distributed dislocations and solving the resulting nonlinear singular integral equations using Chebyshev polynomials. Numerical simulation is conducted to show the effect of crack orientation, crack interaction and electric boundary condition upon the fracture behaviour of cracked piezoelectric media.     

Closed form solutions for partially debonded circular inclusion in piezoelectric materials

Summary This paper deals with the problem of a partially debonded piezoelectric circular inclusion in a piezoelectric matrix. This boundary value problem is reduced to two Riemann-Hilbert problems through the use of the analytical continuation theory.Closed form solutions are obtained by considering the behavior of the complex field potentials at origin and infinity. The formulae for the electro-elastic field intensity factors of the interfacial crack are derivedexplicitly. Several particular cases are provided to show the effect of the crack angle, the mechanical and electrical properties and the loads on the electroelastic field singularities.     

The interaction of multiple rows of periodical cracks

In this paper, the interaction of multiple rows of periodical cracks contained in an infinite elastic plate with far-field stress loaded is studied. An extremely accurate and efficient numerical method for solving the problem is presented. The method is mainly by means of the crack isolating analysis technique, stress superposition principle, the Chebyshev polynomial expansion of the pseudo-traction as well as the segmental average collocation technique. This method can be used to compute the stress intensity factors of multiple cracks, periodical cracks, and multiple rows of periodical cracks. In the process of dealing with the superposition of interaction of infinite number of periodic cracks, a key series summation technique is used, which aims at numerical results with extremely high accuracy but with less computation work. Many complex computing examples are given in this paper, and, for some typical examples, numerical results are compared with analytic solutions and with previous numerical solutions. For the problem of the one periodical collinear cracks, the accuracy given by this method reaches to 6 significant digits if a/d 0.9 (where a is the half crack length, and d is the half crack spacing). And even if a/d=0.99, the error is still less than 0.5%. The computer results for multiple rows of periodic collinear and echelon cracks show that the interaction effect between two rows rapidly decline with exponential law as the array pitch increases. This method has filled the gaps in the research field on the interaction of multiple rows of general periodical cracks.     

含裂纹压电材料的Cell-Based光滑扩展有限元法

为精确而有效地求解机电耦合作用下含裂纹压电材料的断裂参数,首先,通过将复势函数法、扩展有限元法和光滑梯度技术引入到含裂纹压电材料的断裂机理问题中,提出了含裂纹压电材料的Cell-Based光滑扩展有限元法;然后,对含中心裂纹的压电材料强度因子进行了模拟,并将模拟结果与扩展有限元法和有限元法的计算结果进行了对比。数值算例结果表明:Cell-Based光滑扩展有限元法兼具扩展有限元法和光滑有限元法的特点,不仅单元网格与裂纹面相互独立,且裂尖处单元不需精密划分,与此同时,Cell-Based光滑扩展有限元法还具有形函数简单且不需求导、对网格质量要求低且求解精度高等优点。所得结论表明Cell-Based光滑扩展有限元法是压电材料断裂分析的有效数值方法。      

Identification of multiple cracks based on optimization methods using harmonic elastic waves

A theoretical study of the identification of multiple cracks in an elastic medium based on optimization methods using harmonic elastic waves is presented in the paper. A general interacting crack model is first used to determine the dynamic interaction between arbitrarily located and oriented cracks subjected to plane harmonic waves. The solution of this problem is then implemented into an optimization process for the identification of unknown cracks from known strain components at discrete locations. An optimization scheme based on sensitivity analysis is used to determine the length, the orientation and the position of cracks. Numerical simulation indicates that the sensitivity analysis and the optimization method used are effective in identifying multiple cracks. It is observed that convergent results can be achieved from a set of arbitrarily determined initial values of the crack parameters. Numerical examples are presented to illustrate the determination of the length, orientation, and position of different interacting cracks.     

Interaction of a penny-shaped crack with an elliptic crack under shear loading

The problem of interaction between equal coplanar elliptic cracks embedded in a homogeneous isotropic elastic medium and subjected to shear loading was solved analytically by Saha et al. (1999) International Journal of Solids and Structures 36, 619–637, using an integral equation method. In the present study the same integral equation method has been used to solve the title problem. Analytical expression for the two tangential crack opening displacement potentials have been obtained as series in terms of the crack separation parameter _i up to the order _i⁵,(i=1,2) for both the elliptic as well as penny-shaped crack. Expressions for modes II and III stress intensity factors have been given for both the cracks. The present solution may be treated as benchmark to solutions of similar problems obtained by various numerical methods developed recently. The analytical results may be used to obtain solutions for interaction between macro elliptic crack and micro penny-shaped crack or vice-versa when the cracks are subjected to shear loading and are not too close. Numerical results of the stress-intensity magnification factor has been illustrated graphically for different aspect ratios, crack sizes, crack separations, Poisson ratios and loading angles. Also the present results have been compared with the existing results of Kachanov and Laures (1989) International Journal of Fracture 41, 289–313, for equal penny-shaped cracks and illustrations have been given also for the special case of interaction between unequal penny-shaped cracks subjected to uniform shear loading.     

3D point force solution for a permeable penny-shaped crack embedded in an infinite transversely isotropic piezoelectric medium

This paper considers the non-axisymmetric three-dimensional problem of a penny-shaped crack with permeable electric conditions imposed on the crack surfaces, subjected to a pair of point normal forces applied symmetrically with respect to the crack plane. The crack is embedded in an infinite transversely isotropic piezoelectric body with the crack face perpendicular to the axis of material symmetry. Applying the symmetry of the problem under consideration then leads to a mixed–mixed boundary value problem of a half-space, for which potential theory method is employed for the purpose of analysis. The cases of equal eigenvalues are also discussed. Although the treatment differs from that for an impermeable crack reported in literature, the resulting governing equation still has a familiar structure. For the case of a point force, exact expressions for the full-space electro-elastic field are derived in terms of elementary functions with explicit stress and electric displacement intensity factors presented. The exact solution for a uniform loading is also given.     

Active constrained layer damping of geometrically nonlinear vibrations of smart laminated composite sandwich plates using 1–3 piezoelectric composites

This paper deals with the analysis of active constrained layer damping (ACLD) of geometrically nonlinear vibrations of sandwich plate with orthotropic laminated composite faces separated by a flexible core. The constraining layer of the ACLD treatment is composed of the vertically/obliquely reinforced 1?C3 piezoelectric composites. The Golla?CHughes?CMcTavish method has been implemented to model the constrained viscoelastic layer of the ACLD treatment in time domain. The first-order shear deformation theory and the Von Kármán type nonlinear strain displacement relations are used for analyzing this coupled electro-elastic problem. A three dimensional finite element model of smart laminated composite sandwich plate integrated with ACLD patches has been developed to investigate the performance of these patches for controlling the geometrically nonlinear vibrations of the plates. The numerical results indicate that the ACLD patches significantly improve the damping characteristics of the sandwich plates with laminated cross-ply and angle-ply facings for suppressing their geometrically nonlinear vibrations. Particular emphasis has been placed on investigating the effect of the variation of piezoelectric fiber orientation angle on the performance of the ACLD treatment.     

Electro-elastic stress analysis for a screw dislocation interacting with a coated inclusion in piezoelectric solid

Summary. The electro-elastic stress investigation on the interaction problem of a piezoelectric screw dislocation near a coated inclusion in a piezoelectric material has been carried out. In our study, three dissimilar material phases are involved: the matrix, the inclusion and the coating layer. All the three materials are piezoelectric and with different material constants. Explicit closed-form analytical solutions for the stress and electric displacement fields are obtained by using the complex variable method. The image force acting on the screw dislocation is calculated by using the generalized Peach-Koehler formula. Numerical examples for different material constant combinations are performed. The influences of material properties of the inclusion and the coating layer on the image forces are examined and discussed.     

X-FEM simulation for two-unequal-collinear cracks in 2-D finite piezoelectric specimen

A study for two-unequal-collinear cracks in a 2-D finite piezoelectric specimen is carried out using a new set of six crack-tip enrichment functions proposed here for piezoelectric media in the X-FEM framework. The intensity factors and energy release rate are calculated using interaction integral in conjugation with the near tip behavior given by the Stroh formalism. Effect of finite size of the specimen is analyzed with respect to offset distances of the cracks from the specimen boundaries. ERR variations are investigated with respect to inter- crack space, crack lengths and electrical/mechanical loadings. Hence, two-unequal-collinear cracks in an infinite domain problem is simulated, analyzed and validated using X-FEM. Further, ERR at the crack tips for the asymmetric case of two collinear equal and unequal cracks, is also computed. It is concluded through this investigation that the proposed enrichment functions could be used to handle the problems of fracture mechanics in 2-D piezoelectric media within a good accuracy.     

Thermoelectroelastic analysis of cracks in piezoelectric half-plane by BEM

The Green's function and the boundary element method for analysing fracture behaviour of cracks in piezoelectric half-plane are presented in this paper. By combining Stroh formalism and the concept of perturbation, a general thermoelectroelastic solution for half-plane solid subjected to point heat source and/or temperature discontinuity has been derived. Using the proposed solution and the potential variational principle, a boundary element model (BEM) for 2-D half-plane solid with multiple cracks has been developed and used to calculate the stress intensity factors of the multiple crack problem. The method is available for multiple crack problems in both finite and infinite solids. Numerical results for a two-crack system are presented and compared with those from finite element method (FEM).