Mathematical and Numerical Simulation of Equilibrium of an Elastic Body Reinforced by a Thin Elastic Inclusion

A boundary value problem describing the equilibrium of a two-dimensional linear elastic body with a thin rectilinear elastic inclusion and possible delamination is considered. The stress and strain state of the inclusion is described using the equations of the Euler–Bernoulli beam theory. Delamination means the existence of a crack between the inclusion and the elastic matrix. Nonlinear boundary conditions preventing crack face interpenetration are imposed on the crack faces. As a result, problem with an unknown contact domain is obtained. The problem is solved numerically by applying an iterative algorithm based on the domain decomposition method and an Uzawa-type algorithm for solving variational inequalities. Numerical results illustrating the efficiency of the proposed algorithm are presented.     

Dynamic contact between a spherical inclusion and a matrix upon incidence of an elastic wave

Using the boundary element method, we have studied the dynamic displacements and stresses in an infinite elastic matrix with a spherical elastic inclusion, caused by the propagation of an elastic wave. The original problem has been reduced to a system of boundary integral equations for the contact displacements and tractions on the interface between the inclusion and matrix. Based on the numerical solution of these equations, we have analyzed the influence of the direction of wave propagation and frequency on the important physical parameters, depending on the elastic characteristics of composite constituents.     

三角形弹性夹杂对裂纹的影响

研究了三角形弹性夹杂和裂纹之间的相互影响问题。应用Chau和Wang导出的面力边值问题的边界积分方程为基本方程,用夹杂和基体交界面上的面力和位移的连续性条件为补充方程,从而得到了一组能够解决夹杂和裂纹相互影响问题的方程,最后的方程组用一种新的边界单元法求解。计算了各种不同的夹杂和基体的材料常数以及夹杂和基体之间不同距离情况下裂纹尖端的应力强度因子。文中结果对研究新型复合材料有一定的应用价值。     

Averaging of a finely laminated elastic medium with roughness or adhesion on the contact surfaces of the layers

The governing relations of a laminated elastic medium with non-ideal contact conditions in the interlayer boundaries are obtained by an asymptotic averaging method. The interaction of rough surfaces is described by a non-linear contact condition which simulates the local deformation of the microroughnesses using a certain penetration of the nominal surfaces of the elastic layers. The cohesive forces, caused by the thin adhesive layer, are described within the limits of the Frémond model which includes a differential equation characterizing the change in the cohesion function. A piecewise-linear approximation of the initial positive segment of the Lennard–Jones potential curve is proposed to describe of the adhesive forces between smooth dry surfaces. A comparison is made with the solution obtained within the limits of the Maugis–Dugdale model based on a piecewise-constant approximation. Solutions of the above problems are constructed taking account of the possible opening of interlayer boundaries.     

裂纹与线夹杂垂直接触问题的奇性分析

以短纤维复合等作为工程背景,采用现有文献中单裂纹和单夹杂的基本解,对于限平项（基体）上,裂纹和线夹杂的垂直接触问题从断理解力学的角度作了研究,得到了问题的积分方程,推出了接触点的性线指数,奇性应力及以此表示的接触点附近三个区域内的应力强度因子表达式,并给出一些数值结果,可供工程实际参考。     

Dynamic proper colorings of a graph

We solve an axisymmetric problem of the interaction of harmonic waves with a thin elastic circular inclusion located in an elastic isotropic body (matrix). On both sides of the inclusion, between it and the body (matrix), conditions of smooth contact are realized. The method of solution is based on the representation of displacements in the matrix in terms of discontinuous solutions of Lamé equations for harmonic vibrations. This enables us to reduce the problem to Fredholm integral equations of the second kind for functions related to jumps of normal stress and radial displacement on the inclusion.     

Asymptotics of anisotropic weakly curved inclusions in an elastic body

