The present article addresses the long-standing problem of the polymer brush bilayers under stationary shear flow at non-linear response regime where the system gets a non-Newtonian fluid. The main idea behind this research would be the fact that the immense lubricity of the polymer brush bilayers originates from a global restructuring that takes place at large shear rates. It is shown here that physical quantities like, stress tensor, viscosity tensor, the friction coefficient and the chain extensions could become dependent on the shear rate at non-Newtonian regime. Apparently, the sub-linear scaling of the physical quantities at large shear rates is solely due to the fact that the chains stretch in the shear direction.