Invariant Integrals for the Equilibrium Problem for a Plate with a Crack

We consider the equilibrium problem for a plate with a crack. The equilibrium of a plate is described by the biharmonic equation. Stress free boundary conditions are given on the crack faces. We introduce a perturbation of the domain in order to obtain an invariant Cherepanov–Rice-type integral which gives the energy release rate upon the quasistatic growth of a crack. We obtain a formula for the derivative of the energy functional with respect to the perturbation parameter which is useful in forecasting the development of a crack (for example, in study of local stability of a crack). The derivative of the energy functional is representable as an invariant integral along a sufficiently smooth closed contour. We construct some invariant integrals for the particular perturbations of a domain: translation of the whole cut and local translation along the cut.     

The interface crack in anisotropic bodies. stress singularities and invariant integrals

For the problem of the deformation of a composite anisotropic plate with a crack (in a linear formulation, with no assumption of symmetry), all possible power solutions are listed and general relations between the ordinary and singular solutions are revealed. The asymptotic form of the increment of the potential energy of deformation is computed for the cases of the rectilinear propagation of the crack, deviation of a shoot or branching. The form obtained for the final formula is the same as the classical version of the Griffiths formula and involves two invariant integrals. Two methods of determining the modes of radical singularities of the stress-strain state near the crack tip, associated with the use of force and energy criteria, are proposed.     

Diffraction of elastic waves by an inclined crack in a layer

For the problem of the diffraction of normal modes by an inclined crack in an elastic layer, an integral equation with explicit representation of the Fourier symbol kernel is derived in the form of the product of matrices. The algorithm for calculating the wave fields, based on the analytical representations obtained, enables a rapid parametric analysis to be carried out of the influence of the size and orientation of the crack on the transmission of travelling waves. The influence of the inclination of the crack on the effects of resonance trapping and localization of wave energy, which were established previously for the case of a horizontal crack, is analysed.     

An asymptotic form of the energy functional for an elastic body with a crack and a rigid inclusion. The plane problem

The plane problem in the linear theory of elasticity for a body with a rigid inclusion located within it is investigated. It is assumed that there is a crack on part of the boundary joining the inclusion and the matrix and complete bonding on the remaining part of the boundary. Zero displacements are specified on the outer boundary of the body. The crack surface is free from forces and the stress state in the body is determined by the bulk forces acting on it. The variation in the energy functional in the case of a variation in the rigid inclusion and the crack is investigated. The deviation of the solution of the perturbed problem from the solution of the initial problem is estimated. An expression is obtained for the derivative of the energy functional with respect to a zone perturbation parameter that depends on the solution of the initial problem and the form of the vector function defining the perturbation. Examples of the application of the results obtained are studied.     

Effect of shearing forces on the stressed state of a half space with crack

Using the method of boundary integral equations, we study the stressed state in the neighborhood of a plane crack perpendicular to the boundary of a half space. The crack surfaces are subjected to the action of shearing forces. The problem is reduced to two-dimensional hypersingular integral equations, and their regular kernels, taking into account interaction between the crack and boundary of the half space, are written in explicit form. The dependences of stress intensity factors on the angular coordinate are presented for different loads of the crack. Translated from Matematychni Metody ta Fizyko-Mekhanichni Polya, Vol. 51, No. 1, pp. 112–120, January–March, 2008.     

A singular perturbation approach to integral equations occurring in poroelasticity

Singular perturbation theory is used to solve the integral equationswhich occur when treating finite-length crack problems in porouselastic materials. The method provides the stress intensityfactors which characterize the near crack tip stress and displacementfields for small times. The method also gives the stress andpore pressure fields on the fracture plane for small times relativeto the diffusive time scale. In this paper, the authors treatcrack problems which are unmixed in the pore pressure boundarycondition on the fracture plane. The Abelian result that smalltimes correspond, in Laplace transform space, to large valuesof the transform variable is used to formulate the problemsin terms of a small parameter. Rescaling on this small parameterleads to inner problems which are eigensolutions of the semi-infiniteproblems treated earlier by the authors. The outer solutionsare given by elastic eigensolutions together with appropriatefluid dipole responses. These outer solutions give the completestress and pore pressure fields except in the neighbourhoodof the crack tips; in this region the outer solutions are asymptoticallymatched with inner solutions. The full outer solutions are givenhere as an asymptotic expansion for small times and enable thedevelopment of the outer fields to be followed in real time.A reciprocal theorem in Laplace transform space is used to checkthe small-time solutions. The inner problem is rescaled to asemi-infinite crack problem, so eigensolutions of this semi-infiniteproblem are used together with the known asymptotic behaviourof the real solution to identify the stress intensity factor.The stress intensity factor is then related to an integral involvingthe inner limit of the outer solution together with the eigensolutionof the semi-infinite problem. Using this integral, we recoverthe result for the stress intensity factor found using singularperturbation theory. A ‘nearly’ invariant integralanalogous to the invariant M integral used in elastostaticsis derived. Unfortunately, the poroelastic analogue is not invariant,although it is used to verify the small-time results.     

