In order to study a class of quasilinear p-biharmonic equations with Hardy terms and multi-critical Sobolev-Hardy exponents, the existence theorem of the solutions to the above problem is established by means of the Ekeland variational principle. Firstly, to guarantee the variational functional is bounded from below, it is restricted on a set M_η (usually called Nehari manifold). Secondly, the set M_η is divided into three parts M⁺_η, M⁰_η and M^-_η by using fibering maps. Moreover, the existence of minimum in M⁺_η and M^-_η is proved by studying the properties of the two subsets. Finally, by using implicit function theorem, it is found that there exists a minimizing sequence {u_n} making the (PS)_c conditions hold when the parameters satisfy certain conditions. Therefore, the existence of the solutions to the problem is proved. The conclusions provide a theoretical basis for judging the structure and properties of the solutions.