Consistency and Validity of the Mathematical Models and the Solution Methods for BVPs and IVPs Based on Energy Methods and Principle of Virtual Work for Homogeneous Isotropic and Non-Homogeneous Non-Isotropic Solid Continua

Energy methods and the principle of virtual work are commonly used for obtaining solutions of boundary value problems (BVPs) and initial value problems (IVPs) associated with homogeneous, isotropic and non-homogeneous, non-isotropic matter without using (or in the absence of) the mathematical models of the BVPs and the IVPs. These methods are also used for deriving mathematical models for BVPs and IVPs associated with isotropic, homogeneous as well as non-homogeneous, non-isotropic continuous matter. In energy methods when applied to IVPs, one constructs energy functional (I) consisting of kinetic energy, strain energy and the potential energy of loads. The first variation of this energy functional (δI) set to zero is a necessary condition for an extremum of I. In this approach one could use δI = 0 directly in constructing computational processes such as the finite element method or could derive Euler’s equations (differential or partial differential equations) from δI = 0, which is also satisfied by a solution obtained from δI = 0. The Euler’s equations obtained from δI = 0 indeed are the mathematical model associated with the energy functional I. In case of BVPs we follow the same approach except in this case, the energy functional I consists of strain energy and the potential energy of loads. In using the principle of virtual work for BVPs and the IVPs, we can also accomplish the same as described above using energy methods. In this paper we investigate consistency and validity of the mathematical models for isotropic, homogeneous and non-isotropic, non-homogeneous continuous matter for BVPs that are derived using energy functional consisting of strain energy and the potential energy of loads. Similar investigation is also presented for IVPs using energy functional consisting of kinetic energy, strain energy and the potential energy of loads. The computational approaches for BVPs and the IVPs designed using energy functional and principle of virtual work, their consistency and validity are also investigated. Classical continuum mechanics (CCM) principles i.e. conservation and balance laws of CCM with consistent constitutive theories and the elements of calculus of variations are employed in the investigations presented in this paper.


Literature Review
Even though this paper does not consider specific aspects of formulation and solution of composite mechanics problems, but such materials are non-homogeneous and non-isotropic, and the energy methods are almost exclusively employed in the study of BVPs and IVPs associated with such materials. It is for this reason that we present a review of published work related to composites based on either energy methods or principle of virtual work. Widespread use of synthesized composite materials is incentive to examine the consistency and legitimacy of the mathematical models for composites as well as the techniques of obtaining their solutions. Broadly, we can categorize composite material in two groups: 1) those that consist of a matrix with impregnated particles of a different material. Ceramic composites with dispersed metal inclusions fall into this category; 2) those that are made by stacking laminas. Each lamina consisting of a matrix (resin) may have short fibers dispersed in the matrix or may contain continuous fibers oriented in specific directions in the plane of each lamina. Laminas are glued and cured with the same epoxy as the matrix of the laminas. Upon curing, the interfaces between the laminas are assumed to have the same strength and properties as the matrix of the laminas. Regardless of which composite we consider; the composites are obviously non-homogeneous and non-isotropic. This is the major source of difficulty in deriving the mathematical models for composites that describe their deformation physics. In the following we present a brief literature review consisting of the development of mathematical models using energy methods and principle of virtual work and their numerical solutions primarily using finite element method.
The classical plate theory is based on Kirchhoff hypothesis (kinematic as-sumptions) for the deformation of the entire laminated plate cross section [1] [2]. Some analytical solutions of the laminated plate problems are represented based on series solutions in references [3] [4] [5] [6]. It is well known that due to pronounced shear deformation effects in composites, the theories based on Kirchhoff hypothesis do not work well. The laminated composite models incorporating transverse shear deformation effects either by using an assumed stress field [7] [8] [9] or by considering an assumed displacement field [10] [11] [12] have been reported. The extension of the displacement based mathematical models using higher order functions in thickness coordinate have been reported in references [13]- [23]. A summary of many of these works can be found in reference [24]. Since the advent of p-version of finite element method (earlier works can be found in references [24]- [35]), its application for solving laminated composite problems grew. More recently the p-version hierarchical finite element method has been applied to plates and shells as well as beams including laminated composites in references [36] [37] [38] [39]. In these works, hierarchical p-approximation for displacements is used in the direction transverse to the plane of the laminate. We remark that the mathematical models currently used for composite mechanics are: 1) Mostly derived based on energy methods or principle of virtual work.
2) Only valid for reversible mechanical deformation.
In a recent paper Surana et al. [40] showed that all currently used beam mathematical models for isotropic and homogeneous matter are thermodynamically inconsistent. That is, these mathematical models can neither be derived using the conservation and balance laws of classical continuum mechanics (CCM) nor non-classical continuum mechanics (NCCM). The authors pointed out that the root cause of thermodynamic inconsistency is the use of a priori kinematic assumptions. The authors also presented a kinematic assumptions free formulation for bending of beams that is accurate for slender as well as deep beams and is free of shear correction factors that plague all currently used mathematical models. In another recent paper Surana et al. [41] also showed that all plate/shell mathematical models for isotropic homogeneous matter used currently that are based on kinematic assumption(s) and derived using energy methods also are thermodynamically inconsistent. The authors also presented a thermodynamically consistent plate/shell formulation with complete 3D deformation physics. This new formulation is free of kinematic assumption(s), is accurate for very thin as well as thick plate/shells and does not require any shear correction factors. Thus, we see that introduction of a priori kinematic assumptions followed by the derivation of the beam, plate/shell mathematical models (generally using energy methods) result in thermodynamically inconsistent mathematical models. We remark that in the work presented here, the objective is not to embark on developing better or more sophisticated methods for ob- taining the solutions of the mathematical models for isotropic and non-isotropic matter or development of new approaches for deriving mathematical models based on energy methods. This work focuses on determining consistency and legitimacy of the mathematical models themselves for isotropic as well as non-isotropic matter that are derived based on energy methods and principle of virtual work for BVPs and IVPs as well as the solution methods based on these two approaches.

