How can we generalize factorisation of $(a+b)^n-(a^n+b^n)$ 
How can we generalize factorisation of
$$(a+b)^n-(a^n+b^n)\,?$$
where $n$ is an odd positive integer.

I found the following cases:
$$(a+b)^3-a^3-b^3=3ab(a+b)$$
$$(a+b)^5-a^5-b^5=5 a b (a + b) (a^2 + a b + b^2)$$
$$(a+b)^7-a^7-b^7=7ab(a+b)(a^2+ab+b^2)^2$$
Exapansions and factoring seem terrible. Maybe there is a good way for generalization.
The question is, can we generalize these factorisations?
$$(a+b)^n-(a^n
+b^n)$$
for all odd integer $n$?

There doesn't seem to be any general factored form, but on the few polynomial factors you mention in the question, we can answer exactly when they appear. From symmetry and homogeneity, it is enough to consider the simpler $$p(x)=(1+x)^n-(1+x^n)$$ as with $x = b/a$, we have $a^n p(x)$ is the expression you want to factor.

*

*Now $p(0)=0 \implies x$ is always a factor of $p(x) \implies ab$ is always a factor of the original bivariate expression.  From the Binomial theorem, it is also clear that no higher power of $x$ or $(ab)$ will be a factor.

*Noting $p(-1)= -1-(-1)^n \implies x+1$ is a factor when $n$ is odd, or equivalently $(a+b)$ is a factor for odd $n$.  From $p'(x) = n(1+x)^{n-1}-nx^{n-1}, p'(-1) \neq 0$ so the multiplicity of this factor is always $1$.

*Let $\omega$ be a non-real cube root of unity, i.e. $\omega^2+\omega+1=0$.  Then
$p(\omega) = (-1)^n(\omega)^{2n}-1-\omega^n$.
Case $n=3k$: we have $p(\omega)=(-1)^n$
Case $n=3k+1$: we have $p(\omega)= (-1)^n\omega^2+\omega^2$, which is zero iff $n$ is also odd, i.e. for $n=6k+1$.
Case $n=3k+2$: we have $p(\omega)= (-1)^n\omega+\omega$, which vanishes iff $n$ is also odd, i.e. for $n=6k+5$.

If $p(\omega)=0$ we have $1+x+x^2$ as a factor, or equivalently $a^2+ab+b^2$ is a factor iff $n \equiv \pm1 \pmod 6$. Let us also check multiplicity for these cases using $p'(x) = n(1+x)^{n-1}-nx^{n-1} \implies p'(\omega)=n\left( (-\omega^2)^{n-1} - \omega^{n-1}\right)$:
Case $n=6k+1$: we always have $p'(\omega) = 0$, so always $(a^2+ab+b^2)^2$ is always a factor.  Also $p''(\omega) = n(n-1) \neq 0$, so the multiplicity is exactly $2$.
Case $n=6k+5$: we have $p'(\omega) = n\left(\omega^2-\omega \right)\neq 0$ so here the multiplicity is always $1$.
Summarising all of that, for all odd $n$, we have $ab(a+b)$ as a factor, with no higher power for any of the terms.  In case $n=1 \pmod 6$, we have additionally $(a^2+ab+b^2)^2$ as a factor (with no higher power possible) and in case $n=5 \pmod 6$ we have just $a^2+ab+b^2$ as a factor.  In no other case will $a^2+ab+b^2$ be a factor.
The above few factors with limited multiplicity are clearly not enough to provide required degree for higher order polynomials, so there would necessarily be other polynomial factors.
