PZT-4紧凑拉伸试样的断裂分析

基于线性压电材料的复势理论,通过解析分析,导出了一种分析有限压电板裂纹问题的解析数值方法. 首先,计算了含中心裂纹有限板的断裂参数,与Woo和Wang的解析数值法(Int J Fract, 1993, 62: 203$\sim$218)相比较,表明该方法具有很高的精度和很好的计算效率. 随后,采用该方法和有限元法计算了PZT-4紧凑拉伸试样在绝缘裂纹面边界条件下断裂时的断裂参数,发现各断裂参数的临界值分散性很大,不能作为压电材料的单参数断裂准则. 进而,针对试样真实的裂隙形状,采用有限元法计算了裂隙尖端的应力、电位移场,比较了裂隙内介质的介电性能对裂隙尖端场的影响,计算了带微裂纹的真实裂隙模型的断裂参数并进行了理论分析.      

Applicability of the crack-face electromagnetic boundary conditions for fracture of magnetoelectroelastic materials

This paper discusses the different electromagnetic boundary conditions on the crack-faces in magnetoelectroelastic materials, which possess coupled piezoelectric, piezomagnetic and magnetoelectric effects. A notch of finite thickness in these materials containing air (or vacuum) is also addressed. Four ideal crack-face electromagnetic boundary condition assumptions, that is, (a) electrically and magnetically impermeable crack, (b) electrically impermeable and magnetically permeable crack, (c) electrically permeable and magnetically impermeable crack and (d) electrically and magnetically permeable crack, are investigated separately. The influence of notch thickness on the field intensity factors at notch tips and the electromagnetic field inside the notch are obtained in closed-form. The results are compared with the ideal crack solutions. Applicability of crack-face electromagnetic boundary condition assumptions is discussed.     

A dielectric crack in a functionally graded piezoelectric layer

Considering the dielectric effects inside a crack, the problem of an electrically dielectric crack in a functionally graded piezoelectric layer is addressed in this paper. The energetically consistent crack-face boundary conditions are utilized to analyze the effects of a dielectric of crack interior. Applying the Fourier transform technique, the boundary-value problem is reduced to solving three coupling singular equations. Then a system of non-linear algebraic equations is obtained and the field intensity factors along with the energy release rate are given. Numerical results show the differences of the electric displacement inside a crack, the stress and electric displacement intensity factors and the energy release rate using the permeable, impermeable, semi-permeable and energetically consistent boundary conditions respectively. The effects of the material non-homogeneity, the applied electric field and the discharge field of crack interior on the electrostatic traction acting on the crack faces and the energy release rate are further studied through the energetically consistent boundary conditions.     

2D dynamic analysis of cracks and interface cracks in piezoelectric composites using the SBFEM

The dynamic stress and electric displacement intensity factors of impermeable cracks in homogeneous piezoelectric materials and interface cracks in piezoelectric bimaterials are evaluated by extending the scaled boundary finite element method (SBFEM). In this method, a piezoelectric plate is divided into polygons. Each polygon is treated as a scaled boundary finite element subdomain. Only the boundaries of the subdomains need to be discretized with line elements. The dynamic properties of a subdomain are represented by the high order stiffness and mass matrices obtained from a continued fraction solution, which is able to represent the high frequency response with only 3–4 terms per wavelength. The semi-analytical solutions model singular stress and electric displacement fields in the vicinity of crack tips accurately and efficiently. The dynamic stress and electric displacement intensity factors are evaluated directly from the scaled boundary finite element solutions. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required. Numerical examples are presented to verify the proposed technique with the analytical solutions and the results from the literature. The present results highlight the accuracy, simplicity and efficiency of the proposed technique.     

Investigation of the crack problem in non-local piezoelectric materials under combined electromechanical loadings

The present investigation of the crack problem in piezoelectric materials is performed based on the non-local theory. After some manipulations, the impermeable crack, the permeable crack （the crack gap is full of NaCI solution）, and the semi-permeable crack （the crack gap is full of air or silicon oil） are reduced to a uniform formulation by assuming the normal electric displacement on the crack surfaces to be an unknown variable. Thus, a triple integral equation with the unknown normal electric displacement is established. By using the Newton iterative method and solving the triple integral equation, it is found that the normal electric displacement on the crack surfaces is no longer a constant as determined by previous studies, rather, it depends upon the remote combined electromechanical loadings. Numerical results of the stresses and electric displacement fields show that there are no singularities at the crack tips so that the stresses remain finite. It is of great significance that the concrete electric boundary condition on the crack surfaces exerts significant influence on the near-tip fields and in this way plays an important role in evaluating the crack stability in the non-local piezoelectric materials. More specifically, the impermeable crack model always overestimates the finite stresses at the crack tips, whereas the permeable crack model always underestimates them.     