Scope of Present Work
Based on the principles of classical continuum mechanics (CCM) i.e. the conservation and balance laws, consistent derivation of constitutive theories based on the entropy inequality and elements of the calculus of variations, we establish some guidelines that are helpful in the investigations presented in this paper.
1) Mathematical models consisting of the conservation and balance laws of CCM and non-classical continuum mechanics (NCCM) with consistent constitutive theories based on the entropy inequality ensure thermodynamic equilibrium during the entire evolution. These mathematical models are meritorious, hence preferred if these can be derived for the deformation physics under consideration. This would be the preferred choice for the mathematical models for composites.
2) Calculus of variations considers extremum of functionals and provides a consistent and organized mathematical framework for obtaining solutions of the mathematical models of the initial value problems (IVPs) and boundary value problems (BVPs) regardless of how these models are derived and regardless of the nature of the differential operators in these models [42] [43].
3) We first consider elements of calculus of variations as applied to BVPs. We construct a functional I corresponding to the BVP or corresponding to the physics of the BVP. If I is differentiable in its arguments, then δI is unique and δI = 0 is a necessary condition for an extremum of I. If 2 I δ provides a unique extremum principle (sufficient condition, 0, 0, 0 > = < ), a solution obtained from δI = 0 is unique, hence yields a unique extremum of I. The solution obtained from δI = 0 also satisfies Euler's equation, a differential or partial differential equation(s) extracted from δI = 0 using fundamental lemma or the fourth lemma. Thus, a correspondence between the extremums of functional and solution of BVP (Euler's equation) is established. In structural mechanics for isotropic as well as non-isotropic matter one constructs an energy functional I consisting of strain energy and potential energy of loads. δI = 0 is a necessary condition for an extremum I. At this stage we could use δI = 0: 1) to construct a method of approximation such as finite element method to find an approximate solution that yields extremum of energy functional I and/or 2) we could derive Euler's equation from δI = 0 using fundamental lemma or fourth lemma. This would be the differential form of the mathematical model associated with the energy functional I. In this paper we investigate whether this approach of using an energy functional consisting of strain energy and potential energy to derive the differential form of mathematical model for BVPs is consistent and valid way to derive differential mathematical models for homogeneous, isotropic as well as for non-homogeneous, non-isotropic matter. We also investigate the approach based on principle of virtual work.
4) Next we consider elements of calculus of variations as applied to IVPs. In this case we construct a space-time functional I corresponding to the IVP or corresponding to the physics of the IVP. If I is differentiable in its arguments, then δI is unique and δI = 0 is a necessary condition for an extremum of I. If 2 I δ provides a unique extremum principle (sufficient condition, 0, 0, 0 > = < ), then a solution obtained from δI = 0 also satisfies Euler's equation, a partial differential equation in space and time extracted from δI = 0 using fundamental lemma or the fourth lemma. Thus, we have a correspondence between the extremums of space-time functionals and the solutions of the IVPs (Euler's equation in space and time). In structural mechanics for isotropic, homogeneous as well as for non-isotropic, non-homogeneous matter one constructs energy functional I consisting of the kinetic energy, strain energy and potential energy of loads. δI = 0 is a necessary condition for an extremum of I. At this stage (as in the case of BVPs) we could use δI = 0: 1) to construct a method of approximation such as space-time finite element method to find an approximate solution that yields extremum of the energy functional and/or 2) we could derive Euler's equation from δI = 0 using fundamental or fourth lemma. This would be the differential form of the mathematical model associated with the energy functional of I. In this paper we investigate whether this approach of deriving mathematical models for IVPs in differential form from the energy functional consisting of kinetic energy, strain energy and potential energy of loads yields a consistent and valid mathematical model in the differential form for homogeneous, isotropic as well as for non-homogeneous, non-isotropic matter. Similar investigation is also presented for principle of virtual work. 5) Investigations parallel to 3, and 4 are also presented using principle of virtual work.
6) Homogenization processes, their validity and their impact on energy methods and the mathematical models derived from energy methods are also equally important aspects to consider.

