Abstract: | Since the attenulation of propagating waves through soil/rock is related to the localized material properties as well as the strain developed, the commonly used Rayleigh-type damping model and its variations are not suitable for dynamic finite element analysis of such materials. A linear viscoelastic material model based on the concept of the relaxation spectrum is manipualted in place of the damping model in this paper. The method proposed by Day and Minster11 to transform the convolutional form of the stress–strain relationship to a set of differential operators using the Pade approximant method is generalized to non-scalar waves and implemented for transient finite element analyses. A time-marching scheme is also proposed to incorporate the resultant differential operators into the governing equation of motion. The accuracy related to the Pade approximant method and the time-marching scheme is investigated by critically analysing some scalar wave propagation problems. The proposed technique is further verified using two one-dimensional stress wave propagation problems and a two-dimensional transient propagating wave through an unbounded linear viscoelastic medium. Some encouraging results have been obtained using the proposed technique and guidelines for using this technique are also presented. Comparisons of analytical solutions obtained by Fourier synthesis and numerical results have been provided. |