MULTIPLE PARALLEL SYMMETRIC PERMEABLE MODEL-Ⅲ CRACKS IN A PIEZOELECTRIC/PIEZOMAGNETIC COMPOSITE MATERIAL PLANE

In this paper, the interactions of multiple parallel symmetric and permeable finite length cracks in a piezoelectric/piezomagnetic material plane subjected to anti-plane shear stress loading are studied by the Schmidt method.The problem is formulated through Fourier transform into 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 are directly expanded as a series of Jacobi polynomials.Finally, the relation between the electric field, the magnetic flux field and the stress field near the crack tips is obtained.The results show that the stress, the electric displacement and the magnetic flux intensity factors at the crack tips depend on the length and spacing of the cracks.It is also revealed that the crack shielding effect presents in piezoelectric/piezomagnetic materials.     

Interactions of multiple parallel symmetric permeable mode-III cracks in a piezoelectric material plane

In this paper, the interactions of multiple parallel symmetric and permeable finite length cracks in a piezoelectric material plane subjected to anti-plane shear stress loading were studied by the Schmidt method. The problem was formulated through Fourier transform into dual integral equations, in which the unknown variables are the jumps of displacements across the crack surfaces. To solve the dual integral equations, the jumps of displacements across the crack surfaces were directly expanded as a series of Jacobi polynomials. Finally, the relation between the electric field and the stress field near the crack tips was obtained. The results show that the stress and the electric displacement intensity factors at the crack tips depend on the lengths and spacing of the cracks. It is also revealed that the crack shielding effect presents in piezoelectric materials.     

Coupled field state around three parallel non-symmetric cracks in a piezoelectric/piezomagnetic material plane

In this paper, the behavior of three parallel non-symmetric permeable cracks in a piezoelectric/piezomagnetic material plane subjected to anti-plane shear stress loading was studied by the Schmidt method. The problem was formulated through Fourier transform into three pairs of dual integral equations, in which unknown variables are jumps of displacements across the crack surfaces. To solve the dual integral equations, the jumps of displacements across the crack surfaces were directly expanded as a series of Jacobi polynomials. Finally, the relations among the electric displacement, the magnetic flux and the stress fields near the crack tips can be obtained. The results show that the stress, the electric displacement and the magnetic flux intensity factors at the crack tips depend on the lengths and spacing of cracks. It was also revealed that the crack shielding effect is present in piezoelectric/piezomagnetic materials.     

Basic solution of two parallel non-symmetric permeable cracks in piezoelectric materials

The behavior of two parallel non-symmetric cracks in piezoelectric materials subjected to the anti-plane shear loading was studied by the Schmidt method for the permeable crack electric boundary conditions. Through the Fourier transform, the present problem can be solved with two pairs of dual integral equations ip which the unknown variables are the jumps of displacements across crack surfaces. To solve the dual integral equations, the jumps of displacements across crack surfaces were directly expanded in a series of Jacobi polynomials. Finally, the relations between electric displacement intensity factors and stress intensity factors at crack tips can be obtained. Numerical examples are provided to show the effect of the distance between two cracks upon stress and electric displacement intensity factors at crack tips. Contrary to the impermeable crack surface condition solution, it is found that electric displacement intensity factors for the permeable crack surface conditions are much smaller than those for the impermeable crack surface conditions. At the same time, it can be found that the crack shielding effect is also present in the piezoelectric materials.     

Analysis of stress intensity factors of a ring-shaped interface crack

In this paper, numerical solutions of singular integral equations are discussed in the analysis of axi-symmetric interface cracks under torsion and tension. The problems of a ring-shaped interface crack are formulated in terms of a system of singular integral equations on the basis of the body force method. In the numerical analysis, unknown body force densities are approximated by the products of the fundamental density functions and power series, where the fundamental densities are chosen to express a two-dimensional interface crack exactly. The accuracy of the present analysis is verified by comparing the present results with the results obtained by other researchers for the limiting cases of the geometries. The calculation shows that the present method gives rapidly converging numerical results for those problems as well as for ordinary crack problems in homogeneous material. The stress intensity factors of a ring-shaped interface crack are shown in tables and charts with varying the material combinations and also geometrical conditions.     

Analysis of multiple interfacial cracks in three-dimensional bimaterials using hypersingular integro-differential equation method

By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.     

The dynamic behavior of two parallel symmetric cracks in functionally graded piezoelectric/piezomagnetic materials

The dynamic behavior of two parallel symmetric cracks in functionally graded piezoelectric/piezomagnetic materials subjected to harmonic antiplane shear waves is investigated using the Schmidt method. The present problem can be solved using the Fourier transform and the technique of dual integral equations, in which the unknown variables are jumps of displacements across the crack surfaces, not dislocation density functions. To solve the dual integral equations, the jumps of displacements across the crack surfaces are directly expanded as a series of Jacobi polynomials. Finally, the relations among the electric, magnetic flux, and dynamic stress fields near crack tips can be obtained. Numerical examples are provided to show the effect of the functionally graded parameter, the distance between the two parallel cracks, and the circular frequency of the incident waves upon the stress, electric displacement, and magnetic flux intensity factors at crack tips.     