Under study are the boundary value problems describing the equilibrium of twodimensional elastic bodies with thin anisotropic weakly curved inclusions in presence of separations. The latter implies the existence of a crack between the inclusion and the matrix. Nonlinear boundary conditions in the form of inequalities are imposed on the crack faces that exclude mutual penetration of the crack faces. This leads to the formulation of the problems with unknown contact area. The passage to limits with respect to the rigidity parameters of the thin inclusions is inspected. In particular, we construct the models as the rigidity parameters go to infinity and analyze their properties.     

Thermally Stressed State of Bodies with Two Ellipsoidal Inclusions

Solutions are obtained for the interaction of two ellipsoidal inclusions in an elastic isotropic matrix with polynomial external athermal and temperature fields. Perfect mechanical and temperature contact is assumed at the phase interface. A solution to the problem is constructed. When the perturbations in the temperature field and stresses in the matrix owing to one inclusion are re-expanded in a Taylor series about the center of the second inclusion, and vice versa, and a finite number of expansion terms is retained, one obtains a finite system of linear algebraic equations in the unknown constants. The effect of a force free boundary of the half space on the stressed state of a material with a triaxial ellipsoidal inhomogeneity (inclusion) is investigated for uniform heating. Here it was assumed that the elastic properties of the inclusions and matrix are the same, but the coefficients of thermal expansion of the phases differ. Studies are made of the way the stress perturbations in the matrix increase and the of the deviation from a uniform stressed state inside an inclusion as it approaches the force free boundary.     

A Contact Problem of the Interaction of a Semi-Finite Inclusion with a Plate

A piece wise-homogeneous plane made up of twodifferent materials and reinforced by an elastic unclusion is considered on a semi-finite section where the different materials join. Vertical and horizontal forces are applied to the inclusion which haz a variable thichness and a variable elasticity modulus.Under certain conditions the problem is reduced to integrodifferential equations of third order. The solution is constructed effectively by applying the methods of theory of analytic functions to a boundary value problem of the Carleman type for a strip. Asymptotic estimates of normal contact stress are obtained.     

Dual boundary element method for problems of the theory of thin inclusions

We have proposed to apply the dual boundary element method in problems of the theory of thin inclusions. The contact conditions on the boundary of a thin inclusion are considered as jumps of displacements and stresses in the body on the median surface of this defect. Thus, the relations between the unknown discontinuities and average values of the displacements and stresses are a model of inclusions. For rectilinear boundary elements, we have constructed models of inclusions, taking into account the tension, shear, and bending of a thin inclusion. Examples for rectilinear and curved inclusions have been considered. Comparison of the results obtained by the proposed technique with data based on the direct approach shows the efficiency of the proposed method.     

Seepage refraction in a semicircular lens located at the boundary of two porous massifs

The problem of calculating the two-dimensional seepage field in a structurally inhomogeneous three-component medium in the form of two infinitely porous massifs with a semicircular inclusion in their plane boundary of contact is considered. The distribution of the seepage rate, when two matching conditions along the lines of contact of unlike zones are strictly satisfied, is obtained in closed analytic form by methods of complex analysis. Limiting cases of the conduction of the components of the medium and cases of the degeneration of a three-component medium into a two-component medium are considered.     

Dynamic interaction of an elastic medium with a thin-walled curvilinear piezoelectric inclusion under longitudinal vibrations of a composite

Using the method of matched asymptotic expansions, we obtain models of dynamic interaction of a thin-walled curvilinear piezoelectric inclusion of variable thickness with an elastic isotropic matrix under stationary vibrations of the composite. The elastic system is under conditions of longitudinal shear. Different cases of electric boundary conditions on the surface of the heterogeneity are considered. We propose an algorithm for the construction of boundary layer corrections for refining the behavior of displacements and stresses in the vicinity of the edge of the inclusion for its different shapes.     

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.     

The conjunction problem for thin elastic and rigid inclusions in an elastic body

