Next Article in Journal
Research on and Model Analysis of Flexural Mechanical Properties of Basic Magnesium Sulfate Cement Concrete Beams
Previous Article in Journal
Assessing Microplastics and Nanoparticles in the Surface Seawater of Venice Lagoon—Part I: Methodology of Research
Previous Article in Special Issue
Harnessing a Dielectric/Plasma Photonic Crystal as an Optical Microwave Filter: Role of Defect Layers and External Magnetic Fields
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics in Explicit Gradient Elasticity: Material Frame-Indifference, Boundary Conditions and Consistent Euler–Bernoulli Beam Theory

by
Charalampos Tsakmakis
1,*,
Carsten Broese
1 and
Stergios Alexandros Sideris
2
1
Institute for Mechanics, Civil Engineering, Technical University of Darmstadt, Franziska-Braun-Str. 7, 64287 Darmstadt, Germany
2
Ephorate of Antiquities of Lesvos, 81100 Mytilini, Greece
*
Author to whom correspondence should be addressed.
Materials 2024, 17(8), 1760; https://doi.org/10.3390/ma17081760
Submission received: 1 February 2024 / Revised: 27 March 2024 / Accepted: 3 April 2024 / Published: 11 April 2024

Abstract

:
The paper is concerned with the boundary conditions of explicit gradient elasticity of Mindlin’s type in dynamics. It has been argued in an earlier paper that acceleration terms should not be present in the boundary tractions because of objectivity arguments. This is discussed in the present paper in more detail, and it is supplemented by assuming the validity of the principle of material frame indifference. Furthermore, new examples are discussed in order to illustrate that significant differences exist in the responses predicted by boundary tractions with and without acceleration terms.

1. Introduction

In classical elasticity, the stress tensor at a material point is a function of the strain tensor at that point. Such constitutive theories are called local. If the constitutive relations at a material point account, besides for the values of the state variables at the considered point, also for the values of the state variables in a neighborhood of that point, then the constitutive theory is called nonlocal. A simple possibility to capture nonlocality effects in the material response at a point is to incorporate in the theory the gradients of the state variables at that point. The arising theories are called gradient constitutive theories. Examples of gradient theories in solid mechanics are the gradient elasticity (see the references cited in [1]), the gradient plasticity (see the references cited in [2]) and phase field approaches to fracture (see the references cited in [3]), to name a few. It seems that systematic incorporation of gradient effects in elasticity, the present paper deals with, was initiated by the works of Korteweg in the year 1901 (see [4] and the references cited there) and Cosserat and Cosserat in the year 1909 (see [5] and the references cited there). The latter, accomplished with inertial terms, is nowadays known as micropolar elasticity [5]. The idea in this theory is to extend the notion of the classical continuum by attaching at any (macroscopic) material point a microcontinuum (microstructure), which is allowed to rotate like a rigid body. Whenever the microstructure is postulated to undergo homogenous deformations, the resulting continuum is denoted as micromorphic ([6], p. 5). The micromorphic elasticity introduced by Erringen (see [5], Section 7) and the microstructured elasticity theory introduced by Mindlin [7] are, in essence, the same and represent milestones in the development of nonclassical elasticity theories. Nonlocality effects in these theories are captured by the gradient of the deformation of the microcontinuum. The particular case where the microdeformation is set equal to the macrodeformation has been considered in [7] and is known as Mindlin’s gradient elasticity (the static version of this elasticity has been established by Toupin [8] without reference to a microstructure). Among others, Mindlin’s gradient elasticity can describe so-called length scale effects, which cannot be modeled by classical elasticity. Experimental evidence of such effects may be found in microbending tests of epoxy cantilever beams (see [9]), in vibration tests of nickel cantilever microbeams (see [10]), in pantographic structures (see [11]), in dispersion curves observed in metals, alloys and concrete (see [12]), in human calcaneus bones and in fluid saturated porous materials (cf. the references cited in [13,14]).
Mindlin’s gradient elasticity has been established by using a classical version of Hamilton’s variational principle. It obeys, for dynamic problems, traction boundary conditions, which include acceleration terms. This has been criticized in [15] by applying objectivity arguments. That means, as it has been stated in [15], that the boundary tractions are nonobjective and hence such boundary conditions are inacceptable physically. As a consequence, these authors proposed traction boundary conditions not including acceleration terms. Also, they calculated some uniaxial vibration examples in order to illustrate that the differences in the responses predicted by the theories with and without acceleration terms in the boundary tractions may be significant.
The present paper is concerned with two aspects of the traction boundary conditions in Mindlin’s gradient elasticity, which have not been addressed in [15]. The nonobjectivity of the boundary (contact) forces follows from the transformation rules for a change of observer. (The different types of forces used in classical mechanics and their objectivity properties are sketched in Section 2.2). However, as in the case of inertial forces (see Section 2.2), one could postulate the objectivity of contact forces, even if acceleration terms are present in these forces. This is the first aspect addressed in the paper. It is shown, as one might expect, that if the principle of material frame indifference is assumed for contact forces (and this is assumed in the present work), then to postulate objectivity of contact forces including acceleration terms is equivalent to the violation of this principle. These issues are discussed in Section 3.3. The second aspect addressed in the paper concerns the examples used to demonstrate the differences in predicted responses calculated with the two different types of boundary tractions. Nonvanishing nonclassical boundary tractions have been assumed to apply in the examples calculated in [15]. However, for the time being, it is not yet clear how to realize nonclassical boundary tractions. Therefore, homogenous nonclassical boundary tractions are assumed in the examples of the present paper, which seems to be more realistic. In addition, besides uniaxial loading of a bar, also vibrational loading of consistent Euler–Bernoulli beams, and besides traction-controlled also displacement-controlled loading histories are considered. When formulating boundary conditions, an appropriate way is to employ variational methods. For our purposes, especially, we find it convenient to employ extensions of Hamilton’s variational principle as explained in Section 2.1. (Note that bending of beams for Mindlin’s gradient elasticity in dynamics and within a consistent Euler–Bernoulli framework is studied for the first time in the present work). In summary, the present work completes the analysis provided in [15]. It does not investigate properties of a particular model, but rather provides an examination of the boundary tractions of the whole gradient elasticity of Mindlin’s type.

2. Basic Relations

2.1. Extensions of Hamilton’s Variational Principle

There are various extensions of Hamilton’s principle, which are in common use in solid mechanics. The extensions, which are of interest to our paper, can be discussed in a systematic way in the case of classical elasticity. Then,
δ   t 1 t 2 T Π   i Π   e dt   = 0   ,
is the standard form of the principle, where δ( ) is the variation of ( ) and t is the time. Further, T is the kinetic energy, which for classical elasticity reads
T Τ cl = V   1 2 ρ u ˙ i u ˙ i dV .
Here, ρ is the mass density and V is the range in three-dimensional Euclidean point space, occupied by the material body we consider. The boundary of V is V   and has outward unit normal vector n. Unless explicitly stated otherwise, small deformations are assumed, all tensorial components are referred to a Cartesian coordinate system x i , i = 1 , 2 , 3 , and use is made of the summation convention. Material particles in V are identified with location vectors x = x i e i , where { e i } is the Cartesian coordinate basis. We write   ˙ for the time derivative of ( ) and denote by u i the components of the displacement vector u. The variation δ u x , t is postulated to vanish everywhere at the initial and final times t 1 , t 2 , while δ u vanishes at all times in t 1 , t 2 on the boundary part V u , where displacement boundary conditions are prescribed. Conservative systems are addressed in Equation (1) and Π   i , Π   e are the potential energies of the internal and the external forces, respectively. Omitting volume forces and denoting by φ u , t the potential per unit area for external forces, we have
Π   e = V   φ   dS   δ   t 1 t 2 Π   e   dt   = t 1 t 2 V p   P i δ u i   dS   dt   .
In the integral on the far-right side, P i : = φ u i , V P denotes the part of V where traction boundary conditions are given, and V u V P = V ,   V u V P = . For classical elasticity, there exists a free energy pro unit volume ψ = ψ (ε), so that
Π   i = V   ψ ε   dV   ,
where ε is the strain tensor,
ε ij : = 1 2 i u j + j u i   .
Henceforward, we use the notation i   : =   x i   , x i . Moreover, for classical elasticity, we define the Cauchy stress Σ through
Σ ij : = ψ ε ij   .
Using these relations in Equation (1), after well-known manipulations, we arrive at
t 1 t 2   V   ρ u ¨ j i Σ ij δ u j dV V p   P j n i Σ ij δ u j dS   dt   = 0   .
With the aid of fundamental lemmas of calculus of variations (cf., e.g., [16], Section 2.4.), it can be proved that the balance of linear momentum
i Σ ij = ρ u ¨ j ,
together with the traction boundary conditions
n i Σ ij = P j   on   V P × t 1 , t 2
are necessary and sufficient conditions for Equation (7).
The most frequent extension of the principle (1) is when a part of the external forces is conservative, with potential energy Π   e , and the remaining external forces are not conservative and expend virtual work δ W   e . Then the principle (1) is expressed in the form
δ t 1 t 2 T Π   i Π   e dt + t 1 t 2 δ W   e dt = 0   .
Another extension of the principle (1) arises by introducing, as in d’Alembert’s principle, the inertial force i : =   ρ u ¨ (cf. [17], Section 19), which expends virtual work δ W   inert = V     i j δ u j dV = V     ρ u ¨ j δ u j dV . It can be seen that in this case δ t 1 t 2 T   dt = t 1 t 2 δ W   inert dt . Now, assume that a part of the inertial forces expends virtual work expressible as δ t 1 t 2 T   dt , while the remainder of the inertial forces expend the work t 1 t 2 δ W   inert dt . In this case, the principle (10) is extended as follows
δ t 1 t 2 T Π   i Π   e dt + t 1 t 2 δ W   e +   δ W   inert dt = 0   .
If only δ W   e and only δ W   inert are considered, then the time integration in Equation (11) can be dropped out since
δ W   inert + δ W   e δ Π   i = 0
is the virtual work principle.

2.2. Material Frame-Indifference

Only in this section do the discussions refer to finite deformations. Following Liu [18] and Liu and Sampaio [19] (see also the references cited in these works), we define an observer or frame of reference to be a one-to-one map Φ assigning to a point (event) of space-time, a point on the product space of a three-dimensional Euclidean space and the set of real numbers (time axis). Let x = x X , t be the motion of a material body relative to Φ, where X and x are the location vectors of a material particle in the reference and the actual configuration, respectively. The same motion with respect to another frame Φ* is denoted by x * = x * X * , t and we have
x * = Q t x + c * t   ,
where Q is an orthogonal transformation, c is a relative translation, and t * = t + a , a = const.
Scalars s, vectors v and second-order tensors T are said to be Euclidean objective, or simply objective, if
s * = s ,   v * = Q v ,   T * = Q T Q T   .
In these transformation rules, A T is the transpose of a second-order tensor A and s * , v * , T * are the quantities in the frame Φ* corresponding to s, v, T in the frame Φ, respectively.
In classical mechanics, there are different classifications of forces. On the one hand, forces are divided into those that obey a response law and those that cannot be determined by response laws but are calculated from balance laws and boundary conditions. The existence of the second type of forces arises from imposed geometrical conditions. Examples of the first type of forces are spring forces, obeying an elasticity law, and inertial forces, obeying, e.g., the response law ρ u ¨ with respect to an inertial reference frame. Examples of the second type of forces are reaction forces due to supports, the geometrical constraint being, e.g., no displacement, or the pressure in incompressible continua, the geometrical constraint being isochoric deformations. On the other hand, forces in continuum mechanics are divided into contact, body and inertial forces. A fundamental postulate in classical mechanics is that all forces mentioned above are assumed to be objective.
Now let U be a set of state and kinematical variables and suppose that an objective quantity J obeys a response law J = f Φ U   in the frame Φ. For simplicity, f Φ is supposed here to be a function, but in the general case it may be a functional. Generally, response functions may be observer-dependent, so that J * = f Φ   U * . The response function is said to be observer-independent, or material frame-indifferent (see [6,7,8], Section 19), or simply frame-indifferent, if it has the same form in all frames, i.e., f Φ   = f Φ   . When formulating constitutive laws for ordinary stresses, the response functions are required to be frame-indifferent, which imposes restrictions on the form of the functions that are to be designed. On the contrary, the form of the response functions of inertial forces is supposed to be known with respect to an inertial frame. The forms of the response functions of the inertial forces in other, noninertial frames are derived on the basis of the transformation rules (13) and (14), imposed by the postulated objectivity of inertial forces, and are generally dependent on the observer (cf. [20]). This is a consequence of the unique feature characteristics of inertial forces.