纤维增强复合材料圆柱型界面裂纹分析

以裂纹面上的位错函数为未知量将圆柱型界面裂纹问题化成一组奇异积分方程的求解问题．应用Ｍｕｓｋｈｅｌｉｓｈｖｉｌｉ的奇异积分方程理论，分析了圆柱型界面裂纹尖端应力场．针对裂纹尖端分别存在和不存在接触区两种情况，确定了裂纹尖端应力场的奇异性．利用数值方法计算了圆柱型界面裂纹尖端接触区尺寸对剪应力强度因子的影响．     

Analytical solution of multiple moving cracks in functionally graded piezoelectric strip

The dynamic behaviors of several moving cracks in a functionally graded piezoelectric (FGP) strip subjected to anti-plane mechanical loading and in-plane electrical loading are investigated. For the first time, the distributed dislocation technique is used to construct the integral equations for FGP materials, in which the unknown variables are the dislocation densities. With the dislocation densities, the field intensity factors are determined. Moreover, the effects of the speed of the crack propagation on the field intensity factors are studied. Several examples are solved, and the numerical results for the stress intensity factor and the electric displacement intensity factor are presented graphically finally.     

Nonlocal theory solution of two collinear cracks in the functionally graded materials

In this paper, the interaction of two collinear cracks in functionally graded materials subjected to a uniform anti-plane shear loading is investigated by means of nonlocal theory. The traditional concepts of the nonlocal theory are extended to solve the fracture problem of functionally graded materials. To make the analysis tractable, it is assumed that the shear modulus varies exponentially with the coordinate vertical to the crack. By use of the Fourier transform, the problem can be solved with the help of a pair of triple integral equations, in which the unknown variable is the displacement on the crack surfaces. To solve the triple integral equations, the displacement on the crack surfaces is expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present near the crack tips. The nonlocal elastic solutions yield a finite hoop stress at the crack tip, thus allowing us to use the maximum stress as a fracture criterion in functionally graded materials. The magnitude of the finite stress field depends on the crack length, the distance between two cracks, the parameter describing the functionally graded materials and the lattice parameter of the materials.     

The dynamic behavior of two collinear interface cracks in magneto-electro-elastic materials

In this paper, the dynamic behavior of two collinear symmetric interface cracks between two dissimilar magneto-electro-elastic material half planes under the harmonic anti-plane shear waves loading is investigated by Schmidt method. By using the Fourier transform, the problem can be solved with a set of triple integral equations in which the unknown variable is the jump of the displacements across the crack surfaces. To solve the triple integral equations, the jump of the displacements across the crack surface is expanded in a series of Jacobi polynomials. Numerical solutions of the stress intensity factor, the electric displacement intensity factor and the magnetic flux intensity factor are given. The relations among the electric filed, the magnetic flux field and the stress field are obtained.     

Solutions to a limited-permeable crack or two limited-permeable collinear cracks in piezoelectric/piezomagnetic materials

The solutions of a limited-permeable crack (case I) or two collinear limited-permeable cracks (case II) in piezoelectric/piezomagnetic materials subjected to a uniform tension loading were investigated in this paper using the generalized Almansi’s theorem. At the same time, the electric permittivity and the magnetic permeability of air in crack were firstly considered. Through the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations, in which the unknown variables were jumps of displacements across crack surfaces, not the dislocation density functions or the complex variable functions. To solve the dual integral equations, the jumps of displacements across crack surfaces were directly expanded in a series of Jacobi polynomials to obtain the relations among electric displacement intensity factors, magnetic flux intensity factors and stress intensity factors at crack tips.     

反平面弹性中边缘内分叉裂纹的奇异积分方程解法

利用复变函数和奇异积分方程方法,求解反平面弹性中半平面边缘内分叉裂纹问题。提出了满足半平面边界自由的由分布位错密度表示的半平面中单裂纹的基本解,此基本解由主要部分和辅助部分组成。将半平面边缘内分叉裂纹问题看作是许多单裂纹问题的叠加,建立了以分布位错密度为未知函数的Cauchy型奇异积分方程组。然后,利用半开型积分法则求解奇异积分方程,得到了裂纹端处的应力强度因子。文中给出两个数值算例的计算结果。     

Interaction between a circular inclusion and a symmetrically branched crack

The interaction problem between a circular inclusion and a symmetrically branched crack embedded in an infinite elastic medium is solved. The branched crack is modeled as three straight cracks which intersect at a common point and each crack is treated as a continuous contribution of edge dislocations. Green's functions are used to reduce the problem into a system of singular equations consisting of the distributions of Burger's dislocation vectors as unknown functions through the superposition technique. The resulting integral equations are solved numerically by the method of Gauss-Chebychev quadrature. The proposed integral equation approach is first verified for two limiting cases against the literature. More effort is paid on the effect of inclusion on both the Mode I and Mode lI stress intensity factors at the branch tips. The effect of inclusion on the branching path is also investigated.     

