In this paper we study the convergence of nonlinear Dirichlet problems for systems of variational elliptic PDEs defined on randomly perforated domains of $$\mathbb {R}^n$$
R
n
. Under the assumption that the perforations are small balls whose centres and radii are generated by a stationary short-range marked point process, we obtain in the critical-scaling limit an averaged nonlinear analogue of the extra term obtained in the classical work of Cioranescu and Murat (Res Notes Math III, 1982). In analogy to the random setting recently introduced by Giunti, Höfer and Velázquez (Commun Part Differ Equ 43(9):1377–1412, 2018) to study the Poisson equation, we only require that the random radii have finite $$(n-q)$$
(
n
-
q
)
-moment, where $$1<q<n$$
1
<
q
<
n
is the growth-exponent of the associated energy functionals. This assumption on the one hand ensures that the expectation of the nonlinear q-capacity of the spherical holes is finite, and hence that the limit problem is well defined. On the other hand, it does not exclude the presence of balls with large radii, that can cluster up. We show however that the critical rescaling of the perforations is sufficient to ensure that no percolating-like structures appear in the limit.