Solution of the first boundary value problem of plane elasticity theory using conservation laws

封面

如何引用文章

详细

A huge number of works are devoted to solving boundary value problems for the equations of plane elasticity theory. The largest number of studies in this area are based on the formula found by G. V. Kolosov. He was the first to express the general solution to the problem of plane elastic deformation by finding two independent functions of a complex variable. This made it possible to apply a well-developed theory of analytic functions to solving problems of elasticity theory. Later, the solution method based on Kolosov's formula was developed by his student N. I. Muskhelishvili. But the described method also has significant limitations. It is applicable only to those areas that can be conformally mapped onto a circle. Therefore, other methods for solving elasticity theory problems are also needed, since a large number of practically important problems are solved for areas that do not satisfy this condition. The method developed in the work is based on the use of conservation laws that are constructed for equations describing a plane deformable state. The assumptions made in the work make it possible to construct a solution to the first boundary value problem for arbitrary plane areas bounded by a piecewise smooth contour. In this case, finding the components of the stress tensor is reduced to calculating contour integrals along the boundary of the region under consideration. As in the case considered by G. V. Kolosov, the solution to the problem is based on two exact solutions of the Cauchy – Riemann equations, which have singularities at an arbitrary point in the region under consideration.

全文:

Introduction

The solution to boundary value problems for equations of elasticity theory in a plane stationary case is described in a huge number of articles and monographs. The classic work in this direction is a book written by G. V. Kolosov’s student [1]. Despite the long history of solving such problems, the interest in their solution has not waned. This is due to the fact that the classical formulae by G. V. Kolosov allow solving equations of elasticity theory not for all boundary value problems arising in science and technology. The main limitation is caused by the smoothness of the boundary and some features of complex variable functions. The other methods involved in decomposing the desired functions into series for different types of special functions also have natural limitations, which are related to the similarity of the series being used as well as the complexity of the results obtained.

Let us describe some results of the recent studies on the theory of elasticity. The paper [2] presents a brief historical review of the studies devoted to the bending theory of elastic plates. In [3] the author considers the problem of detecting and identifying elastic inclusions in an isotropic linearly elastic plane. The article [4] considers a complicated version of the well-known Lame problem posed in 1852 describing the solution to the static equilibrium of a parallelepiped with free lateral surfaces exposed to the action of opposite end forces, and also for the case of impact effects of end forces. In paper [5] the construction of fundamental solutions for harmonic oscillations equations in elastic anisotropic elastic media theory is carried out; the fundamental solution to oscillations equations for isotropic medium in closed form is constructed. The work [6] presents a general solution to the problems of the elastic theory for anisotropic half-planes and bands with arbitrary holes and cracks using complex potentials of the plane problem of the elastic theory of an anisotropic body. The article [7] is devoted to the study of the functions of tensions, which allow satisfying exactly the equations of equilibrium of the classical theory of elasticity and to obtain a solution in tensions. In order to obtain dependencies between voltages and voltage functions, a mathematical apparatus of general relativity is used. In [8] terms of complex-significant displacements, the system of equations of the axisymmetric theory of elasticity is written; the fundamental solution to this system of equations is a general representation of the field of displacements in an axisymmetric case, similar to the formulae by Kolosov-Muskhelishvili in a plane problem. The basic equations of linear moment elasticity theory are presented in [9]. The determining ratios are written for the case of general anisotropy in the form of linear equations. Some simplified options are considered, in particular with constricted rotation and plane deformation when only shear stresses are present. The work [10] is devoted to the analysis of a boundary value problem with an unknown area of contact, describing the equilibrium of two-dimensional elastic bodies with a thin slightly curved web. The paper [11] explores the evolution of the wave pattern in a multi-modular elastic half-space with a non-stationary one-axis piecewise-linear motion of its boundary in the mode “tension – compression – stop”.

In this work the laws of conservation of differential equations of elasticity are being used. The laws, which allow reducing the presence of tensor component of the tensions at the point to the contour integral along the boundary of the area being considered, have been presented. In these conditions the boundary of an area requires only partial smoothness. Note that earlier some conservation laws were presented in works [12; 13] but they were not used to solve any problems.

Statement of the problem

Let us consider the equations describing plane elastic deformation.

The relations between strain tensor components and displacement vector components in the case of small deformations have the following form

εx=w1x,εy=w2y,εxy=w1y+w2x.                                                        (1)

The Hooke’s will be written as follows:

εx=1E(σxνσy),εy=1E(σyνσx),εxy=2(1+ν)Eτ.                                             (2)

Strain compatibility conditions:

2εx2y+2εy2x=2εxyxy.                                                                                        (3)

By putting (2) into (3), we obtain

2(σxνσy)2y+2(σyνσx)2x=2(1+ν)2τxy                                                  (4)

Equilibrium equations:

σxx+τy=0,τx+σyy=0.

