When sliding a hard cylinder along the surface of glass, periodic surface cracks appear on the flat surface due to tensile stresses induced by the slider. These cracks propagate into the substrate and will affect the fracture properties of a body. Crack spacings and the directions of crack propagation into glass were calculated numerically by applying the finite element method and linear elastic fracture mechanics. The calculated crack spacings were in the range of the experimental results. Stress intensity factors and crack extension angles depended on the radius of slider and the load, and from these two factors the possible directions of crack propagation were calculated. The calculated propagation directions were in good agreement with real crack propagation.