Application of complex path-independent integrals to the solution of the problem of a straight crack in a finite plane isotropic elastic medium

The classical problem of a straight crack in a finite, plane, isotropic, elastic medium of arbitrary shape is reconsidered by the well-known method of Muskhelishvili for such a crack (but in an infinite medium). Both the crack and the boundary of the medium are assumed loaded in an arbitrary way. It is shown that this problem can be completely solved if the numerical values of the first complex potential (z) of Muskhelishvili are known along a closed contour surrounding the crack, probably along the boundary of the medium. To this end, complex path-independent integrals associated with (z) and Chebyshev polynomials have been used. Numerical results for the stress intensity factors are displayed in an application. Generalizations of the method are also proposed and the second fundamental crack problem, the problem of a crack in an anisotropic medium and the problem of an interface crack between two isotropic media are considered in some detail.     

The use of complex valued functions for the solution of three-dimensional elasticity problems

For the treatment of plane elasticity problems the use of complex functions has turned out to be an elegant and effective method. The complex formulation of stresses and displacements resulted from the introduction of a real stress function which has to satisfy the 2-dimensional biharmonic equation. It can be expressed therefore with the aid of complex functions. In this paper the fundamental idea of characterizing the elasticity problem in the case of zero body forces by a biharmonic stress function represented by complex valued functions is extended to 3-dimensional problems. The complex formulas are derived in such a way that the Muskhelishvili formulation for plane strain is included as a special case. As in the plane case, arbitrary complex valued functions can be used to ensure the satisfaction of the governing equations. Within the solution of an analytical example some advantages of the presented method are illustrated.     

The degenerate scale problem for the Laplace equation and plane elasticity in a multiply connected region with an outer circular boundary

This paper investigates the degenerate scale problem for the Laplace equation and plane elasticity in a multiply connected region with an outer circular boundary. Inside the boundary, there are many voids with arbitrary configurations. The problem is analyzed with a relevant homogenous BIE (boundary integral equation). It is assumed that all the inner void boundary tractions are equal to zero, and tractions on the outer circular boundary are constant. Therefore, all the integrations in BIE are performed on the outer circular boundary only. By using the relation z * conjg(z) = a * a, or conjg(z) = a * a/z on the circular boundary with radius a, all integrals can be reduced to an integral for complex variable and they can be integrated in closed form. The degenerate scale a = 1 is found in the Laplace equation and in plane elasticity regardless of the void configuration.     

Complex variable formulation for non-slipping plane strain contact of two elastic solids in the presence of interface mismatch eigenstrain

In this paper, the problems of non-slipping contact, non-slipping adhesive contact, and non-slipping adhesive contact with a stretched substrate are sequentially studied under the plane strain theory. The main results are obtained as follows:(i) The explicit solutions for a kind of singular integrals frequently encountered in contact mechanics (and fracture mechanics) are derived, which enables a comprehensive analysis of non-slipping contacts. (ii) The non-slipping contact problems are formulated in terms of the Kolosov–Muskhelishvili complex potential formulae and their exact solutions are obtained in closed or explicit forms. The relative tangential displacement within a non-slipping contact is found in a compact form. (iii) The spatial derivative of this relative displacement will be referred to in this study as the interface mismatch eigenstrain. Taking into account the interface mismatch eigenstrain, a new non-slipping adhesive contact model is proposed and its solution is obtained. It is shown that the pull-off force and the half-width of the non-slipping adhesive contact are smaller than the corresponding solutions of the JKR model (Johnson et al., 1971). The maximum difference can reach 9% for pull-off force and 17% for pull-off width, respectively. In contrast, the new model may be more accurate in modeling the non-slipping adhesion. (iv) The non-slipping adhesions with a stretch strain (S-strain) imposed to one of contact counterparts are re-examined and the analytical solutions are obtained. The accurate analysis shows that under small values of the S-strain both the natural adhesive contact half-width and the pull-off force may be augmented, but for the larger S-strain values they are always reduced. It is also found that Dundurs’ parameter β may exert a considerable effect on the solution of the pull-off problem under the S-strain.These solutions may be used to study contacts at macro-, micro-, and nano-scales.     