2.3. A Simple Model of Explicit Gradient Elasticity

Consider material bodies which, for omitting body forces, are described by balance laws of linear and angular momentum, so that at any point in V the field equations
j Σ jk + i k = 0   ,
Σ jk = Σ kj  
apply. As in Equation (8), Σ is the Cauchy stress tensor, while i is the inertial force. In Mindlin’s gradient elasticity (see [7]), i is decomposed into classical and nonclassical parts,
i = i cl + i noncl   ,  
with
i k cl : =   ρ   u ¨ k   .
For simplicity,
i k noncl : = γ p p   u ¨ k =   γ   Δ   u ¨ k
is chosen in the present paper, where Δ is the Laplacian operator and γ = const. is a material parameter. The inertial law (19) corresponds to the “isotropic case” in Mindlin’s work (see p. 70 in [7]).
A simple gradient elasticity model arises by assuming the free energy per unit volume ψ to be given by
ψ = ψ ε , k = 1 2 ε jk K jkmn ε mn + 1 2 l 2 k ijk K jkmn   k imn   ,
where
k ijk : = i ε jk
is the gradient of ε. Further, K is an anisotropic fourth-order elasticity tensor exhibiting the symmetries K jkmn = K kjmn = K mnjk and l = const. is an internal material length. According to the gradient elasticity Form II of Mindlin’s theory (see [7]), Σ obeys the constitutive law
Σ jk : = τ jk i μ ijk   ,
τ jk = τ kj : = ψ ε , k ε jk = K jkmn   ε mn   ,  
μ ijk = μ ikj : = ψ ε , k k ijk = l 2 K jkmn i ε mn   ,
where τ is a classical, second-order stress tensor of Cauchy type and μ is a nonclassical, third-order stress tensor. From Equations (22)–(24),
Σ jk = K jkmn ε mn l 2 K jkmn Δ ε mn   .
To our knowledge, this constitutive law, with K being an isotropic elasticity tensor, has been proposed for the first time in Altan and Aifantis [21]. In Broese et al. [22], the gradient elasticity law (25) has been interpreted to represent the gradient elasticity analog of the Kelvin viscoelastic solid. In the following, we shall denote it as Kelvin-Gradient elasticity Model (KG-Model).
To accomplish the theory, concomitant boundary conditions for the field equations remain to be formulated. As mentioned in the introduction, an appropriate way to establish these is to invoke variational principles.

3. Gradient Elasticity in the Setting of Hamilton’s Principle

3.1. Variational Formulation of the Field Equations

After multiplying Equation (15) by the virtual displacement δ u k , integrating over V and using partial integration, we obtain
V   j Σ jk δ u k dV V   Σ jk δ ε jk   dV + V   i k δ u k dV = 0   .
Next, we replace Σ jk in the second integral with the aid of the constitutive law (22),
V   j Σ jk δ u k dV V   τ jk δ ε jk   dV + V   ( m   μ mjk )   δ ε jk   dV + V   i k δ u k dV = 0   ,
or equivalently
V   j Σ jk δ u k + μ jik i δ u k dV V   [ τ jk δ ε jk + μ mjk   δ m ε jk ]   dV + V   i k δ u k dV = 0   .
We recall from Equations (20), (23) and (24), that the second integral is the virtual work of internal forces, i.e.,
V   j Σ jk δ u k + μ jik i δ u k dV δ Π   i + V   i k δ u k dV = 0   ,
with (cf. Equation (4))
Π   i : = V   ψ ε ,   k   dV   .
In order to recast the virtual work expended by the inertial force i, we recall from Equations (17)–(19) that
V   i k δ u k dV   = V     ρ u ¨ k   δ u k   dV   + V   i k noncl   δ u k   dV   = d dt   V   ρ u ˙ k   δ u k   dV   + δ V   1 2 ρ u ˙ k u ˙ k   dV   + V   i k noncl   δ u k   dV = d dt V   ρ u ˙ k   δ u k   dV + δ T cl + V   i k noncl   δ u k   dV ,
where T cl is the classical kinetic energy defined in Equation (2). Thus, by taking the time integral between t 1 and t 2 of Equation (29), and using the result (31), we find that
t 1 t 2 V     j Σ jk δ u k + μ jik i δ u k dVdt δ t 1 t 2 Π   i   dt   + δ t 1 t 2 T cl   dt   + t 1 t 2 V   i k noncl   δ u k   dVdt = 0   .
Now, there are two ways to recast the volume integral of the last term. One is according to Mindlin [7] and leads to boundary conditions involving acceleration terms, whereas the second way, according to Broese et al. [15], leads to boundary conditions in which acceleration terms are not present.

3.2. Gradient Elasticity with Acceleration Terms Present in the Boundary Conditions

