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.     

Strip-coalesced interior zone model for two unequal collinear cracks weakening piezoelectric media

In this paper, a mathematical strip-saturation model is proposed for a poled transversely isotropic piezoelectric plate weakened by two impermeable unequal-collinear hairline straight cracks. Remotely applied in-plane unidirectional electromechanical loads open the cracks in mode-I such that the saturation zone developed at the interior tips of cracks gets coalesced. The developed saturation zones are arrested by distributing over their rims in-plane normal cohesive electrical displacement. The problem is solved using the Stroh formalism and the complex variable technique. The expressions are derived for the stress intensity factors （SIFs）, the lengths of the saturation zones developed, the crack opening displacement （COD）, and the energy release rate. An illustrative numerical case study is presented for the poled PZT-5H ceramic to investigate the effect of prescribed electromechanical loads on parameters affecting crack arrest. Also, the effect of different lengths of cracks on the SIFs and the local energy release rate （LERR） has been studied. The results obtained are graphically presented and analyzed.     

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.     

SELF-SIMILAR CRACK EXPANSION METHOD FOR THE STRESS INTENSITY FACTOR ANALYSIS

The Self-Similar Crack Expansion (SSCE) method is proposed to evaluate stress intensi-ty factors at crack tips, whereby stress intensity factors of a crack can be determined by the crackopening displacement over the crack, not just by the local displacement around the crack tip. The crackexpansion rate is estimated by taking advantage of the crack self-similarity. Therefore, the accuracy ofthe calculation is improved. The singular integrals on crack tip elements are also analyzed and are pre-cisely evaluated in terms of a special integral analysis. Combination of these two techniques greatly in-creases the accuracy in estimating the stress distribution around the crack tip. A variety of two-dimen-sional cracks, such as subsurface cracks, edge cracks, and their interactions are calculated in terms ofthe self-similar expansion rate. Solutions are satisfied with errors less than 0.5% as compared with theanalytical solutions. Based on the calculations of the crack interactions, a theory for crack interactionsis proposed such that for a group of aligned cracks the summation of the square of SIFs at the right tipsof cracks is always equal to that at the left tips of cracks. This theory was proved by the mehtod ofSelf-Similar Crack Expansion in this paper.     

Thermal stress intensity factors for a crack in an anisotropic half plane

On the basis of the two-dimensional theory of anisotropic thermoelasticity, a solution is given for the thermal stress intensity factors due to the obstruction of a uniform heat flux by an insulated line crack in a generally anisotropic half plane. The crack is replaced by continuous distributions of sources of temperature discontinuity and dislocations. First, the particular thermoelastic dislocation solutions for the half plane are obtained; then the corresponding isothermal solutions are superposed to satisfy the traction-free conditions on the crack surfaces. The dislocation solutions are applied to calculate the thermal stress intensity factors, which are validated by the exact solutions. The effects of the uniform heat flux, the ply angle and the crack length are investigated.     

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.     

THERMAL FRACTURE OF FUNCTIONALLY GRADED PLATE WITH PARALLEL SURFACE CRACKS

This work examines the fracture behavior of a functionally graded material （FGM） plate containing parallel surface cracks with alternating lengths subjected to a thermal shock. The thermal stress intensity factors （TSIFs） at the tips of long and short cracks are calculated using a singular integral equation technique. The critical thermal shock △Tc that causes crack initiation is calculated using a stress intensity factor criterion. Numerical examples of TSIFs and △Tc for an Al2O3/Si3N4 FGM plate are presented to illustrate the effects of thermal property gradation, crack spacing and crack length ratio on the TSIFs and △Tc. It is found that for a given crack length ratio, the TSIFs at the tips of both long and short cracks can be reduced significantly and △Tc can be enhanced by introducing appropriate material gradation. The TSIFs also decrease dramatically with a decrease in crack spacing. The TSIF at the tips of short cracks may be higher than that for the long cracks under certain crack geometry conditions. Hence, the short cracks instead of long cracks may first start to grow under the thermal shock loading.     

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.     

压电体中孔边Ⅲ型界面裂纹的动应力强度因子

采用Green函数法研究界面上含圆孔边界径向有限长度裂纹的两半无限压电材料对SH波的散射和裂纹尖端动应力强度因子问题.首先构造出具有半圆型凹陷半空间的位移Green函数和电场Green函数,然后采用裂纹"切割"方法构造孔边裂纹,并根据契合思想和界面上的连接条件建立起求解问题的定解积分方程.最后作为算例,给出了孔边界面裂纹尖端动应力强度因子的计算结果图并进行了讨论.     

界面下主微裂纹间干涉的研究

涉及两相正交各向异性体界面干涉问题的研究,多裂纹问题被分解为只含单裂纹的子问题,利用位错理论和裂面应力自由条件,列出一组可数值求解位错密度函数的奇异积分方程,从耐 注得应力强度因子。     

基于权函数法的核主泵飞轮应力强度因子研究

