Two Positive Solutions for Kirchhoff Type Problems with Hardy-Sobolev Critical Exponent and Singular Nonlinearities

We consider the following singular Kirchhoff type equation with Hardy-Sobolev critical exponent \[ \begin{cases} \displaystyle -\left( a + b \int_{\Omega} |\nabla u|^2 \, dx \right) \Delta u = \frac{u^{3}}{|x|} + \frac{\lambda}{|x|^{\beta} u^{\gamma}}, & x \in \Omega, \\ u > 0, & x \in \Omega, \\ u = 0, & x \in \partial \Omega, \end{cases} \] where $\Omega \subset \mathbb{R}^{3}$ is a bounded domain with smooth boundary $\partial \Omega$, $0 \in \Omega$, $a,b,\lambda \gt 0$, $0 \lt \gamma \lt 1$, and $0 \leq \beta \lt (5+\gamma)/2$. Combining with the variational method and perturbation method, two positive solutions of the equation are obtained.