CC BY
6
УДК 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
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
Россия
В 'работе рассмотрены законы сохранения специального вида для систем дифференциальных уравнений первого порядка, зависящие от двух зависимых и независимых переменных. Показано как законы сохранения могут быть использованы для решения систем уравнений гиперболического и эллиптического типов, которые встречаются в теории пластичности. Приведены примеры эффективного применения описанной методики. С помощью законов сохранения найдена упругопла-стическая граница в задаче о напряженно деформированном состоянии пластины с отверстиями произвольной формы.
Ключевые слова: законы сохранения, упругопластическая граница.
