User:Strettle/Stress, Strain and Displacement Functions

Bold text

Stress, Strain and Displacement Functions in Linear Elasticity new article content ...

SUMMARY

Stress, strain and displacement (SSD) functions, so called because of their representation of the respective field quantities in three-dimensional linear elastic continua undergoing deformation, could be regarded as a little anachronistic in the digital age where numerical solutions to complex problems of elastic continua are routinely solved, often with just a passing understanding of the underlying mechanics. These functions are now regarded as vestiges of a century of effort that went into unlocking the equations of linear elasticity such that engineers and others could come to understand the stresses, strains and displacements that relate to particular boundary value problems of the real world. Much of the research into SSD functions for engineering purposes has now been scaled back with the advent of the emerging digital technologies. Of course mathematicians take a different view, examining the theory of linear elasticity in the broader and deeper context of the mechanics of solids which continues to the present day, in a notation (tensor indicial) that few non-mathematicians could hope to understand in any meaningful way. Although an excursion into the applications of SSD functions, in pragmatic terms, may be of limited value, they offer an illuminating pathway to an intellectual grasp of the underlying mechanics.

As mentioned, the vast majority of interested readers of this article (which includes students of engineering) would be unfamiliar with indicial notation, therefore, the multi-dimensional extension of vector notation and algebra, namely dyadic notation and algebra, whilst acknowledging its limitations, is employed for greater transparency for those readers.

INTRODUCTION

In their purest form SSD functions represent a class of functions independent of any constitutive (stress-strain) relation. Ultimately however they serve to resolve the complexity of the coupled and inter-dependent sets of equations of equilibrium and compatibility which emerge from the formulations of linear elasticity.

In determining a useful characterisation of the field quantities for three-dimensional homogeneous isotropic elastic continua it is desirable to seek an explicit function or functions in which the scalar components can be defined over the elastic region by a tractable independent set of governing equations for each of the scalar functions. That is, for the function or functions to be of any use, their ability to decouple the sets of interrelated equations of equilibrium or compatibility leading to one or more independent equations is essential. A simple example of this concept is that of the Galerkin vector.[1] which is employed later to illustrate the point.

As a second theme to this article it is shown that the known SSD functions of linear elasticity can be unified into a common format using dyadic representation. The Beltrami-Schaeffer [2][3] stress functions are recognised extensions to those of Maxwell [4] and Morera [5] for non-self-equilibrated stress fields in the linear three-dimensional elastic theory. What has generally not appeared in the literature is the notion of unification that these functions enable in respect of many other explicit functions such as Galerkin, Papkovich-Neuber[6][7] and others. The Helmholtz characterization of the displacement vector[8] for three-dimensional continua also emerges as a significant part of the unified formulation. The explicit functions can be shown to be solutions to the Beltrami-Schaeffer formulation and a conceptual hierarchical model is presented here. It is postulated that any derivative which may be described as an SSD function would likely conform to this model. Accordingly, this article is a practical attempt to describe the nature of SSD functions in the context of a unified theory within homogeneous isotropic linear elasticity.

Firstly the fundamental equations of elasticity which characterize the field quantities are required.

NOTATION AND THE BASIC EQUATIONS OF LINEAR ELASTICITY

The dyadic representation of stress and strain at a point in three-dimensional elastic space is given respectively

        

where , and are the usual unit vectors and , , ......... are unit dyads resulting from the extension into second-order formulation.

Note that we assume both tensors are symmetric. Intrinsic material properties relate these two tensors as shown in the next section.

The Constitutive Equations of Linear Elasticity

The relationship between the stress and strain tensors for an isotropic linear elastic material can be expressed in the dyadic form

where

is the trace or first invariant operator
is the idemfactor or unit tensor
in practical terms is Poisson's ratio, and
is Kirchoff's modulus (shear modulus)

The dyadic form above resolves to the six scalar equations

As an adjunct the alternative form of the constitutive equation above can be expressed

Both forms of the constitutive equation are required in the development of the SSD functions below.

The strain tensor can further be defined in terms of the displacement vector as follows

where

is the displacement vector ()
is the nabla or del operator , and
represents the deformator operator

In algebraic form the above strain-displacement relation is

The immediate corollary of defining the strain tensor via a displacement vector is that the stress tensor can also be characterized in terms of the displacement vector giving the important stress-displacement relation. In doing so a first-order reduction is achieved in the number of scalar components from six to three (strains to displacements) in which a potential solution to the boundary value problems becomes more tractable. The relation becomes

