What probability measure(s) maximize the quantity $\iiint\_{\mathbb{S}^1}|(x-z)\times(y-z)|d\mu(x)d\mu(y)d\mu(z)$? 


The answer appears to be uniform measure, since informally it appears better to have more triangles in the support of $\mu$ which the function $|(x-z)\times(y-z)|$ computes the area of. 


Is there an argument that shows this is true generally for $\mathbb{S}^{d-1}$?


**Edit:** Thanks to @fedja for the answer to the main question and the request for clarification. While it could be possible to replace the area of triangles with volumes of simplices (and so on), what was meant by generalization was most immediately the question of whether uniform measure maximizes $$\iiint\_{\mathbb{S}^2}|(x-z)\times(y-z)|d\mu(x)d\mu(y)d\mu(z).$$



Maximizing $\iiint|(x-z)\times(y-z)|d\mu d\mu d\mu$ over probability measures on the unit circle
Yes, it is true for the circle (though the reason is not quite the one you suggested). We shall consider the discrete version of the problem, which is to put some odd number $n\ge 3$ of points (think of them as of being assigned the probability of $1/n$ each) on the circle so that the sum of triangle areas is maximized. Then an optimal configuration exists by compactness. We just want to show that it should be an equispaced distribution. Since you can obtain any measure on the circle as a weak limit of such discrete measures (and the integrand is continuous, so a small displacement of the measure does not change the integral much), you'll get the required optimality of the uniform measure as a limiting statement.


Take any of the points $e$ and enumerate all other points by the numbers $1,\dots,n-1$ clockwise starting from $e$ (so the whole system is $e,e\_1,\dots,e\_{n-1}$. Fix all points except $e$ and think a bit of how the total double area depends on the position of $e$. You'll realize quite soon that it is just some constant plus $v\_e\times e$ where 
$$v\_e=\sum\_{1\le k< m \le n-1}(e\_m-e\_k).$$


If you now consider for fixed $\ell<\frac n2$ all vectors $e\_m-e\_k$ with [$(\ell\le k< m\le n-\ell$) and ($k=\ell$ or $m=n-\ell$)], you'll see that $v\_e=\sum\_{\ell<\frac n2}A\_\ell(e\_{n-\ell}-e\_\ell)$ with some $A\_\ell>0$. You will be able to increase the sum of (even oriented) areas if $e$ is not perpendicular to $v\_e$, so in the optimal configuration we must always have $\langle e,v\_e\rangle=0$ for every choice of $e$.


Now recall the famous "drive around the circle without running out of gas" problem (more precisely, its particular case when all gas stations have the same amount of fuel sufficient for driving $\frac 1n$ of the entire circle). The conclusion of it is that you can choose $e$ so that the arc distance to $e\_k$ is at most $\frac k n$ of the whole circle for every $k=1,\dots,n-1$, which makes all scalar products $\langle e,e\_{n-\ell}-e\_\ell\rangle$ with $\ell<\frac n2$ non-positive with the only chance of the equality $\langle e,v\_e\rangle=0$ when the arc distance from $e$ to $e\_k$ is exactly $\frac kn$ for all $k=1,\dots,n-1$.


As to higher dimensions, can you, please, specify first what you mean by "this"? (there are several ways to try to generalize the setup, so the word seems rather ambiguous to me :-) )


