Abstract Body force modeling is studied in the Generalized Newtonian (GN) fluid flow simulation using a single relaxation time lattice Boltzmann (LB) method. First, in a shear thickening Poiseuille flow, the necessity for studying body force modeling in the LB method is explained. Then, a parametric unified framework is constructed for the first time which is composed of a parametric LB model and its associated macroscopic dual equations in both steady state and transient simulations. This unified framework is used to compare the macroscopic behavior of different forcing models. Besides, using this unified framework, a new forcing model for steady state simulations is devised. Finally, by solving a number of test cases it is shown that numerical results confirm the theoretical arguments presented in this paper.