Abstract: Current calculation processes on dynamic derivatives of projectiles and the corresponding stability analysis are usually conducted with the influences of the motion coupling effect ignored. However, the pitching, rolling and coning motions exist simultaneously in most cases. In this paper, a method based on the Euler rotation theorem and spherical sliding meshing technique has been proposed to deal with the simulation of complex angular motion. The modification value of the projectile's angular velocity at each time step is interpolated out and the sliding meshing zone is allocated with the application of the Rodriguez transformation matrix in the method. The influences of the projectile's rolling and coning motion on pitching combination dynamic derivatives and lift coefficient are analyzed after solving and identifying the unsteady aerodynamic parameters. The results indicate that the proposed method can effectively eliminate the accumulative errors in calculation of attitude angles and realize the accurate simulation under arbitrary given angular motion. For the coupling angle movement, the hysteresis loop of aerodynamic loads are greatly oscillated and offset, and the pitching combination dynamic derivatives and lift coefficient change significantly and nonlinearly with the increase of roll frequency and taper frequency.