From hence we obtain

2σy2y+2σx2x=22τxy.                                                                             (5)

From (4) and (5) we obtain

Here are strain tensor components;  are voltage tensor components;  are displacement vector components;  are elastic constants.

The system has the following final form:

σxx+τy=0,τx+σyy=0,    Δ(σx+σy)=0.                                       (6)

Let us put the first boundary value problem for the system (6):

σxn1+τn2|L=X(x,y),τn1+σyn2|L=Y(x,y).                                               (7)

Here  are components of an outer normal vector to a piecewise smooth contour limiting the end area

Let us find the solution to the problems (6) and (7) in the form

σx+σy=pconst, p ≠0.                                                                              (8)

Let us introduce new variables:

u=σxp,v=τp,σy'=σyp,f=Xp,g=Yp.                                            (9)

Then the problem (6), (7) is written in the following form

F1=ux+vy=0,F2=uyvx=0,                                                                 (10)

un1+vn2=f,vn1un2=gn2,                                                                 (11)

here and then the index at the bottom means a derivative of this variable.

Thus, it is necessary to solve the boundary value problem (11) for the system of equations (10) with the use of conservation laws.

Conservation laws of the system of equations (10)

Definition. The conservation law for the system of equations (10) will be called the following expression

Ax+By=ω1F1+ω2F2,                                                                          (12)

where  are linear differential operators that are simultaneously nonzero identical,

A=α1u+β1v+γ1,B=α2u+β2v+γ2,                                                   (13)

 are some smooth functions that only depend on x,y.

Remark. A more general definition of the conservation law suitable for the arbitrary systems of equations can be found in [14; 15].

From (12) taking into account (13) we obtain

αx1u+α1ux+β1xv+β1vx+γ1x+αy2u+α2uy+βy2v+β2vy+γy2=ω1(ux+vy)+ω2(uyvx)=0.      (14)

From (14) it follows

αx1+αy2=0,βx1+βy2=0,α1=ω1,β1=ω2,α2=ω2,β2=ω1,γ1+γ2=0.

From hence we obtain

α1=β2,α2=β1.                                                                                       (15)

Therefore

αx1βy1=0,αy1+βx1=0.                                                                              (16)

From the given formulae it follows that the system of equations (10) allows infinitely many laws of conservation. Only those that can solve the problem will be listed below.

Therefore, the conserved current is:

A=α1u+β1v+γ1,B=β1u+α1v+γ2.

From (13) we obtain

S(Ax+By)dxdy=LAdy+Bdx=0,                                                 (17)

where S is the area bounded by the L curve.

Solution to the problem (10), (11)

In order to find the  values within the S area, it is necessary to construct solutions to the Cauchy – Riemann system (16) having singularities at the arbitrary point (x0,y0)S.

The first of these solutions has the folloing form:

α1=xx0(xx0)2+(yy0)2,β1=yy0(xx0)2+(yy0)2,γ1=γ2=0.     (18)

Remark. If mass forces are included in the equilibrium equation, then γ1,γ2 will no longer be equal to zero.

At the point (x0,y0)S the functions α1,β1 have singularities, thus let us circle this point with the circumference ε:(xx0)2+(yy0)2=ε2.

Then we obtain from the formula (17)

LAdy+Bdx+εAdy+Bdx=0,                                                      (19)

Let us calculate the second integral in the formula (19). We have

Ady+Bdx=ε(u(xx0)(xx0)2+(xx0)2v(yy0)(xx0)2+(xx0)2)dy+

+u(yy0)(xx0)2+(yy0)2v(xx0)(xx0)2+(yy0)2dx.

Let us introduce the new coordinates xx0=εcosφ,yy0=εsinφ, we obtain

εAdy+Bdx=02π[(ucosφ+vsinφ)cosφ(usinφ+vcosφ)sinφ]dφ==02πudφ=2πu(x0,y0).              (20)

The last equality is obtained by the mean-value theorem at ε0.

For the final construction of the solution, let us find the values u,v on the L boundary. From the formula (11) we obtain

u=fn1gn2+n22,v=fn2+gn1+n1n2.                                                            (21)

Let us put (21) into (20); and taking into account (19) we obtain

2πu(x0,y0)=2πσx(x0,y0)/p==L((fn1gn2+n12)(xx0)(xx0)2+(yy0)2(fn2+gn1n1n2)(yy0)(xx0)2+(yy0)2)dy+(fn1gn2+n12)(yy0)(xx0)2+(yy0)2+(fn2+gn1n1n2)(xx0)(xx0)2+(yy0)2dx.

The second solution to the system of equations (16) is taken as

α1=yy0(xx0)2+(yy0)2,β1=xx0(xx0)2+(yy0)2,                                               (22)

Having made calculations similar to those made with the solution (18), we obtain

