Axisymmetrical computation of the stress distribution in orthotropic rotating discs