Mechanics Based Design of Structures and Machines, 2025 (SCI-Expanded, Scopus)
This study aims to investigate the buckling behavior of axially functionally graded (AFG) columns with continuous elastic restraint within the framework of Euler–Bernoulli beam theory (EBT). The effects of elastic restraint stiffness, flexural rigidity variation, and boundary conditions on the buckling response of both homogeneous isotropic and AFG columns are examined. The variation of flexural rigidity along the column axis is modeled using a polynomial function. The influence of boundary conditions is investigated by examining both classical cases and a range of generalized combinations. The governing differential equations for buckling are reformulated as a first-order system suitable for the Complementary Functions Method (CFM), a numerical approach for boundary value problems (BVPs). The resulting initial value problems (IVPs) are solved using the fifth-order Runge–Kutta (RK5) method through a custom-developed Python code. The findings confirm the accuracy, efficiency, and applicability of the proposed approach for evaluating the buckling behavior of AFG columns with continuous elastic restraint under generalized boundary conditions.