Abstract and Applied Analysis

Volume 2013, Article ID 702132, 8 pages

http://dx.doi.org/10.1155/2013/702132

## Conservation Laws of Three-Dimensional Perfect Plasticity Equations under von Mises Yield Criterion

^{1}Siberian State Aerospace University, 31 Krasnoyarsky Rabochy Avenue, Krasnoyarsk 660014, Russia^{2}Departamento de Matemáticas, CUCEI, Universidad de Guadalajara, Boulevard Marcelino García Barragán 1421, 44430 Guadalajara, JAL, Mexico

Received 4 July 2013; Revised 22 September 2013; Accepted 26 September 2013

Academic Editor: A. F. Cheviakov

Copyright © 2013 S. I. Senashov and A. Yakhno. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

For the first time the conservation laws for von Mises plasticity equations in three dimensions as well as for plane stress equations are given. In the plane case conservation laws are used to construct characteristics for the Cauchy problem. For the system of the plane strain, the conservation laws are used to solve the free boundary problem for any convex smooth contour loaded with constant normal and zero tangential stresses.

#### 1. Introduction

Conservation laws are becoming one of the most important tools for studying and solving differential equations. Such significance was gained after the article by Noether [1]. In this paper it was shown that there is a close relationship between the so-called Noether’s symmetries and conservation laws. Namely, the conservation law corresponds to each such symmetry.

Later it was shown that such a relationship exists only for the equations derived from the variational principle. For other equations the so-called operator of universal linearization should be constructed. And then it is necessary to consider the kernel of its (formally) adjoint operator [2, 3]. In particular, it permits consistently to construct the conservation laws of plasticity equations and moreover to use them to solve the main boundary value problems of the plane theory of perfect plasticity [4, 5].

Let us give some basic definitions of the theory of conservation laws due to [2]. Let us consider a system of differential equations for unknown functions and its derivatives , up to the order with respect to independent variables ().

A *conserved current* is an -dimensional vector function with components (*fluxes*) depending on , , , which satisfy the following relation for any solution of :
where are some scalar differential operators. The above relation (2) is called *conservation law* of system . Let us note that all conserved currents are considered to be equivalence classes with respect to *trivial* currents (for which for any form of ). Hereafter the summation over repeating indices is assumed.

There follows (see [4] and cited bibliography) the following theorem.

Theorem 1. *Let system of differential equations (1) be derived from a variational principle and admit the following point Noether’s symmetry (; ):
**
Then theth component of conserved current has the following form:
**
where is the density of the corresponding Lagrangian, .*

Let us recall [2] that a functional
is called a *Lagrangian* (or an *action* or a *variational functional*), and function is called the *density* of .

#### 2. Conservation Laws of Three-Dimensional Plasticity Equations

Let us consider [6] the system of perfect isotropic plasticity in the three-dimensional stationary case which describes the state of the plastic flow of an incompressible plastic medium and consists of three equilibrium equations: the von Mises yield criterion which is defined when the medium achieves a plastic state the incompressibility equation and the relation between components of the deviatoric stress and components of the symmetric strain-rate tensor: where is a Cartesian coordinate system, is the hydrostatic pressure,are the components of the velocity vector,is the yield point, is a positive function defined by yield criterion, are components of stress tensor, is the Kronecker delta, and , . Eliminating and from (6)–(9) we come to four nonlinear equations for functions , , , only: where .

System (10) admits Lie algebra of point transformations spanned by the following generators [7] (): Generators , , , and are obtained from , by cyclic permutation of indices. The group of point transformations is generated by the following monoparametric subgroups: where , , , and , are group parameters.

Transformations (12) are translations with respect to , , ; (13) are scaling on , ; (14) are three rotations in the planes , , and and there are three rigid displacements (15).

System (10) can be deduced from the Lagrangian (): where the density of the Lagrangian has the form and it represents the velocity of dissipation of mechanical energy in elemental volume.

It is easy to verify that (10) are the Euler-Lagrange equations corresponding to which look as follows: where .

Let us note that the von Mises yield criterion (7) was obtained from a variation principle: actual stresses in plastic state maximize the velocity of surface forces work [8] for an incompressible medium. For such a medium quantity is called the intensity of the share strain rate [6].

For system (10) one can construct the conservation laws with the help of Noether’s symmetries. It is well known that Noether’s symmetries are a subalgebra of Lie algebra of point symmetries admitted by . It is sufficient from symmetries (11) to determine the ones that do not change the density (17). To verify if the given point symmetry is a Noether’s one, let us use the known condition of Lagrangian density invariance [9]
Here is the first prolongation of of the form (3), where are the coefficients of , are the coefficients of , and is the coefficient of ; is the *total derivative* operator along :

The direct verification of (19) for symmetries (11) shows that , , , and are point Noether’s symmetries of (17). Using Theorem 1 and taking into account that one can obtain the fluxes of conserved current, corresponding to generator () and the conservation law represents the conservation of generalized momentum along the -axis and reflects the homogeneity of the medium.

For generator , the components of the conserved current look as follows: and represent the conservation law of angular momentum.

