Analysis of a crack in a half-plane piezoelectric solid with traction-induction free boundary

In this paper, the problem of a crack embedded in a half-plane piezoelectric solid with traction-induction free boundary is analyzed. A system of singular integral equations is formulated for the materials with general anisotropic piezoelectric properties and for the crack with arbitrary orientation. The kernel functions developed are in complex form for general anisotropic piezoelectric materials and are then specialized to the case of transversely isotropic piezoelectric materials which are in real form. The obtained coupled mechanical and electric real kernel functions may be reduced to those kernel functions for purely elastic problems when the electric effects disappear. The system of singular integral equations is solved numerically and the coupling effects of the mechanical and electric phenomena are presented by the generalized stress intensity factors for transversely isotropic piezoelectric materials.     

Crack energy density in arbitrary direction for piezoelectrics and its characteristic features in electromechanical mixed mode crack

The concepts of crack energy density (CED) and its derivatives in arbitrary direction were established for piezoelectric material and, keeping their application to mixed mode fracture in mind, the characteristic features of them as fracture parameters were investigated based on the approximate equations for CED and its derivatives. That is, CED and its derivatives in arbitrary direction are defined first and separation into their each mode contribution is made. Subsequently, path independent integral expressions of them are derived, and then using them, approximate equations of each mode contribution of CED are obtained concretely for the case where linear singular solution is known. The resulting equations are then used to investigate the effects of electric field and electrical boundary condition on CED and its derivatives. An infinite piezoelectric plane with a crack inclined with respect to the poling direction is considered as a numerical example. Mode I contribution of mechanical CED is mainly employed as a possible fracture parameter for the study and it was shown that applied electric field significantly influences on fracture parameters especially for the impermeable crack perpendicular to the poling direction. The effect of electric field has the tendency to decrease as crack inclination angle increases. It was also found that, even for the impermeable crack perpendicular to the poling direction, crack propagation could be deviated from self-similar direction under a strong negative electric field, and this fact is qualitatively consistent with an existing experimental observation. For the ideally sharp crack with no width, impermeable and Hao and Shen type boundary conditions are admissible showing qualitative agreement with experimental results, but exact boundary condition is not suitable and finally consistent with permeable boundary condition.     

Interface crack problems for mode Ⅱ of double dissimilar orthotropic composite materials

The fracture problems near the interface crack tip for mode Ⅱ of double dissimilar orthotropic composite materials are studied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized bi-harmonic equations,the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions,a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about himaterial engineering parameters. According to the uniqueness theorem of limit,both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same,the stress singularity exponents,stress intensity factors and stresses for mode Ⅱ crack of the orthotropic single material are obtained.     

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)     

Thermal response of a partially insulated interface crack in a graded coating-substrate structure under thermo-mechanical disturbance: Energy release and density

Consider the thermal fracture problem of a functionally graded coating-substrate structure of finite thickness with a partially insulated interface crack subjected to thermal-mechanical supply. A new model is proposed that the heat conduction through the crack region occurs and the temperature drop across the crack surfaces is the result of the thermal resistance. For the first time, real fundamental solutions are derived for the fracture analysis of functionally graded materials. The complicated mixed boundary problems of equations of heat conduction and elasticity are converted analytically into singular integral equations, which are solved numerically. The asymptotic expressions with higher order terms for the singular integral kernels are considered to improve the accuracy and efficiency of the numerical integration. Explicit expressions of various failure modes including stress intensity factors, energy release rate and strain energy density, are provided. Numerical results are presented to illustrate the effects of non-homogeneity parameters and the dimensionless thermal resistance on the temperature distribution along the crack surfaces and extended crack line, the thermal stress intensity factors and minimum strain energy density.     

Interface crack problems for mode II of double dissimilar orthotropic composite materials

The fracture problems near the similar orthotropic composite materials are interface crack tip for mode Ⅱ of double disstudied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized hi-harmonic equations, the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions, a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about bimaterial engineering parameters. According to the uniqueness theorem of limit, both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same, the stress singularity exponents, stress intensity factors and stresses for mode II crack of the orthotropic single material are obtained.     

压电材料中椭圆片状裂纹问题精确解的一种方法

在线性压电陶瓷本构关系和裂纹边界绝缘的框架下,用超奇异积分方程的方法对椭圆类片状裂纹问题进行了重新研究.超奇异积分方程中的未知位移间断和电势间断近似地表示为基本密度函数与多项式之积,其中基本密度函数反映了椭圆片状裂纹前沿电弹性场的奇异性,而多项式在均布载荷作用下可用一个常数来表达.引入椭球坐标系后,得到了均布载荷作用下未知位移间断和电势间断的解析解.使用这些解析解和电弹性场强度的定义,得到了裂纹前沿Ⅰ型、Ⅱ型和Ⅲ型应力强度因子以及电位移强度因子的精确表达式.法向均布载荷作用下的结果与现有精确解完全一致,切向均布载荷作用下的结果则尚未见有其它报道.     

