...
Here we note that the principal normal stresses σ1 = σ2 = σr = σθ since the out-of-plane shear stresses, τrz = τθz = 0 and σ3 = σz. And we can deduce that τmax = |τ1|=|τ2|=|(σ1-σ2) / 2|. The effective stress (using the Von-Mises criterion) along the y-axis can be computed as the following:
Lastly, we also confirm that the applied load at the top vertex of the sphere matches our numerical contact pressure, integrated along the interface.
...