Here is a more visual approach: Fill a square $n\times n$ grid with the values of $i+j$:
$$
\begin{matrix}
1+1&1+2&1+3&\cdots&1+n\\
2+1&2+2&2+3&\cdots&2+n\\
3+1&3+2&3+3&\cdots&3+n\\
\vdots&\vdots&\vdots&\ddots&\vdots\\
n+1&n+2&n+3&\cdots&n+n
\end{matrix}
$$
We are after the sum of all the elements above the main diagonal (the diagonal that goes from $1+1$ to $n+n$). Note that this is equal to the sum of all the elements below the main diagonal. So if we can find the sum of the entire table, and the sum of the main diagonal, then a simple subtraction followed by a division by 2 will get us there.
The sum of the main diagonal is $n(n+1)$.
For the sum of the entire table, collect each $j$ column-by column, and all the $i$'s row-by-row. The $j$'s in the first column sum to $n$. The $j$'s in the second column sum to $2n$, and so on. And similarly for the $i$'s. Thus the sum of all the elements in the table is
$$
\underbrace{n+2n+\cdots+n^2}_{\text{the }j\text{'s}}+\underbrace{n+2n+\cdots+n^2}_{\text{the }i\text{'s}}=2n(1+2+\cdots+n)\\
=n^2(n+1)
$$
Our final answer is therefore
$$
\frac{n^2(n+1)-n(n+1)}2=\frac12n(n^2-1)
$$