基于修正权函数的多裂纹无网格模拟

本文采用了一种基于不连续场修正权函数的无网格方法来处理二维平面多裂纹问题。相较于传统的无网格断裂不连续场和奇异场模拟方法,修正权函数法算法简便易实现。采用修正权函数处理多裂纹时,只需要对每一段裂纹周围节点的权函数进行修正,就能同时模拟多裂纹不连续位移场和多裂尖奇异场。本文采用基于不连续场修正权函数的无单元Galerkin方法（EFGM）,对Y型裂纹板、十字型裂纹板和孔边双裂纹板进行了分析。数值结果表明,在不引入扩展基函数情况下,通过修正权函数法能够得到精度较高的应力强度因子解,能较好地拟合多裂纹的裂尖奇异场。     

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

In this paper, the dynamic behavior of two collinear cracks in the anisotropic elasticity material plane subjected to the harmonic anti-plane shear waves is investigated by use of the nonlocal theory. To overcome the mathematical difficulties, a one-dimensional nonlocal kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress field near the crack tips. By use of the Fourier transform, the problem can be solved with the help of a pair of triple integral equations, in which the unknown variable is the displacement on the crack surfaces. To solve the triple integral equations, the displacement on the crack surfaces is expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present near crack tips. The nonlocal elasticity solutions yield a finite hoop stress at the crack tips, thus allowing us to using the maximum stress as a fracture criterion. The magnitude of the finite stress field not only depends on the crack length but also on the frequency of the incident waves and the lattice parameter of the materials.     

Asymptotic elastic solution for two straight cracks of arbitrary length and location

An infinite elastic plane containing two straight cracks of arbitrary length and location is analyzed within the framework of elastostatics. The mathematical formulation is based on the stress solution for a single crack and leads to a system of singular integral equations that govern the crack surface displacement densities. The solution series in terms of the reciprocal of the crack centre distance is not suitable for cracks that are spaced too closely. It is shown by way of examples that the method of asymptotic solution is convenient for developing approximation expressions of the stress and displacement field with certain characteristics. The formulas for the stress intensity factors and crack opening are given for the case of a constant tensile load. Graphical results are given for the variations of the stress intensity factors with parameters depending on the relative positions of the cracks.     

COMPUTATION OF STRESS INTENSITY FACTORS BY THE SUB-REGION MIXED FINITE ELEMENT METHOD OF LINES

Based on the sub-region generalized variational principle,a sub-region mixed ver- sion of the newly-developed semi-analytical‘finite element method of lines’(FEMOL)is pro- posed in this paper for accurate and efficient computation of stress intensity factors(SIFs)of two-dimensional notches/cracks.The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used,with the sought SIFs being among the unknown coefficients.The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements.A mixed system of ordinary differential equations(ODEs) and al- gebraic equations is derived via the sub-region generalized variational principle.A singularity removal technique that eliminates the stress parameters from the mixed equation system even- tually yields a standard FEMOL ODE system,the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver.A number of numerical examples,including bi-material notches/cracks in anti-plane and plane elasticity,are given to show the generally excellent performance of the proposed method.     

含圆形嵌体弹性平面中径向裂纹问题的超奇异积分方程方法

根据含圆形嵌体平面问题在极坐标下的弹性力学基本解,使用Betti互换定理,在有限部积分意义下将问题归结为两个以裂纹岸位移间断为基本未知量、对于Ⅰ型和Ⅱ型问题相互独立的超奇异积分方程,对含圆形嵌体弹性平面中的径向裂纹问题进行了研究.根据有限部积分原理,建立了问题的数值算法.计算结果表明,嵌体半径、裂纹位置及材料剪切弹性模量等都对裂纹应力强度因子具有较为明显的影响.     

The nonlocal solution of two parallel cracks in functionally graded materials subjected to harmonic anti-plane shear waves

In this paper, the dynamic interaction of two parallel cracks in functionally graded materials (FGMs) is investigated by means of the non-local theory. To make the analysis tractable, the shear modulus and the material density are assumed to vary exponentially with the coordinate vertical to the crack. To reduce mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the dynamic problem to obtain stress fields near the crack tips. By use 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 displacements across the crack surfaces. To solve the dual integral equations, the jumps of displacements across the crack surfaces are expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present at the crack tips. The non-local elastic solutions yield a finite hoop stress at the crack tips. The present result provides theoretical references helpful for evaluating relevant strength and preventing material failure of FGMs with initial cracks. The magnitude of the finite stress field depends on relevant parameters, such as the crack length, the distance between two parallel cracks, the parameter describing the FGMs, the frequency of the incident waves and the lattice parameter of materials. The project supported by the National Natural Science Foundation of China (90405016, 10572044) and the Specialized Research Fund for the Doctoral Program of Higher Education (20040213034). The English text was polished by Yunming Chen.