Fredholm integral equation approach for multiple crack problem in antiplane shear and inplane electric field of piezoelectric materials

In this paper, the basic presentation in antiplane shear and inplane electric field of piezoelectric materials is refreshed. In order that the functions used in the formulation can be distinguished by their usage, four analytic functions, or four complex potentials, are introduced. A multiple crack problem for piezoelectric materials is studied. After taking the traction or the electric displacement on the crack face as unknown functions, one can naturally obtain a Fredholm integral equation for the multiple crack problem. It is found that the Fredholm integral equation approach is effective for solving the multiple crack problem. Finally, numerical examples are given.     

Higher order fields for damaged nonlinear antiplane shear notch, crack and inclusion problems

Damaged nonlinear antiplane shear problems with a variety of singularities are studied analytically. A deformation plasticity theory coupled with damage is employed in analysis. The effect of microscopic damage is considered in terms of continuum damage mechanics approach. An exact solution for the general damaged nonlinear singular antiplane shear problem is derived in the stress plane by means of a hodograph transformation, then corresponding higher order asymptotic solutions are obtained by reversing the stress plane solution to the physical plane. As example, traction free sharp notch and crack, rigid sharp wedge and flat inclusion, and mixed boundary sharp notch problems are investigated, respectively. Consequently, higher order fields are obtained, in which analytical expressions of the dominant and second order singularity exponents and angular distribution functions of the near tip fields are derived. Effects of the damage and hardening exponents of materials and the geometric angle of notch/wedge on the near tip quantities are discussed in detail. It is found that damage leads to a weaker dominant singularity of stress, but to little stronger singularities of the dominant and second order terms of strain compared to that for undamaged material. It is also seen that damage has important effect on the angular distribution functions of the near tip stress and strain fields. As special cases, higher order analytical solutions of the crack and rigid flat inclusion tip fields are obtained, respectively, by reducing the notch/wedge tip solutions. Effects of damage and hardening exponents on the dominant and second order terms in the solutions of the crack and inclusion tip fields are discussed.     

Interfacial crack problem in layered soft ferromagnetic materials in uniform magnetic field

The magnetoelastic plane strain problem of an interfacial Griffith crack between two dissimilar soft ferromagnetic elastic materials subjected to a uniform magnetostatic field is considered within the framework of linear magnetoelasticity. By making use of the Fourier integral transform technique, the mixed boundary problem is then reduced to a pair of singular integral equations of the second kind. Solutions of the singular integral equations are obtained numerically by means of a Jacobi polynomial expansion method. Effects of the magnetic field, the combinations of the magnetic properties of materials and the geometric parameters on the magnetoelastic stress intensity factors in the vicinity of crack tip are shown graphically.     

A Crack with an Electric Displacement Saturation Zone in an Electrostrictive Material

A crack with an electric displacement saturation zone in an electrostrictive material under purely electric loading is analyzed. A strip saturation model is here employed to investigate the effect of the electrical polarization saturation on electric fields and elastic fields. A closed form solution of electric fields and elastic fields for the crack with the strip saturation zone is obtained by using the complex function theory. It is found that the K _I -dominant region is very small compared to the strip saturation zone. The generalized Dugdale zone model is also employed in order to investigate the effect of the saturation zone shape on the stress intensity factor. Using the body force analogy, the stress intensity factor for the asymptotic problem of a crack with an elliptical saturation zone is evaluated numerically.     

Crack problem under shear loading in cubic quasicrystal

IntroductionQuasicrystalasanewstructureofsolidmatter[1,2 ]bringsprofoundnewideastothetraditionalcondensedmatterphysicsandencouragesconsiderabletheoreticalandexperimentalstudiesonthephysicalandmechanicalpropertiesofthematerial,includingtheelasticitytheoryofthequasicrystal,manyvaluableresultsweregiven[3～ 5 ].Defectsinthematerialwereobservedsoonafterthediscoveryofthequasicrystal[6 ,7].Cracksareonetypeofdefects,theirexistencegreatlyinfluencesthephysicalandmechanicalpropertiesofthequasicrystalinem…     

Moving crack for functionally grated material in an infinite length strip under antiplane shear

Considered is a Yoffe crack in an infinite strip of functionally grated material (FGM) subjected to antiplane shear. The shear moduli in two directions of FGM are assumed to be of exponential form. The dynamic stress intensity factor and strain energy density factor at the crack tip are obtained by using integral transforms and dual-integral equations. The numerical results show that the decrease of the strain energy density factor varies with the shear moduli gradient, and the increase of the strain energy density factor varies with the increase of the moving crack speed. The ratio of shear moduli in material vertical orientation has a great influence on the strain energy density factor.     

Stress intensity factor in a contact problem for a plane crack under an antiplane shear wave

The frictional contact interaction of the finite edges of a plane crack under the action of a normally incident harmonic shear wave that produces antiplane deformation is studied. The influence of the forces of contact interaction on the stress intensity factor is analyzed Published in Prikladnaya Mekhanika, Vol. 43, No. 9, pp. 115–119, September 2007.     