We defer further analysis of this for now.

The Equations of Equilibrium

Following the dyadic form the equations of equilibrium can be expressed as

where

is the body force vector

The above embodies the following scalar equilibrium equations in the three-dimensional domain

Substitution of the constitutive relation for stress into the stress-equilibrium equation yields (homogeneity of G not assumed)

In tensorial form we have the equivalent of three independent (but coupled) governing scalar equations embodying six scalar components of strain.

With further inclusion of the strain-displacement relation this leads to the important displacement-equilibrium equation

For homogeneous isotropic linear elastic domains the above resolves to Navier's [9]equation

with the constant material property absorbed into

The Equations of Compatibility

The strain equations of compatibility are generally referred to as the compatibility conditions and may be written in the dyadic form

In scalar form these may be written

The strain-compatibility equation is immediately satisfied by substitution of the strain-displacement relation

The curl of the grad of the displacement vector immediately resolves the above by definition. This brings us to the notion that the displacement vector as defined is in fact a special strain function which, together with a deformator operator, identically satisfies the compatibility conditions. The displacement vector is thus one of the few functions of elastic theory having a physical as well as mathematical property. This is a precursor to still higher-order functions. For example, the substitution of the strain-displacement relation into the formulation enables a reduction in the number of unknown variables from six (strains) to three (displacements). By extension, higher-order functions can further reduce the number of unknowns from three to one (stress or strain function).

If the constitutive relation is now introduced into the strain-compatibility equation without assuming homogeneity of the elastic domain we get the stress-compatibility equation

Noting that the stress-compatibilty equation is an identity then it's trace must also be zero, that is (ignoring constants)

Since , this gives the interesting result that the trace of the stress tensor is harmonic, that is .

Employing the general identity (work it out long hand for proof!)

then, again since , and assuming is symmetric, the stress-compatibility equation, when homogeneity is imposed, reduces to

In this form the equation is known as the Beltrami-Michell[10] compatibility equation. See for example,Nadeau[11]

For completeness the tensorial form of the Beltrami-Michell equation resolves in scalar form to

In scalar form it is clear that managing these equations is impractical and confirms the need for the tensorial or dyadic form to gain any traction in finding solutions to them. Also, at this juncture, it is assumed that the inclusion of the body force vector leads to unnecessary complexity in furthering the underlying theory of SSD functions and it is ignored in the subsequent formulations.

Coupling and Dependency in the Equilibrium and Compatibility Equations

Before delving further into possible solutions to Navier's displacement-equilibrium equations or the Beltrami-Michell stress-compatibility equations it is worth visiting the notions of coupling and dependency within the equations which lead us to require higher-order functions to enable us to find tractability. That is, coupling and dependency within the equations force us to seek higher-order functions to enable a resolution of boundary value problem formulations.

A distinction is made between coupling and dependency. For example the six scalar components of the Beltami-Michell equation which define the six components of stress over the elastic domain exhibit both coupling and dependency. By contrast there are only three scalar components of Navier's displacement-equilibrium equations. These equations exhibit coupling but not dependency.

Coupling between equations is demonstrated when one scalar equation contains more than one scalar component of stress or displacement. This results in the 'simultaneous' nature of the equations. Naviers equations though coupled are independent since each represents a unique condition to be satisfied. With the Beltami-Michell equations from any two a third can be derived. This is the result of dependency. No such manipulation can be applied to Navier's equations.

For the Beltrami-Michell equations there are effectively the equivalent of three independent equations only to define the six scalar components of stress. The addition of the stress-equilibrium equations to these gives us a total of six coupled but mutually independent equations which are necessary for the complete description of the six components of stress. However if one is to find a solution to the elastic boundary value problem based on for example the Beltrami-Michell equations only (a sufficient condition) then the equivalent of three scalar equations only, together with the six components of stress, requires that higher-order functions be introduced to reduce the number of independent scalar variables to three. This means that the stress components must be characterized through higher order functions. This is the principle justification for the introduction of the SSD functions as a means of obtaining the general solutions to the basic equations of elasticity.

The simplest example of this concept is briefly examined below in relation to the so-called Galerkin vector which is able to decouple Navier's displacement-equilibrium equations.

Westergaard[12] after Galerkin showed that the displacement vector could be represented as follows

where is an explicit displacement vector function otherwise know as the Galerkin vector

Substituting the above into Navier's equation yields a tractable set of independent equations in biharmonic form for each of the three scalar components of , namely

