...
Note here that the stress in x direction is not constant as assumed in the analytical method. The point load at the end of the bar means the stress will be greater on the top and bottom of the bar. There are also higher stresses in the area close to the wall.
Sigma_y:
...
The analytical method assumed a long bar, therefore by definition, the stresses in the y direction are assumed to be zero. From the simulation, it is clear to see this is obviously not the case.
Tau_xy:
Von Mises:
...