Periodic system of collinear cracks in an elastic porous medium

We consider a problem for cracks in a specific porous elastic material described by the Cowin-Goodman-Nunziato model. For an arbitrary set of collinear cracks, the problem is reduced to an integral equation over the crack surface. In the special case of a periodic crack arrangement, the kernel of the integral equation is written in the form of a Fourier series. By explicitly summing the leading part of the kernel, we can show that it exhibits hypersingular behavior as the argument tends to zero. Further, we solve the main integral equation by using the direct numerical method proposed in our preceding paper. In conclusion, we study the stress concentration at the crack tips and compare this case with the case of a single isolated crack.     

On dynamic stress analysis for cracks in elastic materials with voids

The present paper is concerned with the development of a semi-analytical approach to the dynamic problem of the concentration of stresses near the edges of a crack located in a porous elastic space (two-dimensional problem) and subjected to a normal oscillating load applied to the crack faces. Our analysis is made in the context of the Goodman–Cowin–Nunziato (G–C–N) theory for porous media. In previous work we studied static crack problems for such materials; now we introduce an analysis of the relevant dynamic aspects. By using the Fourier transform, the problem is reduced in explicit form to a hyper-singular integral equation with a convolution kernel valid over the crack length. Then, we apply a collocation technique developed in our previous work to solve this equation, and study the stress intensity factor. The principal goal is to compare the stress intensity factor for the static and dynamic cases. We also compare our results with the case of an ordinary linear elastic medium.     

Fracture analysis of mode-II crack perpendicular to imperfect bimaterial interface

The problem of a mode-II crack close to and perpendicular to an imperfect interface of two bonded dissimilar materials is investigated.The imperfect interface is modelled by a linear spring with the vanishing thickness.The Fourier transform is used to solve the boundary-value problem and to derive a singular integral equation with the Cauchy kernel.The stress intensity factors near the left and right crack tips are evaluated by numerically solving the resulting equation.Several special cases of the mode-II crack problem with an imperfect interface are studied in detail.The effects of the interfacial imperfection on the stress intensity factors for a bimaterial system of aluminum and steel are shown graphically.The obtained observation reveals that the stress intensity factors are dependent on the interface parameters and vary between those with a fully debonded interface and those with a perfect interface.     

Orthotropic semi-infinite medium with a crack under thermal shock

A cracked orthotropic semi-infinite plate under thermal shock is investigated. The thermal stresses are generated due to sudden cooling of the boundary by ramp function temperature change. The superposition technique is used to solve the problem. The crack problem is formulated by applying the thermal stresses obtained from the uncracked plate with opposite sign to be the only external loads on the crack surfaces as the crack surface tractions. The Fourier transform technique is used to solve the problem leading to a singular equation of the Cauchy type. The singular integral equation is solved numerically using the expansion method. The influence of the material orthotropy on the stress intensity factors is shown by comparing the results obtained for different orthotropic materials and isotropic materials in the case of plane stress. The numerical results of the stress intensity factors are demonstrated as a function of time, crack length, location of the crack and the duration of the cooling rate.     

Interfacial fracture analysis of bonded dissimilar strips with a functionally graded interlayer under antiplane deformation

This paper provides the solution to the problem of dissimilar, homogeneous semi-infinite strips bonded through a functionally graded interlayer and weakened by an embedded or edge interfacial crack. The bonded system is assumed to be under antiplane deformation, subjected to either traction-free or clamped boundary conditions along its bounding planes. Based on the Fourier integral transform, the problem is formulated in terms of a singular integral equation which has a simple Cauchy kernel for the embedded crack and a generalized Cauchy kernel for the edge crack. In the numerical results, the effects of geometric and material parameters of the bonded system on the crack-tip stress intensity factors are presented in order to quantify the interfacial fracture behavior in the presence of the graded interlayer.     

Modelling of dynamical crack propagation using time-domain boundary integral equations

We study dynamic antiplane cracks in the time domain by the boundary integral equation method (BIEM) based on the integral equation for displacement discontinuity (or crack opening displacement, COD) as a function of stress on the crack. This displacement discontinuity formulation presents the advantage, with respect to methods developed by Das and others in seismology, that it has to be solved only inside the crack. This BIEM is, however, difficult to implement numerically because of the hypersingularity of the kernel of the integral equation. Hence it is rewritten into a weakly singular form using a regularization technique proposed by Bonnet. The first step, following a method due to Sladek and Sladek, consists in converting the hypersingular integral equation for the displacement discontinuity into an integral equation for the displacement discontinuity and its tangential derivatives (dislocation density distribution); the latter involves a Cauchy type singular kernel. The second step is based on the observation that the hypersingularity is related to the static component of the kernel; the static singularity is then isolated and can be expressed in terms of weakly singular integrals using a result due to Bonnet. Although numerical applications discussed in this paper are all for the antiplane problem, the technique can be applied as well to in-plane crack dynamics.

