Abstract: A new direct integration method based on Gurtin variational principle of displaement model is presented. The interpolation function with non-time-step-parameter is introduced to approximate the bnodal displacement in time domain An unconditional sfable formula is given Examples indicate that the algorithm is simple and accurate. It can be epplied to many enginering种areas.