Abstract
Viscoelastic solids may be effectively treated by the boundary element method (BEM) in the Laplace domain. However, calculation of transient response via the Laplace domain requires the inverse transform. Since all numerical inversion formulas depend heavily on a proper choice of their parameters, a direct evaluation in the time domain seems to be preferable. On the other hand, direct calculation of viscoelastic solids in the time domain requires the knowledge of viscoelastic fundamental solutions. Such solutions are simply obtained in the Laplace domain with the elastic–viscoelastic correspondence principle, but not in the time domain. Due to this, a quadrature rule for convolution integrals, the ‘convolution quadrature method' proposed by Lubich, is applied. This numerical quadrature formula determines their integration weights from the Laplace transformed fundamental solution and a linear multistep method. Finally, a boundary element formulation in the time domain using all the advantages of the Laplace domain formulation is obtained. Even materials with complex Poisson ratio, leading to time-dependent integral free terms in the boundary integral equation, can be treated by this formulation. Two numerical examples, a 3D rod and an elastic concrete slab resting on a viscoelastic half-space, are presented in order to assess the accuracy and the parameter choice of the proposed method. Copyright © 1999 John Wiley & Sons, Ltd.
Highlights
The boundary element method (BEM) has become a widely used numerical tool in statics and dynamics
A review about fractional calculus applied to dynamic problems has recently been given by Rossikhin and Shitikova.[6]
Viscoelastic boundary element formulations are mostly published for the quasi-static case, or in dynamics using a frequency or Laplace domain representation of the governing integral equation
Summary
Viscoelastic solids may be eectively treated by the boundary element method (BEM) in the Laplace domain. Direct calculation of viscoelastic solids in the time domain requires the knowledge of viscoelastic fundamental solutions. Such solutions are obtained in the Laplace domain with the elastic viscoelastic correspondence principle, but not in the time domain. A quadrature rule for convolution integrals, theconvolution quadrature method' proposed by Lubich, is applied. This numerical quadrature formula determines their integration weights from the Laplace transformed fundamental solution and a linear multistep method. KEY WORDS viscoelasticity; boundary integral equations; time domain; transform methods
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have