Explicit Dynamic Analysis

Explicit Dynamic Analysis in OptiStruct is provided through an integration of the Radioss Starter and Radioss Engine via a translator. Explicit integration schemes, are available via Radioss Integration.

Transparent to you, OptiStruct input data is directly translated into Radioss input data. The Radioss Starter and Radioss Engine are then executed, and the results are brought back into the OptiStruct output module to export the different output formats.

Solution Method

The basic concepts of the solution methods to highlight the characteristics of the solution methods and to identify the use of certain parameters to control convergence. The solution utilizes a general Newmark integration scheme.

The following equation of motion shall be solved.(1) M u ¨ + C u ˙ + K u = f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCytaiqahw hagaWaaiabgUcaRiaahoeaceWH1bGbaiaacqGHRaWkcaWHlbGaaCyD aiabg2da9iaahAgaaaa@3F32@
M MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCytaaaa@36CC@
Mass matrix
C MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCytaaaa@36CC@
Damping matrix
K MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCytaaaa@36CC@
Stiffness matrix
These matrices are derived using finite elements.
f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCytaaaa@36CC@
The vector describes the external loads
u MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCytaaaa@36CC@
Displacement vector
The dots describe the derivatives with respect to time.
The equation of motion can be solved using a general Newmark integration scheme. Newmark is a one-step time integration method. All solutions can be derived from it and are formulated in terms of a Time History (Figure 1).


