In this work, the effect of uniaxial stress on the stability of magnetic skyrmion is studied by use of Monte Carlo simulations. The classical Heisenberg model with anisotropic exchanges and Dzyaloshinskii-Moriya interactions is applied to correlate the stress behavior. We calculate magnetization M and skyrmion number Q to estimate the phase transition points, use the spin configurations and their Fourier transform to determine the region of skyrmion state. The h − T phase diagram show the area of skyrmion state is larger with stress than without stress. These results are interpreted qualitatively by the energy competitions between Zeeman, exchanges and Dzyaloshinskii-Moriya interactions. The simulated results suggest that anisotropy induced by uniaxial stress can stabilize the skyrmion phase.