Научная статья на тему 'Solution of boundary value problems of plasticity with the use of conservation laws'

Solution of boundary value problems of plasticity with the use of conservation laws Текст научной статьи по специальности «Математика»

CC BY
51
7
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЗАКОНЫ СОХРАНЕНИЯ / CONSERVATION LAWS / УПРУГОПЛАСТИЧЕСКАЯ ГРАНИЦА / ELASTIC-PLASTIC BOUNDARY

Аннотация научной статьи по математике, автор научной работы — Senashov Sergei I., Savostyanova Irina L., Cherepanova Olga N.

In thepaperconservation lawsofaspecial form for systemsof first-order differentialequationsdepending on two dependent and independent variables are looked at. It is shown how the conservation laws can be used to solve hyperbolic-type and elliptic-type systems of equations that are come across in the theory of plasticity. Examples of an effective use of the described technique are given. With the use of the conservation laws was found the elastoplasticboundaryinaproblemof stress-strain stateofa plate with free-form holes.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Решение краевых задач пластичности с использованием законов сохранения

Вработерассмотрены законы сохранения специального вида для системдифференциальных уравнений первого порядка, зависящие от двух зависимых и независимых переменных. Показано как законы сохранения могут быть использованы для решения систем уравнений гиперболического и эллиптического типов, которые встречаются в теории пластичности. Приведены примеры эффективного применения описанной методики.Спомощью законов сохранения найдена упругопластическая границавзадачеонапряженно деформированномсостоянии пластинысотверстиями произвольной формы.

Текст научной работы на тему «Solution of boundary value problems of plasticity with the use of conservation laws»

УДК 539.374

Solution of Boundary Value Problems of Plasticity with the Use of Conservation Laws

Sergei I. Senashov* Irina L. Savostyanova^

Reshetnev Siberian State University of Science and Technology Krasnoyarsky Rabochy, 31, Krasnoyarsk, 660037

Russia

Olga N.Cherepanova *

Institute of Mathematics and Computer Science Siberian Federal University Svobodny, 79, Krasnoyarsk, 660041

Russia

Received 20.11.2017, received in revised form 26.01.2018, accepted 20.03.2018 In the paper conservation laws of a special form for systems of first-order differential equations depending on two dependent and independent variables are looked at. It is shown how the conservation laws can be used to solve hyperbolic-type and elliptic-type systems of equations that are come across in the theory of plasticity. Examples of an effective use of the described technique are given. With the use of the conservation laws was found the elastoplastic boundary in a problem of stress-strain state of a plate with free-form holes.

Keywords: conservation laws, elastic-plastic boundary DOI: 10.17516/1997-1397-2018-11-3-356-363.

Introduction

Symmetries and conservation laws allow studying of differential equations from different points of view.

Symmetries act on the manifold of solutions of differential equations and that is why allow us to study and find their solutions in two ways. The first way is to find the so-called fixed-points, i.e. such solutions that do not change when some symmetries act. Invariant solutions get found in such a way. This method is realised to the full extent in the papers of L.V.Ovsyannikov [1].

The other approach is frequently declared but is seldom realised. Its main point is that we get symmetries to act on the given solution and we obtain new solutions of the same differential equation system. It was used in the papers [2-5]. There may be an obstruction for realisation of this method which is "poorness" of the admitted group of symmetries or difficulty to interpret the found "multiplied" solutions. We manage to avoid difficulties in interpretation if we act not on the solution itself but on some other objects related to it, for example on characteristics in case of hyperbolic equations. Use of symmetries when solving boundary value problems for

*[email protected] t [email protected] [email protected] © Siberian Federal University. All rights reserved

differential equation systems encounters significant difficulties, though attempts of such a kind have been undertaken, this is evidently explained by local property of symmetries. But on the other hand the technique of using symmetries allows finding and describing fairly wide classes of solutions, and at that fairly effectively.

Conservation laws building technique is in actual fact dual to symmetries building technique and the methods to find them are close to the methods to find symmetries. But the results significantly differ: conservation laws allow solving boundary value problems, and at that fairly effectively. This is illustrated in a series of papers on solving of boundary value problems for plasticity equations and other continuum mechanics equations [6-10].

