On Finsler Geometry and Applications in Mechanics: Review and New Perspectives

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 associatedmetric tensormay 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.


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 positiondependent 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 [1][2][3], geometric field variables describing behavior of a simply connected region of a body depend fundamentally only on referential and spatial coordinate charts {  } and {  } (,  = 1, 2, . . .) 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 [7][8][9][10][11][12][13][14][15].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 2 Advances in Mathematical Physics Riemannian metrics and affine connections include [16,17].Relevant references in geometry and mathematical physics include [18][19][20][21][22][23][24][25][26], 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 [29][30][31][32].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 Kr ö ner 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 [36][37][38][39][40].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,[41][42][43].
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,[35][36][37][38]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.

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.

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  :  → [0, ∞) with the following three properties [31]: (ii) (, ) = (, ) ∀ > 0 (i.e.,  is homogeneous of degree one in ); (iii) the fundamental tensor Restriction of  to a particular tangent space    gives rise to a (local) Minkowski norm which follows from Euler's theorem and the identity Specifically letting   → d  , the length of a differential line element at  depends in general on both  and  as A 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, X = X ( Lowering and raising of indices are enabled via   and its inverse   .Nonlinear connection coefficients on  \ 0 are defined as where   =        .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 ( \ 0); its dual set {  ,   } applies for the cotangent bundle  * ( \ 0).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  \ 0. 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]  (  ) −   ∧    = 0, The first structure equation implies torsion freeness and results in The second leads to the connection coefficients When a Finsler manifold degenerates to a Riemannian manifold,    = 0 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 2.When a Finsler manifold degenerates to a locally Minkowski space ( independent of ), then Γ   =    = 0. Gradients of bases with respect to the Chern-Rund connection and Cartan tensor are As an example of covariant differentiation on a Finsler manifold with Chern-Rund connection ∇, consider a ( Notations (⋅) | and (⋅) ‖ denote respective horizontal and vertical covariant derivatives with respect to ∇.

Geometric Quantities and Identities.
Focusing again on the Chern-Rund connection ∇, curvature 2-forms are with (⋅) the exterior derivative and ∧ the wedge product (no factor of 1/2).HH-, HV-, and VV-curvature tensors of the Chern-Rund connection have respective components VV-curvature vanishes, HV-curvature obeys    =    , and a Bianchi identity for HH-curvature is When a Finsler manifold degenerates to a Riemannian manifold, then    become the components of the usual curvature tensor of Riemannian geometry constructed from    , and    = 0.All curvatures vanish in locally Minkowski spaces.It is not always possible to embed a Finsler Advances in Mathematical Physics space in a Riemannian space without torsion, but it is possible to determine the metric and torsion tensors of a space of dimension 2 − 1 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].

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 {  }, ( = 1, 2, 3) covering the body in its undeformed state.A set of director vectors D () is attached to each , where ( = 1, 2, . . ., ).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  > 1, and no fundamental function is necessarily introduced.Let  denote the spatial location in a deformed body of a point initially at , and let d () denote a deformed director vector.The deformation process is described by ( = 1, 2, 3) The deformation gradient and its inverse are The 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 d 2 = d    d  is a measure of length (  =   for Cartesian coordinates {  }).Fundamental tensors in the spatial frame describing strains of the continuum and directors are Let Γ   and Γ   denote coefficients of linear connections associated with continuum and director fields, related by where The covariant differential of a referential vector field V, where locally V() ∈   , with respect to this connection is the covariant differential of a spatial vector field , where locally () ∈   , is defined as with, for example, the differential of   (, d () ) given by   = (  /  )d  + (  /  () )d  () .Euler-Schouten tensors are Ikeda [35] implies that fundamental variables entering a field theory for directed media should include the set (F, B, E, H, K).Given the fields in (23), the kinematicgeometric theory is fully determined once Γ   are defined.The latter coefficients can be related to defect content in a crystal.For example, setting Γ   = 0 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.

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  ( = 3 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 bases Unlike {/  }, 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 is Let ∇ denote covariant differentiation with respect to a connection on the vector bundles  and , with Consider a horizontal vector field V =   (/  ) and a vertical vector field W =   (/  ).Horizontal and vertical covariant derivatives are defined as Generalization to higher-order tensor fields is given in [30].
In particular, the following coefficients are assigned to the socalled 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  = 1, 2, . . ., .A Lagrangian function L 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  = 0: 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 the usual (horizontal) deformation gradient and its inverse have components It follows from ( 45) that, similar to ( 6) and ( 7), Nonlinear connection coefficients on  and   can be related by Metric tensor components on  and   can be related by Linear connection coefficients on  and   can be related by Using the above transformations, a complete gauge Hconnection 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].

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 elasticdamage 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  = 3 + .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  = 3 + .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 transformations where in general the director deformation map [38]     =     (, ) .
Total deformation gradient F :  →   is defined as In component form, its horizontal and vertical parts are To eliminate further excessive use of Kronecker deltas, let  = 3,   =      , 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], B is associated with lattice distortion and  2 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 Γ   = 0, then F → F, the usual deformation gradient of continuum mechanics.
Restricting attention to the time-independent case, a Lagrangian is posited of the form where Q 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 define Euler-Lagrange equations referred to the reference configuration follow from the variational principle  = 0 analogously to (42).Conjugate forces to variations in kinematic and state variables can be defined as derivatives of L 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 Q and its gradient are considered in particular energy functions, for example, state variables associated with gradients of damage parameters or void volume fractions.

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.

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 D: dimensions of  and  are 3; the dimension of  is 6.Coordinates on  are {  ,   }, ( = 1, 2, 3), 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 Z  is 6.Coordinates on   are {  ,   }, where  ∈   is a point on the base body in its current configuration, and d 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 transformations where the director deformation function is, in general, Introduce the smooth scalar order parameter field  ∈  as in (59), where  = ().Here D is identified with the reference gradient of :  [38], and define the nonholonomic spatial bases as Total deformation gradient F :  →   is defined as In component form, its horizontal and vertical parts are, as in (55) and the theory of [38], As in Section 3.
Specifically for application of the theory in the context of initially homogeneous single crystals, let , (76) The 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.
Horizontal and vertical covariant derivatives in (69) reduce to partial derivatives with respect to   and   , respectively, leading to When  = (), then vertical deformation components F  =    = 0.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.

