Bulk wave propagation

The bulk wave propagation sees the propagation of an elastodynamic wave in the bulk of a 1 x 1 square, with null initial displacements and initial velocity as depicted below (see solid_ic_bc.hpp file):

Simulation is performed under plane strain assumption, with material properties of 7800 kg/m³ of density, 20.9e10 Young modulus and 0.3333 Poisson coefficient, as specified in props.txt file.

The domain is discretised into triangular elements of size 1e-3 m in the interior and 5e-3 m at the boundaries (mesh.geo file).

The generalized alpha-method of solution is employed with 0.75 of rho inf and zero acceleration. The time step id 1e-8 s. Setup can be found in setup.txt file.

Results of horizontal and vertical displacements at x = 0.25 m and y = 0.25 m are compared with the reference solution below:

A comparison the displacements at 1.7e-5 s over the domain follows:

The time evolution of displacements is visualized in the following movies.

Reference

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