Search

For a given material, \emph{controllable deformations} are those deformations that can be maintained in the absence of body forces and by applying only boundary tractions. For a given class of materials, \emph{universal deformations} are those deformations that are controllable for any material within the class. In this paper, we characterize the universal deformations in compressible isotropic implicit elasticity defined by solids whose constitutive equations, in terms of the Cauchy stress $\boldsymbol{\sigma}$ and the left Cauchy-Green strain $\mathbf{b}$, have the implicit form $\boldsymbol{\mathsf{f}}(\boldsymbol{\sigma},\mathbf{b})=\mathbf{0}$. We prove that universal deformations are homogeneous. However, an important observation is that, unlike Cauchy (and Green) elasticity, not every homogeneous deformation is permissible for a given implicit-elastic solid. In other words, the set of universal deformations is material-dependent, yet it remains a subset of homogeneous deformations., Comment: arXiv admin note: text overlap with arXiv:2404.06235

In this paper, we formulate a continuum theory of solidification within the context of finite-strain coupled thermoelasticity. We aim to fill a gap in the existing literature, as the existing studies on solidification typically decouple the thermal problem (the classical Stefan's problem) from the elasticity problem, and often limit themselves to linear elasticity with small strains. Treating solidification as an accretion problem, with the growth velocity correlated with the jump in the heat flux across the boundary, it presents an initial boundary-value problem (IBVP) over a domain whose boundary location is a priori unknown. This IBVP is solved numerically for the specific example of radially inward solidification in a spherical container. Several parametric studies are conducted to compare the numerical results with the rigid cases in the literature and gain insights into the role of elastic deformations in solidification.

Universal displacements are those displacements that can be maintained for any member of a specific class of linear elastic materials in the absence of body forces, solely by applying boundary tractions. For linear hyperelastic (Green elastic) solids, it is known that the space of universal displacements explicitly depends on the symmetry group of the material, and moreover, the larger the symmetry group the larger the set of universal displacements. Linear Cauchy elastic solids, which include linear hyperelastic solids as a special case, do not necessarily have an underlying energy function. Consequently, their elastic constants do not possess the major symmetries. In this paper, we characterize the universal displacements of anisotropic linear Cauchy elasticity. We prove the unexpected result that for each symmetry class, the set of universal displacements of linear Cauchy elasticity is identical to that of linear hyperelasticity.

For a given class of materials, \emph{universal deformations} are those deformations that can be maintained in the absence of body forces and by applying solely boundary tractions. For inhomogeneous bodies, in addition to the universality constraints that determine the universal deformations, there are extra constraints on the form of the material inhomogeneities -- \emph{universal inhomogeneity constraints}. Those inhomogeneities compatible with the universal inhomogeneity constraints are called \emph{universal inhomogeneities}. In a Cauchy elastic solid, stress at a given point and at an instance of time is a function of strain at that point and that exact moment in time, without any dependence on prior history. A Cauchy elastic solid does not necessarily have an energy function, i.e., Cauchy elastic solids are, in general, non-hyperelastic (or non-Green elastic). In this paper we characterize universal deformations in both compressible and incompressible inhomogeneous isotropic Cauchy elasticity. As Cauchy elasticity includes hyperelasticity, one expects the universal deformations of Cauchy elasticity to be a subset of those of hyperelasticity both in the compressible and incompressible cases. It is also expected that the universal inhomogeneity constraints to be more strict than those of hyperelasticity, and hence, the set of universal inhomogeneities to be smaller than that of hyperelasticity. We prove the unexpected result that the sets of universal deformations of isotropic Cauchy elasticity and isotropic hyperelasticity are identical, in both the compressible and incompressible cases. We also prove that their corresponding universal inhomogeneities are identical as well.

It is known that the balance laws of hyperelasticity (Green elasticity), i.e., conservation of mass and balance of linear and angular momenta, can be derived using the first law of thermodynamics by postulating its invariance under rigid body motions of the Euclidean ambient space -- the Green-Naghdi-Rivlin theorem. In the case of a non-Euclidean ambient space, covariance of the energy balance -- its invariance under arbitrary diffeomorphisms of the ambient space -- gives all the balance laws and the Doyle-Ericksen formula -- the Marsden-Hughes theorem. In this note, we show that the constitutive equations as well as the balance laws of hyperelasticity can be derived using the first and second laws of thermodynamics without assuming any invariance.

In this paper we revisit the mathematical foundations of nonlinear viscoelasticity. We study the underlying geometry of viscoelastic deformations, and in particular, the intermediate configuration. Starting from the multiplicative decomposition of deformation gradient into elastic and viscous parts $\mathbf{F}=\Fe\Fv\,$, we point out that $\Fv$ can be either a material tensor ($\Fe$ is a two-point tensor) or a two-point tensor ($\Fe$ is a spatial tensor). We show that based on physical grounds the second choice is unacceptable. It is assumed that the free energy density is the sum of an equilibrium and a non-equilibrium part. The symmetry transformations and their action on the total, elastic, and viscous deformation gradients are carefully discussed. Following a two-potential approach the governing equations of nonlinear viscoelasticity are derived using the Lagrange-d'Alembert principle. We discuss the constitutive and kinetic equations for compressible and incompressible isotropic, transversely isotropic, orthotropic, and monoclinic viscoelastic solids. We finally semi-analytically study creep and relaxation in three examples of universal deformations.

In this paper a geometric field theory of dislocation dynamics and finite plasticity in single crystals is formulated. Starting from the multiplicative decomposition of the deformation gradient into elastic and plastic parts, we use Cartan's moving frames to describe the distorted lattice structure via differential $1$-forms. In this theory the primary fields are the dislocation fields, defined as a collection of differential $2$-forms. The defect content of the lattice structure is then determined by the superposition of the dislocation fields. All these differential forms constitute the internal variables of the system. The evolution equations for the internal variables are derived starting from the kinematics of the dislocation 2-forms, which is expressed using the notions of flow and of Lie derivative. This is then coupled with the rate of change of the lattice structure through Orowan's equation. The governing equations are derived using a two-potential approach to a variational principle of the Lagrange-d'Alembert type. As in the nonlinear setting the lattice structure evolves in time, the dynamics of dislocations on slip systems is formulated by enforcing some constraints in the variational principle. Using the Lagrange multipliers associated with these constraints, one obtains the forces that the lattice exerts on the dislocation fields in order to keep them gliding on some given crystallographic planes. Moreover, the geometric formulation allows one to investigate the integrability -- and hence the existence -- of glide surfaces, and how the glide motion is affected by it. Lastly, a linear theory for small dislocation densities is derived, allowing one to identify the nonlinear effects that do not appear in the linearized setting.

For a given class of materials, universal displacements are those displacements that can be maintained for any member of the class by applying only boundary tractions. In this paper we study universal displacements in compressible anisotropic linear elastic solids reinforced by a family of inextensible fibers. For each symmetry class and for a uniform distribution of straight fibers respecting the corresponding symmetry we characterize the respective universal displacements. A goal of this paper is to investigate how an internal constraint affects the set of universal displacements. We have observed that other than the triclinic and cubic solids in the other five classes (a fiber-reinforced solid with straight fibers cannot be isotropic) the presence of inextensible fibers enlarges the set of universal displacements.

In this paper we formulate a geometric nonlinear theory of the mechanics of accreting-ablating bodies. This is a generalization of the theory of accretion mechanics of Sozio and Yavari (2019). More specifically, we are interested in large deformation analysis of bodies that undergo a continuous and simultaneous accretion and ablation on their boundaries while under external loads. In this formulation the natural configuration of an accreting-ablating body is a time-dependent Riemannian 3-manifold with a metric that is an unknown a priori and is determined after solving the accretion-ablation initial-boundary-value problem. In addition to the time of attachment map, we introduce a time of detachment map that along with the time of attachment map, and the accretion and ablation velocities describes the time-dependent reference configuration of the body. The kinematics, material manifold, material metric, constitutive equations, and the balance laws are discussed in detail. As a concrete example and application of the geometric theory, we analyze a thick hollow circular cylinder made of an arbitrary incompressible isotropic material that is under a finite time-dependent extension while undergoing continuous ablation on its inner cylinder boundary and accretion on its outer cylinder boundary. The state of deformation and stress during the accretion-ablation process, and the residual stretch and stress after the completion of the accretion-ablation process are computed.

In this paper we present a large-deformation formulation of the mechanics of remodeling. Remodeling is anelasticity with an internal constraint -- material evolutions that are mass and volume preserving. In this special class of material evolutions the explicit time dependence of the energy function is via a remodeling tensor (or a set of remodeling tensors) that is (are) the internal variable(s) of the theory. The governing equations of remodeling solids are derived using a two-potential approach and the Lagrange-d'Alembert principle. We consider both isotropic and anisotropic solids and derive their corresponding remodeling equations. We study a particular remodeling of fiber-reinforced solids in which the fiber orientation is time dependent in the reference configuration -- SO(3)-remodeling. We consider an additional remodeling energy, which is motivated by the energy spent in living systems to remodel to enhance stiffness in the direction of loading. We consider the examples of a solid reinforced with either one or two families of reorienting fibers and derive their remodeling equations. This is a generalization of some of the proposed remodeling equations in the literature. We study three examples of material remodeling, namely finite extensions and torsion of solid circular cylinders, which are universal deformation for incompressible isotropic solids and certain anisotropic solids. We consider both displacement and force-control loadings. Detailed parametric studies are provided for the effects of various material and loading parameters on the fiber remodeling.

We present a general formalism for the analysis of mechanical lattices with microstructure using the concept of effective mass. We first revisit a classical case of microstructure being modeled by a spring-interconnected mass-in-mass cell. The frequency-dependent effective mass of the cell is the sum of a static mass and of an added mass, in analogy to that of a swimmer in a fluid. The effective mass is derived using three different methods: momentum equivalence, action equivalence, and dynamic condensation. These methods are generalized to mechanical systems with arbitrary microstructure. As an application, we calculate the effective mass of a $1$D composite lattice with microstructure modeled by a chiral spring-interconnected mass-in-mass cell. A reduced (condensed) model of the full lattice is then obtained by lumping the microstructure into a single effective mass. A dynamic Bloch analysis is then performed using both the full and reduced lattice models, which give the same spectral results. In particular, the frequency bands follow from the full lattice model by solving a linear eigenvalue problem, or from the reduced lattice model by solving a smaller nonlinear eigenvalue problem. The range of frequencies of negative effective mass falls within the bandgaps of the lattice. Localized modes due to defects in the microstructure have frequencies within the bandgaps, inside the negative-mass range. Defects of the outer, or macro stiffness yield localized modes within each bandgap, but outside the negative-mass range. The proposed formalism can be applied to study the odd properties of coupled micro-macro systems, e.g., active matter.

We study the geometric phases of nonlinear elastic $N$-rotors with continuous rotational symmetry. In the Hamiltonian framework, the geometric structure of the phase space is a principal fiber bundle, i.e., a base, or shape manifold~$\mathcal{B}$, and fibers $\mathcal{F}$ along the symmetry direction attached to it. The symplectic structure of the Hamiltonian dynamics determines the connection and curvature forms of the shape manifold. Using Cartan's structural equations with zero torsion we find an intrinsic (pseudo) Riemannian metric for the shape manifold. One has the freedom to define the rotation sign of the total angular momentum of the elastic rotors as either positive or negative, e.g., counterclockwise or clockwise, respectively, or viceversa. This endows the base manifold~$\mathcal{B}$ with two distinct metrics both compatible with the geometric phase. In particular, the metric is pseudo-Riemannian if $\mathsf{A}<0$, and the shape manifold is a $2$D~Robertson-Walker spacetime with positive curvature. For $\mathsf{A}>0$, the shape manifold is the hyperbolic plane $\mathbb{H}^2$ with negative curvature. We then generalize our results to free elastic $N$-rotors. We show that the associated shape manifold~$\mathcal{B}$ is reducible to the product manifold of $(N-1)$ hyperbolic planes $\mathbb{H}^2$~($\mathsf{A}>0$), or $2$D~Robertson-Walker spacetimes~($\mathsf{A}<0$) depending on the convection used to define the rotation sign of the total angular momentum. We then consider elastic $N$-rotors subject to time-dependent self-equilibrated moments. The $N$-dimensional shape manifold of the extended autonomous system has a structure similar to that of the $(N-1)$-dimensional shape manifold of free elastic rotors. The Riemannian structure of the shape manifold provides an intrinsic measure of the closeness of one shape to another in terms of curvature, or induced geometric phase.

An elastic cloak hides a hole or an inhomogeneity from elastic fields. In this paper, a formulation of the optimal design of elastic cloaks based on the adjoint state method, in which the balance of linear momentum is enforced as a constraint, is presented. The design parameters are the elastic moduli of the cloak, while the objective function is a measure of the distance between the solutions in the physical and in the virtual bodies. Both the elastic medium and the cloak are assumed to be made of isotropic linear elastic materials. In order to guarantee smooth inhomogeneous elastic moduli within the cloak a penalization term is added to the objective function. Mixed finite elements are used for discretizing the weak formulation of the optimization problem. Several numerical examples of optimal elastic cloaks designed for both single and multiple loads are presented. We consider different geometries and loading types and observe that in some cases the optimal elastic cloaks for cloaking holes (cavities) are made of auxetic materials., Comment: 28 pages, 12 figures

In this paper we investigate the possibility of elastodynamic transformation cloaking in bodies made of non-centrosymmetric gradient solids. The goal of transformation cloaking is to hide a hole from elastic disturbances in the sense that the mechanical response of a homogeneous and isotropic body with a hole covered by a cloak would be identical to that of the corresponding homogeneous and isotropic body outside the cloak. It is known that in the case of centrosymmetric gradient solids the balance of angular momentum is the obstruction to transformation cloaking. We will show that this is the case for non-centrosymmetric gradient solids as well.

In this paper we formulate the problem of elastodynamic transformation cloaking for Kirchhoff-Love plates and elastic plates with both the in-plane and out-of-plane displacements. A cloaking transformation maps the boundary-value problem of an isotropic and homogeneous elastic plate (virtual problem) to that of an anisotropic and inhomogeneous elastic plate with a hole surrounded by a cloak that is to be designed (physical problem). For Kirchhoff-Love plates, the (out-of-plane) governing equations of the virtual plate is transformed to those of the physical plate up to an unknown scalar field. In doing so, one finds the initial stress and the initial tangential body force for the physical plate, along with a set of constraints that we call cloaking compatibility equations. These constraints involve the cloaking transformation, the unknown scalar field, and the elastic constants of the virtual plate. It is noted that the cloaking map needs to satisfy certain conditions on the outer boundary of the cloak and the surface of the hole. In particular, the cloaking map needs to fix the outer boundary of the cloak up to the third order. Assuming a generic radial cloaking map, we show that cloaking a circular hole in Kirchoff-Love plates is not possible; the cloaking compatibility equations and the boundary conditions that the cloaking map needs to satisfy are the obstruction to cloaking. Next, relaxing the pure bending assumption, the transformation cloaking problem of an elastic plate in the presence of in-plane and out-of-plane displacements is formulated. In this case, there are two sets of governing equations that need to be simultaneously transformed under the cloaking map. We show that cloaking a circular hole is not possible for a general radial cloaking map; the cloaking compatibility equations and the boundary conditions that the cloaking map needs to satisfy obstruct cloaking.

A new family of mixed finite element methods$-$compatible-strain mixed finite element methods (CSFEMs)$-$are introduced for three-dimensional compressible and incompressible nonlinear elasticity. A Hu-Washizu-type functional is extremized in order to obtain a mixed formulation for nonlinear elasticity. The independent fields of the mixed formulations are the displacement, the displacement gradient, and the first Piola-Kirchhoff stress. A pressure-like field is also introduced in the case of incompressible elasticity. We define the displacement in $H^1$, the displacement gradient in $H(curl)$, the stress in $H(div)$, and the pressure-like field in $L^2$. In this setting, for improving the stability of the proposed finite element methods without compromising their consistency, we consider some stabilizing terms in the Hu-Washizu-type functional that vanish at its critical points. Using a conforming interpolation, the solution and the test spaces are approximated with some piecewise polynomial subspaces of them. In three dimensions, this requires using the N\'ed\'elec edge elements for the displacement gradient and the N\'ed\'elec face elements for the stress. This approach results in mixed finite element methods that satisfy the Hadamard jump condition and the continuity of traction on all internal faces of the mesh. This, in particular, makes CSFEMs quite efficient for modeling heterogeneous solids. We assess the performance of CSFEMs by solving several numerical examples, and demonstrate their good performance for bending problems, for bodies with complex geometries, and in the near-incompressible and the incompressible regimes. Using CSFEMs, one can capture very large strains and accurately approximate stresses and the pressure field. Moreover, in our numerical examples, we do not observe any numerical artifacts such as checkerboarding of pressure, hourglass instability, or locking.

We formulate the problems of nonlinear and linear elastodynamics transformation cloaking in a geometric framework. It is noted that a cloaking transformation is neither a spatial nor a referential change of frame (coordinates); a cloaking transformation maps the boundary-value problem of an isotropic and homogeneous elastic body (virtual problem) to that of an anisotropic and inhomogeneous elastic body with a hole surrounded by a cloak that is to be designed (physical problem). The virtual body has a desired mechanical response while the physical body is designed to mimic the same response outside the cloak using a cloaking transformation. We show that nonlinear elastodynamics transformation cloaking is not possible while nonlinear elastostatics cloaking may be possible for special deformations, e.g., radial deformations in a body with either a cylindrical or a spherical cavity. For linear elastodynamics, in line with the previous observations in the literature, we show that the elastic constants in the cloak are not fully symmetric; they do not possess the minor symmetries. We prove that transformation cloaking is not possible regardless of the shape of the hole and the cloak. We next show that linear elastodynamics transformation cloaking cannot be achieved for gradient elastic solids either; similar to classical linear elasticity the balance of angular momentum is the obstruction to transformation cloaking. We finally prove that transformation cloaking is not possible for linear elastic generalized Cosserat solids in 2D for any shape of the hole and the cloak. Particularly, in 2D, transformation cloaking cannot be achieved in linear Cosserat elasticity. We also show that transformation cloaking for a spherical cavity covered by a spherical cloak is not possible in the setting of linear elastic generalized Cosserat solids. We conjecture that this result holds for a cavity of any shape.

In this paper, we present some analytical solutions for the stress fields of nonlinear anisotropic solids with distributed line and point defects. In particular, we determine the stress fields of i) a parallel cylindrically-symmetric distribution of screw dislocations in infinite orthotropic and monoclinic media, ii) a cylindrically-symmetric distribution of parallel wedge disclinations in an infinite orthotropic medium, iii) a distribution of edge dislocations in an orthotropic medium, and iv) a spherically-symmetric distribution of point defects in a transversely isotropic spherical ball.