Doubly-periodic array of cracks in an infinite isotropic medium

The problem of a doubly-periodle array of curvilinear cracks in an infinite isotropic medium under conditions of generalized plane stress or plane strain is reduced, by using the method of complex potentials of Muskhelishvili [1], to a complex Cauchy-type singular integral equation on one of the cracks, which can be further numerically solved by reduction to a system of two real Cauchy-type singular integral equations and application of the Gauss or Lobatto-Chebyshev method of numerical solution of such equations.

---

Zusammenfassung Das Problem einer zwei-periodischer Anordung gekrümmten Rissen in einem unendlichem isotropem Medium im allgemeinen ebener Spannung-bzw Dehnung-Zustand, unter Anwendung der Methode von komplexen Potentialen von Muskhelishvili [1], ist in einer Cauchyscher komplexer singulärer Integral-Gleichung auf einem der Rissen reduziert. Diese Gleichung kann weiter numerisch gelöst werden bei Transformation in einem System von zwei reelen Cauchyschen singulären Integral-Gleichungen und bei Anwendung der Methoden von Gauss oder Lobatto-Chebyshev für die numerische Resolution solcher Gleichungen.

     

Study of the plane problem for a physically nonlinear elastic solid by methods of the theory of functions of one complex Variable

Successful application of methods of complex analysis in linear elasticity problems, initiated by Kolosov, Muskhelishvili, Vekua, and their students, serves as a basis for similar studies in the field of analytical-numerical approximations to solutions of boundary value problems and various nonlinear equations of mathematical physics. In the present paper, we suggest a method for solving plane boundary-value problems for a special class of physically nonlinear elastic solids in the case of small strains. This general method, which can be used for a wide class of domains, is illustrated by the example of a square domain with boundary conditions given in stresses. These methods can also easily be used for boundary conditions of other types.     

三维复杂流动的间断有限元方法模拟

本文基于三维可压缩Euler方程,采用基于Runge-Kutta时间离散的间断有限元方法(RKDG方法),对三维前台阶、三维Riemann问题和球Riemann等问题进行了模拟。结果表明,本文的RKDG方法能够在很少的网格内清晰地捕捉到三维复杂流场中的激波和接触间断;同时,将球Riemann问题中z=0.4平面压强沿到对称轴距离的分布与文献中的近似精确解相比,吻合较好,这也验证了本文的RKDG方法不仅能够进行三维复杂流场的定性描述,也能够应用于三维复杂流场的定量计算。     

Half-planes without coupling under contact loading

Summary Nonrotating half-planes in contact under oblique loading are investigated in this paper. The solution is based on the influence integrals of the Flamant solution. The problem is determined by two integral equations for the normal and tangential stresses, which are uncoupled for special cases, as bodies of similar material in contact. In order to simplify the singular integrals, the method of superposition of flat punches is used. The result for the symmetric case is almost identical with the axisymmetric solution. For polynomial profiles of the form x ^s , the Muskhelishvili potentials can be written in terms of a complex hypergeometric function. The interior stress field is illustrated for an example. Accepted for publication 13 July 1996     

Elastic plates under bending solved by pseudocaustics

The method pf pseudocaustics was applied to the study of out-of-plane bending in elastic plates. It is shown that for bending problems where the loading mode is given, the method determines experimentally the complex potential function at selected points along the boundaries. A conformal mapping of the closed smooth curves of each boundary of the plate on to a unit circle allows the determination of the complex potential ϕ (ζ), expressed in the form of a Laurent series. This in turn yields the complete solution of the bent plate. In order to show the efficiency of the method it was applied to two typical examples of thin infinite plates in cylindrical bending, having either a circular central hole, or a square hole. The experimental results corroborate the theoretical results, thus proving that this combined theoretical and experimental method is suitable for solving elastic problems in applications with high accuracy, where other methods fail to yield satisfactory results.     

On the representation of three-dimensional elasticity solutions with the aid of complex valued functions

