In Finsler geometry, each point of a base manifold can be endowed with coordinates describing its position as well as a set of one or more vectors describing directions, for example. The associated metric tensor may generally depend on direction as well as position, and a number of connections emerge associated with various covariant derivatives involving affine and nonlinear coefficients. Finsler geometry encompasses Riemannian, Euclidean, and Minkowskian geometries as special cases, and thus it affords great generality for describing a number of phenomena in physics. Here, descriptions of finite deformation of continuous media are of primary focus. After a review of necessary mathematical definitions and derivations, prior work involving application of Finsler geometry in continuum mechanics of solids is reviewed. A new theoretical description of continua with microstructure is then outlined, merging concepts from Finsler geometry and phase field theories of materials science.

1. Introduction

Mechanical behavior of homogeneous isotropic elastic solids can be described by constitutive models that depend only on local deformation, for example, some metric or strain tensor that may generally vary with position in a body. Materials with microstructure require more elaborate constitutive models, for example, describing lattice orientation in anisotropic crystals, dislocation mechanisms in elastic-plastic crystals, or cracks or voids in damaged brittle or ductile solids. In conventional continuum mechanics approaches, such models typically assign one or more time- and position-dependent vector(s) or higher-order tensor(s), in addition to total deformation or strain, that describe physical mechanisms associated with evolving internal structure.

Mathematically, in classical continuum physics [13], geometric field variables describing behavior of a simply connected region of a body depend fundamentally only on referential and spatial coordinate charts and    related by a diffeomorphism , with and denoting corresponding points on the spatial and material manifolds covered by corresponding chart(s) and denoting time. State variables entering response functions depend ultimately only on material points and relative changes in their position (e.g., deformation gradients of first order and possibly higher orders for strain gradient-type models [4]). Geometric objects such as metric tensors, connection coefficients, curvature tensors, and anholonomic objects [5] also depend ultimately only on position. This is true in conventional nonlinear elasticity and plasticity theories [1, 6], as well as geometric theories incorporating torsion and/or curvature tensors associated with crystal defects, for example [715]. In these classical theories, the metric tensor is always Riemannian (i.e., essentially dependent only upon or in the spatial or material setting), meaning the length of a differential line element depends only on position; however, torsion, curvature, and/or covariant derivatives of the metric need not always vanish if the material contains various kinds of defects (non-Euclidean geometry). Connections are linear (i.e., affine). Gauge field descriptions in the context of Riemannian metrics and affine connections include [16, 17]. Relevant references in geometry and mathematical physics include [1826], in addition to those already mentioned. Finite deformation director theories of micropolar type are addressed in the context of Riemannian (as opposed to Finslerian) metrics in [1, 27].

Finsler geometry, first attributed to Finsler in 1918 [28], is more general than Riemannian geometry in the sense that the fundamental (metric) tensor generally may depend on additional independent variables, labeled here as and in spatial and material configurations, with corresponding generalized coordinates and . Formal definitions will be given later in this paper; for the present immediate discussion, it suffices to mention that each point can be considered endowed with additional degrees-of-freedom beyond or and that transformation laws among coordinates as well as connection coefficients (i.e., covariant differentials) generally depend on or as well as or . Relevant references in mathematics include [2932]. For descriptions of mechanics of solids, additional degrees-of-freedom can be associated with evolving features of the microstructure of the material, though more general physical interpretations are possible.

The use of Finsler geometry to describe continuum mechanical behavior of solids was perhaps first noted by Krner in 1968 [33] and Eringen in 1971 [3], the latter reference incorporating some basic identities and definitions derived primarily by Cartan [34], though neither developed a Finsler-based framework more specifically directed towards mechanics of continua. The first theory of Finsler geometry applied to continuum mechanics of solids with microstructure appears to be the purely kinematic theory of Ikeda [35], in a generalization of Cosserat-type kinematics whereby additional degrees-of-freedom are director vectors linked to structure. This theory was essentially extended by Bejancu [30] to distinguish among horizontal and vertical distributions of the fiber bundle of a deforming pseudo-Finslerian total space. More complete theories incorporating a Lagrangian functional (leading to physical balance or conservation laws) and couched in terms of Finsler geometry were developed by Saczuk, Stumpf, and colleagues for describing solids undergoing inelastic deformation mechanisms associated with plasticity and/or damage [3640]. To the author’s knowledge, solution of a boundary value problem in solid mechanics using Finsler geometric theory has only been reported once, in [38]. Finsler geometry has been analogously used to generalize fundamental descriptions in other disciplines of physics such as electromagnetism, quantum theory, and gravitation [30, 4143].

