1. Introduction
It is well known that a structure consisting of a linear, isotropic, classical elastic material, subjected to external time-invariant boundary conditions, behaves according to the Navier–Cauchy equilibrium equation [
1,
2,
3], which in terms of displacements
and without the presence of body forces reads:
where
is the gradient operator, while
and
indicate the Lamé constants.
Despite the discrete nature of real materials, the continuum theory of classical elasticity—described by (1) for the elastostatics case—is deduced by considering that the dimensions of the material microstructure are much smaller than a material representative volume element (RVE), which in turn is much smaller than any dimension of the loaded structure. Additionally, the material properties and the generated elastic fields in the RVE are projected, by averaging, around a point
lying at the center of the RVE. That projection imposes the local nature of the classical theory of elasticity and requires that displacements, stresses, and strains vary constantly or linearly throughout the material RVE [
4,
5]. Obviously, the situation becomes problematic when the material inhomogeneity is comparable with the size of the structure and the averaging performed in the RVE requires the consideration of strain gradients in the potential energy density and the introduction of internal length scale parameters, which are able to capture size effect phenomena.
At the beginning of 20th century, the Cosserat brothers [
6] proposed the idea of an enhanced elastic theory in which, except strains and stresses, the gradient of rotations and the dual in energy, couple stresses should be considered. Their idea reached maturity almost fifty years later with the general works of Toupin [
7], Mindlin and Tiersten [
8], Green and Rivlin [
9] and Koiter [
10]. Meanwhile Mindlin [
11], while investigating the influence of couple stresses to stress concentrations mentions that “
Also, it would seem to be desirable to explore the consequences of taking into account the remaining components of the strain gradient and, perhaps, second and higher gradients of the strain”. Indeed, one year later, Mindlin [
12] published his general dynamic theory for elastic materials with microstructure, where the microstructure is considered as an additional micro-continuum embedded at every point of the macro-continuum, thus justifying, the presence of higher order strain gradients in the expressions of potential energy density. Ignoring inertia terms, the couple stress theory of Toupin [
7], Mindlin’s elastic theory with microstructure [
12] and later the virtual power theory of Germain [
13] lead, for a material with microstructural effects, to the same equilibrium equation and boundary conditions. However, the elastostatic Form II version of Mindlin’s theory has the very attractive characteristic of retaining the symmetry of the considered total stress tensor as in the classical elasticity case, and concludes in a simple equilibrium equation of the following form:
where
,
are internal length scale parameters that can facilitate microstructural effects for dilatational and shear deformations, respectively.
Considering uniform microstructural effects for all types of deformation, i.e.,
, Equation (2) can be further simplified obtaining the form:
Equation (3) and the corresponding boundary conditions of Mindlin’s Form II theory consist of an attractive enhanced elastic theory with a microstructure, known as strain gradient elasticity (SGE) or dipolar gradient elasticity (DGE), since it employs only one internal length scale parameter in addition to the two classical Lamé constants, and the most important strains and stresses, appearing in its constitutive equations, are symmetric as in the classical elastic case. During the last thirty years, many authors exploited the simplicity of SGE to solve analytically and/or numerically elastostatic problems with microstructural effects in many fields of linear elastic continuum mechanics, such as fracture and dislocations mechanics [
14,
15,
16,
17,
18,
19,
20] and structural and material response [
21,
22,
23,
24,
25,
26,
27,
28,
29,
30], while interesting remarks on SGE can be found in [
5,
31,
32,
33,
34].
Mindlin [
11] first proposed a solution representation of Equation (2) based on Papkovich–Neuber type vector and scalar potentials
,
[
3,
35], which, in the case of Equation (3), is simplified to [
30]:
Instead of (4), Charalambopoulos and Polyzos [
26] utilized the following representation of the solution of (3):
without providing any correlation of (5) with the solution representation (4). The proof of the decomposition (5) is provided in this paper through the proof of theorem 1 in
Section 2 and via Papkovich–Neuber type potentials in
Appendix A. It should be mentioned at this point that the Papkovich–Neuber type gradient elastic solution (4) agrees with the corresponding one provided by Solyaev et al. [
36,
37], while the solution decomposition (5) is in agreement with the decomposition proposed by Lazar [
38] if one considers that for a classical elastic solution
satisfying Equation (1), the vector function
remains a classical elastic solution. In [
34], Gourgiotis et al. solve a sharp notch problem in microstructured solids utilizing a Knein–Williams technique, and mentioned that their asymptotic solution shows significant departure from those of classical elasticity. This statement is basically the motivation for the present work, which, among others, proposes a solution representation of (3) by comprising the solution of the corresponding classical elastic boundary value problem, which is absent in both representations (4) and (5).
The first attempt at incorporating the classical elastic solution in the solution of (3) is that of Ru and Aifantis [
39]. More specifically, considering the elastic displacement field
that satisfies (1) and the corresponding classical elastic boundary conditions in a domain
confined by a surface
, they proposed as strain gradient elastic solution of (3) the vector
that satisfies the non-homogeneous partial differential equation:
and the extra boundary condition
where
denotes differentiation with respect to the unit normal vector of
.
Comparing (3) with (6), the obvious advantage of this representation is the reduction in the order of the partial differential equation of the problem by two. However, as it is mentioned in Charalambopoulos and Polyzos [
26], the representation (6) and (7) is questionable for the following reasons: (i) the solution of the partial differential equation of second order—as Equation (6)—satisfies a boundary condition of second order as the condition (7) is in contradiction with the mathematically accepted condition that the maximal degree of boundary conditions never exceeds the crucial number
n − 1, where
n stands for the degree of the differential equation. (ii) The boundary condition (7) is arbitrary and not the outcome of a variational process. (iii) The classical solution
satisfies the equilibrium Equation (3), but not Equation (6).
Charalambopoulos et al. [
30] utilized a representation of the solution that satisfies Equation (3) and the corresponding Form II boundary conditions and comprises the corresponding classical elastic solution, without, however, providing any systematic proof of that representation. This is, among others, the goal of the present work. The possible application of the presented here methodology to Equation (1) and to equations describing the behavior of linear pantographic sheets [
40] or obeying to the generalized Hook’s law for isotropic second gradient materials [
41] will be the subject of future work. The structure of the present work is the following: The next section illustrates the Form II SGE theory of Mindlin with only one internal length scale parameter.
Section 3 is entirely devoted to the mathematical establishment of a solution representation of (3), which encompasses the corresponding classical solution and its convergence behavior as the gradient parameter tends to zero. The same solution representation is exploited in
Section 4 to show that the bending stiffness of a Form II SGE Bernoulli–Euler beam depends on the material rigidity
and the internal length scale parameter
and not on
,
plus the area of the cross-section of the beam.
2. Strain Gradient Elastostatics with One Internal Length Scale Parameter
The present section reports a boundary value problem in terms of the simplest possible strain gradient elastostatic theory with one intrinsic parameter. Mindlin [
42] in the second version of his theory considered that the first gradient elastic potential energy density for an elastic body with microstructure is a quadratic form of the strains
and the gradient of strains
. For isotropic materials, this theory provides a potential energy density containing the two Lamé material constants and five constants that normalize the terms of strain gradients, i.e.,
where
denotes spatial differentiation,
are the displacement components,
,
are the well-known Lamé constants having units N/m
2 and
are five constants having units of force, all explicitly provided in [
12].
For the case of
and
,
, the potential energy density
obtains the form
where
is the only internal length scale parameter that correlates the microstructure with macrostructure, having units of length (m).
Strains and gradient of strains are dual in energy with the Cauchy-like stresses and double stresses, respectively, defined as:
and
If the Young modulus and Poisson ratio are used instead of the Lame constants , , then the replacements , should be made.
Considering a material with a microstructure of volume
and external boundary
, the variation of the total potential energy (9) provides, after some algebra [
12,
43,
44], the following relation:
where the vectors
,
represent the traction and double traction vectors, respectively, defined on the boundary
S and written as
and
where
and
represent the tangential and normal gradient operators on
, respectively, and have the form
The vector
in Equation (12) concerns non-smooth boundaries with at least one corner
c in two dimensions or at least one closed edge line
in three dimensions, admitting the form:
where
denotes the difference of
at both sides of the corner
c or the edge
, while
stands for the tangential vector in both sides of a corner or edge.
Equilibrating (12) with the variation of the work performed by external body force
, boundary tractions
, double traction
and jump traction
, we arrive at the following equilibrium equation:
accompanied by the classical essential and natural boundary conditions where the displacement vector
and/or the traction vector
must be defined on the global boundary
, i.e.,
and the non-classical essential and natural boundary conditions where the normal displacement vector
, the double traction vector
or jump traction
are prescribed on
, i.e.,
In terms of the displacement vector
and free of body forces, Equation (17) obtains the form
Theorem 1. The solution of Equation (20) can be written aswithand.
Proof. We denote as
the classical elastic elliptic differential operator
Given that
, we infer that
and
Equations (22) and (24) imply that
We consider the Helmholtz decomposition of the field
. Then, Equation (22) leads to
where
, is a harmonic function. Then, Equations (25) and (26) give
Additionally, it holds that
and therefore
, where
is another harmonic function. One partial solution of this equation is clearly exactly the function
and so the general solution is
, with
. Then, (27) is rewritten as
The first two terms of the decomposition (28) form the classical part while the remaining terms form the component obeying to the homogeneous modified Helmholtz equation. □
3. On Representing Strain Gradient Elastic Solutions via the Solution of the Corresponding Classical Elastic Boundary Value Problem
All the representations of the solution of a strain gradient elastostatic problem appearing in the literature, do not include as a constituent the respective classical elastic solution, while the convergence behavior of these solutions when the gradient parameters fade away is not studied. As in [
45,
46], in the following, we present a general result considering all the above concerns.
Consider a bounded open region
,
, with its boundary
being a Lipschitz surface. Assuming that the body forces are absent, we consider the solution of the boundary value problem consisting of the fourth order partial differential Equation (20) and a set of classical and non-classical boundary conditions as those illustrated in
Section 2. More precisely, we partition the surface
twice, first in two subdomains
and
(with
), where classical conditions are imposed, and secondly, in the subdomains
and
, whose common boundary Γ is of dimension
and represents the corners
or the edges
of the surface
. Then, the general set of mixed-type classical conditions can be formulated as
along with the set of non-classical conditions
where
,
and
.
The solution and the fields , , depend of course on the parameter , but this is omitted for simplicity and will be notified only when it is necessary.
It is assumed that the given functions
,
,
,
and
share all the required regularity for the well-posedness of the traces of the solution of Equation (20)—and its derivatives—on
. To clarify rigorously the last remark, we could additionally invoke the functional theoretic framework of variational problems settled in Sobolev spaces. This approach is described extensively in Wloka [
47] for the general case of elliptic boundary value problems of a higher order and is very profitable since it provides results with crucial influence in the solvability of the problem under consideration with existence, uniqueness, and stability. It is out of the scope of the present work to give an extensive investigation via the aforementioned alternative framework, but it would be very helpful to give some brief concepts, facilitating the comprehension of the inner structure of the boundary value problems under investigation. Then, we recall the Sobolev spaces
and
, which are complete Hilbert spaces built by functions (or distributions) with specific integration behavior over
and
. Among the more recognizable Sobolev spaces, we encounter the space of square integrable measurable functions
and its subspace
, whose elements possess distributional derivatives with the same square integrable behavior. Every Hilbert space
, with
ℕ, has a usual inner product and an induced norm consisting of the
norms of the derivatives up to degree
. Therefore, the inclusion
when
is an obvious relation. For real and positive order
, the norm is defined in a more complicated manner, but generalizes naturally what happens for an integer order. The linear spaces inclusion above still holds. When
, the elements of the space
belong to
, but this is not the case when
and constitute then pure distributions without functional representative.
Although not necessarily classical functions, the elements of have trace on the boundary , which are generalized functions belonging to . The distributional surface normal derivative (of order k) of an element in belongs to in the case that the smoothness of the surface allows the induced differentiability. The space is a subspace of , with elements that, along with all their normal derivatives of order less than , have zero traces on .
Notice at this point that, for the case of the closed surface , the space is the dual space of and so surface terms of the form naturally arise, where and . This term expresses the action of on and it is very reminiscent of the virtual work of a force over a displacement. Only when , we have the reduction in the dual pairing to the usual inner product .
After this brief discussion, we notice that in the framework of a boundary value problem whose differential equation is of order
, the boundary Equations (29) and (30) involve
boundary operators (of possibly mixed type) (when
). These operators are characterized by their own orders, which obey to the rule:
. If, for example, we had exactly the boundary conditions (29) and (30) valid on the whole surface
(with
), the involved boundary operators
would be
(with order
) and
(with order
). This example corresponds to the so called Dirichlet boundary value problem of fourth order. This settlement lies totally in the general framework of the abstract theory of elliptic boundary value problems: The main outcome is the existence and uniqueness of the solution of the problem (20), (29) and (30) with the following regularity general result [
47]
and so, assigning to boundary operators their exact orders, we obtain
We would like to note the double consecutive dimension reduction
, which takes place in the treatment of the fifth term of the r.h.s of the last equation. The regularity of the solutions depends on the regularity of the data in a specific manner. It is not explicitly apparent, but when the parameter
increases, the validity of Equation (31) passes through further assumptions for additional smoothness of the boundary
. The most encountered evocation of the representation (31) is with the selection
(the order of the differential equation). Therefore, for the solution
to have square integrable (in
) derivatives up to fourth order, all the data must belong to “smooth” Sobolev spaces of positive order. Then, the differential equation is satisfied in
-sense. However, this is accompanied with the hypothesis of a smooth
boundary. To permit a boundary with corners or edges and to obtain the broader class of admissible solutions, it is preferable to work with
. Then,
In this case, the data
,
loose regularity, but still belong to the realm of square integrable functions. However, the data
,
, and
pertaining to stresses and jumps of double stresses over corners (edges) might become not square integrable distributions, ready to act—via the mentioned above dual pairings—on their reciprocal fields. In this point, we would like to say that our analysis has been facilitated from the fact that we are in absence of body forces. The differential equation is not any more valid classically, but in the distributional sense. It is noticeable that due to the special coercivity behavior of the bilinear form corresponding to the gradient elasticity operator, the generic constant
appearing in relation (32) cannot present worse behavior than the asymptotic convergence
for
. In addition, (32) implies a fortiori the boundedness
It is noticeable here that the constant is independent of since suppressing downwards the energy bilinear form and keeping only norms of first derivatives involve exclusively the Lamé constants of classical elasticity.
After the brief introduction of the functional theoretic setting, we are in position to present the two main accomplishments of the current work. First, we are going to present the construction of a very useful decomposition of the unique solution of the problem under consideration. In the sequel, we will state the necessary assumptions on the data so that this representation obtains a stable (with respect to ) behavior, incorporating appropriately the classical solution.
Consider the auxiliary second order classical boundary value problem (titled Problem I), satisfied by the solution
:
where
, is the classical surface traction field.
The above defined problem involving Equations (34)–(36) used the classical data from the gradient problem and “ignores” the non-classical ones. We emphasize that in Problem I, the gradient boundary term
offers its place to its classical counterpart
. Applying the classical well-known framework of the preceding analysis concerning this time the traditional second order elliptic boundary value problems, we deduce easily that the unique classical solution satisfies the stability relation
It is worthwhile to mention here that the last equation defines exactly the needed regularity of the data , for the stated stability to be guaranteed. However, the gradient problem has already assigned specific regularity assumptions on the data. Compromising the two groups of requirements, we deduce that as well as .
We are in position to state the main representation theorem of this work. We set first, as induced by the discussion above, the broader possible space in which the data are permitted to belong:
The following theorem holds:
Theorem 2. Let the boundary value problem consist of Equations (20), (29) and (30) where the data . This problem is a well-posed fourth order elliptic boundary value problem with a unique solution. This solution can be represented as followswheresatisfies Problem I, satisfies the classical elastostatic equation and obeys to the modified Helmholtz equation: Proof. As explained above, the problem (20), (29) and (30) disposes a unique solution
, which is stable with respect to the data as Equations (32) and (33) guarantee. We consider the decomposition
The function
satisfies Equation (20) given that both
and the field
obey to this equation too. Based on the Papkovich-type representation [
26,
30], the field
can be written as
where
Working with Equation (41), we have that necessarily
is a harmonic function
. Then,
where
satisfies the modified Helmholtz equation
.
A similar treatment applies to Equation (42) and the result is that the fields of Equations (41) and (42) have the decomposed form
Inserting these representations in Equation (40) leads to the decomposition
where
and
It is evident that
. In addition, it comes out easily that
belongs to the kernel of the operator
. Consequently, we deduce that
from where we infer that
. □
As a conclusion, the decomposition (37) with the differential properties (38) and (39) was derived from the point of view of the underlying differential equations and always can be applied to the unique solution of the fourth order boundary value problem under discussion. The implication of boundary conditions is of course the next step. The representation (37) has the advantage that it disposes additional degrees of freedom since the involved “free” functions satisfy the second order differential equations in which the gradient elasticity law decomposes. In addition, the representation (37) involves the classical solution of the problem in the absence of the microstructure as a cornerstone constituent. It is very interesting to examine whether taking the limit of the expression (37) as leads to the classical solution. What matters of course is primary the construction of the solution of the boundary value problem independently of its convergence behavior. Nevertheless, it is essential to construct sufficient conditions on the data assuring this desirable convergence property.
The next theorem verifies that under specific assumptions on the data, convergence is established. The first requirement is in accordance with the underlying constitutive equations concerning double stresses. Indeed, it is natural that the magnitude of the double stresses and the relevant jump fields imposed on the structure obey to specific order analysis with respect to the microstructure parameter. More precisely, the boundary tensors
and
are selected as follows:
and as stated, this choice has a physical origin. The additional requirements are related to the necessary regularity of the rest of the data. It is necessary, in the convergence setting for the boundary data
, to be a genuine square integrable function, in fact an element of
and that
remains bounded as
varies. In practice, the term
is usually independent of
.
As far as the stated above restriction of is concerned, the concept stems from the treatment of the classical solution. As we will see, under this assumption, the function becomes a square integrable function with a crucial role in the convergence analysis.
The following theorem holds.
Theorem 3. Let the boundary value problem consist of Equations (20), (29) and (30). The dataare elements of B with the restriction that, obey to Equations (48) and (49),andis a bounded function of . Two mutually exclusive alternatives arise:
- (A)
- (B)
In case (B), we impose further regularity on the data by demandingand.
Let be the unique solution of the gradient elasticity boundary value problem constructed in Theorem 2. Then, the following asymptotic analysis holds:rendering Problem I, the asymptotic classical limit of the gradient problem, when the impact of the microstructure fades away. Proof. The classical elastic displacement field satisfies the following stability condition:
as stated before. However, due to the additional regularity
(valid throughout the work) and the extra requirement
introduced in the assumptions of the current theorem, we are in position to invoke the well-known regularity theory for second order elliptic boundary value problems guaranteeing that
Consequently, the classical solution has square integrable derivatives of the second order and the theory of traces on the boundary implies that all the terms
belong to
as well as .
On the basis of Equation (7) and via tensor symbolism—just for condensing the form of the expressions—the following Betti’s form type result, referring to the elastic field
, can be constructed:
Here, we encounter the dual pairings , and between dual spaces that represent surface and curve virtual actions. Only when all the involving fields are regular enough (square integrable functions) these terms give place to the well- known surface inner products.
On the basis of decomposition (45), the boundary conditions satisfied by the two partners of this decomposition and the splitting imposed by Equation (19), we remark that
where
(see Equation (13)).
Handling the non-classical boundary conditions of
leads to
and
Then, Equation (53) obtains the form:
where
Thanks to the introductory discussion of this theorem, pertaining to the regularity of surface terms generated by the classical field, it is clear that all the terms participating in are well defined and belong to . In addition, given that in our elliptic boundary value problem, every surface norm of traces of the classical field and its derivatives is controlled by the volume norm , which is bounded by the data as Equation (52) implies, the surface field is uniformly bounded in with magnitude of order .
The fields and , on surfaces and (or) curves have appropriate norms, which are bounded by the norm of the solution (due to the continuity of the traces with respect to the solution of the boundary value problem), which is, by its turn, bounded uniformly with respect to the gradient parameter , via an estimate of the form (33) (applied to ).
If the case (A) holds, the second dual pairing in r.h.s. of Equation (56) disappears and no need to handle the functional on arises. When we have case (B), we impose further regularity on the displacement and stresses leading to higher regularity on the classical solution. This regularity is optimally selected to guarantee that the crucial field belongs to (since under this choice, the norm is kept bounded).
Finally, the field on equals with the field , which incorporates second order derivatives of . Then, the adequate surface —norm of is again uniformly bounded (due to the continuity of the traces with respect to the solution of the boundary value problem) by , which via (33) is bounded by the norms of the data multiplied by . This coefficient goes to 0 as , since, as mentioned before, the generic constant does not depend on .
In any case, taking the limit as
, in expression (56), all the terms of the right-hand side converge uniformly to 0, therefore
The positivity of the bilinear elastic form guarantees that
, which implies that
is a constant. Given that the trace
, we infer that
from where we obtain immediately the asymptotic behavior (50). □
Much effort was places in finding the necessary regularity of the data assuring the desired convergence when microstructure behavior disappears. It is interesting that this is necessarily valid when the data are smooth analytic functions of their arguments. Indeed, we have the following corollary, which would be the main outcome of this work if the data were considered analytic, but in our opinion, this could not be proven without the herein adopted generalization to abstract functional spaces.
Corollary 1. Consider that the boundary value problem consists of Equations (20), (29) and (30). The data are analytic functions of their arguments with the restriction that,obey to Equations (48) and (49). Letbe the unique solution of the gradient elasticity boundary value problem, expressed on the basis of the decomposition constructed in Theorem 2. Then, the following asymptotic result is obtained: Proof. The assumed analytic smoothness of the data confirms that the assumptions of the alternative (B) of Theorem 3 are always valid and then the convergence outcome holds in any case. □
Remark 1. As far as the field on is concerned, we could easily find that the surface —norm of is again uniformly bounded by . On the basis of (24a), this term is bounded by the norms of the data multiplied by . Therefore, if no extra assumption on data was made, the limit would be ambiguous.
Remark 2. We would like to mention once again here that, despite the convergence, which needs some special settlement of the nature of the imposed data, the decomposition (37) is always valid, and when the conditions are arbitrary, what is exactly stated for the participants of the decomposition is that they satisfy the corresponding differential equations, while their superposition satisfies the boundary value problem.
A number of simple one-dimensional problems can reveal the essence of the previous theorems. Let us consider, for example, the elastostatic response of a finite length string. The considered differential equation is given by , accompanied by the classical boundary conditions , and the non-classical ones and .
The solution is
while the classical elastic solution that satisfies the conditions
and
is
.
The same differential equation with the boundary conditions
,
,
and
has the following solution:
with corresponding classical solution, the function
.
Remark 3. The second example of Remark 2 implies that sometimes it is useful to write Equation (37) in the following formwithas, whereobeys to the same differential regime. This approach consists of a repartition of the partners of the decomposition, which seems to be more flexible in applications. This reordering is realizable since both fieldsandbelong to the kernel of the classical elasticity operator. Remark 4. The result stated by Theorem 3 could be violated if the boundary conditions are set arbitrarily. As an example, consider the same differential equation as in Remark 2, accompanied with the boundary conditions, , and. This problem has the unique solution However, the classical solutionsatisfies the boundary conditions , , but does not satisfy the classical boundary condition. Additionally, the part of the solution containing the displacementdoes not converge to zero asat the endpoint. The critical factor here, is that two classical boundary conditions are prescribed simultaneously at the same part of the boundary.
Remark 5. In the statement of the theorem, we demanded that ). In mathematical terms, this is needed for the uniqueness of the solution. Physically, it reflects the necessity of some anchoring of the structure. However, we would like to say that this condition is not of major importance for the validity of Theorem 3. The modifications needed could be stated mathematically—they include some compatibility conditions for the data—but the essence is very simple: The results are exactly the same modulo motions of a rigid body. Consequently, imposing a-posteriori a condition not allowing rigid motions is sufficient to construct the unique solution of the problem.
4. Revisiting Bending Theories of Strain Gradient Elastic Beams
In this section, the solution representation addressed in
Section 3 is employed to present an answer to contradictions arising by two well-known first-strain gradient elastic models proposed in the literature to describe the strain gradient elastostatic bending behavior of Bernoulli–Euler beams.
During the last two decades, a plethora of papers dealing with the static and dynamic response of Bernoulli–Euler strain gradient elastic beams have appeared in the literature. Most of them are based on variational approaches and, for the elastostatic case, the equilibrium equation they conclude has either the following form [
48,
49,
50,
51,
52]
where
coincides with the neutral axis of the beam,
stands for the material Young modulus,
is the moment of inertia of the beam’s cross-section
,
is the transverse beam deflection,
is the transverse external load, and
is the intrinsic strain gradient elastic parameter, or the form [
53,
54,
55,
56,
57,
58,
59,
60,
61]
The essential difference between these two equations is that the first is derived by considering the strain
and the gradient of strain
in the expression of the potential energy density of the beam, while the second one considers the strain
and the strain gradients
,
with
being the differentiation of
with respect to axis y directed along the thickness of the beam. The result is that the bending stiffness in Equation (60) is the same as that of classical elasticity, while the bending stiffness in (61) depends on the internal length scale parameter
and the cross-section of beam A. The interesting point here is that the same categorization is valid for other works dealing with the bending of strain gradient elastic Timoshenko beams and plates, as well as for experimental and numerical validations on the bending response of strain gradient elastic beams. Here, one can mention the works of Papargyri-Beskou and Beskos [
62], Papargyri-Beskou et al. [
63], Triantafyllou and Giannakopoulos [
64], and Gortsas et al. [
29] for the strain gradient model of Equation (60), and the works of Lazopoulos and Lazopoulos [
65], Khakalo and Niiranen [
57,
58,
60,
66], and Korshunova et al. [
67] for the model of Equation (61).
Since there is a principal difference between the two above-mentioned bending models, the question here is which of them is the correct one. Lurie et al. [
68] and Lurie and Solyaev [
69,
70] proposed elegant answers to that question by proving that Equation (60) is the only correct one for the bending response of a strain gradient elastic Bernoulli–Euler beam. Among others, they mentioned that “…
a formal variational procedure for obtaining the governing equilibrium equations in the beam theories, ignoring boundary conditions on the top and bottom surfaces of the beam leads to an erroneous result of abnormal increasing of the beam normalized bending stiffness with decreasing its thickness. …”. Polizzotto [
31] upholds this argument because the normal derivative of displacements identically vanishes at the beam lateral surface, and thus cannot play any role as a boundary layer. In the present section, we reach the same conclusion under the light of the theorems proved in the previous section.
The starting point of our analysis is the pure bending of a classical elastic, isotropic two-dimensional orthogonal rectangle subjected to pure bending under plane strain conditions, as depicted in
Figure 1. The boundary conditions of the problem are
and
,
and the corresponding solution is provided by Selvadurai [
3] in the following form:
Consider a rectangular 3D plate of length
and cross-section
. The torque at both boundaries of the plate is defined as
with
being the moment of inertia of the cross-section
. Equation (62), valid for the midplane of the plate, obtain the form
while the deflection of the neutral axis is given by
for
, i.e.,
which, for
, is identical to the deflection of a classical elastic Bernoulli–Euler beam subjected to pure bending [
3,
69], i.e.,
Next, we consider the same pure bending problem presented in
Figure 1, however, for a strain gradient elastic material. The classical boundary conditions remain the same as in the classical elastic problem, i.e.,
while in non-classical boundary conditions, we consider the zeroing of double tractions at all the external boundaries, for example:
According to the theorem presented in
Section 3, the strain gradient elastic solution of this problem is written as
Evidently, the solution (69) imposes the following forms for tractions and double tractions, respectively
with
being the unit normal vector of the surface, to which both tractions and double tractions are referred.
Concentrating our attention on the classical part of the solution (69), it is easy to observe that [
30]
and
As it is explained by Charalambopoulos et al. [
30], the two components of the gradient elastic traction and double traction vectors, defined on a surface with unit normal vector
, have the form, respectively,
and
In view of (72)–(75), the classical part of the gradient elastic traction and double traction vectors defined on a surface with unit normal vector
exhibit the following forms, respectively
and
Inserting Equations (76) and (77) into (70) and (71), respectively, and satisfying the boundary conditions (67) and (68) one obtains:
and
Equations (78) and (79) indicate that the parts
,
satisfy a homogeneous system of algebraic equations for arbitrary material properties. Apparently, this leads to the conclusion that
which means that the pure bending of the strain gradient elastic plate presented in
Figure 1 has absolutely the same response as that of the classical elastic one.
Extending this result to the behavior of the neutral axis of the plate for
, we obtain
This result is possible only when Equation (60) is valid. Equation (61) is misleading since the double stress does not contribute to the solution of the problem.