Crack tip field in piezoelectric/piezomagnetic media

This paper considers the magnetoelectroelastic problem of a crack in a medium possessing coupled piezoelectric, piezomagnetic and magnetoelectric effects. Based on the extended Stroh formalism, the general two-dimensional solutions to the magnetoelectroelastic problem are obtained, involving five analytic functions of different variables. The magnetoelectroelastic field around the crack tip is given. It contains five modes of square root singularities. Expressions of the stresses, electric displacements and magnetic inductions in the vicinity of the crack tip are derived and the field intensity factors are provided. The path-independent conservative integral is derived. The energy release rate is written in terms of those field intensity factors. The explicit algebraic results are given for a special case of an anti-plane crack in a magnetoelectroelastic medium.     

Dynamic anti-plane shear of a cracked piezoelectric ceramic

The dynamic theory of antiplane piezoelectricity is applied to solve the problem of a line crack subjected to horizontally polarized shear waves in an arbitrary direction. The problem is formulated by means of integral transforms and reduced to the solution of a Fredholm integral equation of the second kind. The path-independent integral G is extended here to include piezoelectric effects, and is evaluated at the crack tip to obtain the dynamic energy release rate. Numerical calculations are carried out for the dynamic stress intensity factor and energy release rate. The material is piezoelectric ceramic.     

反平面扩展裂纹的J积分与COD

给出一个以任意速率扩展的反平面裂纹与路径无关的Ｊ积分，证明Ｊ积分扩展裂纹尖端的张开位移（动态ＣＯＤ）之间有的简单的关系，Ｊ积分与能量释放率，动应力强度因子之间也有简单关系，利用这些关系，给出了动态ＣＯＤ与动应力强度因子之间的关系式。     

Dynamic energy release rate expression for a spreading circular fracture pattern in a plate

An expression for the dynamic energy release rate of a spreading circular fracture pattern in an uniformly stretched plate is derived from a path-independent energy integral. This is equivalent to the energy dissipated in an incremental growth of a constant velocity circular locus in which the radial tension is zero. The relationship to rapid fracture and multiple crack division becomes apparent as the arc distance separating adjacent crack tips becomes small relative to the radius of the circular locus. The condition at some small fixed distance behind the closely spaced crack tips is then well-stimulated by that of zero radial tension on a circular locus. Experiments could be carried out to establish the accuracy of the analytical dynamic energy release rate expression for multiple crack division.     

Energy release rate and path-independent integral in dynamic fracture of magneto-electro-thermo-elastic solids

The energy release rate and associated energy flux integral in dynamic fracture of magneto-electro-thermo-elastic solids are formulated with the inclusion of multi-field fully coupled effects based on fundamental principles of thermodynamics. The difference between the global and local dynamic contour integrals is caused by unsteady state, mechanical body force, electricity conduction and thermal effect as the closed contour including crack faces is chosen. This formulation successfully captures the crack-tip singularity of coupled fields, offers the right expression for the crack driving force, and resolves the controversial issue on magneto-electro-thermo-elastic fracture criterion. Especially, for steady-state crack propagation in a magneto-electro-elastic solid, the path-independent dynamic contour integral is determined from the asymptotic near-tip field solution based on the Stroh-type formalism and the resulting dynamic energy release rate has an odd dependence on the dynamic magnetic induction intensity factor and the dynamic electric displacement intensity factor.     

矩形压电体中反平面裂纹的电弹性场

基于线性压电理论,本文获得了含有中心反平面裂纹的矩形压电体中的奇异应力和电场。利用Fourier积分变换和Fourier正弦级数将电绝缘型裂纹问题化为对偶积分方程,并进一步归结为易于求解的第二类Fred-holm积分方程。获得了裂纹尖端应力、应变、电位移和电场的解析解,求得了裂纹尖端场的强度因子及能量释放率。分析了压电矩形体的几何尺寸对它们的影响。结果表明,对于电绝缘型裂纹,裂纹尖端附近的各个场变量都具有-1/2阶的奇异性,能量释放率与电荷载的方向及大小有关,并且有可能为负值。     

Arbitrarily oriented crack near interface in piezoelectric bimaterials

Arbitrarily oriented crack near interface in piezoelectric bimaterials is considered. After deriving the fundamental solution for an edge dislocation near the interface, the present problem can be expressed as a system of singular integral equations by modeling the crack as continuously distributed edge dislocations. In the paper, the dislocations are described by a density function defined on the crack line. By solving the singular integral equations numerically, the dislocation density function is determined. Then, the stress intensity factors (SIFs) and the electric displacement intensity factor (EDIF) at the crack tips are evaluated. Subsequently, the influences of the interface on crack tip SIFs, EDIF, and the mechanical strain energy release rate (MSERR) are investigated. The J-integral analysis in piezoelectric bimaterals is also performed. It is found that the path-independent of J₁-integral and the path-dependent of J₂-integral found in no-piezoelectric bimaterials are still valid in piezoelectric bimaterials.     

