The Finite Element Method with High-Order Enrichment Functions for Elastodynamic Analysis

Elastodynamic problems are investigated in this work by employing the enriched finite element method (EFEM) with various enrichment functions. By performing the dispersion analysis, it is confirmed that for elastodynamic analysis, the amount of numerical dispersion, which is closely related to the numerical error from the space domain discretization, can be suppressed to a very low level when quadric polynomial bases are employed to construct the local enrichment functions, while the amount of numerical dispersion from the EFEM with other types of enrichment functions (linear polynomial bases or first order of trigonometric functions) is relatively large. Consequently, the present EFEM with a quadric polynomial enrichment function shows more powerful capacities in elastodynamic analysis than the other considered numerical techniques. More importantly, the attractive monotonic convergence property can be broadly realized by the present approach with the typical two-step Bathe temporal discretization technique. Three representative numerical experiments are conducted in this work to verify the abilities of the present approach in elastodynamic analysis.