Figure 1. Time History
In general Newmark, the state vector is computed:(2) u ˙ t + 1 = u ˙ t + Δ t [ ( 1 γ ) u ¨ t + γ u ¨ t + 1 ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabCyDayaaca WaaSbaaSqaaiaadshacqGHRaWkcaaIXaaabeaakiabg2da9iqahwha gaGaamaaBaaaleaacaWG0baabeaakiabgUcaRiabfs5aejaadshada WadaqaamaabmaabaGaaGymaiabgkHiTiabeo7aNbGaayjkaiaawMca aiqahwhagaWaamaaBaaaleaacaWG0baabeaakiabgUcaRiabeo7aNj qahwhagaWaamaaBaaaleaacaWG0bGaey4kaSIaaGymaaqabaaakiaa wUfacaGLDbaaaaa@4FA4@ (3) u t + 1 = u t + Δ t u ˙ t + ( 1 2 β ) Δ t 2 u ¨ t + β Δ t 2 u ¨ t + 1 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCyDamaaBa aaleaacaWG0bGaey4kaSIaaGymaaqabaGccqGH9aqpcaWH1bWaaSba aSqaaiaadshaaeqaaOGaey4kaSIaaeiLdiaadshaceWH1bGbaiaada WgaaWcbaGaamiDaaqabaGccqGHRaWkdaqadaqaamaalaaabaGaaGym aaqaaiaaikdaaaGaeyOeI0IaeqOSdigacaGLOaGaayzkaaGaaeiLdi aadshadaahaaWcbeqaaiaaikdaaaGcceWH1bGbamaadaWgaaWcbaGa amiDaaqabaGccqGHRaWkcqaHYoGycaqGuoGaamiDamaaCaaaleqaba GaaGOmaaaakiqahwhagaWaamaaBaaaleaacaWG0bGaey4kaSIaaGym aaqabaaaaa@572E@
Then the equation of motion yields:(4) [ M+γΔtC+βΔ t 2 K ] u ¨ t+1 =fC[ u ˙ t +( 1γ )Δt u ¨ t ]K[ u t +Δt u ˙ t +( 1 2 β )Δ t 2 u ¨ t ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaaca WHnbGaey4kaSIaeq4SdCMaaeiLdiaadshacaWHdbGaey4kaSIaeqOS diMaaeiLdiaadshadaahaaWcbeqaaiaaikdaaaGccaWHlbaacaGLBb GaayzxaaGabCyDayaadaWaaSbaaSqaaiaadshacqGHRaWkcaaIXaaa beaakiabg2da9iaahAgacqGHsislcaWHdbWaamWaaeaaceWH1bGbai aadaWgaaWcbaGaamiDaaqabaGccqGHRaWkdaqadaqaaiaaigdacqGH sislcqaHZoWzaiaawIcacaGLPaaacaqGuoGaamiDaiqahwhagaWaam aaBaaaleaacaWG0baabeaaaOGaay5waiaaw2faaiabgkHiTiaahUea daWadaqaaiaahwhadaWgaaWcbaGaamiDaaqabaGccqGHRaWkcaqGuo GaamiDaiqahwhagaGaamaaBaaaleaacaWG0baabeaakiabgUcaRmaa bmaabaWaaSaaaeaacaaIXaaabaGaaGOmaaaacqGHsislcqaHYoGyai aawIcacaGLPaaacaqGuoGaamiDamaaCaaaleqabaGaaGOmaaaakiqa hwhagaWaamaaBaaaleaacaWG0baabeaaaOGaay5waiaaw2faaaaa@70FA@
This can be rewritten into:(5) [ 1 βΔ t 2 M+ γ βΔt C+K ]Δ u t =Δ f ˜ t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada WcaaqaaiaaigdaaeaacqaHYoGycaqGuoGaamiDamaaCaaaleqabaGa aGOmaaaaaaGccaWHnbGaey4kaSYaaSaaaeaacqaHZoWzaeaacqaHYo GycaqGuoGaamiDaaaacaWHdbGaey4kaSIaaC4saaGaay5waiaaw2fa aiaabs5acaWH1bWaaSbaaSqaaiaadshaaeqaaOGaeyypa0JaaeiLdi qahAgagaacamaaBaaaleaacaWG0baabeaaaaa@4E8A@
using:(6) u t+1 = u t +Δ u t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCyDamaaBa aaleaacaWG0bGaey4kaSIaaGymaaqabaGccqGH9aqpcaWH1bWaaSba aSqaaiaadshaaeqaaOGaey4kaSIaaeiLdiaahwhadaWgaaWcbaGaam iDaaqabaaaaa@4113@
In short:(7) AΔu= f ˜ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCyqaiaabs 5acaWH1bGaeyypa0JabCOzayaaiaaaaa@3ADD@

The matrix A MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCyqaaaa@36C1@ is the dynamic stiffness. In nonlinear time-dependent problems, this system becomes nonlinear and its solution requires an additional iteration loop at each time step using a Newton-type method.

A conditionally stable explicit integration scheme can be derived from the Newmark scheme by setting:(8) γ = 1 2 , β = 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4SdCMaey ypa0ZaaSaaaeaacaaIXaaabaGaaGOmaaaacaGGSaGaeqOSdiMaeyyp a0JaaGimaaaa@3E3C@ (9) u ˙ t + 1 = u ˙ t + 1 2 Δ t [ u ¨ t + u ¨ t + 1 ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabCyDayaaca WaaSbaaSqaaiaadshacqGHRaWkcaaIXaaabeaakiabg2da9iqahwha gaGaamaaBaaaleaacaWG0baabeaakiabgUcaRmaalaaabaGaaGymaa qaaiaaikdaaaGaaeiLdiaadshadaWadaqaaiqahwhagaWaamaaBaaa leaacaWG0baabeaakiabgUcaRiqahwhagaWaamaaBaaaleaacaWG0b Gaey4kaSIaaGymaaqabaaakiaawUfacaGLDbaaaaa@4A61@ (10) u t+1 = u t +Δt u ˙ t + 1 2 Δ t 2 u ¨ t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCyDamaaBa aaleaacaWG0bGaey4kaSIaaGymaaqabaGccqGH9aqpcaWH1bWaaSba aSqaaiaadshaaeqaaOGaey4kaSIaaeiLdiaadshaceWH1bGbaiaada WgaaWcbaGaamiDaaqabaGccqGHRaWkdaWcaaqaaiaaigdaaeaacaaI Yaaaaiaabs5acaWG0bWaaWbaaSqabeaacaaIYaaaaOGabCyDayaada WaaSbaaSqaaiaadshaaeqaaaaa@49BB@
From these relationships the central differences explicit integration scheme can be derived.(11) u ˙ t + 1 / 2 = u ˙ t 1 / 2 + Δ t 12 u ¨ t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabCyDayaaca WaaSbaaSqaaiaadshacqGHRaWkcaaIXaGaai4laiaaikdaaeqaaOGa eyypa0JabCyDayaacaWaaSbaaSqaaiaadshacqGHsislcaaIXaGaai 4laiaaikdaaeqaaOGaey4kaSIaaeiLdiaadshadaWgaaWcbaGaaGym aiaaikdaaeqaaOGabCyDayaadaWaaSbaaSqaaiaadshaaeqaaaaa@485B@ (12) u t + 1 = u t + Δ t 2 u ˙ t + 1 / 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCyDamaaBa aaleaacaWG0bGaey4kaSIaaGymaaqabaGccqGH9aqpcaWG1bWaaSba aSqaaiaadshaaeqaaOGaey4kaSIaaeiLdiaadshadaWgaaWcbaGaaG OmaaqabaGcceWH1bGbaiaadaWgaaWcbaGaamiDaiabgUcaRiaaigda caGGVaGaaGOmaaqabaaaaa@460F@ (13) M u ¨ t+1 =fC u ˙ t+1/2 K u t+1 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCytaiqahw hagaWaamaaBaaaleaacaWG0bGaey4kaSIaaGymaaqabaGccqGH9aqp caWHMbGaeyOeI0IaaC4qaiqahwhagaGaamaaBaaaleaacaWG0bGaey 4kaSIaaGymaiaac+cacaaIYaaabeaakiabgkHiTiaahUeacaWH1bWa aSbaaSqaaiaadshacqGHRaWkcaaIXaaabeaaaaa@4911@
Figure 2 illustrates the relationships.


Figure 2. Explicit Integration
Assuming that(14) C u ˙ t+1 C u ˙ t+1/2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaC4qaiqahw hagaGaamaaBaaaleaacaWG0bGaey4kaSIaaGymaaqabaGccqGHijYU caWHdbGabCyDayaacaWaaSbaaSqaaiaadshacqGHRaWkcaaIXaGaai 4laiaaikdaaeqaaaaa@424B@
The equation of motion for the central differences scheme simplifies to:(15) M u ¨ t + 1 = f C u ˙ t + 1 K u t + 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCytaiqahw hagaWaamaaBaaaleaacaWG0bGaey4kaSIaaGymaaqabaGccqGH9aqp caWHMbGaeyOeI0IaaC4qaiqahwhagaGaamaaBaaaleaacaWG0bGaey 4kaSIaaGymaaqabaGccqGHsislcaWHlbGaaCyDamaaBaaaleaacaWG 0bGaey4kaSIaaGymaaqabaaaaa@47A3@ (16) M u ¨ t+1 =f B t+1 t σ t+1 dV MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCytaiqahw hagaWaamaaBaaaleaacaWG0bGaey4kaSIaaGymaaqabaGccqGH9aqp caWHMbGaeyOeI0Yaa8qaaeaacaWHcbWaa0baaSqaaiaadshacqGHRa WkcaaIXaaabaGaamiDaaaaaeqabeqdcqGHRiI8aOGaeq4Wdm3aaSba aSqaaiaadshacqGHRaWkcaaIXaaabeaakiaadsgacaWGwbaaaa@4A57@
These central differences scheme is used for explicit analysis. The time step must always be smaller than the critical time step to ensure stability of the solution. The critical time step depends on the highest frequency in the system and is computed from the corresponding angular frequency ω max as:(17) Δ t cr = 2 ω max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeiLdiaads hadaWgaaWcbaGaam4yaiaadkhaaeqaaOGaeyypa0ZaaSaaaeaacaaI YaaabaGaeqyYdC3aaSbaaSqaaiGac2gacaGGHbGaaiiEaaqabaaaaa aa@40BE@
For a discrete system, the time step must be small enough to excite all frequencies in the finite element mesh. This requires such a short time step that a shock wave does not miss any node when traveling the mesh. Therefore, (18) Δ t l c c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeiLdiaads hacqGHKjYOdaWcaaqaaiaadYgadaWgaaWcbaGaam4yaaqabaaakeaa caWGJbaaaaaa@3CC5@

With l c being the critical length of an element and c is the speed of sound in the given material.

Different ways of time step control are available. The default method is the nodal time step which is computed from the nodal mass m and the equivalent nodal stiffness k such that:(19) Δ t cm = min n 2m k MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeiLdiaads hadaWgaaWcbaGaam4yaiaad2gaaeqaaOGaeyypa0ZaaCbeaeaaciGG TbGaaiyAaiaac6gaaSqaaiaad6gaaeqaaOWaaOaaaeaadaWcaaqaai aaikdacaWGTbaabaGaam4AaaaaaSqabaaaaa@41F1@

The element time step based on the critical length of each element is also available. The choice can be made on the XSTEP Bulk Data Entry.

Problem Setup

Explicit Dynamic Analysis is defined through a SUBCASE.

An XSTEP statement as well as ANALYSIS=EXPDYN must be present in the subcase. To define the termination time a TTERM Subcase Information Entry is mandatory. XSTEP references an XSTEP Bulk Data Entry. Time step control can be defined on the XSTEP Bulk Data Entry.

The definition of a unit system through the DTI, UNITS or UNITS Bulk Data Entry statement is required.

The Explicit Dynamic Analysis loads and boundary conditions are defined in the Bulk Data Entry section of the input deck. They need to be referenced under the SUBCASE using an SPC, NLOAD, LOAD, IC and RWALL statements in a SUBCASE. Each SUBCASE defines one loading condition that is executed separately.

Subcase continuation is available through the use of CNTNLSUB. Any number of explicit and implicit analyses can be linked. However, explicit (ANALYSIS=EXPDYN) analysis subcases cannot yet be linked with small displacement quasi-static nonlinear (ANALYSIS=NLSTAT) analysis subcases and vice versa.

Example: Explicit Analysis

SUBCASE       3
  ANALYSIS = EXPDYN
  SPC = 1
  NLOAD = 2
  XSTEP = 3
  TTERM = 1.0
  DISP = ALL
  STRESS = ALL
BEGIN BULK
XSTEP,3
NLOAD1,2,2,,L,88
TABLED1,88,
+,0.0,0.0,1.0,1.0,ENDT
DTI,UNITS,1,kg,N,m,s

Example: Subcase Continuation

DISP = ALL
  STRESS = ALL
SUBCASE       1
  ANALYSIS = EXPDYN
  IC = 5
  XSTEP = 4
  TTERM = 1.0
SUBCASE       2
  ANALYSIS = EXPDYN
  IC = 5
  XSTEP = 4
  TTERM = 1.1
  CNTNLSUB = 1

User's Considerations

Explicit Dynamic Analysis Properties and Materials

Special element types and nonlinear materials are available for Explicit Dynamic Analysis. As a general rule, property and material definitions that are only applicable in Explicit Dynamic Analysis are defined on extensions to the original property and to a MAT1 material, respectively. The extensions are grouped with the base entry by sharing the same PID or MID. In the case of a subcase that is not an Explicit Dynamic Analysis, these extensions are ignored. Property defaults can be set for shells (XSHLPRM) and solids (XSOLPRM) that may replace the use of property extensions.

Property Example:
PSHELL, 3, 7, 1.0, 7, , 7
PSHELLX, 3, 24, , , 5
Material Example:
MAT1, 102, 60.4, , 0.33, 2.70e-6
MATX02, 102, 0.09026, 0.22313, 0.3746, 100.0, 0.175

Coordinate Systems

In Explicit Dynamic Analysis there are moving and fixed coordinate systems. Rectangular coordinate systems that are defined through grid points (CORD1R, CORD3R) are moving with the deformations of the model. Systems defined in terms of point coordinates (CORD2R, CORD4R) are fixed.

The behavior of loads depends on the coordinate system referenced. If the loads FORCE and MOMENT are desired to be follower forces, a CID that references a moving coordinate system (CORD1R, CORD3R) must be defined. Otherwise these loads are not following the deformation. PLOAD always follows the deformations.

Characteristics

Explicit Dynamic Analysis has the following characteristics:
  • In general, a diagonal mass matrix is used
  • No matrix factorization necessary
  • Equilibrium is always guaranteed
  • Maximum stable time step needs to be respected
  • Small time steps
  • Short-term events

Limitations

The solution will be terminated if unsupported Bulk Data Entries are encountered.
  • The following Bulk Data properties and elements are currently not translated:
    • PBUSHT (partially, KN is translated)
    • PDAMP, CDAMPi
    • PGAP, CGAP, CGAPG (partially, friction is not allowed)
    • PMASS, CMASSi
    • PSHEAR, CSHEAR
    • PVISC, CVISC
  • Additional relevant Bulk Data Entries (except loads) that are currently not translated:
    • CORD1C, CORD1S, CORD2S
    • DMIG
    • MAT2, MAT4, MAT5, MAT8, MAT9, MAT10
    • MATTi, TABLEST
    • MPC, MPCADD
    • RBE1, RROD
  • Relevant loads that are currently not translated:
    • PLOAD1, PLOAD2
    • PLOAD4 (partially, N1, N2, N3 cannot be used)
    • RFORCE (partially, RACC is not supported)
    • TLOAD1, TLOAD2