ATHENA 3D: A finite element code for ultrasonic wave propagation

12th Anglo-French Physical Acoustics Conference

Rose, C.^1 ; Rupin, F.^2 ; Fouquet, T.^1 ; Chassignole, B.^2

EDF RandD SINETICS, Clamart (92), France^1

EDF-RandD MMC, Sites des Renardières, Moret sur Loing (77), France^2

关键词: 3D finite element model;    Anisotropic material;    Complex geometries;    Computing processors;    Elastodynamic equation;    Fictitious domains;    Finite element codes;    Particle velocities;   

Others  :  https://iopscience.iop.org/article/10.1088/1742-6596/498/1/012009/pdf
DOI  :  10.1088/1742-6596/498/1/012009

来源: IOP

【 摘 要 】

The understanding of wave propagation phenomena requires use of robust numerical models. 3D finite element (FE) models are generally prohibitively time consuming. However, advances in computing processor speed and memory allow them to be more and more competitive. In this context, EDF R&D developed the 3D version of the well-validated FE code ATHENA2D. The code is dedicated to the simulation of wave propagation in all kinds of elastic media and in particular, heterogeneous and anisotropic materials like welds. It is based on solving elastodynamic equations in the calculation zone expressed in terms of stress and particle velocities. The particularity of the code relies on the fact that the discretization of the calculation domain uses a Cartesian regular 3D mesh while the defect of complex geometry can be described using a separate (2D) mesh using the fictitious domains method. This allows combining the rapidity of regular meshes computation with the capability of modelling arbitrary shaped defects. Furthermore, the calculation domain is discretized with a quasi-explicit time evolution scheme. Thereby only local linear systems of small size have to be solved. The final step to reduce the computation time relies on the fact that ATHENA3D has been parallelized and adapted to the use of HPC resources. In this paper, the validation of the 3D FE model is discussed. A cross-validation of ATHENA 3D and CIVA is proposed for several inspection configurations. The performances in terms of calculation time are also presented in the cases of both local computer and computation cluster use.