压电介质夹杂功能梯度压电带界面双裂纹反平面问题

利用积分方程方法,本文研究了夹在两个均匀压电半空间的功能梯度压电带界面共线双裂纹的反平面问题。在电渗透型边界条件下,通过Fourier余弦变换将所考虑的问题化为一对偶积分方程,再用Copson方法将该对偶积分方程转化为Fredholm方程进行数值求解,从而给出了裂纹尖端的应力强度因子,电位移强度因子的表达式。分析了裂纹长度,功能梯度非均匀参数以及材料的几何尺寸等对应力强度因子的影响。     

压电拼接电磁复合材料中裂纹对SH波的散射

利用积分变换及奇异积分方程技术研究电磁复合材料底层处裂纹对SH波散射问题。假定裂纹面的边界条件为电渗透性,通过Fourier余弦变化将问题转化为对偶积分方程,并利用Copson方法将对偶积分方程转化为第二类Fredholm积分方程求解。给出标准动应力强度因子表达式;通过数值计算分析裂纹长度、裂纹到界面距离、入射波频率及入射角对标准动应力强度因子影响。     

压电-压磁板条中反平面裂纹的电-磁-弹性分析

基于线性电磁弹性理论,获得了压电-压磁板条中反平面裂纹尖端附近的奇异应力、电场和磁场。假设裂纹位于和板条边界平行的中心位置,并且裂纹是电磁渗透型的。利用Fourier变换,将裂纹面的混合边值问题化为对偶积分方程,即而归结为第二类Fredholm积分方程。通过渐近分析,得到了裂纹尖端附近应力、应变、电位移、电场、磁场和磁感的封闭表达式。结果表明,对于电磁渗透裂纹,电场强度因子和磁场强度因子总为0;板条的宽度对应力强度因子有显著的影响;能量释放率总为正值。     

金属裂纹板复合材料修补结构的超奇异积分方程方法

根据应力强度因子在线弹性范围内具有可叠加性,将金属裂纹板复合材料修补结构进行简化,在表面裂纹线弹簧模型的基础上,建立了基于超奇异积分方程的Line-Spring模型。利用第二类Chebyshev多项式展开的方法,将超奇异积分方程转化为线性方程组,推导出以裂纹面位移表示的应力强度因子表达式,得到了裂纹尖端应力强度因子的数值解,并利用虚拟裂纹闭合法加以验证。参数分析确定了影响对称修补裂纹板应力强度因子的两个主要参数:胶层界面刚度和补片与金属板刚度比,为胶接修补结构的承载能力分析以及改进设计提供理论依据。     

方形截面管横向裂纹的应力强度因子KI

利用裂纹张开能量释放率建立了一个求解方形截面管横向裂纹应力强度因子的一个方法。给出了方形截面管裂纹张开能量释放率的 G*-积分表征,以及和应力强度因子的关系。同时也给出了 G*-积分与载荷、几何参量以及机械性能参数的关系,进而得到方形截面管横向裂纹的应力强度因子。给出的方法不仅适用于一般箱形结构件的裂纹问题,也适用于其它有限边界多边管状结构的三维裂纹问题,过程极为简单。     

加层压电条界面裂纹的稳态扩展

研究当压电条同时与两个不同材料的弹性条粘接在一起,在反平面机械载荷及面内电载荷联合作用下,长度不变的有限Griffith 界面裂纹沿加层压电条界面以常速稳态扩展时裂尖的动态断裂问题。应用Fourier积分变换将问题化为以第二类Fredholm积分方程表示的对偶积分方程,导出了相应的动应力强度因子表达式。给出了动应力强度因子与裂纹传播速度、裂纹长度、压电条及弹性条厚度、电荷载大小及方向的关系曲线。研究结果对结构设计及结构失效的预防具有理论和应用价值。     

多铁性非均匀空心层合柱圆弧界面的反平面断裂分析

讨论反平面载荷作用下多铁性非均匀空心层合柱的圆弧界面裂纹问题,层合柱由梯度铁电层和梯度铁磁层粘接而成,界面处存在圆弧型裂纹。采用分离变量和Cauchy核奇异积分方程方法求解该断裂问题。通过讨论断裂参数的数值解,分析了梯度非均匀参数、几何与材料参数变动等对应力强度因子的影响。     

含任意方向裂纹功能梯度材料的应力分析研究