Non-local theory solution for a Mode I crack in piezoelectric materials

In this paper, the non-local theory of elasticity is applied to obtain the behavior of a Griffith crack in the piezoelectric materials subjected to a uniform tension loading. The permittivity of the air in the crack is considered. By means of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations, in which the unknown variables are the jumps of the displacements across the crack surfaces. To solve the dual integral equations, the jumps of the displacements across the crack surfaces are expanded in a series of Jacobi polynomials. Numerical examples are provided to show the effects of the crack length, the materials constants, the electric boundary conditions and the lattice parameter on the stress and the electric displacement fields near the crack tips. It can be obtained that the effects of the electric boundary conditions on the electric displacement fields are large. Unlike the classical elasticity solutions, it is found that no stress and electric displacement singularities are present at the crack tips. The non-local elastic solutions yield a finite hoop stress at the crack tips, thus allowing us to use the maximum stress as a fracture criterion.     

Dielectric crack problem for a magnetoelectroelastic strip with functionally graded properties

The paper presents a fracture analysis for an electromagnetically dielectric crack in a functionally graded magnetoelectroelastic strip. It is considered that the material properties are varying exponentially along the width direction. Under the assumption of the in-plane magneto-electro-mechanical loadings, the dielectric crack is simulated by using the semi-permeable crack-face boundary conditions. The Fourier transform technique is applied to solve the boundary-value problem and four coupling singular integral equations are determined. A nonlinear system of algebraic equations is further derived and solved numerically to determine the electromagnetic field inside the crack. Then the field intensity factors of stress, electric displacement, and magnetic induction are given. Through the numerical computations, the effects of the material non-homogeneity and the permeability of crack interior on the electric displacement and the magnetic induction at the crack faces are studied. The variations of the intensity factors of stress, electric displacement, and magnetic induction versus the geometry of the crack, the strip width, and the material non-homogeneity are presented in graphics respectively.     

压电材料中裂纹面电边界条件的适用性

先分析理想裂纹面电边界条件(包括导通型裂纹和不导通型裂纹),然后分析缺陷的厚度及其对尖端场的影响,用封闭形式的解表示出了缺陷厚度的影响．结果指出,裂纹长厚比对计算结果是有影响的,除非裂纹极扁,不可通裂纹边界条件可以给出合理的结果,而可通边界条件只有在一些限制条件下才适用。     

Discussion on electromagnetic crack face boundary conditions for the fracture mechanics of magneto-electro-elastic materials

This paper discusses electromagnetic boundary conditions on crack faces in magneto- electroelastic materials, where piezoelectric, piezomagnetic and magnetoelectric effects are coupled. A notch of finite thickness in these materials is also addressed. Four idealized electromagnetic boundary conditions assumed for the crack-faces are separately investigated, i.e. (a) electrically and magnetically impermeable (crack-face), (b) electrically impermeable and magnetically permeable, (c) electrically permeable and magnetically impermeable, and (d) electrically and magnetically permeable. The influence of the notch thickness on important parameters, such as the field intensity factors, the energy release rate at the notch tips and the electromagnetic fields inside the notch, are studied and the results are obtained in closed-form. Results under different idealized electromagnetic boundary conditions on the crack-face are compared, and the applicability of these idealized assumptions is discussed.The project supported by the National Natural Science Foundation of China (10102004) The English text was polished by Yunming Chen.     

ANALYSIS OF 3-D NOTCHED/CRACKED STRUCTURES BY USING EXTENDED BOUNDARY ELEMENT METHOD 1)

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.     

三维切口/裂纹结构的扩展边界元法分析

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.      

改进型扩展比例边界有限元法

