A new beam suspension model for load/unload simulations of magnetic hard disk drive sliders is presented. The model is applicable to the load/unload process as well as to the dynamic flying behavior of the slider on the disk. Suspension and gimbal are modeled with C1 continuous beam elements. The roll degree of freedom of the slider is decoupled from spacing and pitch. The numerical implementation of the model is described. Loading and unloading of a typical subambient pressure type slider is presented.