In this paper we discuss the following Kirchhoff equation
\left\{ \begin{array}{lr} -\left(a+b \int_{\mathbb{R}^3}|\nabla u|^{2} d x\right) \Delta u+V(x)u+\lambda u=\mu|u|^{q-2}u+|u|^{p-2}u \ {\rm in}\ \mathbb{R}^3,&\\ \int_{\mathbb{R}^{3}}u^{2}dx=c^2, \end{array} \right. where a, b, µ and c are positive numbers, λ is unknown and appears as a Lagrange multiplier,
We consider multilinear commutators of singular integrals defined by $T_{\vec{b}}f(x) =\int_{\mathbb{R}^n}\prod^m_{i=1}(b_i(x)-b_i(y))K(x, y)f(y)dy,$
where K is a standard Calder\'{o}n-Zygmund kernel, m is a positive integer and \vec{b} b =(b1,b2,...,bm) is a family of m locally integrable functions. Based on the theory of variable exponent and on generalization of the BMO norm, we prove the boundedness of multilinear commutators T_{\vec{b}} on grand variable Herz spaces. The result is still new even in the special case of m=1.