Abstract: An analytical method for the stability analysis of beams on the layered Gibson soil is developed. The modified Vlasov model is adopted to simulate the foundation behavior. Using the variational principle, the governing equation of the beam and the relation between parameter  and buckling mode shape are derived. The buckling loads and the corresponding mode shapes of the beam can be determined by solving the eigenvalues of the boundary conditions of the beam. Parametric study shows that the soil properties near the surface have larger influence on the vertical elasticity and continuity of the foundation than those near the bottom. The inhomogeneity of the foundation affects the buckling loads and corresponding mode shapes greatly. The contribution of the resting foundation near the surface to the buckling load is more significant than those near the bottom.