We present a comprehensive study on the postbuckling response of nonlocal structures performed by means of a frame-invariant fractional-order continuum theory to model the long-range (nonlocal) interactions. The use of fractional calculus facilitates an energy-based approach to nonlocal elasticity that plays a fundamental role in the present study. The underlying fractional framework enables mathematically, physically, and thermodynamically consistent integral-type constitutive models that, in contrast to the existing integer-order differential approaches, allow the nonlinear buckling and postbifurcation analyses of nonlocal structures. Furthermore, we present the first application of the Koiter’s asymptotic method to investigate postbifurcation branches of nonlocal structures. Finally, the theoretical framework is applied to study the postbuckling behavior of slender nonlocal plates. Both qualitative and quantitative analyses of the influence that long-range interactions bear on postbuckling response are undertaken. Numerical studies are carried out using a 2D fractional-order finite element method (f-FEM) modified to include a combination of the Newton–Raphson and a path-following arc-length iterative methods to solve the system of nonlinear algebraic equations that govern the equilibrium beyond the critical points. The present framework provides a general foundation to investigate the postbuckling response of potentially any type of nonlocal structure.
The study of buckling and postbuckling responses of elastic structures have obvious applications in the domain of structural analysis and design. Seminal work in the study of postbuckling response of continuum structures was conducted by Koiter . Later, Budiansky  extended this study to the postbuckling response of elastic structures. These works generated broad interest in the academic and research communities involved in the structural design . By employing models and methodologies derived, either directly or indirectly, from these seminal works, the postbuckling response of several structures have been successfully studied [4–8]. These studies are based on the classical (local) constitutive theories that involve a point-wise correspondence between the displacement and the stress field generated at the point. Although the localized relations of constitutive theories are still fundamental for most studies on the behavior of solids, the existence of nonlocal interactions between distant points in a solid have been established through experimental investigations. These observations point toward a nonclassical behavior that is not accounted for in classical local elasticity theories. Some common examples of solids exhibiting nonlocal behavior include (microscale and nanoscale) structures where size effects gain prominence [9,10]. In addition, recent experimental observations point to a more pervasive presence of nonlocal interactions even in macroscale applications such as sandwich structures, stiffened panels, and functionally graded materials [11,12].
Owing to the importance of the applications mentioned earlier, several theories have been proposed to model the nonlocal effects in elastic solids. Among the numerous theories existing in the literature, the strain-based integral constitutive relations for the nonlocal solids proposed by Kröner  and Eringen  are likely the most prominent. In these approaches, nonlocal interactions are introduced in the constitutive law by means of an integral stress-strain relation defined over an assigned domain of influence. This integral effectively captures long-range interactions by means of a convolution of the strain tensor weighted by a proper attenuation function. While this approach has certain key advantages, primarily in providing an intuitive representation of the long-range interactions via convolution integrals, it suffers some key mathematical and physical inconsistencies. The salient deficiency of integral models is that it belongs to an ill-posed class of integral equations that do not admit unique solutions [15–17]. To circumvent the implicit formulations associated with the integral approaches, an equivalent gradient-based model of nonlocal elasticity was developed . The modeling advantages introduced by these gradient-based approaches sparked numerous investigations focused on assessing the effect of long-range interactions on the response of structural elements [19–21]. The equivalence between the gradient and integral models is based on certain restrictions like the nonlocal influence domain being unbounded and on certain choices of the attenuation kernel. These conditions place physical restrictions that do not allow the use of gradient models for a general study on nonlocal elasticity. This is because violation of these restrictions would present paradoxical observations, such as the vanishing nonlocal effects, encountered for certain choices of loading and boundary conditions [15,17]. In the context of stability analysis, the nonlocal effects are realized by a modification of the system stiffness terms. When employing integral-type nonlocal constitutive laws, the material stiffness is reduced due to nonlocal interactions [22,23], whereas differential models cause an increase in the geometric stiffness [18,24]. In either case, both models predict a lower critical load in nonlocal structures when compared with their classical local counterpart . Note that most of the literature addressing the effects of nonlocal mechanisms on the critical load has focused on linear buckling. Studies on the critical load in the case of nonlinear buckling and postbuckling response are almost entirely based on differential models. Among the existing studies, nonlocal effects are restricted to carbon nanotubes (CNTs) and graphene. Both classical and modified differential models of nonlocal elasticity were employed to compare the effects of long-range interactions on both the hardening and softening response of low-dimensional structures [25–27]. Additional studies include postbuckling of nonlocal beam and cylinder models for CNTs [28,29] and plate and shell models for graphene sheets [30,31]. Nonlocal effects have also been considered in the study of the postbuckling response of bio-based protein tubules  and in the pull-in instability of the micro- and nano-electromechanical devices [33,34]. To the best of author’s knowledge, all the reported studies are based on the differential models of nonlocal elasticity and conclusively predict a lower value of critical load across boundary conditions followed by a reduction in the postbuckling resistance to deformation.
The lack of similar studies based on integral models of nonlocal elasticity may be attributed to the relative complexity of the integral-type constitutive relations, and the unavailability of a consistent and positive semi-definite definition of the deformation energy density functional within this framework. In fact, it has been argued that such a functional cannot be defined for the nonlocal solids . Since postbuckling analysis is based on energy approaches, this latter point presents a significant technical gap that requires to be addressed. Moreover, it must be noted that the integral and differential-type constitutive laws discussed earlier are thermodynamically inconsistent . This is because of the ad hoc constitutive laws that do not satisfy the thermodynamic balance in a rigorous manner. More specifically, when nonlocal solids are modeled via integral or differential approaches, the first and second law of thermodynamics are only satisfied in a weak manner (over the nonlocal domain) as opposed to the physically consistent localized (strong) form . Thus, a consistent energy formulation may not be developed employing these models. This consideration extends also to the alternative integral-type two-phase definition of constitutive law that includes both local and nonlocal terms in the stress–strain relations. While the two-phase definition admits unique solutions by virtue of the well posedness of Fredholm’s second integral, it does not satisfy the second law of thermodynamics (Clausius–Duhem inequality) in a rigorous manner. This highlights the limitations of the existing constitutive models for nonlocal elasticity and their deficiencies in an analysis over postbuckling studies .
Fractional calculus has demonstrated enormous potential for addressing areas where integer-order models fall short, such as the case of modeling of complex phenomena. Limiting the discussion to the current context of structural models, recent developments include fractional-order constitutive theories of nonlocal elasticity that can potentially address the aforementioned deficiencies [37–41]. Other prominent examples of successful applications of fractional-order theories are the constitutive modeling of memory effects in viscoelastic materials [42,43] and nonlocal effects in a multiscale analysis [38,41,44]. In the case of fractional-order constitutive theories, the differ-integral nature of fractional-order derivatives serves as a suitable alternative to the integro-differential constitutive relations encountered in nonlocal elasticity. Indeed, the fractional-order derivatives appearing in the kinematic strain–displacement relations allow accounting for long-range interactions within the constitutive model of a nonlocal elastic solid. Successful applications of this theory to the development of mathematical and physically consistent models for softening  and stiffening  type of nonlocal interactions have also been proposed. The frame-invariant fractional-order continuum theory based on a physically consistent definition for position-dependent nonlocal length scales is particularly effective in allowing a rigorous application of the thermodynamic principles . This observation is crucial in developing consistent constitutive laws for nonlocal elasticity because it allows defining a positive semi-definite functional corresponding to the deformation energy density of a nonlocal solid. It has been proved that the positive semi-definite form of the deformation energy density functional leads to self-adjoint linear operators in the governing equations [45,46]. By employing weak-form expressions and variational methodologies, finite element techniques suitable for the numerical solution of boundary value problems based on fractional-order equations have also been developed .
Following energy-based methods afforded thanks to the fractional-order model, softening effects due to nonlocal interactions have been accounted for both the linear and nonlinear response of beams and plates [45–48]. The large deformation analysis of the nonlocal solid requiring geometrically nonlinear strain–displacement relations is particularly insightful within the context of stability of fractional-order continuum theory. Unlike the integral and differential models of nonlocal elasticity that modify either the material or the geometric stiffness term alone, physically the nonlocal interactions are expected to be realized on both of these stiffness terms. This is achieved by incorporating the nonlocal interactions at the level of kinematic relations within the constitutive relations for a nonlocal solid. Thus, the effect of the nonlocal interactions over the critical load is decided by the net result of modifications to the material and geometric stiffness terms by the fractional-order derivatives. This observation has been drawn following rigorous parametric studies isolating the nonlocal interactions in each term to identify the net result over the critical load corresponding to linear buckling . It is imperative that a study of the postbuckling response of nonlocal structures, which also depends on both of these stiffness terms, is conducted employing the mathematically, physically, and thermodynamically consistent nonlocal model of the fractional-order continuum theory.
With this objective in mind, we study the critical load corresponding to the nonlinear buckling and the subsequent postbuckling response of nonlocal structures. For this purpose, we build upon the existing framework for the stability analysis of fractional-order structures . As mentioned earlier, this study is possible due to the positive semi-definite deformation energy density functional made available by the fractional-order formulation of nonlocal elasticity [36,45]. The objective of the current work is twofold. First, the nonlinear buckling of fractional-order structures is studied by identifying the singular points of the equilibrium curves. We develop the nonlinear equations required to obtain the critical load and the associated buckling mode following variational principles. Subsequently, the postbifurcation analysis in the immediate vicinity of the singular points is conducted. Second, the developed framework is employed to perform a numerical study (via fractional-order finite element method (f-FEM)) and to assess the effect of nonlocal interactions on the postbuckling response of a slender plate.
This article is organized as follows. First, it provides a brief summary of the fractional-order continuum theory. This is followed by a summary of the energy-based framework employed for the identification of critical points of nonlocal structures corresponding to linear buckling. Then, it presents the derivation of the critical loads corresponding to singular points of nonlinear equilibrium. Thereafter, by means of extending the Koiter’s asymptotic method, the postbuckling response in the immediate vicinity of the critical points is studied. Finally, an analysis of the postbuckling behavior of fractional-order Kirchhoff plates is performed. We provide brief details regarding the 2D f-FEM model developed for this numerical study over the postbuckling response of nonlocal structures.
2 Postbuckling Response of Fractional-Order Nonlocal Structures
In this section, we develop the models required to study the postcritical response of fractional-order nonlocal structures. To provide the necessary mathematical background and notation, we briefly review the basic constitutive relations for the fractional-order continuum theory [41,47] and the stability of equilibrium points for a nonlocal solid evaluated following this theory . We use the boldface notation to describe tensors (including vectors), while subscript indicates the measure number of the component evaluated in Cartesian coordinates.
2.1 Constitutive Model for Fractional-Order Structures.
2.2 Critical Points of Stability for Fractional-Order Structures.
We point to a simultaneous modification of the standard and geometric stiffness of the nonlocal structure (i.e., the numerator and denominator in the above expression, respectively) caused by the long-range interactions. This is because of the fractional-order derivatives being present in the expressions of both stiffness terms in the aforementioned result for critical load of nonlocal structure. Therefore, the effect of the nonlocal interactions on the critical load is not trivial and may either increase or decrease due to long-range effects. This is in contrast with the observations following from existing integer-order nonlocal theories that predict a consistent decrease in the critical load. We clarify this statement by a comparison of the fractional-order model with the standard integer-order Eringen’s integral and differential theories of nonlocal elasticity. Following the integral theory , the effect of the nonlocal interactions is captured only within the material stiffness terms. This causes a reduction of the material stiffness with an increasing degree of nonlocality. However, the geometric stiffness evaluated using this theory is identical to the corresponding result following from classical (local) elasticity . Conversely, the Eringen’s differential theory increases the geometric stiffness terms as a consequence of the nonlocal interactions, while leaving the material stiffness unchanged from classical (local) elasticity predictions . Please refer to Ref.  for a detailed comparison of the critical buckling load evaluated for a nonlocal solid employing the fractional-order and integer-order nonlocal continuum theories. While these observations are drawn for the linear buckling of fractional-order structures, they are relevant for nonlinear buckling studied using the nonlinear eigenvalue problem given in Eq. (9).
2.3 Postbifurcation Analysis of Fractional-Order Structures.
Following the previous study on the identification of singular points of the fundamental solution of the governing equations, our next objective is to analyze the nature of the singularity and develop the models required for the subsequent analysis of the nonlinear response of fractional-order solids beyond these critical points. In the domain of classical elasticity, this study of the postbuckling response is carried out following the energy methodology . Here, we make use of the positive semi-definite deformation energy density (see Eq. (5)) achievable via the fractional-order derivatives for nonlocal solids. We proceed with a straightforward extension of the Koiter’s asymptotic method, originally defined for integer-order models of classical elasticity, to obtain the postbuckling response of fractional-order nonlocal structures.
Increasing the degree of nonlocality would result in lower values for fractional derivatives. Thus, the numerical values of and in Eq. (10) will be reduced for the fractional-order structure. This would result in lower value of for fractional-order solids with the increasing degree of nonlocality. Therefore, we infer from Eq. (13) that nonlocal effects will reduce the difference in energy barrier between the fundamental and the bifurcated paths at the singular points.
3 Postbuckling Analysis of Fractional-Order Kirchhoff Plates
In this section, we apply the aforementioned formulation to study the postbuckling response of a nonlocal Kirchhoff plate. We begin with the derivation of the governing equations of equilibrium for a fractional-order Kirchhoff plate following variational principles. We also develop the numerical model for these nonlinear fractional-order (integro-differential) governing equations employing the path-following method based on arc-length constraints.
3.1 Basic Formulation for Fractional-Order Kirchhoff Plates.
Following the Kirchhoff plate theory, we note that the in-plane stiffness of the plate is much larger than the bending stiffness . This results in transverse displacements, perpendicular to the mid-plane, being considerably larger when compared with the in-plane displacements. Furthermore, transverse shear rigidity of the slender plate is at least an order of magnitude higher than the bending rigidity. Therefore, the contribution of nonzero transverse shear strains to the deformation energy density of the fractional-order Kirchhoff plate may be neglected [45,48].
3.2 Numerical Modeling.
4 Results and Discussion
In this section, we report the numerical results for the models developed in Sec. 3 and provide a quantitative assessment of the effect of nonlocal interactions over the postbuckling response of Kirchhoff plates. More specifically, the effect of the fractional order and of the nonlocal length scales on the nonlinear load–displacement curves will be specifically addressed. In all simulations, we choose an isotropic solid with Young’s modulus E = 1.09 MPa and the Poisson’s ratio ν = 1/3 . The Young’s modulus chosen here applies to a class of soft materials that exhibit large deformations, which is essential for the intended study. However, note that the framework developed here is not specific to this particular choice of material properties.
The constitutive parameters corresponding to fractional-order continuum theory, which are the fractional-order α and the nonlocal length scale lf, are not fixed a priori. They will be varied to identify the effects of these individual parameters on the postbuckling response. Therefore, their value will be provided wherever required. The nonlocal length scales for the isotropic solid are chosen such that for points sufficiently away from the boundaries we have . This symmetry is broken for points close to the boundaries of the plate. As illustrated in Fig. 2, for points whose distance from the boundary (in a given direction) is smaller than lf, we truncate the nonlocal length scale at the physical boundaries of the plate. Concerning the numerical simulations carried out here below, we consider a square plate with in-plane dimensions a = b = 1 m. The thickness of the plate is chosen to be h = a/120 consistently with the assumptions for the Kirchhoff plate displacement theory.
Before proceeding with the analysis of the postbuckling behavior, we inspect the efficacy of the numerical model that was developed. Specifically, we establish the convergence of the 2D f-FEM and then validate the numerical model by comparison with existing results available in the literature.
To establish the convergence of the 2D f-FEM, we compare the postbuckling responses obtained from successive refinements of the finite element mesh. As seen in other studies on nonlocal elasticity, the numerical approximation of the system matrices includes an additional integration . In the case of the fractional-order strain–displacement relations employed here, this is the convolution integral present in the definition of the nonlocal strain (see Eq. (3)). The accuracy of the numerical approximation of this convolution integral is controlled by the “dynamic rate of convergence,” which is defined as . For the 2D f-FEM with convolution integrals defined in both the x1 and x2-directions, we have two parameters in (i = 1, 2) that control the element discretization in the xi-direction. The nonlinear load–displacement curves corresponding to different choices of this parameter are compared for the SSSS-01 plate in Fig. 4(a). From these results, we conclude that an excellent convergence of the postbuckling curves to within 1% difference is achieved for successive refinements beyond . A similar conclusion may be drawn for the convergence of the postbuckling curves for CSCS-01 plate from Fig. 4(b). On the basis of previous convergence results, we use for all the numerical simulations presented in the following section.
Next we validate the postbuckling behavior obtained from the current model when using α = 1.0 (i.e., the local case) with the corresponding results in the literature on classical isotropic plates. The efficacy of the 2D f-FEM as a numerical tool for fractional-order plates was established in our previous studies [46,48], so we validate the numerical method and the efficacy of the arc-length constraints directly on a study of postbuckling response. Note that, for the choice of α = 1.0, the fractional-order continuum theory identically reduces to classical elasticity theory [41,45]. The postbuckling curves for the isotropic square plate with α = 1.0 and subject to SSSS-01 boundary conditions are shown in Fig. 5(a). In this figure, the nonlinear load–displacement curves are compared for different amplitudes of the eccentric load . We validate these curves with the corresponding result available from the classical theory .
A similar exercise is carried on for the validation of the postbuckling curves for the CSCS-01 plate. The results presented in Fig. 5(b) show a comparison of the postbuckling response for the integer-order (α = 1.0) CSCS-01 plate with similar results available in the literature following the classical theory of elasticity . The direct comparison of the numerical results with those available in the literature for both boundary conditions shows a very good agreement. This convergence is clearly evident for load–displacement curves in the postcritical regime of the elastic response with the increasing load factor. In these results, we note that the load–displacement curves evaluated at α = 1.0 converge to the corresponding analytical results from classical elasticity. While this convergence is excellent for points away from the critical (turning) points, slight differences are noted for points in the neighborhood of critical points. This difference can be attributed to the sensitivity of the plate to imperfections resulting in lower values of the critical load for the imperfect structure. Unlike the current study where we consider an eccentric load F3, Yamaki  studied the bifurcation of a perfect structure subject to in-plane loads. Therefore, the case considered in the present study of a critical load of an isotropic plate subject to eccentric load is lower when compared to that of the perfect plate . Moreover, this difference increases with an increase in the amplitude of F3. Therefore, the lower load factor for critical (turning) points of plates subject to eccentric loads is in accordance with the theoretical predictions. Furthermore, as expected, this difference caused by imperfections vanishes for load–displacement curves away from the critical point. In conclusion, we establish a successful convergence of the numerical method and can proceed to use it for the study of the postbuckling response of fractional-order isotropic plates.
4.2 Postbuckling Response of Fractional-Order Plates.
In this section, we study the effect of nonlocal interactions on the postbuckling response of fractional-order plates. Specifically, we intend to analyze the effect of the fractional-order constitutive parameters: fractional-order α, and nonlocal length scale lf on (1) the critical point (turning point or bifurcation point) and (2) on both the postbuckling response and the postcritical stiffness of the structure. This analysis provides an understanding of the effects of the nonlocal behavior on the postbuckling response of thin plates. Note that either an increase in the nonlocal length scale lf or a reduction of the fractional-order α correspond to higher degree of nonlocality.
The nonlinear load–displacement curves of the isotropic nonlocal plate subject to uniaxial compression F1 and an eccentric load F3 for different values of the fractional-order α are compared in Fig. 6. The response for the SSSS-01 plate and CSCS-01 plate are presented in Figs. 6(a) and 6(b), respectively. These figures underline a softening response of the plate, beyond the critical point (turning point), when subject to an increasing degree of nonlocality (i.e., a progressive reduction of the fraction-order α). This observation is valid for both types of boundary conditions being considered. This observed softening in the postbuckling nonlinear response is consistent with the behavior expected from integral-type nonlocal constitutive relations , particularly the fractional-order continuum theories employed [39–41].
Unlike the previous comment on the effect of the nonlocality on the postbuckling response, contrasting observations can be drawn when focusing on the effect on the critical load (bifurcation point of the nonlinear load–displacement curves). In the case of the SSSS-01 beam depicted in Fig. 6(a), a marginal increase in the critical load is observed following the increasing degree of nonlocality. At the same time, a reduction in the critical load is noted with the reducing fractional-order α for the CSCS-01 plate. These observations for the two case studies are in complete agreement with the study on linear buckling of fractional-order structures . The fractional-order strain–displacement relations introduce fractional-order derivatives into the expressions of both the material and the geometric stiffness. This is clear from Eq. (10) for linear buckling where the fractional-order nonlocal terms and are present in both the numerator (material stiffness) and the denominator (geometric stiffness). An increase in the degree of nonlocality reduces the numerical values of the fractional derivatives. Thus, an increasing degree of nonlocality would result in lower material and geometric stiffness terms of the nonlocal structure. While a reduction in the material stiffness would reduce the critical load, a similar reduction in the geometric stiffness leads to higher values of the critical load. A comprehensive account of these effects, including detailed parametric studies on the effect of each stiffness term, is available in Ref. . Thus, we conclude that the effect of the nonlocal interactions on the critical load is dictated by the net result of the contrasting effects caused by the reduction of individual stiffness terms. This observation, drawn in the context of the linear buckling analysis of nonlocal structures, is also applicable to the current discussion on nonlinear buckling. Based on the eigenproblem for nonlinear buckling of fractional-order plates given in Eq. (9), similar conclusions may be drawn over the effect of fractional-order derivatives on the nonlinear material and geometric stiffness terms. Therefore, the effect of the nonlocal interactions on the nonlinear buckling load depends on the net result of the fractional-order terms within the two nonlinear stiffness terms.
It is known that the slope of the load–displacement curves indicate the stiffness of the structure . On the basis of the results presented in Fig. 6, we note a significant reduction in the slope of the nonlinear load–displacement curves for response beyond the critical point. This observation for α = 1.0 is corroborated by the existing literature on postbuckling behavior based on the classical theory of elasticity. In case of fractional values for α, the slope of the load–displacement curves corresponding to nonlocal elasticity is “flatter” when compared with the local analogs. This behavior confirms the lower stiffness value of nonlocal structures beyond the critical point. From the postbuckling nonlinear load–displacement curves, we also note that the influence of nonlocal elasticity increases by moving farther into the postbuckling regime. This is an important observation that points toward reduced resistance to deformation for nonlocal structures in a deep postbuckling regime.
The effect of the nonlocal length scale lf on the postbuckling response of fractional-order isotropic plates is explored in Fig. 7. This behavior is shown in Fig. 7(a) for the SSSS-01 plate and in Fig. 7(b) for the CSCS-01 plate. Similar to the aforementioned discussion regarding the effect of the fractional-order α, a softening behavior in the postbuckled regime is noted for increasing the values of the nonlocal length scale, lf. This observation is consistent for both boundary conditions chosen for this study. Furthermore, in agreement with the aforementioned results, contrasting observations are noted in terms of the effect of nonlocality on the critical load. For the SSSS-01 plate (Fig. 7(a)), the critical load increases with an increasing degree of nonlocality (higher lf). However, the critical load of the CSCS-01 plate (Fig. 7(b)) reduces for higher values of the length scale lf. These contrasting observations are consistent with the previous discussion on the influence of the fractional-order α on the critical load, as well as with previous work on linear buckling load of fractional-order plates .
Finally, we intend to undertake an analysis of the effect of nonlocal interactions on the sensitivity to imperfections of the plates. A detailed analysis can be proposed following the Koiter’s asymptotic approach previously carried out in Sec. 2.3 by including additional deformation energy terms corresponding to (geometric or load) imperfections in Eq. (7). However, for the sake of brevity, we present a quantitative analysis through the case study considered above. As discussed earlier, the transverse load F3 acting on the top surface of the plate amounts to a load imperfection on the isotropic plate subject to uniaxial compressive loads F1 applied the edges of the plate. Note that the role of the eccentric force F3 is to activate the stable branch of the load–displacement equilibrium curve beyond the critical point. This is because the fundamental curve C0 is no longer stable for equilibrium points beyond the critical point (see Sec. 2.3). Another approach in terms of geometric imperfection can be similarly envisioned. The distribution of F3 is motivated by the requirement to activate the transverse displacement mode corresponding to primary mode of buckling . A consequence of the imperfection would be to reduce the critical load corresponding to nonlinear buckling, commonly referred to as the imperfection sensitivity [2,56]. To compare this effect for nonlocal structures, we present the load–displacement curves of the CSCS-01 fractional-order plate for different amplitudes of the transverse force in Fig. 8. Irrespective of the amplitude of the transverse force, the postbuckling curves converge for load factor Λ away from the critical point. However, a lower nondimensionalized critical load (turning point) is noted for the increasing amplitude of F3 over the plate. This points to the imperfection sensitivity of the nonlocal plate. In this figure, we also compare the effect of nonlocal elasticity over the imperfection sensitivity for two different sets of fractional-order constitutive parameters. Any perceptible change in the critical load is not evident from the comparison between Figs. 8(a) and 8(b) for different choices of fractional-order and nonlocal length scales. Therefore, we conclude that postbuckling response of nonlocal plates exhibits imperfection sensitivity, which, however, is not influenced by the nonlocal interactions.
This article presented the theoretical approach and a numerical methodology for the analysis of the postbuckling elastic response of fractional-order structures. The approach is enabled by the thermodynamically consistent and positive semi-definite form of the deformation energy density that can be achieved via the fractional-order continuum theory for nonlocal structures. Further, geometrically nonlinear models of fractional-order strain–displacement relations are integrated in the formulation to enable the stability and postcritical response analyses. This theoretical approach opens the way to using energy methods to identify critical points of nonlinear equilibrium. In addition, it allows considering the influence of nonlocal interactions on both the material and the geometric stiffness terms. The possibility of using energy-based methods, afforded by the fractional formulation, allowed the application of Koiter’s asymptotic approach to ascertain the nature of the singularities exhibited by the nonlinear and nonlocal elastic response. Based on the fractional-order nonlinear stability theory, quantitative and qualitative studies were performed to understand the response of fractional-order (nonlocal) Kirchhoff plates in a postbuckling regime. The fractional-order finite element method (f-FEM) was extended to include an incremental-iterative solver and a path-following method to solve the nonlinear system of algebraic equations necessary to obtain the nonlinear equilibrium load–displacement curves. The effect of the nonlocal interactions on the nonlinear critical load is not trivial due to a simultaneous reduction of the material and geometric stiffness terms caused by the fractional-order kinematic relations. A reduction of the postbuckling stiffness due to the nonlocal interactions is noted, and it is associated with a softening response when compared with the local elastic analogs. Finally, it is established in literature that the frame-invariant fractional-order continuum theory serves as a mathematically, physically, and thermodynamically consistent alternative to the existing integer-order theories of nonlocal elasticity. Moreover, energy methods available for the fractional-order theory make it suitable to developing models that provide a better account of the scale-dependent effects on the stability and postbuckling response of complex structures. Following this, the current study opens an opportunity to perform stability of sandwich composites, layered and porous media, and biological materials like tissues and bones.
The following work was supported by the Defense Advanced Research Project Agency (DARPA) (Grant No. D19AP00052) and the National Science Foundation (NSF) (Grant Nos. DCSD #1825837 and MOMS #1761423). The content and information presented in this manuscript do not necessarily reflect the position or the policy of the government. The material is approved for public release; distribution is unlimited.
Conflict of Interest
There are no conflicts of interest.