Presto! the equilibrium equations are decoupled. The trade-off for this advantage is the requirement to operate with fourth-order functions and to formulate any boundary conditions in terms of the Galerkin functions. However the point is made and we now examine these and other types of higher-order functions in more detail.

THE BELTRAMI-SCHAEFFER FUNCTIONS

The Beltrami-Schaeffer functions are an extension to the generalized form of the complementary sets of Maxwell and Morera functions. The latter are defined in tensorial form as

        

The stress tensor may be characterized respectively by Maxwell's and Morera's functions respectively as

   and   

and in scalar form

Beltrami's generalization of the stress tensor in terms of these functions is

where is the combined Maxwell-Morera stress function tensor. This expression serves as a complete description of for totally self-equilibrated stress fields.

The embodied scalar equations become

Schaeffer introduced a modification to the functions in the form of an additional harmonic vector and arrived at what has become the Beltrami-Schaeffer functions. These may be written

where is a vector which accounts for the non-self-equilibrated component of the stress tensor and is harmonic since .

At this point we introduce a special vectorial operator in shorthand notation as follows . Note that .

The expanded form of the above may be written

and the trace of becomes

Noting that the trace of the constitutive equation for strain, is

then through the equivalence

This expression lends itself to integration which finally defines the displacement vector in terms of the Beltrami-Schaeffer functions as follows

where the vector constant of integration, , has been independently maintained (not absorbed into ) and scaled to better fit the formulation. This is important. As will be demonstrated below, the vector of integration is maintained to ensure that the explicit functions of Galerkin, Papkovich-Neuber and others conform to Beltrami-Schaeffer format. This is particularly so with the Helmholtz functions where directly becomes one of the functions.

To obtain a complete formulation of the field quantities through these functions it remains to deduce the governing tensorial equation for .

Employing the stress-displacement equation

substitution for above yields

If this is now equated to the initial format for the stress tensor, namely

we can generate the underlying tensorial equation which embodies and as follows

A derived solution to this equation for the boundary value problems of linear elasticity requires great insight into the algebra of tensor equations, nevertheless it will be shown that the functions of Galerkin, Papkovich-Neuber and others are clearly explicit solutions which, upon substitution, yield the requisite identities. In this respect the above equation is the generalized form which embodies the characteristics of all equations defining the known functions.

For completeness an independent equation for can be generated by twice applying the divergence operator to the above, yielding (remembering that is harmonic)

divergence

divergence

Accordingly, is neither harmonic nor biharmonic, rendering the function itself quite useless in practical terms. The underlying tensorial equation cannot be manipulated to obtain a simple governing equation for as one can for . What does do is represent, with , a pair of functions that encode all other SSD functions. It is worth noting that the generality of the Beltrami-Schaeffer functions is apparent with the displacement vector being characterized by the divergence of a tensor, , a vector , the curl of a vector, and the grad of a scalar . That is, all classes of variable up to and including second-order tensors are represented in the displacement vector. This may be a little spurious but it allows various classes of variables to be brought in as explicit solutions to the Beltrami-Schaeffer functions.

EXPLICIT REPRESENTATIONS OF THE BELTRAMI-SCHAEFFER FUNCTIONS

Helmholtz Functions

The application of the Helmholtz theorem resolves the displacement vector into an irrotational and a solenoidal component for a displacement field which vanishes at infinity. That is

Accordingly, the displacement vector may be given by

where

is the irrotational component of the displacement vector (, a scalar function), and
is the solenoidal component (, a vector function, the divergence of which is zero, )

The above functions come with a superscript since in the algebra below it will be convenient to scale them with constant material properties such that the above becomes

An elegant proof of Helmholtz theorem for homogeneous linear elastic continua may be demonstrated by considering the stress tensor in terms of an arbitrary admissible tensor function. Let

noting that the coefficient in is chosen in order to obtain the prescribed Helmholtz representation above for . The trace of becomes

Continuing, from the constitutive relation for strain

we get the useful expression

the trace of which becomes

and importantly from which integration yields

where is the vector of integration

We suppose that and also that as before, then we get

Thus the desired outcome is achieved, that is, that the Helmholtz characterization of the displacement vector is confirmed through the linear elastic constitutive equation coupled with the integration of the strain tensor. Again, there is good reason for introducing the scaled functions for and . It happens that the Galerkin vector and other functions below yield quite transparently to unification with this format.

It remains to investigate the Helmholtz functions a little further. We can determine a set of independent and uncoupled biharmonic functions for and by firstly setting up the underlying governing tensorial equation for . Firstly, the Helmholtz displacement vector (in terms of and ) are substituted into the stress-displacement expression

