In this paper, forces and torques on solid, non-spherical, orthotropic particles in Stokes flow are investigated by using a numerical approach on the basis of the Boundary Element Method. Different flow patterns around a particle are considered, taking into account the contributions of uniform, rotational and shear flows, to the force and the torque exerted on the particle. The expressions for the force and the toque are proposed, by introducing translation, rotation and deformation resistance tensors, which capture each of the flow patterns individually. A parametric study is conducted, considering a wide range of non-spherical particles, determined by the parametric superellipsoid surface equation. Using the results of the parametric study, an approximation scheme is derived on the basis of a multivariate polynomial expression. A coefficient matrix for the polynomial model is introduced, which is used as a tunable parameter for a minimization problem, whereby the polynomials are fitted to the data. The developed model is then put to the test by considering a few examples of particles with different shapes, while also being compared to other, currently available solutions. On top of that, the full functionality of the model is demonstrated by considering an example of a pollen grain, as a realistic non-spherical particle. First, a superellipsoid, which best fits the actual particle shape, is found from the considered range. After that, the coefficients of the translation, rotation and deformation resistance tensors are obtained from the present model and compared to the results of other available models. In the conclusion, a superior accuracy of the present model, for the considered range of particles, is established. To the best of the authors knowledge, this is also one of the first models to extend the torque prediction capabilities beyond sphere and prolate particles. At the same time, the model was demonstrated to be simple to implement and very conservative with the computational resources. As such, it is suitable for large scale studies of dispersed two-phase flows, with a large number of particles.