Interaction of harmonic axisymmetric waves with a thin circular absolutely rigid separated inclusion

We study stress concentration near a circular rigid inclusion in an unbounded elastic body (matrix). In the matrix, there are wave motions symmetric with respect to the axis passing through the inclusion center and perpendicular to the inclusion. It is assumed that one of the inclusion sides is completely fixed to the matrix, while the other side is separated and the conditions of smooth contact are realized on that side. The solution method is based on the fact that the displacements caused by waves reflected from the inclusion are represented as a discontinuous solution of the Lamé equations. This permits reducing the original problem to a system of singular integral equations for functions related to the stress and displacement jumps on the inclusion. Its solution is constructed approximately by the collocation method with the use of special quadrature formulas for singular integrals. The approximate solution thus obtained permits numerically studying the stress state in the matrix near the inclusion. Technological defects or constructive elements in the form of thin rigid inclusions contained in machine parts and engineering structure members are stress concentration sources, which may result in structural failure. It is shown that the largest stress concentration is observed near separated inclusions. Static problems for elastic bodies with such inclusions have been studied rather comprehensively [1, 2]. The stress concentration near separated inclusions under dynamic actions on the bodies has been significantly less studied even in the case of harmonic vibrations. The results of these studies can be found in [3, 4], where bodies with a thin separated inclusion were considered, and in [5], where the problem about torsional vibrations of a body with a thin circular separated inclusion was studied. The aim of the present paper is to study stress concentration near such an inclusion in the case of interaction with harmonic waves under axial symmetry conditions.     

Interaction of plane elastic nonstationary waves with an elastic inclusion under complete adhesion

We solve the problem on the interaction of plane elastic nonstationary waves with a thin elastic strip-shaped inclusion. The inclusion is contained in an unbounded body (matrix) which in under conditions of plane strain. It is assumed that the condition of perfect adhesion between the inclusion and the matrix is satisfied. Because of the small thickness of the inclusion we assume that the bending and shear displacements at any inclusion point coincide with the displacements of the corresponding points of its midplane. The displacements on the midplane itself are found from the corresponding equations of the theory of plates. The statement of the boundary conditions for these equations takes into account the forces and moments acting on the inclusion edges from the matrix. The solution method is based on representing the displacements in the space of Laplace transforms as a discontinuous solution of the Lame’ equations for the plane strain with subsequent determining the transforms of the unknown jumps from integral equations. The passage to the original functions is performed numerically by methods based on replacement of the Mellin integral by the Fourier series. As a result, we obtain approximate formulas for calculating the stress intensity factors for the inclusion. These formulas are used to study the time dependence of the stress intensity factors and the influence of the inclusion rigidity on their values. We also study the possibility of treating inclusions of high rigidity as absolutely rigid inclusions.     

Interaction between a compliant disk-shaped inclusion and a crack upon incidence of an elastic wave

The propagation of harmonic elastic wave in an infinite three-dimensional matrix containing an interacting low-rigidity disk-shaped inclusion and a crack. The problem is reduced to a system of boundary integral equations for functions that characterize jumps of displacements on the inclusion and crack. The unknown functions are determined by numerical solution of the system of boundary integral equations. For the symmetric problem, graphs are given of the dynamic stress intensity factors in the vicinity of the circular inclusion and the crack on the wavenumber for different distances between them and different compliance parameters of the inclusion.     

含任意椭圆核各向异性板杂交应力有限元

基于含椭圆核有限大各向异性板弹性问题的复变函数级数解,应用杂交变分原理建立了一种与常规有限元相协调的含任意椭圆核各向异性板杂交应力有限元.单元内的应力场和位移场采用满足平衡方程、几何方程与物理方程的复变函数级数解,假设的复变函数级数解精确满足椭圆核边界处的位移协调条件和应力连续条件,单元外边界上的位移场按常规有限元位移场假设,单元内椭圆核的长轴可以与材料主轴不重合.单元刚度矩阵采用Gauss积分求得,并给出了建立刚度矩阵的主要公式和推倒过程.数值计算结果表明该单元具有计算精度高、计算工作量小等优点.     

Stress concentration near a thin elastic inclusion under interaction with harmonic waves in the case of smooth contact

