This paper is aimed at refinement of the computational model of the turbopump rotor systems associated taking into consideration the effect of rotation of moving parts and compliance of bearing supports elements. The up-to-date approach for investigation of non-linear reactions in rotor’s bearing supports is proposed for turbo-pump units for liquid rocket engines. Five models for modelling contact interaction are investigated, and comparative bearing stiffness characteristics are given. The geometry of the housing and corresponding design scheme are set for each support due to the assembly drawing of the turbopump unit. Rotation of the shaft is taking into account by applying corresponding inertial forces to the inner cage of the bearing. Experimental points of the dependence “load – displacement” as the diagram “F – v” are built by the calculated points as an array of numerical simulation data, obtained by the ANSYS software. As a result of numerical simulation, including loading of the bearing support on the scheme “remote force” in a wide range of rotor speeds, the corresponding displacements are determined. The brand-new approach for evaluation of bearing stiffness coefficients is proposed based on the linear regression procedure. As a result, the obtained values of coefficients are summarized and approximated by the quadratic polynomials.