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 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.     

Debonded arc-shaped interface conducting rigid line inclusions in piezoelectric composites

We consider an arc-shaped conducting rigid line inclusion located at the interface between a circular piezoelectric inhomogeneity and an unbounded piezoelectric matrix subjected to remote uniform anti-plane shear stresses and in-plane electric fields. Moreover, one side of the rigid line inclusion has become fully debonded from the matrix or the inhomogeneity leading to the formation of an insulating crack. After the introduction of two sectionally holomorphic vector functions, the problem is reduced to a vector Riemann–Hilbert problem, which can be decoupled sequentially by repeated application of the orthogonality relations between the eigenvectors for two corresponding generalized eigenvalue problems.     

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.     

Contact stress analysis of an elastic half-plane containing multiple inclusions

This paper presents a two-dimensional contact stress analysis to investigate the effects of multiple inclusions on the contact pressure and subsurface stresses in an elastic half-plane. The boundary element method is used to analyze the contact problem where a set of integral equations is derived on the contact region and the matrix–inclusion interfaces. As the contact region is unknown a priori, an iterative procedure is implemented to determine the actual contact region and the contact pressure, and the tractions and displacements on the matrix–inclusion interfaces are obtained by solving the integral equations numerically. Numerical results show that the inclusions near contact surface could cause significant alterations in the contact pressure distribution. The stiff inclusions could toughen the surrounding material and reduce the internal stresses while the soft inclusions could increase the subsurface stresses.     

PLANE ELASTIC PROBLEM ON RIGID LINES ALONG CIRCULAR INCLUSION

The plane elastic problem of circular-arc rigid line inclusions is considered. The model is subjected to remote general loads and concentrated force which is applied at an arbitrary point inside either the matrix or the circular inclusion. Based on complex variable method, the general solutions of the problem were derived. The closed form expressions of the sectionally holomorphic complex potentials and the stress fields were derived for the case of the interface with a single rigid line. The exact expressions of the singular stress fields at the rigid line tips were calculated which show that they possess a pronounced oscillatory character similar to that for the corresponding crack problem under plane loads. The influence of the rigid line geometry, loading conditions and material mismatch on the stress singularity coefficients is evaluated and discussed for the case of remote uniform load.     

Interacting elliptic inclusions by the method of complex potentials

An accurate series solution has been obtained for a piece-homogeneous elastic plane containing a finite array of non-overlapping elliptic inclusions of arbitrary size, aspect ratio, location and elastic properties. The method combines standard Muskhelishvili’s representation of general solution in terms of complex potentials with the superposition principle and newly derived re-expansion formulae to obtain a complete solution of the many-inclusion problem. By exact satisfaction of all the interface conditions, a primary boundary-value problem stated on a complicated heterogeneous domain has been reduced to an ordinary well-posed set of linear algebraic equations. A properly chosen form of potentials provides a remarkably simple form of solution and thus an efficient computational algorithm. The theory developed is rather general and can be applied to solve a variety of composite mechanics problems. The advanced models of composite involving up to several hundred inclusions and providing an accurate account for the microstructure statistics and fiber–fiber interactions can be considered in this way. The numerical examples are given showing high accuracy and numerical efficiency of the method developed and disclosing the way and extent to which the selected structural parameters influence the stress concentration at the matrix–inclusion interface.     

Stress-strain state of an anisotropic plate with curved cracks and thin rigid inclusions

We solve the bending problem for an anisotropic plate with flaws like smooth curved nonoverlapping through cracks and rigid inclusions. The problem is solved by the method of Lekhnitskii complex potentials specified as Cauchy type integrals over the flaw contours with an unknown integrand density function. We use the Sokhotskii—Plemelj formulas to reduce the boundary-value problem to a system of singular integral equations with the additional conditions that the displacements in the plate are single-valued when going around the cut contours and the equilibrium conditions for stress-free rigid inclusions. After the singular integrals are approximated by the Gauss-Chebyshev quadrature formulas, the problem is reduced to solving a system of linear algebraic equations. We study the local stress distribution near flaw tips. We analyze the mutual influence of flaws on the stress distribution character near their vertices and compare the well-known solutions for isotropic plates with the solutions obtained by passing to the limit in the anisotropy parameters (“weakly anisotropic material”) and by using the method proposed here.     

