A crack of arbitrary shape inside an infinite transversely isotropic cylinder under arbitrary load

Summary  In the first part of the article an infinite circular cylinder is considered, made of transversely isotropic elastic material and weakened by a plane crack perpendicular to its axis O z. The crack is opened by an arbitrary normal stress. The second part is devoted to the same crack loaded by an arbitrary tangential stress. The complete solution in both cases is presented as a sum of the solution of a similar problem of a crack in an infinite space and an integral transform term, the parameters of which are determined from a set of linear algebraic equations derived from the boundary conditions. Governing integral equations with respect to the yet unknown crack displacement discontinuities are obtained. In the case of a circular crack, these equations can be inverted and solved by the method of consecutive interations. Received 30 November 2000; accepted for publication 3 May 2001     

Plastic zone correction in a stretched cylinder with an external circumferential crack

The Dugdale hypothesis is adapted to the problem of an external circumferential crack in a stretched cylinder. The lateral surface of the cylinder is stress free and restrained from radial displacements. An external circumferential edge crack in the cylinder which is considered elastic-perfectly plastic is envisaged with the assumption that the plastic zone forms a very thin in-plane layer surrounding the crack. The solution of the problem is reduced to the solution of dual Dini series which, in turn, is reduced to a Fredholm integral equation of the second kind. Solving this integral equation numerically and using the boundedness of the axial stress, the size of the plastic zone correction is obtained.     

Torsion of a cylinder with a shallow external crack

The torsional problem of a finite elastic cylinder with a circumferential edge crack is studied in this paper. An efficient solution to the problem is achieved by using a new form of regularization applied to dual Dini series equations. Unlike the Srivastav approach, this regularization transforms dual equations into a Fredholm integral equation of the second kind given on the crack surface. Hence, exact asymptotic expansions of the Fredholm equation solution, the stress intensity factor and the torque are derived for the case of a shallow crack. The asymptotic expansions are certain power-logarithmic series of the normalized crack depth. Coefficients of these series are found from recurrent relations. Calculations for a shallow crack manifest that the stress intensity factor exhibits the rather weak dependence upon the cylinder length when the torque is fixed and the triple length is larger than the diameter.     

Internally Pressurized Radially Polarized Piezoelectric Cylinders

The piezoelectric phenomenon has been exploited in science and engineering for decades. Recent advances in smart structures technology have lead to a resurgence of interest in piezoelectricity, and in particular, in the solution of fundamental boundary-value problems. In this paper, we develop an analytic solution to the axisymmetric problem of an infinitely long, radially polarized, radially orthotropic piezoelectric hollow circular cylinder. The cylinder is subjected to uniform internal pressure, or a constant potential difference between its inner and outer surfaces, or both. An analytic solution to the governing equilibrium equations (a coupled system of second-order ordinary differential equations) is obtained. On application of the boundary conditions, the problem is reduced to solving a system of linear algebraic equations. The stress distributions in the cylinder are obtained numerically for two typical piezoceramics of technological interest, namely PZT-4 and BaTiO₃. It is shown that the hoop stresses in a cylinder composed of these materials can be made virtually uniform throughout the cross-section by applying an appropriate set of boundary conditions. This revised version was published online in June 2006 with corrections to the Cover Date.     

A fast 3D dual boundary element method based on hierarchical matrices

In this paper a fast solver for three-dimensional BEM and DBEM is developed. The technique is based on the use of hierarchical matrices for the representation of the collocation matrix and uses a preconditioned GMRES for the solution of the algebraic system of equations. The preconditioner is built exploiting the hierarchical arithmetic and taking full advantage of the hierarchical format. Special algorithms are developed to deal with crack problems within the context of DBEM. The structure of DBEM matrices has been efficiently exploited and it has been demonstrated that, since the cracks form only small parts of the whole structure, the use of hierarchical matrices can be particularly advantageous. Test examples presented show that, with the proposed technique, substantial increase in number of elements over the crack surfaces leads only to moderate increases in memory storage and solution time.     

The effect of an elliptic crack on the stress distribution in a long circular cylinder