On Sanders' energy-release rate integral and conservation laws in finite elastostatics

Sanders showed in 1960, within the framework of two-dimensional elasticity, that in any body a certain integral I around a closed curve containing a crack is path-independent. I is equal to the rate of release of potential energy of the body with respect to crack length. Here we first derive, in a simple way, Sanders' integral I for a loaded elastic body undergoing finite deformations and containing an arbitrary void. The strain energy density need not be homogeneous nor isotropic and there may be body forces. In the absence of body forces, for flat continua, and for special forms of the strain energy density, it is shown that I reduces to the well-known vector and scalar path-independent integrals often denoted by J, L, and M.     

Energy release rate of small cracks in hyperelastic materials

The energy release rate of a small crack in an infinite hyperelastic medium, and subjected to large strain multiaxial loading conditions, is derived by considering the balance of configurational stresses acting on two planes: one cutting the center of the crack face, and the other at an infinite distance in front of the crack tip. The analysis establishes that the energy release rate of a small crack is always proportional to the size of the crack, irrespective of the loading conditions and the crack orientation. The balance of configurational stresses is illustrated for several benchmark cases including simple extension, pure shear and equibiaxial extension, and for perpendicular and inclined cracks.     

APPLICATION OF TREFFTZ BEM TO ANTI-PLANE PIEZOELECTRIC PROBLEM

Anti-plane electroelastic problems are studied by the Trefftz boundary elementmethod (BEM) in this paper. The Trefftz BEM is based on a weighted residual formulation andindirect boundary approach. In particular the point-collocation and Galerkin techniques, in whichthe basic unknowns are the retained expansion coefficients in the system of complete equations,are considered, Furthermore, special Trefftz functions and auxiliary functions which satisfy ex-actly the specified boundary conditions along the slit boundaries are also used to derive a specialpurpose element with local defects. The path-independent integral is evaluated at the tip of acrack to determine the energy release rate for a mode Ⅲ fracture problem. In final, the accuracyand efficiency of the Trefftz boundary element method are illustrated by an example and thecomparison is made with other methods.     

Anti-plane analysis on a finite crack in a one-dimensional hexagonal quasicrystal strip

The anti-plane fracture problem for a finite crack in a one-dimensional hexagonal quasicrystal strip is analyzed. By using Fourier transforms, the mixed boundary value problems are reduced to the dual integral equations. The solution of the dual integral equations is then expressed by the complete elliptic integrals of the first and the third kinds. The expressions for stress, strains, displacements and field intensity factors of the phonon and phason fields near the crack tip are obtained exactly. The path-independent integral derived by a conservation law equals the energy release rate, which can be used as a fracture criterion for a mode III fracture problem.     

Fracture of a rectangular piezoelectromagnetic body

The singular stress, electric fields and magnetic fields in a rectangular piezoelectromagnetic body containing a center Griffith crack under longitudinal shear are obtained. Fourier transforms and Fourier sine series are used to reduce the mixed boundary value problems of the crack, which is assumed to be impermeable, to dual integral equations. The solution of the dual integral equations is then expressed in terms of Fredholm integral equations of the second kind. Expressions for stresses, electric displacements and magnetic inductions in the vicinity of the crack tip are derived. Also obtained are the field intensity factors and the energy release rates. Numerical results obtained show that the geometry of the rectangular body have significant influence on the field intensity factors and the energy release rates.     

Transient response of a rectangular piezoelectric medium with a center crack

The dynamic field intensity factors and energy release rates in a rectangular piezoelectric ceramic medium containing a center crack are obtained for boundary conditions of a permeable and an impermeable crack under electro-mechanical impact loading. An integral transform method is used to reduce the problem to two pairs of dual integral equations, which are then expressed as Fredholm integral equations of the second kind. Numerical values on the dynamic energy release rate are obtained to show the dependences upon the geometry and electric field.     

ELECTROELASTIC FIELD FOR AN IMPERMEABLE ANTI-PLANE SHEAR CRACK IN A PIEZOELECTRIC CERAMICS PLATE

ＩｎｔｒｏｄｕｃｔｉｏｎＤｕｅｔｏｔｈｅｉｎｔｒｉｎｓｉｃｃｏｕｐｌｉｎｇｃｈａｒａｃｔｅｒｉｓｔｉｃｓｂｅｔｗｅｅｎｅｌｅｃｔｒｉｃａｎｄｅｌａｓｔｉｃｂｅｈａｖｉｏｒｓ,ｔｈａｔｉｓ,ａｐｐｌｉｅｄｍｅｃｈａｎｉｃａｌｌｏａｄｉｎｇｐｒｏｄｕｃｅｓｅｌａｓｔｉｃｄｅｆｏｒｍａｔｉｏｎ ,ａｓｗｅｌｌａｓｅｌｅｃｔｒｉｃｆｉｅｌｄ ,ａｎｄｃｏｎｖｅｒｓｅｌｙｅｌｅｃｔｒｉｃｆｉｅｌｄｃａｎｇｉｖｅｒｉｓｅｔｏｅｌａｓｔｉｃｄｅｆｏｒｍａｔｉｏｎ ,ｐｉｅｚｏｅｌｅｃｔｒｉｃｍａｔｅｒｉａｌｓｈａｖｅ…     

