Knowledge of the detailed evolution of the whole charge state distribution of projectile ions colliding with targets is required in several fields of research such as material science, atomic and nuclear physics but also accelerator physics, and in particular in regards of the several foreseen large scale facilities. However, there is a lack for collisions in the non-perturbative energy domain and that involve many-electron projectiles. Starting from the ETACHA model we developed [Rozet et al. Nucl. Instrum. Meth. B 107, 67 (1996)], we present extension of its validity domain towards lower velocities and larger distortions. Moreover, the system of rate equations is able to take into account ions with up to 60 orbital states of electrons. The computed data from the different new versions of the ETACHA code are compared to some test collision systems. The improvements made are clearly illustrated by 28.9 MeV u-1 Pb 56+ ions, and laser-generated carbon ion beams of 0.045 to 0.5 MeV u-1 , passing through carbon or aluminum targets, respectively. Hence, those new developments can efficiently sustain the experimental programs that are currently under progress on the "next generation" accelerators or laser facilities.