Crack energy density of a piezoelectric material under general electromechanical loading

Crack energy density is considered and used as a possible fracture parameter in piezoelectricity under arbitrary electromechanical remote loads. The closed-form solution of a crack in a piezoelectric infinite plate subjected to general static electromechanical loading is obtained through a method alternative to the more common Stroh’s formalism. This analytical method, which is based on the spectral theorem of linear algebra, involves a transformation of similarity induced by the fundamental matrix in order to express the equations governing the problem in terms of complex potentials. The application of the mechanical boundary condition of stress-free crack and of one of the three considered electric boundary conditions (impermeable, permeable or semipermeable) leads then to the formulation of a Hilbert problem whose solution yields the stress and displacement fields. The crack energy density factors for mixed mode are then calculated under different mechanical and electrical loadings, as well as under different electric boundary conditions. The non-singular terms of the stress expressions are retained as well. The definition of the minimum energy density fracture criterion, as proposed by Sih, is given, and the influence of load biaxiality and positive or negative applied electric field on the criterion results is analyzed. The prediction of the incipient branching angle as from the energy density approach is also compared to that arising from the maximum circumferential stress theory for a mixed mode loading condition. Numerical results and graphs are presented and discussed for a PZT-4 piezoelectric ceramic.     

基于Erdogan基本解边界元法计算应力强度因子

引入含裂纹问题基本解(Erdogan基本解),提出了基于Erdogan基本解的样条虚边界 元法,并阐述了该法在实施过程中的特点与具体做法. 采用该方法详细分析了若干 典型裂纹问题,全面考察了方法的计算精度和收敛情况,以及在求解复杂裂纹问题方面 的能力. 结果显示,该方法具有精度高、收敛快、计算能力强等优点,是裂纹问题分析中 一种具有竞争力的通用计算方法.     

功能梯度压电材料反平面裂纹问题

基于三维弹性理论和压电理论，导出了材料系数在横观各向同性平面内梯度分布的压电体的状态方程，进而对材料系数指数函数规律分布的半无限大压电体中的反平面裂纹问题进行了求解，利用Fourier变换给出了半无限大压电体中位移，应力，电势及电位移的解析表达式，并求得了裂纹尖端的应力强度因子和电位移强度因子，分析了不同的非均匀材料系数及几何尺寸对它们的影响。     

Electroelastic analysis of a penny-shaped crack in a piezoelectric ceramic under mode I loading

Following the theory of linear piezoelectricity, we consider the electroelastic problem for a piezoelectric ceramic with a penny-shaped crack under mode I loading. The problem is formulated by means of Hankel transform and the solution is solved exactly. The stress intensity factor, energy release rate and energy density factor for the exact and impermeable crack models are expressed in closed form and compared for a P-7 piezoelectric ceramic. Based on current findings, we suggest that the energy release rate and energy density factor criteria for the exact crack model are superior to fracture criteria for the impermeable crack model.     

The scattering of harmonic elastic anti-plane shear waves by two collinear cracks in the piezoelectric plate by using the non-local theory

In this paper, the dynamic interaction between two collinear cracks in a piezoelectric material plate under anti-plane shear waves is investigated by using the non-local theory for impermeable crack surface conditions. By using the Fourier transform, the problem can be solved with the help of two pairs of triple integral equations. These equations are solved using the Schmidt method. This method is more reasonable and more appropriate. Unlike the classical elasticity solution, it is found that no stress and electric displacement 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 project supported by the Natural Science Foundation of Heilongjiang Province and the National Natural Science Foundation of China(10172030, 50232030)     

Penny-shaped crack in a piezoelectric cylinder under electro-mechanical loads

Summary The fracture behavior of a penny-shaped crack in a axisymmetrical piezoelectric ceramic cylinder of finite radius under mechanical and electrical loads is analyzed under electric continuous boundary conditions on the crack surface. The potential theory and Hankel transform are used to obtain a system of dual integral equations, which is then expressed as a Fredholm integral equation. Singular mechanical and electrical fields and field intensity factors of mode I are obtained. The numerical values of various field intensity factors for PZT-6B piezoelectric ceramics are graphically shown for a uniform load and a ring-shaped load, respectively. The effects of the radius of the cylinder on the field intensity factors are investigated. This work was supported by the Korea Research Foundation Grant (KRF-2001-041-E00057).     

The stress intensity factor history for an advancing crack in a transversely isotropic solid under 3-D loading

The mode I extension of a half plane crack in a transversely isotropic solid under 3-D loading is analyzed. Firstly, the fundamental problem that the crack is subjected to a pair of unit point loads on its faces is considered. Transform methods are used to reduce the boundary value problem to a single integral equation that can be solved by the Wiener–Hopf technique. The Cagniard–de Hoop method is employed to invert the transforms. An exact expression is derived for the mode I stress intensity factor as a function of time and position along the crack edge. Based on the fundamental solution, the stress intensity factor history due to general loading is then obtained. Some features of the solutions are discussed through numerical results.     

