Dynamic stress intensity factors of two collinear mode-III cracks perpendicular to and on the two sides of a bi-FGM weak-discontinuous interface

The mechanical model was established for the anti-plane dynamic fracture problem for two collinear cracks on the two sides of and perpendicular to a weak-discontinuous interface between two materials with smoothly graded elastic properties, as opposed to a sharp interface with discontinuously changing elastic properties. The problem was reduced as a system of Cauchy singular integral equations of the first kind by Laplace and Fourier integral transforms. The integral equations were solved by Erdogan's collocation method and the dynamic stress intensity factors in the time domain were obtained through Laplace numerical inversion proposed by Miller and Guy. The influences of geometrical and physical parameters on the dynamic stress intensity factors were illustrated and discussed, based on which some conclusions were drawn: (a) to increase the thickness of the FGM strip on either side of the interface will be beneficial to reducing the DSIF of a crack perpendicular to a bi-FGM interface and embedded at the center of one of the FGM strips; (b) To increase the rigidity of the FGM strip where the crack is located will increase the DSIF. However, when the material in one side of the interface is more rigid, the DSIF of the interface-perpendicular embedded crack in the other side will be reduced; (c) To decrease the weak-discontinuity of a bi-FGM interface will not necessarily reduce the stress intensity factor of a crack perpendicular to it, which is different from the case of interfacial crack; (d) For two collinear cracks with equal half-length, when the distance between the two inner tips is less than about three times of the half-length, the interaction of them is intensified, however, when the distance is greater than this the interaction becomes weak.     

Impact response of a piezoelectric layered composite plate with a crack

The dynamic response of a central crack in a piezoelectric layered composite plate under normal impact is analyzed. The crack is oriented normally to the interfaces. The Laplace and Fourier transform techniques are used to formulate the problem in terms of a singular integral equation. The order of stress singularity around the tip of the terminated crack is also obtained. Numerical calculations are carried out, and the main results presented are the variations of the dynamic stress intensity factor and the dynamic energy density factor versus time as functions of the geometric parameters and the piezoelectric material properties of the layered composite plate.     

Anti-plane dynamic hole–crack interaction in a functionally graded piezoelectric media

The anti-plane dynamic problem of a functionally graded piezoelectric plane containing a hole–crack system is treated by a non-hypersingular traction-based boundary integral equation method. The material parameters vary exponentially in the same manner in an arbitrary direction. The system is loaded by an incident SH-type wave, and impermeable boundary conditions are assumed. Using a frequency-dependent fundamental solution of the wave equation, the boundary value problem is transformed into a system of integro-differential equations along the boundary of the hole and on the crack line. Its numerical solution yields the dynamic stress intensity factors and stress concentration factors. A parametric study reveals their dependence on the hole–crack scenario and its geometry, characteristics of the dynamic load and magnitude and direction of material inhomogeneity.     

Dynamic analysis of a penny-shaped dielectric crack in a magnetoelectroelastic solid under impacts

The transient response of a magneto-electro-elastic material with a penny-shaped dielectric crack subjected to in-plane magneto-electro-mechanical impacts is made. To simulate an opening crack with a dielectric interior, the crack-face electromagnetic boundary conditions are supposed to depend on the crack opening displacement and the jumps of electric and magnetic potentials across the crack. Four ideal crack-face electromagnetic boundary conditions involving a combination of electrically permeable or impermeable and magnetically permeable or impermeable assumptions can be reduced. The Laplace and Hankel transform techniques are further utilized to solve the mixed initial-boundary-value problem. Three coupling Fredholm integral equations are obtained and solved by the composite Simpson's rule. Dynamic field intensity factors of stress, electric displacement, magnetic induction, crack opening displacement (COD), electric potential and magnetic potential are given in the Laplace transform domain. By means of a numerical inversion of the Laplace transform, numerical results are calculated to show the variations of the physical parameters of concern versus the normalized time in graphics. The effects of applied electric and magnetic loads on the dynamic intensity factors of stress and COD, and the dynamic energy release rate for a BaTiO₃-CoFe₂O₄ composite with a penny-shaped vacuum crack are discussed in detail.     

