Calculations of vertical shear stresses and the distribution of shears in the slabs in post-elastic range obtained by using a two-dimensional nonlinear plate bending finite element programme is presented. The shears at each Gaussian sampling point within each element, in x and y-direction were calculated by using a set of equilibrium equations.