Problem of an elastic half-plane with a crack emerging orthogonally at the boundary

The problem of the half-plane, in which a finite crack emerges orthogonally at the boundary, is studied. On the edges of the crack a self-balancing load is applied. A detailed investigation is carried out for an integral equation with respect to the unknown derivative of the displacement jump, to which the problem can be reduced. The exact solution of the integral equation is constructed with the aid of the Mellin transform and the Riemann boundary value problem for the halfplane. The asymptotic behavior of the solution at both ends of the crack is elucidated. First the asymptotic behavior of the solution at the point of emergence of the crack is obtained and the dependence of this asymptotic behavior on the type of the load is established. For a special form of the load one obtains a simple expression of the stress intensity coefficient. In the case of a general load, the asymptotic behavior is used for the construction of an effective approximate solution on the basis of the method of orthogonal polynomials. As a result, the problem reduces to an infinite algebraic system, solvable by the reduction method.Translated from Dinamicheskie Sistemy, No. 4, pp. 45–51, 1985.     

External circular crack under arbitrary shear loading

An exact closed form solution in terms of elementary functions has been obtained to the governing integral equation of an external circular crack in a transversely isotropic elastic body. The crack is subjected to arbitrary tangential loading applied antisymmetrically to its faces. The recently discovered method of continuity solutions was used here. The solution to the governing integral equation gives the direct relationship between the tangential displacements of the crack faces and the applied loading. Now a complete solution to the problem, with formulae for the field of all stresses and displacements, is possible.     

Driving force acting on a running crack

The asymptotic analysis of the dynamic crack problem for the anti‐plane shear mode is provided. The field near the crack tip is studied in detail for a nonlinear elastic incompressible material whose stored energy behaves asymptotically as a power of the first invariant of the strain tensor at large strains. It is shown that the hardening parameter characterizes fully the singularity degree of the near‐crack‐tip field. Based on the latter knowledge the driving force acting on the crack tip is calculated. Possible scenarios of the crack propagation are discussed.     

On the boundary integral equations for the crack opening displacement of flat cracks

The boundary integral equations for the crack opening displacement in acoustic and elastic scattering problems are discussed in the case of flat cracks by means of the Fourier analysis technique. The pseudo-differential nature of the hypersingular integral operators is shown and their symbols explicited. It is then proved that the variational problems assocaited with these BIE are well-posed in a Sobolev functional framework which is closely linked with the elastic energy. A decomposition of the vector integral equation in the elastic case into scalar integral equations is obtained as a by-product of the variational formulation.     

Boundary value problem for the Laplace equation outside cuts on the plane with different conditions of the third kind on opposite sides of the cuts

We consider a boundary value problem for the Laplace equation outside cuts on a plane. Boundary conditions of the third kind, which are in general different on different sides of each cut, are posed on the cuts. We show that the classical solution of the problem exists and is unique. We obtain an integral representation for the solution of the problem in the form of potentials whose densities are found from a uniquely solvable system of Fredholm integral equations of the second kind.     

A normal crack in an elastic half-space with stress-free surface

The problem of an elastic half-space with stress-free surface and a crack of arbitrary shape with prescribed displacements or tractions is reduced to an equivalent system of integral equations on the crack. For a pressurized crack in a plane perpendicular to the free surface, a scalar integral equation is derived. In properly chosen function spaces, unique solvability of the integral equation and regularity of solutions for regular data are proven.     

A natural order in dynamical systems based on Conley–Markov matrices

We introduce a natural order to study properties of dynamical systems, especially their invariant sets. The new concept is based on the classical Conley index theory and transition probabilities among neighborhoods of different invariant sets when the dynamical systems are perturbed by white noises. The transition probabilities can be determined by the Fokker–Planck equation and they form a matrix called a Markov matrix. In the limiting case when the random perturbation is reduced to zero, the Markov matrix recovers the information given by the Conley connection matrix. The Markov matrix also produces a natural order from the least to the most stable invariant sets for general dynamical systems. In particular, it gives the order among the local extreme points if the dynamical system is a gradient-like flow of an energy functional. Consequently, the natural order can be used to determine the global minima for gradient-like systems. Some numerical examples are given to illustrate the Markov matrix and its properties.     

A Yoffe-type moving crack in one-dimensional hexagonal piezoelectric quasicrystals