结合了扩展有限元法(extended finite elementmethods,XFEM)和比例边界有限元法(scaled boundary finite elementmethods,SBFEM)的主要优点,提出了一种改进型扩展比例边界有限元法(improvedextended scaled boundary finite elementmethods,$i$XSBFEM),为断裂问题模拟提供了一条新的途径.类似XFEM,采用两个正交的水平集函数表征材料内部裂纹面,并基于水平集函数判断单元切割类型;将被裂纹切割的单元作为SBFE的子域处理,采用SBFEM求解单元刚度矩阵,从而避免了XFEM中求解不连续单元刚度矩阵需要进一步进行单元子划分的缺陷;同时,借助XFEM的主要思想,将裂纹与单元边界交点的真实位移作为单元结点的附加自由度考虑,赋予了单元结点附加自由度明确的物理意义,可以直接根据位移求解结果得出裂纹与单元边界交点的位移;对于含有裂尖的单元,选取围绕裂尖单元一圈的若干层单元作为超级单元,并将此超级单元作为SBFE的一个子域求解刚度矩阵,超级单元内部的结点位移可通过SBFE的位移模式求解得到,应力强度因子可基于裂尖处的奇异位移(应力)直接获得,无需借助其他的数值方法.最后,通过若干数值算例验证了建议的$i$XSBFEM的有效性,相比于常规XFEM,$i$XSBFEM的基于位移范数的相对误差收敛性较好;采用$i$XSBFEM通过应力法和位移法直接计算得到的裂尖应力强度因子均与解析解吻合\较好.      

Study of a propagating finite crack in functionally graded piezoelectric materials considering dielectric medium effect

This paper provides a study of the problem of a propagating finite crack under in-plane loading in functionally graded piezoelectric materials (FGPMs). The analytical formulations are developed by Fourier transforms and the resulting singular integral equations are solved by using Chebyshev polynomials. By using a dielectric crack model with deformation-dependent electric boundary condition, numerical simulations are made to show the effects of the dielectric medium, the gradient of material properties and the speed of crack propagation on the fracture parameters, such as the stress, electric displacement and crack opening displacement intensity factors. A critical state for the electromechanical loading applied to the FGPMs is observed, which determines whether the traditionally impermeable (or permeable) crack model serves as the upper or lower bound for the dielectric model. The validity of this dielectric crack model is also examined by comparing the results of different existing crack models.     

Transient response of a magnetoelectroelastic solid with two collinear dielectric cracks under impacts

Dynamic analysis of two collinear electro-magnetically dielectric cracks in a piezoelectromagnetic material is made under in-plane magneto-electro-mechanical impacts. Generalized semi-permeable crack-face boundary conditions are proposed to simulate realistic opening cracks with dielectric. Ideal boundary conditions of a combination of electrically permeable or impermeable and magnetically permeable or impermeable assumptions are several limiting cases of the semi-permeable dielectric crack. Utilizing the Laplace and Fourier transforms, the mixed initial-boundary-value problem is reduced to solving singular integral equations with Cauchy kernel. Dynamic intensity factors of stress, electric displacement, magnetic induction and crack opening displacement (COD) near the inner and outer crack tips are determined in the Laplace transform domain. Numerical results for a special magnetoelectroelastic solid are calculated to show the influences of the dielectric permittivity and magnetic permeability inside the cracks on the crack-face electric displacement and magnetic induction. By means of a numerical inversion of the Laplace transform, the variations of the normalized intensity factors of stress and COD are discussed against applied magnetoelectric impact loadings and the geometry of the cracks for fully impermeable, vacuum, fully permeable cracks and shown in graphics.     

Transient thermal elastic fracture of a piezoelectric cylinder specimen

A finite piezoelectric cylinder with an embedded penny-shaped crack is investigated for a thermal shock load on the outer surface of the cylinder. The theory of linear electro-elasticity is applied to solve the transient temperature field and the associated thermal stresses and electrical displacements without crack. These thermal stresses and electrical displacements are added to the surfaces of the crack to form an electromechanical coupling and mixed mode boundary-value problem. The electrically permeable crack face boundary condition assumption is used, and the thermal stress intensity factor and electrical displacement intensity factor at the crack border are evaluated. The thermal shock resistance of the piezoelectric cylinder is evaluated for the analysis of piezoelectric material failure in practical engineering applications.     

Analysis method of planar cracks of arbitrary shape in the isotropic plane of a three-dimensional transversely isotropic magnetoelectroelastic medium

