Cytoplasmic dynein is a two-headed molecular motor that moves to the minus end of microtubule (MT) using ATP hydrolysis free energy. By employing its two heads (motor domains), cytoplasmic dynein shows various bipedal stepping motions; the inchworm and hand-over-hand motions, as well as non-alternate steps of one head. However, the molecular basis to achieve such diverse stepping manners remains obscure. Here, we propose a kinetic model for bipedal motions of cytoplasmic dynein and performed Gillespie Monte Carlo simulations that reproduces most experimental data obtained to date. The model represents status of each motor domain as five states according to conformations, nucleotide- and MT-binding conditions of the domain. Also, the relative positions of the two domains were approximated by three discrete states. Accompanied by ATP hydrolysis cycles, the model dynein stochastically and processively moved forward in multiple steps via diverse pathways, including inchworm and hand-over-hand motions, same as experimental data. The model reproduced key experimental motility-related parameters including velocity and run-length as functions of ATP concentration and external force. Our model reveals that, in a typical inchworm motion, the leading domain moves via the ATP-dependent power-stroke of the linker coupled with a small change in the stalk angle, whereas the lagging domain moves via diffusion dragged by the leading domain. Moreover, the hand-over-hand motion in the model dynein clearly differs from that of kinesin by the usage of the power-stroke.