Soil-structure interaction (SSI) is present in every problem, wherein a structural element is inrncontact with soil. A concrete beam supported directly by the soil continuum is a very commonrnconstruction form. The response of the beam when it carries external load is influenced byrnthe soil, and the response of the soil is also influenced by the action of the beam under thernload. Thus, developing a subgrade model for soil-structure interaction problem is essential inrnorder to predict the response of both components of the system and arrive at an optimumrndesign. The simplest subgrade model is the single parameter Winkler mechanical model whichrnrepresents the foundation soil by a series of independent springs. Winkler model is widelyrnused and practiced in spite of its deficiency in depicting the continuous behavior of real soils.rnLater, many advanced mechanical subgrade models have been proposed in order to improve onrnthe inherent lack of shear interaction among the individual springs. Moreover, these modelsrnstill have shortcomings with the nature of simplifying assumptions they make to ease thernmathematical relationships and not suggesting ways of estimating the model parameters. Withrnthe objective to improve on such drawbacks, a generalized continuum-based model has beenrnrecently proposed by Worku. This generalized subgrade model satisfies the fundamental elasticrnlaws and all the boundary conditions because it is derived by considering all stress, strain andrndisplacement components. The main objective of this work is to make a comparative studyrnbetween existing single and two parameter selected subgrade models. Microsoft Excel programsrnare written for the computation of deflection, moment and shear force. Finite Element basedrnPlaxis 2D is used as a tool to determine the calibration factor for generalized models. Lastly,rnnumerical illustration is provided using these models in comparison with the Plaxis 2D modelrnand other simplified models for finite and infinite beam subjected to selected symmetricalrnloading conditions. The results of the comparison show that the generalized models especiallyrnthe two-parameter model, are in very good agreement with the FE Plaxis 2D outputs.