A Novel Time‐Domain BIEM for Wave Propagation Analysis in Anisotropic Solids With Cracks

Analysis of elastic wave propagation in anisotropic solids with cracks is of particular interest to quantitative non‐destructive testing and fracture mechanics. For this purpose, a novel time‐domain boundary integral equation method (BIEM) is presented in this paper. A finite crack in an unbounded elastic solid of general anisotropy subjected to transient elastic wave loading is considered. Two‐dimensional plane strain or plane stress condition is assumed. The initial‐boundary value problem is formulated as a set of hypersingular time‐domain traction boundary integral equations (BIEs) with the crack‐opening‐displacements (CODs) as unknown quantities. A time‐stepping scheme is developed for solving the hypersingular time‐domain BIEs. The scheme uses the convolution quadrature formula of Lubich [1] for temporal convolution and a Galerkin method for spatial discretization of the BIEs. An important feature of the present time‐domain BIEM is that it uses the Laplace‐domain instead of the more complicated time‐domain Green's functions. Fourier integral representations of Laplace‐domain Green's functions are applied. No special technique is needed in the present time‐domain BIEM for evaluating hypersingular integrals.     

Transient response of two collinear dielectric cracks in a piezoelectric solid under inplane impacts

The paper is focused on the dynamic analysis of two collinear dielectric cracks in a piezoelectric material under the action of in-plane electromechanical impacts. Considering the dielectric permeability of crack interior, the electric displacements at the crack surfaces are governed by the jumps of electric potential and crack opening displacement across the cracks. The permeable and impermeable crack models are the limiting cases of the general one. The Laplace and Fourier transform techniques are further utilized to solve the mixed initial-boundary-value problem, and then to obtain the singular integral equations with Cauchy kernel, which are solved numerically. Dynamic intensity factors of stress, electric displacement and crack opening displacement are determined in time domain by means of a numerical inversion of the Laplace transform. Numerical results for PZT-5H are calculated to show the effects of the dielectric permeability inside the cracks, applied electric loadings and the geometry of the cracks on the fracture parameters in graphics. The observations reveal that based on the COD intensity factor, a positive electric field enhances the dynamic dielectric crack growth and a negative one impedes the dynamic dielectric crack growth in a piezoelectric solid.     

Transient response of dissimilar piezoelectric layers with multiple interacting interface cracks

In this study, we examine the dynamic behavior of two bonded dissimilar piezoelectric layers containing multiple interfacial cracks subjected to electro-mechanical impact loading. The problem was formulated through Fourier transformation into singular integral equations in which the unknown variables are the jumps of displacement and electric potential across the crack surface in the Laplace transform domain. The resulting integral equations together with the corresponding single-valued conditions are solved numerically for the densities of electro-elastic dislocations on a crack surface. The dynamic field intensity factors and dynamic energy release rate (DERR) history are obtained for both permeable and impermeable crack. The stress field is also obtained for the interface crack under impact loads. The results show that the field intensity factors at the crack tips and dynamic energy release rate depend on the interfacial crack geometry, electromechanical coupling and the electric boundary conditions on the crack surface.     

压电压磁复合材料中界面裂纹对弹性波的散射

利用Schmidt方法分析了压电压磁复合材料中可导通界面裂纹对反平面简谐波的散射问题．经过富里叶变换得到了以裂纹面上的间断位移为未知变量的对偶积分方程A·D2在求解对偶积分方程的过程中,裂纹面上的间断位移被展开成雅可比多项式的形式．数值模拟分析了裂纹长度、波速和入射波频率对应力强度因子、电位移强度因子、磁通量强度因子的影响A·D2从结果中可以看出,压电压磁复合材料中可导通界面裂纹的反平面问题的应力奇异性形式与一般弹性材料中的反平面问题应力奇异性形式相同．     

Transient thermoelastic response in a cracked strip of functionally graded materials via generalized fractional heat conduction

This work is devoted to analyzing a thermal shock problem of an elastic strip made of functionally graded materials containing a crack parallel to the free surface based on a generalized fractional heat conduction theory. The embedded crack is assumed to be insulated. The Fourier transform and the Laplace transform are employed to solve a mixed initial-boundary value problem associated with a time-fractional partial differential equation. Temperature and thermal stresses in the Laplace transform domain are evaluated by solving a system of singular integral equations. Numerical results of the thermoelastic fields in the time domain are given by applying a numerical inversion of the Laplace transform. The temperature jump between the upper and lower crack faces and the thermal stress intensity factors at the crack tips are illustrated graphically, and phase lags of heat flux, fractional orders, and gradient index play different roles in controlling heat transfer process. A comparison of the temperature jump and thermal stress intensity factors between the non-Fourier model and the classical Fourier model is made. Numerical results show that wave-like behavior and memory effects are two significant features of the fractional Cattaneo heat conduction, which does not occur for the classical Fourier heat conduction.     

压电材料中两个非对称平行裂纹的基本解