Consistency of the Methodologies of Deriving Mathematical Models and the Consistency of the Resulting Mathematical Models
We investigate the consistency and validity of following three different approaches of deriving mathematical models, as well as the consistency and validity of resulting mathematical models for initial value problems (IVPs) and boundary value problems (BVPs) for homogeneous, isotropic matter and non-isotropic, non-homogeneous matter.
1) Methodologies based on energy methods.
2) Methods based on conservation and balance laws of CCM and NCCM.
3) Methods based on principle of virtual work. It can be easily shown that methods under 1) and 2) are a subset of calculus of variations, studies in applied mathematics for obtaining extremums of functionals as well as solutions of IVPs and BVPs [42] [43]. Thus, it is perhaps fitting to discuss methods 1) and 2) within the context of calculus of variations, a much broader area of study.

Calculus of Variations, Residual and Energy Functionals: BVPs
To be able to apply elements of calculus of variations to all BVPs rigorously and uniformly for obtaining their solutions or to be able to derive mathematical description of the associated boundary value problem, it is essential to mathematically classify all differential operators appearing in the BVPs into three possible categories: self-adjoint operators (linear and symmetric), non-self-adjoint operators (linear but not symmetric), and non-linear differential operators (neither linear nor symmetric). This mathematical classification allows us to treat all BVPs rigorously regardless of their origin or the field of application. We can use the elements of calculus of variations: 1) to obtain the solutions of BVPs if the mathematical model is known or 2) to derive the mathematical model if the energy functional associated with the BVP of interest is known. Traditionally the calculus of variations begins with 1) Existence of a functional I (generally by construction) 2) If I is differentiable in its arguments then δI is unique and δI = 0 is a necessary condition for an extremum of I.
A solution obtained from δI = 0 also satisfies Euler's equation (PDE(s)) extracted from δI = 0 using fundamental lemma or the fourth lemma. We consider the boundary value problem (1) in which the differential operator A is linear.
with some boundary condition, closed boundary of open domain x Ω .

Residual Functional
The functional I can be the residual functional constructed using (1). If n φ is an approximation of φ over x Ω , then is the residual function and the residual functional I is constructed using E followed by δI and 2 I δ .
( )   (5) Thus, the integral form (4) is variationally consistent (VC) i.e., unconditionally stable. This process can used to obtain solutions of BVP (1). This of course requires that we know the mathematical model (1). It is obvious that in this case the Euler's equation obtained δI = 0 (4) is indeed the BVP (1). A solution obtained from δI = 0 also satisfies Euler's equation, hence is the solution of the BVP. One could show that this method based on residual functional can also be used for nonlinear differential operators [42] with unique extremum principle. Thus, this method based on residual functional for obtaining solutions of BVP is always unconditionally stable (variationally consistent, VC) regardless of the type of differential operator. This method requires that we know the mathematical model of the BVP, thus can always be used for isotropic, homogeneous matter for which the mathematical model can be derived using CCM or NCCM. In case of non-homogeneous and non-isotropic matter, we first need to derive the mathematical model before we can proceed with this residual functional approach described here. It is important to note that the residual functional is not energy functional. This approach of obtaining solutions of BVPs may not allow us to study specific physics of interest that may require concepts of stiffness, mass, damping, etc.

Energy Functional
The functional I can also be used to represent the energy functional consisting of strain energy and potential energy of loads. We consider self-adjoint, non-self-adjoint and non-linear operators A in (1). First let us consider the case when the differential operator A is known i.e. the mathematical model is known.
Then, we can construct the integral form of (1) over x Ω using fundamental lemma of the calculus of variations.