功能梯度材料是在航空航天领域的需求背景下发展起来的,但由于生产技术及工作环境等方面的原因,功能梯度材料内部常常产生各种形式的裂纹并最终导致材料破坏,因此研究含任意方向裂纹功能梯度材料的断裂问题具有重要意义。以含有任意方向裂纹的功能梯度材料为对象,运用积分变换方法,给出了相应材料平面问题的位移场的形式解。通过引入辅助函数并利用相关条件,可将问题转化为求解一组带有Cauchy核的奇异积分方程,继而采用Lobatto-Chebyshev方法对奇异积分方程进行数值求解。最后分析了裂纹方向、材料非均匀指数、载荷条件对混合型应力强度因子的影响。      

非均匀复合材料中反平面裂纹的动态断裂力学研究

对于非均匀复合材料中多个裂纹的动态断裂力学问题, 提出了一种分析方法, 假设复合材料为正交各向异性并含有多个垂直于厚度方向的裂纹, 材料参数沿厚度方向为变化的, 沿该方向将复合材料划分为许多单层, 假设单层材料参数为常数, 应用柔度矩阵/刚度矩阵方法及Fourier变换法, 在L aplace 域内推导出了控制问题的奇异积分方程组, 并用虚位移原理求解, 给出了应力强度因子及能量释放率的表达式, 然后利用Laplace 数值反演, 得出了裂纹尖端的动态应力强度因子和能量释放率。作为算例, 研究了带有两个裂纹的功能梯度结构, 分析了材料参数的优化对降低应力强度因子的意义。      

非均匀材料界面裂纹的Cell-Based光滑有限元法

为提高非均匀材料界面裂纹尖端断裂参数的求解精度,基于非均匀材料界面断裂力学、Cell-Based光滑有限元(Cell-SFEM)和非均匀材料的互交作用积分法,提出了求解非均匀材料界面裂纹尖端断裂参数的CellBased光滑有限元法,推导了基于Cell-Based光滑有限元法的非均匀材料的互交作用积分法,对非均匀材料间的界面裂纹尖端处正则应力强度因子进行了求解,并与参考解进行了比较,讨论了互交积分区域大小和光滑子元个数与正则应力强度因子的关系。数值算例结果表明:本方法具有很高的计算精度,对积分区域大小不敏感,可为设计、制造抗破坏非均匀材料提供依据。     

Scattering of the harmonic anti-plane shear waves by a crack in functionally graded piezoelectric materials

In the present paper, the dynamic behavior of a Griffth crack in the functionally graded piezoelectric material (FGPM) is investigated. It is assumed that the elastic stiffness, piezoelectric constant, dielectric permittivity and mass density of the FGPM vary continuously as an exponential function, and that FGPM is under the anti-plane mechanical loading and in-plane electrical loading. By using the Fourier transform and defining the jumps of displacement and electric potential components across the crack surface as the unknown functions, two pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement and electric potential components across the crack surface are expanded in a series of Jacobi polynomial. Numerical examples are provided to show the effects of material properties on the stress and the electric displacement intensity factors.     

Crack propagating in a functionally graded strip under the plane loading

In the present paper, a finite crack with constant length (Yoffe type crack) propagating in the functionally graded strip under the plane loading is investigated by means of the Schmidt method. By using the Fourier transform and defining the jumps of displacement components across the crack surface as the unknown functions, two pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement components across the crack surface are expanded in a series of Jacobi polynomial. Numerical examples are provided to show the effects of the material properties, the thickness of the functionally graded strip, and speed of the crack propagating upon the dynamic fracture behavior.     

On the Moving Griffith Crack in a Nonhomogeneous Orthotropic Strip

A finite crack with constant length (Yoffe type crack) propagating in the functionally graded orthotropic strip under the plane loading is investigated by means of the Schmidt method. By using the Fourier transform and defining the jumps of displacement components across the crack surface as the unknown functions, two pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement components across the crack surface are expanded in a series of Jacobi polynomial. Numerical examples are provided to show the effects of material properties, the thickness of the functionally graded orthotropic strip and the speed of the crack propagating upon the dynamic fracture behavior.     

Dynamic stress intensity factors of two 3D rectangular permeable cracks in a transversely isotropic piezoelectric material under a time‐harmonic elastic P‐wave

The dynamic behavior of two 3D rectangular permeable cracks in a transversely isotropic piezoelectric material is investigated under an incident harmonic stress wave by using the generalized Almansi's theorem and the Schmidt method. The problem is formulated through double Fourier transform into three pairs of dual integral equations with the displacement jumps across the crack surfaces as the unknown variables. To solve the dual integral equations, the displacement jumps across the crack surfaces are directly expanded as a series of Jacobi polynomials. Finally, the relations among the dynamic stress field and the dynamic electric displacement filed near the crack edges are obtained, and the effects of the shape of the rectangular crack, the characteristics of the harmonic wave, and the distance between two rectangular cracks on the stress and the electric intensity factors in a piezoelectric composite material are analyzed. Copyright © 2015 John Wiley & Sons, Ltd.     