In this paper the representation of three-dimensional displacement fields in linear elasticity in terms of six complex valued functions is considered. The representation includes the complex Muskhelishvili formulation for plane strain as a special case. The completeness of the complex representation for regular solutions is shown and a relationship to the Neuber/Papkovich solutions is given.     

The comparison model method: A new arithmetic approach to the discrete inverse problem of groundwater hydrology

Let the steady-state pressure z(·) of a fluid in a one-dimensional domain be governed by the equation d _x (a d _x z) = f subject to Dirichlet boundary conditions. We consider the identification of the transmissivity a (·), given f(·), and measured pressure z(·) by the comparison model method, a direct method which has been known and applied for some time but lacked theoretical background. With reference to a domain in one spatial dimension, we examine both the infinite-(‘continuous’) and finite-(discrete) dimensional cases. In the former, the method is based on the solution p(·) of an auxiliary flow equation, where f(·) and the two-point boundary conditions are unchanged with respect to the original or z(·) equation, whereas a tentative constant value b is assigned to the auxiliary transmissivity. The ratio of the first derivatives of p(·) and z(·) multiplied by b yields a solution ã(·) to the inverse problem. We examine in detail the nonuniqueness of ã(·) as a function of b, some cases where existence implies uniqueness, the role of positivity constraints, and a special feature: self-identifiability. We then translate all available results into the discrete case, where the good unknowns for the inverse problem are the internode coefficients. Several algebraic and numerical examples are presented.     

On steady plane flows in classical elasticity

A new class of boundary-value problems in mathematical elasticity is proposed, wherein the medium flows steadily relative to a non-embedded surface over which tractions or velocities are prescribed. Such flows are seen in metal forming operations where purely elastic streams enter and leave the working zone. The deformations are assumed here to be plane and isochoric. A general solution is formulated in terms of two complex potentials. Residual stress is accounted for in detail and a uniqueness theorem is proved. Some simple flows are examined, but it remains to develop a systematic procedure for matching the general solution to arbitrary boundary data.     

Circularly cylindrical layered media in plane elasticity

This paper presents an alternative efficient procedure to analyze plane elasticity problems of a circularly cylindrical layered media subject to an arbitrary point force. Based on the method of analytical continuation in conjunction with the alternating technique, the elastic fields of the three-phase media are derived. A rapidly convergent series solution which is expressed in terms of an explicit general term of the complex potential of the corresponding homogeneous problem is obtained in an elegant form. As a numerical illustration, the interfacial stresses are presented for different material combinations and for different positions of the point force.     

The derivation of a thick and thin plate formulation without ad hoc assumptions

For the plate formulation considered in this paper, appropriate three-dimensional elasticity solution representations for isotropic materials are constructed. No a priori assumptions for stress or displacement distributions over the thickness of the plate are made. The strategy used in the derivation is to separate functions of the thickness variable z from functions of the coordinates x and y lying in the midplane of the plate. Real and complex 3-dimensional elasticity solution representations are used to obtain three types of functions of the coordinates x, y and the corresponding differential equations. The separation of the functions of the thickness coordinate can be done by separately considering homogeneous and nonhomogeneous boundary conditions on the upper and lower faces of the plate. One type of the plate solutions derived involves polynomials of the thickness coordinate z. The other two solution forms contain trigonometric and hyperbolic functions of z, respectively. Both bending and stretching (or in-plane) solutions are included in the derivation.     

The state vector methods for space axisymmetric problems in multilayered piezoelectric media

The state vector equations for space axisymmetric problems of transversely isotropic piezoelectric media are established from the basic equations. Using the Hankel transform, the state vector equations are reduced to a system of ordinary differential equations. An analytical solution of the problems in the Hankel transform space is presented in the form of the product of initial state vector and transfer matrix. The transfer matrices are given for the three distinct eigenvalues. Applications of the solutions are discussed. An analytical solution for the transversely isotropic semi-infinite piezoelectric media subjected to concerted point loads on the surface z=0 is presented in the Hankel transform space. Using transfer matrix and the continuity conditions at the layer interfaces, the general solution formulation of N-layered transversely isotropic piezoelectric media is given. A selected set of numerical solutions is presented for a layered semi-infinite piezoelectric solid.     

