(More a comment that is too long for the comment box) ...  Given a sample of size $n$ from a standard Uniform parent, one can show, by induction or otherwise, that the joint pdf of the first $k$ order statistics, is:

 $$f(x_{(1)}, x_{(2)}, \dots, x_{(k)}) = \frac{n!}{(n-k)!}\left(1-x_{(k)}\right)^{n-k} \quad  \quad \text{for }\quad  0 < x_{(1)} < x_{(2)} < \dots <x_{(k)} < 1$$

... which has a neat functional form.

Your problem then reduces to:

> Find the cdf of the sum $X_{(1)} + X_{(2)} + \dots + X_{(k)}$ given joint pdf $f(\centerdot)$ 

... which unfortunately will result in a $k$-part piecewise solution with kinks at $(1, 2, \dots, k)$. 

Using the _mathStatica_ package for  _Mathematica_, I obtained the same solution you provided for $k=2$, and was also able to derive a general exact solution for the sum in the $k = 3$ case (in 3 parts). However, an exact general solution seems to get messy quite rapidly.

For illustration, here is a plot of the cdf of the sum in the $k =2$ case (as $n$ increases from 2 to 10, and then with $n =30$):
[![enter image description here][1]][1]

The intuition, of course, is that when $n$ becomes large, the first two order statistics become very small, and so the cdf of their sum shifts to the left.

Here is the $k=3$ case:

[![enter image description here][2]][2]


  [1]: https://i.sstatic.net/DNmil.png
  [2]: https://i.sstatic.net/2Qg0y.png