Non-local theory solution of a mode-I crack in a piezoelectric/piezomagnetic composite material plane

The non-local theory solution of a mode-I permeable crack in a piezoelectric/piezomagnetic composite material plane was given by using the generalized Almansi’s theorem and the Schmidt method in this paper. The problem was formulated through Fourier transform into two pairs of dual integral equations, in which the unknown variables are the displacement jumps across the crack surfaces. To solve the dual integral equations, the displacement jumps across the crack surfaces were directly expanded as a series of Jacobi polynomials. Numerical examples were provided to show the effects of the crack length and the lattice parameter on the stress field, the electric displacement field and the magnetic flux field near the crack tips. Unlike the classical elasticity solutions, it is found that no stress, electric displacement and magnetic flux singularities are present at the crack tips in piezoelectric/piezomagnetic composite materials. The non-local elastic solution yields a finite hoop stress at the crack tip, thus allowing us to use the maximum stress as a fracture criterion.     

Dynamic behavior of a crack in a functionally graded piezoelectric strip bonded to two dissimilar half piezoelectric material planes

Summary. The dynamic behavior of a crack in a functionally graded piezoelectric material (FGPM) strip bonded to two half dissimilar piezoelectric material planes subjected to combined harmonic anti-plane shear wave and in-plane electrical loading was studied under the limited permeable and permeable electric boundary conditions. It was assumed that the elastic stiffness, piezoelectric constant and dielectric permittivity of the functionally graded piezoelectric layer vary continuously along the thickness of the strip. By using the Fourier transform, the problem can be solved with a set of dual integral equations in which the unknown variables are the jumps of the displacements and the electric potentials across the crack surfaces. In solving the dual integral equations, the jumps of the displacements and the electric potentials across the crack surfaces were expanded in a series of Jacobi polynomials. Numerical results illustrate the effects of the gradient parameter of FGPM, electric loading, wave number, thickness of FGPM strip and electric boundary conditions on the dynamic stress intensity factors (SIFs).     

Basic solutions of two parallel mode-I cracks or four parallel mode-I cracks in the piezoelectric materials

In this paper, the stress and the electric intensity factors of two parallel mode-I cracks or four parallel mode-I cracks in the piezoelectric materials were examined by means of the Schmidt method for the permeable electric boundary conditions. The present problem can be solved by using the Fourier transform and the technique of dual integral equation, in which the unknown variables are the jumps of displacements across the crack surfaces, not the dislocation density functions. To solve the dual integral equations, the displacement jumps are directly expanded in a series of Jacobi polynomials. Finally, the effects of the distance between two parallel cracks and the distance between two collinear cracks on the stress and the electric intensity factors in the piezoelectric materials are analyzed. These results can be used for the strength evaluation of the piezoelectric materials with multi-cracks.     

A time‐dependent formulation of dual boundary element method for 3D dynamic crack problems

In this paper, the dual boundary element method in time domain is developed for three‐dimensional dynamic crack problems. The boundary integral equations for displacement and traction in time domain are presented. By using the displacement equation and traction equation on crack surfaces, the discontinuity displacement on the crack can be determined. The integral equations are solved numerically by a time‐stepping technique with quadratic boundary elements. The dynamic stress intensity factors are calculated from the crack opening displacement. Several examples are presented to demonstrate the accuracy of this method. Copyright © 1999 John Wiley & Sons, Ltd     

Dynamic SIF of interface crack between two dissimilar viscoelastic bodies under impact loading*

In this paper, the transient dynamic stress intensity factor (SIF) is determined for an interface crack between two dissimilar half-infinite isotropic viscoelastic bodies under impact loading. An anti-plane step loading is assumed to act suddenly on the surface of interface crack of finite length. The stress field incurred near the crack tip is analyzed. The integral transformation method and singular integral equation approach are used to get the solution. By virtue of the integral transformation method, the viscoelastic mixed boundary problem is reduced to a set of dual integral equations of crack open displacement function in the transformation domain. The dual integral equations can be further transformed into the first kind of Cauchy-type singular integral equation (SIE) by introduction of crack dislocation density function. A piecewise continuous function approach is adopted to get the numerical solution of SIE. Finally, numerical inverse integral transformation is performed and the dynamic SIF in transformation domain is recovered to that in time domain. The dynamic SIF during a small time-interval is evaluated, and the effects of the viscoelastic material parameters on dynamic SIF are analyzed.