Translation-rotation plastic flow in polycrystals under creep Текст научной статьи по специальности «Физика»

УДК 539.376

Translation-rotation plastic flow in polycrystals under creep

V.E. Egorushkin1 and V.E. Panin1,2

1 Institute of Strength Physics and Materials Science SB RAS, Tomsk, 634055, Russia 2 National Research Tomsk Polytechnic University, Tomsk, 634050, Russia

Grain boundary sliding is the primary process of plastic flow in polycrystals under creep which is accommodated by multiscale translation-rotation modes of intragranular plastic deformation. The steady stage of creep is characterized by constant creep rates which is related to slow rotation of sliding dislocations and material mesofragmentation. The rate of the tertiary creep is being increased which is related to the similarity of grain boundary sliding and plastic flow within near-boundary zones, where high lattice curvature arises. Here we develop a nonlinear theory of micro, meso- and macroscale plastic flow and propose a constitutive equation for its velocity, rotation modes, local irreversible stress, temperature, and heat and mass flux densities with regard to structural transformation. The theory agrees well with experimental data on multiscale translation-rotation in Al polycrystals under creep, suggesting that fracture in tertiary creep is caused by a multilevel increase in lattice curvature in a deformable polycrystal. Fracture is initiated at the interface of grain boundary sliding and near-boundary regions where lattice curvature results in cracks.

Keywords: nonlinear theory of plastic flow, structural transformation, tertiary creep, multiscale rotations

DOI 10.24411/1683-805X-2018-13003

Трансляционно-ротационное пластическое течение поликристаллов

при ползучести

В.Е. Егорушкин1, В.Е. Панин1,2

1 Институт физики прочности и материаловедения СО РАН, Томск, 634055, Россия 2 Национальный исследовательский Томский политехнический университет, Томск, 634050, Россия

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

Ключевые слова: нелинейная теория пластического течения, структурная трансформация, третья стадия ползучести, многоуровневые ротации

1. Introduction

Creep in polycrystalline materials at constant applied stress is a multiscale process behind which lie translationrotation modes of plastic deformation and fracture. The process is guided by grain boundary sliding as a source of primary plastic rotations which are accommodated, accord-

ing to the angular momentum conservation law, by multi-scale intragranular rotations. Creep provides intriguing patterns of grain boundary sliding because grain boundaries lack translation invariance and triple grain junctions form a stochastic system of internal stress concentrators

[1-9].

© Egorushkin V.E., Panin V.E., 2018

It is commonly supposed that high strain rates early in creep result from stress relaxation at triple junctions whose boundaries generate nonlinear defects as sources of intra-granular dislocation deformation. Next comes steady creep in which point defects in grain boundaries create a steady system of positive ion clusters from which dislocations develop to accommodate intragranular strains. This stage of creep features low rates because of easy rotations at grain boundaries [10] and slow accommodating rotations inside grains, which determines actual rotational modes of poly-crystalline materials under creep. Finally, the creep rate increases steeply up to the point of fracture [3, 4, 9]. The question of what mechanism accelerates the creep at its final, tertiary stage is debatable but it is quite clear that all stages of creep involve multiscale translations-rotations and that their effect at the tertiary stage should be strong to eventually culminate in plastic collapse. The foregoing makes it worthy to elucidate what lies behind rotational deformation in polycrystals under creep.

One of the concepts suggests that plastic rotations in solids are much contributed by lattice curvature [3, 4, 10]. This concept, as applied to strength and plasticity, has been discussed in many papers [11-18] but clarifying the relationship between lattice curvature and rotation modes of plastic flow requires a theoretical description and experimental verification. Here we provide such a description and its experimental verification on the example of plastic flow and fracture in polycrystals under creep, supporting the view that any fracture in a deformed solid represents rotation modes and so does tertiary creep.

2. Theoretical approach and verification method

Our theoretical description of nonlinear plastic flow and its translation-rotation modes is based on the following concept. Unlike elastic equilibrium described by elastic wave equations and Hooke's law (equation of state), plastic equilibrium needs equations descriptive of elastic discontinuities, their density and flux [19, 20]. The constitutive equation of plastic state is determined by positive entropy production via defect and energy flows escaping from deformed regions and responsible for local strain irrevers-ibility. Plastic displacements are a sort of compensating fields for elastic strains in a continuous medium, and their dynamics can be described by gauge fields or curved spacetime geometry. The idea of using gauge fields or Riemann geometry to describe defects in a deformed solid dates back to the 1980s [21-23] but it has not been advanced since then. Attempts to explain macroscopic plastic deformation in terms of microscopic or atomic motion are successful only if we consider model specimens with low defect (dislocation) densities and fail if we consider real materials. The escape of a defect from a strain localization zone is just a driving force for entropy production and irreversibi-lity rather than an event responsible for macroscopic deformation.

