The control of polymerization conditions and final properties of polybutadiene produced through solution polymerizations constitute an important industrial challenge, which can be tackled with help of mathematical models. For this reason, the present work focus on the development of a phenomenological mathematical model, and estimation of the respective model parameters, in order to describe the kinetics of 1,3‐butadiene solution polymerizations catalyzed by neodymium versatate. The proposed model can be used to predict some important final properties of polybutadiene resins, including the cis‐content and average molecular weights, and the evolution of operation conditions that are used at plant site to monitor the course of the reaction, including the reaction temperature and pressure. Given the low required computation times, the proposed model can be used for purposes of monitoring and control in actual industrial sites. POLYM. ENG. SCI., 59:2290–2300, 2019. © 2019 Society of Plastics Engineers