The fundamental solutions for the plane problem in piezoelectric media with an elliptic hole or a crack

1.IntroductionItiswell-knownthatthefundame,ltalsolutionsorGreen'sfunctionsplayanimportantroleilllinearelasticity.Forexample,theycanbeusedtoconstructmanyanalyticalsolutionsofpracticalproblems.Itismoreimportantthattheyareusedasthefundamentalsolutionsintheboundaryelementmethod(BEM)tosolvesomecomplicatedproblem.Withthewidely-increasingapplicationofpiezoelectricmaterialsinengineeringproblems,thestudyregardingtheGreen'sfLlnctionsinpiezoelectricsolidshasreceivedmuchinterest.The3DGreen'sfunctionsi…     

Elastic field in an infinite medium of one-dimensional hexagonal quasicrystal with a planar crack

This present work is concerned with planar cracks embedded in an infinite space of one-dimensional hexagonal quasicrystals. The potential theory method together with the general solutions is used to develop the framework of solving the crack problems in question. The mode I problems of three common planar cracks (a penny-shaped crack, an external circular crack and a half-infinite crack) are solved in a systematic manner. The phonon and phason elastic fundamental fields along with some important parameters in crack analysis are explicitly presented in terms of elementary functions. Several examples are given to show the applications of the present fundamental solutions. The validity of the present solutions is discussed both analytically and numerically. The derived analytical solutions of crack will not only play an important role in understanding the phonon–phason coupling behavior in quasicrystals, but also serve as benchmarks for future numerical studies and simplified analyses.     

Analysis on non-oscillatory singularity behaviors of mode Ⅱ interface crack tip in orthotropic bimaterial

The fracture behaviors near the mode Ⅱ interface crack tip for orthotropic bimaterial are studied. The non-oscillatory field, where the stress singularity exponent is a real number, is discussed by the complex function method and the undetermined coefficient method. From the research fracture problems, the stress functions with ten undetermined coefficients and an unknown singularity exponent are introduced when?_1 0 and ?_2 0. By the existence theorem of non-trival solutions for the system of eight homogeneous linear equations, the characteristic equation, the stress singularity exponent, and the discriminating condition of the non-oscillatory singularity are found.By the uniqueness theorem of the solutions for the system of twelve non-homogeneous linear equations with ten unknowns, the ten undermined coefficients in the stress functions are uniquely determined. The definitions of the stress intensity factors are given with the help of one-sided limit, and their theoretical formulae are deduced. The analytic solutions of the stresses near the mode Ⅱ interface crack tip are derived. The classical results for orthotropic material are obtained.     

平行于异材界面的三维椭圆平片裂纹分析

讨论了拉伸载荷作用下平行于两相材料界面的椭圆平片裂纹问题．首先,使用有限部积分概念和两相材料界面完全接合时的点力基本解导出了一组以裂纹表面位移差为未知函数的超奇异积分方程组．该组方程表明,此时三种裂纹模型同时存在;其次,在数值求解该组方程的过程中,未知函数裂纹表面位移差被近似为位移差的基本密度函数与多项式之积．基本密度函数反映了裂纹前沿应力奇性性态;最后,以拉伸载荷为例,讨论了椭圆平片裂纹与界面的距离、裂纹形状比和不同材料组合对应力强度因子的影响,并以图表形式给出。     

Analysis to the stress and strain structures of mode I 3-D crack in elastic-plastic state

In this paper, the stress and strain structures of Mode I 3-D crack in power hardening material are studied by analyzing the fundamental equations of elastic-plastic mechanics. It is shown that three regions, Z1,Z2 and Z3 can be divided in the thickness direction according to the stress characteristic. In region Z1, the stress components in the plane Perpendicular to z axis (thickness direction) can be solved first using the fundamental equations of plane strain state; in region Z3, they can be solved first by the equations of plane stress state. The region Z2 is defined as a transition layer. It is shown that the transition layer is the characteristic of Mode I 3-D crack in elastic-plastic state, and it is significant to the research on 3-D fracture. The crack tip opening displacement CTOD is chosen to describe the amplitude coefficient of the local stress field, and the distribution of CTOD in 3-D state is investigated.The project supported by National Natural Science Foundation of China.     

Penny shaped crack in a piezoelectric cylinder surrounded by an elastic medium subjected to combined in-plane mechanical and electrical loads

The fracture problem of a penny shaped crack in a piezoelectric ceramic cylinder surrounded by an infinite elastic medium under in-plane normal mechanical and electrical loads is considered with the electric continuous boundary conditions on the crack surface. By using the potential theory and Hankel transform, a system of dual integral equations is obtained, and expressed to a Fredholm integral equation of the second kind. The mechanical and electrical field equations and all sorts of field intensity factors of mode I are obtained, and the numerical values of various field intensity factors for PZT-6B piezoelectric ceramic surrounded by several different elastic media are graphically shown for a uniform load and a ring-shaped load, respectively. And the effects of the size of the piezoelectric cylinder and the elastic material properties on various field intensity factors are obtained.