In this paper we prove the existence of small-amplitude quasi-periodic solutions with Sobolev regularity, for the d-dimensional forced Kirchhoff equation with periodic boundary conditions. This is the first result of this type for a quasi-linear equations in high dimension. The proof is based on a Nash-Moser scheme in Sobolev class and a regularization procedure combined with a multiscale analysis in order to solve the linearized problem at any approximate solution.The existence of periodic solutions for the forced Kirchhoff equation in any dimension has been proved by Baldi in [2], while the existence of quasi-periodic solutions in one space dimension under periodic boundary conditions has been proved in [43].Note that equation (1.5) is a quasi-linear PDE and it is well known that the existence of global solutions (even not periodic or quasi-periodic) for quasi-linear PDEs is not guaranteed, see for instance the nonexistence results in [36,39] The existence of periodic solutions for wave-type equations with unbounded nonlinearities has been proved for instance in [46,20,19]. For the water waves equations, which are fully nonlinear PDEs, we mention [32,33,34,1]; see also [3] for fully non-linear Benjamin-Ono equations.The methods developed in the above mentioned papers do not work for proving the existence of quasiperiodic solutions.The existence of quasi-periodic solutions for PDEs with unbounded nonlinearities has been developed by Kuksin [37] for KdV and then . This approach has been improved by 41] to deal with DNLS (Derivative Nonlinear Schrödinger) and Benjamin-Ono equations. These methods apply to dispersive PDEs like KdV, DNLS but not to derivative wave equation (DNLW) which contains first order derivatives in the nonlinearity. KAM theory for DNLW equation has been recently developed by Berti-Biasco-Procesi in [9,10]. Such results are obtained via a KAM-like scheme which is based on the so-called second Melnikov conditions and provides also the linear stability of the solutions.The existence of quasi-periodic solutions can be also proved by imposing only first order Melnikov conditions and the so-called multiscale approach. This method has been developed, for PDEs in higher space dimension, by Bourgain in [17,18,20] for analytic NLS and NLW, extending the result of Craig-Wayne [21] for 1-dimensional wave equation with bounded nonlinearity. Later, this approach has been improved by Berti-Bolle [12, 11] for NLW, NLS with differentiable nonlinearity and by Berti-Corsi-Procesi [14] on compact Lie-groups.This method is especially convenient in higher space dimension since the second order Melnikov conditions are violated, due to the high multiplicity of the eigenvalues. The drawback is that the linear stability is not guaranteed. Indeed there are very few results concerning the existence and linear stability of quasi-periodic solutions in the case of multiple eigenvalues. We mention [22,15] for the case of double eigenvalues and [25,26] in higher space dimension.All the aforementioned results concern semi-lin...