A three-dimensional (3D) full-dynamic damage model for ion implantation in crystalline silicon was proposed to calculate more accurately point defect distributions and ion-implanted concentration profiles during ion implantation process. The developed model was based on the physical monte carlo approach. This model was applied to simulate B and BF2 implantation. We compared our results for damage distributions with those of the analytical kinchin-pease approach. In our result, the point defect distributions obtained by our new model are less than those of kinchin-pease approach, and the vacancy distributions differ from the interstitial distributions. The vacancy concentrations are higher than the interstitial ones before 0.8 . Rp to the silicon surface, and after the 0.8 . Rp to the silicon bulk, the interstitial concentrations are revesrsely higher than the vacancy ones.The fully-dynamic damage model for the accumulative damage during ion implantation follows all of the trajectories of both ions and recoiled silicons and, concurrently, the cumulative damage effect on the ions and the recoiled silicons are considered dynamically by introducing the distributon probability of the point defect. In addition, the self-annealing effect of the vacancy-interstitial recombination during ion implantation at room temperature is considered, which resulted in the saturation level for the damage distribution.