Transverse stress

How is it possible to obtain the transverse shear stresses correctly, not constant, by Reissner Mindlin theory?

Equilibrium equation would be the first solution, but, once the second derivatives of bilinear shape functions are zero, no results is appeared.

Thanks in advance for any help!