We solve the problem on the interaction of plane harmonic waves with a thin elastic plate-shaped inclusion. The ambient medium is assumed to be in plane strain. The smooth contact conditions are satisfied on both sides of the inclusion. The bending displacements of the inclusion are determined from the corresponding differential equation. In the statement of boundary conditions for this equation, one should take into account the transverse forces and bending moments applied to the lateral edges of the inclusion, while the boundary conditions are posed on the midplane of the inclusion. Using the discontinuous solution method, we reduce the problem to a system of two singular integral equations, which are solved numerically by the mechanical quadrature method. We obtain approximate formulas for the stress intensity coefficients near the ends of the inclusion and for the transverse forces and moments applied to the inclusion.     

Interaction of plane nonstationary waves with a thin elastic inclusion under smooth contact conditions

We solve the problem of determining the stress state near a thin elastic inclusion in the form of a strip of finite width in an unbounded elastic body (matrix) with plane nonstationary waves propagating through it and with the forces exerted by the ambient medium taken into account. We assume that the matrix is in the plane strain state, and the smooth contact conditions are realized on both sides of the inclusion. The method for solving this problem consists in using the integral Laplace transform with respect to time and in representing the stress and displacement images in terms of the discontinuous solution of Lamé equations in the case of plane strain. As a result, the initial problem is reduced to a system of singular integral equations for the transforms of the unknown stress and displacement jumps. To invert the Laplace transform, we use a numerical method based on replacing the Mellin integral by the Fourier series. As a result, we obtain approximate formulas for calculating the stress intensity factors (SIF) for the inclusion, which are used to study the SIF time-dependence and its influence on the values of the inclusion rigidity. We also studied the possibility of considering the inclusions of higher rigidity as absolutely rigid inclusions.     

Interaction between rigid-disc inclusion and penny-shaped crack under elastic time-harmonic wave incidence

Influence of a rigid-disc massive inclusion on a neighboring penny-shaped crack induced by the time-harmonic wave propagation in an infinite elastic matrix is investigated by the numerical solution of associated 3D elastodynamic problem. No restrictions on the mutual orientation of interacting objects and direction of wave incidence are assumed. The inclusion is perfectly bonded with a matrix and supposes the translations and rotations, the crack faces are load-free. Frequency-domain problem is reduced to a system of boundary integral equations (BIEs) relative to the interfacial stress jumps (ISJs) on the inclusion and the crack opening displacements (CODs). The subtraction technique in conjunction with mapping technique, under taking into account the structure of solution at the fronts of inclusion and crack, is applied for regularization of BIEs obtained. A discrete analogue of equations is constructed by using the collocation scheme. Numerical calculations are carried out for the grazing incidence of a plane P-wave on the crack, where the interacting inclusion is coplanar and perpendicular to the crack, and has the same radius. The shielding and amplification effects of inclusion are assessed by the analysis of mode-I stress intensity factor (SIF) in the crack vicinity depending on the wave number, incident wave direction, position of the crack front point, inclusion mass, crack-inclusion orientation and distance.     

Interaction of general plane P wave and cylindrical inclusion partially debonded from its viscoelastic matrix

The interaction of a general plane P wave and an elastic cylindrical inclusion of infinite length partially debonded from its surrounding viscoelastic matrix of infinite extension is investigated. The debonded region is modeled as an arc-shaped interface crack between inclusion and matrix with non-contacting faces. With wave functions expansion and singular integral equation technique, the interaction problem is reduced to a set of simultaneous singular integral equations of crack dislocation density function. By analysis of the fundamental solution of the singular integral equation, it is found that dynamic stress field at the crack tip is oscillatory singular, which is related to the frequency of incident wave. The singular integral equations are solved numerically, and the crack open displacement and dynamic stress intensity factor are evaluated for various incident angles and frequencies. The project supported by the National Natural Science Foundation of China (19872002) and Climbing Foundation of Northern Jiaotong University     

Edge crack opening in an elastic plane with a wedge-like cut in contact with a punch

The equilibrium of an elastic plane with a wedge-like cut and an internal or edge crack on the symmetry axis was studied in [1] in the case of punch indentation in the lateral faces of the cut at a distance from the cut tip. In [1], the systemof singular integral equations of the problemwas solved numerically by the mechanical quadrature method. In this paper, the generalized Wiener-Hopf method [2] is used to obtain the analytic solution of a similar problem in the case of an edge crack under punch pressure on parts of the cut lateral faces adjacent to the cut tip. Some special cases of this problem were considered earlier without a crack [3, 4] or a punch [5, 6].     

