Guided wave propagation

For this study, the shape of the propagation domain is chosen as an elongated rectangle with a = 1 m and b = 0.05 m.

The initial conditions are of null displacements and null vertical velocity, while horizontal velocity is depicted below (solid_ic_bc.hpp file):

Material properties are: 7800 kg/m³ of density, 20.9e10 Young modulus, and 0.3333 Poisson coefficient, as specified in props.txt file.

The generalized alpha method is employed with rho inf of 0.75 with zero acceleration, the time step is 1e-8 s, and the grid has 2e-3 m of spacing around the initial velocity patch, and 5e-3 m along the boundaries of the domain. For details, see the setup.txt and mesh.geo files.

The evolution of horizontal and vertical displacements computed with Gales, in comparison with the reference solution follows:

The video shows horizontal and vertical displacements over the space-time domain:

Reference

Chiappa et al. “An analytical benchmark for a 2D problem of elastic wave propagation in a solid” (2021) Engineering Structures