Thermal stresses in an elastic space with a perfectly rigid flat inclusion under perpendicular heat flow

This paper examines the three-dimensional problem of finding thermal stresses due to an insulated rigid sheet-like inclusion (anticrack) in an elastic space under a uniform perpendicular heat flow. By using appropriate harmonic potentials, a general method of solving this problem is presented. The resulting boundary-value problems are reduced to classical mixed problems of potential theory. For the purpose of illustration, a complete solution in terms of elementary functions for a rigid circularly shaped inclusion is given and discussed from the point of view of material failure.     

Interfacial Energies for Incoherent Inclusions

We study a variational problem describing an incoherent interface between a rigid inclusion and a linearly elastic matrix. The elastic material is allowed to slip relative to the inclusion along the interface, and the resulting mismatch is penalized by an interfacial energy term that depends on the surface gradient of the relative displacement. The competition between the elastic and interfacial energies induces a threshold effect when the interfacial energy density is non-smooth: small inclusions are coherent (no mismatch); sufficiently large inclusions are incoherent. We also show that the relaxation of the energy functional can be written as the sum of the bulk elastic energy functional and the tangential quasiconvex envelope of the interfacial energy functional.     

Bond stresses for a composite material reinforced with various arrays of fibers

The problem treated here is that of an isotropic body having a doubly periodic rectangular or triangular array of perfectly bonded circular elastic inclusions. The body is in tension or compression. This simulates a composite material wherein a relatively weak matrix is reinforced by stronger (and more rigid) fibers. Bond stresses for both rectangular and triangular arrays have been calculated using either boundary point matching or boundary point least squares techniques. Numerical results based on a plane strain analysis are given in graphical form.     

Effective thermoelastic properties of graded doublyperiodic particulate matrix composites in varying externalstress fields

We consider a linear elastic composite medium, which consists of a homogeneousmatrix containing aligned ellipsoidal uncoated or coated inclusions arranged in a doubly periodicarray and subjected to inhomogeneous boundary conditions. The hypothesis of effective fieldhomogeneity near the inclusions is used. The general integral equation obtained reduces theanalysis of infinite number of inclusion problems to the analysis of a finite number of inclusions insome representative volume element (RVE) . The integral equation is solved by a modifiedversion of the Neumann series; the fast convergence of this method is demonstrated for concreteexamples. The nonlocal macroscopic constitutive equation relating the cell averages of stress andstrain is derived in explicit iterative form of an integral equation. A doubly periodic inclusion fieldin a finite ply subjected to a stress gradient along the functionally graded direction is considered.The stresses averaged over the cell are explicitly represented as functions of the boundaryconditions. Finally, the employed of proposed explicit relations for numerical simulations oftensors describing the local and nonlocal effective elastic properties of finite inclusion pliescontaining a simple cubic lattice of rigid inclusions and voids are considered. The local andnonlocal parts of average strains are estimated for inclusion plies of different thickness. Theboundary layers and scale effects for effective local and nonlocal effective properties as well as foraverage stresses will be revealed.     

Prediction of the effective elastic properties of spheroplastics by the generalized self-consistent method

The problem of predicting the effective elastic properties of composites with prescribed random location and radius variation in spherical inclusions is solved using the generalized self-consistent method. The problem is reduced to the solution of the averaged boundary-value problem of the theory of elasticity for a single inclusion with an inhomogeneous transition layer in a medium with desired effective elastic properties. A numerical analysis of the effective properties of a composite with rigid spherical inclusions and a composite with spherical pores is carried out. The results are compared with the known solution for the periodic structure and with the solutions obtained by the standard self-consistent methods. Perm’ State Technical University, Perm’ 614600. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 3, pp. 186–190, May–June, 1999.     

A coated rigid elliptical inclusion loaded by a couple in the presence of uniform interfacial and hoop stresses

We consider a confocally coated rigid elliptical inclusion, loaded by a couple and introduced into a remote uniform stress field. We show that uniform interfacial and hoop stresses along the inclusion–coating interface can be achieved when the two remote normal stresses and the remote shear stress each satisfy certain conditions. Our analysis indicates that: (i) the uniform interfacial tangential stress depends only on the area of the inclusion and the moment of the couple; (ii) the rigid-body rotation of the rigid inclusion depends only on the area of the inclusion, the coating thickness, the shear moduli of the composite and the moment of the couple; (iii) for given remote normal stresses and material parameters, the coating thickness and the aspect ratio of the inclusion are required to satisfy a particular relationship; (iv) for prescribed remote shear stress, moment and given material parameters, the coating thickness, the size and aspect ratio of the inclusion are also related. Finally, a harmonic rigid inclusion emerges as a special case if the coating and the matrix have identical elastic properties.     

