Visco-resistive magnetohydrodynamic (MHD) computations with the NIMROD code (Sovinec C R et al 2004 J. Comput. Phys.195 355) are applied to a model tokamak configuration that is subjected to induced vertical displacement. The modeling includes anisotropic thermal conduction within an evolving magnetic topology, and parameters separate the Alfvénic, resistive-wall, and plasma-resistive timescales. Contact with the wall leads to increasingly pervasive kink and tearing dynamics. The computed 3D evolution reproduces distinct thermal-quench and current-quench timescales, a positive bump in plasma current, and net horizontal forcing on the resistive wall. The MHD dynamo effect electric field, is analyzed for understanding the nonlinear effects of the fluctuations on the spreading of parallel current density and the resulting bump in plasma current. Forces on the resistive wall are consistent with Pustovitov's analysis (Pustovitov V D 2015 Nucl. Fusion55 113032); the plasma remains in approximate force-balance with the wall, so net force is accurately computed from integrating stress over the wall's outer surface. Improvements to the modeling that are needed for predictive simulation of asymmetric vertical displacement events are discussed.