1. Fundamental definitions

For simplicity we restrict ourselves to the case of differential equations depending on two independent x,y and two dependent variables u, v. This case is simpler for understanding; best investigated and is frequently come across in mechanics problems. Let us look at the system of the two differential equations

Fi = Fi (x,y.u,v.ux,uy, vx,vy) = 0, F2 = F2 (x,y.u,v.ux,uy,vx,vy) = 0, (1)

where the index below stands for a derivative with respect to the corresponding argument.

Definition 1. As the conservation law for system of equations (1) we will call the expression in the form of

DxA + DyB = fii (Fi) + ^2 (F2), (2)

where Q,i i = 1, 2 some linear differential operators that at the same time are not trivial, Dx, Dy are operators of full differentiation for the corresponding variables Dx = dx + uxdu + vxdv + ..., Dy = dy + uy du + vy dv + ....

Remark. This quite a commonplace definition is reasonably suitable for the presentation of our results. More general approach can be learned in [11].

Definition 2. Vector (A, B) is called conserved current.

In this paper functions and are supposed to be dependent on x, y, u, v only. We can assume that and also depend on the derivatives but such conservation laws for system (2) are difficult to interpret and that is why they will not be considered here. From (1) according to Riemann formula it follows

Ady - Bdx = 0, (3)

where C is closed contour in which Functions and have no singularities for simplicity. Let us show how conservation laws (2) can be used for different types of differential equations.

2. Use of conservation laws to solve hyperbolic systems of differential equations

Assume (1) is a system of hyperbolic differential equations. Its characteristics are written as

Fi : dx = a (u,v,x,y), r : = P(u,v,x,y). (4)

These characteristics proceed from the endpoints of the segment of curve L x\ ^ x ^ x2 with conditions established on it

u\L = uo, v\L = vo.

Assume (xp,yp) is the point of intersection of the characteristics (see Fig. 1).

Fig. 1

Let us find this point. By virtue of (3) we have

We have

(5)

/ Ady - Bdx + Ady - Bdx + Ady - Bdx = 0. Jl jti jt2

i Ady - Bdx = i (A - aB)dy = (A - aB)y\vy1 - i yd(A - aB),

Jfi Jfi P Jfi

