On Some Nonclassical Problems of Fracture Mechanics Taking into Account the Stresses Along Cracks

The paper reports on the results from analysis of some nonclassical fracture problems that allow for stresses along cracks. The results have been obtained by the author and his collaborators at the Institute of Mechanics (Kiev) since 1980     

A Self-Similar Crack Expansion method for three-dimensional cracks in an infinite or semi-infinite medium

The Self-Similar Crack Expansion (SSCE) method is used to calculate stress intensity factors for three-dimensional cracks in an infinite medium or semi-infinite medium by the boundary integral element technique, whereby, the stress intensity factors at crack tips are determined by calculating the crack-opening displacements over the crack surface. For elements on the crack surface, regular integrals and singular integrals are precisely evaluated based on closed form expressions, which improves the accuracy. Examples show that this method yields very accurate results for stress intensity factors of penny-shaped cracks and elliptical cracks in the full space, with errors of less than 1% as compared with analytical solutions. The stress intensity factors of subsurface cracks are in good agreement with other analytical solutions.     

Finite element approximation of the Neumann eigenvalue problem in domains with multiple cracks

^** Email: belhach{at}poncelet.univ-metz.fr^*** Email: bucur{at}math.univ-metz.fr^**** Email: jmse{at}math.univ-metz.fr We study the Neumann–Laplacian eigenvalue problem in domainswith multiple cracks. We derive a mixed variational formulationwhich holds on the whole geometric domain (including the cracks)and implements efficient finite-element discretizations forthe computation of eigenvalues. Optimal error estimates aregiven and several numerical examples are presented, confirmingthe efficiency of the method. As applications, we numericallyinvestigate the behaviour of the low eigenvalues in domainswith a large number of cracks.     

一维六方准晶中带三条不对称裂纹的圆形孔口问题的解析解

利用复变函数方法,通过构造保角映射,研究了一维六方准晶中带不对称三裂纹的圆形孔口的反平面剪切问题,给出了Ⅲ型裂纹问题的应力强度因子,在极限情形下,不仅可以还原为已有的结果,而且求得一维六方准晶中L裂纹问题在裂纹尖端的应力强度因子.     

Transient response of two collinear dielectric cracks in a piezoelectric solid under inplane impacts

The paper is focused on the dynamic analysis of two collinear dielectric cracks in a piezoelectric material under the action of in-plane electromechanical impacts. Considering the dielectric permeability of crack interior, the electric displacements at the crack surfaces are governed by the jumps of electric potential and crack opening displacement across the cracks. The permeable and impermeable crack models are the limiting cases of the general one. The Laplace and Fourier transform techniques are further utilized to solve the mixed initial-boundary-value problem, and then to obtain the singular integral equations with Cauchy kernel, which are solved numerically. Dynamic intensity factors of stress, electric displacement and crack opening displacement are determined in time domain by means of a numerical inversion of the Laplace transform. Numerical results for PZT-5H are calculated to show the effects of the dielectric permeability inside the cracks, applied electric loadings and the geometry of the cracks on the fracture parameters in graphics. The observations reveal that based on the COD intensity factor, a positive electric field enhances the dynamic dielectric crack growth and a negative one impedes the dynamic dielectric crack growth in a piezoelectric solid.     

Reconstruction of cracks of different types from far‐field measurements

In this paper, we deal with the acoustic inverse scattering problem for reconstructing cracks of possibly different types from the far‐field map. The scattering problem models the diffraction of waves by thin two‐sided cylindrical screens. The cracks are characterized by their shapes, the type of boundary conditions and the boundary coefficients (surface impedance). We give explicit formulas of the indicator function of the probe method, which can be used to reconstruct the shape of the cracks, distinguish their types of boundary conditions, the two faces of each of them and reconstruct the possible material coefficients on them by using the far‐field map. To test the validity of these formulas, we present some numerical implementations for a single crack, which show the efficiency of the proposed method for suitably distributed surface impedances. The difficulties for numerically recovering the properties of the crack in the concave side as well as near the tips are presented and some explanations are given. Copyright © 2009 John Wiley & Sons, Ltd.     

Basic solution of four parallel non-symmetric permeable Mode-III cracks in a piezoelectric/piezomagnetic composite plane

The interaction of four parallel non-symmetric permeable cracks in a piezoelectric/piezomagnetic composite plane subjected to anti-plane shear stress loading was studied by the Schmidt method. The problem was formulated through a Fourier transform into four 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 relationships among the electric displacement, magnetic flux and stress fields near the crack tips were 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 composites.     

Dislocation emission from a crack under mixed mode loading studied by molecular statics

The dislocation emission surface in (k _I,?k _II,?k _III) space is calculated by means of atomistic simulations for the {111}?110? crack in Al. For each relevant combination of loading mode, the precise nature of the dislocations and of the emission process are determined. When appropriate, the analytic formulas proposed by Rice are used by calculating the unstable stacking energy including the effect of the mixed mode loading. Quantitative agreement with the full atomistic calculation is found in the case where dislocations glide in the crack plane. This clearly identifies when and how ab initio data can be introduced in the calculation.     

Fatigue Behavior of Composite Sandwich Panels Under Three Point Bending Load

In this paper, the bending fatigue tests of honeycomb sandwich panels are carried out by using an improved three-point bending test fixture, and the S-N curves at different stress ratios are obtained. Through the records of fatigue damage in the experiment, the failure mode of the honeycomb sandwich panels and the source of fatigue damage are determined. At the same time, through the calculation of the shear stress distribution on the honeycomb wall, the reasons for the difference in the failure morphology of the L-direction and W-direction sandwich panels are clarified. Besides, a life prediction method is proposed and its effectiveness in predicting the fatigue life of sandwich panels has been verified.     

Parameter identification for cracks in elastic plate structures based on remote strain fields

A method for the detection of cracks in plate structures is presented. In contrast to most of the common monitoring concepts taking advantage of the reflection of elastic waves at crack faces, the presented approach is based on the strain measured at different locations on the surface of the structure. This allows both the identification of crack position parameters, such as length, location and angles with respect to a reference coordinate system and the calculation of stress intensity factors (SIF). The solution of the direct problem is performed on the basis of the BFM (body force method). The inverse problem is solved applying the particle swarm optimization (PSO) algorithm. The BFM is based on the principle of linear superposition which allows the calculation of the strain field in a cracked body. The strain at an arbitrary point in the structure is replaced by the strain provided by body force doublets in the uncracked structure. The doublets as well as external loads are parameters which have to be determined solving the inverse problem by minimizing a fitness function, which is defined by a square sum of residuals between measured strain distributions and computed ones for an assumed crack. The PSO algorithm applied to the fitness function operates on the basis of a swarm of candidate solutions. Once knowing loading and crack parameters, the SIF can be determined.