We study some differential complexes in continuum mechanics that involve both symmetric and non-symmetric second-order tensors. In particular, we show that the tensorial analogue of the standard grad-curl-div complex can simultaneously describe the kinematics and the kinetics of motions of a continuum. The relation between this complex and the de Rham complex allows one to readily derive the necessary and sufficient conditions for the compatibility of the displacement gradient and the existence of stress functions on non-contractible bodies. We also derive the local compatibility equations in terms of the Green deformation tensor for motions of $2$D and $3$D bodies, and shells in curved ambient spaces with constant curvatures.

A rigorous proof for convergence of the Wolf method for calculating electrostatic energy of a periodic lattice is presented. In particular, we show that for an arbitrary lattice of unit cells, the lattice sum obtained via Wolf method converges to the one obtained via Ewald method.

We impose uniform electric fields both parallel and normal to 180^o ferroelectric domain walls in PbTiO_3 and obtain the equilibrium structures using the method of anharmonic lattice statics. In addition to Ti-centered and Pb-centered perfect domain walls, we also consider Ti-centered domain walls with oxygen vacancies. We observe that electric field can increase the thickness of the domain wall considerably. We also observe that increasing the magnitude of electric field we reach a critical electric field E^c; for E > E^c there is no local equilibrium configuration. Therefore, E^c can be considered as an estimate of the threshold field E_h for domain wall motion. Our numerical results show that Oxygen vacancies decrease the value of E^c. As the defective domain walls are thicker than perfect walls, this result is in agreement with the recent experimental observations and continuum calculations that show thicker domain walls have lower threshold fields.