For we have the following fluxes: and the conservation law is just the first equilibrium equation from (6) which was obtained from conservation of momentum. By analogy, the conservation laws for generators and correspond to the second and the third equilibrium equations.

Generator produces By analogy, one can construct the conserved currents for , , .

It is easy to see that the linear combination is a point Noether’s symmetry of (17). Then, the components of corresponding conserved current from (4) are ():

Let us note that generator is not a point Noether’s symmetry of (17), because relation (19) for does not vanish:
where . In such a case, Noether’s symmetry is called the symmetry *of divergence type* and it produces [9] the conservation law of the form , where and are given by (4). So, for we obtain and which represents the conservation law of the mass given by (8).

For Navier-Stokes equations the complete set of conservation laws was found in [10]. Due to relative analogy of plastic flow system and Navier-Stokes equations it is quite difficult to expect conservation laws other than those specified here.

#### 3. Conservation Laws of Plane Stress Equations

A plane stress state is approximately achieved in a thin lamina deformed under the action of forces which lie in its median plane. Equations for stresses in a Cartesian plane can be written as follows [6]: where the angle is connected with the value of the mean pressure : and is the angle between the first principal direction of stress tensor and the -axis.

To construct the conservation laws let us use the method proposed in [2]. Firstly it is necessary to construct the matrix differential operator of the *universal linearization * (or the Fréchet derivative of [11]). Then considering the formally adjoint operator one needs to solve the equation
restricted to the solutions of (28). Here is called *the generating function* of the conservation law. Let us recall that operator has the form of matrix
where , being a multiindex, , and . And if , then .

After some manipulations (29) takes the following form: where the derivatives, for example, and , should be changed by the corresponding expressions due to (28). It is clear that (31) have an infinite number of solutions. Let us consider one of them, namely, the one of the form , .

Let us note that it is possible to determine other solutions of (29), but as far as the authors know for the conservation laws depending on the derivative there is no application to solve any practical problems.

Relations (31) should be valid for any solution of system (28), and therefore the coefficients of derivatives , are equal to zero identically. It gives two linear equations to determine the components of conserved current : corresponding to generating function with , . It is easy to verify, substituting in (31) and taking into account the above system and its differential consequences. In such a case, the differential operators from (2) are just two functions:

Finally, the following theorem is valid.

Theorem 2. *The system of plane stress state admits infinite series of conservation laws.*

System (28) has two distinct families of real characteristic curves in the hyperbolic region (if ), defined by equations [6]: with relations along characteristics where where is the angle of intersection of characteristic curves. Let us introduce two new unknown functions and (characteristic coordinates or Riemann invariants) in such a way: , .

Then system (28) takes the form

Finally, for the components of conserved current we have the following linear equations:

Let functions be given along the smooth curve (see Figure 1). Curve is not a characteristic and it intersects with each characteristic only once. Let be a point of intersection of characteristics . going from point and . going from point . Then it is necessary to define a solution of the Cauchy problem (37), (39) in curvilinear triangle [12].

Let us determine the coordinates of point . Taking an integral over the closed path which due to the Green theorem is equal to zero we have Then, by integrating parts we obtain Without loss of generality let us take Then we have By analogy if we put Finally, from (40) we have where components of the conserved current , are defined from the linear system (38) and boundary conditions (42), (45).

For the second coordinate of we have Taking we obtain if we take Finally, from (40) one can determine the second coordinate of point where components of conserved current , are solutions of (38) with boundary conditions (48), (50).

The solution of the aforementioned problems for conserved currents gives the solutions of the Cauchy problem for the initial system, because at point one can reconstruct the values of functions and from boundary conditions (39) and, finally, functions and .

In such a way it is possible to determine the families of characteristic curves for the system of the plane stress. These curves do not coincide with the so-called slip lines (lines where the tangent stress achieved its maximum value) as in the case of the plane strain system considered later. But principal directions of the stress tensor bisect the angles between characteristic curves. That is why it is possible to reconstruct slip lines from the known characteristic curves field.

#### 4. Conservation Laws of Plane Strain Equations

Let us consider the system of perfect plane plasticity with Tresca-Saint-Venant-Mises yield condition [6]: where is the hydrostatic pressure, is the angle between the first principal direction of stress tensor and the -axis, and is the plasticity constant. System (52) describes the plain strain state of a plastic deformed medium. In such a state the displacements of particles of the homogeneous and isotropic body are parallel to plane and are independent of the third coordinate .

Quasilinear system (52) is of the hyperbolic type and characteristics are given by the following equations: Along corresponding characteristics the following relations are valid:

The complete set of conservation laws of (52) as a solution of (29) is described in detail in [13]. In [5], the conserved currents depending only on , are applied to solve main boundary problems and are similar to (38) when , . Thus, the solution of problem (38), (42), (45) has a form where function looks like the following: Accordingly, the solution of problem (38), (48), (50) is where is the modified Bessel function of the first kind of a zero order having the following properties:

