Dental implants have substantial significance in modern dentistry, but their osseointegration stage remains the most crucial and time-consuming step for compensating for a lost tooth. One promising approach to improve osseointegration rates is the use of extracorporeal shock wave therapy, which has shown efficacy in treating fractures, bone defects, and bone tissue regeneration in surgical and arthroplasty procedures. To comprehend the potential of shock wave therapy in accelerating implant osseointegration, it is crucial to investigate the influence of mechanical loading on ossification processes of varying scales objectively. This study aims to study numerically the effects of low-energy shock wave therapy at different intensities on the mechanical response of a single osteocyte, the principal bone cell that corresponds to microscale of bone tissue. The investigation employs the method of movable cellular automata for modeling. The computer simulation results and analysis based on mechanobiological principles indicate that low-intensity shock wave loading creates conditions for intramembranous ossification, whereas high-intensity shock wave exposure creates conditions for endochondral ossification.