A three-dimensional (3D) thermoelastic stress analysis is carried out on a single crystal with axisymmetric geometry but with a hexagonal crystallographic symmetry. The crystallographic orientation is off-axis with respect to the cylindrical coordinate system. By applying a Fourier series expansion with respect to the rotational angle φ of the cylindrical coordinates, the 3D boundary value problem is reduced to a sequence of 2D ones on the meridian plane, which are solved by the finite-element method. In our example, the off-axis orientation is towards a direction of high symmetry, and therefore only four of the six stress tensor components are non-zero. In the end, the stress tensor is projected onto the slip system of the crystal.