DYNAMIC STRESS FIELD AROUND THE MODE Ⅲ CRACK TIP IN AN ORTHOTROPIC FUNCTIONALLY GRADED MATERIAL

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ｇｒｅａｔａｔｔｅｎｔｉｏｎｓｈａｖｅｂｅｅｎｐａｉｄｔｏｔｈｅｒｅｓｅａｒｃｈｏｆＦｕｎｃｔｉｏｎａｌｌｙＧｒａｄｅｄＭａｔｅｒｉａｌｓ(ＦＧＭ).Ｆｒｏｍｔｈｅｖｉｅｗｐｏｉｎｔｓｏｆａｐｐｌｉｅｄｍｅｃｈａｎｉｃｓ,ＦＧＭａｒｅｎｏｎ＿ｈｏｍｏｇｅｎｅｏｕｓｓｏｌｉｄｓ.Ｔｈｅｎｏｎ＿ｈｏｍｏｇｅｎｅｉｔｙｏｆＦＧＭｈａｓａｇｒｅａｔｉｎｆｌｕｅｎｃｅｏｎｔｈｅｉｒｍｅｃｈａｎｉｃａｌｂｅｈａｖｉｏｒ,ｅｓｐｅｃｉａｌｌｙｗｈｅｎｔｈｅｃｏｍｐ…     

The mixed-mode fracture mechanics analysis of an embedded arbitrary oriented crack in a two-dimensional functionally graded material plate

Mixed-mode fracture mechanics analysis of an embedded arbitrarily oriented crack in a two-dimensional functionally graded material using plane elasticity theory is considered. The material properties are assumed to vary exponentially in two planar directions. Then, employing Fourier integral transforms with singular integral equation technique, the problem is solved. The stress intensity factors (SIFs) at the crack tips are calculated under in-plane mechanical loads. Finally, the effects of crack orientation, material non-homogeneity, and other parameters are discussed on the value of SIF in mode I and mode II fracture.     

Dynamic analysis of a cracked magnetoelectroelastic medium under antiplane mechanical and inplane electric and magnetic impacts

The transient analysis of a magnetoelectroelastic medium containing a crack is made under antiplane mechanical and inplane electric and magnetic impacts. The crack is assumed to penetrate through the solid along the poling direction. By using the Fourier and Laplace transforms, the associated mixed boundary value problem is reduced to a Fredholm integral equation of the second kind, which is solved numerically. By means of a numerical inversion of the Laplace transform, dynamic field intensity factors are obtained in the time domain. Numerical results are presented graphically to show the effects of the material properties and applied electric and magnetic impacts on the dynamic intensity factors of COD and stress, and dynamic energy density factors. The results indicate that except for the intensity factors of electric displacement and magnetic induction, other field intensity factors exhibit apparent transient feature. Moreover, they depend strongly on mechanical input as well as electric and magnetic impacts.     

SCATTERING OF THE HARMONIC STRESS WAVE BY CRACKS IN FUNCTIONALLY GRADED PIEZOELECTRIC MATERIALS

The present paper considers the scattering of the time harmonic stress wave by a single crack and two collinear cracks in functionally graded piezoelectric material （FGPM）. It is assumed that the properties of the FGPM vary continuously as an exponential function. By using the Fourier transform and defining the jumps of displacements 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 polynomials. Numerical examples are provided to show the influences of material properties on the dynamic stress and the electric displacement intensity factors.     

Dynamic crack tip fields and dynamic crack propagation characteristics of anisotropic material

Based on mechanics of anisotropic material, the dynamic crack propagation problem of I/II mixed mode crack in an infinite anisotropic body is investigated. Expressions of dynamic stress intensity factors for modes I and II crack are obtained. Components of dynamic stress and dynamic displacements around the crack tip are derived. The strain energy density theory is used to predict the dynamic crack extension angle. The critical strain energy density is determined by the strength parameters of anisotropic materials. The obtained dynamic crack tip fields are unified and applicable to the analysis of the crack tip fields of anisotropic material, orthotropic material and isotropic material under dynamic or static load. The obtained results show Crack propagation characteristics are represented by the mechanical properties of anisotropic material, i.e., crack propagation velocity M and fiber direction α. In particular, the fiber direction α and the crack propagation velocity M give greater influence on the variations of the stress fields and displacement fields. Fracture angle is found to depend not only on the crack propagation but also on the anisotropic character of the material.     

Antiplane mechanical and inplane electric time-dependent load applied to two coplanar cracks in piezoelectric ceramic material