采用权函数法确定了含裂纹飞轮在离心力和接触压力作用下的应力强度因子,计算了在同步转速工况下裂纹尖端应力强度因子的值,并与解析法和有限元法计算结果进行了对比。结果表明:权函数法与解析法的误差在3%以内,与有限元法的误差在1%以内,验证了权函数法计算应力强度因子的准确性高;在结构不变的情况下,权函数法可以快速求解不同载荷条件、不同长度裂纹的应力强度因子。通过控制变量法研究了不同参数对应力强度因子的影响,结果表明,飞轮裂纹尖端总应力强度因子随着裂纹尺寸、旋转转速、飞轮尺寸外径与内径比值的增大而增大。     

拉压循环载荷下正交异性钢桥面板肋-面板焊缝裂纹扩展分析

采用四步法计算了考虑循环载荷中压应力影响的正交异性钢桥面板的肋-面板焊缝表面裂纹扩展。第一步是基于正交异性钢桥面板的疲劳分析模型,计算肋-面板焊缝处的应力,第二步是通过肋-面板焊缝的三维局部模型,用Schwartz-Neumann交替法计算焊缝表面裂纹的应力强度因子分布,第三步是用二维断裂力学模型和增量塑性损伤模型,计算循环载荷中的压应力对裂纹扩展的影响,第四步是用第二步中的三维裂纹分析结果和第三步中的二维断裂力学模型得到的裂纹扩展公式,计算钢桥面板的肋-面板焊缝表面裂纹扩展。计算结果表明,对应于正交异性钢桥面板肋-面板焊缝处的循环应力,本文所用模型的裂纹尖端反向塑性区导致裂纹扩展率增加50%以上。研究结果为正交异性钢桥面板肋-面板焊缝裂纹的疲劳寿命分析提供了研究基础。     

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

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

Crack bifurcation predicted for dynamic anti-plane collinear cracks in piezoelectric materials using a non-local theory

A non-local theory of elasticity is applied to obtain the dynamic interaction between two collinear cracks in the piezoelectric materials plane under anti-plane shear waves for the permeable crack surface boundary conditions. Unlike the classical elasticity solution, a lattice parameter enters into the problem that make the stresses and the electric displacements finite at the crack tip. A one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress and electric displacement near the crack tips. By means of the Fourier transform, the problem can be solved with the help of two pairs of triple integral equations in which the unknown variable is the jump of the displacement across the crack surface. The solutions are obtained by means of the Schmidt method. Crack bifurcation is predicted using the strain energy density criterion. Minimum values of the strain energy density functions are assumed to coincide with the possible locations of fracture initiation. Bifurcation angles of ±5° and ±175° are found. The result of possible crack bifurcation was not expected before hand.     

Investigation of the dynamic behavior of two parallel symmetric cracks in piezoelectric materials use of non-local theory

In this paper, the dynamic behavior of two parallel symmetric cracks in piezoelectric materials under harmonic anti-plane shear waves is investigated by use of the non-local theory for permeable crack surface conditions. To overcome the mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the problem to obtain the stress occurs near the crack tips. By means of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations that the unknown variables are the jumps of the displacement along the crack surfaces. These equations are solved using the Schmidt method. Numerical examples are provided. Contrary to the previous results, it is found that no stress and electric displacement singularity is present near the crack tip. The non-local elastic solutions yield a finite hoop stress near 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, the frequency of the incident wave, the distance between two cracks and the lattice parameter of the materials, respectively. Contrary to the impermeable crack surface condition solution, it is found that the dynamic electric displacement for the permeable crack surface conditions is much smaller than the results for the impermeable crack surface conditions. The results show that the dynamic field will impede or enhance crack propagation in the piezoelectric materials at different stages of the dynamic load.     

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.     

The stress field in the neighborhood of a branched crack in an infinite elastic sheet

The problem of a branched crack consisting of a main crack and a straight branch starting from one of its tip located in an infinite elastic sheet is considered under the assumptions of two-dimensional theory of Elasticity. Employing Kolosov-Muskhelishvili representation of the stress function and other well known techniques the problem is reduced to the solution of an integral equation. The nature of the stress singularity at the re-entrant corner, where the two branches of the crack meet, is discussed. Based upon a numerical solution of the integral equation the stress intensity factors at the two tips are computed for two types of prescribed traction at infinity and various geometric configurations of the branched crack.     

INTERFACE CRACK IN BIPIEZOTHERMOELASTIC MEDIA

By using the extended version of Eshelby-Stroh's formulation and the method of analyt-ical continuation,the problems of interface cracks are reduced to a Hilbert problem of vector form.Ageneral explicit closed form solution for the piezothermoelastic interface crack problem is then ob-tained,the whole field solutions of temperature,heat flux,displacements,electric field,stress andelectric induction are given,the explicit expressions for the crack opening displacements and electricpotential are also provided.     

Uniformly moving screw dislocation interacting with interface cracks in anisotropic bimaterials

This paper attempts to investigate the problem for the interaction between a uniformly subsonic moving screw dislocation and interface cracks in two dissimilar anisotropic materials. Using Riemann–Schwarz’s symmetry principle integrated with the analysis singularity of complex functions, we present the general elastic solutions of this problem and the closed form solutions for interface containing one and two cracks. The expressions of stress intensity factors at the crack tips and image force acting on moving dislocation are derived explicitly. The results show that the stress intensity factors at the crack tips decrease with increasing velocity of dislocation, and larger dislocation velocity leads to the equilibrium position of dislocation leaving from crack tips. The presented solutions contain previously known results as the special cases.     

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.