Using the method of anharmonic lattice statics, we calculate the equilibrium structure of steps on 180^o ferroelectric domain walls (DW) in PbTiO_3. We consider three different types of steps: i) Ti-Ti step that joins a Ti-centered DW to a Ti-centered DW, (ii) Pb-Pb step that joins a Pb-centered DW to a Pb-centered DW, and (iii) Pb-Ti step that joins a Pb-centered DW to a Ti-centered DW. We show that atomic distortions due to these steps broaden a DW but are localized, i.e., they are confined to regions with dimensions of a few lattice spacings. We see that a step locally thickens the domain wall; the defective domain wall is two to three times thicker than the perfect domain wall depending on the step type. We also observe that steps distort the polarization distribution in a mixed Bloch-N\'{e}el like way; polarization rotates out of the domain wall plane near the steps. Our calculations show that Pb-Pb steps have the lowest static energy.

In this paper we extend the classical method of lattice dynamics to defective crystals with partial symmetries. We start by a nominal defect configuration and first relax it statically. Having the static equilibrium configuration, we use a quasiharmonic lattice dynamics approach to approximate the free energy. Finally, the defect structure at a finite temperature is obtained by minimizing the approximate Helmholtz free energy. For higher temperatures we take the relaxed configuration at a lower temperature as the reference configuration. This method can be used to semi-analytically study the structure of defects at low but non-zero temperatures, where molecular dynamics cannot be used. As an example, we obtain the finite temperature structure of two 180^o domain walls in a 2-D lattice of interacting dipoles. We dynamically relax both the position and polarization vectors. In particular, we show that increasing temperature the domain wall thicknesses increase.

