We explore the quadratic form of the f(R)=R+bR2 gravitational theory to derive rotating N-dimensions black hole solutions with ai,i≥1 rotation parameters. Here, R is the Ricci scalar and b is the dimensional parameter. We assumed that the N-dimensional spacetime is static and it has flat horizons with a zero curvature boundary. We investigated the physics of black holes by calculating the relations of physical quantities such as the horizon radius and mass. We also demonstrate that, in the four-dimensional case, the higher-order curvature does not contribute to the black hole, i.e., black hole does not depend on the dimensional parameter b, whereas, in the case of N>4, it depends on parameter b, owing to the contribution of the correction R2 term. We analyze the conserved quantities, energy, and angular-momentum, of black hole solutions by applying the relocalization method. Additionally, we calculate the thermodynamic quantities, such as temperature and entropy, and examine the stability of black hole solutions locally and show that they have thermodynamic stability. Moreover, the calculations of entropy put a constraint on the parameter b to be b<116Λ to obtain a positive entropy.