to yield

The divergence of this equation leads us to the equilibrium condition , in terms of the Helmholz functions, that is

Applying the divergence operator and to this expression yields the governing biharmonic equation for . Applying the curl operator independently also demonstrates the biharmonic form of .

    and     

The latter arises by applying the Helmholtz condition . So to be clear, the governing biharmonic equations arise from the equilibrium condition. That is, substitution of the Helmholtz displacement functions into the displacement-equilibrium equation yields the defining equations for and . Since the compatibility condition is identically satisfied by the substitution of the displacement vector for the strain tensor, by extension, all displacement functions identically satisfy compatibility.

Extending the Helmholtz Functions

Now for something quite interesting. Consider the following expansion for an arbitrary symmetric tensor

where, as a reminder,

Applying the same relationship as before we can get an expression for as follows

The latter form is that of the Beltrami-Schaeffer format and it is clear that

    and    

whereupon since then

and the displacement vector in the Beltrami-Schaeffer format, becomes

It is a straitforward assumption that and to yield

as before.

A further interesting observation is that after some manipulation the Helmholtz stress tensor may be written directly from the substitution of the Helmholtz functions into the stress-displacement equation as

or more fully

since component vanishes under the curl curl operation.

As before the Beltrami-Schaeffer format compares to the Helmholtz stress equation above and thus

    and    

Substitution for and into the Beltrami-Schaeffer displacement equation is interesting, giving us

where substitution for the Helmholtz equivalences above yields

If       then    

again, as before.

It is worth noting here that is asymmetric, however the operator terms on the right hand side of the stress equation result in symmetry for the stress tensor as you would expect.

Finally and most interestingly we can eliminate the need for the constant vector of integration in the process of determining a Helmholtz basis for the field quantities.

From above we determined a Beltrami-Schaeffer format of the Helmholtz functions to be

Let us revise this expression such that we introduce a new constant multipier on , giving

    and    

where is a special version of for this and only this configuration of the Helmholtz stress tensor.

Thus

Picking up from the equations above, substitution for and into the Beltrami-Schaeffer displacement equation is again

whereupon substitution for the Helmholtz functions yields

If is carefully chosen, in this case such that then the equation for succumbs to a simple substitution such that the vector of integration is dispensed with, giving for the fourth time

This leads to a relatively important characterization of all of the field quantities of linear elasticity in terms of the Helmholtz functions in the Beltrami-Schaeffer format, which also removes any awkward vectors of integration from the formulations. In summary

       and      

where

       and        

Since both the Beltrami-Schaeffer and Helmholtz functions are fundamental to characterizing the field equations of homogeneous linear elasticity it will be interesting to see how the more explicit functions of Galerkin and others unfold in this format.

Papkovich-Neuber Functions

It is known that the Papkovich_Neuber functions may be simply derived from the Helmholtz functions. Firstly with reference to the Helmholtz stress tensor above

Introducing a new vector function such that

where is harmonic since from above .

then taking the divergence of we get

To proceed we define a position vector

Employing the identity the above can be expressed in the form

A second function is now introduced such that

leading to the harmonic nature of both of the Papkovich-Neuber functions, that is

In reality and are just simple substitutes for the underlying combination of Helmholtz functions.

For completeness the Helmholtz functions become

, and

and thus from the Helmholtz displacement vector , an expression for Papkovich-Neuber displacement vector ensues, namely

From the stress-displacement relation the stress tensor directly becomes

which can be manipulated into the format

which is unsurprising since

Ergo, by observation the Beltrami-Schaeffer stress functions can also encode the Papkovich-Neuber functions by putting

    and    

Thus, along with the Helmholtz functions the Papkovich-Neuber functions are simple variants of the same generalized Beltrami-Schaeffer stress functions.

The Galerkin Vector

The Galerkin vector can also be shown to be a variant of the Beltrami-Schaeffer functions.

From above the Helmholtz displacement vector is

Putting and then the known Galerkin form, or a variant thereof, of the displacement vector ensues

Remembering that

then the stress tensor in the Galerkin form becomes

which upon taking the divergence results in the governing biharmonic equation for as before,

The Beltrami-Schaeffer functions thus take on the Galerkin form

    and    

In the alternative mode using the 'Helmholtz tensor'

which results in the equivalent Galerkin stress function format

A very unique expression can now be found for the Galerkin stress tensor since a null deformator field can be added to the first term and the second can be rationalized as follows