( )
, 0 This integral form represents necessary condition from which the solution n φ is determined. v in (6) is the test function. Thus, at this stage (6) represents first variation of some functional I set to zero i.e. we assume that, there exists a functional I such that ( ) , 0 n φ and v belong to appropriate scalar product space (given or known) on * Γ . We have assumed existence of functional I and have constructed δI or integral form using fundamental lemma, but extremum principle still needs to be established. It is well known [42] that using (7), various methods of approximation can be considered.
For this case we cannot construct functional I using functionals ( ) , B ⋅ ⋅ , and ( ) l ⋅ in (14) or in any other way.
2) Petrov Galerkin Method (PGM) or weighted residual method (WRM): In this method we also consider the integral form (10) but 3) Galerkin method with weak form (GM/WF) In 1) and 2) we have assumed the operator A to be linear. In this section we further assume that adjoint * A of A is the same as A. The assumption of the operator A being linear and that * A A = defines differential operators appearing in small deformation, small strain physics in structural mechanics in which the mechanical deformation is reversible. Such operators contain even order derivatives of the dependent variable with respect to spatial coordinates. We consider the following based on fundamental lemma.

( )
, 0 (14) or We transfer half of the differentiation from n φ to v using integration by parts.
in (16) l v is linear. For this case the functional I can be constructed.  3) We need to ensure that construction of energy functional ( ) n I φ using strain energy and potential energy of loads and the extraction of Euler's equation (the mathematical model) holds for isotropic, homogeneous matter as well as for non-homogeneous and non-isotropic matter (such as composites). 4) We have shown that the energy methods employing energy functional can only be used for linear reversible mechanical deformation in which the mathematical models contain self-adjoint differential operators. This implies small deformation, small strain and reversible mechanical work deformation physics only. Thus, the usefulness of energy methods for designing computational processes (such as finite element method) or for establishing the mathematical model using the energy functional is limited to self-adjoint operators. We present illustrative examples in a subsequent section. 5) Based on (4), it is clear that the energy methods cannot be used to derive mathematical models for finite deformation and finite strain physics even for homogeneous isotropic matter as in this case, the differential operators are non-linear, hence the existence of energy functional cannot be proven or its construction is not possible using calculus of variations. 6) Use of energy functional approach for non-homogeneous and non-isotropic matter in deriving differential mathematical models needs to be investigated.

Calculus of Variation, Residual and Energy Functionals: IVPs
IVPs contain mathematical descriptions in which the dependent variables exhibit simultaneous dependence on space and time and the domain of definition of IVP is naturally a space-time domain. The space-time differential operator can be mathematically classified as: non-self-adjoint (linear but

Residual Functional
The functional ( ) can represent the residual functional constructed using IVP.
with some boundary conditions and initial conditions. We consider the space-time differential operator to be linear (hence it is non-self-adjoint).
Thus, the space-time integral form (22) is space-time variationally consistent (STVF) [43]. This process can be used to obtain solutions of IVP (19). However, this requires that we know the mathematical model (19). It is straight forward to see that in this case the Euler's equation obtained from δI = 0 (23) is indeed the IVP (19). STVC space-time integral form (23) ensures unconditionally stable computational process due to unique extremum principle (23). This method requires that we know the mathematical model (19) of the IVP, thus, can be used for isotropic, homogeneous matter for which the mathematical model (19) can be obtained from CCM or NCCM. In case of non-homogeneous, non-isotropic matter we first need to derive mathematical model before we can proceed with the residual functional approach presented here.

Energy Functional
The functional ( ) can also represent energy functional consisting of kinetic energy, strain energy and the potential energy of loads. This of course can be constructed directly over space time domain xt Ω without the knowledge of the mathematical model.
We recall that the space-time differential operators are either non-self-adjoint or non-linear. Thus, methods of approximation such as, space-time Galerkin method (STGM), STPGM, STGM/WF will yield space-time integral forms that are space-time variationally inconsistent i.e. if we consider (based on funda- and if ( ) ( ) is not valid because Thus, for the space-time differential operators, if the mathematical model (19) is known, the energy functional ( ) , n I t φ cannot be constructed due to the fact that the space-time operator is either non-self-adjoint or non-linear. This raises an extremely important question and a serious concern whether the mathematical models derived using energy functional consisting of kinetic energy, strain energy and the potential energy of loads are even valid for isotropic, homogeneous matter. The validity of the energy functional itself and the resulting mathematical model for non-isotropic and non-homogeneous matter also need to be investigated.

CCM and NCCM: Mathematical Models for Homogeneous Isotropic and Non-Homogeneous, Non-Isotropic Continua (Composites)
The conservation and balance laws of classical continuum mechanics (CCM) in  NCCM permits inclusion of additional physics in the derivation of the conservation and balance laws, but the assumption that the matter is homogeneous and isotropic must hold in order to be able to obtain the balance laws in the differential forms. Derivation of constitutive theory follows from entropy inequality, still based on the assumptions that the matter is isotropic and homogeneous.
Remarks 1) If we examine the conservation and balance laws in the absence of constitutive theory, we can conclude these are independent of the constitutive theories.
The stress tensor σ and heat vector q are assumed to be present in the deforming continua as a consequence of kinematics of deformation and thermal disturbances applied to the volume of matter. This in fact is true.
2) As explained immediately preceding the remarks, that the assumption of homogeneity and isotropy are intrinsic requirements in deriving the differential form of the conservation and balance laws.
3) These assumptions in 2) must hold regardless of the constitutive theories.