Governing Equations: Twinning
Application.Consider a crystal with a single potentially active twin system.Applying (59), let () = 1 ∀ ∈ twinned domains, () = 0 ∀ ∈ parent (original) crystal domains, and () ∈ (0, 1) ∀ ∈ in twin boundary domains.As defined in [53], let  0 denote the twinning eigenshear (a scalar constant), m =     denote the unit normal to the habit plane (i.e., the normal covector to the twin boundary), and s =   (/  ) 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, s ∈  and m ∈  *  are constant fields that obey the orthonormality conditions Twinning deformation is defined as the ( 1 1 ) tensor field Note that det Ξ = 1 +  0 ⟨s, m⟩ = 1.Scalar interpolation function () ∈ [0, 1] monotonically increases between its endpoints with increasing , and it satisfies (0) = 0, (1) = 1, and   (0) =   (1) = 0.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 A: Note that (83) can be considered a version of the Bilby-Kr ö ner 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 A 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 A (assuming    =   ‖ = 0 in (86)),  is the total interfacial energy that includes a double-well function with constant coefficient , and  2 =       is the square of the fundamental Finsler function given in (76).For isotropic twin boundary energy   =  0   , and equilibrium surface energy Γ 0 and regularization width  0 obey  0 = (3/4)Γ 0  0 and  = 12Γ 0 / 0 [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 functional (87) Recall that for solid bodies,  =  = 3.For quasi-static conditions (null kinetic energy), the Lagrangian energy density is L = −.The null first variation of the action integral, appropriate for essential (Dirichlet) boundary conditions on boundary Ω, is where the first variation of potential energy density is defined here by varying  and  within Ω, holding reference coordinates  and reference volume form dΩ fixed: (90)

Advances in Mathematical Physics
The first Piola-Kirchhoff stress tensor is P = /F, 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 ( B = constant).
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.

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.
1,  2 , . . .,   ) , Dual set {  ,    } is a local basis on  *   .Deformation of (, ) to (  ,   ) is dictated by diffeomorphisms in local coordinates: (44))These nonholonomic bases obey simple transformation laws; {/  , /   } serves as a convenient local basis for   where   =   ⊕   .(44) 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./  =   /  and   /  =   /  hold.As in Section 3,    (, ) denote nonlinear connection coefficients on , and the Finsler-type nonholonomic bases 3, a fundamental Finsler function (, ), homogeneous of degree one in , is introduced.Then [53]ication of the divergence theorem, here with vanishing variations of  and  on Ω, leads to the Euler-Lagrange equations[53]