The non-linear dynamics of irrotational dust in General Relativity is studied in synchronous and comoving coordinates. All the equations are written in terms of the metric tensor of spatial sections orthogonal to the flow, which allows an unambiguous expansion in powers of $1/c^2$. To lowest order, the Newtonian approximation in Lagrangian form is derived. At this level the evolution is governed by the Raychaudhuri equation for the Lagrangian-to-Eulerian Jacobian matrix. The Lagrangian spatial metric reduces to that of Euclidean 3-space in time-dependent curvilinear coordinates. A Lagrangian version of the Bernoulli equation for the evolution of the `velocity potential' is also given. At the post-Newtonian (PN) level, an exact and general formula is derived for gravitational-wave emission from non-linear perturbations. It is shown that, in the anisotropic collapse of homogeneous ellipsoids, the ratio of the PN tensor modes to the Newtonian metric tends to diverge like the mass density. It is finally argued that a stochastic gravitational-wave background is produced, with present-day closure density $\Omega_{gw} \sim 10^{-5}$ -- $10^{-6}$ on $1$ - $10$ Mpc scales.

The non--linear dynamics of self--gravitating irrotational dust is analyzed in a general relativistic framework, using synchronous and comoving coordinates. Writing the equations in terms of the metric tensor of the spatial sections orthogonal to the fluid flow allows an unambiguous expansion in inverse powers of the speed of light. The Newtonian and post--Newtonian approximations are derived in Lagrangian form. A general formula for the gravitational waves generated by the non--linear evolution of cosmological perturbations is given. It is argued that a stochastic gravitational--wave background is produced by non--linear cosmic structures, with present--day closure density $\Omega_{gw} \sim 10^{-5}$ -- $10^{-6}$ on Mpc scale.