This paper deals with three-dimensional analysis of stress distribution in a long circular cylinder containing an elliptical crack. The surface of the crack is perpendicular to the axis of the cylinder, and is subjected to a constant pressure. The equations of the classical theory of elasticity are solved in terms of an unknown function which is then shown to be the solution of an integro-differential equation. Numerical solution of the integro-differential equation is obtained. The resultant stress intensity magnification factors for the elliptical crack in the cylinder are presented in graphical forms for various crack aspect ratios, and proximity ratios.     

Application of the angular superposition method to the contact problem on the compression of an elastic cylinder

The angular superposition method is used to construct an approximate solution of the contact problem on the compression of an elastic cylinder by two rigid plates. The solution thus obtained has a closed-form analytic expression and can be used in the entire domain of the cylinder cross-section. We analyze the absolute error, which takes the largest value near the points of contact between the plates and the cylinder, where the boundary conditions are discontinuous. According to the von Mises criterion, when moving into the depth of the cylinder from the contact site along the symmetry axis, the second invariant J ₂ of the stress deviator tensor first decreases and then, after attaining a minimum, increases and attains the largest value at a small depth, which agrees with Johnson’s photoelastic experiments and Dinnik’s computations. We present the graphs of the displacement and normal stress distributions over the contact site, the dependence of the compressing force on the displacements of rigid plates, and the dependence of the invariant J ₂ on the coordinate along the symmetry axis. If 640 computation points are chosen on the cylinder boundary and the Hertz law for the normal pressure on the contact site is used, then the error in the approximate solution near the endpoint of the contact site is approximately 55%, and if the proposed two-parameter normal law is used, then the error is of the order of 4%. On the free lateral surface of the cylinder boundary, we find the critical pointM*, which separates the cylinder contraction and extension parts.The contact problems are the most difficult problems, and their solution is complicated by the discontinuous boundary conditions [1–5]. In [6], the contact problem is solved by the Fourier method, which can be used only for bodies of classical shapes. In such cases, the problem can be reduced to solving coupled integral equations [7]. The interaction between the bandage and a cylindrical body is considered in [2, 6, 7]. In [8], the possibility of using the finite element method is investigated in the case of contact problems for a differential wheel with roughness of the contacting surfaces taken into account. In [9, 10], the method of homogeneous solutions is used to consider contact problems for a finite-dimensional elastic cylinder loaded on its end surfaces. Note that only error estimates are given in the literature cited above; the absolute error over the entire domain of the elastic body is not studied, although this is one of the important characteristics of the obtained approximate solution. A sufficiently complete survey of the literature in the field of contact interactions of elastic bodies is given in [3–5].In what follows, we propose to solve contact problems by the angular superposition method [11]. This method can be used for bodies of nonclassical shapes, which can be multiply connected, and the friction on the contact site can be taken into account. In the present paper, as a first example of applied character, we show how this method can be used in the simplest case. The multiple connectedness and the curvilinearity of the shape of the body, as well as taking into account the friction on the boundary, do not create new essential difficulties in this method.     

Transient thermal elastic fracture of a piezoelectric cylinder specimen

A finite piezoelectric cylinder with an embedded penny-shaped crack is investigated for a thermal shock load on the outer surface of the cylinder. The theory of linear electro-elasticity is applied to solve the transient temperature field and the associated thermal stresses and electrical displacements without crack. These thermal stresses and electrical displacements are added to the surfaces of the crack to form an electromechanical coupling and mixed mode boundary-value problem. The electrically permeable crack face boundary condition assumption is used, and the thermal stress intensity factor and electrical displacement intensity factor at the crack border are evaluated. The thermal shock resistance of the piezoelectric cylinder is evaluated for the analysis of piezoelectric material failure in practical engineering applications.     

The acoustic field in a rigid cylindrical vessel excited by a sphere oscillating by a definite law

The problem on the interaction between a spherical body oscillating by a definite law and a rigid cylinder filled with an ideal compressible liquid is formulated. The geometrical center of the sphere is located on the cylinder axis. The solution is based on the possibility of representing the particular solutions of the Helmholtz equation in cylindrical coordinates in terms of particular solutions in spherical coordinates, and vice versa. As a result of satisfaction of the boundary conditions on the surfaces of the sphere and cylinder, an infinite system of linear algebraic equations is obtained to determine the coefficients of expansion of the potential of liquid velocities into a Fourier series in terms of Legendre polynomials. The use of the reduction technique for solving the infinite system obtained is substantiated. The hydrodynamic characteristics of the liquid filling the cylindrical cavity are determined and compared with the case of a sphere vibrating in an infinite liquid. S. P. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine, Kiev. Translated from Prikladnaya Mekhanika, Vol. 36, No. 6, pp. 88–97, June, 2000.     

