Long-term simulations of the rotational state of small irregular cometary nuclei