The Phenomenological Description of Plastic Flow in Solids

In the framework of the phenomenological scheme a self-consistent description of the transition from the solid state to the plastic flow was presented taking into account the point defects such as interstitials and vacancies. On the basis of the system of synergistic equations the dependencies of the internal stresses and curvature of the velocity profile of the shear displacement on the order parameter as well as the stationary distribution of the concentration of vacancies were found.


Introduction
The processes occurring during the deformation of solids have attracted attention for a long time [1].Moreover it is well known that the elastic deformation vanishes after the removal of the applied stress, and plastic deformation is retained even after removal of all external influences.Exception is the viscous flow of the material, which can continue indefinitely even when exposed to an arbitrarily small external stress (and the flow rate increases with stress).According to [2] it turns that some properties (elasticity and viscosity) are not unique to solids, but they are inherent in continuous media (e.g., granular materials, supercooled liquid polymers, suspensions, some gels, emulsions, etc.).In this regard, the science of the deformation and flow of continuous media (rheology) has received a great development recently [2,3].Particularly the plastic deformation should be noted: on the one hand, it is typical for many natural and technological processes and, on the other hand, it is difficult in terms of the theoretical description.Usually the plastic deformation of crystals is associated with the emergence and movement of the dislocations within the grains and can be implemented by slipping and twinning [1].Although many methods are applied to the numerical simulation of the dynamics of dislocation lines [4][5][6][7][8], the description of the plastic deformation of crystals remains a difficult task.In the case of amorphous materials and continuous media the analytical description of the plastic flow is possible only within the framework of the phenomenological theories, where the material is considered as homogeneous.The purpose of this paper is to construct a self-consistent phenomenological scheme, which allows describing the transition of solids to the plastic flow state.Since the plastic flow always arises at the microscopic level, that is, at the level of elementary point defects of different physical nature, the influence of point defects on the behavior of solids under the action of the given stress will be an actual task.And although at the moment the theories, associated with the movements of linear defects, have received the increasing development, nevertheless, in some cases (e.g., the return of the metal at low temperatures from the state of strain hardening) only interstitials and vacancies will be mobile, in contrast to the dislocation.
In this paper we propose a phenomenological scheme, within framework of which the transition from a solid state to a plastic flow is presented as a self-organization process of the elementary particles of a material taking into account the ensemble of vacancies and the external influence.Section 2 discusses the main variables that characterize the state of the system.The self-consistent system of equations describing the behavior of a solid is presented in Section 3. The analysis of the system in the framework of the adiabatic approximation is carried out in Section 4. Stationary flow regime is studied in Section 5.

Concentration of Vacancies
Let us consider the dynamics of the material in the case of shear strain, when the displacement is given by the function u(r, ).Usually, solid density deviation  can be associated with the expression where  is an average density of the material, || ≪ ;  trace of the strain tensor [1].
The understanding of the basic laws of plastic flow can be achieved in the framework of the hydrodynamic theory [9,10], where the auxiliary variable was introduced The last one meets the local fraction of the free volume.On the other hand, the parameter  for crystalline solids can be interpreted as the concentration of vacancies and is given by the equation where  0 (r, ) is the density of lattice sites, which depends on both the temperature and the strain and (r, ) is the number of particles per unit of material volume (it depends on the strain) [9].In the case of amorphous solids the concept of the "lattice sites" and "vacancies" or "interstitial atom" is not so obvious, and it is necessary to give a more general interpretation of the density  0 (r, ).Further, we will base on the assumption that the amorphous sample must have a set of periodically arranged in the space positions ("sites"), in each of which the atom can be in local equilibrium for a long time.
For crystals  0 could be treated as density of the particles, which would have a system under an infinite compression pressure at a fixed value of the lattice constant.Consequently for the amorphous solids parameter  0 can be defined as the particle density under an infinite compression pressure at a constant distance to the nearest neighbors [10].Of course, it is difficult to realize experimentally the compression at a constant distance to the nearest neighbors, but, the accurate estimate of  0 can be obtained, if the compression is performed at constant temperature and the distance to the nearest neighbor is properly defined.As a result, to determine  0 at a given pressure  and temperature  with reasonable accuracy one can use the approximation [10] where (  , ) is a density of particles at pressure   and temperature  and   (, ) is the distance between nearest neighbors, determined by X-ray or neutron diffraction.It turns out that the value of  will be critical for the description of the plastic deformation of solids.Indeed, in the absence of vacancies ( =  0 ,  = 0) the material is in a solid state when the external stress gives rise to only elastic deformation; in the case of vacancies ( <  0 ,  ≲ 1) we get the state when the external stress leads to plastic flow of the material.

