Theoretical and numerical studies of self-steepening and shock formation of large-amplitude dust ion-acoustic waves in dusty plasmas are presented. A comparison is made between the nondispersive two fluid model, which predicts the formation of large-amplitude compressive and rarefactive dust ion-acoustic shocks, Vlasov simulations, and recent laboratory experiments.