arXiv:2001.08631v2 [cond-mat.soft] 5 Mar 2020 Lerner, Düring, and Wyart, 2012]. The value of the exponent ν is debated [Lerner, Düring, and Wyart, 2012;Gallier et al., 2014a;Mari et al., 2014;Ness and Sun, 2015], the most typical value for spherical particles in the literature being ν = 2 Hermes et al., 2016;Guy et al., 2018]. For long rods, Tapia et al. [2017] measure ν = 1. Shear thickening stems from the fact that φ J , is at a higher value, φ 0 J , for frictionless particles than for frictional ones, at φ 1 J < φ 0 J . Because frictional contacts only exist at large stresses, due to the fact that under small stresses repulsive forces prevent them to form, the suspension has a stress dependent jamming point, φ J (σ), such that η is an increasing function of σ.This scenario is a priori independent of the particle shape. While it has been tested mostly with suspensions of spherical particles, it is expected to be equally valid for suspensions of nonspherical particles, and many such suspensions are known to shear thicken in a qualitatively same way than spherical particles. Cornstarch suspension is the most famous example [Brown and Jaeger, 2009;Fall et al., 2012;Oyarte Gálvez et al., 2017], but suspensions of synthetized particles were also studied [Egres and Wagner, 2005;Brown et al., 2011;Royer et al., 2015;James et al., 2019;Rathee et al., 2019]. Most industrial suspensions known to shear thicken, such as cement paste [Lootens et al., 2004;Papo and Piani, 2004;Feys, Verhoeven, and De Schutter, 2009;Toussaint, Roy, and Jézéquel, 2009;Roussel et al., 2010], suspensions used for mechanical polishing such as fumed silica [Crawford et al., 2012;Amiri, Øye, and Sjöblom, 2012] or quartz powder suspensions [Freundlich and Roder, 1938], fresh paints and coatings [Zupančič, Lapasin, and Žumer, 1997;Khandavalli and Rothstein, 2016], or molten chocolate [Blanco et al., 2019], also contain particles of varied non-spherical shapes. Numerical simulations of suspensions of frictional and repulsive aspherical particles also showed shear thickening [Lorenz et al., 2018].The major difference with spherical particles comes from the different values for φ 0 J and φ 1 J (although other differences exist, e.g. the exponent of the viscosity divergence close to jamming [Tapia et al., 2017]). Indeed, it is known that the jamming point for isotropic random packings generically depends on the shape of the particles [Torquato and Stillinger, 2010;Jiao and Torquato, 2011;Baule and Makse, 2014]. For families of axisymmetric shapes (like axisymmetric ellipsoids, rods, spherocylinders, etc), which can be characterized by one scalar value, the aspect ratio α (the ratio between particle length and width), the behavior is non-monotonic as a function of α. Starting from spheres (α = 1), the jamming volume fraction generally increases up to a maximum reached in the α = 1.2 − 2 range, and then decreases with increasing α for larger aspect ratios [A similar orientation is found in dense suspensions [Egres and Wagner, 2005;Rathee et al., 2019], alt...