2πv(x0,y0)=2πτ(x0,y0)/p==L((fn1gn2+n12)(yy0)(xx0)2+(yy0)2+(fn2+gn1n1n2)(xx0)(xx0)2+(yy0)2)dy+(fn1gn2+n12)(xx0)(xx0)2+(yy0)2+(fn2+gn1n1n2)(yy0)(xx0)2+(yy0)2dx.

Conclusion

This paper proposes a new method for solving the first boundary value problem for the equations of plane elasticity theory in a stationary case. This method makes it possible to find the value of the component of a strain tensor at each point of the area being studied. In this case, the stress calculations are limited only to the calculation of contour integrals along the boundaries of the area.

×

作者简介

Olga Pashkovskaya

Reshetnev Siberian State University of Science and Technology

编辑信件的主要联系方式.
Email: pashkovskaya@sibsau.ru
ORCID iD: 0009-0003-2529-4105

Cand. Sс., associate Professor, Reshetnev

俄罗斯联邦, 31, Krasnoyarskii rabochii prospekt, Krasnoyarsk, 660037

Sergei Lukyanov

Reshetnev Siberian State University of Science and Technology

Email: lukyanovsv@sibsau.ru

postgraduate student

俄罗斯联邦, 31, Krasnoyarskii rabochii prospekt, Krasnoyarsk, 660037

参考

  1. Muskhelishvili N. I. Nekotorye osnovnye zadachi matematicheskoy teorii uprugosti [Some basic problems of the mathematical theory of elasticity]. Moscow, Nauka Publ., 1966, 708 p.
  2. Vasil'ev V. V. [Theory of thin elastic plates – history and current state of the problem]. Izv. RAN. MTT. 2024, No. 2, P. 3–39 (In Russ.).
  3. Kaptsov A. V., Shifrin E. I. [Plane problem of elasticity theory on identification of nodal points of a quadrature inclusion]. Izv. RAN. MTT. 2023, No. 6, P. 47–68 (In Russ.).
  4. Rasulova N. B., Makhmudzade T. M. [Solution of the dynamic Lame problem]. Izv. RAN. MTT. 2023, No 5, P. 131–137 (In Russ.).
  5. Il'yashenko A. V. [Fundamental solutions of the equations of the theory of oscillations for anisotropic elastic media]. Izv. RAN. MTT. 2023, No. 5, P. 138–146 (In Russ.).
  6. Kaloerov S. A., Glushankov E. S., Mironenko A. B. [Solution of elasticity theory problems for multiply connected half-planes and strips]. Izv. RAN. MTT. 2023, No. 4, P. 23–37 (In Russ.).
  7. Vasil'ev V. V., Fedorov L. V. [Stress functions in elasticity theory]. Izv. RAN. MTT. 2022, No. 4, P. 103–113 (In Russ.).
  8. Georgievskiy D. V., Stetsenko N. S. [Aleksandrovich's Complex Representation of Solutions in Displacements in Three-Dimensional Elasticity Theory]. Izv. RAN. MTT. 2022, No. 3, P. 8–15 (In Russ.).
  9. Annin B. D., Ostrosablin N. I., Ugryumov R. I. [Constitutive equations of anisotropic linear moment theory of elasticity and a two-dimensional pure shear problem with constrained rotation]. Sibirskiy zhurnal industrial'noy matematiki. 2023, Vol. 26, No. 1(93), P. 5–19 (In Russ.).
  10. Khludnev A. M. [On the equilibrium of elastic bodies with a slightly curved bridge]. Sibirskiy zhurnal industrial'noy matematiki. 2023, Vol. 26, No. 3(95), P. 154–168 (In Russ.).
  11. Dudko O. V., Lapteva A. A., Ragozina V. E. [Interaction of plane deformation waves in a multimodular elastic half-space at the stage of forced stopping of its boundary after uniaxial tension-compression]. 2023, Vol. 26, No. 4(96), P. 32–48 (In Russ.).
  12. Olver P. Conservation laws in elasticity. I General result. Arch. Rat, Mech Anal. 1984, No. 85, P. 111–129 (In Engl.).
  13. Senashov S. I., Filyushina E. V. [Conservation laws of plane elasticity theory]. Vestnik SibGAU. 2014, Vol. 1 (53), P. 79–81 (In Russ.).
  14. Vinogradov A. M. Local symmetries and conservation laws. Acta Appl. Math. 1984, No. 6, P. 56–64.
  15. Senashov S. I., Vinogradov A. M. Symmetries and Conservation Laws of 2-Dimensional Ideal Plasticity. Proc. of Edinb. Math. Soc. 1988, No. 31, P. 415–439.

补充文件

附件文件
动作
1. JATS XML

版权所有 © Pashkovskaya O.V., Lukyanov S.V., 2025

Creative Commons License
此作品已接受知识共享署名 4.0国际许可协议的许可