The rotary vector reducer (RV reducer) is a key reduction mechanism for industrial robots, and its torsional stiffness directly determines positioning accuracy and dynamic characteristics. Unlike traditional time-varying torsional stiffness analysis methods considering limited factors and statics solution method, a time-varying torsional stiffness analysis method considering output torque variability and dynamic nonlinearity is proposed to assess the time-varying torsional stiffness of the reducer. Firstly, the rigid-flexible coupling multi-body dynamic model was built employing software ADAMS, with crank shaft and cycloidal gears being treated as flexible bodies. Its correctness was verified under standard condition through grey correlation theory and the torsional stiffness under different conditions was obtained, respectively. Secondly, a contact finite element model of the reducer was established, taking contact relationships between different shaft-gear-bearings into account. As a comparison, a rotating-join finite element model was also constructed. And their torsional stiffness under experimental condition was calculated. Finally, the torsional stiffness test was conducted to verify the accuracy of the three models under this condition. The results show that the time-varying torsional stiffness obtained by the proposed method is only 8.83% different from the test results and exhibiting a certain periodic pattern under standard conditions.