4) Based on 1)-3), we have additional restrictions on constitutive relations.
Regardless of whether we are to substitute the constitutive equations in balance laws or if we are to consider them as additional equations augmenting the conservation and balance laws, the constitutive theories must be for the isotropic and homogeneous matter as the conservation and balance laws only hold for isotropic and homogeneous continua. In other words, anisotropic or orthotropic constitutive relations can neither be included as additional equations in the conservation and balance laws of CCM or NCCM nor can these be substituted in them.
In view of the remarks, the phenomenological orthotropic constitutive theories such as in Section 3.1, though possibly verifiable experimentally, but cannot be used in the conservation and balance laws of CCM or NCCM. More specifically, in the stress strain relationship used for non-isotropic matter such as composites, the stresses as a function of strain cannot be substituted in the BLM and FLT (or considered as additional equations) as these balance laws are only valid for isotropic, homogeneous matter.
Thus, it is straightforward to conclude that the conservation and balance laws of CCM (or NCCM) cannot be used as a possible mathematical model for non-isotropic and non-homogeneous matter in conjunction with phenomenological constitutive theory. At this stage, it is worth mentioning that in many works the differential form of BLM (valid for isotropic and homogeneous matter) is used in the process of determining solutions of model problems in composite mechanics using orthotropic constitutive theories. Based on the details presented here, this is not valid and the resulting solutions are not correct solutions of the actual model problems.

Material Properties: Non-Homogeneous, Non-Isotropic, Orthotropic (Laminated Composites)
The process of constructing laminas with preferred orientation of the fibers at the centerline of each lamina and the stacking of lamina in a desired sequence of material direction orientations by gluing them with epoxy (same resin as the matrix of the laminas) and the curing to construct a laminate is a well-known procedure. In this type of construction, each lamina contains fibers assumed to be at the center of the lamina with different material properties in two perpendicular directions in the plane of the lamina. Matrix (i.e. resin) is assumed to be isotropic and homogeneous with mechanical properties significantly different from the fibers. However, at the common boundary and in its vicinity between the glued laminas, the material is essentially resin (or matrix) which is assumed to be isotropic and homogeneous. Thus, our view is that the material properties make a transition from the center of the lamina (properties close to fiber) to the transverse boundaries of the lamina (properties close to resin or matrix).
In the published works based on rules of mixtures and volume fractions and other homogenization considerations, the mechanical properties of the fibers and the matrix are used to obtain equivalent homogenized mechanical properties that are called orthotropic in the two orthogonal directions in the plane of the lamina and the third one orthogonal to the lamina plane. These mechanical properties are assumed to hold throughout the lamina. When laminas with these homogenized properties are used to construct a laminate, naturally some discontinuity of stress is realized at the lamina interfaces if a single common displacement approximation is used for the entire laminate as the case is in classical plate theory (CPT), first-order shear deformation theory (FSDT), higher-order shear deformation theory (HSDT) etc. The homogenization process is somewhat empirical and phenomenological and the assumption of the same homogenized properties to be valid throughout the lamina thickness seems not in accordance with what the construction of the laminate suggests. From the point of view of the work presented here, it is sufficient for us to note that from the actual non-homogeneous, non-isotropic matter we obtain homogenized material properties that still remain non-isotropic.

