Computation Checks

This is a method to check the quality of simulation results. Even if the run does not fail, it needs to be verified that three fundamental conservation laws are respected. The time step variation and the qualitative evaluation of results may help to avoid modeling errors.

Energy Balance

Taking into account the external works, the total energy must remain constant or decrease slightly. The total energy can increase at the end of the computation, during the spring back or at the beginning during the first cycles.

Internal energy + Kinematic energy + Hourglass energy + Contact Energy + … = Variation of the External Work

If under-integrated elements are used, the total hourglass energy must remain lower than 10% of the total energy. If this is not the case, the mesh should be reworked or elements with physical stabilization method should be used.

The contact energy is not really physical. For each subset and for each part the following limitation is recommended:(1)
E h + E c E t 15 % MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca WGfbWaaSbaaSqaaiaadIgaaeqaaOGaey4kaSIaamyramaaBaaaleaa caWGJbaabeaaaOqaaiaadweadaWgaaWcbaGaamiDaaqabaaaaOGaey izImQaaGymaiaaiwdacaGGLaaaaa@408F@
Where,
E h MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaBa aaleaacaWGObaabeaaaaa@37DA@
Hourglass energy
E c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaBa aaleaacaWGJbaabeaaaaa@37D5@
Contact energy
E t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaBa aaleaacaWG0baabeaaaaa@37E6@
Total energy
The energy error is computed by Radioss as:(2)
%Error=100( E k + E k r + E i E k,1 + E k,1 r + E i,1 + E wk E wk,1 1 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiyjaiaadw eacaWGYbGaamOCaiaad+gacaWGYbGaeyypa0JaaGymaiaaicdacaaI WaWaaeWaaeaadaWcaaqaaiaadweadaWgaaWcbaGaam4AaaqabaGccq GHRaWkcaWGfbWaa0baaSqaaiaadUgaaeaacaWGYbaaaOGaey4kaSIa amyramaaBaaaleaacaWGPbaabeaaaOqaaiaadweadaWgaaWcbaGaam 4AaiaacYcacaaIXaaabeaakiabgUcaRiaadweadaqhaaWcbaWaaSba aWqaaiaadUgacaGGSaGaaGymaaqabaaaleaacaWGYbaaaOGaey4kaS IaamyramaaBaaaleaacaWGPbGaaiilaiaaigdaaeqaaOGaey4kaSIa amyramaaBaaaleaacaWG3bGaam4AaaqabaGccqGHsislcaWGfbWaaS baaSqaaiaadEhacaWGRbGaaiilaiaaigdaaeqaaaaakiabgkHiTiaa igdaaiaawIcacaGLPaaaaaa@6058@
Where,
E k
Translational kinetic energy
E k r
Rotational kinetic energy
E i
Internal energy
E w k
External work
E k,1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaBa aaleaacaWGRbGaaiilaiaaigdaaeqaaaaa@3948@
Initial translational kinetic energy
E k,1 r MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaDa aaleaadaWgaaadbaGaam4AaiaacYcacaaIXaaabeaaaSqaaiaadkha aaaaaa@3A78@
Initial rotational kinetic energy
E i,1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaBa aaleaacaWGPbGaaiilaiaaigdaaeqaaaaa@3946@
Initial internal energy
E wk,1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyramaaBa aaleaacaWG3bGaam4AaiaacYcacaaIXaaabeaaaaa@3A44@
Initial external work

The error must be negative and decreasing (except for the first cycle or in spring-back stage). The error must be less than 15% at the end of computation (5% for a good model). If the error increases a little, then decreases may be normal. If the error increases, this means that a problem has occurred. Then, the error may grow to 99.9% with computation failure message.

Mass Balance

If the mass increases, its variation must remain smaller than 1% for each subset and for each part ( d M M < 1 % MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca WGKbGaamytaaqaaiaad2eaaaGaeyipaWJaaGymaiaacwcaaaa@3AFC@ ). If the mass variation is between 1% to 3%, check if the nodes with the added mass are moving or not. If this is the case, the added mass results in an increase in kinetic energy. For more than 3% of variation, the results are probably bad.

Momentum Balance

The dynamic equilibrium of each node is satisfied by the Newton law at the end of each cycle. As Radioss resolves the equilibrium equations at each cycle, normally the momentum balance is satisfied. However in case of a problem, a cross-check between nodal accelerations and the impactor forces (interface, rigid wall, barrier, etc.) can help to better understand problems.(3)
F wall = Monentum variation dt MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOramaaBa aaleaacaWG3bGaamyyaiaadYgacaWGSbaabeaakiabg2da9maalaaa baGaamytaiaad+gacaWGUbGaamyzaiaad6gacaWG0bGaamyDaiaad2 gafaqabeqabaaabaaaaiaadAhacaWGHbGaamOCaiaadMgacaWGHbGa amiDaiaadMgacaWGVbGaamOBaaqaaiaadsgacaWG0baaaaaa@4DB4@

Time Step Evolution

If the time step decreases and then it goes up quickly, there is not a problem. If it varies greatly from one cycle to another, it may be due to the interface stiffness.

If the time step remains low, a problem has occurred. In this case, find the node (or element) controlling the time step and try to understand why the decrease occurred.

Visual Inspection of Deformed Parts

After getting some animation files, the overall deformation of the structure can be compared to the physical behavior. Inspect the deformed shape to verify:
  • There are no flying nodes (parts)
  • The deformed shape is smooth
  • The chord angle between adjacent elements is sufficiently small
  • There are no intersections
Good physical behavior is obtained when:
  • The plastic strain is less than 30%(4)
    E k + E k r E t 1 % MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaaca WGfbWaaSbaaSqaaiaadUgaaeqaaOGaey4kaSIaamyramaaDaaaleaa caWGRbaabaGaamOCaaaaaOqaaiaadweadaWgaaWcbaGaamiDaaqaba aaaOGaeyizImQaaGymaiaacwcaaaa@40D3@

Meaningless Results

Contact forces, von Mises stress, nodal velocities, and accelerations have to be checked carefully. If the values are meaningless (for example: von Mises stress = 1 GPa), the first check may be the unit system consistency. Refer to the Appendix for more information.