The hyper-singular boundary integral equation method of crack analysis in three-dimensional transversely isotropic magnetoelectroelastic media is proposed. Based on the fundamental solutions or Green’s functions of three-dimensional transversely isotropic magnetoelectroelastic media and the corresponding Somigliana identity, the boundary integral equations for a planar crack of arbitrary shape in the plane of isotropy are obtained in terms of the extended displacement discontinuities across crack faces. The extended displacement discontinuities include the displacement discontinuities, the electric potential discontinuity and the magnetic potential discontinuity, and correspondingly the extended tractions on crack face represent the conventional tractions, the electric displacement and the magnetic induction boundary values. The near crack tip fields and the intensity factors in terms of the extended displacement discontinuities are derived by boundary integral equation approach. A solution method is proposed by use of the analogy between the boundary integral equations of the magnetoelectroelastic media and the purely elastic materials. The influence of different electric and magnetic boundary conditions, i.e., electrically and magnetically impermeable and permeable conditions, electrically impermeable and magnetically permeable condition, and electrically permeable and magnetically impermeable condition, on the solutions is studied. The crack opening model is proposed to consider the real crack opening and the electric and magnetic fields in the crack cavity under combined mechanical-electric-magnetic loadings. An iteration approach is presented for the solution of the non-linear model. The exact solution is obtained for the case of uniformly applied loadings on the crack faces. Numerical results for a square crack under different electric and magnetic boundary conditions are displayed to demonstrate the proposed method.     

Exact solutions for the plane problem in piezoelectric materials with an elliptic or a crack

Based on the complex potential approach, the two-dimensional problems in a piezoelectric material containing an elliptic hole subjected to uniform remote loads are studied. The explicit, closed-form solutions satisfying the exact electric boundary condition on the hole surface are given both inside and outside the hole. When the elliptic hole degenerates into a crack, the field intensity factors are obtained. It is shown that the stress intensity factors are the same as that of isotropic material, while the electric displacement intensity factor depends on both the material properties and the mechanical loads, but not on the electric loads. In other words, the uniform electric loads have no influence on the field singularities. It is also shown that the impermeable crack assumption used previously to simply the electric condition is not valid to crack problems in piezoelectric materials.     

Theoretical model of crack branching in magnetoelectric thermoelastic materials

Thermomagnetoelectroelastic crack branching of magnetoelectro thermoelastic materials is theoretically investigated based on Stroh formalism and continuous distribution of dislocation approach. The crack face boundary condition is assumed to be fully thermally, electrically and magnetically impermeable. Explicit Green’s functions for the interaction of a crack and a thermomagnetoelectroelastic dislocation (i.e., a thermal dislocation, a mechanical dislocation, an electric dipole and a magnetic dipole located at a same point) are presented. The problem is reduced to two sets of coupled singular integral equations with the thermal dislocation and magnetoelectroelastic dislocation densities along the branched crack line as the unknown variables. As a result, the formulations for the stress, electric displacement and magnetic induction intensity factors and energy release rate at the branched crack tip are expressed in terms of the dislocation density functions and the branch angle. Numerical results are presented to study the effect of applied thermal flux, electric field and magnetic field on the crack propagation path by using the maximum energy release rate criterion.     

A mixed electric boundary value problem for a two-dimensional piezoelectric crack

In this paper, a mixed electric boundary value problem for a two-dimensional piezoelectric crack problem is presented, in the sense that the crack face is partly conducting and partly impermeable. By the analytical continuation method, the unknown electric charge distributions on the upper and lower conducting crack faces are reduced to two decoupled singular integral equations and then these two equations are converted into algebraic equations to find the full field solution. Though the results suggest that the stress intensity factors at the crack tip are identical to those of conventional piezoelectric materials, but the electric field and electric displacement are related to the electric boundary conditions on the crack faces. The electric field and electric displacement are singular not only at crack tips but also at the junctures between the impermeable part and conducting parts. Numerical results for the variations of the electric field, electric displacement field and J-integral with respect to the normalized impermeable crack length are shown. Some discussions for the energy release rate and the J-integral are made.     

Anisotropic thermopiezoelectric solids with an elliptic inclusion or a hole under uniform heat flow

The two-dimensional problem of a thermopiezoelectric material containing an elliptic inclusion or a hole subjected to a remote uniform heat flow is studied. Based on the extended Lekhnitskii formulation for thermopiezoelectricity, conformal mapping and Laurent series expansion, the explicit and closed-form solutions are obtained both inside and outside the inclusion (or hole). For a hole problem, the exact electric boundary conditions on the hole surface are used. The results show that the electroelastic fields inside the inclusion or the electric field inside the hole are linear functions of the coordinates. When the elliptic hole degenerates into a slit crack, the electroelastic fields and the intensity factors are obtained. The effect of the heat flow direction and the dielectric constant of air inside the crack on the thermal electroelastic fields are discussed. Comparison is made with two special cases of which the closed solutions exist and it is shown that our results are valid.