Constitutive Equations for Orthotropic Materials (Laminated Composites)
CCM and NCCM only provide a thermodynamically consistent framework for deriving constitutive theories for homogeneous and isotropic continua, solid or fluent. For small deformation, small strain, isotropic and homogeneous solid continua with reversible mechanical deformation one could show [44] that the in which ( ) Φ ε is Helmholtz free energy density and ( ) π ε is strain energy density. The frame invariance requirement of the material coefficients in the constitutive theory necessitates that instead of ( ) Φ = Φ ε and ( ) π π = ε we must consider Φ and π functions of the invariants of ε , either principal invariants , , [44]. Since the two sets of invariants are related, we can choose either. Thus, we consider Details of the derivation can be found in reference [44]. A linear constitutive theory resulting from (29) or (30) will contain only two material coefficients: Lame's constants µ and λ or modulus of elasticity E and Poisson's ratio ν .
A constitutive theory with more than two material coefficients will be non-linear constitutive theory in ε and its invariants.
Currently used constitutive theories for anisotropic and orthotropic materials are derived by considering Taylor series expansion of ( ) Φ ε or ( ) π ε in ε about a known configuration and then substituting it in (27)

Remarks
1) The constitutive theory (32) is purely phenomenological. It is calibrated by conducting experiments to determine nine material coefficients.
2) A serious drawback of assuming ( ) Φ = Φ ε and ( ) π π = ε is that material coefficients are now functions of ε in a known configuration, hence are not frame invariant, a strict requirement for the material coefficients. However, if the material coefficients are not dependent on deformation (i.e. are constant), then this is not a handicap, but certainly limits them to be constant throughout the evolution.
3) When the arguments of Φ and π are invariant of ε , the resulting linear constitutive theory for homogeneous, isotropic solid continua requires only two material coefficients. Clearly the derivation based on Taylor series expansion of ( ) Φ ε and ( ) π ε in ε about a known configuration is not supported by continuum mechanics (CCM or NCCM). Thus, this constitutive theory for orthotropic matter is purely phenomenological or empirical. 4) We also remark that the orthotropic (or anisotropic in general) constitutive relations can neither be substituted in the balance laws nor the balance laws can be augmented using these constitutive relations as additional equations due to the fact that differential form of the balance laws only holds for homogeneous and isotropic matter.

Mathematical Models for BVPs Using Energy Functional
We consider methods of deriving mathematical models based on calculus of variations, more specifically energy method for BVPs for homogeneous, isotropic matter as well as for non-homogeneous, non-isotropic matter (such as composites). We recall that the orthotropic constitutive relations used for composites (in particular laminated composites) are also frequently used in other applications as well. As shown earlier, these constitutive relations are neither supported by CCM nor NCCM as these only consider isotropic homogeneous matter.
In this approach of deriving mathematical models for BVPs considered in this section, one constructs an energy functional I consisting of strain energy and potential energy of loads over the spatial domain x Ω of the deforming volume of matter. We consider small deformation and small strain and furthermore we consider the matter to be isotropic. For this case, we can define the energy functional I given in the following.
in which Φ is the Helmholtz free energy density, η is the entropy density, θ in which Thus using (41)  We transfer one order of differentiation from i u to ji σ in the first term of the integrand in (42).
If I is differentiable in its arguments, then I δ is unique and 0 I δ = is necessary condition for an extremum of I.
For simplicity we assume that using fundamental lemma, we obtain the following Euler's equation from (45) We remark that since the volume x Ω contains homogeneous and isotropic matter, each material particle in x Ω is identical. Thus, size of the volume x Ω is irrelevant, hence permitting us to extract (46) from (45) at a material point. When the matter is non-isotropic and non-homogeneous, the first problem is with the strain energy term in (39) which is only valid for isotropic and homogeneous matter. Even if we could find a more justifiable expression for strain energy for non-isotropic matter, Euler's equation similar to (46) cannot be extracted from δI = 0 due to the fact that volume Ω contains non-isotropic matter, hence each material point could potentially be different in volume x Ω . This forces us to maintain the volume integral form (45) instead of differential form (46) so that we know the specific material point in it. We note that (46) is balance of linear momenta for homogeneous, isotropic matter (for BVPs). Thus, when the deforming matter is isotropic and the mechanical deformation is reversible, the energy functional for BVPs only contains strain energy and potential energy of load in which case the Euler's equations derived from the energy functional are indeed balance of linear momenta as expected. In this case, the differential operator (in displacements) is linear and its adjoint is the same as the operator itself. An integral from constructed using (46) based on fundamental lemma and GM/WF will yield a weak form CCM yield the same mathematical model. This is only possible for a linear system in which the mechanical deformation is reversible. Strains and deformation are small and the matter is isotropic and homogeneous.

