Surface cracking of a piezoelectric strip bonded to an elastic substrate (Mode I crack problem)

Summary A piezoelectric material layer bonded to an elastic substrate is investigated. The piezoelectric layer contains an edge crack that is perpendicular to the surface of medium. The poling axis of the piezoelectric layer is parallel to its edge. The elastic layer can be an ideal insulator or an ideal conductor. The Fourier transform technique is used to reduce the problem to a solution of singular integral equations. Both impermeable crack and permeable crack assumptions are considered. Stress and electric displacement intensity factors are investigated for different dimensions of the medium. A double-edge cracked symmetric piezoelectric laminate under symmetric electro-mechanical load is also investigated. BLW would like to thank the National Science Foundation of China (#10102004) and the City University of Hong Kong for the support of this work (DAG #7100219). YGS also thanks the Multidiscipline Scientific Research Foundation Project (HIT. MD 2001. 39) of the Harbin Institute of Technology and the SRF for ROCS, SEM.     

Anti-plane crack moving along the interface of dissimilar piezoelectric materials

The problem of an anti-plane Griffith crack moving along the interface of dissimilar piezoelectric materials is solved by using the integral transform technique. It is shown from the result that the intensity factors of anti-plane stress and electric displacement around the crack tip are dependent on the speed of the Griffith crack as well as the material coefficients. When the two piezoelectric materials are identical, the present result will be reduced to the result for the problem of an anti-plane moving Griffith crack in homogeneous piezoelectric materials. Supported by the National Natural Science Foundation and the National Post-doctoral Science Foundation of China.     

The behavior of a crack in functionally graded piezoelectric/piezomagnetic materials under anti-plane shear loading

Summary In this paper, the behavior of a crack in functionally graded piezoelectric/piezomagnetic materials subjected to an anti-plane shear loading is investigated. To make the analysis tractable, it is assumed that the material properties vary exponentially with the coordinate parallel to the crack. By using a Fourier transform, the problem can be solved with the help of a pair of dual integral equations in which the unknown variable is the jump of the displacements across the crack surfaces. These equations are solved using the Schmidt method. The relations among the electric displacement, the magnetic flux and the stress field near the crack tips are obtained. Numerical examples are provided to show the effect of the functionally graded parameter on the stress intensity factors of the crack.The authors are grateful for financial support from the Natural Science Foundation of Hei Long Jiang Province (A0301), the National Natural Science Foundation of China (50232030, 10172030), the Natural Science Foundation with Excellent Young Investigators of Hei Long Jiang Province(JC04-08) and the National Science Foundation with Excellent Young Investigators (10325208).     

Near-tip fields for a crack growing along an elastic perfectly plastic interface

The stress and deformation fields near the tip of an anti-plane crack growing quasi-statically along an interface of elastic perfectly plastic materials are given in this paper. A family of solutions for the growing crack fields is found covering all admissible crack line shear stress ratios. The project supported by the National Natural Science Foundation of China     

Circular crack in a transversely isotropic piezoelectric space under point forces and point charges

In this paper, two kinds of circular crack including external circular crack and penny-shaped crack in a transversely isotropic piezoelectric space are considered. Firstly, we obtain the solution to the problem of an external circular crack in a transversely isotropic piezoelectric space subjected to antisymmetric normal point forces and point charges. Based on this, the solution of one-sided loading of an external circular crack is constructed. Secondly, the real shape of an external circular crack and the opening displacement of a penny-shaped crack under an arbitrary point force and point charge are further obtained. At last, the results are presented in a graphical form. The project supported by the National Natural Science Foundation of China (19872060 and 69982009) and the Postdoctoral Foundation of China     

The scattering of harmonic elastic anti-plane shear waves by a finite crack in infinitely long strip using the non-local theory

In this paper, the scattering of harmonic anti-plane shear waves by a finite crack in infinitely long strip is studied using the non-local theory. The Fourier transform is applied and a mixed boundary value problem is formulated. Then a set of dual integral equations is solved using the Schmidt method instead of the first or the second integral equation method. A one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress occurring at the crack tips. Contraty to the classical elasticity solution, it is found that no stress singularity is present at the crack tip. The non-local dynamic elastic solutions yield a finite hoop stress at the crack tip, thus allowing for a fracture criterion based on the maximum dynamic stress hypothesis. The finite hoop stress at the crack tip depends on the crack length, the width of the strip and the lattice parameter. Supported by the Post Doctoral Science Foundation of Heilongjiang Province, the Natural Science Foundation of Heilongjiang Province and the National Foundation for Excellent Young Investigators.     

Dynamic behavior of two parallel symmetric permeable cracks in a piezoelectric strip

The dynamic behavior of two parallel symmetric cracks in a piezoelectric strip under harmonic anti-plane shear waves is studied using the Schmidt method for permeable crack surface conditions. The cracks are parallel to the edge of the strip. By means of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations. These equations are solved using the schmidt method. The results show that the stress and the electric displacement intensity factors depend on the geometry of the cracks, the frequency of incident waves, the distance between cracks and the thickness of the strip. It is also found that the electric displacement intensity factors for the permeable crack surface conditions are much smaller than those for the impermeable crack surface conditions. Project supported by the Post Doctoral Science Foundation of Heilongjiang Province, the Natural Science Foundation of Heilongjiang Province, the National Science Foundation with the Excellent Young Investigator Award (No. 19725209) and the Scientific Research Foundation of Harbin Institute of Technology (HIT.2000.30).     

The singularity behavior of an interface crack under a dynamic load

The singularity behavior of a crack on the interface of two different media under dynamic load is investigated. By introducing a small region in which the crack faces make frictionless contact and making use of a kind of integral equations with moving boundaries, it is proved that there are only square-root singularities near the interface crack tips in case that a dynamic load acts on it. Numerical results show that the normal stress in the contact region remains negative. The results of the stress intensity factor and the length of the crack face contact region are given to illustrate the dynamic behavior of the interface crack.This work is supported by the National Natural Science Foundation of China.     

Dynamic interaction between a sub-interface crack and the interface in a bi-material: time-domain BEM analysis

Summary Transient response of a sub-interface crack in a bi-material is studied with emphasis on the dynamic interaction between the crack and the interface, by combining the traditional time-domain displacement boundary element method (BEM) and the non-hypersingular traction BEM. Computations are performed for an unbounded bi-material with a crack subjected to impact tensile loading on its faces or incident impact waves and a bounded rectangular bi-material plate under remote impact tensile loading. Numerical results of the dynamic stress intensity factors (DSIFs) and dynamic interface tractions are presented for various material combinations and crack locations. It is shown that pronounced increases in DSIFs and the interface tractions may be caused in some cases because of the dynamic interaction between the crack and the interface.This work was initialized during the second author's stay at Institute of Mechanics, TU Darmstadt, Germany under the support of the Alexander von Humboldt Foundation. Discussion on the BEM formulation with Dr. Seelig is gratefully acknowledged. The first two authors are also grateful for the partial support by the China National Natural Science Foundation under Grant No. 10025211 and the NJTU Scientific Paper Fund (PD195).     

Elastic-plastic crack growth on plane strain bimaterial interface

Finite element computation are carried out to simulate plane strain crack growth on a bimaterial interface under the assumption of small scale yielding. The modified Gurson constitutive equation and the element vanish technique introduced by Tvergaard et al. are used to model the final formation of an open crack. It is found from the calculation that the critical fracture toughness for crack growth is much lower in bimaterials than that in homogeneous material. The critical fracture toughness is strongly dependent on material properties of the bimaterial pair and the mixed mode of remote loads. The interface crack grows in the more compliant (lower hardening) material or in the weaker (lower yield strength) material. In Mode-I loading, the crack grows zigzag along the interface. Project supported by Fok Ying-Tung Education Foundation and National Natural Science Foundation of China.     

Asymmetrical dynamic propagation problems on mode III interface crack

By the application of the theory of complex functions, asymmetrical dynamic propagation problems on modeⅢinterface crack are studied. The universal representations of analytical solutions are obtained by the approaches of self-similar function. The problems researched can be facilely transformed into Riemann-Hilbert problems and analytical solution to an asymmetrical propagation crack under the condition of point loads and unit-step loads, respectively, is acquired. After those solutions were used by superposition theorem, the solutions of arbitrarily complex problems could be attained.     

Scattering of love waves by an interface crack between a piezoelectric layer and an elastic substrate

The scattering of Love waves by an interface crack between a piezoelectric layer and an elastic substrate is investigated by using the integral transform and singular integral equation techniques. The dynamic stress intensity factors of the left and the right crack tips are determined. It is found from numerical calculation that the dynamic response of the system depends significantly on the crack configuration, the material combination and the propagating direction of the incident wave. It is expected that specifying an appropriate material combination may retard the growth of the crack for a certain crack configuration. Project supported by the National Natural Science Foundation of China (No. 19891180), the Fundamental Research Foundation of Tsinghua University (JZ 2000.007) and the Fund of the Education Ministry of China.     

Scattering of SH-waves by an interacting interface linear crack and a circular cavity near bimaterial interface

An analytical method is developed for scattering of SH-waves and dynamic stress concentration by an interacting interface crack and a circular cavity near bimaterial interface. A suitable Green‘s function is contructed, which is the fundamental solution of the displacement field for an elastic half space with a circular cavity impacted by an out-plane harmonic line source loading at the horizontal surface. First, the bimaterial media is divided into two parts along the horizontal interface, one is an elastic half space with a circular cavity and the other is a complete half space. Then the problem is solved according to the procedure of combination and by the Green‘s function method. The horizontal surfaces of the two half spaces are loaded with undetermined anti-plane forces in order to satisfy continuity conditions at the linking section, or with some forces to recover cracks by means of crack-division technique. A series of Fredholm integral equations of first kind for determining the unknown forces can be set up through continuity conditions as expressed in terms of the Green‘s function. Moreover, some expressions are given in this paper, such as dynamic stress intensity factor (DSIF) at the tip of the interface crack and dynamic stress concentration factor (DSCF) around the circular cavity edge. Numerical examples are provided to show the influences of the wave numbers, the geometrical location of the interface crack and the circular cavity, and parameter combinations of different media upon DSIF and DSCF.     

Toughening analysis of mode III crack in transformation toughened ceramics

The fracture toughness enhancement of mode III crack in transformation toughened ceramics such as Ce-TZP and Mg-PSZ is predicted by using a pressure sensitive transformation criterion and the weight function method. The theoretical analyses of the toughening effect for both stationary and steady-state growing crack are given, respectively. The results show that the toughening is only associated with crack growth and there is no toughening effect for stationary crack.The project was supported by the National Science Foundation of China     

Griffith crack moving in a piezoelectric strip

Summary  The dynamic problem of an impermeable crack of constant length 2a propagating along a piezoelectric ceramic strip is considered under the action of uniform anti-plane shear stress and uniform electric field. The integral transform technique is employed to reduce the mixed-boundary-value problem to a singular integral equation. For the case of a crack moving in the mid-plane, explicit analytic expressions for the electroelastic field and the field intensity factors are obtained, while for an eccentric crack moving along a piezoelectric strip, numerical results are determined via the Lobatto–Chebyshev collocation method for solving a resulting singular integral equation. The results reveal that the electric-displacement intensity factor is independent of the crack velocity, while other field intensity factors depend on the crack velocity when referred to the moving coordinate system. If the crack velocity vanishes, the present results reduce to those for a stationary crack in a piezoelectric strip. In contrast to the results for a stationary crack, applied stress gives rise to a singular electric field and applied electric field results in a singular stress for a moving crack in a piezoelectric strip. Received 14 August 2001; accepted for publication 24 September 2002 The author is indebted to the AAM Reviewers for their helpful suggestions for improving this paper. The work was supported by the National Natural Science Foundation of China under Grant 70272043.     

Study for 2D moving contact elastic body with closed crack using BEM

Using a sub-regional boundary element method, an algorithm for the two-dimensional elastic bodies with a closed crack loaded by a moving contact elastic body is proposed. Since the extent and status of the contact surface of two elastic bodies and the crack within the body are all not known in advance, a double iterative contact algorithm is used. The BEM program for solving the closed crack problems is developed, some numerical examples are calculated, and the results of the center crack cases are shown to be in good agreement with the analytical solution in the classical fracture mechanics. In the condition of friction and non-friction, some coupling computational results of the SIF for the closed crack, with different angles and loaded by a moving contact elastic body, are also obtained by a numerical computation. The project supported by the National Natural Science Foundation of China (10172053) and NJTU Foundation of China (PD-157)     

A MODE Ⅱ CRACK IN A TWO-DIMENSIONAL OCTAGONAL QUASICRYSTALS

The present study develops the fracture theory for a two-dimensional octagonal quasicrystals. The exact analytic solution of a Mode Ⅱ Griffith crack in the material was obtained by using the Fourier transform and dual integral equations theory, then the displacement and stress fields, stress intensity factor and strain energy release rate were determined, the physical sense of the results relative to phason and the difference between mechanical behaviors of the crack problem in crystal and quasicrystal were figured out. These provide important information for studying the deformation and fracture of the new solid phase.     

On anti-plane shear behavior of a Griffith permeable crack in piezoelectric materials by use of the non-local theory

In this paper, the non-local theory of elasticity is applied to obtain the behavior of a Griffith crack in the piezoelectric materials under anti-plane shear loading for permeable crack surface conditions. By means of the Fourier transform the problem can be solved with the help of a pair of dual integral equations with the unknown variable being the jump of the displacement across the crack surfaces. These equations are solved by the Schmidt method. Numerical examples are provided. Unlike the classical elasticity solutions, it is found that no stress and electric displacement singularity is present at the crack tip. The non-local elastic solutions yield a finite hoop stress at the crack tip, thus allowing for a fracture criterion based on the maximum stress hypothesis. The finite hoop stress at the crack tip depends on the crack length and the lattice parameter of the materials, respectively. The project supported by the National Natural Science Foundation of China (50232030 and 10172030)     

Mode I and mode II crack tip asymptotic fields with strain gradient effects

The strain gradient effect becomes significant when the size of fracture process zone around a crack tip is comparable to the intrinsic material lengthl, typically of the order of microns. Using the new strain gradient deformation theory given by Chen and Wang, the asymptotic fields near a crack tip in an elastic-plastic material with strain gradient effects are investigated. It is established that the dominant strain field is irrotational. For mode I plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist simultaneously. In the stress dominated asymptotic field, the angular distributions of stresses are consistent with the classical plane stress HRR field; In the couple stress dominated asymptotic field, the angular distributions of couple stresses are consistent with that obtained by Huang et al. For mode II plane stress and plane strain crack tip asymptotic fields, only the stress-dominated asymptotic fields exist. The couple stress asymptotic field is less singular than the stress asymptotic fields. The stress asymptotic fields are the same as mode II plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode I and mode II, because the present theory is based only on the rotational gradient of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient. The project supported by the National Natural Science Foundation of China (19704100), National Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20), CAS K.C. Wong Post-doctoral Research Award Fund and Post-doctoral Science Fund of China     

The anti-plane shear field in an infinite slab of elasto-damaged material with a semi-infinite crack

This paper deals with an infinite slab with a semi-infinite crack, which is subjected to the anti-plane sheark _III field at infinity. The slab is made of an elasto-damaged material. Analytical solution is obtained by use of conformal mapping. The shape of damaged-zone, the dissipative energy, the shear opening displacement on the crack surface and several stress distribution curves are given. The far field condition is checked, The asymptotic behavior near the crack-tip is given. The project supported by National Natural Science Foundation of China