ON A DIFFERENTIAL CONSTRAINT IN THE CONTINUUM THEORY OF GROWING SOLIDS Текст научной статьи по специальности «Математика»

Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki

[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2019, vol. 23, no. 4, pp. 646-656 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1696

Mechanics of Solids

MSC: 74A20, 74A60

On a differential constraint in the continuum theory of growing solids

E. V. Murashkin, Yu. N. Radayev

A. Ishlinsky Institite for Problems in Mechanics, Russian Academy of Sciences, 101, pr. Vernadskogo, Moscow, 119526, Russian Federation.

Abstract

The present paper is devoted to the problem of boundary conditions formulation for asymmetric problems in the mechanics of growing solids (MGS). The boundary conditions on the propagating growing surface (PGS) is the fundamental problem of this branch of mechanics. Results from the algebra of rational invariants are used for deriving constitutive equations on PGS. Geometrically and mechanically consistent differential constraints are obtained on PGS. Those are valid for a wide range of materials and metamaterials. A number of constitutive equations on PGS of different complexity levels are proposed. The boundary conditions simultaneously can be treated as differential constraints within the frameworks of variational formulations. The differential constraints imply an experimental identification of constitutive functions. For this reason, the obtained results furnish a general ground in applied problems of the MGS.

Keywords: 3D printing, surface growth, stress, constitutive equation, rational invariant, differential constraint, complete system.

Received: 30th April, 2019 / Revised: 12th August, 2019 / Accepted: 16th September, 2019 / First online: 18th November, 2019

Research Article

3 ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:

Murashkin E. V. , Radayev Yu. N. On a differential constraint in the continuum theory of growing solids, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2019, vol. 23, no. 4, pp. 646-656. doi: 10.14498/vsgtu1696.

Authors' Details:

Evgenii V. Murashkin https://orcid.org/0000-0002-3267-4742

Cand. Phys. & Math. Sci., PhD, MD; Senior Researcher; Lab. of Modeling in Solid Mechanics; e-mail: [email protected]

Yuri N. Radayev © https://orcid.org/0000-0002-0866-2151

D.Sc. (Phys. & Math. Sci.), Ph.D., M.Sc., Professor; Leading Researcher; Lab. of Modeling in Solid Mechanics; e-mail: [email protected], [email protected]

1. Introduction

Traditional methods of manufacturing products of complex shape imply a variety of technological processes. Such processes are coupled with the removal of material and based on the synthesis of products by successively adding material to a surface of complex shape. Manufacturing products by adding new material is widely used in modern engineering [1]. Such additive manufacturing methods include: laser stereolithography [2], selective laser sintering [3], electron beam melting [4], fusion simulation, multi-jet modeling method [5], manufacturing by lamination, 3D printing [6-8], computer axial lithography [9], layer-by-layer concreting [10].

3D printing (3DP) is similar to selective laser sintering technology, but melting isn't used here. An object is formed from powder material by gluing, using inkjet printing to apply liquid glue. 3D printing technology allows color modeling by adding dyes to the adhesive (directly during printing), or by using multiple printheads with color glue [6-8].

Layer-by-layer concreting of constructions also belongs to the methods of additive technologies [10]. Layer-by-layer concreting can be divided into two types, depending on the time between layer fills. In the case of a "hot seam," the break between the layers is less than 12 hours. The second type of layer-by-layer filling is a "cold seam." In this case, it is necessary to wait for the complete hardening of the previous layer to exclude the possible cracking of the non-hardened part under the action of the newly added material.

The additive manufacturing techniques described above basically use well-known ideas from natural processes: accretion of space objects, the formation of avalanches and glaciers, crystal growth processes, the growth of atherosclerotic plaques [11,12]. All these phenomena are characterized by presence of PGS. The growth of atherosclerotic plaque can be described as the process of initial infiltration of blood plasma components into a thin surface layer of the inner wall of the artery. The growth of a crystal nucleus occurs by adding individual atoms or their groups to crystal surface.

A solution of applied problem of growing solids mechanics is sometimes a sophisticated and time-consuming procedure [13-19]. An substantial feature of the formulation of boundary value problems of the MGS is the formulation of boundary conditions on the interface between the source material and the added part [20-22]. In this paper, several variants of the constitutive equations on the growing surface are discussed, starting from the simplest relations (see the book by G.I. Bykovtsev: [22, Pp. 288-292]) to some significant generalizations of the theory. Throughout the paper the terminology and notations adopted in the publications [16-19,22,24] are used.

2. Rate equations of mechanics of growing solids (MGS)

The governing equations in MGS are often conveniently furnished as the rate equations. In this case, the equilibrium equation formulated in an arbitrary spatial system of curvilinear coordinates xk (k = 1, 2, 3) for the asymmetric stress tensor can be presented as partial differential equation

V (d.aji) + d.X1 = 0,

where is the Cauchy force stress tensor, X% is the volume forces, Vj is the covariant differential operator with respect to a given spatial metrics, d. is the time derivative take with fixed spatial coordinates xk.

In order to close the system of differential equations, it is required to specify constitutive equations according to the constitutive law which is an intrinsic characteristics of a particular solid/process. A model of growing solid should also be supplemented by boundary conditions on PGS.

Growing surface £ propagating in space is determined by an implicit form equation:

t = t (xi).

In this case, the unit normal vector n on the surface £ directed towards its propagation is given by the spatial gradient as follows

nl = cdiT, c = |Vr |-1, (1)

where c is the linear velocity of PGS in the normal direction n, which is defined according to

1^1 c = lim —-—.

St^Q St

Here iPP'| is the length of normal vector directed from an arbitrary point P on PGS at time t to point P' which is intersection of normal vector n and PGS at time t + St.

To find out the force stress tensor by the given rates of stresses, one can employ the formulae given in [22]:

aij = f [d.aij(xk, t')]dt' + Iji + aij(xk), (2)

jt+q C *+Q

Iij = [d.aij (xk ,t' )]dt',

Jt-Q

where — stress jump related integral, a^(xk) = a^(xk,t)\t_*^xS)-Q are the

force stress tensor components at the time t = t(xs) — 0 right before an element

starting formation of the main solid. The time t = r(xs) + 0 corresponds to moment right after an element has been attached to PGS. In order to simplify script of formulas we retain t = *(xs) instead of a more correct t = * (xs) +0. After substituting of actual stresses (2) into the equilibrium equation

V aji + Xi = 0 (3)

one can readily come to

V,{jf [d.aji(xk,t')]dt' + Iij + aji(xkXi = 0 (t ^ * + 0). (4)

After transformations according to the rule of differentiation of the integral depending on a parameter, the equation (4) is turned to

J* Vj[d.aji(xk, t')]dt' + VjIji - (Vj*)[d^aji(xk, i)j

* + t=r(xs)

+ Vj aji(xk )+ Xi = 0. (5)

Upon substituting the Eq. (1) in Eq. (5)

J Vj[d.aji(xk, t')]dt' + VjIji - c-1nj [d^aji(xk, i)]

t=r(xs )

+

+ Vj aJt(xk)+ X* = 0. (6)

The following equation can be obtained by integrating (6) and taking account of (3), then changing the order of derivatives

-X1 + X '

t=r(xs)

+ Vj Iji - c-1nj [rlaji(xk ,i)] * + Vj aji(xk )+ Xi = 0. (7

t=T(xs)

After rearrangement of terms Eq. (7) reads

Vj aji(xk ) + Vj Iji +

t=T(xs)

- nj [d.aji(xk,i)]

t=T(xs)

Hereinafter the notation ( ■ )l », „. will be dropped out for values calculated on

v y H=t(Xs)

PGS.

Finally, the constitutive boundary condition on PGS can be formulated in form

c[Vjaji(xk) + VjIji + Xi] - nj[d.aji(xk, i)j = 0 (t ^ * + 0). (9)

The Eq. (9) should be considered as a differential constraint for stresses on

*

PGS. If it is possible to express the stresses a in terms of the actual stresses a\t=T(Xs), then from the equation (8) we can derive the differential equation for

stresses. Alternative approaches to the boundary conditions derivation are discussed in detail, for example, in [20,21].

*

In the general case, the stresses a are to be expressed in terms of the actual stresses on PGS by a tensor constitutive equation

a = n,...).

(10)

Omitted arguments of the function F in the Eq. (10) stand as additional parameters characterizing the growth process. Those can be related to multiphysics phenomena. In a simple model, list of additional parameters may be empty. In particular, the function F may depend on the microstructural directors and the thermophysical hidden variables associated with PGS. The physical sense of the

1The spatial coordinates xk are the parameters, and the lower limit of integration depends

0

c

additional directors [23-26] may be associated with the characteristic directions of the fibers laying in woven composite materials, reinforcement in concrete structures, winding of threads in a bobbin, etc. The function F, in fact, should depend on combinations of arguments which are invariant under coordinate-frame rotations around of the director n.

3. A simple formulation of differential constraint on PGS

Consider in further details a geometric interpretation of the case discussed earlier in [22]. Let traction vector t be known or prescribed on PGS:

t = n ■ a.

* *

In this case, we introduce the following notation a for the reduced tensor a

2d

to the two-dimensional plane element T and assume a simple form of the Eq. (10) as follows

a = F(t, n).

2d 2d

*

Tensors a and F are represented in the Cartesian coordinate frame by 2x2 2d 2d

matrices

* i &11 0"12 A F = ( Fll F12

2d \ ^21 &22 J 2d \ F21 F22

Hereafter, it is convenient to use the decomposition of the traction vector t into the parts

t = t± + t||.

Here t^ is the projection of the traction vector on the tangent plane T to the instantaneous PGS, t| denotes the projection on the normal direction.

As for arguments of the tensor function F, we choose the joint rational invariants of the second-rank tensor a and the unit vector n, which are unchanged under coordinate frame rotations around the director n. In this case, we write the rationally independent system of invariants [27] in the following form

|t|2 = t ■ t, |t||| = |t ■ n|, |t±|2. (11)

There is an obvious rational syzygy in the system of invariants (11)

|t|2 = |t|||2 + |t±|2.

After eliminating the first invariant from (11) the system of independent rational invariants takes the following form

|t|||, |t±|2. (12)

The constitutive equation on PGS, taking account of (12), can be furnished in the form

a = F(|t|||, |t±|2). (13)

2d 2d "

The equation (13) has a clear mechanical sense. We choose a Cartesian coordinate system so that the unit vector k is directed along the normal n to PGS,

and for the projection of the traction vector t in the tangent plane T to PGS, the following equation is valid

|t±|2 = |t±- i|2 + |t±-j|2. (14)

Then the projections of the vector t onto the unit vectors of the specified coordinate system can be expressed in terms of actual stresses by

|ty| = k^L |t±'*!2 = ^231), |t±-3\2 = °232).

As a result, the constitutive equation (13) is rewritten in the form

= 2(|CT<33>|,ct231> + a232)). (15)

Note that the equation (15) is similar to those discussed in [22]. However,

contrary to the work of [22], the tensor F is not isotropic, i.e. four constitutive

2d

functions are to be determined on PGS.2

Substituting the Eq. (15) into the equation (8), keeping in mind that n = k and introducing notation (...) for indices in the specific coordinate frame, we obtain the differential constraints on PGS in coordinate form

c[d{]>Fm + X{i>] - dM{3i)(x{k>,t) = 0 ((i), (j) = 1, 2, 3), (16)

In Eq. (16) values Xand d^G{3i> must be considered on PGS. Furthermore, d{j> denotes the directional derivatives:

d{3> = n ■ V, d{i> = I ■ V, d{2> = J ■ v. (17)

In Eq. (17) index (3) corresponds to the direction of the normal, and indices (1) and (2) to the tangent directions to PGS. Indepth discussion of directional operators d{j> apparatus see in monograph [28].

Note that the system of independent joint rational invariants (12) is not complete. It does not take account of joint invariants containing the squares of the force stress tensor a.

4. A full invariant formulation of differential constraint on PGS

A complete system of joint rational invariants of the second-rank tensor a and the vector n, in addition to the invariants (12), includes invariants

|t||| = |t ■ n|, |t±|2, |t±■ t±|2. (18)

In the equation (18), the vector t is defined according to

t = n • a2, t = t • T, ts = n-iG^Gis = tlGis. 2 2 2

2There will be three constitutive functions in the symmetric case. Only one constitutive function remains for the isotropic symmetric case.

Following to the discussion in the previous section, for the vector t we adopt

the decomposition

t = tii +1±. 2 2| 2

As (w is the angle determined by directions t^ and t±)

|t± ■ t±|2 = |t±|2|t±|2cos2 W

then one can replace the invariant |t^ ■ t^|2 in (18) by cos2 w.

The complete system of joint rational invariants of the stress tensor a and the vector n takes the form

|t|||, |t±|2, |t|||, |t±|2, |t± ■ t± |2. (19)

The constitutive equation (10) on PGS in terms of the complete system of joint rational invariants (19) of the asymmetric second-order tensor a and the vector n for reduced tensors onto the two-dimensional plane element T takes the form

£ = F(|t|||, |tx|2, |2||, |2±|2, |t± ■ t,|2). (20)

2d 2d 2 2 2

We choose, as in the previous case, the Cartesian rectangular coordinate system so that the unit vector k is directed along the normal n to PGS. For the projection of the traction vector t in the tangent plane T, we assume (14) to PGS, and for the vector t^ obvious equality is satisfied

|t±|2 = |t±- i|2 + (21)

The invariants (12) and (19) and projection lengths (14) and (21) are calculated via the actual stresses a on the surface £ according to

|t±-*|2 = |t±-j|2 = a2{з2), |t|| = k^L

= |ct(31) CT(11) + ct(32)ct21 + ст(33)ст(31)|2,

|t± 'J|2 = |ct(31) ct(12) + ct(32)ct(22) + ct(33)ct(32) |2,

2 2 (22)

|2||| = |ct(31)ct(13) + ct(32) ct(23) + ^(33)^

|t± ■ t±|2 = k231)CT(11) + ct(31)ct(32)ct21 + ct231)ct(33) +

+ ct(32)ct(31)ct(12) + ct232)CT(22) + ct(32)CT(33)|2.

The latter equation containing cubics in stresses should be classified as never discussed and used in the mechanics of solids.

The constitutive equation on PGS (20) taking account of the expressions (22) after introducing following notation

2 2 2

1 = F^h 11 = <31) + CT(32), 111 = |CT(31)CT(13) + ct(32)ct(23) + ct<33) |,

IV = |CT(31)CT(11) + ct(32)ct(21) + + |ct(31)ct(12) + ct(32) ct(22) + ct(33) ct(32) 12,

22

V = |CT231)CT(11) + ct(31)ct(32)ct(21) + ^31)^(33) + ct(32)ct(31)ct(12) +

+ ct232) CT(22) + ct232) CT(33)12.

can be rewritten in functional form

a = F(I, II, III, IV, V)

2d 2d

or in a simplified variant

a = F(I, II, III, IV).

2d 2d

The boundary conditions in the form of differential constraints on PGS (8), in the case when the added material has microstructural features, can be generalized by introducing in arguments of the function (10) the additional microstructural directors associated with the characteristic directions of the laying of the material in the processes of winding threads or the production of woven composites.

5. Conclusion

1. In present study, geometrically and mechanically consistent boundary conditions on PGS in the form of differential constraints have been obtained and discussed.

2. A general form of the mentioned constraints for the asymmetric force stress tensor has been obtained. It is valid for a wide range of materials and metamaterials.

3. The arguments of the constitutive tensor function on PGS have been determined by a set of invariants that are constant with respect to the rotations of the coordinate frame.

4. A geometric interpretation of the simplest variant of the differential constraint has been considered. A full invariant formulation has been proposed.

5. The developed approach must involve the experimental identification of the constitutive tensor functions on PGS.

6. The obtained results afford a general ground in applied studies on the MGS with an asymmetric force stress tensor.

Competing interests. We declare that we have no competing interests. Authors' contributions and responsibilities. Each author has participated in the article concept development and in the manuscript writing. The authors are absolutely responsible for submitting the final manuscript in print. Each author has approved the final version of manuscript.

Funding. This work was supported by the Russian Science Foundation (project no. 17-19-01257).

Acknowledgments. The authors are grateful to the reviewer for careful reading of the paper and valuable suggestions and comments.

References

1. Berman B. 3-D printing: The new industrial revolution, Business Horizons, 2012, vol. 55, no. 2, pp. 155-162. doi: 10.1016/j.bushor.2011.11.003.

2. Mankovich N. J., Cheeseman A. M., Stoker N. G. The display of three-dimensional anatomy with stereolithographic models, J. Digit. Imaging, 1990, vol.3, no. 3, pp. 200-203. doi: 10. 1007/bf03167610.

3. Stampfl J., Baudis S., Heller C., et al. Photopolymers with tunable mechanical properties processed by laser-based high-resolution stereolithography, J. Micromech. Microeng., 2008, vol.18, no. 12, 125014. doi: 10.1088/0960-1317/18/12/125014.

4. Murr L. E., Gaytan S. M., Ceylan A., et al. Characterization of titanium aluminide alloy components fabricated by additive manufacturing using electron beam melting, Acta Materialia, 2010, vol.58, no. 5, pp. 1887-1894. doi: 10.1016/j.actamat.2009.11.032.

5. Bauer W., Knitter R., Emde A., et al. Replication techniques for ceramic microcomponents with high aspect ratios, Microsystem Technologies, 2002, vol. 9, no. 1-2, pp. 81-86. doi: 10. 1007/s00542-002-0200-z.

6. Rengier F., Mehndiratta A., von Tengg-Kobligk H., et al. 3D printing based on imaging data: review of medical applications, Int. J. Comput. Assist. Radiol. Surg., 2010, vol.5, no. 4, pp. 335-341. doi: 10.1007/s11548-010-0476-x.

7. Lipson H., Kurman M. Fabricated: The new world of 3D printing. Indiana, John Wiley & Sons, 2013.

8. Ventola C. L. Medical applications for 3D printing: current and projected uses, Pharmacy and Therapeutics, 2014, vol. 39, no. 10, pp. 704-711.

9. Ozel T., Bourret G. R., Mirkin C. A. Coaxial lithography, Nature Nanotechnology, 2015, vol.10, no. 4, pp. 319-324. doi: 10.1038/nnano.2015.33.

10. Panda B., Paul S. C., Hui L. J., et al. Additive manufacturing of geopolymer for sustainable built environment, Journal of Cleaner Production, 2017, vol.167, pp. 281-288. doi: 10. 1016/j.jclepro.2017.08.165.

11. Stadnik N. E., Dats E. P. Continuum mathematical modelling of pathological growth of blood vessels, Journal of Physics: Conference Series, 2018, vol.991, 012075. doi: 10.1088/ 1742-6596/991/1/012075.

12. Stadnik N. E., Murashkin E. V., Dats E. P. Residual stresses in blood vessel wall during atherosclerosis, AIP Conference Proceedings, 2019, vol.2116, no. 1, 380013. doi: 10.1063/ 1.5114394.

13. Southwell R. V. An introduction to the theory of elasticity. For engineers and physicists, Oxford Engineering Science Series. London, Oxford Univ. Press, 1936.

14. Rashba E. I. Stresses computation in massive construction under their own weight taking into account the construction sequence, Proc. Inst. Struct. Mech. Acad. Sci. Ukrainian SSR, 1953, no. 18, pp. 23-27 (In Russian).

15. Charlab V. D. Linear creep theory of the build-up body. Mechanics of rod systems and solid mediums, In: The Proceedings of the Leningrad Construction Institute, vol.49. Leningrad, Leningrad Construction Institute, 1966, pp. 93-119 (In Russian).

16. Arutyunyan N. Kh., Naumov V. E., Radayev Yu. N. Dynamic expansion of an elastic layer. Part 1. Motion of a flow of precipitated particles at a variable rate, Izv. Akad. Nauk. Mekh. Tverd. Tela, 1992, no. 5, pp. 6-24 (In Russian).

17. Arutyunyan N. Kh., Naumov V. E., Radayev Yu. N. Dynamical expansion of an elastic layer. Part 2. The case of drop of accreted particles at a constant rate, Izv. Akad. Nauk. Mekh. Tverd. Tela, 1992, no. 6, pp. 99-112 (In Russian).

18. Naumov V. E., Radayev Yu. N. Thermomechanical model of an growing solids: Variational formulation, Preprint no. 527. Moscow, IPMech RAS, 1993, 39 pp. (In Russian)

19. Dmitrieva A. M., Naumov V. E., Radayev Yu. N. Growth of thermoelastic spherical layer: Application of the variational approach, Preprint no. 528. Moscow, IPMech RAS, 1993, 64 pp. (In Russian)

20. Arutyunyan N. Kh., Naumov V. E. The boundary value problem of the theory of viscoelastic plasticity of a growing body subject to aging, J. Appl. Math. Mech., 1984, vol.48, no. 1, pp. 1-10. doi: 10.1016/0021-8928(84)90099-6.

21. Trincher V. K. On the formulation of the problem of stresses calculation in the gravitational state of a growing solid, Izv. Akad. Nauk SSSR. Mekh. Tverd. Tela, 1984, no. 2, pp. 119-124 (In Russian).

22. Bykovtsev G. I. Izbrannye problemnye voprosy mekhaniki deformiruemykh sred [Selected Problems from Solid Mechanics. Collection of articles]. Vladivostok, Dal'nauka, 2002 (In Russian).

23. Kovalev V. A., Radayev Yu. N. Mathematical models and contemporary theories of physical fields, Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2009, vol.9, no. 4(2), pp. 41-94 (In Russian). doi: 10.18500/1816-9791-2009-9-4-2-41-94.

24. Kovalev V. A., Radayev Yu. N. Volnovye zadachi teorii polia i termomekhanika [Wave Problems of Field Theory and Thermomechanics]. Saratov, Saratov State Univ., 2010 (In Russian).

25. Kovalev V. A., Radayev Yu. N. On precisely conserved quantities of coupled micropolar thermoelastic field, Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2012, vol. 12, no. 4, pp. 71-79 (In Russian). doi: 10.18500/1816-9791-2012-12-4-71-79.

26. Kovalev V. A., Radayev Yu. N. Covariant field equations and d-tensors of hyperbolic thermoelastic continuum with fine microstructure, Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2013, vol.13, no. 2(1), pp. 60-68 (In Russian). doi: 10.18500/ 1816-9791-2013-13-2-1-60-68.

27. Gurevich G. B. Foundations of the theory of algebraic invariants. Groningen, P. Noordhoff, 1964.

28. Radayev Yu. N. Prostranstvennaia zadacha matematicheskoi teorii plastichnosti [Three-dimensional Problem of the Mathematical Theory of Plasticity]. Samara, Samara State Univ., 2007 (In Russian).

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2019. Т. 23, № 4. С. 646—656

ISSN: 2310-7081 (online), 1991-8615 (print)

d https://doi.org/10.14498/vsgtu1696

УДК 539.319

Об одном дифференциальном ограничении в континуальной механике растущих тел

Е. В. Мурашкин, Ю. Н. Радаев

Институт проблем механики им. А. Ю. Ишлинского РАН,

Россия, 119526, Москва, просп. Вернадского, 101, корп. 1.

Аннотация

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

Ключевые слова: 3Р-печать, поверхностный рост, напряжение, определяющее уравнение, рациональный инвариант, дифференциальное ограничение, полная система.

Получение: 30 апреля 2019 г. / Исправление: 12 августа 2019 г. / Принятие: 16 сентября 2019 г. / Публикация онлайн: 18 ноября 2019 г.

Конкурирующие интересы. Заявляем, что в отношении авторства и публикации этой статьи конфликта интересов не имеем.

Авторский вклад и ответственность. Все авторы принимали участие в разработке концепции статьи и в написании рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами.

Финансирование. Работа выполнена при поддержке Российского научного фонда (проект № 17-19-01257).

Благодарности. Авторы благодарят рецензента за внимательное прочтение статьи, ценные предложения и комментарии.

Научная статья

Q ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования

Murashkin E.V., Radayev Yu. N. On a differential constraint in the continuum theory of growing solids, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2019, vol. 23, no. 4, pp. 646-656. doi: 10.14498/vsgtu1696.

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

Евгений Валерьевич Мурашкин& https://orcid.org/0000-0002-3267-4742 к.ф.-м.н.; старший научный сотрудник; лаб. моделирования в механике деформируемого твердого тела; e-mail: evmurashkin@google. com

Юрий Николаевич Радаев © https://orcid.org/0000-0002-0866-2151

доктор физико-математических наук, профессор; ведущий научный сотрудник; лаб. моделирования в механике деформируемого твердого тела; e-mail: [email protected], [email protected]