Humanoid robots are expected to operate in unstructured environments where foot-soil interaction strongly affects mobility. This work develops a simulation framework that couples multibody dynamics with deformable terrain models to study locomotion under such conditions. We compare two approaches: the Discrete Element Method (DEM), which provides particle-level accuracy at high computational cost, and the Soil Contact Model (SCM), an empirical formulation calibrated through DEM-based virtual bevameter tests. Through a progression of experiments–from single-foot validation to full-robot walking, climbing, running, and jumping–we show that SCM predicts foot-soil forces and joint loads with accuracy sufficient for motion planning, while achieving near real-time efficiency. DEM, while more computationally demanding, remains essential for analyzing extreme maneuvers involving rapid soil deformation. Together, these results highlight a pathway for scalable, high-fidelity simulation of humanoid locomotion in granular environments, with direct implications for planetary exploration, disaster response, and other field robotics applications.