采用Schmidt方法分析压电材料中非对称平行的双可导通裂纹的断裂性能．利用Fourier变换使问题的求解转换为求解两对以裂纹面位移之差为未知变量的对偶积分方程．为了求解对偶积分方程,直接把裂纹面位移差函数展开成Jacobi多项式形式．最终得到了裂纹的应力强度因子与电位移强度因子之间的关系．数值结果表明,应力强度因子和电位移强度因子与裂纹间的距离、裂纹的几何尺寸有关;与不可导通裂纹有关结果相比,可导通裂纹的电位移强度因子远小于相应问题不可导通裂纹的电位移强度因子．同时可以发现裂纹间的“屏蔽”效应也在压电材料中出现．     

利用Schmidt方法研究压电材料Ⅰ-型界面裂纹问题

在一定的假设条件下,即不考虑界面裂纹尖端处裂纹面的相互叠入现象,研究了压电材料Ⅰ-型界面裂纹问题．利用Fourier变换使问题的求解转换为求解两对对偶积分方程．进而把裂纹表面位移差展开成Jacobi多项式形式来求解对偶积分方程．结果表明裂纹尖端应力场和电位移场的奇异性与均匀材料裂纹问题的奇异性相同．当上下半平面材料相同时,解可以退化而得到其精确解．     

横观各向同性电磁弹性介质中裂纹对SH波的散射

研究横观各向同性电磁弹性介质中裂纹和反平面剪切波之间的相互作用．根据电磁弹性介质的平衡运动微分方程、电位移和磁感应强度微分方程,得到SH波传播的控制场方程．引入线性变换,将控制场方程简化为Helmholtz方程和两个Laplace方程A·D2通过Fourier变换,并采用非电磁渗透型裂面边界条件,得到了柯西奇异积分方程组．利用Chebyshev多项式求解积分方程,得到应力场、电场和磁场以及动应力强度因子的表达,并给出了数值算例．     

Non-local theory solution of two collinear mode-I cracks in piezoelectric materials

In this paper, the basic solution of two collinear cracks in a piezoelectric material plane subjected to a uniform tension loading is investigated by means of the non-local theory. Through the Fourier transform, the problem is solved with the help of two pairs of integral equations, in which the unknown variables are the jumps of displacements across the crack surfaces. To solve the integral equations, the jumps of displacements across the crack surfaces are directly expanded in a series of Jacobi polynomials. Numerical examples are provided to show the effects of the interaction of two cracks, the materials constants and the lattice parameter on the stress field and the electric displacement field near crack tips. Unlike the classical elasticity solution, it is found that no stress and electric displacement singularities are present at crack tips. The non-local elastic solutions yield a finite hoop stress at the crack tip, thus allowing us to using the maximum stress as a fracture criterion in piezoelectric materials.     

An interface crack moving between magnetoelectroelastic and functionally graded elastic layers

A constant crack moving along the interface of magnetoelectroelastic and functionally graded elastic layers under anti-plane shear and in-plane electric and magnetic loading is investigated by the integral transform method. Fourier transforms are applied to reduce the mixed boundary value problem of the crack to dual integral equations, which are expressed in terms of Fredholm integral equations of the second kind. The singular stress, electric displacement and magnetic induction near the crack tip are obtained asymptotically and the corresponding field intensity factors are defined. Numerical results show that the stress intensity factors are influenced by the crack moving velocity, the material properties, the functionally graded parameter and the geometric size ratios. The propagation of the moving crack may bring about crack kinking, depending on the crack moving velocity and the material properties across the interface.     

压电材料中两平行不相等界面裂纹的动态特性研究

利用Schmidt方法,研究了压电材料中两个平行不相等的可导通界面裂纹对简谐反平面剪切波的散射问题．利用Fourier变换,使问题的求解转换为对两对以裂纹面张开位移为未知变量的对偶积分方程的求解．数值计算结果表明,动态应力强度因子及电位移强度因子受裂纹的几何参数、入射波频率的影响．在特殊情况下,与已有结果进行了比较分析．同时,电位移强度因子远小于不可导通电边界条件下相应问题的结果．     

压电压磁复合材料中一对平行裂纹对弹性波的散射

利用Schmidt方法对压电压磁复合材料中一对平行对称裂纹对反平面简谐波的散射问题进行了分析,借助富里叶变换得到了以裂纹面上的间断位移为未知变量的对偶积分方程．在求解对偶积分方程的过程中,裂纹面上的间断位移被展开成雅可比多项式的形式,最终获得了应力强度因子、电位移强度因子、磁通量强度因子三者之间的关系．结果表明,压电压磁复合材料中平行裂纹动态反平面断裂问题的应力奇异性与一般弹性材料中的动态反平面断裂问题的应力奇异性相同,同时讨论了裂纹间的屏蔽效应．     

Mathematical model for crack arrest of a transversely cracked piezoelectromagnetic strip – Part I

