The simulation of QCD on dynamical (Nf=2) anisotropic lattices is described. A method for nonperturbative renormalisation of the parameters in the anisotropic gauge and quark actions is presented. The precision with which this tuning can be carried out is determined in dynamical simulations on 8^3x48 and 8^3x80 lattices.