We prove a stable blowup result for solutions to the Einstein-scalar field and Einstein-stiff fluid systems. Our result applies to small perturbations of the spatially flat Friedmann–Lemaitre–Robertson–Walker (FLRW) solution with topology $$(0,\infty ) \times \mathbb {T}^3$$ . The FLRW solution models a spatially uniform scalar-field/stiff fluid evolving in a spacetime that expands towards the future and that has a “Big Bang” singularity at $$\lbrace 0 \rbrace \times \mathbb {T}^3$$ , where various curvature invariants blow up. We place “initial” data on a Cauchy hypersurface $$\Sigma _1'$$ that are close, as measured by a Sobolev norm, to the FLRW data induced on $$\lbrace 1 \rbrace \times \mathbb {T}^3$$ . We then study the asymptotic behavior of the perturbed solution in the collapsing direction and prove that its basic qualitative and quantitative features closely resemble those of the FLRW solution. In particular, for the perturbed solution, we construct constant mean curvature-transported spatial coordinates covering $$(t,x) \in (0,1] \times \mathbb {T}^3$$ and show that it also has a Big Bang at $$\lbrace 0 \rbrace \times \mathbb {T}^3$$ , where its curvature blows up. The blowup confirms Penrose’s Strong Cosmic Censorship hypothesis for the “past-half” of near-FLRW solutions. Furthermore, we show that the equations are dominated by kinetic (that is, time-derivative-involving) terms that induce approximately monotonic behavior near the Big Bang. As a consequence of the monotonicity, we also show that various time-rescaled components of the solution converge to regular functions of x as $$t \downarrow 0$$ . The most difficult aspect of the proof is showing that the solution exists for $$(t,x) \in (0,1] \times \mathbb {T}^3$$ , and to this end, we derive a hierarchy of energy estimates that allow for the possibility of mild energy blowup as $$t \downarrow 0$$ . To close these estimates, it is essential that we are able to rule out more singular energy blowup, a step that is in turn tied to the most important ingredient in our analysis: a collection of integral identities that, when combined in the right proportions, yield an $$L^2$$ -type approximate monotonicity inequality, a key point being that the error terms are controllable up to the singularity for near-FLRW solutions. In a companion article, we derived similar approximate monotonicity inequalities for linearized versions of the Einstein-scalar field equations and used them to prove linear stability results for a family of spatially homogeneous background solutions. The present article shows that the linear stability of the FLRW background solution can be upgraded to a full proof of the nonlinear stability of its singularity.