The second-order stiffness matrix and corresponding loading vector of a prismatic beam–column subjected to a constant axial load and supported on a uniformly distributed elastic foundation (Winkler type) along its span with its ends connected to elastic supports are derived in a classical manner. The stiffness coefficients are expressed in terms of the ballast coefficient of the elastic foundation, applied axial load, support conditions, bending, and shear deformations. These individual parameters may be dropped when the appropriate effect is not considered; therefore, the proposed model captures all the different models of beams and beam–columns including those based on the theories of Bernoulli–Euler, Timoshenko, Rayleigh, and bending and shear.The expressions developed for the load vector are also general for any type or combinations of transverse loads including concentrated and partially nonuniform distributed loads. In addition, the transfer equations necessary to determine the transverse deflections, rotations, shear, and bending moments along the member are also developed and presented.