A moving mode-III crack in functionally graded piezoelectric material: permeable problem

In this paper a moving mode-III crack in functionally graded piezoelectric materials (FGPM) is studied. The crack surfaces are assumed to be permeable. The governing equations for FGPM are solved by means of Fourier cosine transform. The mathematical formulation for the permeable crack condition is derived as a set of dual integral equations, which, in turn, are reduced to a Fredholm integral equation of the second kind. The results obtained indicate that the stress intensity factor of moving crack in FGPM depends only on the mechanical loading. The gradient parameter of the FGPM and the moving velocity of the crack do have significant influence on the dynamic stress intensity factor.     

功能梯度压电材料反平面裂纹问题

基于三维弹性理论和压电理论，导出了材料系数在横观各向同性平面内梯度分布的压电体的状态方程，进而对材料系数指数函数规律分布的半无限大压电体中的反平面裂纹问题进行了求解，利用Fourier变换给出了半无限大压电体中位移，应力，电势及电位移的解析表达式，并求得了裂纹尖端的应力强度因子和电位移强度因子，分析了不同的非均匀材料系数及几何尺寸对它们的影响。     

含不可渗透边裂纹压电材料反平面问题的解

研究了含边缘裂纹的矩形截面压电材料在平面内电场和反平面荷载作用下的问题。得到了满足拉普拉斯方程、裂纹面边界条件的位移函数解和电势函数解及电弹场的基本解。最后,用边界配置法计算了能量释放率。本文提出的这种半解析半数值的方法计算简便,而且具有广泛的应用性。     

Crack problem for superconducting strip with finite thickness

The inclined crack problems are considered for a thin strip and a strip with finite thickness in a perpendicular magnetic field. The critical current density is assumed to be a constant. The crack orientation is varied and the effect of crack on the magnetic field distribution is neglected. Based on the analytical results and variational inequality, the field and current distributions are computed for both thin strip and strip with finite thickness cases, respectively. Then, the stress intensity factors at the crack tip are determined using the finite element method for magnetic field loads. The numerical results are presented for different inclined crack angles, magnetization processes and geometry parameters of the strip. The results show that the fracture behavior of the strip with finite thickness is more complicated than that of the thin strip. With the numerical results, we can predict the largest possibility of cracking as the strip is in an external field.     

Mode III fracture of an arbitrary oriented crack in two dimensional functionally graded material

In this paper, a two dimensional functionally graded material (2D-FGM) under an anti-plane load with an internal crack is considered. The crack is oriented in an arbitrary direction. The material properties are assumed to vary exponentially in two planar directions. The problem is analyzed and solved by two different methods namely Fourier integral transforms with singular integral equation technique, and then by the finite element method. The effects of crack orientation, material non-homogeneity, and other parameters on the value of stress intensity factor (SIF) are studied. Finally, the obtained results for Mode III stress intensity factor of different methods are compared.     

Dynamic fracture mechanics analysis for composite material with material non-homogeneity in thickness direction

The problem considered here is the response of a non-homogeneous composite material containing some cracks subjected to dynamic loading. It is assumed that the composite material is orthotropic and all the material properties depend only on the coordinatey (along the thickness direction). In the analysis, the elastic region is divided into a number of plies of infinite length. The material properties are taken to be constants for each ply. By utilizing the Laplace transform and Fourier transform technique, the general solutions for plies are derived. The singular integral equations of the entire elastic region are obtained and solved by the virtual displacement principle. Attention is focused on the time-dependent full field solutions of stress intensity factor(SIF) and strain energy release rate. As a numerical illustration, the dynamic stress intensity factor of a substrate/functionally graded film structure with two cracks under suddenly applied forces on cracks face are presented for various material non-homogeneity parameters.     

Antiplane shear fields near a crack-tip in a brittle damaged material

In this paper, a simplified brittle damage model is proposed according to the Mazarz-Lemaitre damage model for concrete. A closed-form solution for a mode III crack is obtained based on the simplified model under small scale damage conditions, which allows for discontinuities of displacement-gradient and tangential stress on the damage boundary. It is pointed out that the discontinuities of field variables near the tip region exist for the brittle damaged material induced by the softening effect of the material. The preoject supported by the National Natural Science Foundation of China     

Elastic stress distributions for hyperbolic and parabolic notches in round shafts under torsion and uniform antiplane shear loadings

Closed-form solutions are developed for the stress fields induced by circumferential hyperbolic and parabolic notches in axisymmetric shafts under torsion and uniform antiplane shear loading. The boundary value problem is formulated by using complex potential functions and two different coordinate systems, providing two classes of solutions. It is also demonstrated that some solutions of linear elastic fracture and notch mechanics reported in the literature can be derived as special cases of the general solutions proposed herein.Finally the analytical frame is used to link the Mode III notch stress intensity factor to the maximum shear stress at the notch tip, as well as to give closed-form expressions for the strain energy averaged over a finite size volume surrounding the notch root.     

