New phenomenological models are proposed to describe the effect of an ordered lattice structure of crystalline targets on the as-implanted doping profiles of low-energy heavy ions. The models account for the channeling kinetics and clarify the effect of bi-directional transitions of ions between random-like and channeled modes of motion on the target depth dependencies of dopant concentration. They also incorporate a simple model of the target radiation damaging effect on doping profiles. The presented results of model validation against the experimental and Monte Carlo computer simulation data and the comparative analysis of the capabilities of the proposed and the existing models show that the application of a more physically grounded approach allows us to improve the quality of doping profile description. The theoretical models developed are useful for obtaining physical parameters of low-energy ion channeling kinetics from the experimental data.