A numerical scheme that incorporates a self-consistent cosmic-ray (CR, hereafter) injection model into the combined gas dynamics and CR diffusion-convection code has been developed. The hydro/CR code can follow in a very cos-effective way the evolution of CR modified shocks by adopting subzone shock-tracking and multi-level Adaptive Mesh Refinement techniques. The injection model is based on interactions of the suprathermal particles with self-generated MHD waves in quasi-parallel shocks. The particle injection is followed numerically by filtering the diffusive flux of suprathermal particles across the shock to upstream region according to a velocity-dependent transparency function, which represents the fraction of leaking suprathermal particles. In the strong shock limit of Mach numbers ≥ 20, significant physical processes such as the injection and acceleration seem to become independent of M, while they are sensitively dependent on M for M < 10. Although some particles injected early in the evolution continue to be accelerated to higher energies, the postshock CR pressure reaches a time asymptotic value due to balance between acceleration and diffusion of the CR particles.