A Yoffe-type moving crack in one-dimensional hexagonal piezoelectric quasicrystals is considered. The Fourier transform technique is used to solve a moving crack problem under the action of antiplane shear and inplane electric field. Full elastic stresses of phonon and phason fields and electric fields are derived for a crack running with constant speed in the periodic plane. Obtained results show that the coupled elastic fields inside piezoelectric quasicrystals depend on the speed of crack propagation, and exhibit the usual square-root singularity at the moving crack tip. Electric field and phason stresses do not have singularity and electric displacement and phonon stresses have the inverse square-root singularity at the crack tip for a permeable crack. The field intensity factors and energy release rates are obtained in closed form. The crack velocity does not affect the field intensity factors, but alters the dynamic energy release rate. Bifurcation angle of a moving crack in a 1D hexagonal piezoelectric quasicrystal is evaluated from the viewpoint of energy balance. Obtained results are helpful to better understanding crack advance in piezoelectric quasicrystals.     

Shape sensitivity of a plane crack front

The 3D‐elasticity model of a solid with a plane crack under the stress‐free boundary conditions at the crack is considered. We investigate variations of a solution and of energy functionals with respect to perturbations of the crack front in the plane. The corresponding expansions at least up to the second‐order terms are obtained. The strong derivatives of the solution are constructed as an iterative solution of the same elasticity problem with specified right‐hand sides. Using the expansion of the potential and surface energy, we consider an approximate quadratic form for local shape optimization of the crack front defined by the Griffith criterion. To specify its properties, a procedure of discrete optimization is proposed, which reduces to a matrix variational inequality. At least for a small load we prove its solvability and find a quasi‐static model of the crack growth depending on the loading parameter. Copyright © 2003 John Wiley & Sons, Ltd.     

An augmented galerkin procedure for the boundary integral method applied to two-dimensional screen and crack problems

Here we present a new solution procedure for Helm-holtz and Laplacian Dirichlet screen and crack problems in IR² via boundary integral equations of the first kind having as an unknown the jump of the normal derivative across the screen or a crack curve T. Under the assumption of local finite energy we show the equivalence of the integral equations and the original boundary value problem. Via the method of local Mellin transform in [5]-[lo] and the calculus of pseudodifferential operators we derive existence, uniqueness and regularity results for the solution of our boundary integral equations together with its explicit behaviour near the screen or crack tips.With our integral equations we set up a Galerkin scheme on T and obtain high quasi-optimal convergence rates by using special singular elements besides regular splines as test and trial functions.     

On the existence of invariant tori in nearly-integrable Hamiltonian systems with finitely differentiable perturbations

We prove the existence of invariant tori in Hamiltonian systems, which are analytic and integrable except a 2n-times continuously differentiable perturbation (n denotes the number of the degrees of freedom), provided that the moduli of continuity of the 2n-th partial derivatives of the perturbation satisfy a condition of finiteness (condition on an integral), which is more general than a Hölder condition. So far the existence of invariant tori could be proven only under the condition that the 2n-th partial derivatives of the perturbation are Hölder continuous.     

On fracture of viscoelastic bodies

On the basis of a thermodynamic variational principle, a criterion for quasistatic crack growth in a viscoelastic solid is deduced: the change in the sum of the scattering and the rate of decrease of elastic energy equals the increment in surface dissipation. The appropriate invariant contour integral is found. The theory proposed is suitable for cracks of any shape for any loading path on viscoelastic solids. Examples are considered, including crack growth with localized viscous dissipation.     

Dual boundary integral formulation for 2-D crack problems

An efficient integral equation formulation for two-dimensional crack problems is proposed with the displacement equation being used on the outer boundary and the traction equation being used on one of the crack faces. Discontinuous quarter point elements are used to correctly model the displacement in the vicinity of crack tips. Using this formulation a general crack problem can be solved in a single-region formulation, and only one of the crack faces needs to be discretised. Once the relative displacements of the cracks are solved numerically, physical quantities of interest, such as crack tip stress intensity factors can be easily obtained. Numerical examples are provided to demonstrate the accuracy and efficiency of the present formulation.     

Diffraction in a cylindrically orthotropic elastic solid containing a stress free crack

A cylindrically orthotropic elastic solid is excited by a point impulsive body force. The solid contains a semi-infinite stress free crack. The resulting anti-plane wave motion problem has been solved in the form of a finite series representing the incident and reflected pulses plus an integral representing the diffraction pulse. The series part of the solution has been previously treated. In the present investigation the diffraction integral is integrated when λ (which measures the anisotropy of the solid) is an odd integer number. The diffraction integral is also integrated when λ is half an odd integer, for the special case in which the source lies in the plane of the crack and parallel to the crack edge. The displacement jump across the circular diffraction wave front is given for unrestricted (positive) values of λ.