The BIEM is implemented numerically using continuous linear space-time base functions to model the COD on the crack. In the present scheme the COD gradient interpolation is discontinuous at the element nodes while the integral equations are collocated at the element midpoints. This leads to an overdetermined discrete problem which is solved by standard least-squares methods. We use the dynamic BIEM to study a set of problems that appear in earthquake source dynamics, including the spontaneous dynamic crack propagation for a very simple rupture criterion. The numerical results compare favorably with the few exact solutions that are available. Then we demonstrate that difficulties experienced with finite difference simulations of spontaneous crack dynamics can be removed with the use of BIEM. The results are improved by the use of singular crack tip elements.     

Stress intensity factors for an oblique edge crack in a coating/substrate system with a graded interfacial zone under antiplane shear

This paper investigates the edge crack problem for a coating/substrate system with a functionally graded interfacial zone under the condition of antiplane deformation. With the interfacial zone being modeled by a nonhomogeneous interlayer having the continuously varying shear modulus between the dissimilar, homogeneous phases of the coated medium, the coating is assumed to contain an edge crack at an arbitrary angle to the interfacial zone. The Fourier integral transform method is used in conjunction with the coordinate transformations of basic field variables. Formulation of the proposed crack problem is then reduced to solving a singular integral equation with a generalized Cauchy kernel. The mode III stress intensity factors are defined and evaluated in terms of the solution to the integral equation. In the numerical results, the values of the stress intensity factors are plotted, illustrating the effects of the crack orientation angle for various material and geometric combinations of the coating/substrate system with the graded interfacial zone.     

Contact Problem for Porous Elastic Half-Plane

The paper is concerned with a static contact problem about a rigid punch on the free surface of a linear porous elastic half-plane. With the use of the Fourier transform the problem is reduced to a singular integral equation holding over the contact zone. This integral representation permits consideration of the Flamant problem (a line load on the half-plane) to be explicitly reduced to some quadratures. It is shown that in the classical linear elasticity limit the main integral equation has a Cauchy-type kernel, so distribution of the contact pressure is like in the Sadowsky punch-problem. For arbitrary porosity a numerical co-location technique is applied that allows one to analyze in detail the distribution of the contact pressure versus porosity. Both in the Flamant and Sadowsky problems we demonstrate a higher compliance of the porous foundation, with respect to the classical linear elastic results. This revised version was published online in July 2006 with corrections to the Cover Date.     

Solution of multiple crack problem in a finite plate using coupled integral equations

This paper studies a numerical solution of multiple crack problem in a finite plate using coupled integral equations. After using the principle of superposition, the multiple crack problem in a finite plate can be converted into two problems: (a) the multiple crack problem in an infinite plate and (b) a usual boundary value problem for the finite plate. For the former problem, the Fredholm integral equation is used. For the latter problem, a BIE based on complex variable is suggested in which a Cauchy singular kernel exists. For the proposed BIE, after using the inverse matrix technique, the dependence of the traction at a domain point from the boundary tractions is formulated indirectly. This is a particular advantage of the present study. Several numerical examples are provided and the computed results for stress intensity factor and T-stress at crack tips are given.     

Stress intensity factor for an edge crack in two bonded dissimilar materials under convective cooling

The transient thermal stress crack problem for two bonded dissimilar materials subjected to a convective cooling on the surface containing an edge crack perpendicular to the interface is considered. The problem is solved using the principle of superposition and the uncoupled quasi-static thermoelasticity. The crack problem is formulated by applying the transient thermal stresses obtained from the uncracked medium with opposite sign on the crack surfaces to be the only external loads. Fourier integral transform is used to solve the perturbation problem resulting in a singular integral equation of Cauchy type in which the derivative of the crack surface displacement is the unknown function. The numerical results of the stress intensity factors are calculated for both the edge crack and the crack terminating at the interface using two different composite materials and illustrated as a function of time, crack length, coefficient of heat transfer, and the thickness ratio.     

Crack problem for an inhomogeneous plane bonded by two different inhomogeneous half-planes