This paper is organized as follows. In Section 2, requisite mathematical background on Finsler geometry (sometimes called Riemann-Finsler geometry [31]) is summarized. In Section 3, the aforementioned theories from continuum physics of solids [30, 3538, 40] are reviewed and compared. In Section 4, aspects of a new theory, with a primary intention of description of structural transformation processes in real materials, are proposed and evaluated. Conclusions follow in Section 5.

2. Finsler Geometry: Background

Notation used in the present section applies to a referential description, that is, the initial state; analogous formulae apply for a spatial description, that is, a deformed body.

2.1. Coordinates and Fundamentals

Denote by an -dimensional manifold. Each element (of support) of is of the form , where and , with the tangent bundle of . A Finsler structure of is a function with the following three properties [31]:(i)The fundamental function is on ;(ii)   (i.e., is homogeneous of degree one in );(iii)the fundamental tensor is positive definite at every point of . Restriction of to a particular tangent space gives rise to a (local) Minkowski normwhich follows from Euler’s theorem and the identitySpecifically letting , the length of a differential line element at depends in general on both and asA Finsler manifold reduces to a Minkowskian manifold when does not depend on and to a Riemannian manifold when does not depend on . In the latter case, a Riemannian metric tensor is . Cartan’s tensor, with the following fully symmetric covariant components, is defined for use later:

Consider now a coordinate transformation to another chart on ; for example,From the chain rule, holonomic basis vectors on then transform as [30, 31]

2.2. Connections and Differentiation

Christoffel symbols of the second kind derived from the symmetric fundamental tensor areLowering and raising of indices are enabled via and its inverse . Nonlinear connection coefficients on are defined aswhere . The following nonholonomic bases are then introduced:It can be shown that unlike (6), these nonholonomic bases obey simple transformation laws like (7). The set serves as a convenient local basis for ; its dual set applies for the cotangent bundle . A natural Riemannian metric can then be introduced, called a Sasaki metric [31]:The horizontal subspace spanned by is orthogonal to the vertical subspace spanned by with respect to this metric. Covariant derivative , or collectively connection 1-forms , define a linear connection on pulled-back bundle over . Letting denote an arbitrary direction,A number of linear connections have been introduced in the Finsler literature [30, 31]. The Chern-Rund connection [29, 44] is used most frequently in applications related to the present paper. It is a unique linear connection on characterized by the structural equations [31]The first structure equation implies torsion freeness and results inThe second leads to the connection coefficientsWhen a Finsler manifold degenerates to a Riemannian manifold, and . Cartan’s connection 1-forms are defined by where correspond to (14); its coordinate formulae and properties are listed in [3]. It has been shown [45] how components of Cartan’s connection on a Finsler manifold can be obtained as the induced connection of an enveloping space (with torsion) of dimension . When a Finsler manifold degenerates to a locally Minkowski space ( independent of ), then . Gradients of bases with respect to the Chern-Rund connection and Cartan tensor areAs an example of covariant differentiation on a Finsler manifold with Chern-Rund connection , consider a tensor field on the manifold . The covariant differential of isNotations and denote respective horizontal and vertical covariant derivatives with respect to .

2.3. Geometric Quantities and Identities

Focusing again on the Chern-Rund connection , curvature 2-forms arewith the exterior derivative and the wedge product (no factor of ). HH-, HV-, and VV-curvature tensors of the Chern-Rund connection have respective componentsVV-curvature vanishes, HV-curvature obeys , and a Bianchi identity for HH-curvature isWhen a Finsler manifold degenerates to a Riemannian manifold, then become the components of the usual curvature tensor of Riemannian geometry constructed from , and . All curvatures vanish in locally Minkowski spaces. It is not always possible to embed a Finsler space in a Riemannian space without torsion, but it is possible to determine the metric and torsion tensors of a space of dimension in such a way that any -dimensional Finsler space is a nonholonomic subspace of such a space with torsion [46].

