A　seismic　response　analysis　of　layered,　liquefiable　sites　plays　an　important　role　in　the　seismic　design　of　both　aboveground　and　underground　structures.　This　study　presents　a　detailed　dynamic　site　response　analysis　procedure　with　advanced　nonlinear　soil　constitutive　models　for　non-liquefiable　and　liquefiable　soils　in　the　OpenSees　computational　platform.　The　stress　ratio　controlled,　bounding　surface　plasticity　constitutive　model,　PM4Sand,　is　used　to　simulate　the　nonlinear　response　of　the　liquefiable　soil　layers　subjected　to　two　seismic　ground　motions　with　different　characteristics.　The　nonlinear　hysteretic　behavior　of　the　non-liquefiable　soil　under　earthquake　excitations　is　captured　by　the　Pressure　Independent　Multi　Yield　kinematic　plasticity　model　with　a　von　Mises　multi-yield　surface.　The　soil　elements　are　modelled　utilizing　the　solid-fluid　fully　coupled　plane-strain　u-p　elements.　The　seismic　response　of　the　layered　liquefiable　site　in　terms　of　the　development　of　excess　pore　water　pressure,　acceleration,　ground　surface　settlement,　and　stress-strain　and　effective　stress　path　time　histories　under　two　representative　earthquake　excitations　are　investigated　in　this　study.　The　numerical　results　indicate　that　both　the　characteristics　of　ground　motions　and　the　site　profile　have　a　significant　influence　on　the　dynamic　response　of　the　layered　liquefiable　site.　Under　the　same　intensity　of　ground　motion,　the　loose　sand　layer　with　a　35%　relative　density　is　more　prone　to　liquefaction　and　contractive　deformation,　which　causes　irreversible　residual　deformation　and　vertical　settlement.　The　saturated　soil　layer　can　effectively　filter　the　high-frequency　components　and　amplify　the　low-frequency　components　of　ground　motions.　Moreover,　the　liquified　site　produces　a　40%　post-earthquake　consolidation　settlement　after　the　excess　pore　pressure　dissipation.