On the method of reciprocal theorem to find solutions of the plane problems of elasticity

In this paper the method of reciprocal theorem is extended to find solutions of plane problems of elasticity of the rectangular plates with various edge conditions. First we give the basic solution of the plane problem of the rectangular plate with four edges built-in as the basic system and then find displacement expressions of the actual system by using the reciprocal theorem between the basic system and actual system with various edge conditions. When only displacement edge conditions exist, obtaining displacement expressions by means of the method of reciprocal theorem is actual. But in other conditions, when static force edge conditions or mixed ones exist, the obtained displacements are admissible. In order to find actual displacement, the minimum potential energy theorem must be applied. Calculations show that the method of reciprocal theorem is a simple, convenient and general one for the solution of plane problems of elasticity of the rectangular plates with various edge conditions. Evidently, it is a new method.     

The penny-shaped crack at a bonded plane with localized elastic non-homogeneity

This paper examines the axisymmetric problem pertaining to a penny-shaped crack which is located at the bonded plane of two similar elastic halfspace regions which exhibit localized axial variations in the linear elastic shear modulus, which has the form G(z)=G₁+G₂e^±ζz. The equations of elasticity governing this type of non-homogeneity are solved by employing a Hankel transform technique. The resulting mixed boundary value problem associated with the penny-shaped crack is reduced to a Fredholm integral equation of the second kind which is solved in a numerical fashion to generate the crack opening mode stress intensity factor at the tip.     

Influence of couple-stresses on stress concentrations around the cavity

ＩｎｔｒｏｄｕｃｔｉｏｎＳｔｒｅｓｓｃｏｎｃｅｎｔｒａｔｉｏｎｉｓｏｎｅｉｍｐｏｒｔａｎｔｐｒｏｂｌｅｍｏｆｍｅｃｈａｎｉｃｓｒｅｓｅａｒｃｈｄｏｍａｉｎ .Ｉｎｔｈｅｍｉｃｒｏｐｏｌａｒｅｌａｓｔｉｃｉｔｙｔｈｅｏｒｙ ,ｉｔｉｓｍｏｒｅａｂｓｏｒｐｔｉｖｅ .Ｉｔｉｓｗｅｌｌ＿ｋｎｏｗｎｔｈａｔｔｈｅｍｉｃｒｏｐｏｌａｒｅｌａｓｔｉｃｉｔｙｔｈｅｏｒｙｕｓｕａｌｌｙｇｉｖｅｓａｌｏｗｃｏｅｆｆｉｃｉｅｎｔａｎｄｃｈａｎｇｅｓｔｈｅｓｔｒａｎｇｅｎｅｓｓａｔｔｈｅｔｉｐｏｆｔｈｅｃｒａｃｋ .Ｔｈｅｓ…     

The stress state of a plate subjected to a piecewise homogeneous load

We consider the stress-strain state of a plate having a doubly connected domain S bounded from the outside by a circle of radius R and from the inside by an ellipse with two rectilinear cuts. The cuts lie symmetrically on the x-axis. The plate is subjected to various forces: the hole contour (the ellipse) is under the action of uniformly distributed forces of intensity q, and the cut shores are free of loads; at the points ±ib of the imaginary axis, the plate is under the action of a lumped force P.The solution of the problem is reduced to determining two analytic functions φ(z) and ψ(z) satisfying certain boundary conditions (depending on the type of the acting loads).We use the Kolosov-Muskhelishvili method to reduce the problem to a system of linear algebraic equations for the coefficients in the expansions of the functions φ(z) and ψ(z). The solution thus obtained is illustrated by numerical examples.     

接触问题应力分析的混合解法

接触问题是一个极其复杂的非线性问题，单独使用数值方法或实验方法求解应力都有一定的困难．有限元计算与平面光弹性实验相结合的混合法是对接触问题进行应力分析简单而有效的途径．