which further becomes the interesting expression

or the most compact form of the Galerkin representation for stress

where, not to be confused with , Kirchoff's modulus

Clearly,     and    ,   whilst the elemental function within , that is , remains biharmonic as before.

The associated displacement vector takes the form

As a reminder, Navier's equation is similar

The resulting substitution is not very helpful unless is brought into the equation via however.

Generalised Cylinder Stress Function

A new stress function is currently being researched which has a particular property of being applicable to uni-directional non-homogeneous continua. That is for an elastic field where G, Kirchoff's modulus, is functionally defined along one axis. This will be added to the article as the work is completed.

References

edit
  1. ^ Galerkin, B. (1930). Contribution a la solution generale du probleme de la theorie de l’elasticite’ dans le cas de trois dimensions. Compt. Rend., vol. 190, pp.1047-1048.
  2. ^ Beltrami, E (1892). Observaxionii sulla Nota precedente. Alt. Accad. Nazl. Lincei., vol. pp. 141-142.
  3. ^ Schaeffer, H. (1953). Die Spannungsfunktionen des dreidimensionalen Kontinuums und des elastischen Korpers. Z. angew Math. Mech., vol. 33, pp. 356-362.
  4. ^ Maxwell, J. (1870). On the reciprocal figures, frames and diagrams of forces. Trans. Roy. Soc. Edinburgh, vol.26 pp. 1-40.
  5. ^ Morera, G. (1892). Soluzione generale delle equazioni indefinite dell'equilibrio di un corpo continuo. Atti Accad. Nazl. Lincei. vol 1, pp. 137-141
  6. ^ Papkovich, P.F. (1932) Solution generale des equations differentielles fondamentals d'elasticite exprimee par troois functions harmoniques. Comp. Rend., vol 195, pp.513-515.
  7. ^ Neuber, H. (1934). Ein neuer Anzatz zur Losung raumlicher Probleme de Elastizitatsheorie. Z.angew Math. Mech., vol. 14, pp. 203-212.
  8. ^ Morse, P.M. and Feshbach, H. (1953). Methods of Theoretical Physics, McGraw-Hill, New York.
  9. ^ Navier, C.L.M.H. (1827). Memoire sur les lois de l'equilibre et du movement des corps solides elastiques. Mem. Acad. Sci. Inst. France. vol. 7, pp.375-393
  10. ^ Michell, J.H. (1899). The uniform torsion and flexure of incomplete tores, with application to helical springs,Proc. Lond. Math. Soc. 31, 130-146.
  11. ^ Nadeau G (1964), Introduction to Elasticity, Holt, Rinehart and Winston, New York.
  12. ^ Westergaard, H.M. (1952). Theory of Elasticity and Plasticity, Wiley, New York.
  13. ^ Airy, G. (1863). On the strains in the interior of beams. Phil. Trans. Roy. Soc. London, vol. 153, pp 49-80.
  14. ^ Little, R.W. (1973). Elasticity. Prentice-Hall, New Jersey.

Additional reading

edit
  • Morse, P.M. and Feshbach, H. (1953). Methods of Theoretical Physics. McGraw-Hill, New York.
  • Youngdahl, C.K. (1969). On the completeness of a set of stress functions appropriate to the solution of elasticity problems in general cylindrical coordinates. Int. J. Enging. Sci., vol. 7, pp. 61-79.
  • Gurtin, M.E. (1972). The linear theory of elasticity. S. Flugge, ed., Handbuch der Physik, col. VI1/2, Springer-Verlag, Berlin.
  • Lure, A.I. (1964). Three-Dimensional Problems of the Theory of Elasticity. Interscience, New York.
  • Beasley, A.J. (1986). A New Pair of Galerkin Based Stress Functions in Homogeneous Linear Elasticity. Univ. of Tas., Dept of Civ and Mech. Engng, Research Report No CM 86/4.
  • Beasley, A.J. (1995). An interim unified model of stress functions for homogeneous linear elasticity, Proc. Fourteenth Australasian Conference on the Mechanics of Structures and Materials, Hobart.
  • Sternberg, E. (1960). On some recent developments in the linear theory of elasticity, Structural Mechanics, Pergamon, New York.
  • Sternberg, E., Eubanks, R.A. and Sadowsky, W.A. (1951). On the stress-function approaches of Boussinesq and Timpe to the axisymmetric problem of elasticity theory, J. App. Phys. vol. 22, pp 1121-1124.
edit