Stress intensity factors for a moving cylindrical crack in a nonhomogeneous cylindrical layer in composite materials

Stresses are determined for a finite cylindrical crack that is propagating with a constant velocity in a nonhomogeneous cylindrical elastic layer, sandwiched between an infinite elastic medium and a circular elastic cylinder made from another material. The Galilean transformation is employed to express the wave equations in terms of coordinates that are attached to the moving crack. An internal gas pressure is then applied to the crack surfaces. The solution is derived by dividing the nonhomogeneous interfacial layer into several homogeneous cylindrical layers with different material properties. The boundary conditions are reduced to two pairs of dual integral equations. These equations are solved by expanding the differences in the crack surface displacements into a series of functions that are equal to zero outside the crack. The Schmidt method is then used to solve for the unknown coefficients in the series. Numerical calculations for the stress intensity factors were performed for speeds and composite material combinations.     

三维有限体中平片裂纹的超奇异积分方程方法——Ⅰ、理论分析

本文利用单个平片裂纹的基本解,将三维有限体中的平片裂纹问题,归为解一组超奇异积分方程,然后使用主部分析方法,对这组方程的求解作了理论分析,其结果在本文的第Ⅰ部分给出,关于这组方程的数值法求解,则给出于本文的第Ⅱ部分。     

The line spring model with arbitrary loads on crack surfaces and its applications in thermal shock analysis

In this paper the line spring model taking account of arbitrary loads on crack surfaces, and the corresponding constitutive relations, are proposed. The general expressions of the additional outfield loads, which are equivalent to the distributed loads on crack surfaces, are derived. The model is used to compute stress intensity factors in a hollow cylinder with an axial surface crack subjected to thermal shock. Several results of calculations are presented and discussed.     

Interaction of crack-tip and notch-tip stress singularities for circular cylinder in torsion

Singular integral equations are used to formulate the torsion problem of a circular cylinder containing a polygonal opening and a line crack. The formulation is based on degenerating a system of connecting line cracks to that of a polygon, the sides of which coincides with the cracks. Considered, in particular, is the torsion of a circular cylinder with a rectangular hole and a nearby slanted line crack. Mode III stress intensity factors are computed at both ends of the crack to reflect their relative position to the rectangular hole in addition to change in the dimensions of the crack relative to the other geometric variables. Recognizing that the singular behavior of the stresses near a reentrant's corner differs from that of the crack tip, intensification of the local stresses at the corners of the rectangular hole is also examined. The results show the influence of the crack size and position.     

Axisymmetric displacement boundary value problem for a penny-shaped crack

Summary A solution is derived from equations of equilibrium in an infinite isotropic elastic solid containing a penny-shaped crack where displacements are given. Abel transforms of the second kind stress and displacement components at an arbitrary point of the solid are known in the literature in terms of jumps of stress and displacement components at a crack plane. Limiting values of these expressions at the crack plane together with the boundary conditions lead to Abel-type integral equations, which admit a closed form solution. Explicit expressions for stress and displacement components on the crack plane are obtained in terms of prescribed face displacements of crack surfaces. Some special cases of the crack surface shape functions have been given in the paper.     

Two classes of vibrations of a short circular cylinder

This paper gives a solution of the stationary dynamic problem of elasticity which describes two classes of natural nonaxisymmetric vibrations of a finite circular cylinder. In the particular case of axial symmetry, the resulting solution describes two well-known classes of axisymmetric vibrations: vibrations of the first class become longitudinal-transverse vibrations and vibrations of the second class become torsional vibrations. The existence of two classes of nonaxisymmetric vibrations is due to the boundary conditions at the ends. It is shown that as the length (height) of the cylinder increases, the effect of the boundary conditions at the ends on the frequency spectrum reduces, and the vibration frequencies of the two classes become similar and then identical.     

Axisymmetric circumferential internal crack problem of a thick-walled cylinder with inner and outer claddings

