  The detailed formulation for loop quantum cosmology (LQC) in the Bianchi I
model with a massless scalar field was recently constructed. In this paper, its
effective dynamics with the LQC discreteness corrections is studied and the
equations of motion are analytically solved, showing that the big bang
singularity is replaced by the big bounces, which take place up to three times,
once in each diagonal direction, whenever each of the area scale factors
approaches its critical value in the Planck regime measured by the reference of
the scalar field momentum.
