Seismic liquefaction, a phenomenon which occurs during earthquakes and causes enormous damage and failures of foundations and structures. Based on documented records of liquefaction case histories, different procedures have been developed to assess the liquefaction potential and established suitable liquefaction countermeasures. Site-specific geologic details, ground motion characteristics and structural configurations cannot be accounted for effectively in the empirical and semi-empirical approaches of liquefaction potential assessment. This paper presents the use of finite element analysis (FEA) as a computational tool using advanced constitutive relations to assess the liquefaction triggering and its associated consequences. Using data from the cyclic direct simple shear (CDSS) tests of the foundation soil, the UBC3D-PLM model is calibrated and further validated using the results of centrifuge test available in the literature. A plane strain condition in PLAXIS 2D is used to model the earthen embankment resting on a liquefiable soil deposit. The UBC3D-PLM model is used to predict the triggering of liquefaction and subsequent increase in excess pore pressure (EPP) during earthquake loading. Emphasis is laid on simulating the different factors affecting the liquefaction such as the effect of initial effective confining stress, load amplitude and initial static shear stress. A flow failure of the embankment is also simulated. Because of liquefaction triggering, there is a sudden change in the ground motion frequency and amplitude. This effect is captured by adopting the tied degrees of freedom (TDOF) as lateral boundary conditions in the FEA. It is concluded that the finite element method can simulate the complex phenomenon of liquefaction triggering and model its consequences more effectively. These simulation results can be used to implement liquefaction countermeasures in the field.