Considering structured media on meso- and macroscales in continuum mechanics and electrodynamics, we are to average microscale equations, and in dislocation physics, such averaging leads to dislocation polarization [22, 23]. Electrodynamics considers a medium and a magnetic field as two different objects for which a field-averaged problem is formulated. Plasticity physics and mechanics treat a medium and defects as a single object and needs averaging over its scale. One of the averaging procedures uses scale transformation and e extension [24], making it possible to average microscale equations and to derive their macroscale analogues with extensions descriptive of nonlinear plastic waves in a deformed solid [16, 19, 20]. Such a nonlinear wave is basic in our analysis of deformation processes, allows their description from a single standpoint, and provides the way to better understand the origin of translationrotation plastic flow.

The method of experimental verification consists in analyzing the processes occurring in high-purity (99.999 at %) A999 Al polycrystals under creep at ambient temperature. When creep is investigated at 50°C, this type of polycrystals can experience highly pronounced grain boundary sliding and attendant strong plastic rotations on meso- and macroscales. According to the angular momentum conservation law, the structure of the material is involved in intragranular deformation which is accommodating process by multiscale rotations of opposite sign.

For our verification, we used flat dumbbell-shaped specimens with a gage section of 1.2 x 8 x 37 mm cut from rolled sheet and annealed at 513 K for 0.5 h such that the average grain size was d = 400 ^m; for structural analysis, the specimen surface was electrolytically polished.

The specimens were tested for creep under uniaxial tension at a stress of 11, 13, and 18 MPa and temperature of 323 K (0.35Tm). The elongation was measured with a dial indicator accurate to ±1 ^m. The structure of the material at different creep stages was analyzed by optical microscopy (Axiovert 25CA), interference microscopy (NewView 6200), transmission electron microscopy (JEOL 2100), and scanning electron microscopy (Quanta 2003D).

3. Results

3.1. Plastic strain waves

Let us analyze the dynamics of deformation in an elasto-plastic medium using a macroscale plastic equilibrium equation which relates the flux and density of defects (discontinuities of displacement vector fields), elastic stresses, and potential external sources of plastic strains [14, 15, 20]:

1 cPy

^«5 = 4VM--fpY, (1)

c E

where e^x5 is the Levi-Civita symbol, dx is the derivative with respect to space variables, d0 is the time derivative with summation over reiterating indices, a5 is the linear defect density tensor (for dislocations if the discontinuity

is a multiple of the lattice constant), jR is the plastic flux or defect flux density, aj = E ln uv are local elastic stress concentrators determined by Hooke's law, cRV are elastic moduli, E is Young's modulus, c is the velocity of sound, and pvy is plastic distortion. The plastic distortion pV is of dual nature: in combination with a and J, it determines the current strain, and in combination with aj, it can represent external strain-stimulating sources.

Generally, it is impossible to discard p^ from Eq. (1) in a way similar to that in electrodynamics because gauge conditions are specified by specific problem conditions.

For the elastic sources determined by the scalar x, t) rather than by the displacement vector (volume fraction of deformed material), Eq. (1) takes the form

_ 1 j- y v

£j*8d*a8 - 2 J j + aj - Cjvpy '

(2) c

where a„ = cYdY ln^(x, t), cjj = £ cRV E. For certain load-

R,v

ing conditions, Eq. (1) can also be given in vector form.

Deriving macroscale relations requires scale averaging of Eqs. (1), (2) for which we can use scale transformation and £ extension that gives us the plastic equilibrium equation for macro- and mesoscales [20]:

Vj = 6^*08-^ p;. (3)

Here,

O = -c2Jadt, V = -J, p,v = c2JpYvdt, cjv = £<cRvy>/e,

