Notes on Global Stress and Hyper-Stress Theories
Reuven Segev

TL;DR
This paper introduces a geometric framework for global stress and hyper-stress theories in continuum mechanics, extending classical concepts to infinite-dimensional configuration spaces and providing a unified mathematical foundation.
Contribution
It develops an infinite-dimensional geometric formulation of stress and hyper-stress theories, generalizing traditional continuum mechanics using manifold and bundle structures.
Findings
Forces are represented as linear functionals on tangent vectors.
Configuration space topology allows stress objects to be generalized.
Framework unifies finite and infinite-dimensional continuum mechanics.
Abstract
The fundamental ideas and tools of the global geometric formulation of stress and hyper-stress theory of continuum mechanics are introduced. The proposed framework is the infinite dimensional counterpart of statics of systems having finite number of degrees of freedom, as viewed in the geometric approach to analytical mechanics. For continuum mechanics, the configuration space is the manifold of embeddings of a body manifold into the space manifold. Generalized velocity fields are viewed as elements of the tangent bundle of the configuration space and forces are continuous linear functionals defined on tangent vectors, elements of the cotangent bundle. It is shown, in particular, that a natural choice of topology on the configuration space, implies that force functionals may be represented by objects that generalize the stresses of traditional continuum mechanics.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsElasticity and Material Modeling · Dynamics and Control of Mechanical Systems · Gas Dynamics and Kinetic Theory
Notes on Global Stress and Hyper-Stress Theories
Reuven Segev
Reuven Segev
Department of Mechanical Engineering
Ben-Gurion University of the Negev
Beer-Sheva, Israel
Abstract.
The fundamental ideas and tools of the global geometric formulation of stress and hyper-stress theory of continuum mechanics are introduced. The proposed framework is the infinite dimensional counterpart of statics of systems having finite number of degrees of freedom, as viewed in the geometric approach to analytical mechanics. For continuum mechanics, the configuration space is the manifold of embeddings of a body manifold into the space manifold. Generalized velocity fields are viewed as elements of the tangent bundle of the configuration space and forces are continuous linear functionals defined on tangent vectors, elements of the cotangent bundle. It is shown, in particular, that a natural choice of topology on the configuration space, implies that force functionals may be represented by objects that generalize the stresses of traditional continuum mechanics.
Key words and phrases:
Continuum mechanics; Differentiable manifold; Stress; Hyper-stress; Global analysis; Manifold of mappings; de Rham currents.
2000 Mathematics Subject Classification:
74A10; 58Z05; 58A25; 53Z05; 57N35
March 12, 2024
Contents
-
3 Banachable Spaces of Sections of Vector Bundles over Compact Manifolds
-
5.3 Open neighborhoods for using vector bundle neighborhoods
-
6.1 The case of a trivial fiber bundle—manifolds of mappings
-
8.1 Spaces of differentiable sections over a manifold without boundary and linear functionals
-
8.2 Localization of sections and linear functionals for manifolds without boundaries
-
8.3 Localization of sections and linear functionals for manifolds with corners
-
8.4 Supported sections, static indeterminacy and body forces
-
10.3.1 The inner product of a vector valued current and a vector field
-
10.3.2 The tensor product of a current and a co-vector field
-
10.3.4 The exterior product of a vector valued current and a multi-vector field
-
10.3.5 The contraction of a vector valued current and a form
-
11 The Representation of Forces by Hyper-Stresses and Non-Holonomic Stresses
1. Introduction
These notes provide an introduction to the fundamentals of global analytic continuum mechanics as developed in [ES80, Seg81, MH94, Seg86b, Seg86a, Seg16]. The terminology “global analytic” is used to imply that the formulation is based on the notion of a configuration space of the mechanical system as in analytic classical mechanics. As such, this review is complementary to that of [Seg13], which describes continuum mechanics on differentiable manifolds using a generalization of the Cauchy approach to flux and stress theory.
The setting for the basics of kinematics and statics is quite simple and provides an elegant geometric picture of mechanics. Consider the configuration space containing all admissible configuration of the system. Then, construct a differentiable manifold structure on the configuration space, define generalized (or virtual) velocities as tangent vectors, elements of , and define generalized forces as linear functions defined on the space of generalized velocities, elements of . The result of the action of a generalized force on a generalized velocity is interpreted as mechanical power. Thus, such a structure may be used to encompass both classical mechanics of mass particles and rigid bodies as well as continuum mechanics. The difference is that the configuration space for continuum mechanics and other field theories is infinite dimensional.
It is well known that the transition from the mechanics of mass particles and rigid bodies to continuum mechanics is not straightforward and requires the introduction of new notions and assumptions. The global analytic formulation explains this observation as follows. Linear functions, and forces in particular, are identically continuous when defined on a finite dimensional space. However, in the infinite dimensional situation, one has to specify exactly the topology on the infinite dimensional space of generalized velocities with respect to which forces should be continuous. Then, the properties of force functionals are deduced from the continuity requirement through a representation theorem. In other words, the properties of forces follow directly from the kinematics of the theory.
For continuum mechanics of a body in space , the basic kinematic assumption is traditionally referred to as the axiom of material impenetrability. A configuration of the body in space is specified by a mapping which is assumed to be injective and of full rank at each point—an embedding. Hence, the configuration space for continuum mechanics should be the collection of embeddings of the body manifold into the space manifold. It turns out that the -topology is the natural one to use in order to endow the collection of embeddings with a differentiable structure of a Banach manifold. The -topologies for are admissible also.
It follows that forces are continuous linear functionals on the space of vector fields over the body of class , , equipped with the -topology with a special role for the case . A standard procedure based on the Hahn-Banach theorem leads to a representation theorem for a force functional in terms of vector valued measures.
The measures representing a force generalize the stress and hyper-stress objects of continuum mechanics. On the one hand, as expected, a stress measure is not determined uniquely by a force. This is in accordance with the inherent static indeterminacy of continuum mechanics and it follows directly from the representation procedure. While the case leads to continuum mechanics of order one, the cases are extensions of higher order continuum mechanics. Thus, an existence theorem for hyper-stresses follows naturally. The relation between a force and a representing stress object is a generalization of the principle of virtual work in continuum mechanics and so it is analogous to the equilibrium equations.
The representation of forces by stress measures is significant for two reasons. First, the existence of the stress object as well as the corresponding equilibrium condition are obtained for stress distributions that may be as singular as Radon measures. In addition, while force functionals cannot be restricted to subsets of a body, measures may be restricted to subsets. This reflects a fundamental feature of stress distributions—they induce force systems on bodies. It is emphasized that in no further assumptions of mathematical or physical nature are made.
The framework described above applies to continuum mechanics on general differentiable manifolds without any additional structure such as a Riemannian metric or a connection. The body manifold is assumed here to be a compact manifold with corners. However, as described in [Mic20], it is now possible to extend the applicability of this framework to a wider class of geometric object—Whitney manifold germs.
Starting with the introduction of notation used in the manuscript in Section 2, we continue with the construction of the manifold structure on the space of embeddings. Thus, Section 3 describes the Banachable vector spaces used to construct the infinite dimensional manifold structure on the configuration space and Section 4 is concerned with the Banach manifold structure on the set of -sections of a fiber bundle . This includes, as a special case, the space of -mappings of the body into space and also provides a natural extension to continuum mechanics of generalized media. After describing the topology in in Section 5, we show in Section 6, that the set of embeddings is open in , . As such, it is a Banach manifold also and the tangent bundle is inherited from that of . In Section 7 we outline the framework for the suggested force and stress theory as described roughly above. Sections 8, 9 and 10 introduce relevant spaces of linear functionals on manifolds , and present some of their properties. These include some standard classes of functionals such as de Rham currents and Schwartz distributions on manifolds. The representation theorem of forces by stress measures in considered in Section 11. Section 12 discusses the natural situation of simple forces and stress, that is, the case .
Concluding remarks and references to further studies are made in Section 13.
2. Notation and Preliminaries
2.1. General notation
A collection of indices , will be represented as a multi-index and we will write , the length of the multi-index. In general, multi-indices will be denoted by upper-case roman letters and the associated indices will be denoted by the corresponding lower case letters. Thus, a generic element in a -multilinear mapping is given in terms of the array , . In what follows, we will use the summation convention for repeated indices as well as repeated multi-indices. Whenever the syntax is violated, e.g., when a multi-index appears more than twice in a term, it is understood that summation does not apply.
A multi-index induces a sequence in which is the number of times the index appears in the sequence . Thus, . Multi-indices may be concatenated naturally such that .
In case an array is symmetric, the independent components of the array may be listed as with . A non-decreasing multi-index, that is, amulti-index that satisfies the condition , will be denoted by boldface, upper-case roman characters so that a symmetric tensor is represented by the components , . In particular, for a function , a particular partial derivative of order is written in the form
[TABLE]
where is a non-decreasing multi-index with .
The notation , will be used for both the partial derivatives in and for the elements of the basis of the tangent space of a manifold at a point . The corresponding dual basis for will be denoted by .
Greek letters, , , will be used for strictly increasing multi-indices used in the representation of alternating tensors and forms, e.g.,
[TABLE]
Given a strictly increasing multi-index with , we will denote the strictly increasing -multi-index that complements to by . In this context, , , etc. will indicate generic increasing multi-indices. The Levi-Civita symbol will be denoted as or , so that for example , where we also set
[TABLE]
(Note that we view and as two distinct indices so summation is not implied in a term such as .)
The following identifications will be implied for tensor products of vector spaces and vector bundles
[TABLE]
For vector bundles and over a manifold , let be a section of and a section of . The notation is used for the section of given by
[TABLE]
For two manifolds and , will denote the collection of -mappings from to . If is a fiber bundle, is the space of -sections .
2.2. Manifolds with corners
Our basic object will be a fiber bundle where is assumed to be an oriented manifold with corners. We recall (e.g., [DH73, Mic80, Mel96, Lee02, MRD08]) that an -dimensional manifold with corners is a manifold whose charts assume values in the -quadrant of , that is, in
[TABLE]
In the construction of the manifold structure, it is understood that a function defined on a quadrant is said to be differentiable if it is the restriction to the quadrant of a differentiable function defined on . If is an -dimensional manifold with corners, a subset is defined to be a -dimensional, , submanifold with corners of if for any there is a chart , , such that .
With an appropriate natural definition of the integral of an -form over the boundary of a manifold with corners, Stokes’s theorem holds for manifolds with corners (see [Lee02, pp. 363–370]).
Relevant to the subject at hand is the following result (see [DH73, Mic80, Mel96, Mic20]). Every -dimensional manifold with corners is a submanifold with corners of a manifold without boundary of the same dimension. In addition, if is compact, it can be embedded as a submanifold with corners in a compact manifold without boundary of the same dimension [Mel96, pp. I.24–26]. Furthermore, -forms defined on , may be extended continuously and linearly to forms defined on . Such a manifold is referred to as* an extension* of . Each smooth vector bundle over extends to a smooth vector bundle over . Each immersion (embedding) of into a smooth manifold without boundary is the restriction of an immersion (embedding) of into .
It is emphasized that manifold with corners do not model some basic geometric shapes such as a pyramid with a rectangular base or a cone. However, much of material presented in this review is valid for a class of much more general objects,* Whitney manifold germs *as presented in [Mic20].
2.3. Bundles, jets and iterated jets
We will consider a fiber bundle , where is -dimensional and the typical fiber is -dimensional. The projection is represented locally by , , . Let
[TABLE]
be the tangent mapping represented locally by
[TABLE]
The vertical sub-bundle of is the kernel of . An element is represented locally as . With some abuse of notation, we will write both and . For with and , we may view as an element of . In other words, elements of the vertical sub-bundle are tangent vectors to that are tangent to the fibers.
Let be a section and let
[TABLE]
be the pullback of the vertical sub-bundle. Then, we may identify with the restriction of the vertical bundle to .
2.3.1. Jets
We will denote by the corresponding -jet bundle of . When no ambiguity may occur, we will often use the simpler notation and refer to a section of as a section of . One has the additional natural projections , , and in particular [Sau89]. The jet extension mapping associates with a -section, , of , a continuous section of the jet bundle .
Let be a section of which is represented locally by
[TABLE]
, . Then, denoting the -th derivative of by , a local representative of is of the form
[TABLE]
Accordingly, an element is represented locally by the coordinates
[TABLE]
2.3.2. Iterated (non-holonomic) jets
Completely non-holonomic jets for the fiber bundle are defined inductively as follows. Firstly, one defines the fiber bundles
[TABLE]
and projections
[TABLE]
Then, we define the *iterated -jet bundle *as
[TABLE]
with projection
[TABLE]
where,
[TABLE]
By induction, is a well defined fiber bundle..
When the projections are used inductively -times, we obtain a projection
[TABLE]
Let be a -section of . The iterated jet extension mapping
[TABLE]
is naturally defined by
[TABLE]
Note that we use here as a generic jet extension mapping, omitting the indication of the domain.
There is a natural inclusion
[TABLE]
Let be a vector bundle, then is a vector bundle. Continuing inductively,
[TABLE]
is a vector bundle. In this case, the inclusion is linear. Naturally, elements in the image of are referred to as holonomic.
2.3.3. Local representation of iterated jets
The local representatives of iterated jets are also constructed inductively. Hence, at each step, , to which we refer as* generation, *the number of arrays is multiplied. Hence powers of two are naturally used below. Thus, it is proposed to use multi-indices of the form , where , , etc. are binary numbers that enumerate the various arrays included in the representation. For example, a typical element of , in the form
[TABLE]
is written as
[TABLE]
and for short
[TABLE]
Here, is the generation where the -th array appears and it is given by
[TABLE]
where denotes the integer part of . In (2.23,2.24) the generations are separates by semicolons. As indicated in the example above, with each we associate a multi-index as follows. For each binary digit 1 in there is an index , . Thus, the total number of digits 1 in , which is denoted by , is the total number of indices, , in . In other words, the length, , of the induced multi-index satisfies
[TABLE]
Note also that the expression , is not a multi-index since we use upper-case letters to denote multi-indices. Here, the subscript serves for the enumeration of the indices. If no ambiguity may arise, we will often make the notation somewhat shorter and write for . Continuing by induction, let a section of be represented locally by , . Then, its -jet extension, a section of , is of the form
[TABLE]
or equivalently,
[TABLE]
where is the binary representation of . It is noted that the array contains the derivatives of the array , and that . Thus indeed, the number of digits that appear in , i.e., , determine the length of the index .
It follows that an element of may be represented in the form
[TABLE]
or
[TABLE]
That is, for each with , we have an index such that if and if .
A similar line of reasoning leads to the expression for the local representatives of the iterated jet extension mapping. For a section , the iterated jet extension , a section of , the local representation , , , satisfies
[TABLE]
Indeed, if , , with represent , then, is represented locally by
[TABLE]
Thus, by induction, any with and , may be written as , , so that .
Let an element be represented by , then is represented by .
2.4. Contraction
The right and left contractions of a -form and a -vector are given respectively by
[TABLE]
for every -vector . We will use the notation
[TABLE]
and
[TABLE]
for the mappings satisfying
[TABLE]
respectively. The left and right contraction differ by a factor of .
For the case , ; as the mappings and are injective, they are invertible. Specifically, consider the mappings
[TABLE]
and
[TABLE]
given by
[TABLE]
respectively. One can easily verify that these mappings are isomorphisms, and in fact, they are the inverses of the the contraction mappings defined above.
For example,
[TABLE]
where we view s an element of the double dual. Thus,
[TABLE]
3. Banachable Spaces of Sections of Vector Bundles over Compact Manifolds
For a compact manifold , the infinite dimensional Banach manifold of mappings to a manifold and the manifold of sections of the fiber bundle , are modeled by Banachable spaces of sections of vector bundles over , as will be described in the next section. In this section we describe the Banachable structure of such a space of differentiable vector bundle sections and make some related observations. Thus, we consider in this section a vector bundle , where is a smooth compact -dimensional manifold with corners and the typical fiber of is an -dimensional vector space. The space of -sections , , will be denoted by or by if no confusion may arise. A natural real vector space structure is induced on by setting and , .
3.1. Precompact atlases
Let , , be a finite collection of compact subsets whose interiors cover such that for each , is a subset of a domain of a chart on and
[TABLE]
is some given vector bundle chart on . Such a covering may always be found by the compactness of (using coordinate balls as, for example, in [Lee02, p. 16] or [Pal68, p. 10]). We will refer to such a structure as a precompact atlas. The same terminology will apply for the case of a fiber bundle.
3.2. The -topology on
For a section of and each , let
[TABLE]
satisfying
[TABLE]
be a local representative of .
Such a choice of a vector bundle atlas and subsets makes it possible to define, for a section ,
[TABLE]
Palais [Pal68, in particular, Chapter 4] shows that is indeed a norm endowing with a Banach space structure. The dependence of this norm on the particular choice of atlas and sets , makes the resulting space Banachable, rather than a Banach space. Other choices will correspond to different norms. However, norms induced by different choices will induce equivalent topological vector space structures on .
3.3. The jet extension mapping
Next, one observes that the foregoing may be applied, in particular, to the vector space of continuous sections of the -jet bundle of . As a continuous section of is locally of the form
[TABLE]
the analogous expression for the norm induced by a choice of a precompact vector bundle atlas is
[TABLE]
Once, the topologies of and have been defined, one may consider the jet extension mapping
[TABLE]
For a section , with local representatives , is represented by a section , the local representatives of which satisfy,
[TABLE]
Clearly, the mapping is injective and linear. Furthermore, it follows that
[TABLE]
(Note that since we take the supremum of all partial derivatives, we could replace the non-decreasing multi-index by a regular multi-index .) Thus, in view of Equation (3.4),
[TABLE]
and we conclude that is a linear embedding of into . Evidently, is not surjective as a section of need not be compatible, i.e., it need not satisfy (3.8), for some section of . As a result of the above observations, has a continuous right inverse
[TABLE]
3.4. The iterated jet extension mapping
In analogy, we now consider the iterated jet extension mapping
[TABLE]
Specializing Equation (2.31) for the case of the non-holonomic -jet bundle
[TABLE]
a section of is represented locally in the form
[TABLE]
Thus, the induced norm on is given by
[TABLE]
where the supremum is taken over all , , , with , and with .
Specializing (2.32) for the case of a vector bundle, it follows that if the section of , satisfies , its local representatives satisfy
[TABLE]
It follows that in
[TABLE]
(where the supremum is taken over all , , , with , and with ), it is sufficient to take simply all derivatives , for . Hence,
[TABLE]
where the supremum is taken over all , , , and with It is therefore concluded that
[TABLE]
In other words, one has a sequence of linear embeddings
[TABLE]
where, is the natural inclusion (2.21) and defined as , for every continuous section of , is the inclusion of sections. These embeddings are not surjective. In particular, sections of need not have the symmetry properties that hold for sections of .
4. The Construction of Charts for the Manifold of Sections
In this section, we outline the construction of charts for the Banach manifold structure on the collection of sections as in [Pal68]. (See a detailed presentation of the subject in this volume [Mic20].)
Let be a -section of . Similarly to the construction of tubular neighborhoods, the basic idea is to identify points in a neighborhood of with vectors at which are tangent to the fibers. This is achieved by defining a second order differential equation, so that a neighboring point in the same fiber as is represented through the solution of the differential equations with the initial condition by . In other words, is the image of under the exponential mapping.
To ensure that the image of the exponential mapping is located on the same fiber, , the spray inducing the second order differential equation is a vector field
[TABLE]
which is again tangent to the fiber in the sense that for
[TABLE]
This condition, together with the analog of the standard condition for a second order differential equation, namely,
[TABLE]
imply that is represented locally in the form
[TABLE]
Finally, is a *bundle spray *so that
[TABLE]
Bundle sprays can always be defined on compact manifolds using partitions of unity and the induced exponential mappings have the required properties.
The resulting structure makes it possible to identify an open neighborhood —*a vector bundle neighborhood—*of in with
[TABLE]
(We note that a rescaling is needed if is to be identified with the whole of . Otherwise, only an open neighborhood of the zero section of will be used to parametrize .)
Once the identification of with is available, the collection of sections may be identified with , . Thus, a chart into a Banachable space is constructed, where is identified with the zero section.
The construction of charts on the manifold of sections, implies that curves in in a neighborhood of are represented locally by curves in the Banachable space . Thus, tangent vectors may be identified with elements of . We therefore make the identification
[TABLE]
5. The -Topology on the Space of Sections of a Fiber
Bundle
The topology on the space of sections of fiber bundles is conveniently described in terms of filters of neighborhoods (e.g., [Trè67]).
5.1. Local representatives of sections
We consider a fiber bundle , where is assumed to be a compact manifold with corners and the typical fiber is a manifold without a boundary. Let , , and , be a precompact (as in Section 3.1) fiber bundle trivialization on . That is, the interiors of cover , and . Let , , be an atlas on so that cover .
Consider a -section . For any , we can set
[TABLE]
Let
[TABLE]
so that , and so, the local representatives of are
[TABLE]
Thus, re-enumerating the subsets and we may assume that we have a precompact trivialization , , , on , and an atlas on such that . The local representatives of relative to these atlases are
[TABLE]
5.2. Neighborhoods for and the -topology
Let be given. Consider sets of sections of the form induced by the collection of local representatives as above and numbers in the form
[TABLE]
where the supremum is taken over all
[TABLE]
The -topology on uses all such sets as a basis of neighborhoods. Using the transformation rules for the various variables, it may be shown that other choices of precompact trivialization and atlas will lead to equivalent topologies. It is noted that we use here the compactness of which implies that the weak and strong -topologies (see [Hir76, p. 35]) become identical.
Remark 5.1*.*
The collection of neighborhoods for the various values of generate a basis of neighborhoods for the topology of . If one keeps the value of , fixed, then the collection of sections
[TABLE]
where the supremum is taken over all
[TABLE]
is a neighborhood as it contains the open neighborhood . In fact, since is assumed to be compact, the collection of sets of the form is a sub-basis of neighborhoods of for the topology on .
5.3. Open neighborhoods for using vector bundle
neighborhoods
In order to specialize the preceding constructions for the case where a vector bundle neighborhood is used, we first consider local representations of sections.
Let be a section and let
[TABLE]
be the vector bundle identified with an open subbundle of . (We will use the two aspects of the vector bundle neighborhood, interchangeably.) Since the typical fiber of is , one may choose a precompact vector bundle atlas , , , on , such that
[TABLE]
Thus, if we identify all open subsets in Section 5.1 above with the typical fiber , the representatives of a section are of the form
[TABLE]
A basic neighborhood of is given by Equation (5.5). However, from the point of view of a vector bundle neighborhood, is represented by the zero section and each is viewed as a section of the vector bundle , which we may denote by . Thus
[TABLE]
In other words, using the structure of a vector bundle neighborhood we have
[TABLE]
so that is identified with a ball of radius in at the zero section.
It is concluded that the charts on induced by the vector bundle neighborhood are compatible with the -topology on .
6. The Space of Embeddings
The kinematic aspect of the Lagrangian formulation of continuum mechanics is founded on the notion of a configuration, an embeddings of a body manifold into the space manifold . The restriction of configurations to be embeddings, rather than generic -mappings of the body into space follows from the traditional principle of material impenetrability which requires that configurations be injective and that infinitesimal volume elements are not mapped into elements of zero volume.
It is noted that any configuration may be viewed as a section of the trivial fiber bundle . Thus, the constructions described above apply immediately to configurations in continuum mechanics. In this particular case, we will write for the collection of all -mappings. Our objective in this section is to describe how the set of embeddings constitutes an open subset of for . In particular, it will follow that that at each configuration , . Since the -topologies, for , are finer than the -topology, it is sufficient to prove that is open in . This brings to light the special role that the case plays in continuum mechanics.
6.1. The case of a trivial fiber bundle—manifolds of mappings
It is observed that the definitions of Sections 5.1 and 5.2 hold with natural simplifications for the case of the trivial bundle. Thus, we use a precompact atlas , , and , in (the interiors, , cover ). Given , we can find an atlas on such that . The local representatives of are of the form
[TABLE]
For the case , Equation (5.5) reduces to
[TABLE]
where the supremum is taken over all
[TABLE]
Remark 6.1*.*
It is noted that in analogy with Remark 5.1, for a fixed , a subset of mappings of the form
[TABLE]
where the supremum is taken over all
[TABLE]
is a neighborhood of as it contains a neighborhood as defined above. The collection of such sets for various values of and form a subbasis of neighborhoods for the topology on .
6.2. The space of immersions
Let be an immersion, so that is injective for every . We show that there is a neighborhood of such that all are immersions.
Note first, that since the evaluation of determinants of matrices is a continuous mapping, the collection of matrices for which all minors vanish is a closed set. Hence, the collection of all injective matrices is open in . Let be an immersion with representatives as above. For each , the derivative mapping
[TABLE]
is continuous, hence, is a compact set of injective linear mappings. Choosing any norm in , one can cover by a finite number of open balls all containing only injective mappings. In particular, setting
[TABLE]
let be the least radius of balls in this covering. Thus, we are guaranteed that any linear mapping , such that for some , is injective. Specifically, for any , if
[TABLE]
is injective everywhere in . Letting , any configuration in as in (6.2) is an immersion.
6.3. Open neighborhoods of local embeddings
Let and It is shown below that if is injective, then there is a neighborhood of mappings of such that every is injective in some fixed neighborhood of . Specifically, there is a neighborhood of , and a neighborhood of such that for each , is injective.
Let and be coordinate neighborhoods of and , respectively, such that . Let and be the local representative of and relative to these charts. Thus, we are guaranteed that
[TABLE]
By a standard corollary of the inverse function theorem, due to the injectivity of , we can choose to be small enough so that the restrictions of and to and its image under , respectively, are injective. Next, let be a neighborhood of such that is convex and its closure, , is a compact subset of . Thus, define the neighborhood whose elements, , satisfy the conditions
[TABLE]
By the definition of neighborhoods in in (6.2), contains a neighborhood of , hence, it is also a neighborhood.
Next, it is shown that the fact that the values of the derivatives of elements of are close to the injective everywhere in , implies that these mappings are close to the linear approximation using , which in turn, implies injectivity in of these elements. Specifically, for and , since
[TABLE]
the triangle inequality implies that
[TABLE]
Using the mean value theorem, there is a point such that
[TABLE]
Hence,
[TABLE]
It follows that
[TABLE]
which proves the injectivity.
6.4. Open neighborhoods of embeddings
Finally, it is shown how every has a neighborhood consisting of embeddings only. It will follow that is an open subset of . This has far-reaching consequences in continuum mechanics and it explains the special role played by the -topology in continuum mechanics.
Let be a given embedding. Using the foregoing result, for each there is an open neighborhood of and a neighborhood of , such that for each , is injective. The collection of neighborhoods , , is an open cover of and by compactness, it has a finite sub-cover. Denote the finite number of open sets of the form as above by , , so that is a compact subset of , . For each , we have a neighborhood of such that each satisfies the condition that is injective. Let so that for each , is injective for all . Let be a neighborhood of which contains only immersions as in Section 6.2. Thus, contains immersions which are locally injective.
Let be an embedding and as above. If there is no neighborhood of that contains only injective mappings, then, for each , there is a , , and points , , such that . As is a neighborhood of , we may assume that for all . By the compactness of and , we can extract a converging subsequence from the sequence in We keep the same notation for the converging subsequences and let
[TABLE]
We first exclude the possibility that . Assume , for some . Then, for any neighborhood of and any neighborhood of , there is a configuration such that is not injective. This contradicts the construction of local injectivity above.
Thus, one should consider the situation for which . Assume and for . By the definition of , the local representatives of and converge uniformly to the local representatives of and , respectively. This implies that
[TABLE]
However, since for each , , it follows that , which contradicts the assumption that is an embedding.
It is finally noted that the set of Lipschitz embeddings equipped with the Lipschitz topology may be shown to be open in the manifold of all Lipschitz mappings . See [FN05] and an application in continuum mechanics in [FS15].
7. The General Framework for Global Analytic Stress Theory
The preceding section implied that for the case where the kinematics of a material body is described by its embeddings in a physical space , the collection of configurations—the configuration space
[TABLE]
—is an open subset of the manifold of mappings , for . As a result, the configuration space is a Banach manifold in its own right and
[TABLE]
where is the natural projection of the trivial fiber bundle.
In view (4.7), , where now
[TABLE]
Hence, one may make the identifications
[TABLE]
and
[TABLE]
A section of is of the form
[TABLE]
and may be viewed as a vector field along , i.e., a mapping
[TABLE]
Thus, a tangent vector to the configuration space at the configuration may be viewed as a -vector field along . This is a straightforward generalization of the standard notion of a virtual velocity field and we summarize these observations by
[TABLE]
In the case of generalized continua, where need not be a trivial vector bundle, this simplification does not apply of course. However, the foregoing discussion motivates the definition of the configuration space for a general continuum mechanical system specified by the fiber bundle as
[TABLE]
We note that the condition that configurations are embeddings is meaningless in the case of generalized continua.
The general framework for global analytic stress theory adopts the geometric structure for the statics of systems having a finite number of degrees of freedom. Once a configuration manifold is specified, generalized or virtual velocities are defined to be elements of the tangent bundle, and generalized forces are defined to be elements of the cotangent bundle . The action of a force on a virtual velocity is interpreted as virtual power and as such, the notion of power has a fundamental role in this formulation.
The foregoing discussion, implies that a force at a configuration is an element of —a continuous and linear functional on the Banachable space of -section of the vector bundle. Thus, in the following sections we consider the properties of linear functionals on the space of -sections of a vector bundle . Of particular interest is the fact that our base manifold, or body manifold, is a manifold with corners rather than a manifold without boundary. The relation between such functionals, on the one hand, and Schwartz distribution and de Rham currents, on the other hand, is described. In Section 11 we show that the notions of stresses and hyperstresses emerge from a representation theorem for such functionals and in Section 12 we study further the properties of stresses.
8. Duals to Spaces of Differentiable Sections of a Vector Bundle: Localization
of Sections and Functionals
As follows from the foregoing discussion, generalized forces are modeled mathematically as elements of the dual space of the space of -sections of a vector bundle . This section reviews the basic notions corresponding to continuous linear functional in the dual space with particular attention to localization properties. While we assume that our base manifold is compact with corners, we want to relate the nature of functionals defined on sections over with analogous settings where is a manifold without boundary. Thus, one can make a connection of the properties of generalized forces and objects like distributions, de Rham currents and generalized sections on manifolds. (See, in particular, Section 8.5.) As an additional motivation for considering sections over manifolds without boundaries, it is observed that in both the Eulerian formulation of continuum mechanics and in classical field theories, the base manifold, either space or space-time, is usually taken as manifold without boundary. We start with the case where is a manifolds without a boundary and continue with the case where bodies are modeled by compact manifolds with corners.
8.1. Spaces of differentiable sections over a manifold without boundary
and linear functionals
A comprehensive introduction to the subject considered here is available in the Ph.D thesis [Ste00] and the corresponding [GKOS01, Chapter 3]. See also [GS77, Chapter VI] and [Kor08].
Consider the space of -sections of a vector bundle , for . For manifolds without boundary that are not necessarily compact, the setting of Section 3.2 will not give a norm on the space of sections. Thus, one extends the settings used for Schwartz distributions and de Rham currents to sections of a vector bundle (see also [Die72, Chapter XVII]). Specifically, we turn our attention to , the space of test sections—-sections of having compact supports in .
Let , be a vector bundle atlas so that
[TABLE]
and let be a compact subset of . Consider the vector subspace of sections, the supports of which are contained in . Let , indicate a finite collection of charts such that cover , and for each let , , be a fundamental sequence of compact sets, i.e., , covering . Then, for a section , the collection of semi-norms
[TABLE]
induce a Fréchet space structure for . Since for each compact subset , one has the inclusion mapping , one may define the topology on as the inductive limit topology generated by these inclusions, i.e., the strongest topology on for which all the inclusions are continuous. A sequence of sections in converges to zero, if there is a compact subset such that the supports of all sections in the sequence are contained in and the -jets of the sections converge uniformly to zero in .
A linear functional is continuous when it satisfies the following condition. Let be a sequence of sections of all of which are supported in a compact subset for some . In addition, assume that the local representatives of and their derivatives of all orders converge uniformly to zero in . Then,
[TABLE]
Functionals in for a finite value of are referred to as functionals of order .
For a linear functional , the support, is defined as follows. An open set is termed a null set of if for any section of with . The union of all null sets, is an open set which is a null set also. Thus, one defines
[TABLE]
8.2. Localization of sections and linear functionals for manifolds without
boundaries
Let be a locally finite vector bundle atlas on and consider
[TABLE]
the natural zero extension of sections supported in a compact subsets of to the space of sections that are compactly supported in . This is evidently a linear and continuous injection of the subspace. On its image, the subspace of sections with we have a left inverse, the natural restriction
[TABLE]
a surjective mapping. However, it is well known (e.g., [Sch63, Trè67, pp. 245–246]) that the inverse is not continuous.
The dual,
[TABLE]
is the restriction of functionals on to sections supported on , and as is not continuous, is not surjective (loc. cit.). We will write
[TABLE]
We also note that the restrictions satisfy the condition
[TABLE]
for any section supported in .
Consider the mapping
[TABLE]
given by
[TABLE]
Due to the overlapping between domains of definition, the mapping is not injective. However, is surjective because using a partition of unity, , which subordinate to this atlas, for each section, , is a compactly supported in and . Hence, the dual mapping,
[TABLE]
given by,
[TABLE]
is injective. In other words, a functional is determined uniquely by the collection of its restrictions. Note that no compatibility condition is imposed above on the local sections .
Since satisfy the compatibility condition (8.9), is not surjective. However, it is easy to see that is exactly the subspace of containing the compatible collections of local functionals. For let be local functionals that satisfy (8.9) and a partition of unity. Consider the functional given by
[TABLE]
If is supported in for , then
[TABLE]
Thus, is a well defined functional on and it is uniquely determined by the collection —its restrictions, independently of the partition of unity chosen.
As mentioned above, a partition of unity induces an injective right inverse to in the form
[TABLE]
that evidently satisfies . It is noted that is not a left inverse. In particular, for a section supported in , with for all ,
[TABLE]
which need not be equal to . Thus, depends on the partition of unity.
For the surjective dual mapping
[TABLE]
[TABLE]
we note that while , in general. The surjectivity of implies that every functional may be represented by a non-unique collection in the form
[TABLE]
which depends on the partition of unity. Here, denotes the functional defined by .
Nevertheless, we may restrict to the subspace of compatible local functionals, , i.e., those satisfying (8.9). Thus, the restriction
[TABLE]
is an isomorphism (which depends on the partition of unity). It follows that
[TABLE]
(For additional details, see [Die72, p. 244–245], which is restricted to the case of de Rham currents, and [GKOS01, pp. 234–235].)
8.3. Localization of sections and linear functionals for manifolds with
corners
In analogy with Section 8.2, we consider the various aspects of localization relevant to the case of compact manifolds with corners. Thus, the base manifold for the vector bundle is assumed to be a manifold with corners and we are concerned with elements of acting on sections that need not necessarily vanish together with their first jets on the boundary of .
In [Pal68, pp 10–11], Palais proves what he refers to as the “Mayer-Vietoris Theorem”. Adapting the notation and specializing the theorem to the -topology, the theorem may be stated as follows.
Theorem 8.1**.**
Let be a compact smooth manifold and let be compact -dimensional submanifolds of whose interiors cover such as in a precompact atlas. Given the vector bundle , set
[TABLE]
and define
[TABLE]
Then, is an isomorphism of Banach spaces.
We will refer to the condition in (8.23) as the compatibility condition for local representatives of sections. The most significant part of the proof is the construction of . Thus, one has to construct a field when a collection , satisfying the compatibility condition, is given. This is done using a partition of unity which is subordinate to the interiors of .
It is noted that the situation may be viewed as “dual” to that described in Section 8.2. For functionals on spaces of sections with compact supports defined on a manifold without boundary, there is a natural restriction of functionals, , and the images of a functional under the restrictions satisfy the compatibility condition (8.9). The collection of restrictions determine uniquely. Here, it follows from Theorem 8.1 that we have a natural restriction of sections, and the restricted sections satisfy the compatibility condition (8.23). The restrictions also determine the global section , uniquely.
In Section 8.2, we observed that sections with compact supports on cannot be “restricted” naturally to sections with compact supports on the various . Such restrictions depend on the chosen partition of unity. The analogous situation for functionals on manifolds with corners is described below.
Corollary 8.2**.**
Let , then may be represented (non-uniquely) by , , in the form
[TABLE]
Indeed, as in Theorem 8.1 is an embedding of into a subspace of , one has a surjective
[TABLE]
given by
[TABLE]
8.4. Supported sections, static
indeterminacy and body forces
The foregoing observations are indicative of the fundamental problem of continuum mechanics—that of static indeterminacy. Given a force on a body as an element of for some vector bundle , and a sub-body , there is no unique restriction of to a force on in . This problem is evident for standard continuum mechanics in Euclidean spaces and continues all the way to continuum mechanics of higher order on differentiable manifolds.
Adopting the notation of [Mel96], denote by the space of sections of , the -jet extensions of which vanish on all the components of the boundary . Let be a manifold without boundary extending and let
[TABLE]
be an extension of . Then, we may use zero extension to obtain an isomorphism
[TABLE]
If is a sub-body of , then, one has the inclusion
[TABLE]
The dual to the space of sections supported in is the space of extendable functionals. From [Mel96, Proposition 3.3.1] it follows that the restriction
[TABLE]
is surjective and its kernel is the space of functionals on supported in .
Thus, if we interpret as a force, is interpreted as the corresponding body force. For a sub-body , using the dual of (8.30), one has
[TABLE]
We conclude that even in this very general settings, body forces of any order may be restricted naturally to sub-bodies.
8.5. Supported functionals
Distributions on closed subsets of have been considered by Glaeser [Gla58], Malgrange [Mal66, Chapter 7] and Oksak [Oks76]. The basic tool in the analysis of distributions on closed sets is Whitney’s extension theorem [Whi34] (see also [See64, Hör90]) which guarantees that a differentiable function on a closed subset of may be extended to a compactly supported smooth function on . The extension mapping between the corresponding function spaces is continuous. The extension theorem implies that restriction of functions is surjective and so, the dual of the restriction mapping associates a unique distribution in an open subset of with a linear functional defined on the given closed set. Distributions and functionals on manifold with corners have been considered by Melrose [Mel96, Chapter 3], whom we follow below.
Thus, let and let be an extension of the vector bundle , where is a manifold without a boundary. The Whitney-Seeley extension
[TABLE]
is a continuous injection. It follows that the natural restriction
[TABLE]
its left inverse satisfying , is surjective and the inclusion
[TABLE]
is injective. In other words, each functional , determines uniquely a functional satisfying
[TABLE]
The last equation implies also that for any section supported in . Hence, is supported in .
Conversely, every , with represents a functional , i.e., . This may be deduced as follows. For any such , consider . One needs to show that . Let , then,
[TABLE]
It is observed that vanishes on so that . Since , for any section supported in . However, approximating the section , supported in the closure, , by sections supported in , one concludes that also.
Due to this construction, Melrose [Mel96, Chapter 3] refers to such functionals (distributions) as supported. It is noted that such functionals of compact support are of a finite order .
8.6. Density dual and smooth functionals
A simple example for functionals on spaces of sections of a vector bundle is provided by smooth functionals. Consider, in analogy with the dual of a vector bundle, the vector bundle of linear mappings into another one-dimensional vector vector bundle, that of -alternating tensors. Thus, for a given vector bundle, , we use the notation (see Atiyah and Bott [AB67])
[TABLE]
Let be a vector bundle morphism over . Then, in analogy with the dual mapping, one may consider
[TABLE]
It is also noted that we have
[TABLE]
and as is isomorphic with , one has a natural isomorphism
[TABLE]
For the vector bundles , ,
[TABLE]
We will refer to as the density-dual bundle and to as the density-dual mapping.
As an example, for the case , we have an isomorphism (see Section 2.4),
[TABLE]
given by
[TABLE]
Smooth functionals may be induced by smooth sections of . For a section of , and a section of , let be the -form
[TABLE]
The smooth functional induced by is defined by
[TABLE]
8.7. Generalized sections and distributions
Let be a vector bundle and consider the case where the vector bundle above is set to be
[TABLE]
Thus, the corresponding functionals on sections of are elements of
[TABLE]
In this case, smooth functionals are represented by smooth sections of
[TABLE]
One concludes that smooth functionals in are represented by sections of . It is natural therefore to refer to elements of
[TABLE]
as generalized sections of (see [AB67, GS77, GKOS01], [Kor08, p. 676]).
In the particular case where is the natural line bundle, smooth functionals are represented by real valued functions on . Consequently, elements of
[TABLE]
are referred to as generalized functions.
The apparent complication in the definition of generalized sections using the density dual is justified in the sense that each element in may be approximated by a sequence of smooth functionals induced by sections of . (See [GKOS01, p. 241].)
In the literature, the term section distributions is used in different ways in this context. For example, in [GKOS01] and [Kor08, p. 676], -valued distributions are defined as elements of , i.e., what are referred to here as generalized sections of . (In [AB67, AS68] they are referred to as distributional sections.) On the other hand, in [GS77], distributions are defined as generalized sections of —elements of . See further comments on this issue and the corresponding terms *section distributional densities *and generalized densities in [GS77, GKOS01, Kor08].
9. de Rham currents
For a manifold without boundary , de Rham currents (see [Rha84, Sch73, Fed69]) are functionals corresponding to the case of the vector bundle
[TABLE]
so that test sections are smooth -forms having compact supports. Thus, a -current of order on is a continuous linear functional on .
A particular type of -currents, smooth currents, are induced by differential -forms. Such an form, , induces the currents and by
[TABLE]
Another simple -current, is induced by an oriented -dimensional submanifold . It is naturally defined by
[TABLE]
These two examples illustrate the two points of views on currents. On the one hand, the example of the current suggests that a current in is viewed as a generalized -form. With this point of view in mind, elements of are referred to as currents of degree . Consequently, the space of -currents on is occasionally denoted by
[TABLE]
On the other hand, the example of the current induced by a -dimensional manifolds , suggests that currents be viewed as a geometric object of dimension . Thus, an element of is referred to as a p-dimensional current.
9.1. Basic operations with currents
The contraction operations of a -current and a -form , yields the -currents defined by
[TABLE]
so that
[TABLE]
Note that our notation is different from that of [Rha84] and different in sign form that of [Fed69]. In particular, given a -current , any -form induces naturally a zero-current
[TABLE]
The -current defined above can be expressed using contraction in the form
[TABLE]
For a -current and a -multi-vector field , the -currents and are defined by
[TABLE]
for a -form . Using, \xi\raisebox{0.4pt}{\mbox{\lrcorner}}\,\psi=(-1)^{qp}\psi\raisebox{0.4pt}{\mbox{\llcorner}}\,\xi, one has
[TABLE]
in analogy with the corresponding expression for multi-vectors. Thus, the wedge product of a -current and an -multi-vector is an -current. Note that a real valued function defined on may be viewed both as a zero-form and as a zero-multi-vector. Hence, we may write for any of the four operations defined above so that .
The boundary operator
[TABLE]
defined by
[TABLE]
is a linear and continuous operator. In other words, the boundary of a -current is a -current. In particular, for a smooth current, represented by the -form , one has
[TABLE]
Hence,
[TABLE]
Similarly,
[TABLE]
In order to strengthen further the point of view that a -current is a generalized -form, the exterior derivative of a -current is defined by
[TABLE]
Thus, in the smooth case,
[TABLE]
In addition, Stokes’s theorem implies that for the -current induced by the -dimensional submanifold with boundary , the boundary, a -current, is given by
[TABLE]
It is quite evident, therefore, that the notion of a boundary generalizes and unites both the exterior derivative of forms and the boundaries of manifolds.
9.2. Local representation of currents
We consider next the local representation of de Rham currents in coordinate neighborhoods.
9.2.1. Representation by [math]-currents
Let be the restriction of a -current to forms supported in a particular coordinate neighborhood—a local representative of . Writing
[TABLE]
(we could have used R\raisebox{0.4pt}{\mbox{\llcorner}}\,\textrm{d}X^{\lambda} just the same as the are real valued functions), one notes that locally
[TABLE]
Using the exterior product of a multi-vector field and a current in (9.9), we may write
[TABLE]
and so a current may be represented locally in the form
[TABLE]
This representation suggests that be interpreted as a generalized multi-vector field (cf. [Whi57]).
In the sequel, when we refer to local representative of a current , we will often keep the same notation, , and it will be implied that we consider the restriction of to forms (or sections, in general) supported in a generic coordinate neighborhood.
9.2.2. Representation by -currents
Alternatively (cf. [Rha84, p. 36]), for a -current defined in a coordinate neighborhood and with , consider the -currents
[TABLE]
Then, for every -form ,
[TABLE]
where we used
[TABLE]
Also,
[TABLE]
implies
[TABLE]
and so,
[TABLE]
as expected. Hence,
[TABLE]
and we conclude that may be represented by the -currents
[TABLE]
with
[TABLE]
This representation suggests again that a -current be interpreted as a generalized -form. In particular, an -current is a generalized function and is often referred to as a distribution on the manifold (e.g., [Mel96, Chapter 3]).
Remark 9.1*.*
It is noted that one may set
[TABLE]
Using (9.10) for the -multi-vector
[TABLE]
In addition, by (9.30) and (9.6), for the -current and the -form ,
[TABLE]
and
[TABLE]
10. Vector-valued currents
A natural extension of the notions of generalized sections and de Rham currents yields vector valued vector valued currents that will be used to model stresses. Vector valued currents and their local representations will be considered in this section.
10.1. Vector valued forms
Let be a given vector bundle whose typical fiber is -dimensional. We will refer to sections of
[TABLE]
as vector valued -forms, which is short for the more appropriate vector bundle valued -form (cf. [Sch73, p. 340])*. *Thus in particular, sections of the density dual, are vector valued forms. In the mechanical context, we will also be concerned with co-vector valued forms, that is, sections of
[TABLE]
The terminology follows from the observation that using the isomorphism induced by transposition, i.e., , a co-vector valued form may be viewed as a section of
[TABLE]
Given a co-vector valued -form, , and a vector valued -form, , one can define the bilinear action f\dot{\mbox{\wedge}}\chi by setting
[TABLE]
for sections , , , of , , , , respectively. Thus, \dot{\mbox{\wedge}} induces a bilinear mapping
[TABLE]
or a linear
[TABLE]
The mapping \dot{\mbox{\wedge}}_{\raisebox{0.4pt}{\mbox{\lrcorner}}\,} gives rise to an extension of the isomorphism e_{\raisebox{0.4pt}{\mbox{\lrcorner}}\,} considered above to an isomorphism (we keep the same notation)
[TABLE]
Let be a vector bundle trivialization for the vector bundle so that
[TABLE]
where is the -dimensional typical fiber. Given a basis in , let and be the local bases and dual bases induced by on . Then, a co-vector valued form and a vector valued -form are represented locally in the forms
[TABLE]
respectively.
10.2. Vector valued currents
We now substitute the vector bundle for the vector bundle in definition (8.50) of generalized sections. Thus,
[TABLE]
Using the isomorphism e_{\raisebox{0.4pt}{\mbox{\lrcorner}}\,} as defined above, it is concluded that we may make the identifications
[TABLE]
(see [Sch73, p. 340]). Comparing the last equation to (10.1) we may refer to elements of these spaces as *generalized vector valued -forms *or as vector valued -currents.
A smooth vector valued -current may be represented by a valued -form—a smooth section of by
[TABLE]
where it is noted that is an -form. Locally, for ,
[TABLE]
Alternatively, a smooth element of is induced by a section of in the form
[TABLE]
for every -section of . Locally, for , ,
[TABLE]
Comparing the last two expressions for the resulting densities, one concludes that
[TABLE]
Globally, it follows that
[TABLE]
where
[TABLE]
is induced by the right contraction \theta\raisebox{0.4pt}{\mbox{\llcorner}}\,\eta_{1}(\eta_{2})=\theta(\eta_{2}\wedge\eta_{1}). What determined the direction of the contraction is the choice of action of in (10.7) as in Remark 9.1.
10.3. Local representation of vector valued currents
We now consider the local representation of the restriction of a vector valued -current to vector valued forms supported in some given vector bundle chart. We introduce first some basic operations.
10.3.1. The inner product of a vector valued current and a vector field
Given a vector valued current in and a -section of , we define the (scalar) -current by
[TABLE]
For local representation, one may consider the -currents
[TABLE]
Thus, in analogy with (9.19,9.20) we have
[TABLE]
10.3.2. The tensor product of a current and a co-vector field
A (scalar) -current, , and a -section of , , induce a vector valued current g\otimes T\in C^{r}\bigl{(}W\otimes{\textstyle\bigwedge}^{p}T^{*}\mathcal{X}\bigr{)}^{*} by setting
[TABLE]
In particular, locally,
[TABLE]
Utilizing this definition, one may write for local representatives
[TABLE]
and so, complementing (10.21), one has
[TABLE]
10.3.3. Representation by [math]-currents
Proceeding as in Section 9.2, the -current may be represented by the [math]-currents
[TABLE]
Using (9.9) and (9.22), we finally have
[TABLE]
In the case where the [math]-currents are represented locally by smooth -forms , one has
[TABLE]
in accordance with (10.13).
10.3.4. The exterior product of a vector valued current and a multi-vector
field
Next, in analogy with Section 9.2, for a vector valued -current and a -multi-vector , , consider the vector valued -current defined by
[TABLE]
In particular, for multi-indices , , we define locally the vector valued -currents
[TABLE]
so that
[TABLE]
10.3.5. The contraction of a vector valued current and a form
Also, for a vector valued -current and a -form , , define the vector valued -currents \psi\raisebox{0.4pt}{\mbox{\lrcorner}}\,T and T\raisebox{0.4pt}{\mbox{\llcorner}}\,\psi as
[TABLE]
and
[TABLE]
so that T\raisebox{0.4pt}{\mbox{\llcorner}}\,\psi=(-1)^{pq}\psi\raisebox{0.4pt}{\mbox{\lrcorner}}\,T. In the case where , we obtain an element \omega\raisebox{0.4pt}{\mbox{\lrcorner}}\,T\in C^{r}(W)^{*}, a vector valued [math]-current, satisfying
[TABLE]
Locally, one may consider the functionals—vector valued [math]-currents,
[TABLE]
Hence,
[TABLE]
is a local representation of the action of using vector valued [math]-currents. It is implied by the identity , that
[TABLE]
10.3.6. Representation by -currents
Next, for a local basis of , using (10.30) and (10.34),
[TABLE]
Following the same procedure as that leading to (9.30) and (9.31), one concludes that the vector valued -current may be represented locally by the vector valued -currents in the form
[TABLE]
and
[TABLE]
Using (10.19) and (10.25), we may define the (scalar) -currents
[TABLE]
and
[TABLE]
Considering the -currents in (10.20), the local components are defined by , hence,
[TABLE]
and we conclude that
[TABLE]
Thus, Equations (10.39) and (10.40) may be rewritten as
[TABLE]
and
[TABLE]
Comparing the last equation with (10.26) we arrive at
[TABLE]
In the smooth case, the -currents are represented by functions that make up the vector valued -form
[TABLE]
Remark 10.1*.*
In summary, the representation by zero currents (e.g., (10.28) corresponds to viewing the vector valued current as an element of . On the other hand, the representation as in (10.28,10.48) corresponds to the point of view that by the isomorphism of with ,
[TABLE]
Remark 10.2*.*
In the foregoing discussion we have made special choices and used, for example, the definitions, and T^{\lambda}:=\textrm{d}X^{\lambda}\raisebox{0.4pt}{\mbox{\lrcorner}}\,T rather than and T^{\prime\lambda}:=T\raisebox{0.4pt}{\mbox{\llcorner}}\,\textrm{d}X^{\lambda}, respectively. The correspondence between the two schemes is a natural extension of Remark 9.1. In particular, will be replaced by .
11. The Representation of Forces by Hyper-Stresses and Non-Holonomic
Stresses
11.1. Stresses and non-holonomic stresses
We recall that the tangent space to the Banach manifold of -sections of the fiber bundle at the section may be identified with the Banachable space of sections of the pullback vector bundle . Elements of the tangent space at to the configuration manifold represent generalized velocities of the continuous mechanical system. Consequently, a generalized force is modeled mathematically by and element . The central message of this section is that although such functionals cannot be restricted naturally to subbodies of , as discussed in Section 8.3, forces may be represented, non-uniquely, by stress objects that enable restriction of forces to sub-bodies. In order to simplify the notation, we will consider a general vector bundle , as in Section 3, and the notation introduced there will be used throughout. The construction is analogous to representation theorem for distributions of finite order (e.g., [Sch73, p. 91] or [Trè67, p. 259]).
Consider the jet extension linear mapping
[TABLE]
as in Section 3.3. As noted, is an embedding and under the norm induced by an atlas, it is even isometric. Evidently, due to the compatibility constraint, is a proper subset of and its complement is open. Hence, the inverse
[TABLE]
is a well defined linear homeomorphism. Given a force , the linear functional
[TABLE]
is a continuous and linear functional on . Hence, by the Hahn-Banach theorem, it may be extended to a linear functional . In other words, the linear mapping
[TABLE]
is surjective.
By the definition of the dual mapping, represents a force , i.e.,
[TABLE]
if and only if,
[TABLE]
for all -virtual velocity fields . The object is interpreted as a generalization of the notion of hyper-stress in higher-order continuum mechanics and will be so referred to. For , is a generalization of the standard stress tensor. The condition (11.6), resulting from the representation theorem, is a generalization of the principle of virtual work as it states that the power expended by the force for a virtual velocity field is equal to the power expended by the hyper-stress for —containing the first derivatives of the velocity field. Accordingly, Equation (11.5) is a generalization of the equilibrium equation of continuum mechanics.
It is noted that is not unique. The non-uniqueness originates from the fact that the image of the jet extension mapping, containing the compatible jet fields, is not a dense subset of . Thus, the static indeterminacy of continuum mechanics follows naturally from the representation theorem.
In view of (3.19), the same procedure applies if we use the non-holonomic jet extension . A force may then be represented by a non-unique, non-holonomic stress in the form
[TABLE]
The mapping of Section 3.4 is an embedding. Hence, a hyper-stress may be represented by some non-unique, non-holonomic stress in the form
[TABLE]
and in the following commutative diagram all mappings are surjective.
[TABLE]
11.2. Smooth stresses
In view of the discussion in Section 8.7, hyper-stresses are elements of
[TABLE]
and so they may be approximated by smooth sections of , i.e., -forms valued in the dual of the -jet bundle.
Similarly, non-holonomic stresses are elements of
[TABLE]
and smooth non-holonomic stresses are -forms valued in the dual of the -iterated jet bundle.
11.3. Stress measures
Analytically, stresses are vector valued zero-currents that are representable by integration. (See [Fed69, Section 4.1], for the scalar case).
As noted in Section 8.2, given a vector bundle atlas , a linear functional is uniquely determined by its restrictions to sections supported in the various domains —its local representatives. In particular, for the case of an -dimensional vector bundle , and the space of functionals , a typical local representative is an element . Thus, each component is a Radon measure or a distribution representable by integration. We will use the same notation for the measure. Consequently, for a section compactly supported in , we may write
[TABLE]
Given a partition of unity subordinate to the atlas, one has
[TABLE]
For the case of stresses, one has to replace by , by , and by , . In addition, as is a manifold with corners, representing measures may be viewed as measures on the extension which are supported in . Thus, for a section of , represented locally by ,
[TABLE]
We note that the components have the same symmetry under permutations of as sections of the jet bundle. If is a section of a vector bundle , then,
[TABLE]
The same reasoning applies to the representation by non-holonomic stresses, only here we consider sections of the iterated jet bundle represented locally by , . The local non-holonomic stress measures have components and
[TABLE]
[TABLE]
where summation is implied on all values of , , for all values of such that .
It is concluded that for a given force , there is some non-unique vector valued hyper-stress measure and a non-holonomic stress measure , such that
[TABLE]
11.4. Force system induced by stresses
It was noted in Section 8.3 that given a force on a body , a manifold with corners, there is no unique way to restrict it to an -dimensional submanifold with corners, a sub-body . We view this as the fundamental problem of continuum mechanics—static indeterminacy.
Stress, though not determined uniquely by a force, provide means for inducing a force system, the assignment of a force to each sub-body . Indeed, once a stress measure is given, be it a hyper-stress or a non-holonomic stress, integration theory makes it possible to consider the force system given by
[TABLE]
for any section of .
Further details on the relation between hyper-stresses and force systems are available in [SD91]. It is our opinion that the foregoing line of reasoning captures the essence of stress theory in continuum mechanics accurately and elegantly.
12. Simple Forces and Stresses
We restrict ourselves now to the most natural setting for continuum mechanics, the case —the first value for which the set of -embeddings is open in the manifold of mappings. (See [FS15] for consideration of configurations modeled as Lipschitz mappings.) Evidently, hyper-stresses and non-holonomic stresses become identical now, and therefore, it is natural in this case to use the terminology simple forces and stresses.
12.1. Simple stresses
A simple stress on a body is an element of
[TABLE]
which implies that smooth stress distributions are sections of
[TABLE]
Following the discussion in Section 8.5, may be viewed as a generalized section of , which is supported in , where we use the extension of the vector bundle to a vector bundle over a compact manifold without boundary .
A typical local representative of a section of the jet bundle is of the form
[TABLE]
so that locally,
[TABLE]
Here, and are [math]-currents defined by
[TABLE]
In the smooth case, is represented by a section of ( in the form
[TABLE]
Locally, such a vector valued form is represented as
[TABLE]
so that, for the domain of a chart, , and a section with ,
[TABLE]
12.2. The vertical projection
The vertical sub-bundle
[TABLE]
is kernel of the natural projection
[TABLE]
In other words, elements of the vertical sub-bundle at a point are jets of sections that vanish at . Thus, if a typical element of is represented locally in the form ), an element of the vertical sub-bundle, has the form in any adapted coordinate system. The vertical sub-bundle may be identified with the vector bundle . Denoting the natural inclusion by
[TABLE]
one has the induced inclusion
[TABLE]
Clearly, is injective and a homeomorphism onto its image. Hence, its dual
[TABLE]
is a well defined surjection. Simply put, is the restriction of the stress to sections of the vertical sub-bundle. Accordingly, we will refer to an element of as a vertical stress and to as the vertical projection.
In case the stress is represented locally by the [math]-currents as in (12.4), then is represented by .
Let be a vertical stress and let . Then, defined by
[TABLE]
is a 1-current. This is an indication of the fact that may be viewed as a vector valued -current.
We may use the local representation of currents as in Section 10.3 to represent by the scalar -currents given as
[TABLE]
so that
[TABLE]
Similarly, the -current may be represented locally as in Section 9.2 by [math]-currents given by
[TABLE]
in the form
[TABLE]
Evidently,
[TABLE]
and so
[TABLE]
It is concluded that
[TABLE]
[TABLE]
where it is recalled that in case , then, .
In the smooth case, the vertical projection of the stress is represented by a section of so that
[TABLE]
where
[TABLE]
If is represented by a section of then, is represented locally by
[TABLE]
For a vertical stress which is represented by as above and a field , the -current is given locally by
[TABLE]
for a -form supported in . In other words, if the vertical stress is represented by the section of , the -current is represented by the density , a section of given by
[TABLE]
12.3. Traction stresses
Using the transposition , one has a mapping on the space of vertical stresses
[TABLE]
We define traction stress *distributions *to be elements of
[TABLE]
Thus, it is noted that a traction stress is not much different than a vertical stress distribution but transposition enables its representation as a vector valued current. Using local representation in accordance with Section 10.3.6, a traction stress is represented locally by -currents , in the form
[TABLE]
where,
[TABLE]
Hence,
[TABLE]
Introducing the notation
[TABLE]
a simple stress distribution induces a traction stress distribution by
[TABLE]
Let , then, comparing the last equation with (12.23), it is concluded that, in accordance with (10.47), locally,
[TABLE]
for any function , and so
[TABLE]
Remark 12.1*.*
Continuing Remark 10.2, it is observed that one may consider
[TABLE]
so that
[TABLE]
Thus,
[TABLE]
and comparison with (12.23) implies that
[TABLE]
where it is observed that as are zero currents, the order of the contraction and wedge product in the last equation is immaterial.
12.4. Smooth traction stresses
In the smooth case, we adapt (10.15) to the current context. The traction stress is represented by a section of so that
[TABLE]
Locally,
[TABLE]
so that the local components of are the -currents—functions—that represent locally.
Let be a smooth stress represented by the vector valued form , a section of as in (12.8, 12.9) and let be its vertical component as in (12.26). In view of (10.17, 10.16), is represented by the section of with
[TABLE]
The terminology, traction stress, originates from the fact that a traction stress represented by a smooth section of , induces the analog of a traction field on oriented hyper-surfaces in the body as follows. (See [Seg13] for further details.) We consider, for any given section of and a field , the -form given by
[TABLE]
for sections of . Consider an -dimensional oriented smooth submanifold . Let
[TABLE]
be the natural inclusion and
[TABLE]
the corresponding restriction of -forms. Combining the above, one may define a linear mapping
[TABLE]
whereby
[TABLE]
A section of is interpreted as a *surface traction distribution *on the hyper-surface . Its action is interpreted as the power density of the corresponding surface force, and may be integrated over . In particular, the condition
[TABLE]
is a generalization of Cauchy’s formula for the relation between traction fields and stresses.
Remark 12.2*.*
The factor that appears in (12.50) above, and is absent in [Seg13], originates from our use of the choice to use exterior multiplication on the left as in Remarks 9.1 and 10.2. Evidently, if we represented by the smooth vector valued form such that
[TABLE]
instead of (12.42), the factor would not appear in the analogous computation and
[TABLE]
In addition, the second of Equations (12.44), may be rewritten as
[TABLE]
12.5. The generalized divergence of the stress
Let be a stress distribution, , and . We compute using (12.33) a local expression for the boundary of the -current as
[TABLE]
Here, for a [math]-current , is the “partial boundary” operator or the dual to the partial derivative, a [math]-current defined by
[TABLE]
Since and are -currents, definition (9.16) implies that the exterior derivatives satisfy
[TABLE]
and
[TABLE]
The computations above imply that there are “dual” linear differential operators
[TABLE]
and
[TABLE]
such that
[TABLE]
Consequently, we define the generalized divergence, a differential operator
[TABLE]
by
[TABLE]
Here, we view the various terms as vector valued [math]-currents, elements of , so that each may be contracted with to give an element of Thus,
[TABLE]
The local expression for the generalized divergence in a coordinate neighborhood is obtained using (12.54). For a smooth function having a compact support in ,
[TABLE]
so that
[TABLE]
In the smooth case
[TABLE]
where we have used
[TABLE]
as is compactly supported in . Thus, noting that in the smooth case are represented by the -forms , we conclude that locally
[TABLE]
In other words, locally, the vector valued [math]-current is represented by the vector valued form
[TABLE]
12.6. The balance equation
We define now the body force current and the boundary force current corresponding to the stress , elements of , by
[TABLE]
From the definition of the divergence in (12.63) we deduce
[TABLE]
The last equation is yet another generalization of the principle of virtual work in continuum mechanics.
For the smooth case, is represented by the smooth vector valued form as in (12.42) and we can compute, for any differentiable function defined on ,
[TABLE]
where Stokes’s theorem was utilized in the fourth line and (12.50) was used in the fifth line. Thus, in the smooth case, contains, upon appropriate choices of , information regarding the action of the surface force.
12.7. Application to non-holonomic stresses
In spite of numerous attempts (see [Seg86a, Seg17, SŚ18b, SŚ18a]), for the general geometry of manifolds, we were not able to extend the foregoing analysis to hyper-stresses, even for the case of stresses represented by smooth densities. Yet, the introduction of non-holonomic stresses makes it possible to carry out one step of the reduction.
Let be a vector bundle over and consider forces in . Using the representation by non-holonomic stresses as in (11.7) in Section 11.1, let
[TABLE]
Then, a force is represented by an element
[TABLE]
in the form
[TABLE]
Thus, one may apply the foregoing analysis of simple stresses to the study the action for elements
[TABLE]
In other words, the analysis of simple stresses is used where we substitute and for and above respectively. In particular, the balance equations for this reduction will yield
[TABLE]
where
[TABLE]
are interpreted as hyper surface traction and hyper body force distributions, respectively.
13. Concluding Remarks
The forgoing text is meant to serve as an introduction to global geometric stress and hyper-stress theory. We used a simple geometric model of a mechanical system in which forces are modeled as elements of the cotangent bundle of the configuration space and outlined the necessary steps needed in order to use it in the infinite dimensional case of continuum mechanics. The traditional choice of configurations as embeddings of a body in space, led us to the natural -topology which determined the properties of forces as linear functionals. In particular, the stress object emerges from a representation theorem for force functionals.
The general stress object we obtain preserves the basic feature of the stress tensor—it induces a force system on the body and its sub-bodies as described in Section 11.4. Further details of the relation between hyper-stresses and force systems are presented in [SD91] for the general case where stresses are as irregular as measures.
Generalizing continuum mechanics to differentiable manifolds implies that derivatives can no longer be decomposed invariantly from the values of vector fields and jets, combining the values of the field and its derivatives, are used. As a result, simple stresses mix both components dual to the values of the velocity fields, , and components dual to the derivatives, . This distinction form the classical stress tensor may be treated if additional mathematical structure is introduced. It is noted that no conditions of equilibrium, which are equivalent to invariance of the virtual power under the action of the Euclidean group, were imposed. In the general case, one may assume the action of a Lie group on the space manifold and obtain corresponding balance laws (see [Seg94]).
Another subject that has been omitted here is that of constitutive relations. Constitutive relations, in particular the notion of locality have been considered from the global point of view in [Seg88]. Roughly speaking, it is shown in 11.4 that a local constitutive relation, viewed form the global point of view as a mapping that assigns a stress distribution to a configuration, which is continuous relative to the -topology is a constitutive relation for a material of grade . Thus, the notion of locality is tied in with that of continuity.
A framework for the dynamics of a continuous body, for the geometry of differentiable manifolds, was suggested in [KOS17]. The dynamics of the system is specified using a Riemannian metric on the infinite dimensional configuration space.
Acknowledgement*.*
This work has been partially supported by H. Greenhill Chair for Theoretical and Applied Mechanics and the Pearlstone Center for Aeronautical Engineering Studies at Ben-Gurion University.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[AB 67] M. Atiyah and R. Bott. A Lefschetz fixed point formula for elliptic complexes: I. Annals of Mathematics , 86:374–407, 1967.
- 2[AS 68] M. Atiyah and I. Singer. The index of elliptic operators: I. Annals of Mathematics , 87:484–530, 1968.
- 3[DH 73] A. Douady and L. Hérault. Arrondissement des variétés à coins in Appendix of A. Borel J-P. Serre: Corners and arithmetic groups. Commentarii Mathematici Helvetici , 48:484–489, 1973.
- 4[Die 72] J. Dieudonné. Treatise on Analysis , volume III. Academic Press, 1972.
- 5[ES 80] M. Epstein and R. Segev. Differentiable Manifolds and the Principle of Virtual Work in Continuum Mechanics. Journal of Mathematical Physics , 21:1243–1245, 1980.
- 6[Fed 69] H. Federer. Geometric Measure Theory . Springer, 1969.
- 7[FN 05] K. Fukui and T. Nakamura. A topological property of Lipschitz mappings. Topology and its Applications , 148:143–152, 2005.
- 8[FS 15] L. Falach and R. Segev. The configuration space and principle of virtual power for rough bodies. Mathematics and Mechanics of Solids , 20:1049–1072, 2015.