Thermal stress analysis of a three-dimensional anticrack in a transversely isotropic solid

A three-dimensional analysis is performed for an infinite transversely isotropic elastic body containing an insulated rigid sheet-like inclusion (an anticrack) in the isotropy plane under a remote perpendicularly uniform heat flow. A general solution scheme is presented for the resulting boundary-value problems. Accurate results are obtained by constructing suitable potential solutions and reducing the thermal problem to a mechanical analog for the corresponding isotropic problem. The governing boundary integral equation for a planar anticrack of arbitrary shape is obtained in terms of a normal stress discontinuity. As an illustration, a complete solution for a rigid circular inclusion is obtained in terms of elementary functions and analyzed. This solution is compared with that corresponding to a penny-shaped crack problem.     

On the singularity induced by certain mixed boundary conditions in linearized and nonlinear elastostatics

This study concerns the local character of the elastostatic field in plane strain near a point that separates a free from an adjoining fixed segment of a rectilinear boundary-component. The well-known singular field behavior predicted by the linear theory, as such a point is approached, exhibits oscillatory deformations and stresses. It is shown here by means of an asymptotic analysis that the foregoing anomalous behavior does not occur within the nonlinear theory of harmonic elastic materials. In preparation for this task certain general aspects of the latter theory are reviewed. The results obtained in the nonlinear asymptotic treatment of the class of mixed boundary-value problems considered are discussed in detail with particular attention to the problem of a bonded flat-ended rigid punch.     

Interaction of a plane harmonic wave with a thin rigid inclusion of the shape of a cylindrical shell

The paper presents the solution of the problem of determining the stress state in an elastic matrix containing a rigid inclusion of the shape of a thin cylindrical shell. It is assumed that harmonic vibrations occur in the matrix under the conditions of axial symmetry (the symmetry axis is the inclusion axis) and the conditions of full adhesion between the inclusion and the matrix are satisfied. The vibrations are caused by the propagation of a plane wave whose front is perpendicular to the inclusion axis. The solution method is based on representing the displacements in the matrix as discontinuous solutions of the equations of axisymmetric oscillations of an elastic medium with unknown stress jumps on the inclusion surface. The realization of the boundary conditions for these jumps leads to a system of integral equations. Its solution is constructed numerically by the mechanical quadrature method with the use of special quadrature formulas for specific integrals. It is numerically investigated how the ratio of the inclusion geometric dimensions and the propagating wave frequency affect the stress concentration near the inclusion.     

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.     

Antiplane problem of circular ARC interfacial rigid line inclusions

The antiplane problem of circular arc rigid line inclusions under antiplane concentrated force and longitudinal shear loading was dealt with. By using RiemannSchwaz‘s symmetry principle integrated with the singularity analysis of complex functions, the general solution of the problem and the closed form solutions for some important practical problems were presented. The stress distribution in the immediate vicinity of circular arc rigid line end was examined in detail. The results show that the singular stress fields near the rigid inclusion tip possess a square-root singularity similar to that for the corresponding crack problem under antiplane shear loading, but no oscillatory character. Furthermore, the stresses are found to depend on geometrical dimension, loading conditions and materials parameters. Some practical results concluded are in agreement with the previous solutions.     

On multiple circular inclusions in plane magnetoelasticity

A general series solution to the magnetoelastic problem of interacting circular inclusions in plane magnetoelasticity is provided in this paper. By the use of complex variable theory and Laurent series expansion method, the general expression of the magnetic and the magnetoelastic complex potentials for the circular inclusion problem is derived. Expanding the definition of the Airy’s stress function of pure elastic field into the magnetoelastic field and applying the superposition method, the general expression then can be reduced to a set of linear algebraic equations and solved in a series form. An approximate closed form solution for the case of two arbitrarily located inclusions is also provided. For illustrating the effect of the pertinent parameters, the numerical results of the interfacial magnetoelastic stresses are displayed in graphic form.