Mathematical Models for IVPs Using Energy Functional
The space-time differential operators appearing in IVPs are either non-self-adjoint or non-linear [43]. We can use STGM, STPGM, STWRM and STGM/WF for constructing space-time integral forms based on the fundamental lemma of the calculus of variations for space-time functionals (Surana et al. [43]). The integral form resulting from any of the methods for both types of operators are space-time variationally inconsistent. That is, corresponding to these integral forms, representing first variation (δI) of some functional I, a unique extremum principle (a sufficient condition) does not exist. In this case, the resulting computations based on the integral form (δI = 0 for some functional I) are not ensured to be unconditionally stable. The bilinear functional (for linear space-time differential operators) in the space-time integral form is not symmetric. Hence, based on the calculus of variations [42] [43], it is not possible to construct a space-time functional using the functionals appearing in the space-time integral form. That is, corresponding to a valid IVP, a valid energy functional cannot be constructed. Conversely, there does not exist a valid energy functional such that the Euler's equation extracted from its first variation set to zero represents a valid IVP. In the present work we show that if we construct an energy functional consisting of kinetic energy, strain energy and potential energy of loads, then the Euler's equations extracted from its first variation of this functional set to zero cannot yield a valid mathematical model. This is obviously due to the fact that for a valid mathematical model for an IVP, it is not possible to construct a corresponding valid space-time functional due to the fact that for space-time linear differential operators * A A ≠ . Thus, even if we know the PDEs describing IVP, then the functional I corresponding to these cannot be constructed. On the other hand, if we construct energy functional I directly without using space-time differential operator, the Euler's equations resulting from the first variation of this functional cannot be the valid mathematical model (PDEs). First, we consider a simple IVP to illustrate the main source of problem in extracting mathematical models from energy functional for IVPs.

Model Problem
Consider the following IVP representing 1D balance of linear momenta in which all coefficients in the dimensionless form are unity. We have considered small deformation, small strain reversible mechanical deformation physics with linear constitutive theory for the Cauchy stress tensor and the matter is isotropic and homogeneous.
( ) ( ) in which the space-time differential operator A is given by be the scalar product space of functions admissible in (49).

Linearity of A:
Clearly operator A is linear as ( ) ( ) We use BCs (50) and ICs (51) and the following to simplify the integrals over closed boundary Γ in (55) and obtain the following.
in which the concomitant is given by , Au v is given by [43] ( )( ) ( )( ) In (62) , Au v is the concomitant. We can write (61) as Clearly the differential operator * A (adjoint of A) from the right hand side of (63) is not the same as the differential operator A.
Remarks 1) Thus, in the IVP (for isotropic matter) the space-time differential operator is linear but its adjoint * A is not the same as the operator A itself. This is pri-

Integral Form for Model Problem Using Fundamental Lemma and STGM/WF
Let v u δ = , be the test function, then based on fundamental lemma we can write the following using (49) over the space-time domain We transfer one order of differentiation from u to v for each term in the integrand of (65) ( ) Simplification of the boundary integrals (see Equation (57) and (59) which is not the same as (69)

Derivation of Mathematical Model for IVP from the Energy Functional
In this section, we construct an energy functional I consisting of kinetic energy, strain energy and potential energy of loads that is believed to correspond to the initial value problem (49). We consider isotropic and homogeneous matter. Since (49) is in dimensionless form, we can construct an energy functional as follows We take first variation of I and set it to zero, necessary condition for an extremum of I. 11 1 1 d 0 We transfer differentiation from u δ to u t ∂ ∂ and 11 σ using integration by it has been shown using (82) in (81) In the absence of the scalar product term over 4 Γ , the Euler's equation would have been (when the matter is homogeneous and isotropic)

Energy Functional in  3 and the Mathematical Model for IVP
The energy functional I consisting of kinetic energy, strain energy and potential energy of loads in 3 R can be written as (not nondimensionalized) follows for homogeneous, isotropic matter.
we note that (84) which is balance of linear momenta. However, since the integrals over 4 Γ are not zero, Euler's equations (91) cannot be derived from (90). Remarks 1) In structural mechanics for homogeneous, isotropic matter as well as non-homogeneous, non-isotropic matter such as composites, the mathematical model (91) is used and is derived using the energy functional (86). The derivation presented here clearly shows this involves approximation of neglecting the boundary integrals over 4 Γ (open boundary at t τ = ) in (90).
2) Another important point to note is that, in the derivation of the Euler's equations that they are valid at a material point and for every material point in the volume, there is the inherent assumption that the continua are homogeneous and isotropic. Only then the choice of volume is arbitrary as every material point in the volume is identically the same, hence we can set the integrand to zero which gives us the differential model at a material point.
3) We finally conclude that the derivation presented is only possible for isotropic matter and that (91) cannot be derived from the functional I in (86) as used commonly in the published works for both isotropic and non-isotropic solid continua.

