Abstract: Vibration of pile was studied and a new theory was developed, in which surrounding soil was described as multi-components visco-elastic models. In detail, one soil model is that two Voigt bodies are connected in series (called as Model A), the other is that two Maxwell bodies are connected in parallel (called as Model B). Then, a series of analytical responses in the frequency or time domain were obtained by means of Laplace transform and invert Fourier transform. Based on this, several figures were prepared to illustrate influences of soil model parameters on the responses, and some important conclusions were drawn for pile integrity test. With Model A, it is shown that the vibration characteristics of pile are influenced mainly by one Voigt body with smaller parameters, and secondarily by the other. It is also shown that fluctuation of curves of amplitude versus frequency will decrease and reflection from pile toe in velocity response of pile top in time-domain will reduce too with damping coefficient of Voigt body being increased. On the other hand, with Model B, some different regulars are found that spring or damper in a Maxwell body won't act alone, that is to say, one have to depend on the other. Out of the ordinary, the influence of increasing of damper or spring coefficient would be equivalent to that of increasing of the other to some extent from traditional angle. Either Model A or Model B could be simplified into the traditional soil model of single Voigt body. So, the new theory is more perfect than the current one.