Under consideration is the conjunction problem for a thin elastic and a thin rigid inclusions that are in contact at one point and placed in an elastic body. Depending on what kind of conjunction conditions are set at the contact point of inclusions, we consider the two cases: the case of no fracture, where, as the conjunction conditions, we take the matching of displacements at the contact point and preservation of the angle between the inclusions, and the case with a fracture in which only the matching of displacements is assumed. At the point of conjunction, we obtain the boundary conditions for the differential formulation of the problem. On the positive face of the rigid inclusion, there is delamination. On the crack faces, some nonlinear boundary conditions of the type of inequalities are set, that prevent mutual penetration of the faces. The existence and uniqueness theorems for the solution of the equilibrium problem are proved in both cases.     

Numerical solution of an equilibrium problem for an elastic body with a thin delaminated rigid inclusion

Under consideration is a 2D-problem of elasticity theory for a body with a thin rigid inclusion. It is assumed that there is a delamination crack between the rigid inclusion and the elastic matrix. At the crack faces, the boundary conditions are set in the form of inequalities providing mutual nonpenetration of the crack faces. Some numerical method is proposed for solving the problem, based on domain decomposition and the Uzawa algorithm for solving variational inequalities.We give an example of numerical calculation by the finite element method.     

含椭圆形夹杂的压电体平面应变问题

利用复变函数理论,对在无限远处均匀应力和电位移载荷作用下的含有椭圆形弹性夹杂的横观各向同性压电材料,作了力电分析．在该文有限元结果和前人相关理论解的基础上,提了出一个可接受的认为弹性夹杂体内的应力场为常应力场的假设．在采用了不导通电边界条件之后,获得了以复势形式表示的压电基体的和弹性夹杂体内部的应力场解．     

Instability of thermoelastic interaction between a half-space and a rigid base through a thin liquid layer

We investigate the instability of thermoelastic interaction between elastic and rigid half-spaces through a liquid interlayer under the conditions of heat transfer across the interfaces. Due to the small thickness of the liquid layer, its influence on the temperature field is taken into account by the thermal resistance of the contact between the bodies, which depends on the normal displacement of the boundary of the elastic body. The pressure inside the liquid is equal to the external pressure applied to the bodies. We determined the critical value of the external heat flow for which the instability becomes possible in such a system and studied the dependence of this value on the parameters of the elastic half-space, the thickness of the liquid layer, and its thermal conduction. Institute of Mathematics, Ukrainian Academy of Sciences, Kiev. Translated from Matematychni Metody ta Fizyko-Mekhanichni Polya, Vol. 41, No. 2, pp. 76–82, April–June, 1998.     

The Inclusion Interval of Basic Coneigenvalues of a Matrix

In this paper,a compound-type inclusion interval of basic coneigenvalues of(com- plex)matrix is obtained.The corresponding boundary theorem and isolating theorem are given.     

In-plane loading of a cracked elastic solid by a disc inclusion with a Mindlin-type constraint

The paper examines the problem related to the axisymmetric interaction between an external circular crack and a centrally placed penny-shaped rigid inclusion located in the plane of the crack. The interface between the inclusion and the elastic medium exhibits a Mindlin-type imperfect bi-lateral contact. Analytical results presented in the paper illustrate the manner in which the lateral translational stiffness of the inclusion and the stress intensity factor at the boundary of the external circular crack are influenced by the inclusion/crack radii ratio.     

Analysis of a quasistatic contact problem for piezoelectric materials

We consider a mathematical model which describes the frictional contact between a piezoelectric body and an electrically conductive foundation. The process is quasistatic, the material behavior is modeled with an electro-viscoelastic constitutive law and the contact is described with subdifferential boundary conditions. We derive the variational formulation of the problem which is in the form of a system involving two history-dependent hemivariational inequalities in which the unknowns are the velocity and electric potential field. Then we prove the existence of a unique weak solution to the model. The proof is based on a recent result on history-dependent hemivariational inequalities obtained in Migórski et al. (submitted for publication) [16].