A semi-analytical model of the Yarkovsky-O''Keefe-Radzievskii-Paddack (YORP) effect on an asteroid spin in a non-principal axis rotation state is developed. The model describes the spin-state evolution in Deprit-Elipe variables, first-order averaged with respect to rotation and Keplerian orbital motion.
Assuming zero conductivity, the YORP torque is represented by spherical harmonic series with vectorial coefficients, allowing us to use any degree and order of approximation. Within the quadrupole approximation of the illumination function we find the same first integrals involving rotational momentum, obliquity and dynamical inertia that were obtained by Cicaló & Scheeres.
The integrals do not exist when higher degree terms of the illumination function are included, and then the asymptotic states known from Vokrouhlický et al. appear. This resolves an apparent contradiction between earlier results.
Averaged equations of motion admit stable and unstable limit cycle solutions that were not previously detected. Non-averaged numerical integration by the Taylor series method for an exemplary shape of 3103 Eger is in good agreement with the semi-analytical theory.