We introduce a Newtonian model for the deformations of a compressible, autogravitating, and continuously stratified neutron star. The present framework can be applied to a number of astrophysical scenarios as it allows us to account for a great variety of loading forces. In this first analysis, the model is used to study the impact of a frozen adiabatic index in the estimate of rotation-induced deformations: we assume a polytropic equation of state for the matter at equilibrium but, since chemical reactions may be slow, the perturbations with respect to the unstressed configuration are modelled by using a different adiabatic index. We quantify the impact of a departure of the adiabatic index from its equilibrium value on the stressed stellar configuration and we find that a small perturbation can cause large variations both in displacements and strains. As a first practical application, we estimate the strain developed between two large glitches in the Vela pulsar showing that, starting from an initial unstressed configuration, it is not possible to reach the breaking threshold of the crust, namely to trigger a starquake. In this sense, the hypothesis that starquakes could trigger the unpinning of superfluid vortices is challenged and, for the quake to be a possible trigger, the solid crust must never fully relax after a glitch, making the sequence of starquakes in a neutron star an history-dependent process.

We study a class of Newtonian models for the deformations of non-magnetized neutron stars during their spin-down. The models have all an analytical solution, and thus allow to understand easily the dependence of the strain on the star's main physical quantities, such as radius, mass and crust thickness. In the first ``historical'' model the star is assumed to be comprised of a fluid core and an elastic crust with the same density. We compare the response of stars with different masses and equations of state to a decreasing centrifugal force, finding smaller deformations for heavier stars: the strain angle is peaked at the equator and turns out to be a decreasing function of the mass.We introduce a second, more refined, model in which the core and the crust have different densities and the gravitational potential of the deformed body is self-consistently accounted for. Also in this case the strain angle is a decreasing function of the stellar mass, but its maximum value is at the poles and is always larger than the corresponding one in the one-density model by a factor of two. Finally, within the present analytic approach, it is possible to estimate easily the impact of the Cowling approximation: neglecting the perturbations of the gravitational potential, the strain angle is 40\% of the one obtained with the complete model.