Nonholonomicity (i.e., nonintegrability) of the horizontal distribution is measured by [32]where is the Lie bracket and can be interpreted as components of a torsion tensor [30]. For the Chern-Rund connection [31],Since Lie bracket (21) is strictly vertical, the horizontal distribution spanned by is not involutive [31].

3. Applications in Solid Mechanics: 1973–2003

3.1. Early Director Theory

The first application of Finsler geometry to finite deformation continuum mechanics is credited to Ikeda [35], who developed a director theory in the context of (pseudo-) Finslerian manifolds. A slightly earlier work [47] considered a generalized space (not necessarily Finslerian) comprised of finitely deforming physical and geometrical fields. Paper [35] is focused on kinematics and geometry: descriptions of deformations of the continuum and the director vector fields and their possible interactions, metric tensors (i.e., fundamental tensors), and gradients of motions. Covariant differentials are defined that can be used in field theories of Finsler space [41, 48]. Essential concepts from [35] are reviewed and analyzed next.

Let denote a pseudo-Finsler manifold in the sense of [35], representative of a material body with microstructure. Let denote a material point with coordinate chart , covering the body in its undeformed state. A set of director vectors is attached to each , where . In component form, directors are written . In the context of notation in Section 2, is similar to a Finsler manifold with , though the director theory involves more degrees-of-freedom when , and no fundamental function is necessarily introduced. Let denote the spatial location in a deformed body of a point initially at , and let denote a deformed director vector. The deformation process is described by The deformation gradient and its inverse areThe following decoupled transformations are posited:with . Differentiating the second of (25),where denotes the total covariant derivative [20, 26]. Differential line and director elements can be related by

Let denote the metric in the reference configuration, such that is a measure of length ( for Cartesian coordinates ). Fundamental tensors in the spatial frame describing strains of the continuum and directors areLet and denote coefficients of linear connections associated with continuum and director fields, related bywhere and . The covariant differential of a referential vector field , where locally , with respect to this connection isthe covariant differential of a spatial vector field , where locally , is defined aswith, for example, the differential of given by . Euler-Schouten tensors areIkeda [35] implies that fundamental variables entering a field theory for directed media should include the set . Given the fields in (23), the kinematic-geometric theory is fully determined once are defined. The latter coefficients can be related to defect content in a crystal. For example, setting results in distant parallelism associated with dislocation theory [7], in which case (29) gives the negative of the wryness tensor; more general theory is needed, however, to represent disclination defects [49] or other general sources of incompatibility [6, 50] requiring a nonvanishing curvature tensor. The directors themselves can be related to lattice directions, slip vectors in crystal plasticity [51, 52], preferred directions for twinning [6, 53], or planes intrinsically prone to cleavage fracture [54], for example. Applications of multifield theory [47] towards multiscale descriptions of crystal plasticity [50, 55] are also possible.

3.2. Kinematics and Gauge Theory on a Fiber Bundle

The second known application of Finsler geometry towards finite deformation of solid bodies appears in Chapter 8 of the book of Bejancu [30]. Content in [30] extends and formalizes the description of Ikeda [35] using concepts of tensor calculus on the fiber bundle of a (generalized pseudo-) Finsler manifold. Geometric quantities appropriate for use in gauge-invariant Lagrangian functions are derived. Relevant features of the theory in [30] are reviewed and analyzed in what follows next.

Define as a fiber bundle of total space , where is the projection to base manifold and is the fiber. Dimensions of and are ( for a solid volume) and , respectively; the dimension of is . Coordinates on are , where is a point on the base body in its reference configuration, and is a director of dimension that essentially replaces multiple directors of dimension 3 considered in Section 3.1. The natural basis on is the field of frames . Let denote nonlinear connection coefficients on , and introduce the nonholonomic basesUnlike , these nonholonomic bases obey simple transformation laws; serves as a convenient local basis for adapted to a decomposition into horizontal and vertical distributions:Dual set is a local basis on . A fundamental tensor (i.e., metric) for the undeformed state isLet denote covariant differentiation with respect to a connection on the vector bundles and , withConsider a horizontal vector field and a vertical vector field . Horizontal and vertical covariant derivatives are defined asGeneralization to higher-order tensor fields is given in [30]. In particular, the following coefficients are assigned to the so-called gauge H-connection on :Comparing with the formal theory of Finsler geometry outlined in Section 2, coefficients are analogous to those of the Chern-Rund connection, and (36) is analogous to (16). The generalized pseudo-Finslerian description of [30] reduces to Finsler geometry of [31] when and a fundamental function exists from which metric tensors and nonlinear connection coefficients can be derived.

