by L. Idoux, EDF R&D / AMA
The energy balance can be a very useful information. First, it provides information on the quality of results in the same way that a residue of equilibrium in the Newton algorithm, and secondly it helps to identify the physical phenomena involved and so meet the following questions: what is the energy brought to my structure and what becomes of it?
Up to now, the calculation of energy is available in Code_Aster, but not for every case. For example, it is not possible to know the energy dissipated by damping. In addition, the access to the energy balance is only possible afterwards by post-processing, which raises some difficulties.
With this in mind, a new calculation option was developed. It calculates the energy balance of the full model during the calculation. This option is activated by the keyword ENERGIE in the commands STAT_NON_LINE U4.51.03, DYNA_NON_LINE U4.53.01 and DYNA_VIBRA U4.53.02 ; for the latter, in the case of calculations on physical basis. The presence of this keyword triggers the calculation of the energy balance as follows:
The results are displayed in the message file and saved in a table at the end of the calculation for recovery. For theoretical precisions, please see the new documentation R4.09.01 (automatic translation).
To illustrate this development, the calculation of the energy balance is enabled in the test case sdlv120a (automatic translation). It consists of the absorption of a compression wave in an elastic bar, with absorbing boundary elements. The blue curve represents the energy supplied from the outside, as an imposed displacement at the end of the bar. This energy propagates along the bar in the form of kinetic energy and total strain energy, and is fully dissipated at the other end of the bar through the absorbing boundary elements.
Available in the 11.1.7 version and future releases