Flexural vibrations of a thin circular elastic inclusion in an unbounded body under the action of a plane harmonic wave

The interaction of a plane harmonic longitudinal wave with a thin circular elastic inclusion is considered. The wave front is assumed to be parallel to the inclusion plane. Since the inclusion is thin, the matrix-inclusion interface conditions (perfect bonding) are formulated on the mid-plane of the inclusion. The bending displacements of the inclusion are determined from the bending equation for a thin plate. The problem is solved using discontinuous Lamé solutions for harmonic vibrations. Therefore, the problem can be reduced to the Fredholm equation of the second kind for a function associated with the discontinuity of normal stresses on the inclusion. The equation obtained is solved by the method of mechanical quadratures using Gaussian quadrature formulas. Approximate formulas for the stress intensity factors are derived. Results from a numerical analysis of the dependence of the SIFs on the dimensionless wave number and the stiffness of the inclusion are presented __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 5, pp. 16–21, May 2008.     

Interaction of plane harmonic waves with a thin elastic inclusion of zero flexural rigidity

We solve the problem on the interaction of plane elastic harmonic waves with a thin elastic strip-shaped inclusion. The inclusion is contained in an unbounded body (matrix) that is under plane strain conditions. The normal forces applied by the medium to the inclusion side edges are taken into account. Because of the small thickness of the inclusion, we assume that its flexural rigidity is zero and that the shear displacements at any of its points coincide with the displacements of the corresponding points of its midplane. The displacements on the midplane itself can be found from the corresponding equation of the theory of plates. The solution method consists in representing the displacements as discontinuous solutions of the Lamé equations and then determining the unknown jump from a singular integral equation. This equation is solved numerically by the collocation method, and formulas for the approximate calculation of the stress intensity factors near the inclusion ends are obtained.     

Influence of sliding interfaces on the response of a layered viscoelastic medium under a moving load

This article presents a method to compute the response of a viscoelastic layered half-space to a moving load when interlayer slip is considered. The Navier equations of equilibrium are solved for each layer in the frequency domain. The solution in the spatial coordinate system is subsequently obtained by means of Fast Fourier Transform and quadrature rules applied to integrable singularities. Following the global solution technique, the developed method compiles all the interface and the boundary conditions within a global matrix and it solves a unique linear system per couple of wave numbers. This method proves to be effective and is validated in an elastic case by comparison with the ALIZE-LCPC software that implements the Burmister axisymmetric solution. The influence of the interface sliding condition on the response of a layered viscoelastic medium is studied through an application to pavement structures. In this application, the effect of the load speed on vertical and horizontal profiles of the longitudinal strain and the normal stress is analyzed. It is shown, inter alia, that the maximum extension in the medium is not systematically observed at the location of an interface and that, as expected, low speeds and interlayer slip are more damaging to the structure when either a strain or a stress criterion is considered.     

Interaction between a circular inclusion and a symmetrically branched crack

The interaction problem between a circular inclusion and a symmetrically branched crack embedded in an infinite elastic medium is solved. The branched crack is modeled as three straight cracks which intersect at a common point and each crack is treated as a continuous contribution of edge dislocations. Green's functions are used to reduce the problem into a system of singular equations consisting of the distributions of Burger's dislocation vectors as unknown functions through the superposition technique. The resulting integral equations are solved numerically by the method of Gauss-Chebychev quadrature. The proposed integral equation approach is first verified for two limiting cases against the literature. More effort is paid on the effect of inclusion on both the Mode I and Mode lI stress intensity factors at the branch tips. The effect of inclusion on the branching path is also investigated.     

A dynamic photoelastic study of stress-wave propagation through an inclusion

A photoelastic study of the elastodynamic-stress fields around a circular, elastic inclusion (Solithane 113) embedded in an elastic plate (Hysol 4485) is presented. The edge of the plate was loaded by an explosive charge, which produced a plane, compressional stress wave of triangular shape. Isochromatic-fringe patterns were obtained, which give the maximum shear stresses, both inside the inclusion and in the surrounding medium. The principal stresses on the axis of symmetry were determined through the use of the oblique-incidence method. It was found that small tensile stresses are generated at the interface on the shadow side of the inclusion. The focusing effect inside the inclusion predicted by ray theory was not observed. Finally, the shape of the wavefront as the wave passes the inclusion was determined.     