Sudden deceleration or acceleration of an intersonic shear crack

Sudden jumps in the crack tip velocity were revealed by numerical simulation (in both continuum/cohesive element and molecular dynamics approaches) and experiments for rapid shear cracking. The cracking velocity may accelerate from a sub-Rayleigh speed to the intersonic range, or from an intersonic speed to a higher one, when the reflected impact wave reloads the crack tip. On the other hand, the cracking velocity may decelerate from an intersonic speed to a lower one or recede to the sub-Rayleigh range when the fracture driving force declines. The velocity change encountered during intersonic cracking plays a different role from that in the acceleration or deceleration of a subsonic crack. A crack propagating at an intersonic speed would leave a shear wave trailing behind. When the crack decelerates or accelerates, the effect of the trailing wave will lead to a transition period from one steady-state solution of crack tip singularity to another. This investigation aims at quantifying these processes. The full field solution of an intersonic mode II crack whose speed changed suddenly from one velocity (intersonic or subsonic) to another (intersonic or subsonic) is given in closed form. The solution is facilitated via superposing a series of propagating crack problems that are loaded by dislocations to seal the unwanted crack-face sliding or by concentrated forces moving at various speeds to negate the crack-face traction. In contrast to the subsonic solution, the results in the intersonic case indicate that the elastic fields around the crack tip depend on the deceleration or acceleration history that is traced back over a long time. Singularity matching dictates the jump that may actually take place.     

A rigid line in a confocal elliptic inhomogeneity embedded in an infinite medium under antiplane shear

An effective method is developed and used to investigate the antiplane problem of a rigid line in a confocal elliptic inhomogeneity embedded in an infinite medium. The analytical solution is obtained. The proposed method is based upon the use of conformal mapping and the theorem of analytic continuation. Special solutions which are verified by comparison with existing ones are provided. Finally, the characteristics of stress singularity at the tip of the rigid line inhomogeneity are analyzed and the extension forces for the crack and the rigid line inhomogeneity are derived. The project supported by the National Natural Science Foundation of China, the State Education Commission Foundation and the Failure Mechanics Lab of SEdC.     

Experimental caustics analysis in fracture mechanics of anisotropic materials

The originally developed reflection method of caustics is presented for application to cracks in mechanically anisotropic materials such as fiber-reinforced composites. The derived solutions for the combination of crack-opening modes I and II show that the size of the dark spot depends on the load intensity, whereas its shape depends strongly on the mechanical properties of the material, the orientation of the crack tip and the mixed-mode mixture. The evaluation of optical effects is possible using the diameter-measuring method or the advanced multipoint overdeterministic data reduction method. To find the exact position of caustics, the experimental images are analyzed by the simple boundary value method and a more sophisticated differential method, which is accomplished by shifting the real image onto the gradient image. The standard experimental testing procedure is performed for cracks oriented 0 deg, 45 deg and 90 deg to the material axes in carbon-fiber-reinforced polymer composites.     

ANTI-PLANE FRACTURE ANALYSIS OF FUNCTIONALLY GRADIENT MATERIAL INFINITE STRIP WITH FINITE WIDTH

The special case of a crack under mode III conditions was treated, lying parallel to the edges of an infinite strip with finite width and with the shear modulus varying exponentially perpendicular to the edges. By using Fourier transforms the problem was formulated in terms of a singular integral equation. It was numerically solved by representing the unknown dislocation density by a truncated series of Chebyshev polynomials leading to a linear system of equations. The stress intensity factor (SIF) results were discussed with respect to the influences of different geometric parameters and the strength of the non-homogeneity. It was indicated that the SIF increases with the increase of the crack length and decreases with the increase of the rigidity of the material in the vicinity of crack. The SIF of narrow strip is very sensitive to the change of the non-homogeneity parameter and its variation is complicated. With the increase of the non-homogeneity parameter, the stress intensity factor may increase, decrease or keep constant, which is mainly determined by the strip width and the relative crack location. If the crack is located at the midline of the strip or if the strip is wide, the stress intensity factor is not sensitive to the material non-homogeneity parameter.     

Crack driving force and energy-momentum tensor in electroelastodynamic fracture

The energy flux integral and the energy-momentum tensor for studying the crack driving force in electroelastodynamic fracture are formulated within the framework of the nonlinear theory of coupled electric, thermal and mechanical fields based on fundamental principles of thermodynamics. This formulation lays a foundation for in-depth understanding of the fracture behavior of piezoelectric materials. Remarkably, the dynamic energy release rate thus obtained has an odd dependence on the electric displacement intensity factor for steady-state propagation of a conventional (unelectroded) crack with exact, electrically permeable, semi-permeable, or impermeable crack surface condition, which is in agreement with experimental evidence.