This work is concerned with the dynamic response of two coplanar cracks in a piezoelectric ceramic under antiplane mechanical and inplane electric time-dependent load. The cracks are assumed to act either as an insulator or as a conductor. Laplace and Fourier transforms are used to reduce the mixed boundary value problems to Cauchy-type singular integral equations in Laplace transform domain. A numerical Laplace inversion algorithm is used to determine the dynamic stress and electric displacement factors that depend on time and geometry. A normalized equivalent parameter describing the ratio of the equivalent magnitude of electric load to that of mechanical load is introduced in the numerical computation of the dynamic stress intensity factor (DSIF) which has a similar trend as that for the pure elastic material. The results show that the dynamic electric field will impede or enhance crack propagation in a piezoelectric ceramic material at different stages of the dynamic electromechanical load. Moreover, the electromechanical response is greatly affected by the ratio of the crack length to the ligament between the cracks. The stress and electric displacement intensity factor can be combined by the energy density factor or function to address the fracture of piezoelectric materials under the combined influence of electromechanical loading.     

Analysis of a mode-I crack perpendicular to an imperfect interface

The elastostatic problem of a mode-I crack embedded in a bimaterial with an imperfect interface is investigated. The crack is in proximity to and perpendicular to the imperfect interface, which is governed by linear spring-like relations. The Fourier transform is applied to reduce the associated mixed-boundary value problem to a singular integral equation with Cauchy kernel. By numerically solving the resulting equation, stress intensity factors near both crack tips are evaluated. Obtained results reveal that the stress intensity factors in the presence of the imperfect interface vary between that with a perfect interface and that with a completely debonding interface. Moreover, an increase in the interface parameters decreases the stress intensity factors. In particular, when crack approaches to the weakened interface closer, the stress intensity factors become larger for a sliding interface, and become larger or smaller for a Winkler interface, depending on the crack lying in a stiffer or softer material. The influences of the imperfection of the interface on the stress intensity factors for a bimaterial composed of aluminum and steel are presented graphically.     

复合材料桥连的断裂动力学模型

复合材料产生裂纹后，其纤维处形成“桥连”，这是一个不可避免的现象。由于桥连问题很复杂．在数学方法的处理上有很大困难，至今人们研究大多是桥连的静力学问题．而对其动力学问题研究得很少。为了便于分析复合材料的问题，将桥连处用载荷代替，当裂纹高速扩展时．其纤维也连续地断裂。只有建立复合材料的桥连动力学模型，才能更好地研究复合材料的断裂动力学问题。通过复变函数论的方法，将所讨论的问题转化为Riemann—Hilbert问题。利用建立的动态模型和自相似方法，得到了正交异性体中扩展裂纹受运动的集中力P及阶跃载荷作用下位移、应力和动态应力强度因子的解析解，并通过叠加原理，最终求得了该模型的解。     

Effects of material gradients on transient dynamic mode-III stress intensity factors in a FGM

This paper presents a transient dynamic crack analysis for a functionally graded material (FGM) by using a hypersingular time-domain boundary integral equation method. The spatial variations of the material parameters of the FGM are described by an exponential law. A numerical solution procedure is developed for solving the hypersingular time-domain traction BIE. To avoid the use of time-dependent Green’s functions which are not available for general FGM, a convolution quadrature formula is adopted for approximating the temporal convolution, while a Galerkin method is applied for the spatial discretization of the hypersingular time-domain traction BIE. Numerical results for the transient dynamic stress intensity factors for a finite crack in an infinite and linear elastic FGM subjected to an impact anti-plane crack-face loading are presented and discussed. The effects of the material gradients of the FGM on the transient dynamic stress intensity factors and their dynamic overshoot over the corresponding static stress intensity factors are analyzed.     

Transient anti-plane crack problem for two bonded functionally graded piezoelectric materials

In this paper, the dynamic anti-plane crack problem for two bonded functionally graded piezoelectric materials is considered. The crack is perpendicular to the interface and assumed to be electrically impermeable or permeable. Integral transforms are employed to reduce the problem to Cauchy singular equations that can be solved numerically. The effects of the loading parameter λ, material constants and the geometry parameters on the stress intensity factor and the energy density factor are studied. It is found that for the impermeable case, the normalized dynamic stress intensity factor may increase or decrease in different time domains determined by the sign and magnitude of λ.     

An axisymmetric problem of an embedded crack in a graded layer bonded to a homogeneous half-space

In an attempt to simulate non-uniform coating delamination, the elasto-static problem of a penny shaped axisymmetric crack embedded in a functionally graded coating bonded to a homogeneous substrate subjected to crack surface tractions is considered. The coating’s material gradient is parallel to the axisymmetric direction and is orthogonal to the crack plane. The graded coating is modeled as a non-homogeneous medium with an isotropic constitutive law. Using Hankel transform, the governing equations are converted into coupled singular integral equations, which are solved numerically to yield the crack tip stress intensity factors. The Finite Element Method was additionally used to model the crack problem. The main objective of this paper is to study the influence of the material non-homogeneity and the crack position on the stress intensity factors for the purpose of gaining better understanding on the behavior of graded coatings.     

