Time domain simulation of Shapiro steps in Josephson junctions