Synergetic System of Equations
The dynamic equations for the density and momentum density in the deterministic case can be represented in a standard way [8,11] where k = u/ is the rate of change of the shear displacement,  is the elastic stresses in the sample,  0 is the dynamic viscosity of the material, and the operator ∇ = /r.Equations ( 5)-( 6) have the same form as the hydrodynamic equations for one-component fluids [12].
Using the definition of , the relationships for the velocity k and stress with the free energy functional  = {, u} ( 0 is a kinetic coefficient), and the identities one can obtain [13] the equation of motion for the parameter instead of (5).
Parameter  is convenient to be chosen as an order parameter, which distinguishes the solid state and the plastic flow state.Then, to find the variational derivative in (9) we write the free energy of the system in the form of the Landau expansion for the second order phase transition where ,  are the Landau expansion coefficients.For further consideration we restrict the quadratic term in (10), taking  > 0.
It should also be noted that in [13] the consideration of the corresponding transition is limited by the analysis of ( 6) and ( 9) supplemented by the stochastic components.In this case, the dependence of the shear modulus  on the parameter  is fundamental, when two approximations are using the expansion in powers of  and the inverse exponential function.As a result, near the glass transition with decreasing parameter  the shear modulus drops sharply from the final value  0 to zero (value, which is characteristic for the liquid state).However, for the construction of an adequate analytical scheme according to the Ruelle-Takens theorem [14], a self-consistent behavior of three degrees of freedom should be considered.It is known that a simple threeparameter scheme, describing the self-organizing system is the Lorenz system [15,16].As the study of the selforganization processes, occurring in the ensemble of particles [17] or in a solid [18], shows, it is necessary to take into account the characteristics of each system and the relationships of its basic parameters, except the use of the main approaches for the construction of Lorentz mathematical model.
Based on the above it is natural to assume that the third parameter characterizing the system under study is the field of elastic stresses .As a result the system of synergistic equations, describing the transition to a plastic flow state and taking into account dissipative terms, becomes In contrast to the original relation (9) for the rate of change of the concentration of vacancies  in (11) we have used two additional terms.The first of these is a microscopic dissipation channel with a characteristic relaxation time   and the last one is the hydrodynamic contribution ]∇ 2 k, taking into account the spatial variation of the velocity of the shear displacement (  > 0 is the coupling coefficient, ] =  0 / is a kinematic viscosity).In (12) (comparing with (6)) the additional term is also present.It takes into account the microscopic dissipation due to the mutual action of shear stresses and the number of vacancies on the change in the rate of shear displacement ( V -a positive constant).The velocity of shear displacement plays the role of the field conjugated to the order parameter .The role of the control parameter in the system of synergetic ( 11)-( 13) is played by the stress , for which the microscopic relaxation during time   to the value   , given by the external influences, is taken into account and the relaxation due to the redistribution of vacancies by the spatial variation of the velocity of the shear displacement (  is a positive coupling constant) is considered.
The system of ( 11)-( 13) is very hard to research, because it has a large set of constants: the relaxation times   ,   , material parameters ,  0 ,  0 (] =  0 /), the constant expansion of Landau , voltage   , and coupling constants   ,  = , V, .Therefore, to analyze the system (11)-( 13) it is convenient to use dimensionless variables for the time coordinate (for simplicity restrict dimensional case), the order parameter, the conjugate field, and control parameter to the appropriate scale: where we have introduced the hydrodynamic time scale  =  0 /.As a result of the behavior of solids under external loading with the ensemble of vacancies there is a dimensionless system of equations Here the prime denotes the coordinate, and the additional parameters characterizing the ratio of the constants and parameters of the material are given by

