By means of molecular dynamics simulations we demonstrate power laws for macroscopic transport properties of strongly compressed polymer-brush bilayers to stationary shear motion beyond the Newtonian response. The corresponding exponents are derived from a recently developed scaling theory, where the interpenetration between the brushes is taken as the relevant length scale. This allows to predict the dependence of the critical shear rate, which separates linear and non-linear behavior, on compression and molecular parameters of the bilayer. We present scaling plots for chain extension (R), viscosity (η) , and shear force (F over a wide range of Weissenberg numbers, W . In agreement with our theory, the simulation reveals simple power laws, R ∼ W (0.53), η ∼ W (-0.46), and F ∼ W (0.54), for the non-Newtonian regime.