The elastostatic axisymmetric problem for a long thick-walled cylinder containing an axisymmetric circumferential internal crack with two claddings is considered. The claddings having different elastic properties than the hollow cylinder are assumed to be bonded to inner and outer wall of the hollow cylinder. The problem is formulated in terms of a singular integral equation of a well known type, the derivative of the crack surface displacement being the density function, using the standard transform technique. By using appropriate quadrature formulas, the integral equation is reduced to a system of linear algebraic equations. This system is solved numerically and the related stress-intensity factors are calculated for the cases of hollow cylinder with two claddings bonded to inner and outer wall of the cylinder, a cladding bonded to inner wall of the cylinder, a cladding bonded to outer wall of the cylinder and no cladding under axial tensile load. The influence of the geometrical configuration, the claddings and internal crack length on the stress-intensity factors is shown graphically.     

Cracked elastic layer under compressive mechanical loads

We consider boundary value problem in which an elastic layer containing a finite length crack is under compressive loading. The crack is parallel to the layer surfaces and the contact between crack surfaces are either frictionless or with adhesive friction or Coulomb friction.Based on fourier integral transformation techniques the solution of the formulated problems is reduced to the solution of a singular integral equation, then, using Chebyshev’s orthogonal polynomials, to an infinite system of linear algebraic equations. The regularity of these equations is established. The expressions for stress and displacement components in the elastic layer are presented. Based on the developed analytical algorithm, extensive numerical investigations have been conducted.The results of these investigations are illustrated graphically, exposing some novel qualitative and quantitative knowledge about the stress field in the cracked layer and their dependence on geometric and applied loading parameters. It can be seen from this study that the crack tip stress field has a mode II type singularity.     

Thermal stress intensity factors for a crack in an anisotropic half plane

On the basis of the two-dimensional theory of anisotropic thermoelasticity, a solution is given for the thermal stress intensity factors due to the obstruction of a uniform heat flux by an insulated line crack in a generally anisotropic half plane. The crack is replaced by continuous distributions of sources of temperature discontinuity and dislocations. First, the particular thermoelastic dislocation solutions for the half plane are obtained; then the corresponding isothermal solutions are superposed to satisfy the traction-free conditions on the crack surfaces. The dislocation solutions are applied to calculate the thermal stress intensity factors, which are validated by the exact solutions. The effects of the uniform heat flux, the ply angle and the crack length are investigated.     

JKR adhesion in cylindrical contacts

Planar JKR adhesive solutions use the half-plane assumption and do not permit calculation of indenter approach or visualization of adhesive force–displacement curves unless the contact is periodic. By considering a conforming cylindrical contact and using an arc crack analogy, we obtain closed-form indenter approach and load–contact size relations for a planar adhesive problem. The contact pressure distribution is also obtained in closed-form. The solutions reduce to known cases in both the adhesion-free and small-contact solution (Barquins, 1988) limits. The cylindrical system shows two distinct regimes of adhesive behavior; in particular, contact sizes exceeding the critical (maximum) size seen in adhesionless contacts are possible. The effects of contact confinement on adhesive behavior are investigated. Some special cases are considered, including contact with an initial neat-fit and the detachment of a rubbery cylinder from a rigid cradle. A comparison of the cylindrical solution with the half-plane adhesive solution is carried out, and it indicates that the latter typically underestimates the adherence force. The cylindrical adhesive system is novel in that it possesses stable contact states that may not be attained even on applying an infinite load in the absence of adhesion.     

Periodic heat conduction with relaxation time in cylindrical geometry

Steady-periodic heat conduction with relaxation time in an infinitely long hollow cylinder is considered. Four boundary value problems, with boundary conditions of the first and of the second kind, are solved analytically. The solution for a solid cylinder with a sinusoidally varying surface temperature is obtained as a special case of a solution found for the hollow cylinder. The effects of the relaxation time on the steady-periodic temperature field are analysed, in details, for a solid cylinder with a sinusoidally varying surface temperature and for a hollow cylinder with a sinusoidally varying heat flux at the inner surface and with a constant temperature at the outer surface. The results show that thermal resonances may occur and suggest that accurate measurements of the relaxation time could be obtained by means of experiments on steady-periodic heat conduction in cylindrical geometry. Received on 15 April 1997