( Ady - Bdx = f (A - ¡3B)dy = (A - [lB)y\yi - ( yd(A - pB).

jf2 -T2 jf2

Let us assume that and is the solution of equation (2) with the following boundary conditions (A - aB) r = const, (A - ftB) \r2 = const.

In this case from (5) we can obtain coordinate yp. In the same way coordinate xp shall be obtained. In this way we manage to build the characteristics of System of equations (1), and therefore to solve the set problem. Moreover for many of mechanics problems building of the characteristics often give more information to the researcher than finding of an explicit solution. The described method is effective, with the use of which were solved Cauchy and Riemann problems for perfect plasticity equations, and also some problems for other equations of mechanics [12].

3. Use of conservation laws to solve elliptic of differential equation systems

Assume system (1) is of an elliptic type. We will be finding the solution of equation (2) that has a singularity in point (x0, y0). Assume that on closed curve L the functions u\L = u0,v\L = v0 are defined.

Let us encircle point (x0,y0) with a circumference with a radius e : (x — x0)2 + (y — y0)2 = e2. Let us connect it as shown on Fig. 2, with a cut with contour L. We have

/ Ady - Bdx + Ady - Bdx + Ady - Bdx + Ady - Bdx = 0. 'l J ri jF2 Je

Fig. 2

From (6) we obtain

Ady - Bdx = - Ady - Bdx.

Let us set e ^ 0. We obtain a relation relating u (x0,y0) ,v (x0,y0) and J Ady — Bdx. Having chosen another conserved current with the same conditions we get another relation relating u (xo,yo) ,v (x0,y0) and J Ady — Bdx. These two relations allow finding values

u (x0,y0) ,v (x0,y0) and because point (x0,y0) is arbitrary, then therefore the solution for the set problem has been built. This method is realised in papers [8-10]. It allowed calculating of the boundary between plastic strain range and elastic strain range for rolled section rods being twisted. The results of the calculations are given in the Fig. 3.

Fig. 3

L

4. Solution of a plane elastic-plastic problem

Within the plastic strain range in a bidimensional case with a general yield criterion the equations will be written as

dax dr dr da,, , %

—x +--=0,--1--y- = 0, (7)

dx dy dx dy

(ax,a- ,t ) = 0. (8)

Within the elastic strain range relations (7) and compatibility condition of deformations are performed that is written as

dax dr dr day

—x +--=0,--1--y = 0,

dx dy dx dy

d2(°x + ay) + d2 (ax + ay ) = 0

+ O 0.

(9)

dx2 dy

5. Problem setting

Let us consider a finite rectangular plate with dimensions of a x b with holes of an arbitrary shape limited by contours ri, r2,..., rm. Let us assume that the plate is under different stresses in x,y directions, and the contour lines of the holes are under the action of the set stresses. On the external boundaries of r0 we get the following conditions

ax|x=a = ax|x=0 = qi, a-\-=b = a- -=0 = q2. (10)

Other stress components on these boundaries are equal to zero. On the boundary of contours r1, r2,..., rm the condition is fulfilled.

axni + Tmi\ri = pi, Tni + a- mi |r = p2, (11)

where qi,pj are constants, (ni,mi) is vector components of the normal to contour

We assume that stresses ax, ay,t on contours r1, r2,..., rm are such that plasticity condition (8) is fulfilled.

From equation (9) we obtain

ax + a- = F (x,y), (12)

where F (x, y) some solution of Laplace equation (9). From (10) and (11) it follows that in order to find F (x, y) we need to solve Laplace equation with the following boundary conditions

F \ x=a = F |x=0 = qi, F\y=b = F-=0 = qs, F |r = di. (13)

Value di shall be determined from the solution of the three equations (11) and (8). Let us consider the solution of problem (12), (13) as known because this is a classical problem for Laplace equation. In this case in order to find ax,T we will get the system of equations

dax dT , dT dax dF . .

= ~sT + â"=0, $2 ^^ — ^x + 1T=0- (14)

dx dy dx dy dy

For system (14) let us write conservation laws in the form of (2)

DxA + D-B = fii ($i)+fi2 ($2) . Let us find conserved current in the form of:

A = ai (x, y) ax + fîi (x, y) t + 7i (x, y),

B = «2 (x, y) ax + P2 (x, y) t + 72 (x, y) .

As a result we have

dai - ßßi =0 dßi + dai =0 dji + dj2 = ^ 8F_

dx dy ' dx dy ' dx dy dy

Finally conserved current shall be written as

A = ai (x, y) ax + ßi (x, y) t + Yi (x, y), B = -ßi (x, y) ax + ai (x, y) t + Y2 (x, y),

(16)

(17)

coefficients of which are related by relations (16).

To simplify further calculations let us suppose that we have only one hole limited by contour r, (see Fig. 4) according to Green formula we get

Ady - Bdx + f Ady - Bdx = 0. (18)

ro Jr 1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Fig. 4

Let us look at the two singular solutions of equations (16)

, l

x - xo

i (x - xo)2 + (y - yof li = - J 9y "2dy, y1 = 0.

ßi = -

y - yo

(x - xo f + (y - yo)

(19)

2

y - yo

1 (x - xof + (y - yo)2'

2 dF 2

Y2 = - J 9y "2dy, Y i = 0.

ß2 =

x - xo

(x - xo)2 + (y - yo)2 '

From (19) and (20) according to the above shown formulas we finally obtain

2nax(xo, yo) = -

Ü

y - yo

x - xo

ro V (x - xo) +(y - yo)2

(x - xo) +(y - yo)

jt dy-

y - yo

(x - xo)2 + (y - yo)

J&x +

x - xo

(20)

(x - xo)2 + (y - yo)

2 t + Yi

) dx - J Yldy,

2nT(xo, yo) = -

hi, (

x - xo

Tax +

y - yo

IrA (x - xo)2 + (y - yo)2 x (x - xo)2 + (y - yo)2

t dy-

x - xo

— I —

, \2 t \2 x

(x - xo) +(y - yo)

ax +

y - yo

(x - xo)2 + (y - yo)

2 t + Yi

) dx - J Y2dy,

a„ —

Oy(xo,yo) = F(xo,yo) - ax(xo,yo).

We insert the obtained stresses into yield condition (8). In those points where f (ox, oy, t) < 0,

it will be elastic range, and other points of the plate belong to plastic range.

Conclusion

The results of the use of conservation laws to solve boundary value problems have shown that

this method is effective for the problems of the theory of plasticity and elastic plasticity.

References

[1] L.V.Ovsyannikov, Group Analysis of Differential Equations, NewYork, Academic Press, 1982.

[2] S.I.Senashov, A.N.Yakhno, Reproduction of solutions for bi-dimensional ideal plasticity, Journal of Non-Linear Mechanics, 42(2007), 500-503.

[3] S.I.Senashov, A.N.Yakhno, Deformation of characteristic curves of the plane ideal plasticity equations by point symmetries, Nonlinear analysis , 71(2009), 1274-1284.

[4] S.I.Senashov, O.N.Cherepanova, New classes of solution of minimal surfaces, J. Siberian Federal Univ., Math. and Physics., 3(2010), no. 2, 248-255 (in Russian).

[5] S.I.Senashov, Conservation Laws. Hodograph Transformation and Boundary Value Problems of Plane Plasticity, SIGMA, 8(2012), 071.

[6] S.I.Senashov, A.N.Yakhno, Some symmetry group aspects of a perfect plane plasticity system, J. Phys. A: Math. Theor, 46(2013), 355202.

[7] S.I.Senashov, A.N.Yakhno, Conservation Laws of Three-Dimensional Perfect Plasticity Equations under von Mises Yield Criterion, Abstract and Applied Analysis, 2013(2013), 702-732.

[8] S.I.Senashov, A.V.Kondrin, O.N.Cherepanova, On Elastoplastic Torsion of a Rod with Multiply Connected Cross-Section, J. Siberian Federal Univ., Math. and Physics., 7(2015), no. 1, 343-351.

[9] S.I.Senashov, O.N.Cherepanova, A.V.Kondrin, Elasto-plastic Bending of Beam, J. Siberian Federal Univ., Math. and Physics., 7(2014), no. 2, 218-223.

[10] S.I.Senashov, E.V.Filyushina, O.V.Gomonova, Construction of elasto-plastic boundaries using conservation laws, Vestnik SibGAU , 16(2015), no. 2, 343-359.

[11] S.I.Senashov, A.M.Vinogradov, Symmetries and Conservation Laws of 2-Dimensional Ideal Plasticity, Proc. of Edinb. Math. Soc, 31(1988), 415-439.

[12] S.I.Senashov, P.P.Kiriakov, A.N.Yakhno, Application of symmetries and conservation laws of solution of Differential Equation, Novosibirsk, Nauka, 2001 (in Russian).

Решение краевых задач пластичности с использованием законов сохранения

Сергей И. Сенашов Ирина Л. Савостьянова

Сибирский государственный университет науки и технологий им. М. Ф. Решетнева

Красноярский рабочий, 31, Красноярск, 660037

Россия

Ольга Н. Черепанова

Институт математики и фундаментальной информатики Сибирский федеральный университет Свободный, 79, Красноярск, 660041

Россия

В 'работе рассмотрены законы сохранения специального вида для систем дифференциальных уравнений первого порядка, зависящие от двух зависимых и независимых переменных. Показано как законы сохранения могут быть использованы для решения систем уравнений гиперболического и эллиптического типов, которые встречаются в теории пластичности. Приведены примеры эффективного применения описанной методики. С помощью законов сохранения найдена упругопла-стическая граница в задаче о напряженно деформированном состоянии пластины с отверстиями произвольной формы.

Ключевые слова: законы сохранения, упругопластическая граница.

i Надоели баннеры? Вы всегда можете отключить рекламу.