Investigation on thermomechanical fracture in the framework of configurational forces

A configurational force approach is developed for providing a fresh look onto classical aspects of thermomechanical fracture. The theoretical framework is based on the finite deformation and makes no restrictions on the material response. The integral form of configurational force balance at the crack tip is constructed, and the concentrated configurational body force is decomposed into the inertial and internal parts. The energy release rate is evaluated through the generalized second law of thermodynamics applicable to configurational force system. The theoretical investigation shows that the negative of the projection of the internal configurational force concentrated at the crack tip along the direction of crack propagation plays the role of the energy release rate and acts directly in response to crack propagation. This finding enables us to deal with the thermomechanical fracture problems in material space.     

Angle cracks in two bonded functionally graded piezoelectric material under anti-plane shear

Consider two bonded functionally graded piezoelectric material (FGPM) with finite height. Each material contains an arbitrary oriented crack. The material properties are assumed in exponential forms in the direction normal to the interface. The crack surface condition is assumed to be electrically impermeable or permeable. Using the Fourier transform technique, the problem can be reduced to a system of singular integral equations, which are then solved numerically by applying the Gauss-Chebyshev integration formula to obtain the stress intensity factors at the crack tips. Numerical calculations are carried out to obtain the energy density factor S and the energy release rate G. In impermeable case, the energy release rate has been shown to be negative as the electric loads are applied. The positive definite characteristic of the energy density factor makes it possible for predicting the fracture behavior of the cracked structure. The influences of the non-homogeneous parameters and crack orientation on the energy density factors at the crack tips are discussed in detail. The results show that the energy density factor at the crack tip will be increased when the crack tip is located within the softer material.     

Interaction of three interfacial Griffith cracks between bonded dissimilar orthotropic half planes

The paper deals with the interaction of a pair of outer cracks on a central crack situated at the interface of two dissimilar orthotropic half-planes. The mixed boundary value problem is reduced to solving a pair of simultaneous singular integral equations which have finally been solved numerically by using Jacobi polynomials. The analytical expressions for stress intensity factors at the central crack tip and the expression of the strain energy release rate have been derived for general loading. Numerical values of the interaction effects of the outer cracks on the central crack have been calculated through stress magnification factors. It is seen that the interaction effects are either shielding or amplification depending on the size of the outer cracks and their spacing from the central crack.     

Anti-plane stress intensity, energy release and energy density at crack tips in a functionally graded strip with linearly varying properties

The fracture problem of a crack in a functionally graded strip with its properties varying in a linear form along the strip thickness under an anti-plane load is considered. The embedded anti-plane crack is located in the middle of strip half way through the thickness. The third mode stress intensity factor is derived using two different methods. In the first method, by employing Fourier integral transforms, the governing equation is converted to a singular integral equation, which is subsequently solved numerically by the collocation method based on Chebyshev polynomials. Then, the problem is solved by means of finite element method in which quadrilateral 8-node singular elements around each crack tip are used. After inspecting the validity of the solution technique, effects of crack geometry and non-homogeneous material parameter on the stress intensity, energy release and energy density are studied and the results of analytical and FEM solutions are compared.     

陶瓷基复合材料基体唇形裂纹失效规律研究

建立横向拉伸载荷下的唇形裂纹数学模型,采用复变函数的方法,通过保角映射,推导了唇形裂纹尖端应力场和位移场的解析解,建立了唇形裂纹的应力强度因子准则和最大能量释放率准则,结合算例分析陶瓷基复合材料基体唇形裂纹的几何参数、外载荷和纤维分布对失效准则的影响规律.结果 表明,(1)裂纹尖端应力场和位移场的解析解与有限元计算结果...     

A contour integral method for dynamic stress intensity factors

The contour integral method previously used to determine static stress intensity factors is applied to dynamic crack problems. The required derivatives of the traction in the reference problem are obtained numerically by the displacement discontinuity method. Stress intensity factors are determined by an integral around a contour which contains a crack tip. If the contour is chosen as the outer boundary of the body, the stress intensity factor is obtained from the boundary values of traction and displacement. The advantage of this path-independent integral is that it yields directly both the opening-mode and sliding-mode stress intensity factors for a straight crack. For dynamic problems, Laplace transforms are used and the dynamic stress intensity factors in the time domain are determined by Durbin's inversion method. An indirect boundary element method, incorporating both displacement discontinuity and fictitious load techniques, is used to determine the boundary or contour values of traction and displacement numerically.