Next Article in Journal
(L,M)-Fuzzy k-Pseudo Metric Space
Next Article in Special Issue
Analysis of Industrial Engineering Students’ Perception after a Multiple Integrals-Based Activity with a Fourth-Year Student
Previous Article in Journal
Optimized Power Supply Rejection Ratio Modeling Technique for Simulation of Automotive Low-Dropout Linear Voltage Regulators
Previous Article in Special Issue
Novel Photovoltaic Empirical Mathematical Model Based on Function Representation of Captured Figures from Commercial Panels Datasheet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Representing Strain Gradient Elastic Solutions of Boundary Value Problems by Encompassing the Classical Elastic Solution

by
Antonios Charalambopoulos
1,
Theodore Gortsas
2 and
Demosthenes Polyzos
2,*
1
School of Applied Mathematics and Physical Sciences, National Technical University of Athens, 15780 Athens, Greece
2
Department of Mechanical Engineering and Aeronautics, University of Patras, 26504 Patras, Greece
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(7), 1152; https://doi.org/10.3390/math10071152
Submission received: 24 February 2022 / Revised: 29 March 2022 / Accepted: 31 March 2022 / Published: 2 April 2022
(This article belongs to the Special Issue Mathematics and Its Applications in Science and Engineering)

Abstract

:
The present work aims to primarily provide a general representation of the solution of the simplified elastostatics version of Mindlin’s Form II first-strain gradient elastic theory, which converges to the solution of the corresponding classical elastic boundary value problem as the intrinsic gradient parameters become zero. Through functional theory considerations, a solution representation of the one-intrinsic-parameter strain gradient elastostatic equation that comprises the classical elastic solution of the corresponding boundary value problem is rigorously provided for the first time. Next, that solution representation is employed to give 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.

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 u c l a s s i c a l and without the presence of body forces reads:
μ 2 u c l a s s i c a l + ( λ + μ ) u c l a s s i c a l = 0 ,
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 x 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:
( λ + 2 μ ) ( 1 l 1 2 2 ) u g r a d i e n t μ ( 1 l 2 2 2 ) × × u g r a d i e n t = 0 ,
where l 1 2 , l 2 2 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., l 1 2 = l 2 2 = g , Equation (2) can be further simplified obtaining the form:
( 1 g 2 2 ) [ μ 2 u g r a d i e n t + ( λ + μ ) u g r a d i e n t ] = 0
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 B , B 0 [3,35], which, in the case of Equation (3), is simplified to [30]:
u g r a d i e n t = B λ + μ 2 ( λ + 2 μ ) [ r ( 1 g 2 2 ) B + B 0 ] ( 1 g 2 2 ) 2 B = 0 ( 1 g 2 2 ) 2 B 0 = 0 .
Instead of (4), Charalambopoulos and Polyzos [26] utilized the following representation of the solution of (3):
u g r a d i e n t = u e + u g μ 2 u e + ( λ + μ ) u e = 0 ( 1 g 2 2 ) u g = 0 .
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 u e satisfying Equation (1), the vector function ( 1 g 2 2 ) u e 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 u c l a s s i c a l 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 u g r a d i e n t that satisfies the non-homogeneous partial differential equation:
( 1 g 2 2 ) u g r a d i e n t = u c l a s s i c a l ,
and the extra boundary condition
2 n 2 u g r a d i e n t = 0 ,
where / n 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 u c l a s s i c a l 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 E I and the internal length scale parameter g and not on E I , g 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 ε i j and the gradient of strains κ i j k . 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.,
U = 1 2 λ ε i i ε j j + μ ε i j ε i j + α ^ 1 κ i i k κ k j j + α ^ 2 κ i j j κ i k k + α ^ 3 κ i i k κ j j k + α ^ 4 κ i j k κ i j k + α ^ 5 κ i j k κ k j i ε i j = 1 2 ( i u j + j u i ) = ε j i κ i j k = i ε j k = 1 2 ( i j u k + i k u j ) = κ i k j ,
where i denotes spatial differentiation, u i are the displacement components, λ , μ are the well-known Lamé constants having units N/m2 and α ^ 1 ÷ α ^ 5 are five constants having units of force, all explicitly provided in [12].
For the case of α ^ 1 = α ^ 3 = α ^ 5 = 0 and α ^ 2 = λ g 2 , α ^ 4 = μ g 2 , the potential energy density U obtains the form
U = 1 2 λ ε i i ε j j + μ ε i j ε i j + 1 2 λ g 2 κ i j j κ i k k    + μ g 2 κ i j k κ i j k ,
where g 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:
τ i j = W ε i j = 2 μ ε i j + λ ε k k δ i j ,
and
μ i j k = W κ i j k = g 2 i τ j k
If the Young modulus E and Poisson ratio v are used instead of the Lame constants λ , μ , then the replacements λ = E v ( 1 + v ) ( 1 2 v ) , μ = E 2 ( 1 + v ) should be made.
Considering a material with a microstructure of volume V and external boundary S , the variation of the total potential energy (9) provides, after some algebra [12,43,44], the following relation:
V δ U d V = V [ j ( τ j k i μ i j k ) ] δ u k d V + S R k D δ u k d S + S 1 S 2 p k δ u k d S + E k δ u k ,
where the vectors p k , R k represent the traction and double traction vectors, respectively, defined on the boundary S and written as
p k = n j ( τ j k i μ i j k ) D j ( n i μ i j k ) + ( D m n m ) n i n j μ i j k ,
and
R k = n i n j μ i j k ,
where D j and D represent the tangential and normal gradient operators on S , respectively, and have the form
D j = ( δ j m n j n m ) m D = n m m . .
The vector E k 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:
E k δ u k = { n i m j μ i j k c o r n e r c δ u k f o r 2 D n i m j μ i j k e d g e δ u k d f o r 3 D ,
where denotes the difference of at both sides of the corner c or the edge , while m j 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 F k , boundary tractions p ¯ k , double traction R ¯ k and jump traction E ¯ k , we arrive at the following equilibrium equation:
j ( τ j k i μ i j k ) + F k = 0 ,
accompanied by the classical essential and natural boundary conditions where the displacement vector u k and/or the traction vector p k must be defined on the global boundary S S 1 S 2 , i.e.,
u k ( x ) = u ¯ k ( x ) , x S 1 p k ( x ) = p ¯ k ( x ) , x S 2 ,
and the non-classical essential and natural boundary conditions where the normal displacement vector q k = D u k , the double traction vector R k or jump traction E k are prescribed on S S 3 S 4 , i.e.,
q k ( x ) = q ¯ k ( x ) , x S 3 R k ( x ) = R ¯ k ( x ) , x S 4 E k ( x ) = E ¯ k ( x ) , x corner or edge .
In terms of the displacement vector u ( x ) and free of body forces, Equation (17) obtains the form
( 1 g 2 2 ) [ μ 2 u ( x ) + ( λ + μ ) u ( x ) ] = 0 .
Theorem 1.
The solution of Equation (20) can be written as
u u g r a d i e n t = u c l a s s i c a l + u g ,
with u c l a s s i c a l ker ( μ 2 + ( λ + μ ) ) and u g ker ( 1 g 2 2 ) .
Proof. 
We denote as Δ * the classical elastic elliptic differential operator μ Δ + ( λ + μ ) ( ) Given that ( 1 g 2 Δ ) Δ * u g r a d i e n t = Δ * ( 1 g 2 Δ ) u g r a d i e n t = 0 , we infer that
Δ * u g r a d i e n t = w g ker ( 1 g 2 Δ ) ,
and
( 1 g 2 Δ ) u g r a d i e n t = w e ker ( Δ * ) .
Consequently,
( 1 g 2 Δ ) u g r a d i e n t = w e μ Δ u g r a d i e n t = μ g 2 u g r a d i e n t μ g 2 w e .
Equations (22) and (24) imply that
μ g 2 u g r a d i e n t μ g 2 w e + ( λ + μ ) ( u g r a d i e n t ) = w g u g r a d i e n t = λ + μ μ g 2 ( u g r a d i e n t ) + w e + g 2 μ w g .
We consider the Helmholtz decomposition of the field w g : w g = H + × K . Then, Equation (22) leads to
[ μ Δ u g r a d i e n t + ( λ + μ ) ( u g r a d i e n t ) ] = ( H + × K ) Δ [ ( λ + 2 μ ) u g r a d i e n t H ] = 0 ( λ + 2 μ ) u g r a d i e n t = H + B 0 ,
where B 0 , is a harmonic function. Then, Equations (25) and (26) give
u g r a d i e n t = λ + μ μ ( λ + 2 μ ) g 2 ( H + B 0 ) + w e + g 2 μ w g .
Additionally, it holds that Δ H = w g ker ( 1 g 2 Δ ) and therefore Δ [ ( 1 g 2 Δ ) H ] = 0 ( 1 g 2 Δ ) H = B 1 , where B 1 is another harmonic function. One partial solution of this equation is clearly exactly the function B 1 and so the general solution is H = u g + B 1 , with u g ker ( 1 g 2 Δ ) . Then, (27) is rewritten as
u g r a d i e n t = λ + μ μ ( λ + 2 μ ) g 2 ( B 0 + B 1 ) + w e + g 2 μ w g λ + μ μ ( λ + 2 μ ) g 2 u g .
The first two terms of the decomposition (28) form the classical part u e ker Δ * while the remaining terms form the component u g 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 Ω R d , d = 2 , 3 , 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 Ω D and Ω N (with m e a s ( Ω D ) > 0 ), where classical conditions are imposed, and secondly, in the subdomains Ω Q and Ω R , whose common boundary Γ is of dimension ( d 2 ) and represents the corners c or the edges l of the surface Ω . Then, the general set of mixed-type classical conditions can be formulated as
u ( x ) = f ( x ) , x   Ω D   ( a ) P ( x ) = g ( x ) , x   Ω N   ( b ) ,
along with the set of non-classical conditions
u ( x ) n = h ( x ; g ) ,   x   Ω Q   ( a ) R ( x ) = r ( x ; g ) ,   x   Ω R   ( b ) E ( x ) = s ( x ; g ) ,   x Γ   ( c ) ,
where P = P i x ^ i , R = R i x ^ i and E = E i x ^ i .
The solution u ( x ) and the fields P ( x ) , R ( x ) , E ( x ) depend of course on the parameter g , but this is omitted for simplicity and will be notified only when it is necessary.
It is assumed that the given functions f , g , h , r and s 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 H s ( Ω ) and H s ( Ω ) , 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 H 0 ( Ω ) = L 2 ( Ω ) and its subspace H 1 ( Ω ) , whose elements possess distributional derivatives with the same square integrable behavior. Every Hilbert space H m ( Ω ) , with m ℕ, has a usual inner product and an induced norm consisting of the L 2 norms of the derivatives up to degree m . Therefore, the inclusion H m 1 ( Ω ) H m 2 ( Ω ) when m 1 > m 2 is an obvious relation. For real and positive order s , 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 s 0 , the elements of the space H s ( Ω ) belong to L 2 ( Ω ) , but this is not the case when s < 0 and constitute then pure distributions without functional representative.
Although not necessarily classical functions, the elements of H s ( Ω ) have trace on the boundary Ω , which are generalized functions belonging to H s 1 2 ( Ω ) . The distributional surface normal derivative (of order k) of an element in H s ( Ω ) belongs to H s k 1 2 ( Ω ) in the case that the smoothness of the surface allows the induced differentiability. The space H 0 s ( Ω ) is a subspace of H s ( Ω ) , with elements that, along with all their normal derivatives of order less than s , have zero traces on Ω .
Notice at this point that, for the case of the closed surface Ω , the space H s ( Ω ) is the dual space of H s ( Ω ) and so surface terms of the form h , f | Ω naturally arise, where h H s ( Ω ) and f H s ( Ω ) . This term expresses the action of h on f and it is very reminiscent of the virtual work of a force over a displacement. Only when s = 0 , we have the reduction in the dual pairing h , f | Ω to the usual inner product Ω h ( x ) f ( x ) d S x .
After this brief discussion, we notice that in the framework of a boundary value problem whose differential equation is of order 2 m = 4 , the boundary Equations (29) and (30) involve m ( = 2 ) boundary operators (of possibly mixed type) (when Γ = ). These operators are characterized by their own orders, which obey to the rule: 0 m j 2 m 1 . If, for example, we had exactly the boundary conditions (29) and (30) valid on the whole surface Ω (with Γ = ), the involved boundary operators B j , j = 1 , 2 would be B 1 = I (with order m 1 = 0 ) and B 2 = / n (with order m 2 = 1 ). 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]
u H S ( Ω ) C g ( f H s m 1 1 2 ( Ω D ) + g H s m 2 1 2 ( Ω N ) + h H s m 3 1 2 ( Ω Q )        + r H s m 4 1 2 ( Ω R ) + s H ( s m 4 1 2 ) 1 2 ( Γ ) ) ,
and so, assigning to boundary operators their exact orders, we obtain
u H S ( Ω ) C g ( f H s 1 2 ( Ω D ) + g H s 7 2 ( Ω N ) + h H s 3 2 ( Ω Q ) + r H s 5 2 ( Ω R ) + s H s 3 ( Γ ) ) .
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 s 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 s = 2 m = 4 (the order of the differential equation). Therefore, for the solution u 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 L 2 -sense. However, this is accompanied with the hypothesis of a smooth C 1 , 1 boundary. To permit a boundary with corners or edges and to obtain the broader class of admissible solutions, it is preferable to work with m = 2 . Then,
u H 2 ( Ω ) C g ( f H 3 2 ( Ω D ) + g H 3 2 ( Ω N )             + h H 1 2 ( Ω Q ) + r H 1 2 ( Ω R ) + s H 1 ( Γ ) ) .
In this case, the data f , h loose regularity, but still belong to the realm of square integrable functions. However, the data g , r , and s 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 C g appearing in relation (32) cannot present worse behavior than the asymptotic convergence O ( g 2 ) for g 0 . In addition, (32) implies a fortiori the boundedness
u H 1 ( Ω ) C ( f H 3 2 ( Ω D ) + g H 3 2 ( Ω N )             + h H 1 2 ( Ω Q ) + r H 1 2 ( Ω R ) + s H 1 ( Γ ) ) .
It is noticeable here that the constant C is independent of g 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 g ) behavior, incorporating appropriately the classical solution.
Consider the auxiliary second order classical boundary value problem (titled Problem I), satisfied by the solution u c l a s s i c a l ( x ) :
Δ * u c l a s s i c a l ( x ) = 0 , x   Ω Δ * μ 2 + ( λ + μ ) ,
u c l a s s i c a l ( x ) = f ( x ) , x   Ω D ,
t c l a s s i c a l ( x ) = g ( x ) , x   Ω N ,
where t c l a s s i c a l = τ i j x ^ i x ^ j n ^ , 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 P ( x ) offers its place to its classical counterpart t c l a s s i c a l ( x ) . 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
u c l a s s i c a l H 1 ( Ω ) C ( f H 1 2 ( Ω D ) + g H 1 2 ( Ω N ) ) .
It is worthwhile to mention here that the last equation defines exactly the needed regularity of the data f , g 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 f H 3 2 ( Ω D ) H 1 2 ( Ω D ) as well as g H 1 2 ( Ω Ν ) H 3 2 ( Ω Ν ) .
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:
( f , g , h , r , s ) B = H 3 2 ( Ω D ) × H 1 2 ( Ω N ) × H 1 2 ( Ω Q ) × H 1 2 ( Ω R ) × H 1 ( Γ ) .
The following theorem holds:
Theorem 2.
Let the boundary value problem consist of Equations (20), (29) and (30) where the data ( f , g , h , r , s ) B . This problem is a well-posed fourth order elliptic boundary value problem with a unique solution. This solution can be represented as follows
u ( x ; g ) = u c l a s s i c a l ( x ) + w ( x ; g ) , x Ω w ( x ; g ) = u N ( x ; g ) + g 2 u G ( x ; g ) ,  
where u c l a s s i c a l satisfies Problem I,  u N satisfies the classical elastostatic equation and  u G obeys to the modified Helmholtz equation:
μ 2 u N ( x ; g ) + ( λ + μ ) u N ( x ; g ) = 0 , x Ω ,  
( 1 g 2 2 ) u g ( x ; g ) = 0 , x Ω .
Proof. 
As explained above, the problem (20), (29) and (30) disposes a unique solution u ( x , g ) , which is stable with respect to the data as Equations (32) and (33) guarantee. We consider the decomposition
u ( x ; g ) = u c l a s s i c a l ( x ) + w ( x ; g ) .
The function w ( x ; g ) satisfies Equation (20) given that both u c l a s s i c a l ( x ) and the field u ( x ; g ) obey to this equation too. Based on the Papkovich-type representation [26,30], the field w ( x ; g ) can be written as
w ( x ; g ) = B λ + μ 2 ( λ + 2 μ ) { ( 1 g 2 2 ) B +      [ ( 1 g 2 2 ) B ] r + B 0 } ,
where
( 1 g 2 2 ) 2 B ( x ; g ) = 0 ,
( 1 g 2 2 ) 2 B 0 ( x ; g ) = 0 .
Working with Equation (41), we have that necessarily ( 1 g 2 2 ) B ( x ; g ) is a harmonic function B L ( x ; g ) . Then,
B ( x ; g ) = ( 1 g 2 2 ) B ( x ; g ) + g 2 2 B ( x ; g ) = B L ( x ; g ) + g 2 B G ( x ; g ) ,
where B G ( x ; g ) satisfies the modified Helmholtz equation ( 1 g 2 2 ) B G ( x ; g ) = 0 .
A similar treatment applies to Equation (42) and the result is that the fields of Equations (41) and (42) have the decomposed form
B ( x ; g ) = B L ( x ; g ) + g 2 B G ( x ; g ) B 0 ( x ; g ) = B 0 L ( x ; g ) + g 2 B 0 G ( x ; g ) ,
2 B L ( x ; g ) = 0 , ( 1 g 2 2 ) B G ( x ; g ) = 0 2 B 0 L ( x ; g ) = 0 , ( 1 g 2 2 ) B 0 G ( x ; g ) = 0 .
Inserting these representations in Equation (40) leads to the decomposition
w ( x ; g ) = u N ( x ; g ) + g 2 u G ( x ; g ) ,
where
u N = λ + 3 μ 2 ( λ + 2 μ ) B L λ + μ 2 ( λ + 2 μ ) B L r λ + μ 2 ( λ + 2 μ ) B 0 L ,
and
u G = B G λ + μ 2 ( λ + 2 μ ) B 0 G .
It is evident that u G ker ( 1 g 2 2 ) . In addition, it comes out easily that u N belongs to the kernel of the operator 2 . Consequently, we deduce that
0 = ( 1 g 2 2 ) [ μ 2 w + ( λ + μ ) w ]   = μ 2 u N + ( λ + μ ) u N + ( 1 g 2 2 ) [ μ 2 u G + ( λ + μ ) u G ]   = μ 2 u N + ( λ + μ ) u N ,
from where we infer that u N ker [ μ 2 + ( λ + μ ) ] . □
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 g = 0 as a cornerstone constituent. It is very interesting to examine whether taking the limit of the expression (37) as g 0 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 R ( x ) and E ( x ) are selected as follows:
r ( x ; s ) = g 2 r ˜ ( x ) , x   Ω R ,
s ( x ; g ) = g 2 s ˜ ( x ) , x   Γ ,
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 g , to be a genuine square integrable function, in fact an element of H 1 2 ( Ω Ν ) and that h ( , g ) H 1 2 ( Ω Q ) remains bounded as g varies. In practice, the term h ( x , g ) is usually independent of g .
As far as the stated above restriction of g is concerned, the concept stems from the treatment of the classical solution. As we will see, under this assumption, the function n u c l a s s i c a l ( x ) 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 data ( f , g , h , r , s ) are elements of B with the restriction that r , s obey to Equations (48) and (49), g H 1 2 ( Ω Ν ) H 1 2 ( Ω Ν ) and h ( , g ) H 1 2 ( Ω Q ) is a bounded function of g . Two mutually exclusive alternatives arise:
(A)
h ( x ; g ) = h ( x ) = n u c l a s s i c a l ( x ) ,   x Ω Q   a . e .
(B)
h ( x ; g ) n u c l a s s i c a l ( x ) 0 ,   x Ω Q
In case (B), we impose further regularity on the data by demanding f H 5 2 ( Ω D ) and g H 3 2 ( Ω N ) .
Let  u ( x ; g ) be the unique solution of the gradient elasticity boundary value problem constructed in Theorem 2. Then, the following asymptotic analysis holds:
u G = B G λ + μ 2 ( λ + 2 μ ) B 0 G ,  
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:
u c l a s s i c a l H 1 ( Ω ) C ( f H 1 2 ( Ω D ) + g H 1 2 ( Ω N ) ) ,
as stated before. However, due to the additional regularity f H 3 2 ( Ω D ) H 1 2 ( Ω D ) (valid throughout the work) and the extra requirement g H 1 2 ( Ω Ν ) H 1 2 ( Ω Ν ) 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
u c l a s s i c a l H 2 C ( f H 3 2 ( Ω D ) + g H 1 2 ( Ω N ) ) .
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 ι j u c l a s s i c a l | Ω belong to H 1 2 ( Ω ) as well as i j k u c l a s s i c a l | Ω H 3 2 ( Ω ) .
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 w , can be constructed:
0 = Ω ( ( τ ˜ w μ ˜ w ) ) w d x = R w , D w | Ω + P w , w | Ω + Ε w , w | Γ Ω ( τ ˜ w : ε ˜ w + μ ˜ w ε ˜ w ) d x .
Here, we encounter the dual pairings R w , D w | Ω = R w , D w | H 1 2 ( Ω ) × H 1 2 ( Ω ) , P w , w | Ω = P w , w | H 3 2 ( Ω ) × H 3 2 ( Ω ) and E w , w | Γ = E w , w | H 1 ( Γ ) × H 1 ( Γ ) 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 L 2 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
w ( x ; g ) = 0 ,   x Ω D ,
P w ( x ; g ) = g ( x ) P u c l a s s i c a l ( x , g ) = ( ( D n ^ ( x ) ) n ^ ( x ) n ^ ( x ) D n ^ ( x ) ) : μ ˜ u c l a s s i c a l ( x ; g ) + n ^ ( x ) n ^ ( x ) : D μ ˜ u c l a s s i c a l ( x ; g ) + n ^ ( x ) ( D μ ˜ u c l a s s i c a l ( x ; g ) ) + n ^ ( x ) ( D μ ˜ u c l a s s i c a l 213 ( x ; g ) ) ,   x Ω Ν P w ( x ; g ) = g 2 ( D n ^ ( x ) ( D n ^ ( x ) ) n ^ ( x ) n ^ ( x ) ) : τ ˜ u c l a s s i c a l       + g 2 n ^ ( x ) n ^ ( x ) : D τ ˜ u c l a s s i c a l       + g 2 n ^ ( x ) ( D τ ˜ u c l a s s i c a l ) + g 2 n ^ ( x ) ( D τ ˜ u c l a s s i c a l 213 ) ,
where D = x ^ i D i (see Equation (13)).
Handling the non-classical boundary conditions of w ( x ; g ) leads to
D w ( x ; g ) = n u ( x ; g ) n u c l a s s i c a l ( x ) = h ( x ; g ) n u c l a s s i c a l ( x ) ,   x Ω Q ,
and
R w ( x ; g ) = g 2 r ˜ ( x ) R u c l a s s i c a l ( x ; g ) = g 2 ( r ˜ ( x ) n ^ ( x ) n ^ ( x ) : τ ˜ u c l a s s i c a l ) ,   x Ω R .
Then, Equation (53) obtains the form:
Ω ( τ ˜ w : ε ˜ w + μ ˜ w   ε ˜ w ) d x = g 2 k ( ) , w ( ; g ) | Ω Ν + + R w ( x ; g ) , h ( x ; g ) n u c l a s s i c a l ( x ) | Ω Q + g 2 ( r ˜ n ^ ( x ) n ^ ( x ) : τ ˜ u c l a s s i c a l ) , D w | Ω R + g 2 s ˜ ( x ) , w | Γ g 2 s u c l a s s i c a l ( x ) , w | Γ ,
where
k ( x ) = { ( D n ^ ( x ) ( D n ^ ( x ) ) n ^ ( x ) n ^ ( x ) ) : τ ˜ u c l a s s i c a l + n ^ ( x ) n ^ ( x ) : D τ ˜ u c l a s s i c a l + n ^ ( x ) ( D τ ˜ u c l a s s i c a l ) + n ^ ( x ) ( D τ ˜ u c l a s s i c a l 213 ) } .
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 k ( x ) are well defined and belong to H 3 2 ( Ω Ν ) . 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 u c l a s s i c a l H 2 ( Ω ) , which is bounded by the data as Equation (52) implies, the surface field g 2 k ( x ) is uniformly bounded in H 3 2 ( Ω Ν ) with magnitude of order g 2 .
The fields w and D w , on surfaces and (or) curves have appropriate norms, which are bounded by the norm of the solution w H 1 ( Ω ) (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 g , via an estimate of the form (33) (applied to w ).
If the case (A) holds, the second dual pairing in r.h.s. of Equation (56) disappears and no need to handle the functional R w ( x ; g ) on Ω D 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 n u c l a s s i c a l ( x ) belongs to H 3 2 ( Ω ) (since under this choice, the norm u c l a s s i c a l H 3 ( Ω ) is kept bounded).
Finally, the field R w ( x ; g ) on Ω Q equals with the field g 2 n ^ ( x ) n ^ ( x ) : τ w , which incorporates second order derivatives of w . Then, the adequate surface H 3 2 —norm of R w ( x ; g ) is again uniformly bounded (due to the continuity of the traces with respect to the solution of the boundary value problem) by g 2 w H 1 ( Ω ) , which via (33) is bounded by the norms of the data multiplied by g 2 C . This coefficient goes to 0 as g 0 , since, as mentioned before, the generic constant does not depend on g .
In any case, taking the limit as g 0 , in expression (56), all the terms of the right-hand side converge uniformly to 0, therefore
( τ ˜ w g = 0 : ε ^ w g = 0 + μ ˜ w g = 0   ε ^ w g = 0 ) d x = 0 .
The positivity of the bilinear elastic form guarantees that ε ˜ w g = 0 = 0 ˜ , which implies that w ( x , 0 ) is a constant. Given that the trace w ( x , 0 ) | Ω D = 0 , we infer that
lim g 0 w ( x , g ) = 0 ,   x Ω ,
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  ( f , g , h , r , s ) are analytic functions of their arguments with the restriction that r , s obey to Equations (48) and (49). Let u ( x ; g ) be 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:
u ( x ; g ) u c l a s s i c a l ( x )   as   g 0 ,   x Ω .
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 R w ( x ; g ) on Ω Q is concerned, we could easily find that the surface H 1 2 —norm of R w ( x ; g ) is again uniformly bounded by g 2 w H 2 ( Ω ) . On the basis of (24a), this term is bounded by the norms of the data multiplied by g 2 C g . Therefore, if no extra assumption on data was made, the limit lim g 0 g 2 C g 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 u x x g 2 u x x x x = 0 ,   x ( 0 , 1 ) accompanied by the classical boundary conditions u ( 0 ) = 1 , u x ( 1 ) g 2 u x x x ( 1 ) = 2 and the non-classical ones g 2 u x x x ( 0 ) = 3 g 2 and u x ( 1 ) = 0 .
The solution is
u ( x ) = 2 x + 1 e g 2 + [ ( 3 g 2 + 2 g e 1 g ) e x g + g ( 3 g e 1 g 2 ) e x 1 g ] ( 1 + e 2 g ) 1 ,
while the classical elastic solution that satisfies the conditions u ( 0 ) = 1 and u x ( 1 ) = 2 is u c l ( x ) = 2 x + 1 .
The same differential equation with the boundary conditions u ( 0 ) = 0 , u ( 1 ) = 1 , u x ( 0 ) = 2 and u x x x ( 1 ) = 0 has the following solution:
u ( x ) = x + g [ ( 1 e 2 g ) ( x 1 ) + e x g e x 2 g ] ( g 1 ( g + 1 ) e 2 g ) 1
with corresponding classical solution, the function u c l . ( x ) = x .
Remark 3.
The second example of Remark 2 implies that sometimes it is useful to write Equation (37) in the following form
u ( x , g ) = ( 1 + δ ( g ) ) u c l a s s i c a l ( x ) + u N ( x , g ) + g 2 u G ( x , g )
with δ ( g ) 0 as g 0 , where u N ( x , ξ ) obeys 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 fields u c l a s s i c a l ( x ) and u N ( x ) belong 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 u ( 0 ) = 0 , u ( 1 ) = 1 , u x ( 0 ) g 2 u x x x ( 0 ) = 2 and u x x ( 1 ) = 0 . This problem has the unique solution u ( x ) = 2 x 1 + [ e x g e 1 g e x 1 g ] ( 1 e 2 g ) 1 . However, the classical solution 2 x 1 satisfies the boundary conditions  u ( 1 ) = 1 , u x ( 0 ) = 2 , but does not satisfy the classical boundary condition u ( 0 ) = 0 . Additionally, the part of the solution containing the displacement u G does not converge to zero as g 0 at the endpoint x = 0 . 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 m e a s ( Ω D ) > 0 ). 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]
E I u E I g 2 u + q ( x ) = 0
where x coincides with the neutral axis of the beam, E stands for the material Young modulus, I is the moment of inertia of the beam’s cross-section A , u ( x ) is the transverse beam deflection, q ( x ) is the transverse external load, and g 2 is the intrinsic strain gradient elastic parameter, or the form [53,54,55,56,57,58,59,60,61]
( E I + g 2 A ) u E I g 2 u + q ( x ) = 0 .
The essential difference between these two equations is that the first is derived by considering the strain e x x and the gradient of strain e x x x in the expression of the potential energy density of the beam, while the second one considers the strain e x x and the strain gradients e x x x , e y x x with e y x x being the differentiation of e x x 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 g 2 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 p k ( x 1 , ± a ) = 0 and p 2 ( 0 , x 2 ) = p 2 ( L , x 2 ) = 0 , p 1 ( 0 , x 2 ) = p 1 ( L , x 2 ) = P 1 x 2 / a and the corresponding solution is provided by Selvadurai [3] in the following form:
u 1 c l a s s i c a l = P 1 ( 1 ν 2 ) E a x 1 x 2 u 2 c l a s s i c a l = P 1 2 E a [ ν ( 1 + ν ) x 2 2 + ( 1 ν 2 ) x 1 2 ] .
Consider a rectangular 3D plate of length L and cross-section 2 a × b . The torque at both boundaries of the plate is defined as
M = a a x 2 ( P 1 x 2 a ) d S = a a x 2 ( P 1 x 2 a ) b d x 2 = P 1 I a
with I = 2 b a 3 / 3 being the moment of inertia of the cross-section 2 a × b . Equation (62), valid for the midplane of the plate, obtain the form
u 1 c l a s s i c a l = M ( 1 ν 2 ) E I x 1 x 2 u 2 c l a s s i c a l = M 2 E I [ ν ( 1 + ν ) x 2 2 + ( 1 ν 2 ) x 1 2 ]
while the deflection of the neutral axis is given by u 2 c l a s s i c a l for x 2 = 0 , i.e.,
u 2 c l a s s i c a l = M 2 E I ( 1 ν 2 ) x 1 2
which, for ν = 0 , is identical to the deflection of a classical elastic Bernoulli–Euler beam subjected to pure bending [3,69], i.e.,
u B e r n o u l l i E u l e r c l a s s i c a l = M 2 E I x 1 2 .
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.,
p k ( x 1 , ± a ) = 0 p 2 ( 0 , x 2 ) = p 2 ( L , x 2 ) = 0 p 1 ( 0 , x 2 ) = p 1 ( L , x 2 ) = M I x 2
while in non-classical boundary conditions, we consider the zeroing of double tractions at all the external boundaries, for example:
R k ( x 1 , ± a ) = 0 R k ( 0 , x 2 ) = R k ( 0 , x 2 ) = 0 .
According to the theorem presented in Section 3, the strain gradient elastic solution of this problem is written as
u 1 ( x 1 , x 2 ; g ) = u 1 c l a s s i c a l ( x 1 , x 2 ) + u 1 N ( x 1 , x 2 ; g ) + g 2 u 1 G ( x 1 , x 2 ; g ) u 2 ( x 1 , x 2 ; g ) = u 2 c l a s s i c a l ( x 1 , x 2 ) + u 2 N ( x 1 , x 2 ; g ) + g 2 u 2 G ( x 1 , x 2 ; g ) .
Evidently, the solution (69) imposes the following forms for tractions and double tractions, respectively
p 1 ( n ^ , x 1 , x 2 ; g ) = p 1 c l a s s i c a l ( n ^ , x 1 , x 2 ; g ) + p 1 N ( n ^ , x 1 , x 2 ; g )          + g 2 p 1 G ( n ^ , x 1 , x 2 ; g ) p 2 ( n ^ , x 1 , x 2 ; g ) = p 2 c l a s s i c a l ( n ^ , x 1 , x 2 ; g ) + p 2 N ( n ^ , x 1 , x 2 ; g )          + g 2 p 2 G ( n ^ , x 1 , x 2 ; g )
R 1 ( n ^ , x 1 , x 2 ; g ) = R 1 c l a s s i c a l ( n ^ , x 1 , x 2 ; g ) + R 1 N ( n ^ , x 1 , x 2 ; g )          + g 2 R 1 G ( n ^ , x 1 , x 2 ; g ) R 2 ( n ^ , x 1 , x 2 ; g ) = R 2 c l a s s i c a l ( n ^ , x 1 , x 2 ; g ) + R 2 N ( n ^ , x 1 , x 2 ; g )          + g 2 R 2 G ( n ^ , x 1 , x 2 ; g )
with n ^ 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]
τ 11 c l a s s i c a l = ( λ + 2 μ ) 1 u 1 c l a s s i c a l + λ 2 u 2 c l a s s i c a l = M I x 2 τ 22 c l a s s i c a l = λ 1 u 1 c l a s s i c a l + ( λ + 2 μ ) 2 u 2 c l a s s i c a l = 0 τ 12 c l a s s i c a l = τ 21 c l a s s i c a l = μ ( 2 u 1 c l a s s i c a l + 1 u 2 c l a s s i c a l ) = 0
and
μ 111 c l a s s i c a l = ( λ + 2 μ ) g 2 1 2 u 1 c l a s s i c a l + λ g 2 1 2 u 2 c l a s s i c a l = 0 μ 222 c l a s s i c a l = λ g 2 1 2 u 1 c l a s s i c a l + ( λ + 2 μ ) g 2 2 2 u 2 c l a s s i c a l = 0 μ 112 c l a s s i c a l = μ 121 c l a s s i c a l = μ g 2 ( 1 2 u 1 c l a s s i c a l + 1 2 u 2 c l a s s i c a l ) = 0 μ 122 c l a s s i c a l = λ g 2 1 2 u 1 c l a s s i c a l + ( λ + 2 μ ) g 2 1 2 u 2 c l a s s i c a l = 0 μ 211 c l a s s i c a l = ( λ + 2 μ ) g 2 1 2 u 1 c l a s s i c a l + λ g 2 2 2 u 2 c l a s s i c a l = g 2 M I μ 212 c l a s s i c a l = μ 221 c l a s s i c a l = μ g 2 ( 2 2 u 1 c l a s s i c a l + 1 2 u 2 c l a s s i c a l ) = 0 .
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 n ^ ( n 1 , n 2 ) , have the form, respectively,
p 1 = n 1 τ 11 + n 2 τ 21 + n 1 ( n 1 2 2 ) 1 μ 111 + n 2 ( n 2 2 2 ) 2 μ 221 + n 2 ( n 1 2 1 ) ( 1 μ 121 + 1 μ 211 ) + n 1 ( n 2 2 1 ) ( 2 μ 121 + 2 μ 211 ) + n 2 n 1 2 2 μ 111 + n 1 n 2 2 1 μ 221 p 2 = n 2 τ 22 + n 1 τ 12 + n 1 ( n 1 2 2 ) 1 μ 112 + n 2 ( n 2 2 2 ) 2 μ 222 + n 2 ( n 1 2 1 ) ( 1 μ 212 + 1 μ 122 ) + n 1 ( n 2 2 1 ) ( 2 μ 122 + 2 μ 212 ) + n 2 n 1 2 2 μ 112 + n 1 n 2 2 1 μ 222
and
R 1 = n 1 2 μ 111 + n 1 n 2 μ 121 + n 1 n 2 μ 211 + n 2 2 μ 221 R 2 = n 2 2 μ 222 + n 1 2 μ 112 + n 1 n 2 μ 122 + n 1 n 2 μ 212
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 n ^ ( n 1 , n 2 ) exhibit the following forms, respectively
R 1 = n 1 2 μ 111 + n 1 n 2 μ 121 + n 1 n 2 μ 211 + n 2 2 μ 221 R 2 = n 2 2 μ 222 + n 1 2 μ 112 + n 1 n 2 μ 122 + n 1 n 2 μ 212
and
R 1 c l a s s i c a l ( n ^ , x 1 , x 2 ; g ) = n 1 n 2 μ 211 c l a s s i c a l = n 1 n 2 M I R 2 c l a s s i c a l ( n ^ , x 1 , x 2 ; g ) = 0 . .
Inserting Equations (76) and (77) into (70) and (71), respectively, and satisfying the boundary conditions (67) and (68) one obtains:
p k N ( n 2 , x 1 , ± a ; g ) + g 2 p k G ( n 2 , x 1 , ± a ; g ) = 0 , k = 1 , 2 p k N ( n 1 , 0 , x 2 ; g ) + g 2 p k G ( n 1 , 0 , x 2 ; g ) = 0 , k = 1 , 2 p k N ( n 1 , L , x 2 ; g ) + g 2 p k G ( n 1 , L , x 2 ; g ) = 0 , k = 1 , 2
and
R k N ( n 2 , x 1 , ± a ; g ) + g 2 R k G ( n 2 , x 1 , ± a ; g ) = 0 , k = 1 , 2 R k N ( n 1 , 0 , x 2 ; g ) + g 2 R k G ( n 1 , 0 , x 2 ; g ) = 0 , k = 1 , 2 R k N ( n 1 , L , x 2 ; g ) + g 2 R k G ( n 1 , L , x 2 ; g ) = 0 , k = 1 , 2 .
Equations (78) and (79) indicate that the parts u k N ( x 1 , x 2 ; g ) , u k G ( x 1 , x 2 ; g ) satisfy a homogeneous system of algebraic equations for arbitrary material properties. Apparently, this leads to the conclusion that
u k N ( x 1 , x 2 ; g ) = u k G ( x 1 , x 2 ; g ) = 0
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 ν = 0 , we obtain
u B e r n o u l l i E u l e r c l a s s i c a l u B e r n o u l l i E u l e r g r a d i e n t = M 2 E I x 1 2 .
This result is possible only when Equation (60) is valid. Equation (61) is misleading since the double stress μ y x x does not contribute to the solution of the problem.

