A boundary element methodology for viscoelastic analysis: Part I with cells

In this study Kelvin and Boltzmann viscoelastic models are implemented in a two-dimensional boundary element atmosphere. This general methodology is based on differential constitutive relations for viscoelasticity, avoiding the use of relaxation functions. Part I describes a methodology using internal cells. This methodology makes it possible to consider viscous parameters, which are not proportional to elastic tensor. From the kinematical relation between material and strain velocities at the approximation level a simple time marching process is achieved. At the end of Part I, numerical examples are provided to validate the methodology. The BEM viscoelastic formulation without using cells is carefully described in Part II.