Let denote a set of differentiable state variables, where . A Lagrangian function of the following form is considered on :Let be a compact domain of , and define the functional (action integral)Euler-Lagrange equations referred to the reference configuration follow from the variational principle :These can be rewritten as invariant conservation laws involving horizontal and vertical covariant derivatives with respect to the gauge-H connection [30].

Let be the deformed image of fiber bundle , representative of deformed geometry of the body, for example. Dimensions of and are and , respectively; the dimension of is . Coordinates on are , where is a point on the base body in its current configuration, and is a director of dimension . The natural basis on is the field of frames . Let denote nonlinear connection coefficients on , and introduce the nonholonomic basesThese nonholonomic bases obey simple transformation laws; serves as a convenient local basis for whereDual set is a local basis on . Deformation of to is dictated by diffeomorphisms in local coordinates:Let the usual (horizontal) deformation gradient and its inverse have componentsIt follows from (45) that, similar to (6) and (7),Nonlinear connection coefficients on and can be related byMetric tensor components on and can be related byLinear connection coefficients on and can be related byUsing the above transformations, a complete gauge H-connection can be obtained for from reference quantities on if the deformation functions in (45) are known. A Lagrangian can then be constructed analogously to (40), and Euler-Lagrange equations for the current configuration of the body can be derived from a variational principle where the action integral is taken over the deformed space. In application of pseudo-Finslerian fiber bundle theory similar to that outlined above, Fu et al. [37] associate with the director of an oriented area element that may be degraded in strength due to damage processes such as fracture or void growth in the material. See also related work in [39].

3.3. Recent Theories in Damage Mechanics and Finite Plasticity

Saczuk et al. [36, 38, 40] adapted a generalized version of pseudo-Finsler geometry similar to the fiber bundle approach of [30] and Section 3.2 to describe mechanics of solids with microstructure undergoing finite elastic-plastic or elastic-damage deformations. Key new contributions of these works include definitions of total deformation gradients consisting of horizontal and vertical components and Lagrangian functions with corresponding energy functionals dependent on total deformations and possibly other state variables. Constitutive relations and balance equations are then derived from variations of such functionals. Essential details are compared and analyzed in the following discussion. Notation usually follows that of Section 3.2, with a few generalizations defined as they appear.

Define as a fiber bundle of total space , where is the projection to base manifold and is the fiber. Dimensions of and are 3 and , respectively; the dimension of is . Coordinates on are , where is a point on the base body in its reference configuration, and is a vector of dimension . Let be the deformed image of fiber bundle ; dimensions of and are 3 and , respectively; the dimension of is . Coordinates on are , where is a point on the base body in its current configuration, and is a vector of dimension . Tangent bundles can be expressed as direct sums of horizontal and vertical distributions:Deformation of to is locally represented by the smooth and invertible coordinate transformationswhere in general the director deformation map [38]Total deformation gradient is defined asIn component form, its horizontal and vertical parts areTo eliminate further excessive use of Kronecker deltas, let , , and . Introducing a fundamental Finsler function with properties described in Section 2.1 (), the following definitions hold [38]:Notice that correspond to the Chern-Rund connection and to the Cartan tensor. In [36], is presumed stationary () and is associated with a residual plastic disturbance, is remarked to be associated with dislocation density, and the HV-curvature tensor of (e.g., of (19)) is remarked to be associated with disclination density. In [38], is associated with lattice distortion and with residual strain energy density of the dislocation density [6, 56]. In [40], a fundamental function and connection coefficients are not defined explicitly, leaving the theory open to generalization. Note also that (52) is more general than (45). When the latter holds and , then , the usual deformation gradient of continuum mechanics.