5. Conclusions

A material with microstructural effects obeys the simplified elastostatic version of Mindlin’s Form II first-strain gradient elastic theory, and its displacement field u ( x ) satisfies the fourth-order partial differential Equation (3) and the relevant classical and non-classical boundary conditions. In the present work, it has been rigorously proved that the solution of Equation (3) admits the following representation:
u ( x , g ) = ( 1 + δ ( g ) ) u c l a s s i c a l ( x ) + u N ( x , g ) + g 2 u G ( x , g ) ,
which has the following convenient advantages:
  • Incorporates the solution u c l a s s i c a l ( x ) of the respective classical elastic boundary value problem, that satisfies the same classical boundary conditions with the strain gradient elastic problem.
  • Converges to the classical elastic solution as g 0 .
  • Comprises two displacements fields u N ( x , g ) , u G ( x , g ) , which satisfy the simpler equations μ 2 u N + ( λ + μ ) u N = 0 and ( 1 g 2 2 ) u G = 0 , respectively.
The representation of the solution presented above was employed to prove that a strain gradient elastic Bernoulli–Euler beam subjected to pure bending does not present microstructural effects and its behavior is identical to that of a classical elastic Bernoulli–Euler beam. This result is in full agreement with the corresponding conclusions provided by Lurie and Solyaev [69] on the same subject.

Author Contributions

Conceptualization, D.P.; Data curation, A.C.; Formal analysis, A.C. and D.P.; Funding acquisition, D.P.; Investigation, A.C. and T.G.; Writing—original draft, A.C., T.G. and D.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Hellenic Foundation for Research and Innovation (HFRI) under the “First Call for HFRI. Research Projects to support Faculty members and Researchers and the procurement of high-cost research equipment grant” (Project Number: 2060).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this Appendix, the proof of decomposition (5) through the solution representation (4) is provided.
It is well known that the solution of Equation (1), via Papkovich–Neuber potentials, has the form [3,64]
u c l a s s i c a l = B c l a s s i c a l 1 4 ( 1 v ) ( r B c l a s s i c a l + B 0 c l a s s i c a l ) ,
with
2 B c l a s s i c a l = 0 2 B 0 c l a s s i c a l = 0 .
On the other hand, the solution representation (4) reads
u g r a d i e n t = B g r a d i e n t 1 4 ( 1 v ) [ r ( 1 g 2 2 ) B g r a d i e n t + B 0 g r a d i e n t ] ,
with
( 1 g 2 2 ) 2 B g r a d i e n t = 0 ( 1 g 2 2 ) 2 B 0 g r a d i e n t = 0 .
Comparing (A2) with (A4), it is apparent that ( 1 g 2 2 ) B g r a d i e n t B c l a s s i c a l , and subsequently (A3) can be written as
u g r a d i e n t = B g r a d i e n t 1 4 ( 1 v ) [ r B c l a s s i c a l + B 0 g r a d i e n t ] .
By adding and subtracting B c l a s s i c a l and B 0 c l a s s i c a l in (A5), we obtain
u g r a d i e n t = B g r a d i e n t + B c l a s s i c a l B c l a s s i c a l 1 4 ( 1 v ) [ r B c l a s s i c a l + B 0 g r a d i e n t + B 0 c l a s s i c a l B 0 c l a s s i c a l ] u g r a d i e n t = B c l a s s i c a l 1 4 ( 1 v ) [ r B c l a s s i c a l + B 0 c l a s s i c a l ] + B g r a d i e n t B c l a s s i c a l 1 4 ( 1 v ) [ B 0 g r a d i e n t B 0 c l a s s i c a l ] ,
and because of (A1), it is apparent that
u g r a d i e n t = u c l a s s i c a l + B g r a d i e n t B c l a s s i c a l 1 4 ( 1 v ) [ B 0 g r a d i e n t B 0 c l a s s i c a l ] .
However,
( 1 g 2 2 ) ( B g r a d i e n t B c l a s s i c a l ) ( 1 g 2 2 ) ( B 0 g r a d i e n t B 0 c l a s s i c a l ) } = ( A . 4 )   B c l a s s i c a l B c l a s s i c a l B 0 c l a s s i c a l B 0 c l a s s i c a l } = 0 .
Equations (A6) and (A7) easily imply that
u g r a d i e n t = u c l a s s i c a l + u g u c l a s s i c a l ker ( μ 2 ( λ + μ ) ) u g ker ( 1 g 2 2 )
which confirms (5).

References

  1. Timoshenko, S.P.; Goodier, J.N. Theory of Elasticity, 3rd ed.; McGraw-Hill: New York, NY, USA, 1970. [Google Scholar]
  2. Gurtin, M.E. An Introduction to Continuum Mechanics, 1st ed.; Academic Press: Cambridge, UK, 1981. [Google Scholar]
  3. Selvadurai, A.P.S. Partial Differential Equations in Mechanics 2, 1st ed.; Springer: Heidelberg, Germany, 2000. [Google Scholar]
  4. Exadaktylos, G.E.; Vardoulakis, I. Microstructure in linear elasticity and scale effects: A reconsideration of basic rock mechanics and rock fracture mechanics. Tectonophysics 2001, 335, 81–109. [Google Scholar] [CrossRef]
  5. Polyzos, D.; Fotiadis, D.I. Derivation of Mindlin’s gradient elastic theory via simple lattice and continuum models. Int. J. Solids Struct. 2012, 49, 470–480. [Google Scholar] [CrossRef] [Green Version]
  6. Cosserat, E.; Cosserat, F. Theorie des Corps Deformables; A. Hermann et Fils: Paris, France, 1909; Volume 81. [Google Scholar]
  7. Toupin, R.A. Εlastic materials with couple stresses. Arch. Ration. Mech. Anal. 1962, 11, 385–414. [Google Scholar] [CrossRef] [Green Version]
  8. Mindlin, R.D.; Tiersten, H.F. Effects of couple-stresses in linear elasticity. Arch. Ration. Mech. Anal. 1962, 11, 415–448. [Google Scholar] [CrossRef]
  9. Green, A.E.; Rivlin, R.S. Multipolar continuum mechanics. Arch. Ration Mech. Anal. 1964, 17, 113–147. [Google Scholar] [CrossRef]
  10. Koiter, W.T. Couple-stresses in the theory of elasticity, I&II. Philos. Trans. R. Soc. Lond. B 1964, 67, 17–44. [Google Scholar]
  11. Mindlin, R.D. Influence of couple stresses on stress concentrations. Exp. Mech. 1963, 3, 1–7. [Google Scholar] [CrossRef]
  12. Mindlin, R.D. Micro-structure in linear elasticity. Arch. Rat. Mech. Anal. 1964, 16, 51–78. [Google Scholar] [CrossRef]
  13. Germain, P. The method of virtual power in the mechanics of continuous media, I: Second-gradient theory. Math. Mech. Complex Syst. 2020, 8, 153–190. [Google Scholar] [CrossRef]
  14. Vardoulakis, I.; Exadaktylos, G.; Aifantis, E. Gradient elasticity with surface energy: Mode-III crack problem. Int. J. Solids Struct. 1996, 33, 4531–4559. [Google Scholar] [CrossRef]
  15. Exadaktylos, G. Gradient elasticity with surface energy: Mode-I crack problem. Int. J. Solids Struct. 1998, 35, 421–456. [Google Scholar] [CrossRef]
  16. Aravas, N.; Giannakopoulos, A.E. Plane asymptotic crack-tip solutions in gradient elasticity. Int. J. Solids Struct. 2009, 46, 4478–4503. [Google Scholar] [CrossRef] [Green Version]
  17. Gourgiotis, P.A.; Georgiadis, H.G. Plane-strain crack problems in microstructured solids governed by dipolar gradient elasticity. J. Mech. Phys. 2009, 57, 1898–1920. [Google Scholar] [CrossRef] [Green Version]
  18. Karlis, G.; Tsinopoulos, S.V.; Polyzos, D.; Beskos, D.E. Boundary Element Analysis of Mode I and Mixed Mode (I and II) Crack Problem of 2-D Gradient Elasticity. Comput. Methods Appl. Mech. Eng. 2007, 196, 5092–5103. [Google Scholar] [CrossRef]
  19. Karlis, G.; Tsinopoulos, S.V.; Polyzos, D.; Beskos, D.E. 2D and 3D Boundary Element Analysis of Mode-I Cracks in Gradient Elasticity. CMES 2008, 26, 189–208. [Google Scholar]
  20. Lazar, M.; Maugin, G.A.; Aifantis, E.C. Dislocations in second strain gradient elasticity. Int. J. Solids Struct. 2006, 43, 1787–1817. [Google Scholar] [CrossRef] [Green Version]
  21. Gao, X.L.; Park, S. Variational formulation of a simplified strain gradient elasticity theory and its application to a pressurized thick-walled cylinder problem. Int. J. Solids Struct. 2007, 44, 7486–7499. [Google Scholar] [CrossRef]
  22. Georgiadis, H.G.; Anagnostou, D.S. Problems of Flamant–Boussinesq and Kelvin type in dipolar gradient elasticity. J. Elast. 2008, 90, 71–98. [Google Scholar] [CrossRef]
  23. Gao, X.L.; Ma, H.M. Solution of Eshelby’s inclusion problem with a bounded domain and Eshelby’s tensor for a spherical inclusion in a finite spherical matrix based on a simplified strain gradient elasticity theory. J. Mech. Phys. Solids 2010, 58, 779–797. [Google Scholar] [CrossRef]
  24. Karlis, G.F.; Charalambopoulos, A.; Polyzos, D. An advanced boundary element method for solving 2D and 3D static problems in Mindlin’s strain-gradient theory of elasticity. Int. J. Numer. Methods Eng. 2010, 83, 1407–1427. [Google Scholar] [CrossRef]
  25. Papargyri-Beskou, S.; Tsinopoulos, S.V. Lame’s strain potential method for plain gradient elasticity problems. Arch. Appl. Mech. 2015, 9–10, 1399–1419. [Google Scholar] [CrossRef]
  26. Charalambopoulos, A.; Polyzos, D. Analytical solutions for a plane strain gradient elastic rectangle in tension. Arch. Appl. Mech. 2015, 85, 1421–1438. [Google Scholar] [CrossRef]
  27. Anagnostou, D.S.; Gourgiotis, P.A.; Georgiadis, H.G. The Cerruti problem in dipolar gradient elasticity. Math. Mech. Solids 2015, 20, 1088–1106. [Google Scholar] [CrossRef]
  28. Khakalo, S.; Niiranen, J. Gradient-elastic stress analysis near cylindrical holes in a plane under bi-axial tension fields. Int. J. Solids Struct. 2017, 110–111, 351–366. [Google Scholar] [CrossRef]
  29. Gortsas, T.; Tsinopoulos, S.V.; Rodopoulos, D.; Polyzos, D. Strain gradient elasticity and size effects in the bending of fiber composite plates. Int. J. Solids Struct. 2018, 143, 103–112. [Google Scholar] [CrossRef]
  30. Charalambopoulos, A.; Tsinopoulos, S.V.; Polyzos, D. Plane strain gradient elastic rectangle in bending. Arch. Appl. Mech. 2020, 90, 967–986. [Google Scholar] [CrossRef]
  31. Polizzotto, C. A hierarchy of simplified constitutive models within isotropic strain gradient elasticity. Eur. J. Mech. A Solids 2017, 61, 92–109. [Google Scholar] [CrossRef]
  32. Zhou, S.; Li, A.; Wang, B. A reformulation of constitutive relations in the strain gradient elasticity theory for isotropic materials. Int. J. Solids Struct. 2016, 80, 28–37. [Google Scholar] [CrossRef]
  33. Lazar, M.; Po, G. On Mindlin’s isotropic strain gradient elasticity: Green tensors, regularization, and operator-split. J. Micromechanics Mol. Phys. 2018, 3, 1840008. [Google Scholar] [CrossRef]
  34. Gourgiotis, P.A.; Sifnaiou, M.D.; Georgiadis, H.G. The problem of sharp notch in microstructured solids governed by dipolar gradient elasticity. Int. J. Fract. 2010, 166, 179–201. [Google Scholar] [CrossRef] [Green Version]
  35. Andreou, E.; Dassios, G.; Polyzos, D. Korn’s constant for a spherical shell. Q. Appl. Math. 1988, XLVI, 583–591. [Google Scholar] [CrossRef] [Green Version]
  36. Solyaev, Y.; Lurie, S.; Korolenko, V. Three-phase model of particulate composites in second gradient elasticity. Eur. J. Mech. A Solids 2019, 78, 103853. [Google Scholar] [CrossRef]
  37. Solyaev, Y.; Lurie, S.; Altenbach, H.; dell’Isola, F. On the elastic wedge problem within simplified and incomplete strain gradient elasticity theories. Int. J. Solids Struct. 2022, 239–240, 111433. [Google Scholar] [CrossRef]
  38. Lazar, M. On gradient field theories: Gradient magnetostatics and gradient elasticity. Philos. Mag. 2014, 94, 2840–2874. [Google Scholar] [CrossRef] [Green Version]
  39. Ru, C.Q.; Aifantis, E.C. A simple approach to solve boundary value problems in gradient elasticity. Acta Mech. 1993, 101, 59–68. [Google Scholar] [CrossRef]
  40. Eremeyev, V.A.; Dell’Isola, F.; Boutin, C.; Steigmann, D. Linear pantographic sheets: Existence and uniqueness of weak solutions. J. Elast. 2018, 132, 175–196. [Google Scholar] [CrossRef] [Green Version]
  41. Dell’Isola, F.; Sciarra, G.; Vidoli, S. Generalized Hook’s law for isotropic second gradient materials. Proc. R. Soc. A Math. Phys. Eng. Sci. 2009, 465, 2177–2196. [Google Scholar]
  42. Mindlin, R.D. Second gradient of strain and surface-tension in linear elasticity. Int. J. Solids Struct. 1965, 1, 417–438. [Google Scholar] [CrossRef]
  43. Polyzos, D.; Tsepoura, K.G.; Tsinopoulos, S.V.; Beskos, D.E. A Boundary Element Method for Solving 2D and 3D Static Gradient Elastic Problems, Part 1: Integral formulation. Comput. Methods Appl. Mech. Eng. 2003, 192, 2845–2873. [Google Scholar] [CrossRef]
  44. Polyzos, D. A Boundary Element Formulation for Solving Dipolar Gradient Elastic Problems. Comput. Mech. 2005, 35, 292–304. [Google Scholar] [CrossRef]
  45. Marin, M.; Othman, M.I.A.; Seadawy, A.R.; Carstea, C. A domain of influence in the Moore-Gibson-Thompson theory of dipolar bodies. J. Taibah Univ. Sci. 2020, 14, 653–660. [Google Scholar] [CrossRef]
  46. Ghita, C.; Pop, N.; Cioban, H. Quasi-static behavior as a limit process of a dynamical one for an anisotropic hardening material. Comput. Mater. Sci. 2012, 52, 217–225. [Google Scholar] [CrossRef]
  47. Wloka, J. Partial Differential Equations; Cambridge University Press: Cambridge, UK, 1987. [Google Scholar]
  48. Papargyri-Beskou, S.; Tsepoura, K.G.; Polyzos, D.; Beskos, D.E. Bending and stability analysis of gradient elastic beams. Int. J. Solids Struct. 2003, 40, 385–400, Erratum in Int. J. Solids Struct. 2005, 42, 4911–4912. [Google Scholar] [CrossRef]
  49. Challamel, N.; Wang, C.M. The small length scale effect for a non-local cantilever beam: A paradox solved. Nanotechnology 2008, 19, 345703. [Google Scholar] [CrossRef] [PubMed]
  50. Giannakopoulos, A.E.; Stamoulis, K. Structural analysis of gradient elastic components. Int. J. Solids Struct. 2007, 44, 3440–3451. [Google Scholar] [CrossRef] [Green Version]
  51. Challamel, N. Variational formulation of gradient or/and nonlocal higher-order shear elasticity beams. Compos. Struct. 2013, 105, 351–368. [Google Scholar] [CrossRef]
  52. Polizzotto, C. Stress gradient versus strain gradient constitutive models within elasticity. Int. J. Solids Struct. 2014, 51, 1809–1818. [Google Scholar] [CrossRef] [Green Version]
  53. 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]
  54. Shokrieh, M.M.; Zibaei, I. Determination of the Appropriate Gradient Elasticity Theory for Bending Analysis of Nano-beams by Considering Boundary Conditions Effect. Lat. Am. J. Solids Struct. 2015, 12, 2208–2230. [Google Scholar] [CrossRef] [Green Version]
  55. Akgoz, B.; Civalek, O. Strain gradient elasticity and modified couple stress models for buckling analysis of axially loaded micro-scaled beams. Int. J. Eng. Sci. 2011, 49, 1268–1280. [Google Scholar] [CrossRef]
  56. Akgoz, B.; Civalek, O. A size-dependent shear deformation beam model based on the strain gradient elasticity theory. Int. J. Eng. Sci. 2013, 70, 1–14. [Google Scholar] [CrossRef]
  57. Khakalo, S.; Niiranen, J. Form II of Mindlin’s second strain gradient theory of elasticity with a simplification: For materials and structures from nano- to macro-scales. Eur. J. Mech. A Solids 2018, 71, 292–319. [Google Scholar] [CrossRef]
  58. Khakalo, S.; Balobanov, V.; Niiranen, J. Modelling size-dependent bending, buckling and vibrations of 2D triangular lattices by strain gradient elasticity models: Applications to sandwich beams and auxetics. Int. J. Eng. Sci. 2018, 127, 33–52. [Google Scholar] [CrossRef] [Green Version]
  59. Zhao, B.; Liu, T.; Chen, J.; Peng, X.; Song, Z. A new Bernoulli–Euler beam model based on modified gradient elasticity. Arch. Appl. Mech. 2019, 89, 277–289. [Google Scholar] [CrossRef]
  60. Niiranen, J.; Balobanov, V.; Kiendl, J.; Hosseini, S.B. Variational formulations, model comparisons and numerical methods for Euler–Bernoulli micro- and nano-beam models. Math. Mech. Solids 2019, 24, 312–335. [Google Scholar] [CrossRef] [Green Version]
  61. Fu, G.; Zhou, S.; Qi, L. On the strain gradient elasticity theory for isotropic materials. Int. J. Eng. Sci. 2020, 154, 103348. [Google Scholar] [CrossRef]
  62. Papargyri-Beskou, S.; Beskos, D.E. Static Analysis of Gradient Elastic Bars, Beams, Plates and Shells. Open Mech. J. 2010, 4, 65–73. [Google Scholar]
  63. Papargyri-Beskou, S.; Giannakopoulos, A.E.; Beskos, D.E. Variational analysis of gradient elastic flexural plates under static loading. Int. J. Solids Struct. 2010, 47, 2755–2766. [Google Scholar] [CrossRef] [Green Version]
  64. Triantafyllou, A.; Giannakopoulos, A.E. Structural analysis using a dipolar elastic Timoshenko beam. Eur. J. Mech. A Solids 2013, 39, 218–228. [Google Scholar] [CrossRef]
  65. Lazopoulos, K.A.; Lazopoulos, A.K. On a strain gradient elastic Timoshenko beam model. ZAMM Z. Angew. Math. Mech. 2011, 91, 875–882. [Google Scholar] [CrossRef]
  66. Khakalo, S.; Niiranen, J. Lattice structures as thermoelastic strain gradient metamaterials: Evidence from full-field simulations and applications to functionally step-wise-graded beams. Compos. Part B 2019, 177, 107224. [Google Scholar] [CrossRef]
  67. Korshunova, N.; Alaimo, G.; Hosseini, S.B.; Carraturo, M.; Reali, A.; Niiranen, J.; Auricchio, F.; Rank, E.; Kollmannsberger, S. Bending behavior of octet-truss lattice structures: Modelling options, numerical characterization and experimental validation. Mater. Des. 2021, 205, 109693. [Google Scholar] [CrossRef]
  68. Lurie, S.A.; Volkov-Bogorodsky, D.B.; Belov, P.; Lykosova, E. Do nanosized rods have abnormal mechanical properties? on some fallacious ideas and direct errors related to the use of the gradient theories for simulation of scale-dependent rods. Nanosci. Technol. Int. J. 2016, 7, 261–295. [Google Scholar] [CrossRef]
  69. Lurie, S.; Solyaev, Y. Revisiting bending theories of elastic gradient beams. Int. J. Eng. Sci. 2018, 126, 1–21. [Google Scholar] [CrossRef]
  70. Lurie, S.; Solyaev, Y. On the formulation of elastic and electroelastic gradient beam theories. Contin. Mech. Thermodyn. 2019, 31, 1601–1613. [Google Scholar] [CrossRef]
Figure 1. Pure bending of an elastic beam.
Figure 1. Pure bending of an elastic beam.
Mathematics 10 01152 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Charalambopoulos, A.; Gortsas, T.; Polyzos, D. On Representing Strain Gradient Elastic Solutions of Boundary Value Problems by Encompassing the Classical Elastic Solution. Mathematics 2022, 10, 1152. https://doi.org/10.3390/math10071152

AMA Style

Charalambopoulos A, Gortsas T, Polyzos D. On Representing Strain Gradient Elastic Solutions of Boundary Value Problems by Encompassing the Classical Elastic Solution. Mathematics. 2022; 10(7):1152. https://doi.org/10.3390/math10071152

Chicago/Turabian Style

Charalambopoulos, Antonios, Theodore Gortsas, and Demosthenes Polyzos. 2022. "On Representing Strain Gradient Elastic Solutions of Boundary Value Problems by Encompassing the Classical Elastic Solution" Mathematics 10, no. 7: 1152. https://doi.org/10.3390/math10071152

APA Style

Charalambopoulos, A., Gortsas, T., & Polyzos, D. (2022). On Representing Strain Gradient Elastic Solutions of Boundary Value Problems by Encompassing the Classical Elastic Solution. Mathematics, 10(7), 1152. https://doi.org/10.3390/math10071152

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