In this paper we first obtain the order of stress singularity for a dynamically propagating self-affine fractal crack. We then show that there is always an upper bound to roughness, i.e. a propagating fractal crack reaches a terminal roughness. We then study the phenomenon of reaching a terminal velocity. Assuming that propagation of a fractal crack is discrete, we predict its terminal velocity using an asymptotic energy balance argument. In particular, we show that the limiting crack speed is a material-dependent fraction of the corresponding Rayleigh wave speed.

In this letter we obtain the finite-temperature structure of 180^o domain walls in PbTiO3 using a quasi-harmonic lattice dynamics approach. We obtain the temperature dependence of the atomic structure of domain walls from 0K up to room temperature. We also show that both Pb-centered and Ti-centered 180^o domain walls are thicker at room temperature; domain wall thickness at T=300K is about three times larger than that of T=0K. Our calculations show that Ti-centered domain walls have a lower free energy than Pb-centered domain walls and hence are more likely to be seen at finite temperatures.

In this paper we study the effect of normal and shear strains and oxygen vacancies on the structure of 180^o ferroelectric domain walls in PbTiO_3. It is known that oxygen vacancies move to the domain walls and pin them. Hence, we assume a periodic arrangement of oxygen vacancies on both Pb-centered and Ti-centered domain walls in PbTiO_3. We use a semi-analytic anharmonic lattice statics method for obtaining the relaxed configurations using a shell potential. In agreement with recent ab initio calculations, we observe that a Pb-centered domain wall with oxygen vacancies is not stable even under strain. Our semi-analytic calculations for PbTiO_3 show that oxygen vacancies affect the structure of 180^o domain walls significantly but do not have a considerable effect on the thickness of domain walls; they broaden the domain walls by about fifty percent. We also study the effect of normal and shear strains on both perfect and defective 180^o domain walls. We observe that normal and shear strains affect the structure but do not change the domain wall thickness.