where (cRV> is the space average derived elsewhere [20] and the time average is over an observation period. The opposite signs of O and V indicate that these are mass parameters rather than defect characteristics.

The averaged equations provide the following valuable data. First, meso- and macroscale plastic deformation, like microscale deformation [25], develops for no change of volume because such changes vanish during scale averaging (volume conservation law). Second, we have an elastic stress concentrator after averaging, which also provides similarity between macro- and microscale equations, in particular (1), except only for variations in averaged physical parameters (elastic moduli, density, etc.). Third, strain-induced defects are responsible for instantaneous plastic flow velocities. Thus, scale similarity holds for plastic deformation.

The above conclusions follow from nonlinear equilibrium and from any local solutions of its equation. Because plastic deformation is actually the evolution of local displacement fields and their discontinuities, we need solutions of (2), (3) for a local region rather than for the whole coordinate space to trace the behavior of a plastic strain wave as a certain space curve in this region. The theory of space curves is based on the Frenet equation in a local coordinate system, allowing us to obtain nonlinear equations of motion solvable by the inverse scattering method [20].

The equations for localized plastic strain (2), (3) show that deformation develops as nonlinear spiral waves the characteristics of which are determined by the geometry of a deformed region (bending-torsion), its size (length, lateral dimension), and physical parameters (elastic modulus, density) [20]. Reasoning that plastic strain waves are a fundamental phenomenon of macro-, meso-, and microscales, let us consider these equations and their solutions at greater length as they allow us to better understand how such waves contribute to plastic and other types of deformation, to analyze differences in plastic flow on different scales (up to nano), and to link the physics of defects and the mechanics of plastic deformation.

Following the method proposed elsewhere [20], we can represent the term cjv p^ as the sum of the rotor of a certain vector O and the scalar gradient f:

cjv pV=£j^8d*°'8+dj f (4)

and take the rotor from both parts of Eq. (3). Then, (3) assumes the form

£nkj£j*8d kd x°8 = Wn, (5)

where Wn = £nkjdkVj is the plastic flow vorticity and O = O + O' contains the rotation component of external plastic distortion. Solving (5) in a local coordinate system (normal n, binormal b, tangent t) together with the Frenet equation gives the plastic wave velocity proportional to the curvature k(S , t) of a deformed region and directed along the binormals b(s, t):

V = — x(s, t) b(s, t) A(L, d), 4n

(6)

where L is the length of a deformed region, d is its lateral dimension, s is the current length, and b is the Burgers vector. The scale factor is A = ln(2L/d) -1 for L >> d and A = 0.2 for L ~ d.

Passing from the motion equation for plastic flow to a nonlinear equation allows us to calculate the curvature k(S, t) = 4R sech(2R)(S - at), where a and R are the torsion and maximum curvature of a deformed region: R determines the plastic flux and a gives the rate of its change along S.

If we select a global coordinate system such that the z axis is directed along L, and X and E are varied in the range of lateral dimensions, we have the following expressions for the strain and rotation tensors, £jV and rojV, in the Cartesian coordinate system:

£ „ = 1 —

1 + v2

-sech2£,

(7.1)

£ z* (®z*) =

1

abxt

(

1+ V2 U(x2 + y2)

1

1

cosh2£ cosh2n ( sin(v£) sinh £ v cos(v£) ^

cosh2£

cosh £

(7.2)

Fig. 1. Plastic wave form and velocity (arrows) in a deformed region

ezy (ro^ ) =

abyt

1 + v2 U(x2 + y2) cos(v£) sinh £

cosh2£ v sin(v£)

cosh2n

(7.3)

cosh2£ cosh £

The subscripts in brackets in the left-hand side of (7) correspond to the rotations , £ = 2Pz + 2a(3b/nAt, n = 2Pz, and v = -a/p.

The quantity abA/ n represents the plastic wave velocity in the global coordinate system. In the local coordinate system, the wave velocity is determined solely by the torsion a.

It should be noted that the scale factor A was initially used to determine the time t in a common coordinate system. This means that the time, being a parameter of plastic deformation, has no absolute meaning and the wave velocity is controlled only by torsion.

Figure 1 shows the form of a plastic strain wave in a deformed region according to (6), (7), with arrows indicating its velocities polarized along the binormals b perpendicular to the propagation direction.

3.2. Constitutive equation

Any deformed region involved in plastic strain localization generates defect and energy flows with attendant heat and mass transfer. In this case, the work function determines the specific entropy production S (dissipative function) and the entropy flux I:

S = (Q - Jy)d T - JdjV dTvk + a^vk

I

T T

= roikvk -Q- J

T T T '

(8)

where roik = eijkaj on microscale, a—strain-induced defect, dkUi - dUk on macroscale or antisymmetric component of plastic distortion (rotation mode), v is the plastic strain rate, a is the local stress, T is the temperature at local equilibrium, eijk is the Levi-Civita symbol, Qf = -k9T/p is the specific heat flux, k is the heat conductivity, p is the material density, dT are the temperature gradient components, J = -Ddic(z)/p is the specific mass flux, D is the diffusion coefficient, dic are the concentration gradient

components, and y = - y2 is the difference between the chemical potentials of material and diffusing component.

The first and second terms in the right-hand side of (8) for S represent the entropy production due to heat release and mass transfer; the other two, due to plastic deformation [20]. The former two terms are of second-order infinitesimal with respect to the gradients 5fT5fT, 5fT5iC, 5iC5iy, and the thermoplastic term is of first-order infinitesimal.

The second law of thermodynamics gives us a constitutive equation being the equation of plastic state for a deformed system:

ak vk - - J*-^ <% vk * 0. (9)

It is now evident that the plastic flow supporting this state becomes irreversible. For the second-order terms being small with respect to the gradients, relation (9) is given by

\IVk > 0. (10)

In the elastic case, the strain and the stress are related by Hooke's law. Likewise, in the plastic case, expressions (9) and (10) relate the macroscale flux, density (translation and rotation modes), and local stress. The deviator ak in (9), (10) is the irreversible stress equal to the difference between the total local stress and its value at the yield point (elastic limit).

Inequalities (9), (10) represent plasticity with work hardening; for perfect plasticity at vk ^ 0, an equality is possible. With no temperature gradient, Eq. (10) coincides with the well-known equation of state in plasticity mechanics [26]. Let us dwell on some particular aspects of irreversibility.

Relations (9), (10) suggest that the increment in the work function is positive (Drucker postulate). For unloading (a < 0), equality (10) holds only at vk = 0, implying that unloading proceeds via elastic deformation. Inequality (10) for unloading holds only at vk < 0: reverse plastic deformation occurs at total stresses lower than the yield stress (elastic limit), i.e., the process involves softening (Bauschinger effect).

The perfect thermoplastic state with no work hardening, i.e., with (10) equal to zero, corresponds to steady conditions (9S/dvk = 0) and zero total thermodynamic strength:

ak-T ®ik = 0. (11)

This implies that every event of plastic deformation provides conditions for stress relaxation via rotations at gradient temperature.

Creep at ambient temperature results in steady plastic flow when absorbed heat and mass transfer (migration and diffusion) transform into mechanical energy for no change in stress. The constitutive equation of state (11) for these

0.04 "" 006 Distance, mm

Fig. 2. Denuded grain boundary MN of grains A and B in A999 aluminum under creep at T = 323 K, o = 11 MPa: interference image (a), surface profile normal to MN (NewView 6200) (b)

conditions represents plastic rotation modes at different creep stages. Let us show it experimentally.

3.3. Experimental verification of theoretical predictions

When analyzing three-stage creep in polycrystals, one should allow for a wide range of potential stress sources Vf. The high rate of creep at its first stage is due to stochastically distributed stress concentrators at triple junctions. At the stage of prefracture, constrained rotations induced by grain boundary sliding create highly nonequilibrium conditions and increase the stress in locally thinned zones (analogue of quasi-necking in cylindrical specimens under uniaxial tension). These potential stress sources are ignored in most creep theories, and the expressions derived above are correct only for the second or steady creep stage.

From relation (11) it follows that the macroscale rotation mik in steady creep must increase with increasing applied stresses ok, which is confirmed by experiments.

Figures 2, 3 show grain boundary sliding zones in high-purity A999 Al polycrystals under creep at oJ =11 MPa (Fig. 2) and o2 = 18 MPa (Fig. 3). It is seen that although the denuded grain boundary at oJ =11 MPa measures 12 ^m, the accommodating rotation mode of near-boundary deformation is very local and homogeneous. This is quite clear because such conditions provide intragranular

rotations via microscale formation and misorientation of a cellular dislocation substructure [3, 4]. At o2 = 18 MPa, the denuded grain boundary measures a mere few micrometers but the material involved in plastic rotations near the boundary is extruded through mesoscale displacements of coarse bands. It is this scaling effect of rotation modes that Eq. (11) predicts for creep at increasing applied stresses.

The most severe rotations in A999 Al polycrystals develop at the third creep stage (Fig. 4). Early in tertiary creep, grain boundary sliding causes multiscale fragmentation in near-boundary grain zones. As can be seen in Fig. 4, a, grain C is involved in sliding with a highly pronounced tangential component along its denuded boundary (shown by an arrow), which results in complex near-boundary fragmentation. Outer near-boundary layer AB is severely fragmented such that the initial structure degrades. Inner layers 1-3 are displaced relative to each other with the formation of denuded subboundaries ab, cd, ef and cross steps. Thus, the material is fragmented into blocks preserving their translation-invariant structure. It is this invariance that provides dS/dvk = 0 in a fragmented material. Noteworthy is that while grain C as a whole rotates clockwise due upward displacements of its zone B (Fig. 4, b), layers 1-3 are displaced anticlockwise with the formation of blocks. In other words, grain boundary sliding and fragmentation show self-consistency obeying the angular momentum conservation

1.5

S 0.0'

Cg

ie

ffi -1.5-

-3.0-

A B

0.00

0.02

0.04 0.06 Distance, mm

0.08

0.10

Fig. 3. Extrusion bands along boundary MN of grains A and B in A999 aluminum under creep at T = 323 K, o = 18 MPa: interference image (a), surface profile normal to MN (NewView 6200) (b)

-i-1-1-1-1-1-1-r

0.00 0.05

Distance, mm

Fig. 4. Extruded mesofragments along boundary AB of grains C and D in A999 aluminum under creep at T = 323 K, a = = 18 MPa: interference image (a), profile of fragmented meso-band 1 (NewView 6200) (b)

law. Cleary, such structural degradation in tertiary creep cannot be described by relations (8)—(11).

Late in the tertiary creep at a2 = 18 MPa, the surface structure of A999 aluminum reveals a highly nonequilibrium state and multiscale local cross-sectional thinning (quasi-necking) indicative of prefracture. It is seen from Fig. 5 that the surface layer in the prefracture zone is involved in rotations via thin mesobands of localized plastic flow. At higher magnification (Fig. 5, b), one can identify that such mesobands propagate stepwise with step heights of about 2 ^m. Mesoband AB, which causes quasi-necking before fracture, consists of thin transverse shear bands with steps of ~0.5 ^m. Thus, the material in the quasi-neck zone undergoes submicron fragmentation. Such multiscale fragmentation induced by transverse sliding within mesobands can provide any rotations in a material involved in tertiary creep.

The potential stress sources responsible for quasi-necking before fracture greatly increase the rate of tertiary creep. Because macromechanics is yet unable to account for stochastic potential stress sources at the stages of primary and

2.00

i 1.25

ie 0.50-

0.25

1.00

b

■y . B

0.00

0.01 0.02

Distance, mm

0.03

0.04

Fig. 5. Transverse bands inside strain localization mesoband AB (quasi-neck region) in A999 aluminum under creep near the point of fracture at T = 323 K, a = 18 MPa: interference image (a); profile of fragmented mesoband AB (NewView 6200) (b)

tertiary creep, Eqs. (8)-(11) refer to nonlinear mesomecha-nics of multiscale steady creep.

As has been shown [3, 4], polycrystals in tertiary creep undergo delamination in near-boundary intragranular sliding zones (Fig. 6), and this is associated with high lattice curvature and with bifurcations arising at interstices. It is logical to expect that the latter two factors will increase the role of rotation modes in near-boundary intragranular deformation. Let us analyze such rotations in the context of concentration fluctuations [27, 28].

3.4. Concentration fluctuations in translation-rotation plastic flow

Let us consider concentration fluctuations 8c in a deformed region at plastic equilibrium when O ^ 0 (see Eq. (3)) and write the free energy for them [27]:

Fig. 6. Noncrystallographic sliding within near-boundary zone of fracture surface as a consequence of grain boundary sliding along boundary AB between grain C and D; tertiary creep of Al polycrystal; scanning electron microscopy

F = j{-a| 8c|2 + ||8c|4 +

+ envk<^vrk |8c

+ 02 \ dv.

(12)

Let 8c be maximal only along the z axis, i.e., dc/dx, dc/dy = dc/dz and c = c(z). In the general case, reasoning that 8c represents short-wave fluctuations, we have the coefficient a = a(T) + a'(q - q0)2 [27] for \q - q01 << q0. The vector q 0 in the reciprocal space corresponds to the vector l in the straight space, which limits the scale of plastic mis-orientation O such that q01 = 2n. The q-dependent term in the expression for a is responsible for segregation (weak crystallization) in the region of concentration fluctuations. Because the effect of segregation on superplasticity is slight, we can omit a'| q - q01 and assume the coefficient a to depend only on temperature a(T) = a(T - T ), where T is the temperature of stability loss for a high-temperature phase.

Minimization of (12) with respect to 8c gives d28c d8c

a2-

- + Sc -|Sc|2 8c = 0,

(13)

"z dt

where 38c/dt = v d8c/dz, and v = 2uJ(al2) represents the velocity of fluctuation motion along the z axis and depends on plastic displacements u^ = e^vk Ov , temperature a(T), and misorientation scale.

For deriving a relation for concentration fluctuations and plastic flow, we minimize (12) with respect to uu for which the parameter 8C = |8C|ce(r) is expressed as the product of its absolute value and argument such that the total velocity of superplastic flow is determined by the phase of nonuniform fluctuations.

The resulting equation together with (3) gives the system

X

du0 = ^Vu +

8

|2 ^vfcdv<£k,

J

e^vkk = Vn + c¿vPp ,

(14)

(15)

where X = b/a, £ = al2/l, |8c| * 0.

The above two equations describe the interaction and energy exchange between long-wave plastic fluctuations and short-wave concentration fluctuations. The incompatibility of Eqs. (14) and (15) provides low rates of steady creep. High-rate rotations in grain boundary sliding cannot be immediately accommodated by microscale rotations of intragranular dislocation glide, as is required by the angular momentum conservation law. Because dislocations are generated by structural transformation flows in the planar subsystem of a solid (surface layers and all internal interfaces, including grain boundaries) [14, 15], there should be a tangential displacement component in grain boundary sliding and a sinusoidal normal stress distribution at grain boundaries [29]. In the grain boundary zones of tensile

normal stress, clusters of positive ions arise and generate dislocation cores in the zones of intragranular sliding. The rotation mode of dislocation deformation develops slowly, and this greatly decreases the grain boundary sliding velocity.

Equations (14) and (15) show compatibility for tertiary creep when high lattice curvature and high concentrations of interstitial states arise in near-boundary grain zones. In this case, the concentration fluctuations described by (14) occur both at grain boundaries and in near-boundary intra-granular sliding zones. The velocities vu and V take very high values of the same order. The rotation modes of grain boundary sliding are fully accommodated by the concentration fluctuations in near-boundary intragranular sliding zones where the material structure is delaminated and degenerated. Such an effect was observed earlier [3, 4] but had no explanation at that time (Fig. 6).

The energy exchange between grain boundary sliding and concentration fluctuations in near-boundary plastic flow is most efficient under resonance conditions, with vu and V having the same order. In this case, collective motion provides displacements of an area equal in size to a concentration kink, i.e., Eq. (13) has the solution

8c: ±th

(

1 z - vt

(16)

where v and t are determined after solving (13).

Combination of (3) and (14) gives the velocity of plastic flow for tertiary creep:

X

X

cP Pv

c|av pp .

IS 12 ^ - ^ 15- 12 H (17)

|8c| |8c|

The first term in the right-hand side of (17) represents the velocity of plastic wave flow (6) increased by a factor of X/|8c|2.

The quantity Vu was calculated with regard to the rotation component of plastic distortion O' and plastic displacement uu. Thus, the term £vu in the right-hand side of (17) is determined by the gradient of external sources, and the next term by delayed polarization when Vu reaches its maximum value. At cUv ^ 0 in a deformed region, the last term of (17) is responsible only for near-boundary intra-granular deformation. The rate of tertiary creep calculated by (17) is in good agreement with experimental results.

4. Discussion

The necessity of accounting for lattice curvature in a deformed solid is pointed out in many papers [11-18]. For example, high lattice curvature is responsible for shear banding and plastic distortion in surface-hydrogenated Ti and commercial Al with a thin A999 Al film under alternate cyclic bending [17] and for dynamic rotations in metal surface layers during shock wave propagation [10] and opening mode cracking [14, 15]. These effects of lattice curvature suggest that its zones develop translation-rotation de-

Fig. 7. TEM image showing extinction contours in near-boundary zones with high lattice curvature and structural degradation in polycrystalline aluminum; tertiary creep at T = 50°C

Fig. 8. TEM image of subboundary KLM formed during submicron fragmentation in polycrystalline aluminum; tertiary creep at T = 50°C

fects, which are not considered by conventional theories in which the lattice is taken as translation-invariant.

The lattice curvature x induces a counter stress field azz and is expressed as [30]:

X = 4|3sech[2P(s + 4vt)], P = (18)

dE

where d is the width of a curved plane, s is the current coordinate, t is time, and v is the velocity of local structural transformation in the curved plane in a local coordinate system. From expression (18) it follows that translationrotation defects should arise in a deformed solid, providing relaxation of the stress field azz and transformation of the elastic curvature into plastic rotation.

When the pair interaction potential has a single minimum, the particle system potential can take bifurcation minima due to increasing interatomic spacing [12, 13] such that atoms can pass from lattice sites to interstices, forming "translation-rotation" strain-induced defects and providing noncrystallographic plastic flow. Such strain-induced defects arise in all deformed solids but their identification is easiest in polycrystals under creep which reveal highly pronounced accommodation of their primary grain boundary sliding by translation-rotation plastic flow in near-boundary intragranular regions. Let us consider some experimental evidences for the generation of rotation defects in tertiary creep (Figs. 7, 8).

Figure 7 shows TEM image with extinction contours demonstrating the presence of high lattice curvature which is equal to 25° |xm-1 in a near-boundary region. The sliding velocities V and v^, in this case, are of the same order, and the result is an abrupt increase in the rate of tertiary creep.

Very important data on the formation mechanism of subboundaries in grains at the stage of tertiary creep is represented in Fig. 8. Early in the formation, the structure of subboundaries points to their shear band nature having a

lot of the short local noncrystallographic shears within the high lattice curvature region (Fig. 8, fragment LM). At the second stage, the shear band makes a sharp turn and form a fiber structure (Fig. 8, fragment KL) which greatly contributes to low-angle misorientations of subgrains in tertiary creep. This type of strain-induced defects should be introduced to describe nonlinear plastic deformation and fracture in tertiary creep [10, 17].

5. Conclusions

Thus, the developed nonlinear theory of micro-, meso-and macroscale plastic flow gives a constitutive equation for the velocity of plastic flow, rotation modes, local irreversible stress, temperature, heat and mass flux densities with regard to structural transformations. The theory agrees well with experimental data on multiscale translation-rotation modes in Al polycrystals under creep. The results suggest that such modes of intragranular plastic flow are initiated by grain boundary sliding as a primary process at all creep stages. Slow intragranular rotations provide constant rates of steady creep, and very high rates of tertiary creep are provided by high lattice curvature and numerous non-equilibrium interstices that arise in near-boundary regions due to grain boundary sliding. The structure of a material thus transforms through noncrystallographic sliding in local high curvature zones up to the point of plastic collapse and fracture.

Acknowledgments

The work was performed under Fundamental Research Program of the State Academies of Sciences of Russia for 2013-2020 and supported by the International Research Staff Exchange Scheme under 7th EU Framework Program (grant No. 612552), Russian Foundation for Basic Research

(project No. N170100691), and President of the Russian

Federation (grant for support of leading scientific schools No. NSh-10186.2016.1).

References

1. Kassner M.E. Fundamentals of Creep in Metals and Alloys. - Oxford: Elsevier Science, 2009.

2. Gifkins R.C. Grain-boundary sliding and its accommodation during creep and superplasticity // Met. Trans. A. - 1976. - V. 7. - No. 8. -P. 1225-1232.

3. Panin V.E., Elsukova T.F., Surikova N.S., Popkova Yu.F., BorisyukD.V. The role of rotation modes of plastic deformation in failure of high-purity Al at ambient temperature creep // Deform. Razrush. Mater. -2016. - No. 12. - P. 2-9.

4. Panin V.E., Egorushkin V.E., Panin A.V., Chernyavskii A.G. Plastic distortion as a fundamental mechanism in the nonlinear mesomecha-nics of plastic deformation and fracture // Phys. Mesomech. - 2016. -V. 19. - No. 3. - P. 255-268.

5. Wadsworth J., Ruano O.A., Sherby O.D. Denuded zones, diffusional creep and grain boundary sliding // Metall. Mater. Trans. A. - 2002. -V. 33. - P. 219-229.

6. Matsunaga T., Kameyama T., Sato E. Grain boundary sliding induced by lattice dislocation activity during ambient temperature creep in hcp metals // IOP Conf. Mater. Sci. Eng. - 2009. - V 3. - Article 012014.

7. Matsunaga T., Kameyama T., Ueda S., Sato E. Grain boundary sliding during ambient-temperature creep in hexagonal close-packed metals // Phil. Mag. - 2010. - V. 90(30). - P. 4041-4054.

8. Ruano O.A., Sherby O.D., Wadsworth J., Wolfenstine J. Diffusional creep and diffusion-controlled dislocation creep and their relation to denuded zones in Mg-ZrH2 materials // Scripta Mater. - 1998. -V. 38. - P. 1307-1314.

9. Wilshire B., BurtH. Tertiary creep of metals and alloys // Int. J. Mater. Res. Adv. Tech. - 2005. - V. 96. - No. 6. - P. 552-557.

10. Panin V.E., Egorushkin V.E., Elsukova T.F., Surikova N.S., Pochiva-lov Yu.I., Panin A.V. Multiscale Translation-Rotation Plastic Flow in Polycrystals. Handbook of Mechanics of Materials / Ed. by C.-H. Hsuch et al. - Singapore: Springer Nature, 2018. - doi 10.1007/978-981-10-6855-3_77-1.

11. Myasnikov V.P., Guzev M.A. A geometrical model of the defect structure of an elastoplastic continuous medium // J. Appl. Mech. Tech. Phys. - 1999. - V. 40. - No. 2. - P. 331-340.

12. Guzev M.A. Structure of displacement and stress fields in a non-Euclidean continuum model // Solid Mechanics: Nizhny Novgorod State Univ. Vestnik. - 2011. - V. 4. - P. 1461-1462.

13. Guzev M.A., Dmitriev A.A. Bifurcational behaviour of potential energy in a particle system // Phys. Mesomech. - 2013. - V. 16. - No. 4.-P. 287-293.

14. Moiseenko D.D., Panin S. V., Maksimov P. V., Panin V.E., Babich D.S., Berto F. Computer simulation of material behaviour at the notch tip: Effect of microrotations on elastic energy release // AIP Conf. Proc. -2016. - V. 1783. - P. 020157. - doi 10.1063/1.4966450.

15. Panin V.E., Grinyaev Yu.V, Egorushkin V.E. Foundations of physical mesomechanics of structurally inhomogeneous media // Mech. Solids. - 2010. - V. 45. - No. 4. - P. 501-518.

16. Panin V.E., Egorushkin V.E., Panin A.V. Nonlinear wave processes in a deformable solids as a multiscale hierarchically organized system // Phys. Usp. - 2012. - V. 55(12). - P. 1260-1267.

17. Panin V.E., Egorushkin V.E., Surikova N.S., Pochivalov Yu.I. Shear bands as translation-rotation mode of plastic deformation in solids under alternate bending // Mater. Sci. Eng. A. - 2017. - V. 703. -P. 451-460.

18. Tyumentsev A.N., Ditenberg I.A., Korotaev A.D., Denisov K.I. Lattice curvature evolution in metal materials on meso- and nanostructural scale of plastic deformation // Phys. Mesomech. - 2013. - V. 16. -No. 4. - P. 319-334.

19. Egorushkin V.E. Gauge dynamic theory of defects in nonuniformly deformed media with structure. Interface behavior // Russ. Phys. J. -1990. - V. 33. - P. 135-149.

20. Egorushkin V.E. Dynamics of plastic deformation: Waves of localized plastic deformation in solids // Russ. Phys. J. - 1992. - V. 35. -P. 316-334.

21. Kroner E. Gauge Field Theories of Defects in Solids. - Stuttgart: Max. Planck. Inst., 1982.

22. Kadic A., Edelen D.G.B. Gauge Theory of Dislocations and Disclina-tions: Lecture Notes. - Heidelberg: Springer-Verlag, 1983.

23. Edelen D.G.B., Lagoudas D.C. Gauge Theory and Defects in Solids. -Amsterdam: North Holland Publ., 1988.

24. Sanchez-Palencia E. Non-Homogeneous Media and Vibration Theory: Lecture Notes in Physics. - Berlin: Springer, 1982.

25. Kosevich A.V. Dynamical theory of dislocations // Sov. Phys. Usp. -1965. - V. 7. - P. 579-609.

26. Ziegler H. Some Externum Principles in Irreversible Thermodynamics with Application to Continuum Mechanics // Progress in Solid Mechanics. - Amsterdam: North Holland, 1963.

27.Lebedev V.V. Fluctuation Effects in Macrophysics. Moscow: MTsNMO, 2004.

28. Patashinsky A.Z., Pokrovsky V.L. Fluctuation Theory of Phase Transitions. - Moscow: Nauka, 1982.

29. Cherepanov G.P. On the theory of thermal stresses in thin bounding layer // J. Appl. Phys. - 1995. - V. 78. - P. 6826-6832.

30. Landau L.D., Lifshits E.M. Theory of Elasticity. - Moscow: Nauka, 1987.

Поступила в редакцию 23.03.2018 г.

Сведения об авторах

Егорушкин Валерий Ефимович, д.ф.-м.н., проф., вне ИФПМ СО РАН, root@ispms.tomsk.ru

Панин Виктор Евгеньевич, д.ф.-м.н., проф., акад. РАН, зав. лаб. ИФПМ СО РАН, проф. ТПУ, paninve@ispms.tsc.ru