Ek is the translational kinetic energy at current time Ekr is the rotational kinetic energy at current time Ei is the internal energy Ewk is the work of external forces (energy brought to the system) E,1 is the energy at beginning of the RUN (not at time t=0) The Hourglass energy is not counted in this energy balance, so that a negative energy error generally occurs (except if using QEPH or BATOZ shells, and fully integrated solids, for which there is no Hourglass). It is bounded to 99%. The energy error is reset after each RESTART. It is possible to stop a job and rerun it by using a control file containing /CHKPT. In this case, error and energy values will restart from their last value. |
The Energy Error computed by RADIOSS is a percentage. If the error is negative, it means that some energy has been dissipated. In case of under integrated elements (Belytschko shells, solids with 1 integration point), the Hourglass energy can also explain a negative Energy Error since it is not counted in the energy balance. The normal amount of Hourglass energy is about -10% to -15%. If the error is positive, there is an energy creation. In case of using QEPH shell formulation or fully integrated elements, the Energy Error can be slightly positive since there is no Hourglass energy and the computation is much more accurate. An error of +1% or +2% is acceptable. If the positive energy error is greater than 2%, the source of this energy has to be identified. Incompatible kinematic conditions can lead to such a situation. An increasing Energy Error that reaches 99% can indicate the simulation has diverged. However, in certain situations high energy error can be acceptable:
|
The added mass can be due to Interface Type 2, Spotflag =1. In this case, the added mass is totally made at time t=0. It can also be due to options for constant time step (/DT/NODA/CST or /DT/INTER/CST). In case of added mass in the model, it is necessary to check if it is not too important with respect to the total mass of the model (see the DM/M value in the last column in the RADIOSS Engine listing file (Runname_nnnn.out)). It is also important to post-process this added mass in order to check that it is not too large locally, since this could mean false results (for checking this, the Animations written with /ANIM/NODA/DMAS have to be visualized). Also note that /DT/NODA/CST can lead to added inertia, which is much more difficult to appreciate from an engineering standpoint and can not be post-processed prior to RADIOSS V9.0. Added mass due to /DT/NODA/CST is often the cause of the divergence of computation. If interfaces are badly defined (too small gaps, initial penetrations, and so on), /DT/NODA/CST can not be the solution to increase time step. The model has to be improved. It is recommended not to use either /DT/NODA/CST or /DT/INTER/CST. It is possible to first make a short run with added mass by using /DT/NODA/CST (one cycle is sufficient) and then to check the added mass by post-processing the Animation written at time t=0. The computation can then be continued without using /DT/NODA/CST, in order to limit the possible problems it can generate while getting a better time step. |