Dynamic stress intensity factors of multiple cracks in a functionally graded orthotropic half-plane

In the present paper dynamic stress intensity factor and strain energy density factor of multiple cracks in the functionally graded orthotropic half-plane under time-harmonic loading are investigated. By utilizing the Fourier transformation technique the stress fields are obtained for a functionally graded orthotropic half-plane containing a Volterra screw dislocation. The variations of the material properties are assumed to be exponential forms which the equilibrium has an analytical solution. The dislocation solution is utilized to formulate integral equation for the half-plane weakened by multiple smooth cracks under anti-plane deformation. The integral equations are of Cauchy singular type at the location of dislocation which are solved numerically to obtain the dislocation density on the faces of the cracks. The dislocation densities are employed to determined stress intensity factor and strain energy density factors (SEDFs) for multiple smooth cracks under anti-plane deformation. Numerical examples are provided to show the effects of material properties and the crack configuration on the dynamic stress intensity factors and SEDFs of the functionally graded orthotropic half-plane with multiple curved cracks.     

Transient response of a piezoelectric ceramic strip with an eccentric crack under electromechanical impacts

The transient response of a piezoelectric strip with an eccentric crack normal to the strip boundaries under applied electromechanical impacts is considered. By using the Laplace transform, the mixed initial-boundary-value problem is reduced to triple series equations, then to a singular integral equation of the first kind by introducing an auxiliary function. The Lobatto–Chebyshev collocation technique is adopted to solve numerically the resulting singular integral equation. Dynamic field intensity factors and energy release rate are obtained for both a permeable crack and an impermeable crack. The effects of the crack position and the material properties on the dynamic stress intensity factor are examined and numerical results are presented graphically.     

Electromechanical impact of an impermeable parallel crack in a functionally graded piezoelectric strip

The dynamic fracture problem for a functionally graded piezoelectric material (FGPM) strip containing a crack parallel to the free boundaries is considered in this study. It is assumed that the electroelastic properties of the strip vary continuously along the thickness direction of the strip, and that the strip is under the in-plane mechanical and electric impact. Integral transform techniques and dislocation density functions are employed to reduce the problem to the solutions of a system of singular integral equations. The dynamic stress and electric displacement intensity factors versus time are presented for various values of dimensionless parameters representing the crack size, the crack location, the material nonhomogeneity and the loading combination.     

Finite element analysis for steady-state hydride-induced fracture in metals by composite model

Delayed hydride cracking, which is observed in hydride-forming metals, due to the precipitation of hydrides near the crack tip, is investigated under conditions of constant temperature and crack velocity, plane strain and small-scale hydride-precipitation. The coupling of the operating physical processes of hydrogen-diffusion, hydride precipitation and material deformation is taken into account. The material is assumed to be an elastic composite made of hydrides and solid solution, with properties depending locally on the volume fraction of the hydrides. In the present analysis, the composite elastic properties have been derived by a generalized self consistent model for particulate composites. With respect to hydride-precipitation, two cases have been considered: (i) precipitation in a homogeneous medium with elastic properties, equal to the effective properties of the composite and (ii) precipitation in an inhomogeneous medium, where the expanding hydride has different elastic properties than those of the surrounding solid solution. The differences between the near-tip field distributions, produced by the two precipitation models, are relatively small. The effect of the hydrogen concentration far from the crack tip, on the near-tip field is also studied. It is shown that for small crack growth velocities, near the threshold stress intensity factor, the remote hydrogen concentration weakly affects the normalized stress distribution in the hydride-precipitation zone, which is controlled by the thermodynamically required hydrostatic stress, under hydrogen chemical equilibrium. However, for values of the applied stress intensity factor and the crack tip velocity, away from the threshold stress intensity factor and crack arrest, the effect of remote hydrogen concentration on the normalized near-tip stress field is strong. Reduction of the remote hydrogen concentration generally leads to reduction of the hydride-precipitation zone and increase of the near-tip stresses. Also reduction of the remote hydrogen concentration leads to distributions closer to those under hydrogen chemical equilibrium.     

用光弹性法测双层复合板的应力强度因子

本文用光弹性法分析了双层复合材料结构的裂纹问题，研究了两种不同材料的存在对应力强度因子的影响，以及讨论了由于裂纹不断加深，应力强度因子的变化规律，并与计算的结果作了比较，在裂纹不接近界面情况下，两者结果吻合的较好。