Let us consider the problem of a free boundary. Let the constant normal and null tangent components of stresses be given along a smooth contourof some convex cavity (see Figure 1), situated in an infinite plastic deformed medium: It is necessary to determine the free boundary along which The relation between functions , of (52) and normal components , looks as follows [6]: where is the angle between the normal to the contour and -axis.

Substituting (59) in (61), we have along () Let us draw two characteristics where , are the values of at points , respectively.

Applying the above method we define the coordinates of where we have On the other hand, given (60) we obtain Equating the right-hand sides of (64), (65) we have

As a particular example, let us consider contour in the form of an ellipse (Figure 1). Equality (66) is valid, if . Changing the coordinates of points , , one can find curve , which will be the boundary free of stresses. In Figure 2 some points of the free boundary for the contour are given for .

#### 5. Conclusions

For all point symmetries of a three-dimensional perfect plasticity system with the von Mises yield criterion the corresponding conservation laws are calculated. Some of them (conservation of the mass and of the impulse) are in the basis of this system; however, there are new conservation laws.

Recently a rigid plasticity constitutive model with the linear kinematic hardening has been analyzed in [14]. Focusing on planar simple shear the constitutive functions were classified according to continuous symmetry groups. It will be interesting to consider the conservation laws of the obtained systems.

How two linear systems for the components of conserved currents for the system of plane stress state can be used to determine its characteristic curves is shown.

For the perfect strain plane plasticity with Tresca-Saint-Venant-Mises yield condition, the problem of a free boundary for arbitrary convex cavity situated in the infinite plastic deformed medium is solved with the use of the corresponding conservation laws.

#### Acknowledgments

The work was supported by the Ministry of Education and Science of the Russian Federation (Project 1.3720.2011) to S. I. Senashov and by PRO-SNI 2013-UDG to A. Yakhno.

#### References

- E. Noether, “Invariante variations probleme,”
*Nachrichten der Königlichen Gesellschaft der Wissenschaften zu Göttingen*, pp. 235–257, 1918. View at Google Scholar - I. S. Krasil'shchik and A. M. Vinogradov, Eds.,
*Symmetries and Conservation Laws for Differential Equations of Mathematical Physics*, vol. 182 of*Translations of Mathematical Monographs*, American Mathematical Society, Providence, RI, USA, 1999. View at MathSciNet - A. M. Vinogradov, “Local symmetries and conservation laws,”
*Acta Applicandae Mathematicae*, vol. 2, no. 1, pp. 21–78, 1984. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. I. Senashov and A. Yakhno,
*Aplicación de Simetrías y Leyes de Conservación a la Resolución de Ecuaciones Diferenciales de Mecánica*, Universidad de Guadalajara, 2008, (Spanish). - S. I. Senashov and A. Yakhno, “Conservation laws, hodograph transformation and boundary value problems of plane plasticity,”
*SIGMA*, vol. 8, no. 071, 2012. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. M. Kachanov,
*Fundamentals of the Theory of Plasticity*, Dover Books on Engineering, Dover, 2004. - S. I. Senashov, “Invariant solutions of a three-dimensional ideal plasticity problem,”
*Journal of Applied Mechanics and Technical Physics*, vol. 21, no. 3, pp. 417–420, 1980. View at Google Scholar - A. M. Freudenthal and H. Geiringer, “The mathematical theories of the inelastic continuum,” in
*Elasticity and Plasticity*, S. Flügge, Ed., vol. 3–6 of*Encyclopedia of Physics*, pp. 229–433, Springer, Berlin, Germany, 1958. View at Google Scholar · View at MathSciNet - L. V. Ovsiannikov,
*Group Analysis of Differential Equations*, Academic Press, New York, NY, USA, 1982. View at MathSciNet - V. N. Gusyatnikova and V. A. Yumaguzhin, “Symmetries and conservation laws of Navier-Stokes equations,”
*Acta Applicandae Mathematicae*, vol. 15, no. 1-2, pp. 65–81, 1989. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. J. Olver,
*Applications of Lie Groups to Differential Equations*, vol. 107 of*Graduate Texts in Mathematics*, Springer, New York, NY, USA, 1986. View at Publisher · View at Google Scholar · View at MathSciNet - B. L. Roždestvenskiĭ and N. N. Janenko,
*Systems of Quasilinear Equations and Their Applications to Gas Dynamics*, vol. 55 of*Translations of Mathematical Monographs*, American Mathematical Society, Providence, RI, USA, 1983, Translated from the Second Russian Edition by J. R. Schulenberger. View at MathSciNet - S. I. Senashov and A. M. Vinogradov, “Symmetries and conservation laws of 2-dimensional ideal plasticity,”
*Proceedings of the Edinburgh Mathematical Society*, vol. 31, no. 3, pp. 415–439, 1988. View at Publisher · View at Google Scholar · View at MathSciNet - A. F. Cheviakov, J. F. Ganghoffer, and R. Rahouadj, “Finite strain plasticity models revealed by symmetries and integrating factors: the case of Dafalias spin model,”
*International Journal of Plasticity*, vol. 44, pp. 47–67, 2013. View at Publisher · View at Google Scholar