Restricting attention to the time-independent case, a Lagrangian is posited of the formwhere is a generic vector of state variables and denotes a generic gradient that may include partial, horizontal, and/or vertical covariant derivatives in the reference configuration as physically and mathematically appropriate. Let be a compact domain of , and defineEuler-Lagrange equations referred to the reference configuration follow from the variational principle analogously to (42). Conjugate forces to variations in kinematic and state variables can be defined as derivatives of with respect to these variables. Time dependence, dissipation, first and second laws of thermodynamics, and temperature effects are also considered in [38, 40]; details are beyond the scope of this review. In the only known application of Finsler geometry to solve a boundary value problem in the context of mechanics of solids with microstructure, Stumpf and Saczuk [38] use the theory outlined above to study localization of plastic slip in a bar loaded in tension, with specifying a preferred material direction for slip. In [40], various choices of and its gradient are considered in particular energy functions, for example, state variables associated with gradients of damage parameters or void volume fractions.

4. Towards a New Theory of Structured Media

4.1. Background and Scope

In Section 4 it is shown how Finsler geometry can be applied to describe physical problems in deformable continua with evolving microstructures in a manner somewhat analogous to the phase field method. Phase field theory [57] encompasses various diffuse interface models wherein the boundary between two (or more) phases or states of material is distinguished by the gradient of a scalar field called an order parameter. The order parameter, denoted herein by , typically varies continuously between values of zero and unity in phases one and two, with intermediate values in phase boundaries. Mathematically,Physically, phases might correspond to liquid and solid in melting-solidification problems, austenite and martensite in structure transformations, vacuum and intact solid in fracture mechanics, or twin and parent crystal in twinning descriptions. Similarities between phase field theory and gradient-type theories of continuum mechanics are described in [58]; both classes of theory benefit from regularization associated with a length scale dependence of solutions that can render numerical solutions to governing equations mesh independent. Recent phase field theories incorporating finite deformation kinematics include [59] for martensitic transformations, [54] for fracture, [60] for amorphization, and [53] for twinning. The latter (i.e., deformation twinning) is the focus of a more specific description that follows later in Section 4.3 of this paper, though it is anticipated that the Finsler-type description could be adapted straightforwardly to describe other deformation physics.

Deformation twinning involves shearing and lattice rotation/reflection induced by mechanical stress in a solid crystal. The usual elastic driving force is a resolved shear stress on the habit plane, in the direction of twinning shear. Twinning can be reversible or irreversible depending on material and loading protocol; the physics of deformation twinning is described more fully in [61] and Chapter 8 of [6]. Theory that follows in Section 4.2 is thought to be the first to recognize analogies between phase field theory and Finsler geometry and that in Section 4.3 the first to apply Finsler geometric concepts to model deformation twinning. Notation follows that of Section 3.2 and Section 3.3 with possible exceptions highlighted as they appear.

4.2. Finsler Geometry and Kinematics

As in Section 3.3, define as a fiber bundle of total space with . Considered is a three-dimensional solid body with state vector : dimensions of and are 3; the dimension of is 6. Coordinates on are , (), where is a point on the base body in its reference configuration. Let be the deformed image of fiber bundle ; dimensions of and are 3, and that of is 6. Coordinates on are , where is a point on the base body in its current configuration, and is the updated state vector. Tangent bundles are expressed as direct sums of horizontal and vertical distributions as in (51):Deformation of to is locally represented by the smooth and invertible coordinate transformationswhere the director deformation function is, in general,Introduce the smooth scalar order parameter field as in (59), where . Here is identified with the reference gradient of :For simplicity, here is taken as a Cartesian coordinate chart on such that , and so forth. Similarly, is chosen Cartesian. Generalization to curvilinear coordinates is straightforward but involves additional notation. Define the partial deformation gradientBy definition and then from the chain rule,leading toFrom (63) and (65), local integrability (null curl) conditions and hold. As in Section 3, denote nonlinear connection coefficients on , and the Finsler-type nonholonomic basesLet denote nonlinear connection coefficients on [38], and define the nonholonomic spatial bases asTotal deformation gradient is defined asIn component form, its horizontal and vertical parts are, as in (55) and the theory of [38],As in Section 3.3, a fundamental Finsler function , homogeneous of degree one in , is introduced. ThenSpecifically for application of the theory in the context of initially homogeneous single crystals, letThe metric tensor with Cartesian components is identified with the gradient energy contribution to the surface energy term in phase field theory [53], as will be made explicit later in Section 4.3. From (76) and (77), reference geometry is now Minkowskian since the fundamental Finsler function is independent of . In this case, (72)–(75) and (67)-(68) reduce toHorizontal and vertical covariant derivatives in (69) reduce to partial derivatives with respect to and , respectively, leading toWhen , then vertical deformation components . In a more general version of (76) applicable to heterogeneous material properties, and is homogeneous of degree zero with respect to , and the above simplifications (i.e., vanishing connection coefficients) need not apply.