In this paper we formulate a geometric theory of thermal stresses. Given a temperature distribution, we associate a Riemannian material manifold to the body, with a metric that explicitly depends on the temperature distribution. A change of temperature corresponds to a change of the material metric. In this sense, a temperature change is a concrete example of the so-called referential evolutions. We also make a concrete connection between our geometric point of view and the multiplicative decomposition of deformation gradient into thermal and elastic parts. We study the stress-free temperature distributions of the finite-deformation theory using curvature tensor of the material manifold. We find the zero-stress temperature distributions in nonlinear elasticity. Given an equilibrium configuration, we show that a change of the material manifold, i.e. a change of the material metric will change the equilibrium configuration. In the case of a temperature change, this means that given an equilibrium configuration for a given temperature distribution, a change of temperature will change the equilibrium configuration. We obtain the explicit form of the governing partial differential equations for this equilibrium change. We also show that geometric linearization of the present nonlinear theory leads to governing equations that are identical to those of the classical linear theory of thermal stresses.

In this paper we formulate a geometric theory of the mechanics of growing solids. Bulk growth is modeled by a material manifold with an evolving metric. Time dependence of metric represents the evolution of the stress-free (natural) configuration of the body in response to changes in mass density and "shape". We show that time dependency of material metric will affect the energy balance and the entropy production inequality; both the energy balance and the entropy production inequality have to be modified. We then obtain the governing equations covariantly by postulating invariance of energy balance under time-dependent spatial diffeomorphisms. We use the principle of maximum entropy production in deriving an evolution equation for the material metric. In the case of isotropic growth, we find those growth distributions that do not result in residual stresses. We then look at Lagrangian field theory of growing elastic solids. We will use the Lagrange-d'Alembert's principle with Rayleigh's dissipation functions to derive all the governing equations. We make an explicit connection between our geometric theory and the conventional multiplicative decomposition of deformation gradient $\mathbf{F}=\mathbf{F}_e\mathbf{F}_g$ into growth and elastic parts. We linearize the nonlinear theory and derive a linearized theory of growth mechanics. Finally, we obtain the stress-free growth distributions in the linearized theory.