In this work, we explore the spin ice model under uniaxial pressure using the Monte Carlo simulation method. For the known spin ices, the interaction correction (δ) introduced by the uniaxial pressure varies in quite a wide range from positive to negative. When δ is positive, the ground state characterized by the ferromagnetic spin chains is quite unstable, and in real materials it serves as intermediate state connecting the ice state and the long range ordered dipolar spin ice ground state. In the case of negative δ, the system relaxes from highly degenerate ice state to ordered ferromagnetic state via a first order phase transition. Furthermore, the domain walls in such ferromagnetic state are the hotbed of the excitations of magnetic monopoles, thus indicating that the uniaxial pressure can greatly increase the monopole density.