A new, easy-to-implement method for the free vibration analysis of Euler-Bernoulli beam-columns embedded in a two-parameter elastic soil is presented. The non-homogeneity of the soil was accounted for by assuming a linear distribution of the modulus of subgrade reaction. The governing differential equation (GDE) of the proposed structural element is solved using the differential transform method (DTM). This method is implemented since finding an analytical solution by means of other traditional approaches is a complex task and becomes even more challenging when the non-homogeneity of the soil is introduced in the formulation. Examples are presented to validate the proposed formulation. The results show that the introduction of a shear layer and soil homogeneity has a strong influence on the natural frequencies of rigid and semirigid beam-columns. Comparisons are performed to demonstrate the simplicity and accuracy of the proposed method.