The Adiabatic Approximation
Despite a significant simplification of the system (15) compared with ( 11)-( 13), it is still quite difficult to analyze, so you need to use the approximation.The latter two equations of the system are given by the time ,   , and   .Take into account that the corresponding relaxation time of the vacancy concentration is sufficiently large and is far superior to the hydrodynamic time scale  and the stress relaxation   .As a result of dissipative systems for real picture of the transition from a solid state to a plastic state it is realized when Since the rate of change of the dimensionless ṁ, V , σ is of the same order of magnitude, then (17) allows us to neglect the left side of the last two equations of (15).Moreover, by analyzing the expression (16) in the adiabatic approximation (17) we can put ,  ≪ 1.This leads to the simplified system of equations from which we can obtain the dependence of the conjugate field and control parameter of the order parameter Thus, as shown in Figure 1(a) an increase in the vacancy concentration in the range bounded by the maximum value increases the curvature of the velocity of the shear displacement (dashed lines correspond to nonphysical area  > ).In the case of internal stress  increase in concentration  in the same interval leads to gradual relaxation of the value   , given by external shocks to the value   /2 (see Figure 1(b)).
Further, substituting the first relation (21) in the first equation (18), we obtain the Ginzburg-Landau-Khalatnikov which is determined by the energy of vacancy As shown in Figure 2, with   <  2 energy (23) is a monotonically increasing form (curve 1), with a minimum of  0 = 0, corresponding to solid state.If the voltage   exceeds the critical level  2 (2, 3 curves in Figure 2), you have at least corresponding to the state of plastic flow (see Figure 3).In this case, substituting the value of (24) to (21), we obtain the final values for the curvature of the velocity of the shear displacement V  0 = V  ( 0 ) =  0 and internal stress  0 = ( 0 ) =  2 , corresponding to a given state.

The Stationary Solution
Consider the solution of (22) in the stationary case, when ṁ = 0.The first integral in this case is given by where  0 is integration constant.To determine the latter we consider an ordered phase (state of plastic flow) when the concentration of  =  0 and   = 0, respectively.As a result, constant  0 represents the energy of the ordered state It should be noted that at the boundary of the material ( → ∞), which implemented the disordered phase (solid state), the concentration of  ∞ = 0 and thus the vacancy formation energy  = 0, according to (25), lead to a non-zero gradient where  is characteristic scale, the value of which diverges when the elastic stresses   critical value  2 .The presence of the gradient   ∞ ̸ = 0 does not mean the appearance of the flow at infinity, since the gradient is due to the disappearance of vacancies ( ∞ = 0).
Using the approximation   ≈  2 and given that the rms fluctuations are reduced along the axis , the first derivative we obtain is The final expression for the concentration of vacancies in the definition (27) for the correlation length and the steadystate value (24) takes the form shown in Figure 4 (  is the constant of integration).
In this case the constant   ≥  determines the width of the border area, in which the concentration of vacancies is reduced to steady state (24).Given the dependence of the elastic stress on the concentration of vacancies (21) and expression (29) to find the spatial distribution of the elastic stresses (see Figure 5), we have As follows from Figure 5, in the transition region elastic stresses increase monotonically from the critical value  =  2 to the value   , which is determined by an external influence, which in turn corresponds to increasing the coefficient of friction  = / values from the hydrodynamic  2 0 /( 0   ) to the value   /, corresponding to dry friction.
However, in determining the spatial dependence of the shear displacement difficulties, for example, the first equation (21) we can find the first integral

Conclusion
Our analysis shows that the use of the system ( 11)-( 13) allows us to take into account the effect of point defects (vacancies and interstitials) on the self-consistent picture of the transition from a solid state to a state of plastic flow under the action of external loads.In this case, the order parameter characterizing the state of the material is the concentration of vacancies for crystalline solid or the proportion of free local volume for amorphous solids.
In the framework of the adiabatic approximation, it is shown that an increase in the curvature of the vacancy rate of the shear displacement increases and the internal voltage drops below the level fixed by external factors.Obtained from the Ginzburg-Landau-Khalatnikov vacancy formation energy becomes non-zero minimum  0 , typical of the state of plastic flow, provided that the external voltage   exceeds a critical level  2 , determining the parameters of the medium.In addition, it is shown that the steady-state spatial distribution of the concentration of vacancies is a falling form while shearing stress component increases monotonically.

Figure 3 :
Figure 3: The dependence of the steady-state value of the system parameters  0 .
At the boundary of the material, when  → −∞, the value of the concentration of vacancies When  → +∞, the expression for the vacancy concentration and rate of shear displacement, respectively, take the form  ≈ √   −  2 exp (−2  + which has no analytical solution.But the form of the dependence V() can be analyzed asymptotically.