Principle of Virtual Work
Principle of virtual work has been used in CCM as well as NCCM. The basic concepts remain the same in both cases. Thus, it suffices to consider CCM only in the following (in Langrangian description). Definition: When a deformed body (volume of matter) is in stable equilibrium due to actions and reactions, then the work done by all forces due to virtual displacements is zero.
The virtual displacements at the material points are hypothetical or imaginary displacements such that due to the application of these, the boundary conditions and the loads as well as their points of application remain unaltered, hence ensuring that the equilibrium of the deformed body is not affected due to virtual displacements.
We consider a solid continuum with deformed volume V bounded by V ∂ .
We consider small deformation, small strain, hence the undefined and defined coordinates as approximately same ( x x ≈ ), thus no change in volume, implying that (V V ≈ ) and ( V V ∂ ≈ ∂ ). Let the volume V be subjected to body forces, inertial forces and tractions P on V ∂ . Then, be virtual displacements at the material point x or x then the total virtual work for the volume V is (no sum over i, but sum over j) is given by the following and is zero.
where i v are test functions, then we can write (93) as we note that which in fact are the balance of linear momenta. Remarks 1) We see that starting with the statement of virtual work (92), we are able to obtain the corresponding differential form of the mathematical model using the fundamental lemma of the calculus of variations.
2) It is obvious that choice of the work terms in (92) is crucial. In this example we considered time dependent deformation of a simple solid for which the inertial forces, body forces and those due to tractions P on V ∂ are well understood and known. This may not be the case in more complex deformation physics.
3) Equivalence of the virtual work statement (92) and the integral form (98) resulting from it using the fundamental lemma is important to note. This confirms that (92)

Summary and Conclusion
In the following we summarize the work presented in this paper and draw some conclusions.
1) The objective of this work is to investigate the consistency and legitimacy of the mathematical models for BVPs and IVPs for homogeneous and isotropic matter as well as for non-homogeneous and non-isotropic matter that are derived using energy functional approach. For BVPs, the energy functional consists of strain energy and the potential energy of loads. For IVPs, the energy functional consists of kinetic energy, strain energy and the potential energy of loads.
2) We consider BVPs in the following: a) For linear, reversible mechanical deformation with small strain and small deformation, calculus of variations, energy methods and CCM, all yield the same mathematical model for isotropic, homogeneous matter. This of course is due to the fact that the differential operators in this case are linear and their adjoint is the same as the operator itself.
b) For non-homogeneous and non-isotropic matter with orthotropic or anisotropic constitutive relations, there is no basis for using ji ij σ ε as strain energy.
We note that work conjugate pair ji ij σ ε is due to entropy inequality containing ji ij σ ε as rate of work conjugate pair. When the matter is not isotropic, there is no basis to assume ji ij σ ε as work conjugate pair. Even if we assume this to be true, the volume x Ω in this case contains non-isotropic matter, hence the differential model (Euler's equation) cannot be derived as we must maintain the volume integral containing the different material points of non-isotropic matter.
Thus, for non-isotropic matter regardless of whether the problem is linear or non-linear (but reversible mechanical deformation), a valid energy functional needs to be established using valid strain energy expression. Even if this is possible to accomplish, the derivation of the differential form of the mathematical model is not possible due to lack of isotropy. and third issue is that, in this case we cannot set the integrand to zero in the 0 I δ = due to the fact that the volume x Ω now contains different material points due to lack of isotropy. d) Thus for IVPs, derivation of a valid differential form of mathematical model from the energy functional consisting of kinetic energy, strain energy and potential energy of loads is not possible for isotropic as well as non-isotropic continua (including composites). This is primarily due to the fact that the space-time differential operators are not self-adjoint (may be linear but * A A ≠ ).
4) Based on work presented in this paper, we can conclude that derivation of mathematical models from the energy methods only remains valid for isotropic, homogeneous matter and that a valid mathematical model using energy methods can only be derived for BVPs with linear, isotropic reversible mechanical work with small strain, small deformation physics.

5)
We also remark that use of the energy functionals directly in processes such as finite element method employing integral forms over element domain e x Ω also requires caution: a) In using ji ij σ ε as strain energy when the matter is non-isotropic as there is no basis for it based on entropy inequality for isotropic matter. Thus, in view of the thermodynamically consistent mathematical models based on CCM or NCCM and calculus of variations, there is virtually no benefit to