A magnetic, electric and mechanical yield model is proposed for a cracked piezoelectromagnetic ceramic narrow strip. The strip is subjected to anti-plane mechanical and in-plane electric and magnetic loads, consequently the crack opens in self-similar fashion forming a magnetic, a saturation and a slide zone ahead each tip. These in turn are arrested by prescribing a magnetic, electric and mechanical load, respectively. Employing Fourier integral transform the problem reduces to the solution of three dual integral equations. The solution of dual integral equations is then expressed in terms of Fredholm integral equation of second kind. Expressions are derived for yield induction zone, slide-yield zone and saturation zone lengths, energy release rate. A case study is carried for BaTiO₃–CoFe₂O₄ and results are presented graphically. It is shown that proposed model is capable of crack opening arrest under small-scale-yielding.     

Boundary integral equation formulation for generalized thermoelastic diffusion – Analytical aspects

The present work deals with the formulation of the boundary integral equations for the solution of equations under linear theory of generalized thermoelastic diffusion in a three-dimensional Euclidean space. A mixed initial-boundary value problem is considered in the present context and the fundamental solutions of the corresponding coupled differential equations are obtained in the Laplace transform domain by employing the treatment of scalar and vector potential theory. A reciprocal relation of Betti type is established. Then we formulate the boundary integral equations for generalized thermoelastic diffusion on the basis of these fundamental solutions and the reciprocal relation.     

不同功能梯度压电压磁层状介质中共线界面裂纹的动态性能分挝

对不同功能梯度压电压磁层状介质中,共线界面裂纹对简谐应力波作用下的动态问题,进行了分析.经Fourier变换,使问题的求解转换为求解以裂纹面上位移间断为未知量的三重对偶积分方程,三重对偶方程可以采用Schmidt方法来求解,进而分析了功能梯度参数、入射波频率和层状介质厚度对应力、电位移和磁通量强度因子的影响.     

平面十次对称准晶中Ⅱ型Griffith裂纹的求解

应用应力函数法,求解了二维十次对称准晶中的Ⅱ型Griffith裂纹问题。特点是把二维准晶的弹性力学问题分解成一个平面应变问题与一个反平面问题的叠加,通过引入应力函数,把平面应变问题的十八个弹性力学基本方程简化成一个八阶偏微分方程,并且求出了其在Ⅱ型Griffith裂纹情况的混合边值问题的解,所有的应力分量和位移分量都用初等函数表示出来,并且由此得出了准晶中Ⅱ型Griffith裂纹问题的应力强度因子和能量释放率。     

On the convergence of a numerical method for a hypersingular integral equation on a closed surface

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises in the solution of the Neumann boundary value problem for the Laplace equation with a representation of a solution in the form of a double-layer potential. We consider the case in which the interior or exterior boundary value problem is solved in a domain; whose boundary is a smooth closed surface, and an integral equation is written out on that surface. For the integral operator in that equation, we suggest quadrature formulas like the method of vortical frames with a regularization, which provides its approximation on the entire surface for the use of a nonstructured partition. We construct a numerical scheme for the integral equation on the basis of suggested quadrature formulas, prove an estimate for the norm of the inverse matrix of the related system of linear equations and the uniform convergence of numerical solutions to the exact solution of the hypersingular integral equation on the grid.     

剪切载荷作用下含损伤胶接材料界面动应力强度因子的研究

主要针对剪切载荷作用下,胶接材料接合区域界面裂纹尖端动态应力强度因子进行了分析,其中考虑了裂尖区域的损伤．通过积分变换,引入位错密度函数,奇异积分方程被简化为代数方程,并采用配点法求解;最后经过Laplace逆变换,得到动态应力强度因子的时间响应．Ⅱ型动应力强度因子随着黏弹性胶层的剪切松弛参量、弹性基底的剪切模量和Poisson比的增加而增大;随膨胀松弛参量的增加而减小．损伤屏蔽发生在裂纹扩展的起始阶段．裂纹尖端的奇异性指数(-0.5)是与材料参数、损伤程度和时间无关的,而振荡指数由黏弹性材料参数控制．     

压电材料中两平行对称可导通裂纹断裂性能分析

采用Schmidt研究了压电材料中对称平行的双可导通裂纹的断裂性能,利用富里叶变换使问题的求解转换为求解两对以裂纹面位移之差为未知变量的对偶积分方程,并采用Schmidt方法来对这两对对偶积分程进行数值求解。结果表明应力强度因子和电位移强度因子与裂纹的几何尺寸有关。与不可导通裂纹有关结果相比,可导通裂纹的电位移强度因子远小于相应问题不可导通裂纹的电位移强度因子。     

Application of the discrete Fourier transform to solving the Cauchy integral equation

The uniquely solvable system of the Cauchy integral equation of the first kind and index 1 and an additional integral condition is treated. Such a system arises, for example, when solving the skew derivative problem for the Laplace equation outside an open arc in a plane. This problem models the electric current from a thin electrode in a semiconductor film placed in a magnetic field. A fast and accurate numerical method based on the discrete Fourier transform is proposed. Some computational tests are given. It is shown that the convergence is close to exponential.