The aim of Mindlin [7] was to bring Equation (32) to a form corresponding to Equation (1), or to the more general Equation (10). To achieve this, the last term in Equation (32) is rewritten using the definition (19):
t 1 t 2 V   i k noncl   δ u k   dVdt   = t 1 t 2 V   γ j j u ¨ k δ u k   dVdt   = t 1 t 2 V   j γ j u ¨ k δ u k dVdt   t 1   t 2   V   γ j u ¨ k   δ ( j u k )   dVdt   = t 1 t 2 V   j γ j u ¨ k δ u k dVdt   t 1 t 2 d dt   V   γ j u · k   δ ( j u k )   dV   dt   + t 1 t 2 V   γ   j u · k   δ j u · k   dV   dt   .
On the far-right side, the second integral vanishes since j δ u k t 1 = j δ u k t 2 0 in V, while the last integral can be represented by the variation of nonclassical kinetic energy
T ¯ noncl : = V   1 2 γ   j u · k   j u · k dV   .
Hence, Equation (33) implies
t 1 t 2 V   i k noncl   δ u k   dVdt = t 1 t 2 V   j γ j   u ¨ k δ u k dVdt + δ t 1 t 2 T ¯ noncl   dt   ,
and by inserting into Equation (32), we obtain
δ t 1 t 2 T cl + T ¯ noncl Π   i dt + t 1 t 2 V     j Σ jk + γ j   u ¨ k δ u k + μ jik i δ u k dVdt   = 0   .  
This has been interpreted by Mindlin to suggest defining a total Cauchy stress tensor Σ t involving acceleration terms,
Σ jk t : = Σ jk + γ j u ¨ k   ,
a virtual work of external forces
δ W ¯ e : = V   n j Σ jk t δ u k + μ jik i δ u k dS   ,  
and the total kinetic energy
T ¯ : =   T cl + T ¯ noncl   .
Then, Equation (36) becomes
δ t 1 t 2 T ¯ Π   i dt   + t 1 t 2 δ W ¯ e dt   = 0   ,  
which is of the form (10).
To accomplish Mindlin’s approach, the surface integral δ W ¯ e must be resolved further, since the gradient j δ u k is not independent of δ u k on V . After lengthy and elaborate algebraic manipulations, it can be proved (see [7,15]) that
δ W ¯ e = V   P ¯ k   δ u k + R ¯ k   D δ u k dS   ,
where P ¯ and R ¯ are classical and nonclassical traction vectors, respectively, defined by
P ¯ k : = n j Σ jk t D j n i μ ijk + D l n l n i n j μ ijk = n j Σ jk + γ j u ¨ k n j D j n i μ ijk + D l n l n i n j μ ijk   ,
R ¯ k : = n i n j μ ijk   .  
For a function f(x,t), the normal derivative Df(x,t) and the surface derivative D i f x , t are defined through
Df : = n l l f   ,   D i f : = i f   n i Df   .
This way, δ W ¯ e is interpreted as the virtual work expended by the tractions P ¯ and R ¯ , and Equation (40) is assumed to be the appropriate form of Hamilton’s principle for the considered material. Since now δu and Dδu are independent of each other on V , the adjoint boundary conditions suggested by Equation (41) are
either   P ¯ k   or   u k   and  
either   R ¯ k   or   Du k  
have to be prescribed on V .
Altogether, for the gradient elasticity based on the KG-Model, Equations (15)–(19) are the governing equations of motion and Equations (45) and (46) are the adjoint boundary conditions proposed by Mindlin. The remarkable feature in the traction boundary condition (45) is the presence of acceleration terms in P ¯ (see Equation (42). Moreover, Σ t cannot be a proper classical Cauchy stress tensor, for it is generally not symmetric (because of the presence of the acceleration term in Equation (37)).

3.3. Gradient Elasticity without Acceleration Terms Present in the Boundary Conditions

As mentioned in Section 2.2, contact, body and inertial forces are postulated to be objective and, in addition, constitutive functions of ordinary Cauchy stresses are required to be frame-indifferent. It can be proved, that the presence of acceleration terms in Equation (42) renders the traction vector P to be non objective. If someone would get the idea to postulate P as objective, as in the case of inertial forces, then one would have to conclude that the constitutive function of stress Σ t in Equation (37) will not be frame-indifferent. Either way, we believe that the presence of acceleration terms in boundary conditions for contact forces is physically not acceptable and this is, in principle, the criticism made by Broese et al. [15]. These authors proposed an alternative but not equivalent definition for the virtual work expended by external forces. Considering inertial forces to have the nature of conventional body forces (see [17], Section 21), they proposed to interpret/rewrite Equation (32) as follows.
Let δ W ^ inert , noncl denote the virtual work of the inertial force i noncl ,
δ W ^ inert , noncl : = V   i k noncl δ u k dV = V   γ Δ u ¨ k   δ u k dV ,
and reformulate Equation (32) in the form
  δ t 1 t 2 T cl   Π   i   dt   + t 1 t 2 δ W ^ e + δ W ^ inert , noncl   dt   = 0 ,
where
δ W ^ e : = V   j Σ jk δ u k + μ jik i δ u k   dV .
We view Equation (48) to be of the form (11), and recognize δ W ^ e as the proper virtual work of the external forces and Σ jk = Σ kj as the proper Cauchy stress tensor, which should enter into the boundary tractions. By performing similar calculations as for the transition from Equation (38) to Equations (41)–(43), we can introduce proper boundary tractions P ^ and R ^ , so that
δ W ^ e = V   P ^ k δ u k + R ^ k δ Du k dS ,
where (cf. [15,16])
P ^ k : = n j Σ jk D j n i μ ijk + D l n l n i n j μ ijk ,
R ^ k : = n i n j μ ijk .
Accordingly, the proper boundary conditions read
either   P ^ k   or   u k   and  
either   R ^ k   or   Du k  
have to be prescribed on V . No acceleration terms are now involved in the boundary condition for classical traction, but otherwise, the governing equations of motion are the same as in Equations (15)–(19).

4. Consistent Euler–Bernoulli Beam Theory in Dynamics

4.1. Main Assumptions

The traditional Euler–Bernoulli beam theory relies upon two fundamental assumptions (see [23], Sections 5.1 and 5.4.2 and [13], p. 90): (1) The material response is isotropic elastic. (2) Plane cross sections of the beam remain plane and perpendicular to the deformed beam axis and the shape of the cross sections does not change (no deformation of cross sections). It is obvious, that no deformation of cross sections is not consistent with isotropic elasticity. This inconsistency is reflected by the fact that the elasticity law is satisfied in local form, but, e.g., the equilibrium equations in statics are not satisfied in local form (see [24] and the references cited there).
To overcome this problem, Sideris and Tsakmakis [24] (see also [25]) proposed to drop the assumption of isotropy and instead suppose transversal isotropic material behavior subject to geometrical constraints. This way, they obtained consistent Euler–Bernoulli beam theories. The main kinematical and constitutive equations of this approach for the KG-Model can be summarized as follows (see [24]).
Consider the beam in Figure 1, which is of length L, and constant cross section A, with width 2b and height 2c. The origin of the Cartesian coordinate system { x i } is located on the left boundary plane, the x 1 x 3 plane is a symmetry plane, and the x 1 axis is the centroidal axis of the beam. The beam might be subject to a transverse load, which acts in the x 1 x 3 plane at x 3 = ±   c , and to axial loads, which act on the boundary planes x 1 = 0 and x 1 = L . In addition, problem-specific boundary/reaction forces will apply. The material response is assumed to be transverse isotropic with vanishing in-plane Poisson ratios and subject to internal constraints. All assumptions together cause conditions of plane strain and plane stress with u, ε, Σ being functions only of x 1 , x 3 and
u 2 = ε 2 i = Σ 2 i = 0 .
The remaining components of the displacement vector and the strain tensor are given by
u 1 = U x 1 , t w x 1 , t x 3 ,   u 2 0 ,   u 3 = w x 1 , t ,
ε 11 = U   w x 3 ,   ε 33 = ε 13 = 0 ,
where U is axial displacement, w x 1 , t is the deflection curve and f x 1 , t is the derivative of the function f x 1 , t with respect to x 1 . The assumed geometrical constraints do not allow to determine every stress component from the elasticity laws (22)–(24). In fact, we have (see [24,25])
Σ 11 = E ε 11 μ 111 =   E   U w x 3 l 2   E   U w x 3   ,
τ 11 = E ε 11 = E U w   x 3 ,
μ 111 = l 2 E ε 11 = l 2 E U w   x 3 ,
μ 311 = l 2 E     3   ε 11 = l 2 E   w ,
τ 33 , τ 13 , μ i 33 , Σ 33 , Σ 13 : not   determinable   by   constitutive   law ,
while ψ in Equation (20) becomes
ψ = 1 2 E ε 11 2 + l 2 2 E   1 ε 11 2 + l 2 2 E 3 ε 11 2 .
Every other stress component is vanishing, while E is Young’s modulus in x 1 -direction. Only statics has been addressed in Sideris and Tsakmakis [24] and Broese et al. [25], and therefore only equilibrium equations have come into play in these works. However, the present paper is concerned with dynamics and hence the equations of motion (15)–(19) must be used. For the plane stress and plane strain conditions set up above, these reduce to the two equations
1 Σ 11 + 3 Σ 13 = ρ u ¨ 1   γ u ¨ 1   ,
1 Σ 13 + 3 Σ 33 = ρ u ¨ 3   γ u ¨ 3   .
In the remainder of the paper, we shall use frequently the definitions
f 1 x 1 : = EA U l 2 U ,   f 2 x 1 : = l 2 EA U ,
f 3 x 1 : = EI + l 2 EA w + l 2 EI w   ,   f 4 x 1 : = l 2 EI w .
Furthermore, we denote by I the moment of inertia and we notice that, with respect to the chosen coordinate system, we have
I = A   x 3 2 dS   ,   A   x 3 dS = 0 .
When solving equations of the form (63) and (64), the effort can be reduced considerably by simplifying the problem using section and resultant forces. The resulting equations correspond to a one-dimensional continuum, which is bounded by the points x 1 = 0 , x 1 = L . Points x 1 0 , L are interior points and traction boundary conditions for the three-dimensional body on the planes x 3 = ± c are accounted for at any x 1 of the one-dimensional beam continuum as resultants of these tractions. These resultants may be viewed as body forces for the one-dimensional continuum acting on interior points x 1 . Thus, the forces which act on the one-dimensional beam, or any sub-body of it, are section forces on the boundaries and resultant (body) forces distributed along x 1 . A convenient way to introduce these forces is to invoke variational methods.

4.2. Variational Methods of Gradient Elastic Euler-Bernoulli Beams

At the beginning of Section 3.3, it has been argued that the presence of acceleration terms in the traction boundary conditions is physically unacceptable. Nevertheless, for comparison, we shall elaborate variational formulations of gradient elastic Euler–Bernoulli beams for both cases of the KG-Model with and without acceleration terms in the boundary tractions.

4.2.1. Hamilton’s Principle for the Case Where Acceleration Terms Are Present in the Traction Boundary Conditions

Approach Based on Hamilton’s Principle Equation (40) 

It is instructive first to evaluate Hamilton’s principle (40) by specializing δ W ¯ e , Π   i , T ¯ to the assumed Euler–Bernoulli kinematics. Several algebraic manipulations are similar to those elaborated in Sideris and Tsakmakis [24] and Broese et al. [25], and therefore they will be only sketched briefly here.
Using steps quite similar to those in Sideris and Tsakmakis [24] (see also [25]), it can be seen that δ W ¯ e in Equations (41)–(43) can be expressed in terms of section forces N ¯ ,   H ¯ , V ¯ , M ¯ , m ¯ and resultant forces p ¯ , q ¯ as follows:
δ W ¯ e = V   P ¯ k δ u k + R ¯ k D δ u k dS   = N ¯ δ U x 1 = 0 x 1 = L + H ¯ δ U x 1 = 0 x 1 = L + V ¯ δ w x 1 = 0 x 1 = L + M ¯ δ w x 1 = 0 x 1 = L + m ¯ δ w x 1 = 0 x 1 = L + 0 L p ¯ δ U   dx 1 + 0 L q ¯ δ w   dx 1 .
The potential of the internal forces can be determined from the free energy in Equation (62),
Π   i = V   1 2 E ε 11 2 + l 2 1 ε 11 2 + l 2 3 ε 11 2 dV ,
or equivalently, by virtue of Equation (57),
Π   i = V   1 2 E U w x 3 2 + l 2 U w x 3 2 + l 2 w 2 dV .
If U   0 , then the potential (70) reduces to the one, which was essentially the starting point in the work of Lazopoulos and Lazopoulos [26]. After lengthy manipulations and repeatedly applying partial integration, we find that
δ Π   i = f 1 δ U x 1 = 0 x 1 = L + f 2 δ   U x 1 = 0 x 1 = L + f 3 δ w x 1 = 0 x 1 = L + f 3   δ w x 1 = 0 x 1 = L   + f 4   δ w x 1 = 0 x 1 = L 0 L f 1   δ U   dx 1 0 L f 3   δ w   dx 1
The parts T cl and T ¯ noncl , which make up the total kinetic energy T ¯ in Equation (39), are defined in Equations (2) and (34), and can be calculated with the aid of Equation (56):
Τ cl = V   1 2 ρ u ˙ k u ˙ k dV = 0 L 1 2 ρ A U ˙ 2 +   I w ˙ 2 + A   w ˙ 2     dx 1 ,
T ¯ noncl = V   1 2 γ j u ˙ k j u ˙ k dV = 0 L 1 2 γ A   U ˙ 2 + I w ˙ 2 + 2 A   w ˙ 2   dx 1 .
After repeated use of partial integration and the divergence theorem, we arrive at
t 1 t 2 δ T ¯ dt   = t 1 t 2 δ T cl + δ T ¯ noncl   dt   = t 1 t 2 { γ A U ¨ δ U x 1 = 0 x 1 = L ρ I + 2 γ A w ¨ γ I w ¨ δ w x 1 = 0 x 1 = L + γ I w ¨ δ w x 1 = 0 x 1 = L   + 0 L ρ A U ¨ + γ A U ¨ δ U + ρ A w ¨ + ρ I + 2 γ A w ¨ γ I w ¨ δ w dx 1 } dt .
The next step is to substitute formulas (68), (71) and (74) into Hamilton’s principle (40) and to apply the fundamental lemmas of calculus of variations. This gives the sectional constitutive laws
N ¯ = f 1 + γ A U ¨ = EA   U l 2 U + γ A U ¨ ,
H ¯ = f 2 = l 2 EA U ,
M ¯ = f 3 γ I w ¨ = EI   + l 2 EA w + l 2 EIw γ I w ¨ ,
m ¯ = f 4 = l 2 EI w ,
the sectional balance law for V ¯
V ¯ = f 3 γ Ι w ¨ + ρ Ι + 2 γ A w ¨ = EI + l 2 EA w + l 2 EIw γ I w ¨ + ρ I + 2 γ A w ¨  
V ¯ M ¯ = ρ I + 2 γ A w ¨ ,
the governing equations of motion for the axial displacement U
p ¯ + f 1 ρ A U ¨ + γ A U ¨ = 0  
ρ A U ¨ + γ A U ¨ + E A U l 2 E A U + p ¯ = 0 ,
the governing equation of motion for the deflection w
ρ A w ¨ + ρ I + 2 γ A w ¨ γ I w ¨ + f 3 + q ¯ = 0
ρ A w ¨ + ρ I + 2 γ A w ¨ γ I w ¨ EI + l 2 EA w + l 2 EIw + q ¯ = 0
and the boundary conditions
either N ¯ or   U ,   either H ¯ or U ,   either V ¯ or   w ,
either M ¯ or w and   either m ¯ or w
have to be prescribed at x 1 = 0   and   x 1 = L .
Now, it will be shown that the variational formulation (40), with δ W ¯ e , Π   i and T ¯ = T cl + T ¯ noncl being given in Equations (68), (69), (72) and (73), respectively, can be converted into an equivalent reduced form.

Approach Based on the Balance of Linear Momentum for the Beam 

Alternatively, one might establish a formulation of Hamilton’s principle on the basis of the balance equations of linear momentum for the beam in (63) and (64), which can be rewritten with the help of the stress components Σ ij t defined in Equation (37):
1 Σ 11 t + 3 Σ 31 t = ρ u ¨ 1 ,
1 Σ 13 t + 3 Σ 33 t = ρ u ¨ 3 ,
with
Σ 11 t = Σ 11 + γ U ¨ w ¨ x 3 ,   Σ 31 t = Σ 31 γ w ¨ ,
Σ 13 t = Σ 13 + γ w ¨   ,   Σ 33 t = Σ 33 .
Scalar multiplication of Equations (87) and (88) with the virtual displacement vector δu (see Equation (56)), and integration over V , yields
V   1 Σ 11 t + 3 Σ 31 t δ u 1 + 1 Σ 13 t + 3 Σ 33 t δ u 3 dV V   ρ ρ u ¨ 1 δ u 1 + ρ u ¨ 3 δ u 3 dV = 0 .
In order to replace the stress components by sectional and resultant forces, we introduce the following definitions
N ¯ r : = A   Σ 11 t dS   , H ¯ r   : = A   μ 111 dS , V ¯ r : = A   3 Σ 31 t   dS   ,
M ¯ r : = A   Σ 11 t x 3 dS   , m ¯ r   : = A   μ 111 x 3 dS ,
p ¯ r : = 2 b Σ 31 t x 3 = c x 3 = c , q ¯ r : = 2 b Σ 33 + 1 Σ 13 x 3 x 3 = c x 3 = c + γ A w ¨ .
Using these forces, it is shown in Appendix A that Equation (91) is equivalent to the following virtual work statement for the beam, which is of the form (12):
δ W ¯ r inert + δ W ¯ r e δ Π r i = 0 ,
with
δ W ¯ r inert : = 0 L ρ A U ¨ δ U + ρ I w ¨ δ w + ρ A w ¨ δ w dx 1 0 L γ A U ¨ δ U + γ I w ¨ δ w dx 1 ,
δ W ¯ r e   : = N ¯ r   δ U x 1 = 0 x 1 = L + H ¯ r δ U x 1 = 0 x 1 = L + V ¯ r δ w x 1 = 0 x 1 = L + M ¯ r   δ w x 1 = 0 x 1 = L + m ¯ r   δ w x 1 = 0 x 1 = L + 0 L p ¯ r   δ U   dx 1 + 0 L q ¯ r   δ w   dx 1 ,
Π r i : = V   1 2 E ε 11 2 + l 2 1 ε 11 2 dV .
It is readily seen, by using partial integration, that the time integral of δ W ¯ r inert between t 1 and t 2 can be represented in terms of kinetic energy T ¯ r ,
T ¯ r : = T cl + 0 L 1 2 γ A U ˙ 2 I w ˙ 2   dx 1
with T cl as given in Equation (72), i.e.,
t 1 t 2 δ W ¯ r inert   dt   =   δ t 1 t 2 T ¯ r   dt .
Thus, by taking the time integral of (95), we arrive at the variational formulation
δ t 1 t 2 T ¯ r Π r i   dt   + t 1 t 2 δ W ¯ r e   dt = 0 ,
which is of the form (10).
We would like to draw attention to the fact that Π r i in Equation (98) includes one term less than Π i in Equation (66). That is why we consider (101) as a reduced form of the principle (40).
Now, using partial integration repeatedly, the variation of Π r i may be expressed in terms of displacement components as follows:
δ Π r i = f 1 δ U x 1 = 0 x 1 = L + f 2 δ U x 1 = 0 x 1 = L + f 3 + l 2 EA w δ w x 1 = 0 x 1 = L + f 3 + l 2 EA w δ w x 1 = 0 x 1 = L + f 4   δ w x 1 = 0 x 1 = L 0 L f 1   δ U   dx 1 0 L f 3 + l 2 EAw δ wdx 1 .
Similarly, from Equation (96) we obtain
δ W ¯ r inert : = γ A U ¨ δ U x 1 = 0 x 1 = L ρ I w ¨ γ I w ¨ δ w x 1 = 0 x 1 = L γ I w ¨   δ w x 1 = 0 x 1 = L + 0 L ρ A U ¨ + γ A U ¨ δ Udx 1 + 0 L ρ A w ¨ + ρ I w ¨ + γ I w ¨ δ w   dx 1 .
Substitution of (97), (102) and (103) into (95) and application of the fundamental lemmas of calculus of variations leads to the sectional constitutive laws
N ¯ r = f 1 + γ A U ¨ = EA U l 2 U + γ A U ¨ ,
H ¯ r = f 2 = l 2 EA U ,
M ¯ r = f 3 + l 2 EA w γ I w ¨ = EI w l 2 w γ I w ¨   ,
m ¯ r = f 4 = l 2 EI w ,
the sectional balance law for V ¯ r
V ¯ r = f 3 + l 2 EA w + ρ I w ¨ γ I w ¨ = EI w + l 2 EIw + ρ I w ¨ γ I w ¨
V ¯ r M ¯ r = ρ I w ¨ ,
the governing equation of motion for the axial displacement U
p ¯ r + f 1 ρ A U ¨ + γ A U ¨ = 0  
ρ A U ¨ + γ A U ¨ + EA U l 2 EAU + p ¯ r = 0 ,
the governing equation of motion for the deflection w
q ¯ r + f 3 + l 2 EAw ρ A w ¨ + ρ I w ¨ γ I w ¨ = 0
ρ A w ¨ + ρ I w ¨ γ I w ¨ EIw + l 2 EIw + q ¯ r = 0
and the boundary conditions
either   N ¯ r   or   U ,   either H ¯ r   or U ,   either V ¯ r   or   w ,
either   M ¯ r   or w and   either m ¯ r   or w
have to be prescribed at x 1 = 0 and x 1 = L .
It should be noted that the sectional constitutive laws (104)–(107) can alternatively be derived from definitions (92) and (93). Also, a bending theory based on Π r i in Equation (98) has been formulated for the first time in Papargyri-Beskou et al. [27].
The two beam approaches considered in this section are equivalent to each other. To see this, it suffices to show that the second approach implies the first one. Let p 0 be a given loading function and assume that p ¯ = p ¯ r = p 0 . Then, the two equations of motion (82) and (111) are identical. Further, assume q ¯ r to have the form q ¯ r = q 0 l 2 EAw + 2 γ A w ¨ with q 0 being a known external body force. Then, the equation of motion (113) implies the equation of motion (84) with q ¯ = q 0 . The displacement boundary conditions and the traction boundary conditions for N ¯ , H ¯ , m , ¯ and N ¯ r , H ¯ r , m ¯ r are identical for both approaches. By comparing Equation (79) with (108) and Equation (77) with (106), we recognize that V ¯ = V ¯ r l 2 EA w + 2 γ A w ¨ and M ¯ = M ¯ r l 2 EA w . Thus, assume that at the boundaries x 1 = 0 ,   L , V ¯ r and M ¯ r have the forms V ¯ r = V 0 + l 2 EA w 2 γ A w ¨ , M ¯ r = M 0 + l 2 EA w ,   with   V 0 , M 0 being given. Then, to the boundary conditions for V ¯ r and M ¯ r correspond the boundary conditions V ¯ = V 0 and M ¯ = M 0 . Note that such boundary conditions include acceleration terms, but this is characteristic of Mindlin’s approach. Finally, by adding the time integral of
δ V   1 2 l 2 E 3 ε 11 2 dV + l 2 EA w δ w x 1 = 0 x 1 = L + l 2 EA w δ w x 1 = 0 x 1 = L 0 L l 2 EAw δ w   dx 1 = 0  
and the identity
δ t 1 t 2 0 L γ A w ˙ 2 dx 1 dt + t 1 t 2 0 L 2 γ A w ¨ δ wdx 1 dt δ t 1 t 2 2 γ A w ¨ δ w x 1 = 0 x 1 = L dt = 0
to Hamilton’s principle (101), and replacing V ¯ r l 2 EA w + 2 γ A w ¨ by V ¯ and M ¯ r l 2 EA w by M ¯ , we can deduce that (101) implies Hamilton’s principle (40) with δ W ¯ e , Π   i and T ¯ = T cl + T ¯ noncl as given in Equation (68), (70), (72) and (73), respectively.

4.2.2. Hamilton’s Principle for the Case Where Acceleration Terms Are Not Present in the Traction Boundary Conditions

Similar to Section 4.2.1, again, there are two equivalent methods for deriving the equations of motion. The first one starts with Hamilton’s principle (48), assuming δ W ^ e and δ W ^ inert ,   noncl are as defined in Equations (49) and (47), respectively. The main feature of this method is that it deals with the potential for the internal forces Π   i in Equation (69). The second one starts with the equations of linear momentum (63) and (64) and leads to the reduced form of the potential of internal forces Π r i in Equation (98). The calculations are quite similar to those in Section 4.2.1. In the examples below, we shall employ only the equations of motion derived according to the second method. Therefore, only this method will be discussed here. Accordingly, some quantities are designated with the subscript r.
In the last section, the balance of linear momentum of the beam (Equations (87) and (88)), expressed in terms of the total Cauchy stress tensor Σ t , was the start point. Opposite to this, we now start with the balance of linear momentum in Equations (63) and (64), which is expressed in terms of the proper Cauchy stress Σ:
Σ 11 + 3 Σ 13 + i 1 = 0 ,
Σ 13 + 3 Σ 33 + i 3 = 0 ,
with the components of the inertial force i being given by
i 1 = i 1 cl + i 1 noncl ,   i 1 cl = ρ u ¨ 1 ,   i 1 noncl = γ u ¨ 1 ,
i 3 = i 3 cl + i 3 noncl ,   i 3 cl = ρ u ¨ 3 ,   i 1 noncl = γ u ¨ 3 .
By scalar multiplication of (118) and (119) with virtual displacement (see Equation (56)) and integration over V, we find that
V   Σ 11 + 3 Σ 13 δ u 1 + Σ 13 + 3 Σ 33 δ u 3 dV + δ W r inert = 0 .
In this equation, δ W r inert is the virtual work of the inertial force composed of a classical and a non-classical term,
δ W r inert = δ W inert , cl + δ W r inert , noncl ,
δ W r inert , cl = V   i i cl δ u i dV = V   ρ u ¨ 1 δ u 1 + ρ u ¨ 3 δ u 3 dV = 0 L ( ρ A U ¨ δ U + ρ I w ¨ δ w + ρ A w ¨ δ w )   dx 1 = ρ I w ¨ δ w x 1 = 0 x 1 = L + 0 L [ ρ A U ¨ δ U ρ A w ¨ ρ I w ¨ δ w ]   dx 1 ,
δ W r inert , noncl = V   i i noncl δ u i dV = V   γ u ¨ 1 δ u 1 + γ u ¨ 3 δ u 3 dV = 0 L ( γ A U ¨ δ U + γ I w ¨ δ w )   dx 1 = γ I w ¨ δ w x 1 = 0 x 1 = L + 0 L γ A U ¨ δ U γ I w ¨ δ w dx 1 .
The appropriate definitions of sectional and resultant forces are now
N r : = A   Σ 11 dS ,   H r : = A   μ 111 dS ,   V r : = A   3 Σ 13 x 3 dS ,
M r : = A   Σ 11 x 3 dS ,   m r : = A   μ 111 x 3 dS ,
p r : = 2 b   Σ 13 x 3 = c x 3 = c   ,   q r : = 2 b   Σ 33 + 1 Σ 13 x 3 x 3 = c x 3 = c .
It is worth noting that these definitions apply equally in statics (cf. [25]). With the help of these definitions, it can be shown (see Appendix B), that Equation (122) implies
δ W r inert + δ W r e δ Π r i = 0 ,
where
δ W r e = N r   δ U x 1 = 0 x 1 = L + H r δ U x 1 = 0 x 1 = L + V r δ w x 1 = 0 x 1 = L + M r δ w x 1 = 0 x 1 = L + m r δ w x 1 = 0 x 1 = L + 0 L p r δ Udx 1 + 0 L q r δ wdx 1 .
Equation (129) is a virtual work principle of the form (12). We recall the well-known identity (cf. Equation (125))
t 1 t 2 δ W inert , cl dt = δ t 1 t 2 1 2 [ ρ A U ˙ 2 + ρ I ( w ˙ ) 2 + ρ A w ˙ 2 ] dt = δ t 1 t 2 T cl dt ,
and that Equation (123) applies. Thus, time integration of Equation (129) furnishes
δ t 1 t 2 T cl Π   i dt + t 1 t 2 δ W r e + δ W r inert , noncl dt = 0 ,
which is a variant of Hamilton’s principle of the form (11).
On the other hand, we can invoke in Equation (129) the result (102) for δ Π r i and the result (123), (124)3 and (125)3 for δ W r inert , and by applying the fundamental lemmas of calculus for variations to infer the sectional constitutive laws
N r = f 1 = EA U l 2 U ,
H r = f 2 = l 2 EA U ,
M r = f 3 + l 2 EA w = EI w l 2 w ,
m r = f 4 = l 2 EI w ,
the sectional balance law for V r
V r = f 3 + l 2 EA w + ρ I w ¨ γ I w ¨ = Ε Ι w + l 2 EI   w + ρ I w ¨ γ I w ¨  
V r M r = ρ I w ¨ γ I w ¨ ,
the governing equation of motion for the axial displacement U
p r + f 1 ρ A U ¨ + γ A U ¨
ρ A U ¨ +   γ A U ¨ + EA U l 2 EAU + p r = 0 ,
the governing equation of motion for the deflection w
q r + f 3 + l 2 EA   w ρ A w ¨ + ρ I w ¨ γ I w ¨ = 0
ρ A w ¨ + ρ I w ¨ γ I w ¨ EIw + l 2 EI   w + q r = 0  
and the boundary conditions
either   N r   or   U ,   either   H r   or U ,   either   V r   or   w ,
either   M r   or   w and   either   m r   or   w
have to be prescribed at x 1 = 0   and x 1 = L   .
We see from Equations (104)–(113) and (133)–(142), that differences between the two approaches only exist in the boundary conditions for the sectional forces; some contain acceleration terms while others do not.

5. Examples

The proper aim of the examples is to demonstrate with reference to bending loading that the presence of acceleration terms in the boundary conditions may have significant qualitative and quantitative effect on the predicted responses. Note, however, that Euler–Bernoulli beam theories are especially tempting for they render the equations of motion one-dimensional. That naturally arouses interest to compare with each other predicted responses according to one-dimensional tension/compression and bending loadings. For our purposes, it suffices to perform the comparison only with respect to the beam approaches, which rely upon the reduced form of the potential of internal forces. Additionally, we shall concentrate on the following three versions of the theory.
Version 1: γ = 0, i.e., nonclassical acceleration terms are omitted.
Version 2: γ 0 , nonclassical acceleration terms are present in both the equations of motion and the traction boundary conditions.
Version 3: γ 0 , nonclassical acceleration terms are present only in the equations of motion.
In the case of Version 1, the equations of motion follow from (140) and (142) by setting γ = 0, while the boundary conditions are the same as for Version 3. Furthermore, the definitions of sectional and resultant distributed forces for Version 1 are the same as for Version 3. The predicted responses will be presented in dimensionless form by employing the definitions
x ˜ : = x 1 L   , t ˜ : = c * L   t   ,   x ˜ =   ,   x ˜   ,     t ˜ =   ,   t ˜   , U ˜ : = U L   , w ˜ : = w L   , N ˜ : = N EA ,
V ˜ : = V EL 2   , M ˜ : = M EL 3   , γ ˜ : = γ c * 2 EL 2   ,   c * = E ρ   , A ˜ : = A L 2 , I ˜ : = I L 4 ,   l ˜ : = l L , ω ˜ : = ω L c * .
Here, N =   N r , V =   V r , M =   M r for Versions 1, 3, while N = N ¯ r , V = V ¯ r   ,   M   = M ¯ r for Version 2. We shall compare the three versions with each other for the case of a cantilever beam and for harmonically with time-varying loading conditions. In all calculations, the values of A ˜ = 15   · 10 4 and I ˜ = 2.813   · 10 9 have been chosen.

5.1. Uniaxial Tension/Compression Loading

5.1.1. Governing Equations

Consider a cantilever beam, subject only to axial load, so that w = 0, and let p ¯ r = p r = 0 . Then, from Equation (111) or (140), we find that
  ρ U ¨ + γ U ¨ + E U l 2 E   U = 0 ,
or in dimensionless form
U ˜ , t ˜ t ˜ + γ ˜ U ˜ , x ˜ x ˜ t ˜ t ˜ + U ˜ , x ˜ x ˜ l ˜ 2 U ˜ , x ˜ x ˜ x ˜ x ˜ = 0 .
If γ ˜ = 0 ,   then (148) represents the equation of motion for Version 1 and if γ ˜ 0 , then (148) is the equation of motion for versions 2, 3.
A comprehensive discussion of Equation (148), with respect to size effects and the convergence behavior for l ˜ 0 , has been provided in Broese et al. [15], but with nonclassical boundary conditions different from those we shall assume in the present paper. Especially, the interest here is focused on homogenous nonclassical traction boundary tractions. The important loading condition is at x 1 = L and has the form B A e i ω t , with ω being an operating frequency, B A being a displacement-or force-like amplitude and i being the imaginary unit. This kind of loading suggests assuming for the solution of Equation (148) the form
U ˜ x ˜ , t ˜ = U ˜ 0 x ˜   e i ω ˜ t ˜   .  
After substitution of this into Equation (148), and elimination of the factor e i ω ˜ t ˜ , we obtain
ω ˜ U ˜ 0 + 1 γ ˜ ω ˜ 2 ( U ˜ 0 ) , x ˜ x ˜ l ˜ 2 ( U ˜ 0 ) , x ˜ x ˜ x ˜ x ˜ = 0 ,
which, along with a set of boundary conditions, can be solved by employing standard methods. Having available solution of displacement (149), it is straightforward to establish solutions of sectional forces from the corresponding formulae. We will now discuss force and displacement-controlled loadings for a cantilever beam.

5.1.2. Force Controlled Loading

The boundary conditions for Versions 1, 3 are U 0 ,   t = 0 ,   H r 0 , t = H r L , t = 0 and N L , t = F A e i ω t , where F A = constant is a force amplitude. From these, we can gain boundary conditions for the differential Equation (150) by taking Equations (133), (134) and (149) into account, eliminating the factor e i ω t , and using dimensionless variables. In particular, the dimensionless expression of the boundary condition for N becomes N ˜ 1 , t ˜ = F ˜ A e i ω ˜ t ˜ , with F ˜ A = F A EA . Moreover, it is readily seen from Equation (133), that, for Versions 1, 3, N ˜ = N ˜ 0 x ˜   e i ω ˜ t ˜ , with N ˜ 0 = ( U ˜ 0 ) , x ˜ l ˜ 2 ( U ˜ 0 ) , X ˜ X ˜ X ˜ . With the help of Equations (104) and (105) and using similar manipulations as for Versions 1, 3, it is straightforward to establish corresponding boundary conditions and distributions of N ˜ 0 for Version 2, where again N ˜ = N ˜ 0 x ˜ e i ω ˜ t ˜ . Altogether, we have the following boundary conditions and distributions of N ˜ 0 .
Versions 1, 3
B . C . : U ˜ 0 x ˜ = 0 = [ ( U ˜ 0 ) , x ˜ x ˜ ] x ˜ = 0 = [ ( U ˜ 0 ) , x ˜ x ˜ ] x ˜ = 1 = 0 ,
[ ( U ˜ 0 ) , x ˜ l ˜ 2 ( U ˜ 0 ) , x ˜ x ˜ x ˜ ] x ˜ = 1 = F ˜ A .
Solution   N ˜ 0 : N ˜ 0 x ˜ = ( U ˜ 0 ) , x ˜ l ˜ 2 ( U ˜ 0 ) , x ˜ x ˜ x ˜
Version 2
B . C . : U ˜ 0 x ˜ = 0 = [ ( U ˜ 0 ) , x ˜ x ˜ ] x ˜ = 0 = [ ( U ˜ 0 ) , x ˜ x ˜ ] x ˜ = 1 = 0 ,
[ 1 γ ˜ ω ˜ 2 ( U ˜ 0 ) , x ˜ l ˜ 2 ( U ˜ 0 ) , x ˜ x ˜ x ˜ ] x ˜ = 1 = F ˜ A .
Solution   N ˜ 0 : N ˜ 0 x ˜ = 1 γ ˜ ω ˜ 2 ( U ˜ 0 ) , x ˜ l ˜ 2 ( U ˜ 0 ) , x ˜ x ˜ x ˜ .
Resulting distributions of U ˜ 0 and N ˜ 0 are illustrated in Figure 2 and Figure 3. In addition, corresponding distributions predicted by classical elasticity are shown in these figures.
It is well known that gradient elasticity includes parameters which control length-scale effects captured by the constitutive theory. In the case of Version 1 ( γ ˜ = 0 ), the only material parameter responsible for length scale effects is the internal material length l ˜ . Figure 2 illustrates the effect of l ˜ on the predicted model responses for the case of Version 1. It can be seen that for small values of frequencies as, e.g., ω ˜ = 1.5 , all distributions of U ˜ 0 are monotonically increasing, do not intersect for x ˜ > 0   and indicate the gradient stiffening effect in comparison to the classical solution. The stiffening effect is increasing with increasing values of l ˜ . All corresponding N ˜ 0 -distributions, shown in Figure 2b, are monotonically decreasing, do not intersect for x ˜ < 1 and are below the classical one. Thus, the distributions of U ˜ 0 and N ˜ 0 indicate a common intersection point, respectively.
No regular tendencies in the distributions of U ˜ 0 and N ˜ 0 can be recognized, or even the opposite may happen, for sufficiently large values of ω ˜ . In particular, new intersection points in the graphs may occur and positions of intersection points can change depending on the applied frequency. Since the relationships for large values of ω ˜ are similar to those reported in Broese et al. [15], they will not be further discussed here.
In the case of Versions 2 and 3 ( γ ˜   0 ), length scale effects can be controlled besides by l ˜ , also by γ ˜ (see Figure 3). Apparently, for both versions, and for the sufficiently small frequency ω ˜ = 1.5 , the U ˜ 0 - and N ˜ 0 -distributions, for l ˜ = constant, look similar to the ones for Version 1 in Figure 2. Moreover, the U ˜ 0 -distributions increase with increasing values of γ ˜ , and can exceed the one predicted by classical elasticity. This, in turn, indicates that, for l ˜ = constant, the nonclassical acceleration terms controlled by γ ˜ , cause (dynamical) softening in the material behavior. However, it must be emphasized that the amounts of these distributions are considerably smaller for Version 3. This behavior carries over to the N ˜ 0 -distributions as well. Altogether, the predicted responses by Version 2 and Version 3 are qualitatively similar to each other but, depending on the γ ˜ -values, significant quantitative differences can occur.

5.1.3. Displacement Controlled Loading

All the boundary conditions and the solutions for N ˜ 0 are the same as in the last section, except for the boundary condition at x ˜ = 1 for the sectional force N ˜ , which is now replaced by the displacement boundary condition U ˜ 1 , t ˜ = U ˜ A e i ω ˜ t ˜ ,   U ˜ A = constant.
Versions 1, 3
B . C . : U ˜ 0 x ˜ = 0 = [ ( U ˜ 0 ) , x ˜ x ˜ ] x ˜ = 0 = [ ( U ˜ 0 ) , x ˜ x ˜ ] x ˜ = 1 = 0   ,   U ˜ 0 x ˜ = 1 = U ˜ A
Solution   N ˜ 0 : N ˜ 0 x ˜ = ( U ˜ 0 ) , x ˜ l ˜ 2 ( U ˜ 0 ) , x ˜ x ˜ x ˜ .
Version 2
B . C . : U ˜ 0 x ˜ = 0 = [ ( U ˜ 0 ) , x ˜ x ˜ ] x ˜ = 0 = [ ( U ˜ 0 ) , x ˜ x ˜ ] x ˜ = 1 = 0   ,   U ˜ 0 x ˜ = 1 = U ˜ A .
Solution   N ˜ 0 : N ˜ 0 x ˜ = 1 γ ˜ ω ˜ 2 ( U ˜ 0 ) , x ˜ l ˜ 2 ( U ˜ 0 ) , x ˜ x ˜ x ˜ .
Predicted distributions for Version 1 ( γ ˜ = 0) and classical elasticity for the above displacement-controlled boundary conditions are depicted in Figure 4 for frequency ω ˜ and U ˜ A = 5 · 10 3 . The general observations concerning U ˜ 0 -responses are similar to those for force-controlled loading. Clearly, due to the imposed boundary conditions, the U ˜ 0 -distributions intersect now at x ˜ = 1 as well. It can be recognized from Figure 4a, that for the sufficiently small value ω ˜ = 1.5 , only small quantitative differences are visible, which could be expected because of the assumed displacement boundary conditions. In the related N ˜ 0 -distributions only small quantitative differences are visible as well (cf. Figure 4b). However, there is the remarkable qualitative difference that, now, these distributions do intersect for some 0 < x ˜ < 1 . For sufficiently large values of ω ˜ no regular tendencies in the predicted responses can be stated. As the differential Equation (148) for U ˜ 0 and the associated boundary conditions (155) and (157) are identical for Versions 2 and 3, the predicted U ˜ 0 -distributions according to both versions are identical as well. The graphs of U ˜ 0 -distributions in Figure 5 attest, for the small frequency ω ˜ = 1.5 , only small quantitative differences for various values of γ ˜ . After magnification of Figure 5, it becomes clear that these distributions increase with increasing values of γ ˜ but are always below the one predicted by classical elasticity. On the other hand, observing closely the N ˜ 0 -distributions in Figure 6, it becomes clear that these distributions for Version 2 decrease with increasing values of γ ˜ and intersect the one predicted by classical elasticity for some x ˜ 0 , 1 . Comparison with N ˜ 0 -distributions for Version 3, depicted in Figure 6b, reveal for Version 2 smaller values of N ˜ 0 than for Version 3. In addition, by increasing the value of γ ˜ , these distributions in Figure 6b are increasing. This observation is the opposite of the one made in Figure 6a.
Concluding this section, we can state that there are significant quantitative and qualitative differences for the considered tension/compression loading conditions.

5.2. Cantilever Beam under Dynamical Transverse Load

5.2.1. Governing Equations

Consider, now, the cantilever beam to be loaded only transversely so that U 0 , and suppose that, q ¯ r = q r 0 . For this case, we find from Equation (113) (or from Equation (142)) that
ρ A w ¨ + ρ I w ¨ γ I   w ¨ EIw + l 2 EI   w = 0 .
The dimensionless form of it reads
A ˜ w ˜ , t ˜ t ˜ + I ˜ w ˜ , t ˜ t ˜ x ˜ x ˜ γ ˜ I ˜ w ˜ , t ˜ t ˜ x ˜ x ˜ x ˜ x ˜ I ˜ w ˜ , x ˜ x ˜ x ˜ x ˜ + l ˜ 2 I ˜ w ˜ , x ˜ x ˜ x ˜ x ˜ x ˜ x ˜ = 0 .
Obviously, for γ ˜ = 0 , Equation (160) represents the equation of motion for Version 1, and for γ ˜ 0 it represents the equation of motion for Versions 2, 3. We assume that, at x 1 = L , the beam is subject to transverse load, force or deflection controlled, of the form Be ω it , where again ω is an operating frequency, B is a force-or deflection-like amplitude and i is the imaginary unit. This suggests for the solution of Equation (160) to make the Ansatz
w ˜ x ˜ ,   t ˜ = w ˜ 0 x ˜ e i ω ˜ t ˜ .
After substitution of (161) into (160), and elimination of the factor e i ω ˜ t ˜ , we arrive at
A ˜ ω ˜ 2 w ˜ 0 ω ˜ 2 Ι ˜ ( w ˜ 0 ) , x ˜ x ˜ + Ι ˜ γ ˜ ω ˜ 2 1 ( w ˜ 0 ) , x ˜ x ˜ x ˜ x ˜ + l ˜ 2 Ι ˜ ( w ˜ 0 ) , x ˜ x ˜ x ˜ x ˜ x ˜ x ˜ = 0 ,
which, along with a set of boundary conditions, can be solved by employing standard methods. Distribution of sectional forces can be established by inserting into corresponding formulae the solutions w ˜ x ˜ ,   t ˜ . In the subsequent sections, we will discuss force and deflection controlled bending of the cantilever beam.

5.2.2. Force Controlled Bending

The manipulations in the current and the next section are similar to those in Section 5.1.2 and Section 5.1.3, respectively. We suppose for Versions 1, 3 the boundary conditions w 0 , t = 0 , w 0 ,   t = 0 , m r 0 ,   t = m r L ,   t = 0 , M r L ,   t = 0 and V L , t = V r L , t = F A e i ω t , where again F A = constant is a force amplitude. From these boundary conditions, we can deduce boundary conditions for the differential Equation (162) by taking into account Equations (135)–(138) and (161), and eliminating the factor e i ω ˜ t ˜ . In particular, the boundary condition for V ˜ becomes V ˜ 1 ,   t ˜ = F ˜ A e i ω ˜ t ˜ , where now F ˜ A : = F A E L 2 . Also, we can deduce from Equations (135) and (137) that, for Versions 1, 3, V ˜ = V ˜ 0 x ˜ e i ω ˜ t ˜ , M ˜ = M ˜ 0 x ˜ e i ω ˜ t ˜ , with V ˜ 0 x ˜ = ω ˜ 2 Ι ˜ w ˜ 0 , x ˜ + I ˜ γ ˜ ω ˜ 2 1 w ˜ 0 , x ˜ x ˜ x ˜ + l ˜ 2 I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ x ˜ and M ˜ 0 x ˜ = Ι ˜ w ˜ 0 , x ˜ x ˜ + l ˜ 2   I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ .
With the help of Equation (106)–(108), and using similar manipulations as for Versions 1, 3, we can deduce corresponding boundary conditions and distributions V ˜ 0 , M ˜ 0 for Version 2. Altogether, we have the following.
Versions 1, 3
B . C . : w ˜ 0 x ˜ = 0 = ( w ˜ 0 , x ˜ ] x ˜ = 0 = ( w ˜ 0 , x ˜ x ˜ x ˜ ] x ˜ = 0 = ( w ˜ 0 , x ˜ x ˜ x ˜ ] x ˜ = 1 = ( w ˜ 0 , x ˜ x ˜ l ˜ 2 ( w ˜ 0 ) , x ˜ x ˜ x ˜ x ˜ ] x ˜ = 1 = 0 ,
ω ˜ 2 Ι ˜ w ˜ 0 , x ˜ + Ι ˜ γ ˜ ω ˜ 2 1 w ˜ 0 , x ˜ x ˜ x ˜ + l ˜ 2 I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ x ˜ x ˜ = 1 = F ˜ A .
Solution   V ˜ 0 : V ˜ 0 x ˜ = ω ˜ 2 Ι ˜ w ˜ 0 , x ˜ + Ι ˜ γ ˜ ω ˜ 2 1 w ˜ 0 , x ˜ x ˜ x ˜ + l ˜ 2 I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ x ˜ .
Solution   M ˜ 0 : M ˜ 0 x ˜ = Ι ˜ w ˜ 0 , x ˜ x ˜ + l ˜ 2 I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ .
Version 2
B . C . : w ˜ 0 x ˜ = 0 = ( w ˜ 0 , x ˜ ] x ˜ = 0 = ( w ˜ 0 , x ˜ x ˜ x ˜ ] x ˜ = 0 = ( w ˜ 0 , x ˜ x ˜ x ˜ ] x ˜ = 1 = [ Ι ˜ γ ˜ ω ˜ 2 1 ( w ˜ 0 ) , x ˜ x ˜ + l ˜ 2 ( w ˜ 0 ) , x ˜ x ˜ x ˜ ] x ˜ = 1 = 0 ,
ω ˜ 2 Ι ˜ w ˜ 0 , x ˜ + Ι ˜ γ ˜ ω ˜ 2 1 w ˜ 0 , x ˜ x ˜ x ˜ + l ˜ 2 I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ x ˜ x ˜ = 1 = F ˜ A .
Solution   V ˜ 0 : V ˜ 0 x ˜ = ω ˜ 2 Ι ˜ w ˜ 0 , x ˜ + Ι ˜ γ ˜ ω ˜ 2 1 w ˜ 0 , x ˜ x ˜ x ˜ + l ˜ 2 I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ x ˜ .
Solution   M ˜ 0 : M ˜ 0 x ˜ = Ι ˜ 1 γ ˜ ω ˜ 2 w ˜ 0 , x ˜ x ˜ + l ˜ 2 I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ .
By comparing Figure 7 with Figure 2, it can be recognized that the w ˜ 0 - and V ˜ 0 -distributions are essentially similar to the U ˜ 0 - and N ˜ 0 -distributions, respectively, the only difference being that the w ˜ 0 -distributions are convex whereas the U ˜ 0 -distributions are concave.
Also, similar to the uniaxial loading cases in Figure 3a,c, for the small frequency ω ˜ = 0.05 , and for keeping l ˜ constant, by increasing values of γ ˜ , the w ˜ 0 -distributions predicted by Versions 2, 3 are increasing, which indicates that dynamic gradient softening can occur (see Figure 8a,c). However, the amounts of the w ˜ 0 -distributions predicted by Version 3 are smaller than those predicted by Version 2. Furthermore, for forced controlled loadings, in the case of Version 3, there is a noticeable qualitative difference between the U ˜ 0 -distributions in Figure 3c and the w ˜ 0 -distributions in Figure 8c: All w ˜ 0 -distributions are bounded from above by the classical one, in opposite to the U ˜ 0 -distributions.

5.2.3. Deflection Controlled Bending

All the boundary conditions and the solutions V ˜ 0 , M ˜ 0 are the same as in the last section, except for the boundary conditions at x ˜ = 1 for the sectional forces V ˜ , which are now replaced by the deflection boundary conditions w ˜ 1 ,   t ˜ = w ˜ A e i ω ˜ t ˜ , w ˜ A = constant .
Versions 1, 3
B . C . : w ˜ 0 x ˜ = 0 = ( w ˜ 0 , x ˜ ] x ˜ = 0 = ( w ˜ 0 , x ˜ x ˜ x ˜ ] x ˜ = 0 = ( w ˜ 0 , x ˜ x ˜ x ˜ ] x ˜ = 1 = ( w ˜ 0 , x ˜ x ˜ l ˜ 2 ( w ˜ 0 ) , x ˜ x ˜ x ˜ x ˜ ] x ˜ = 1 = 0 ,   w ˜ 0 x ˜ = 1 = w ˜ A .
Solution   V ˜ 0 : V ˜ 0 x ˜ = ω ˜ 2 Ι ˜ w ˜ 0 , x ˜ + Ι ˜ γ ˜ ω ˜ 2 1 w ˜ 0 , x ˜ x ˜ x ˜ + l ˜ 2 I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ x ˜ .
Solution   M ˜ 0 : M ˜ 0 x ˜ = Ι ˜ w ˜ 0 , x ˜ x ˜ + l ˜ 2 I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ .
Version 2
B . C . : w ˜ 0 x ˜ = 0 = ( w ˜ 0 , x ˜ ] x ˜ = 0 = ( w ˜ 0 , x ˜ x ˜ x ˜ ] x ˜ = 0 = [ ( w ˜ 0 ) , x ˜ x ˜ x ˜ ] x ˜ = 1 = [( Ι ˜ γ ˜ ω ˜ 2 1 ( w ˜ 0 ) , x ˜ x ˜ + l ˜ 2 ( w ˜ 0 ) , x ˜ x ˜ x ˜ x ˜ ] x ˜ = 1 ,   w ˜ 0 x ˜ = 1 = w ˜ A .
Solution   V ˜ 0 : V ˜ 0 x ˜ = ω ˜ 2 Ι ˜ w ˜ 0 , x ˜ + Ι ˜ γ ˜ ω ˜ 2 1 w ˜ 0 , x ˜ x ˜ x ˜ + l ˜ 2 I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ x ˜ .
Solution   M ˜ 0 : M ˜ 0 x ˜ = Ι ˜ 1 γ ˜ ω ˜ 2 + l ˜ 2 I ˜ w ˜ 0 , x ˜ x ˜ x ˜ x ˜ .
Figure 9 illustrates predicted w ˜ 0 - and V ˜ 0 -distributions according to Version 1 ( γ ˜ = 0 ) and classical elasticity for w ˜ = 0.05 and w ˜ A = 5 · 10 3 . Except for the differences in the curvatures, the w ˜ 0 – distributions in Figure 9a look similar to the U ˜ 0 -distributions in Figure 4a. However, there are significant qualitative and quantitative differences between the V ˜ 0 -distributions in Figure 9b and the N ˜ 0 -distributions in Figure 4b. The most important one is that there are no intersections between the graphs of the curves in Figure 9b. The remarks concerning the predicted responses according to Versions 2, 3 in Figure 10 and Figure 11 are similar to those for the uniaxial loadings in Figure 5 and Figure 6.
The main conclusion which can be drawn from the calculated responses in the present and the last section is that, depending on the imposed boundary conditions, both similarities and differences in the predicted distributions due to uniaxial and bending loadings may occur. The most important observation, however, is that if γ 0, then significant qualitative and quantitative differences in the predicted responses can occur, depending on whether acceleration terms are present in the boundary conditions or not.

6. Concluding Remarks

The explicit gradient elasticity proposed by Mindlin [7] was an important step toward modeling nonlocalities and length scale effects in the material response. The original formulation of this theory is characterized by the presence of acceleration terms in the boundary tractions. It was argued in [15] that such boundary tractions are physically unacceptable as they are not objective. The argumentation against boundary tractions including acceleration terms has been extended in the present paper within the context of the principle of material frame indifference. It is shown that if this principle is assumed to hold, then boundary tractions including acceleration terms are not acceptable since either these tractions will not be objective, or the response functions of the associated stresses will not satisfy the principle of material frame indifference. The differences in the model responses, according to the two different forms of boundary tractions, are not negligible. In fact, it has been demonstrated with reference to uniaxial and bending loadings, that significant qualitative and quantitative differences may occur between the corresponding model responses. Different to [15], homogenous nonclassical boundary tractions are assumed here in the calculated examples, which, for the time being, seems to be physically more realistic. Besides force controlled loading histories, displacement controlled ones have also been considered and compared with each other. Thus, the present paper, compared with the analysis given [9,10], offers a more comprehensive discussion of the boundary conditions, and in addition, for the first time, it provides a consistent Euler–Bernoulli beam theory for bending in dynamics.

Author Contributions

Conceptualization, C.T., C.B. and S.A.S.; Methodology C.T., C.B. and S.A.S.; Software, C.T., C.B. and S.A.S.; Validation, C.T., C.B. and S.A.S.; Formal Analysis, C.T., C.B. and S.A.S.; Investigation, C.T., C.B. and S.A.S.; Resources, C.T., C.B. and S.A.S.; Data Curation, C.T., C.B. and S.A.S.; Writing—Original Draft Preparation, C.T., C.B. and S.A.S.; Writing—Review & Editing, C.T., C.B. and S.A.S.; Visualization, C.T., C.B. and S.A.S.; Supervision, C.T., C.B. and S.A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Substitution of Equation (56) into Equations (89) and (91) gives
V   1 Σ 11 t δ UdV + V   1 Σ 11 t x 3 δ   w dV + V   3 Σ 31 t δ UdV V   1 Σ 11 t δ UdV + V   1 Σ 11 t x 3 δ w dV + V   3 Σ 31 t δ UdV V   ρ u ¨ 1 δ u 1 + u ¨ 3 δ u 3 dV = 0 .
The first term can be calculated with the aid of the constitutive law (58), (59) and (89) and the definitions (92):
V   1 Σ 11 t δ UdV = V   1 Σ 11 t δ U dV V   Σ 11 t δ U dV = 0 L 1 N ¯ r δ U dx 1 V   Σ 11 δ U dV V   γ U ¨ w ¨ x 3 δ U dV = N ¯ r   δ U x 1 = 0 x 1 = L V   E ε 11 δ U dV + V   μ 111 δ U dV 0 L γ A U ¨ δ U dx 1 = N ¯ r   δ U x 1 = 0 x 1 = L V   E ε 11 δ U dV + V   1 A   μ 111 dS   δ U dx 1 V   μ 111 δ U dV 0 L γ A U ¨ δ U dx 1 = N ¯ r   δ U x 1 = 0 x 1 = L + H ¯ r   δ U x 1 = 0 x 1 = L V   E ε 11 δ U + μ 111 δ U dV 0 L γ A U ¨ δ U dx 1 .
The second term in Equation (A1) can be evaluated in a similar fashion:
V   1 Σ 11 t x 3 δ w dV = V   1 Σ 11 t x 3 δ w dV V   Σ 11 t x 3 δ w dV = 0 L 1 M ¯ r δ w dx 1 V   E ε 11 x 3 δ w dV + V   μ 111 x 3 δ w dV V   γ U ¨ w ¨ x 3 x 3 δ w dV = M ¯ r   δ w x 1 = 0 x 1 = L + m ¯ r   δ w x 1 = 0 x 1 = L V   E ε 11 x 3 δ w + μ 111 δ w dV 0 L γ I w ¨ δ w dx 1 .  
The third term in Equation (A1) gives
V   3 Σ 31 t δ UdV = 0 L 2 b Σ 31 t x 3 = c x 3 = c δ Udx 1 = 0 L p ¯ r δ Udx 1 .
For the other integrals in Equation (A1) including stress components, we have
V   3 Σ 33 t δ wdV + V   3 Σ 13 t x 3 δ w dV + V   1 Σ 13 t δ wdV 0 L 2 b Σ 33 x 3 = c x 3 = c   δ wdx 1 0 L 1   A   3 Σ 31 t x 3 dS   δ w dx 1 + 0 L 1 A   3 Σ 31 t x 3 dS δ wdx 1 + 0 L A   1 Σ 13 t δ w   dS δ x 1 = 0 L 2 b Σ 33 x 3 = c x 3 = c   δ wdx 1 + V ¯ r   δ w   dx 1 x 1 = 0 x 1 = L + 0 L 1 A   3 Σ 13 γ w ¨ x 3 dS δ wdx 1 + 0 L 1 A   Σ 13 + γ w ¨ dS δ wdx 1 = 0 L 2 b Σ 33 x 3 = c x 3 = c   δ wdx 1 + V ¯ r   δ w   dx 1 x 1 = 0 x 1 = L + 0 L 1 A   3   Σ 13 x 3 dS δ wdx 1 + 0 L 1 A   Σ 13 dS δ wdx 1 + 0 L γ A   w ¨ δ wdx 1 = V ¯ r δ w x 1 = 0 x 1 = L + 0 L 2 b Σ 33 x 3 = c x 3 = c + 2 b 1 Σ 13 x 3 = c x 3 = c + γ A   w ¨ δ wdx 1 = V ¯ r δ w x 1 = 0 x 1 = L + 0 L q ¯ r δ wdx 1 .
Keeping in mind the relations (56) for the displacement components, the last integral in Equation (A1), including acceleration terms, becomes
V   ρ u ¨ 1 δ u 1 + u ¨ 3 δ u 3 dV = 0 L ρ A U ¨ δ U + ρ I w ¨ δ w + ρ A w ¨ δ w dx 1 ,  
where, once more, partial integration has been applied.
Summing the results (A2)–(A6), and collecting terms, we may rewrite Equation (A1) in the following form:
N ¯ r δ U x 1 = 0 x 1 = L + H ¯ r δ U x 1 = 0 x 1 = L + V ¯ r δ w x 1 = 0 x 1 = L + M ¯ r δ w x 1 = 0 x 1 = L + m ¯ r δ w x 1 = 0 x 1 = L + 0 L p ¯ r   δ Udx 1 + 0 L q ¯ r   δ wdx 1 V   E ε 11 δ ε 11 + μ 111 δ   ε 11 dV 0 L ρ A U ¨ δ U + ρ I w ¨ δ w + ρ A w ¨ δ w dx 1 0 L γ A U ¨ δ U + γ I w ¨ δ w dx 1 = 0 .
If μ 111 is replaced by the constitutive law (60), then this equation implies Equation (95).

Appendix B

We observe that Equation 120 can be recast with the aid of Equation (56) as follows
V   1 Σ 11 δ UdV + V   1 Σ 11 x 3 δ w dV + V   3 Σ 13 δ UdV + V   3 Σ 13 x 3 δ w dV + V   1 Σ 13 δ wdV + V   3 Σ 33 δ wdV +   δ W r in = 0 .  
The calculation of the integrals containing stress components is based on the definitions (126)–(128) and is quite similar to corresponding calculations in the last section (cf. also [25]).
The results are
V   1 Σ 11 δ UdV = N r δ U x 1 = 0 x 1 = L + H r δ U x 1 = 0 x 1 = L V   E ε 11 δ U + μ 111 δ U dV V   μ 111 δ U dV ,
V   1 Σ 11 x 3 δ w dV = M r δ w x 1 = 0 x 1 = L + m r δ w x 1 = 0 x 1 = L V   [ E ε 11 x 3 δ w + μ 111 x 3 δ w dV ,
V   3 Σ 13 δ UdV = 0 L p r   δ Udx 1
V   3 Σ 33 δ wdV + V   3 Σ 13 x 3 δ w dV + V   1 Σ 13 δ wdV = V r δ w x 1 = 0 x 1 = L + 0 L q r   δ wdx 1 .  
Substitution of these results into Equation (A8), and collecting terms, leads to
N r δ U x 1 = 0 x 1 = L + H r δ U x 1 = 0 x 1 = L + V r δ w x 1 = 0 x 1 = L + M r δ w x 1 = 0 x 1 = L + m r δ w x 1 = 0 x 1 = L + 0 L p r   δ Udx 1 + 0 L q r   δ wdx 1 + δ W r in = V   E ε 11 δ ε 11 + μ 111 δ ε 11 dV = δ Π i ,
with Π r i being defined in Equation (98).

References

  1. Askes, H.; Aifantis, E.C. Gradient elasticity in statics and dynamics: An overview of formulations, length scale identification procedures, finite element implementations and new results. Int. J. Solids Struct. 2011, 48, 1962–1990. [Google Scholar] [CrossRef]
  2. Russo, R.; Girot Mata, F.A.; Forest, S.; Jacquin, D. A review on strain gradient plasticity approaches in simulation of manufacturing processes. J. Manuf. Mater. Process. 2020, 4, 87. [Google Scholar] [CrossRef]
  3. Zhuang, X.; Zhou, S.; Huynh, G.D.; Areias, P.; Rabczuk, T. Phase field modeling and computer implementation: A review. Eng. Fract. Mech. 2022, 262, 108234. [Google Scholar] [CrossRef]
  4. Dunn, J.; Serrin, J. On the thermodynamics of interstitial working. Arch. Ration. Mech. Anal. 1988, 88, 95–133. [Google Scholar] [CrossRef]
  5. Eringen, A.C. Theory of Micropolar Elasticity; Springer: New York, NY, USA, 1999. [Google Scholar]
  6. Eringen, A.C. Microcontinuum Field Theories: I. Foundations and Solids; Springer Science & Business Media: New York, NY, USA, 2012. [Google Scholar]
  7. Mindlin, R.D. Microstructure in linear elasticity. Arch. Ration. Mech. Anal. 1964, 16, 51–78. [Google Scholar] [CrossRef]
  8. Toupin, R. Elastic materials with couple-stresses. Arch. Ration. Mech. Anal. 1962, 11, 385–414. [Google Scholar] [CrossRef]
  9. Lam, D.C.; Yang, F.; Chong, A.C.; Wang, J.; Tong, P. Experiments and theory in strain gradient elasticity. J. Mech. Phys. Solids 2003, 51, 1477–1508. [Google Scholar] [CrossRef]
  10. Lei, J.; He, Y.; Guo, S.; Li, Z.; Liu, D. Size-dependent vibration of nickel cantilever microbeams: Experiment and gradient elasticity. AIP Adv. 2016, 6, 105202. [Google Scholar] [CrossRef]
  11. dell’Isola, F.; Seppecher, P.; Spagnuolo, M.; Barchiesi, E.; Hild, F.; Lekszycki, T.; Giorgio, I.; Placidi, L.; Andreaus, U.; Cuomo, M.; et al. Advances in pantographic structures: Design, manufacturing, models, experiments and image analyses. Contin. Mech. Thermodyn. 2019, 31, 1231–1282. [Google Scholar] [CrossRef]
  12. De Domenico, D.; Askes, H.; Aifantis, E.C. Capturing wave dispersion in heterogeneous and microstructured materials through a three-length-scale gradient elasticity formulation. J. Mech. Behav. Mater. 2018, 27, 20182002. [Google Scholar] [CrossRef]
  13. Chakraborty, A. Prediction of negative dispersion by a nonlocal poroelastic theory. J. Acoust. Soc. Am. 2008, 123, 56–67. [Google Scholar] [CrossRef] [PubMed]
  14. Tong, L.; Yu, Y.; Hu, W.; Shi, Y.; Xu, C. On wave propagation characteristics in fluid saturated porous materials by a nonlocal Biot theory. J. Sound Vib. 2016, 379, 106–118. [Google Scholar] [CrossRef]
  15. Broese, C.; Tsakmakis, C.; Beskos, D. Mindlin’s micro-structural and gradient elasticity theories and their thermodynamics. J. Elast. 2016, 125, 87–132, Erratum in J. Elast. 2016, 125, 133–136. [Google Scholar] [CrossRef]
  16. Bedford, A. Hamilton’s Principle in Continuum Mechanics; Pitman Advanced Publishing Program; Springer Nature: Boston, MA, USA; London, UK; Melbourne, Australia, 1985. [Google Scholar]
  17. Gurtin, M.E.; Fried, E.; Anand, L. The Mechanics and Thermodynamics of Continua; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  18. Liu, I.S. Constitutive theories: Basic principles. In Encyclopedia of Life Support Systems (EOLSS); Continuum Mechanics; United Nations Educational, Scientific and Cultural Organization: London, UK, 2009; Volume 1, p. 198. [Google Scholar]
  19. Liu, I.S.; Sampaio, R. On objectivity and the principle of material frame-indifference. Mec. Comput. 2012, 31, 1553–1569. [Google Scholar] [CrossRef]
  20. Leigh, D.C. Nonlinear Continuum Mechanics: An Introduction to the Continuum Physics and Mathematical Theory of the Nonlinear Mechanical Behavior of Materials; McGraw-Hill: New York, NY, USA, 1968. [Google Scholar]
  21. Altan, B.S.; Aifantis, E.C. On some aspects in the special theory of gradient elasticity. J. Mech. Behav. Mater. 1997, 8, 231–282. [Google Scholar] [CrossRef]
  22. Broese, C.; Tsakmakis, C.; Beskos, D. Gradient elasticity based on Laplacians of stress and strain. J. Elast. 2018, 131, 39–74. [Google Scholar] [CrossRef]
  23. Bauchau, O.A.; Craig, J.I. Structural Analysis: With Applications to Aerospace Structures; Springer Science & Business Media: New York, NY, USA, 2009. [Google Scholar]
  24. Sideris, S.A.; Tsakmakis, C. Consistent Euler–Bernoulli beam theories in statics for classical and explicit gradient elasticities. Compos. Struct. 2022, 282, 115026. [Google Scholar] [CrossRef]
  25. Bröse, C.; Sideris, S.A.; Tsakmakis, C.; Üngör, Ö. Equivalent Formulations of Euler–Bernoulli Beam Theory for a Simple Gradient Elasticity Law. J. Eng. Mech. 2023, 149, 04022081. [Google Scholar] [CrossRef]
  26. Lazopoulos, K.A.; Lazopoulos, A.K. Bending and buckling of thin strain gradient elastic beams. Eur. J. Mech. A Solids 2010, 29, 837–843. [Google Scholar] [CrossRef]
  27. Papargyri-Beskou, S.; Tsepoura, K.G.; Polyzos, D.; Beskos, D. Bending and stability analysis of gradient elastic beams. Int. J. Solids Struct. 2003, 40, 385–400. [Google Scholar] [CrossRef]
Figure 1. Rectangular beam of length L, width 2b and height 2c.
Figure 1. Rectangular beam of length L, width 2b and height 2c.
Materials 17 01760 g001
Figure 2. Force controlled loading. Model responses according to Version 1 for various values of l ˜ with w ˜ = 1.5, F ˜ A = 5 · 10 4   ; distributions of (a) U ˜ 0 and (b) N ˜ 0 .
Figure 2. Force controlled loading. Model responses according to Version 1 for various values of l ˜ with w ˜ = 1.5, F ˜ A = 5 · 10 4   ; distributions of (a) U ˜ 0 and (b) N ˜ 0 .
Materials 17 01760 g002
Figure 3. Force controlled loading. Model responses according to Version 2 (top) and Version 3 (bottom) for various values of γ ˜ with w ˜ = 1.5, F ˜ A = 5 · 10 4 and l ˜ 2 = 0.05 ; (a,c) distributions of U ˜ 0 and (b,d) distributions of N ˜ 0 .
Figure 3. Force controlled loading. Model responses according to Version 2 (top) and Version 3 (bottom) for various values of γ ˜ with w ˜ = 1.5, F ˜ A = 5 · 10 4 and l ˜ 2 = 0.05 ; (a,c) distributions of U ˜ 0 and (b,d) distributions of N ˜ 0 .
Materials 17 01760 g003
Figure 4. Displacement controlled loading. Model responses according to Version 1 for various values of l ˜ , with w ˜ = 1.5 , U ˜ A = 5 · 10 3 ; (a) distributions of U ˜ 0 and (b) distributions of N ˜ 0 .
Figure 4. Displacement controlled loading. Model responses according to Version 1 for various values of l ˜ , with w ˜ = 1.5 , U ˜ A = 5 · 10 3 ; (a) distributions of U ˜ 0 and (b) distributions of N ˜ 0 .
Materials 17 01760 g004
Figure 5. Displacement controlled loading. Identical U ˜ 0 -distributions for Versions 2 and 3; w ˜ = 1.5, U ˜ A = 5 · 10 3 , l ˜ 2 = 0.04 .
Figure 5. Displacement controlled loading. Identical U ˜ 0 -distributions for Versions 2 and 3; w ˜ = 1.5, U ˜ A = 5 · 10 3 , l ˜ 2 = 0.04 .
Materials 17 01760 g005
Figure 6. Displacement controlled loading. Predicted N ˜ 0 -distributions for Version 2 (a) and Version 3 (b); w ˜ = 1.5, U ˜ A = 5 · 10 3 , l ˜ 2 = 0.04 .
Figure 6. Displacement controlled loading. Predicted N ˜ 0 -distributions for Version 2 (a) and Version 3 (b); w ˜ = 1.5, U ˜ A = 5 · 10 3 , l ˜ 2 = 0.04 .
Materials 17 01760 g006
Figure 7. Force control loading. Model responses according to Version 1 for various values of l ˜ and w ˜ = 0.05, F ˜ A = 5 · 10 9 . Distributions of (a) w ˜ 0 and (b) V ˜ 0 .
Figure 7. Force control loading. Model responses according to Version 1 for various values of l ˜ and w ˜ = 0.05, F ˜ A = 5 · 10 9 . Distributions of (a) w ˜ 0 and (b) V ˜ 0 .
Materials 17 01760 g007
Figure 8. Force controlled loading. Model responses according to Version 2 (top) and Version 3 (bottom) for various values of γ ˜ and w ˜ = 0.05, F ˜ A = 5 · 10 9 , l ˜ = 0.05 ; (a,c) distributions of w ˜ 0 and (b,d) distributions of V ˜ 0 .
Figure 8. Force controlled loading. Model responses according to Version 2 (top) and Version 3 (bottom) for various values of γ ˜ and w ˜ = 0.05, F ˜ A = 5 · 10 9 , l ˜ = 0.05 ; (a,c) distributions of w ˜ 0 and (b,d) distributions of V ˜ 0 .
Materials 17 01760 g008
Figure 9. Deflection controlled loading. Model responses according to Version 1 for various values of l ˜ and w ˜ = 0.05, w ˜ A = 5 · 10 3 ; distributions of (a) w ˜ 0 and (b) V ˜ 0 .
Figure 9. Deflection controlled loading. Model responses according to Version 1 for various values of l ˜ and w ˜ = 0.05, w ˜ A = 5 · 10 3 ; distributions of (a) w ˜ 0 and (b) V ˜ 0 .
Materials 17 01760 g009
Figure 10. Displacement controlled loading. Identical w ˜ 0 -distributions for Versions 2 and 3; w ˜ = 0.05, U ˜ A = 5 · 10 3 , l ˜ 2 = 0.04 .
Figure 10. Displacement controlled loading. Identical w ˜ 0 -distributions for Versions 2 and 3; w ˜ = 0.05, U ˜ A = 5 · 10 3 , l ˜ 2 = 0.04 .
Materials 17 01760 g010
Figure 11. Deflection controlled loading. Predicted V ˜ 0 -distributions for Version 2 (a) and Version 3 (b); w ˜ = 0.05, U ˜ A = 5 · 10 3 , l ˜ 2 = 0.04 .
Figure 11. Deflection controlled loading. Predicted V ˜ 0 -distributions for Version 2 (a) and Version 3 (b); w ˜ = 0.05, U ˜ A = 5 · 10 3 , l ˜ 2 = 0.04 .
Materials 17 01760 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tsakmakis, C.; Broese, C.; Sideris, S.A. Dynamics in Explicit Gradient Elasticity: Material Frame-Indifference, Boundary Conditions and Consistent Euler–Bernoulli Beam Theory. Materials 2024, 17, 1760. https://doi.org/10.3390/ma17081760

AMA Style

Tsakmakis C, Broese C, Sideris SA. Dynamics in Explicit Gradient Elasticity: Material Frame-Indifference, Boundary Conditions and Consistent Euler–Bernoulli Beam Theory. Materials. 2024; 17(8):1760. https://doi.org/10.3390/ma17081760

Chicago/Turabian Style

Tsakmakis, Charalampos, Carsten Broese, and Stergios Alexandros Sideris. 2024. "Dynamics in Explicit Gradient Elasticity: Material Frame-Indifference, Boundary Conditions and Consistent Euler–Bernoulli Beam Theory" Materials 17, no. 8: 1760. https://doi.org/10.3390/ma17081760

APA Style

Tsakmakis, C., Broese, C., & Sideris, S. A. (2024). Dynamics in Explicit Gradient Elasticity: Material Frame-Indifference, Boundary Conditions and Consistent Euler–Bernoulli Beam Theory. Materials, 17(8), 1760. https://doi.org/10.3390/ma17081760

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop