Euler–Bernoulli elastic beam models of Eringen’s differential nonlocal type revisited within a C0-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf{C }^{0}-$$\end{document}continuous displacement framework

A theory of the Erigen’s differential nonlocal beams of (isotropic) elastic material is prospected independent of the original integral formulation. The beam problem is addressed within a C(0)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C^{(0)}-$$\end{document}continuous displacement framework admitting slope discontinuities of the deflected beam axis with the formation of bending hinges at every cross section where a transverse concentrated external force is applied, either a load or a reaction. Concepts sparsely known from the literature are in this paper used within a more general context, in which the beam is envisioned as a macro-beam whose microstructure is able to take on a size dependent initial curvature dictated by the loading and constraint conditions. Indeed, initial curvature seems to be an effective analytical tool to inject size effects into micro- and nano-beams. The proposed theory is applied to a set of benchmark beam problems showing that a softening behaviour is always predicted without the appearance of paradoxical situations. Comparisons with other theories are also presented.


Introduction
Eringen [1] proposed a method to solve integrodifferential elasticity problems by means of a differential equation whose Green function coincides with the kernel of the integral equation. Eringen et al. [1][2][3][4][5][6][7] mainly applied the proposed differential method to problems with unbounded domains (as for instance wave propagation, crack tip singularities, dislocation analysis problems) whereby the asymptotic conditions at infinity, featured by evanescent values of the response functions, make the differential and the integral methods lead to the same solution. In the case of bounded domains, the differential method-based solution may coincide with the integral method-based one but only if the former solution satisfies-in addition to the standard boundary conditions-some extra nonlocality boundary conditions dictated by the integral method-based solution. This latter assessment is the result of recent research work within beam mechanics [8,9], though roots of it were already known within the mathematical world of integral equation theories [10].
Always with reference to bounded domains, let us recall that, on the one hand, the Eringen's differential method applied to a (fully) nonlocal Euler-Bernoulli (EB) beam model leads to a fourth order governing differential equation in the beam deflection w as [11] w 0000 ðxÞ ¼ 1 D pðxÞ À ' 2 p 00 ðxÞ Â Ã ð0\x\LÞ ð1Þ here the primes denote derivatives with respect to the abscissa x, whereas D ¼ EI ¼ bending stiffness, pðxÞ ¼ distributed load, ' ¼ nonlocality parameter. Equation (1), together with the inherent four standard boundary conditions, leads to a unique solution of the beam problem, but this solution is generally different from the solution of the nonlocal integral problem due to the impossibility to accommodate the mentioned extra boundary conditions [8,9] (a sixth order differential equation would be needed to this purpose). This means that the Eringen's differential nonlocal method, though different from the nonlocal integral method, possesses a capacity of its own to predict size effects. These effects, arising from the existence of a length scale parameter ð'Þ, make the response of the nonlocal beam exhibit either an increased stiffness feature (stiffness hardening, or simply stiffening), or instead a decreased stiffness feature (stiffness softening, or simply softening), the more the larger is '. Therefore the mentioned Eringen's method constitutes a method apart in which size effects are carried in by the applied distributed load, but the properties of this method are perhaps not well understood yet. On the other hand, the mentioned nonlocality boundary conditions read as [8] Mð0Þ À 'M 0 ð0Þ ¼ MðLÞ þ 'M 0 ðLÞ ¼ 0 These, being in evident contrast with the boundary equilibrium conditions, cannot be satisfied and therefore the integral method likely does not admit a solution; in other words, the governing integral equation constitutes a Fredholm integral equation of the first kind, known to lead to a ill-posed boundaryvalue problem with a multiple solution, or no solution at all [10]. Peddieson et al. [11] first used the Eringen's differential method to address nonlocal beams simulating sensor and actuator devices within nanotechnologies, using for this purpose (1) in association with the standard boundary conditions. It was found that Eq. (1) leads to a (unique) solution which predicts size effects of different types, either softening or stiffening, with increasing ', apparently without a precise rule. It was also found the existence of paradoxical beam cases, like the cantilever beam under a tip concentrated load, in which the obtained solution coincides with the classical solution, namely no size effects are predicted by the obtained solution. Another paradoxical case is the cantilever beam under uniform load which is usually referred to as a stiffening beam case, but surprisingly no mention is given from the literature about the beam's tendency to deform raising up against the applied load at higher values of the nonlocality parameter. Several other paradoxical beam cases were reported in the literature [9,17,25,39].
The work by Peddieson et al. [11] stimulated a large amount of research which essentially developed along two streams. A first stream includes a huge amount of works in which the differential method is used to solve bending, buckling and vibration problems for nonlocal beams and plates, for which we make reference to the review papers [12][13][14]. The other stream includes works in which the inconsistencies of the nonlocal integral and differential theories are in some way overcome, often by modifying the constitutive model. The two-phase local/nonlocal constitutive model previously prospected by Eringen [4] and subsequently discussed by Polizzotto [15], was used either within the integral formulation with which the governing equation turns out to be a Fredholm integral equation of the second kind [16][17][18][19][20], or within the differential formulation with which the governing equation is a sixth order differential equation and thus all the boundary conditions (including the extra nonlocality ones) were accomodated [21][22][23][24]. A hybrid constitutive model in which the nonlocal integral model is mixed with the strain gradient one-often also called ''nonlocal strain gradient'' model-was advanced by [25], then developed by [26][27][28][29][30], whereby a constitutive integral equation was used, having as driving variables the strain and strain gradient. A third group of works operates with the Eringen's differential constitutive model itself and reports solutions in which the beam deflection has slope discontinuities at points of application of concentrated loads with size effects in all beam cases, including the so-called paradox cases [31][32][33]. In addition, a bending hinge was introduced by [33] at every clamped end(s) of the beam as well as within the beam axis, with which bending hinge lattice based solutions were provided.
The last quoted papers [31][32][33] provide solutions of the differential nonlocal beam problem out of the usual continuity framework of elasticity theory, with the incorporation of bending hinges; all this by remaining within the original Eringen's nonlocal differential constitutive model. The obtained beam deflections turn out to be size dependent in all cases, including the paradoxical cases. Indeed, this result seems to be the best one for its capacity to better understand the Eringen's differential nonlocal beam model and to solve the paradox beam cases. However, a constitutive-like relationship between the bending hinge mechanism and the related concentrated load, valid also for hinges located at the beam constrained end(s), is still lacking, hence the paradoxical beam cases seem to be not completely solved yet.
In the following, we intend to reconsider and rediscuss the above issues of the literature, in the purpose to build a coherent theory for nonlocal beams governed by Eq. (1), but useful to predict size effects without paradoxes, nor other drawbacks.

Objectives of the present work
In the present paper our concern will be the differential Eq. (1) with the appropriate four standard boundary conditions. The solution of Eq. (1) is the response of the nonlocal beam to the load p augmented by an additional fictitious size dependent load p ð0Þ ¼ À' 2 p 00 ðxÞ. However, by the existence of paradox cases we know that the simulation of size effects by the extra load p ð0Þ is not always effective, at least within the classical framework of continuous displacement solutions. Since p' 2 =D is dimensionally an inverse length, we may interpret this latter term to be proportional to a size dependent initial (inelastic) curvature, say v ð0Þ , so that (1) may be rewritten in the form In this way the solution turns out to be the response of the beam to the load p and to the initial curvature v ð0Þ carring in size effects. The passage from (1) to (3) is more than a formal transformation, because the constitutive-like relation between v ð0Þ and p holds true even if the load is a concentrated force P, no matter if active or reactive, in which case v ð0Þ transforms into the relative rotatiom H of a bending hinge, with a consequent wider class of state variables to describe the beam strain/stress states. All this can be reasonably cast within the mechanics framework of a beam whose microstructure undergoes an initial (inelastic) curvature v ð0Þ leading to solution displacements belonging to the class of C 0 Àcontinuous functions. Furthermore, the so obtained C 0 solutions always resort to a softening behavior of the beam model, without paradoxes, nor other behavioral shortcomings, which means that the Eringen's differential nonlocal method correspondingly gains the effectiveness of a consistent size effect analysis method.
The theory proposed here above recalls an analogy principle known from the literature [34]. According to this principle, the response of a nonlocal EB beam subjected to a distributed load p(x) computed through the Eringen's differential nonlocal method-that is, solving (1)-coincides with the response of the same beam considered of local type and subjected-in addition to the load p(x)-to the inelastic curvature v in ðxÞ v ð0Þ ðxÞ. The validity of this principle is here confirmed and extended to a wider context whereby concentrated forces, either active or reactive, are also considered among the external actions, whereas bending hinges are correspondingly considered as a particular manifestation of the inelastic deformation of the beam.
The outline of the subsequent developments within this paper are as follows. In Sect. 2, the beam model with microstructure is presented together with the concept of initial (inelastic) curvature taken on by the microstructure, along with the concept of bending hinge forming up at any cross section where a concentrated load is applied (but no discontinuity of the deflection curve is generated by a concentrated couple). In Sect. 3, the solution method is presented. Section 4 is dedicated to the application of the proposed theory, as well as to comments and comparisons. Conclusions are drawn in Sect. 5.
The notation will be defined in the text at the first appearance of the related symbols.

EB beam model with microstructure
As explained in the preceding section, the Eringen's differential nonlocal EB beam model is based on a Helmholtz differential equation whose Green function coincides with the kernel function of the related integral constitutive equation. This basic Helmholtz equation may be cast as where vðxÞ denotes the beam's curvature and k is a (dimensionless) constant (0\k\1). Accounting for the equilibrium equation, M 00 ðxÞ þ pðxÞ ¼ 0, (4) loses its differential character and takes on an algebraic form as where M ð0Þ ðxÞ :¼ k' 2 pðxÞ and thus the size dependence of (4) is saved through the initial bending moment M ð0Þ carrying in size effects. A further transformation of (5) is obtained by introducing the quantity v ð0Þ ðxÞ : which coincides with the initial (inelastic) curvature mentioned before. Then, (5) can be rewritten in the form Since the magnitude of v ð0Þ increases with k, this constant characterizes the relative deformation of the microstructure with respect to the continuum, thus k is referred to as the compliance coefficient of the microstructure. For k ! 0 it is v ð0Þ ðxÞ 0, that is, no relative deformation of the microstructure is allowed to occur and Eq. (7) takes on the form of classical elasticity correspondingly, i.e. MðxÞ ¼ DvðxÞ.
Equation (7) is the constitutive equation of an EB beam model with microstructure, whereby the bending moment M is proportional to the net (elastic) curvature v ðeÞ :¼ v À v ð0Þ , difference between the total (compatible) curvature v ¼ Àw 00 and the initial (inelastic) curvature v ð0Þ . Equation (7) recalls an analogous relation given by [34] within an analogy principle, but here (7) has a wider application domain involving concepts as concentrated forces, either active or reactive, and concomitant bending hinges with consequent slope jumps in the deflection curve w(x).
Also, the bending hinge herein introduced appeals to the plastic hinge of limit analysis of beam structures in bending, but these two types of hinges are cenceptually different from each other. In fact, a plastic hinge obeys a threshold law and is irreversible, whereas in contrast a bending hinge obeys a one-toone reversible load-curvature law (6) which is comparable to a thermal-like law, but not to a Hooke law whereby deformation is related to internal forces, as in (7).
As a rule, we shall operate taking k somewhere within the interval (0, 1), but a precise criterion to calibrate the value of k is lacking for the moment. Heuristically, we may just conjecture that the microstructure's relative deformation with respect to the macro-beam, and thus the k coefficient as well, be the smaller the higher is the hyperstatic degree of the beam, such that k may be fixed, for instance, as k ' 1=2 for a statically determinate beam, and k ' 1=4 for a fully clamped beam.

Microstructure
The beam's microstructure is conceived as a continuous set of micro-cell elements as in [35,36]. Here the cells are one-dimensional beam-like (local-type) elements of equal length, say 2d, every element being capable to deform with respect to the beam matrix. The generic element is appended to the beam matrix, for instance having the centroid points of the end cross sections simply supported upon the beam matrix (segment a-b in Fig. 1). The term ''macro-beam'' is used here with the same meaning as the term ''continuum'' within 3D, that is, to indicate the beam together with the microstructure; whereas the term ''beam matrix'' indicates what remains of the macrobeam after removing the microstructure.
The load p acts in a two-fold manner, namely: i) Transmitted to the cell elements through the macrobeam, p acts as a thermal-like load inducing an initial (inelastic) curvature v ð0Þ ¼ k' 2 p=D of the cell elements in accord with (6); ii) The load p also inducesthrough the bending moment M-a net (elastic) curvature v ðeÞ of the macro-beam itself in accord with (7).
Assuming pðx þ nÞ ' pðxÞ within every cell, v ð0Þ can there be treated as constant hence the initial deflection of the generic cell-written to within a rigid translation-is given by which describes the parabola a-c-b of Fig. 1.
This curve is the initial pre-strained configuration of the cell which after application of the load p will adhere to the axis of the beam matrix, so losing its initial deformation. The relative (anticlockwise) rotation of the end cross sections of the beam-like cell in the initial deformation state is given by where P ¼ 2dp is the total (downward) load applied upon the cell, whereas U ¼ k' 2 =D is a constant featuring the compliance of the hinge. On letting p and d vary, but keeping P fixed, H remains constant even for d ! 0 and p ! 1, it therefore holds good also in the presence of a concentrated load of intensity P. Notably, (10) applies to all the cell elements of the beam, including those located at the beam ends where there may exist a constraint giving rise to a reacting (transverse) concentrated force. It is worthwhile to observe that in virtue of (10), which expresses H as the jump of the first derivative of the (small) deflection v of the inherent beam-like cell, the hinge relative rotation H must be of the same order of magnitude as twice that of v 0 , OðjHjÞ % 2Oðjv 0 jÞ. In other words, H must satisfy the inequality We shall return to this point next (Sect. 4.5).

Micro-beam
For d ! 0, the continuous beam-like cell system can be interpreted as a micro-beam possessing the same properties of the cell system. Namely, it undergoes an initial inelastic (thermal-like) curvature v ð0Þ ðxÞ ¼ k' 2 pðxÞ=D at points x where a distributed load is applied, but a relative rotation H i ¼ k' 2 P i =D (bending hinge) at points x i where a concentrated load P i is applied. As shown in next subsection, loading couples, either distributed or concentrated, do not cause effects upon the microstructure, their possible presence is thus ignored for the moment. The above is illustrated in Fig. 2, where the macrobeam AB is subjected to concentrated loads P 1 and P 2 at points C 1 , C 2 , respectively, to a uniform loap p distributed within C 1 and C 2 , along with the concentrated (upward) reacting forces V A and V B at the ends A, B, respectively. Whether the end cross sections are clamped or simply supported is irrelevant for the present reasoning. The initial deformation v ð0Þ of the micro-beam is-to within a rigid motion-uniquely described by a (conventional) initial deflection of the beam in the form of a polyline (with a parabolic side under the uniformly distributed load) like the polyline aA 0 C 0 1 c C 0 2 B 0 b of Fig. 2, showing as many hinges as the number of concentrated forces. Note that this polyline is the effect of the initial deformation upon the beam rendered statically determinate. This implies that in general no rigid motion can exist such as to Fig. 1 Geometrical sketch illustrating a beam-like cell element a-b, which is simply supported by the beam matrix. The load p transmits itself through the macro-beam to the micro-structure causing the beam element to undergo an initial pre-strained curvature v ð0Þ ¼ ðk' 2 =DÞp, whereas the macro-beam deforms according to the constitutive law M ¼ Dðv À v ð0Þ Þ make the initial deflection line compatible with the beam's constraints at A and B, except that the beam is statically determinate. Since by global equilibrium V A þ V B À P 1 À P 2 À pL 1 ¼ 0, it results that the relative rotation of the beam ends is zero, that is, the end sides aA 0 and B 0 b of the polyline must be parallel to each other.
V A and V B must be evaluated through equilibrium considerations, as better explained shortly in Sect. 4.

No displacement jump under a concentrated couple
The features of the micro-beam illustrated above do not consider the possible existence of concentrated couples within the loading upon the beam. Here we show that concentrated couples do not produce any discontinuity of the micro-beam initial configuration. For this purpose, let the macro-beam be subjected to two concentrated loads, namely a downward load P at x 0 À e and an upward load P at x 0 þ e, Fig. 3. The initial configuration of the beam is the trilateral line abcd with ab parallel to cd and the bending hinges at the application points of the concentrated loads. On letting e ! 0, P ! 1, but 2eP ¼ C ¼ constant, at the limit the trilateral line abcd tends to the broken line ab 0 c 0 d. The transverse segment b 0 c 0 of length ' 2 C=D Fig. 2 Geometrical sketch showing a beam subjected to two concentrated loads P 1 , P 2 at points C 1 , C 2 , respectively, to a uniform load p distributed within C 1 and C 2 , along with the concentrated reacting forces at the constrained ends A, B where non-zero bending moments are allowed to occur. The initial prestrained deformation includes four bending hinges at points A, C 1 , C 2 , B. Since V A þ V B À P 1 À P 2 À pL 1 ¼ 0, then the relative rotation of the end outer cross sections must be vanishing, therefore the (conventional) initial configuration of the beam is a polyline like aA 0 C 0 1 cC 0 2 B 0 b, with extreme sides aA 0 and B 0 b parallel to each other and a parabolic side under the uniformly distributed load p Fig. 3 Geometrical sketch used to show that no discontinuity is allowed to arise at a point where a concentrated couple is applied, which is a direct consequence of the shear rigidity of the EB beam model represents a displacement jump at x 0 due to the concentrated (anticlockwise) couple C applied at x 0 . Since this latter deformation mechanism is of pure shear nature, and since shear strain is forbidden within the context of the present EB beam model, it follows that the initial configuration of the micro-beam cannot take on the broken form ab 0 c 0 d, but it instead will be some continuous line (as e.g. the dotted line b 1 c 1 of Fig. 3). Therefore, concentrated couples are allowed to play as external actions and thus to have an influence on the beam's response, but they do not produce effects upon the microstructure.

The solution method
For a better understanding of the present method, it may be useful to recall that the boundary ends of a nonlocal beam must be thought of as two boundary layers each of small length, say e, such that the effective length of the macro-beam is L þ 2e and the beam ends are located at x ¼ Àe and x ¼ L þ e. Also, differentiations must be intended in a distribution sense.
The beam model under consideration obeys the constitutive Eq. (7) which, after a double differentiation and recalling that v ¼ Àw 00 , leads to the displacement governing Eq. (3). As an illustrative example we consider a simple beam of length L subjected to a distributed load p(x), to concentrated loads P i ði ¼ 1; :::; mÞ applied at points x i and to a concentrated couple C 0 at x 0 , 0\x 0 \L. The beam has some constraints at the ends giving rise to (upwards) concentrated reaction forces V A , V B . The solution of the beam problem for assigned static loads is derived through Eq. (3) here reported again in the following shape, where the presence of concentrated forces and couples is explicitly indicated through the Dirac delta distribution d, namely, By (7), the bending moment M(x) associated to (12) is MðxÞ ¼ ÀDw 00 ðxÞ À Dv ð0Þ ðxÞ ð13Þ Moreover, as stated before, the initial curvature v ð0Þ is uniquely determined from the external applied forces, either active or reactive. Therefore, we can write v ð0Þ ðxÞ as Equation (14) predicts a bending hinge at every point where a concentrated force, either active or reactive, is applied, therefore the couple C 0 is not there involved. The reactive forces V A and V B are given by and satisfy the global equilibrium equation Next, through the relation v ð0Þ ðxÞ ¼ Àw 00 ð0Þ ðxÞ, Eq. (14) can be integrated to get what can be called conventional initial deflection of the beam, namely, where f(x) is a particular function satisfying the equation f 0000 ðxÞ ¼ pðxÞ 8x 2 ð0; LÞ. Equation (17) gives-to within a rigid motion-the displacement effects produced by the application of the initial curvature upon the beam rendered statically determinate.
Next, substituting (14) into (12) and by integration of (12) we can write: where C 1 ; C 2 ; C 3 ; C 4 are some constants. Substituting (18) into (13) gives and thus, by (15), we have whereby the global equilibrium of the beam is satisfied, namely, At this point it remains to evaluate the four constants C 1 ; C 2 ; C 3 ; C 4 for which the boundary conditions must be used. The resulting procedure will be illustrated in next section dedicated to applications. What has to be pointed out here is that the applied couple C 0 has an influence on the deflection w(x) and the bending moment M(x) of the beam, but not upon the microstructure's deformation (both v ð0Þ and w ð0Þ do not depend on C 0 ). In concluding the present section, we observe that the proposed method can be extended to beams in vibration and buckling, but this extension is left open for the moment. For easy reference, the beam model governed by the above relations will be called shortly as ''C 0 -beam model''.

Applications to benchmark beam cases
A few benchmark beam cases are addressed as illustrative examples. The procedure presented in Sect. 3 is followed hereafter.

Cantilever beam under point load
A cantilever beam is here considered, which is subjected to a concentrated load P at some intermediate point C as sketched in Fig. 4a. There are two concentrated forces, namely, a downward force P applied at point C of abscissa x\L, and an upward reacting force, say V A , at the left end A. Therefore, there arise two bending hinges at points A and C, respectively, such that the related initial inelastic curvature v ð0Þ reads v ð0Þ ðxÞ ¼ By global equilibrium we have V A ¼ P. Equation (22) can be integrated to obtain a conventional initial deflection of the beam w ð0Þ ðxÞ counterpart of (17), that is, which is sketched in Fig. 4b (solid line). As the beam is statically determinate, the initial deflection w ð0Þ ðxÞ is compatible.
The C 0 solution of the beam problem is obtained by integration of the governing differential equation together with (13) and taking into account (22) along with the boundary conditions Written for the interval (0, L), the C 0 deflection proves to be whereas the bending moment is given by Since w ð0Þ ðxÞ is compatible, then we have wðxÞ ¼ w c ðxÞ þ w ð0Þ ðxÞ, whereas the bending moment coincides with its classical form. For x ¼ L (26) and (27) become respectively: MðxÞ ¼ À PðL À xÞ ð0 x LÞ ð 29Þ In Fig. 4c the C 0 solution w(x) of (28) is reported as a function of x/L, normalized with respect to the maximum deflection w c ðLÞ ¼ PL 3 =3D and for different values of k ¼ '=L (= 0., 0.25, 0.5, 0.75, 1.). Analogously, in Fig. 4d the normalized maximum C 0 deflection g :¼ wðLÞ=w c ðLÞ is plotted as a function of k. The value k ¼ 0:5 has been used for the plots of Fig. 4c, d. We remark that the solution (26) differs from the analogous solution by [33], where no bending hinge is considered at point C where the load P is applied. For x ¼ L, the classical paradox beam case is obtained, but without any paradox therein. A softening behavior is exhibited by the beam (gðkÞ [ gð0Þ 8 k 6 ¼ 0).

Cantilever beam under uniform load
A cantilever beam under a uniformly distributed load p is here considered (Fig. 5a), which is usually referred to as the stiffening beam case. There is only one (upward) concentrated force, V A ¼ pL at the clamped end A, hence only one bending hinge at the same point. The initial curvature is expressed as This, by integration, gives the (compatible) conventional initial deflection w ð0Þ ðxÞ, namely, which is reported in Fig. 5b. Proceeding as in the previous case and considering the boundary conditions the C 0 deflection, written for the interval (0, L), reads as which asserts a softening behavior of the beam. Also it is In Fig. 5c the normalized C 0 deflection wðxÞ=w c ðLÞ is plotted as a function of x/L for different values of k ¼ '=L (= 0., 0.25, 0.5, 0.75, 1.), whereas in Fig. 5d the maximum C 0 deflection is plotted as a function gðkÞ. Again the value k ¼ 0:5 has been used. We observe that the C 0 solution (33) coincides with the analogous one by [33], which was derived adding to the classical solution the effects of relaxing the kinematics at the clamped end. Again a softening behavior is exhibited by the beam (gðkÞ [ gð0Þ 8 k 6 ¼ 0).

Clamped-pinned beam under uniform load
Here, a clamped-pinned beam under uniformly distributed load p is considered (Fig. 6a). As there are two concentrated forces, namely the reactions V A and V B at the beam ends, there arise two bending hinges correspondingly, hence the initial curvature reads The related conventional initial deflection proves to be expressed as which is plotted in Fig. 6b (solid line). As the beam under consideration is statically indeterminate, the conventional initial deflection w ð0Þ ðxÞ of (36) is incompatible, since in fact w ð0Þ ðLÞ 6 ¼ 0, though w 0 ð0Þ ð0 À Þ ¼ 0. Next, considering the boundary conditions by integration of the pertaining governing differential equation we get the C 0 deflection as and the bending moment M(x) as both of which hold within the interval (0, L). For k ¼ 0, (38) and (39) provide the classical solution. The hinge relative rotation at x ¼ 0 is found to be expressed as In Fig. 6c  reported as a function gðkÞ. It is shown that the C 0 plot exhibits a softening behavior like that of Peddieson et al. [11], but in a slightly more pronounced proportion; for instance for k ¼ 0:5, g C 0 ð0:5Þ ¼ 2:74 while g Pedd ð0:5Þ ¼ 2:50. The value k ¼ 0:25 has been used.

Clamped-pinned beam under a concentrated load
Here a clamped-pinned beam subjected to a concentrated load P at x ¼ L=2 is considered (Fig. 7a). This time we deal with three concentrated forces, namely the downward force P at midle point C, along with the (upward) reacting forces V A and V B . Therefore there arise three bending hinges at the same points and the inherent initial curvature proves to be to which we can associate the (incompatible) conventional initial deflection as which is reported in Fig. 7b (solid line). Next, considering the boundary conditions and recalling (41), by integration of the governing differential equation the C 0 deflection function w(x) proves to be Both (44) and (45) hold for 0 x L; they recover the respective classic form for k ¼ 0. The hinge relative rotation at A is given by In Fig. 7c the C 0 deflection w(x) is reported as a function of x/L, normalized with respect to the w c ðL=2Þ, for different values of k ¼ '=L (= 0., 0.25, 0.5, 0.75, 1.). Also, in Fig. 7d the normalized maximum C 0 deflection wðL=2Þ=w c ðL=2Þ is plotted as a function gðkÞ. Both plots of Fig. 7c, d show that the C 0 model exhibits a softening behavior. Again k ¼ 0:25 has been assumed.

Remarks on the obtained results
The main result of the applications presented above is that in all beam cases a softening behavior of nonlocal gradient beam models is exhibited, without drawbacks of any sort, and this for all loading and constraint conditions. This result is in strong contrast with the outcome from Peddieson et al. [11], where the displacement solution of analogous beam problems was searched for within the classical continuity framework. The work by Peddieson et al. [11] showed that the solution of differential nonlocal beam problems generally indicates a softening behavior of micro-and nano-beams, but that the continuity restriction therein maintained caused dramatic limitations of the size dependent response in some beam cases. For instance, either size effects are not predicted at all whenever the beam is subjected to concentrated loads (as for a cantilever beam under point load), or the beam subjected to uniform load exhibits stiffening effects increasing with the nonlocal scale parameter, the more the larger is this parameter, which leads the beam to deflect upwards against the (downward) load (as in the case of a cantilever beam under uniform load). As stated before, the examples reported above clearly indicate that the C 0 -continuous solutions of the fourth order differential equation and boundary conditions indicate a softening behavior for all beam cases without paradoxes and no other shortcomings. Indeed, no paradoxical situations are encountered if one accepts the presence of bending hinges within the context of elastic beams, along with the consequent presence of slope discontinuities. However, as mentioned at the end of Sect. 2.1, the presence of bending hinges carries in limitations as (11) to the amplitude of the hinge relative rotations. The relative rotations H A at the left clamped ends of the beam samples worked out before are all of the form where P denotes the total load applied upon the beam, whereas fðkÞ is a case dependent coefficient, that is fðkÞ ¼ This result indicates that 0\fðkÞ 1 8k (and likely for almost all other beam cases). The factor / may be interpreted as the rotation of the end cross section of a cantilever beam of length ffiffi ffi 2 p L and stiffness D ¼ EI, subjected to concentrated load P at the free end. Within the framework of classical linearized elasticity featured by infinitesimal displacements (not too high load P, nor too small stiffness D) we may assume that /\\1. Therefore we may assert that which obviously can be extended to all bending hinges of the beam. Indeed, the amplitude of the hinge relative rotations is infinitesimal within the present theory. On the other hand, a restriction on the boundedness of k seems to be in good agreement with an analogous restriction discending from the requirement that the nonlocal continuum model saves its own effectiveness versus atomic lattice models, to which it in fact tends to resort for k ¼ '=L values larger than unit [1,27,38]. For the above reasons here we have considered k values within the interval (0, 1).

Concluding remarks
We have presented a method to solve EB beam models within the framework of C 0 -continuous displacement functions, in which the beam deflection function is allowed to exhibit slope discontinuities at the application points of concentrated forces, either loads, or reacting forces. Concepts as bending hinge and associated relative rotation, which were sparsely advanced within the literature, are here rationally collected within a consistent continuum mechanics framework in which the beam is endowed with a microstructure undergoing an initial (inelastic) curvature carring in the inherent size dependent effects. The resulting C 0 -beam model leads to a unique solution of the beam problem, featured by a softening behavior for all beam's load and constraint cases, without paradoxical situations or other shortcomings. In contrast with the original Eringen-Peddieson model, the present one may widely be used for size effect analysis of beams in bending, buckling and vibration in a consistent manner.
Funding The authors declare that this study did not receive funding.

Declarations
Conflict of interest The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.