Interaction of plane elastic harmonic waves with a perfectly bonded elastic inclusion

The interaction of plane harmonic waves with a thin elastic inclusion in the form of a strip in an infinite body (matrix) under plane strain conditions is studied. It is assumed that the bending and shear displacements of the inclusion coincide with the displacements of its midplane. The displacements in the midplane are found from the theory of plates. The priblem-solving method represents the displacements as discontinuous solutions of the Lamé equations and finds the unknown discontinuities solving singular integral equations by the numerical collocation method. Approximate formulas for the stress intensity factors at the ends of the inclusion are derived     

Dynamic stress intensity factors of two 3D rectangular cracks in a transversely isotropic elastic material under a time-harmonic elastic P-wave

The dynamic stress intensity factors (DSIFs) of two 3D rectangular cracks in a transversely isotropic elastic material under an incident harmonic stress wave are investigated by generalized Almansi’s theorem and the Schmidt method in the present paper. Using 2D Fourier transform and defining the jumps of displacement components across the crack surface as the unknown functions, three pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement components across the crack surfaces are expanded in a series of Jacobi polynomials. Numerical examples are provided to show the effects of the geometric shape of the rectangular crack, the characteristics of the harmonic wave and the distance between two rectangular cracks on the DSIFs of the transversely isotropic elastic material.     

Anisotropic bodies with inhomogeneities: The case of a set of cracks

We use the Betti theorem to obtain the integral equations of the dynamic theory of elasticity for a multilayer convex body with an arbitrary elastic anisotropy of layers containing plane infinitely thin cracks. The systems of integral equations relate the displacement jumps to the stresses on the crack lips and are stated numerically in terms of Fourier transforms. For the case of plane-parallel layers with a set of plane cracks on the interfaces between the layers, we propose a simple numerical-analytic method for constructing the Fourier symbol, i.e., the matrix of the kernel of the system of integral equations. The method is stable for an arbitrary combination of continuous and discontinuous conditions on the layer boundaries. Numerical examples are given for a packet of four heterogeneous anisotropic layers.     

Analysis of shock wave reflection from fixed and moving boundaries using a stabilized particle method

In the present paper, the efficiency of an enhanced formulation of the stabilized corrective smoothed particle method （CSPM） for simulation of shock wave propagation and reflection from fixed and moving solid boundaries in compressible fluids is investigated. The Lagrangian nature and its accuracy for imposing the boundary conditions are the two main reasons for adoption of CSPM. The governing equations are further modified for imposition of moving solid boundary conditions. In addition to the traditional artificial viscosity, which can remove numerically induced abnormal jumps in the field values, a velocity field smoothing technique is introduced as an efficient method for stabilizing the solution. The method has been implemented for one- and two-dimensional shock wave propagation and reflection from fixed and moving boundaries and the results have been compared with other available solutions. The method has also been adopted for simulation of shock wave propagation and reflection from infinite and finite solid boundaries.     

Jeffery solution for an elastic disk containing a sliding eccentric circular inclusion assembled by interference fit

An analytic solution is presented for stresses induced in an elastic and isotropic disk by an eccentric press-fitted circular inclusion. The disk is also subject to uniform normal stress applied at its outer border. The inclusion is assumed to be of the same material as the annular disk and both elements are in a plane stress or plane strain state. A frictionless contact condition is assumed between the two members. The solution is obtained by using the general expression for a biharmonic stress function in bipolar coordinates. The results show that the maximum of the von Mises effective stress due to the inclusion interference occurs in the ligament for large eccentricity, but it deviates from the symmetry axis for small eccentricity. Moreover, along the border of the circular inclusion the hoop stress locally coincides with the contact pressure, in agreement with a similar classical result valid for a half plane.     

On the load transfer from elastic inclusions to an infinite elastic orthotropic plane with cuts

The present paper deals with the problem of load transfer from elastic inclusions to an infinite elastic orthotropic plane with cuts located on one of the principal orthotropy directions. The constitutive system of equations of this problem is derived under the assumption that the inclusions are in a uniaxial stress state. The obtained system consists of a singular integro-differential equation and a singular integral equation for the jumps of the tangential stresses acting on the inclusion shores and for the derivative of the the cut opening function. The behavior of solutions of the system of constitutive equations at the endpoints of the inclusions and cuts is studied, and the solution of this system is constructed by the numerical-analytic discrete singularity method.