4.3. Governing Equations: Twinning Application

Consider a crystal with a single potentially active twin system. Applying (59), let twinned domains, parent (original) crystal domains, and in twin boundary domains. As defined in [53], let denote the twinning eigenshear (a scalar constant), denote the unit normal to the habit plane (i.e., the normal covector to the twin boundary), and the direction of twinning shear. In the context of the geometric framework of Section 4.2, and with simplifying assumptions (76)–(79) applied throughout the present application, and are constant fields that obey the orthonormality conditionsTwinning deformation is defined as the tensor fieldNote that . Scalar interpolation function monotonically increases between its endpoints with increasing , and it satisfies , , and . A typical example is the cubic polynomial [53, 59]The horizontal part of the deformation gradient in (79) obeys a multiplicative decomposition:The elastic lattice deformation is the two-point tensor :Note that (83) can be considered a version of the Bilby-Krner decomposition proposed for elastic-plastic solids [8, 62] and analyzed at length in [5, 11, 26] from perspectives of differential geometry of anholonomic space (neither nor is necessarily integrable to a vector field). The following energy potentials, measured per unit reference volume, are defined:Here, is the strain energy density that depends on elastic lattice deformation (assuming in (86)), is the total interfacial energy that includes a double-well function with constant coefficient , and is the square of the fundamental Finsler function given in (76). For isotropic twin boundary energy , and equilibrium surface energy and regularization width obey and [53]. The total energy potential per unit volume is . Let be a compact domain of , which can be identified as a region of the material body, and define the total potential energy functionalRecall that for solid bodies, . For quasi-static conditions (null kinetic energy), the Lagrangian energy density is . The null first variation of the action integral, appropriate for essential (Dirichlet) boundary conditions on boundary , iswhere the first variation of potential energy density is defined here by varying and within , holding reference coordinates and reference volume form fixed:Application of the divergence theorem, here with vanishing variations of and on , leads to the Euler-Lagrange equations [53]The first Piola-Kirchhoff stress tensor is , and the elastic driving force for twinning is . These equations, which specify mechanical and phase equilibria, are identical to those derived in [53], but have arrived here via use of Finsler geometry on fiber bundle . In order to achieve such correspondence, simplifications and have been applied, the first reducing the fundamental Finsler function to one of Minkowskian geometry to describe energetics of twinning in an initially homogeneous single crystal body ().

A more general and potentially powerful approach would be to generalize fundamental function and deformed coordinates to allow for all possible degrees-of-freedom. For example, a fundamental function corresponding to nonuniform values of in the vicinity of grain or phase boundaries, wherein properties change rapidly with position , could be used instead of a Minkowskian (position-independent) fundamental function . Such generalization would lead to enriched kinematics and nonvanishing connection coefficients, and it may yield new physical and mathematical insight into equilibrium equations—for example, when expressed in terms of horizontal and vertical covariant derivatives [30]—used to describe mechanics of interfaces and heterogeneities such as inclusions or other defects. Further study, to be pursued in the future, is needed to relate such a general geometric description to physical processes in real heterogeneous materials. An analogous theoretical description could be derived straightforwardly to describe stress-induced amorphization or cleavage fracture in crystalline solids, extending existing phase field models [54, 60] of such phenomena.

5. Conclusions

Finsler geometry and its prior applications towards continuum physics of materials with microstructure have been reviewed. A new theory, in general considering a deformable vector bundle of Finsler character, has been posited, wherein the director vector of Finsler space is associated with a gradient of a scalar order parameter. It has been shown how a particular version of the new theory (Minkowskian geometry) can reproduce governing equations for phase field modeling of twinning in initially homogeneous single crystals. A more general approach allowing the fundamental function to depend explicitly on material coordinates has been posited that would offer enriched description of interfacial mechanics in polycrystals or materials with multiple phases.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.