In this paper the crack problem.for two bonded inhomogeneous half-planes isconsidered.It is assumed that the different materials have the same Poisson ratio v.butgenerally speaking,both Young s moduli vary exponentially with the coordinate x indifferent form.Using the single crack solution of the inhomogeneous plane problem andFourier transform technique.the problem is reduced to a Cauchy-type singular integralequation.Several numerical examples to calculate the stress intensity factors are carriedout.     

A numerical approach for a crack problem by Gauss–Chebyshev quadrature

In this paper, a problem of a crack in an orthotropic strip is studied under plane strain conditions. It is assumed that normal displacements and shear stresses do not act on neither of the boundaries of the strip. Cauchy-type singular integral equation for the crack problem is derived by using the theory of plane elasticity and the Fourier transformation technique. A quadrature collocation approach is adopted for the numerical solutions of the singular integral equation. The effect of relative thickness and mechanical properties of strip on Mode I stress intensity factors (SIFs) are examined under different loading conditions. Some sample results are given for SIFs; also, material orthotropy and geometrical effects are discussed in detail.     

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.     

Electroelastic intensification near anti-plane shear crack in orthotropic piezoelectric ceramic strip

The theory of linear piezoelectricity is applied to solve the antiplane electroelastic problem of an orthotropic piezoelectric ceramic strip with a finite crack, which is situated symmetrically and oriented in a direction normal to the edges of the strip. Fourier transforms are used to reduce the problem to the solution of a pair of dual integral equations. They are then reduced to a Fredholm integral equation of the second kind. Numerical values on the stress intensity factor and the energy release rate for some piezoelectric ceramics are obtained and the results are displayed numerically to exhibit the electroelastic interactions.     

Interface crack between piezoelectric and elastic strips

Summary An interface crack between piezoelectric and elastic strips is analyzed using the theory of linear piezoelectricity. The combined out-of-plane mechanical and in-plane electrical loads are applied to the layered strip. Fourier transforms are used to reduce the problem to a pair of dual integral equations, which is then expressed in terms of a Fredholm integral equation of the second kind. The stress intensity factor is determined, and numerical analysis is performed and discussed. Received 22 September 1999; accepted for publication 3 May 2000     

Layered piezoelectric medium with interface crack under anti-plane shear

The theory of linear piezoelectricity is applied to solve the anti-plane shear problem of a piezoelectric layer sandwiched by two dissimilar homogeneous materials with a crack at the interface. Both mechanical and electrical loads are applied to the piezoelectric laminate. By the use of Fourier transforms, the mixed boundary value problem is reduced to a singular integral equation which is solved numerically to determine the stress intensity factors for several layered piezoelectric media, and the results are presented in graphical form.     

Anti-plane stress intensity, energy release and energy density at crack tips in a functionally graded strip with linearly varying properties

The fracture problem of a crack in a functionally graded strip with its properties varying in a linear form along the strip thickness under an anti-plane load is considered. The embedded anti-plane crack is located in the middle of strip half way through the thickness. The third mode stress intensity factor is derived using two different methods. In the first method, by employing Fourier integral transforms, the governing equation is converted to a singular integral equation, which is subsequently solved numerically by the collocation method based on Chebyshev polynomials. Then, the problem is solved by means of finite element method in which quadrilateral 8-node singular elements around each crack tip are used. After inspecting the validity of the solution technique, effects of crack geometry and non-homogeneous material parameter on the stress intensity, energy release and energy density are studied and the results of analytical and FEM solutions are compared.     

New integral equation for plane elasticity crack problems

In the usual formulation of singular equation approach for crack problems in plane elasticity [1,2], if one changes the right-hand term of the integral equation from tractions to resultant forces, a new integral equation can be obtained and is presented in this paper. The newly obtained integral equation has a log singular kernel. Interpolation equation for the dislocation functions (the undetermined functions in the integral equations) is proposed. Numerical examination is used to demonstrate the efficiency of the present technique, and a number of numerical examples are given.     

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.     

Cylindrical interface crack in a hollow layered functionally graded cylinder under static torsion

A hollow functionally graded composite cylinder under static torsion, which consists of an inner and outer elastic circular tube with a cylindrical interface crack, is studied in this work. By utilizing Fourier integral transform method, the mixed boundary value problem is reduced to a Cauchy singular integral equation, from which the numerical results of the stress intensity factor are obtained by the Lobatto–Chebyshev quadrature technique. Numerical results demonstrate the coupled effects of geometrical, physical, and functionally graded parameters on the interfacial fracture behavior.