Spatial mesh adaptation in semidiscrete finite element analysis of linear elastodynamic problems

A spatial mesh adaptation procedure in semidiscrete finite element analysis of 2D linear elastodynamic problems is presented. The procedure updates, through an automatic remeshing scheme, the spatial mesh when found necessary in order to gain control of the spatial discretization error from time to time. An a posteriori error estimate developed by Zienkiewicz and Zhu (1987) for elliptic problems is extended to dynamic analysis to estimate the spatial discretization error at a certain time, which is found to be reasonable by analyzing an a priori error estimate. Numerical examples are used to demonstrate the performance of the procedure. It is indicated that the extended error estimation and the procedure are capable of monitoring the moving of steep stress regions by updating the spatial mesh according to a prescribed error tolerance, thus providing a reliable finite element solution in an efficient manner.