Nonlinear Transient Analysis

Nonlinear Small Displacement Transient and Nonlinear Large Displacement Transient (LGDISP) Analysis are currently available to solve nonlinear problems which include transient effects.

The predominant difference between Nonlinear Static Analysis and Nonlinear Transient Analysis is the inclusion of inertia and momentum terms in the solution of the equation of motion. This nonlinear transient solution sequence typically supports all nonlinear features supported by Nonlinear Static Analysis, including geometric nonlinearity, material nonlinearity, and contact. Subcase continuation between Nonlinear Transient Analysis and Nonlinear Static Analysis is also supported.

Nonlinearity Sources

Geometric Nonlinearity

In analyses involving geometric nonlinearity, changes in geometry as the structure deforms are considered in formulating the constitutive and equilibrium equations. Many engineering applications require the use of large deformation analysis based on geometric nonlinearity. For example, applications such as metal forming, tire analysis, and medical device analysis. Small deformation analysis based on geometric nonlinearity is required for some applications, like analysis involving cables, arches and shells. Such applications involve small deformation, except finite displacement or rotation.

Material Nonlinearity

Material nonlinearity involves the nonlinear behavior of a material based on current deformation, deformation history, rate of deformation, temperature, pressure, and so on.

Constraint and Contact Nonlinearity

Constraint nonlinearity in a system can occur if kinematic constraints are present in the model. The kinematic degrees-of-freedom of a model can be constrained by imposing restrictions on its movement. In OptiStruct, constraints are enforced with Lagrange multipliers.

In the case of contact, the constraint condition is based on inequalities and such a constraint generally does not allow penetration between any two bodies in contact.

Follower Load

Applied loads can depend upon the deformation of the structure when large deformations are involved. Geometrically, the applied loads (Forces or Pressure) can deviate from their initial direction based on how the model deforms at the location of application of load. In OptiStruct, if the applied load is treated as follower load, the orientation and/or the integrated magnitude of the load will be updated with changing geometry throughout the analysis.

Nonlinear Transient Solution Method

Nonlinear Transient solutions are done by marching through the solution in the time domain. The time-discretized equation is solved using Newton's method (similar to NLSTAT). Two time stepping schemes are provided for nonlinear transient analysis. The TMTD field on the TSTEP entry can be used to switch between the two integration schemes. The Generalized Alpha method (TMTD=1, default) and the Backward Euler method (TMTD=2).

Generalized Alpha (α) Method

In the Generalized Alpha method, the equilibrium equation takes the following form:

f ext t+αh =M 1 α m a t+αh + α m a t +C v t+αh + f int t+αh MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWHMbWaa0baaSqaaiaadwgacaWG4b GaamiDaaqaaiaadshacqGHRaWkcqaHXoqycaWGObaaaOGaeyypa0Ja aCytamaabmaabaWaaeWaaeaacaaIXaGaeyOeI0IaeqySde2aaSbaaS qaaiaad2gaaeqaaaGccaGLOaGaayzkaaGaaCyyamaaCaaaleqabaGa amiDaiabgUcaRiabeg7aHjaadIgaaaGccqGHRaWkcqaHXoqydaWgaa WcbaGaamyBaaqabaGccaWHHbWaaWbaaSqabeaacaWG0baaaaGccaGL OaGaayzkaaGaey4kaSIaaC4qaiaahAhadaahaaWcbeqaaiaadshacq GHRaWkcqaHXoqycaWGObaaaOGaey4kaSIaaCOzamaaDaaaleaacaWG PbGaamOBaiaadshaaeaacaWG0bGaey4kaSIaeqySdeMaamiAaaaaaa a@5F57@

Where,
M
Mass
C
Viscous damping
K
Stiffness matrix
f
Total force

Subscripts, e x t and i n t indicating external and internal force, respectively. Superscripts, t and t + h indicate the time when the variable is calculated, with t being the previous time and t + h being the current time when the displacement increment is being solved for. The Lagrangian coordinate is denoted by x , and the material time derivatives of x are denoted by v (velocity) and a (acceleration).

The velocity and acceleration at current time are computed with Newmark approximation by involving parameters γ and β .

v t + h = v t + h ( ( 1 γ ) a t + γ a t + h ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWH2bWaaWbaaSqabeaacaWG0bGaey 4kaSIaamiAaaaakiabg2da9iaahAhadaahaaWcbeqaaiaadshaaaGc cqGHRaWkcaWGObWaaeWaaeaadaqadaqaaiaaigdacqGHsislcqaHZo WzaiaawIcacaGLPaaacaWHHbWaaWbaaSqabeaacaWG0baaaOGaey4k aSIaeq4SdCMaaCyyamaaCaaaleqabaGaamiDaiabgUcaRiaadIgaaa aakiaawIcacaGLPaaaaaa@49AA@
u t + h = u t + h v t + 1 2 h 2 ( ( 1 2 β ) a t + 2 β a t + h ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWH1bWaaWbaaSqabeaacaWG0bGaey 4kaSIaamiAaaaakiabg2da9iaahwhadaahaaWcbeqaaiaadshaaaGc cqGHRaWkcaWGObGaaCODamaaCaaaleqabaGaamiDaaaakiabgUcaRm aalaaabaGaaGymaaqaaiaaikdaaaGaamiAamaaCaaaleqabaGaaGOm aaaakmaabmaabaWaaeWaaeaacaaIXaGaeyOeI0IaaGOmaiabek7aIb GaayjkaiaawMcaaiaahggadaahaaWcbeqaaiaadshaaaGccqGHRaWk caaIYaGaeqOSdiMaaCyyamaaCaaaleqabaGaamiDaiabgUcaRiaadI gaaaaakiaawIcacaGLPaaaaaa@518C@

The superscript t + α h indicates, for a general quantity, z , that is using parameter α MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3793@ :

z t + α h = ( 1 + α ) z t + h α z t

As the only special case, the acceleration is interpolated using parameter α m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacqaHXoqydaWgaaWcbaGaiqhGd2gaae qaaaaa@3573@ .

a t+ α m h = 1+ α m a t+h α m a t MathType@MTEF@5@5@+= feaahGart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeqabaqaciGa caGaaeqabaqaaeaadaaakeaacaqGHbWaaWbaaSqabeaacaWG0bGaey 4kaSIaeqySde2aaSbaaWqaaiaad2gaaeqaaSGaamiAaaaakiaaykW7 cqGH9aqpcaaMc8+aaeWaaeaacaaIXaGaaGPaVlabgUcaRiaaykW7cq aHXoqydaWgaaWcbaGaamyBaaqabaaakiaawIcacaGLPaaacaaMi8Ua aeyyamaaCaaaleqabaGaamiDaiabgUcaRiaadIgaaaGccaaMc8Uaey OeI0IaaGPaVlabeg7aHnaaBaaaleaacaWGTbaabeaakiaabggadaah aaWcbeqaaiaadshaaaaaaa@53CF@

From the above equations, it can be inferred that the Generalized Alpha method is controlled by four non-dimensional parameters ( α , β , γ , and α m ). With α m = 0.0 or α = α m = 0.0 , the method degenerates into the HHT- α method or Newmark- β method, respectively. The parameters can be selected as:

[ 1 3 < α 0 ] ; [ β 1 4 1 2 ( α m + α ) ] ; [ 1 α m < 1 2 ]

Given α m and α , the default values of γ and β which ensure the method is unconditionally stable (for linear problems), 2nd order accurate and with maximized high frequency dissipation:

[ γ = 1 2 ( α m + α ) ] ; [ β = 1 4 ( 1 α m α ) 2 ]

By default, α m = 0.0 and α = 0.05 . That is, the default scheme is HHT- α method.

The Generalized Alpha method is solved using Newton's method. In each iteration, Δ u t + h , the increment of displacement Δ u t + h = u t + h u h , is solved for. The subscript denotes the iteration count.

K * Δ u j + 1 t + h = f e x t t + h f i n t t + α h M a j t α m h C v j t + h
K * = 1 α m β h 2 M + ( 1 + α ) γ h β C + ( 1 + α ) K j
a j t α m h = α m a t + ( 1 α m ) a j t + h = α m a t + ( 1 α m ) ( 1 β h 2 Δ u j t + h 1 β h v t 1 2 β 2 β a t )
v j t + h = v t + h ( ( 1 γ ) a t + γ a j t + h )

Backward Euler Method

In the Backward Euler method, the equilibrium equation takes the following form:

f e x t t + α h = M a t + h + C v t + h + f i n t t + α h
v t + h = u t + h u t h
a t + h = v t + h v t h

Where,
M
Mass
C
Viscous damping
K
Stiffness matrix
f
Total force

Subscripts, e x t and i n t indicating external and internal force, respectively. Superscripts, t and t + h indicate the time when the variable is calculated, with t being the previous time and t + h being the current time when the displacement increment is being solved for. The Lagrangian coordinate is denoted by x , and the material time derivatives of x are denoted by v (velocity) and a (acceleration).

The Backward Euler method is also solved using Newton's method. In each iteration, Δ u t + h , the increment of displacement Δ u t + h = u t + h u h , is solved for. The subscript denotes the iteration count.

K * Δ u j + 1 t + h = f e x t t + h f i n t t + h M a j t + h C v j t + h
K * = 1 h 2 M + 1 h C + K j
a j t + h = 1 h 2 Δ u j t + h 1 h v t
v j t + h = 1 h Δ u j t + h

The only parameter for Backward Euler method is the time step h .

The Backward Euler method does not require the input of TCi fields on the TSTEP entry. The Alpha and Beta fields introduce subcase-dependent Rayleigh damping, so that the viscous damping matrix C in a particular subcase is C = α M + β K .

Typically, the Generalized Alpha method should be used for most Nonlinear Transient Analysis. In this method, numerical damping can be adjusted through the parameters α , β , γ and α m . In particular, non-zero α and α m introduce damping for high-frequency response components. However, Backward Euler method can be used for Static Analysis such as post-buckling solutions, since this method is dissipative and therefore highly stable.

Automatic Time Stepping

This Nonlinear Transient solution provides automatic time stepping based on Local Truncation Error (LTE).

The HHT- α method is used as an example to derive the method, and the other methods follow similar steps. Using HHT- α method, by Taylor's expansion, the LTE of displacement u from time t to t + h is:

e= 1 6 h 3 u τ t<τ<t+h MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakqaabeqaaiaadwgacqGH9aqpdaWcaaqaai aaigdaaeaacaaI2aaaaiaadIgadaahaaWcbeqaaiaaiodaaaGcceWH 1bGbaqaadaqadaqaaiabes8a0bGaayjkaiaawMcaaaqaaiaadshacq GH8aapcqaHepaDcqGH8aapcaWG0bGaey4kaSIaamiAaaaaaa@4300@

Using acceleration, e can be approximated as:

e= 1 6 h 2 Δa MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWGLbGaeyypa0ZaaSaaaeaacaaIXa aabaGaaGOnaaaacaWGObWaaWbaaSqabeaacaaIYaaaaOWaaqWaaeaa caqGuoGaaCyyaaGaay5bSlaawIa7aaaa@3C39@

For coupled multiple DOF system, the above absolute value is replaced by certain norm of Δ a . Here the mass-weighted norm is used, so that:

e= 1 6 h 2 Δ a T MΔa 1 2 MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWGLbGaeyypa0ZaaSaaaeaacaaIXa aabaGaaGOnaaaacaWGObWaaWbaaSqabeaacaaIYaaaaOWaaeWaaeaa caqGuoGaaCyyamaaCaaaleqabaGaamivaaaakiaah2eacaqGuoGaaC yyaaGaayjkaiaawMcaamaaCaaaleqabaWaaSGaaeaacaaIXaaabaGa aGOmaaaaaaaaaa@403F@

Error estimation in the previous equation requires normalization, as it depends on certain reference displacement, for example, initial displacement, as well as cut-off frequency. The normalized error estimation is:

e ˜ = e K Ω u R

Here, K Ω is the time-averaged Local Truncation Error (LTE) measure of a linear undamped single DOF oscillator with unit initial displacement and oscillating at approximately 10 cycles per second or with cut-off frequency Ω = ω h = 0.6 (where, ω is the natural frequency of the single DOF system). u R is the maximum value of the displacement norm over all the previous time steps.

u R = max 1 j n ( u j ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaadaqbdaqaaiaahwhadaWgaaWcbaGaam OuaaqabaaakiaawMa7caGLkWoacqGH9aqpdaWfqaqaaiGac2gacaGG HbGaaiiEaaWcbaGaaGymaiabgsMiJkaadQgacqGHKjYOcaWGUbaabe aakmaabmaabaWaauWaaeaacaWH1bWaaSbaaSqaaiaadQgaaeqaaaGc caGLjWUaayPcSdaacaGLOaGaayzkaaaaaa@47E4@

Where, u j MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaebbnrfifHhDYfgasaacH8srps0l bbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0R Yxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaabaqaciGa caGaaeqabaqaaeaadaaakeaacaWH1bWaaSbaaSqaaiaadQgaaeqaaa aa@33D4@ is the displacement at time step, j MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOAaaaa@36E6@ .

The norm is defined as:

u j = ( u j ) T M( u j ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaauWaaeaaca WH1bWaaSbaaSqaaiaadQgaaeqaaaGccaGLjWUaayPcSdGaeyypa0Za aOaaaeaadaqadaqaaiaahwhadaWgaaWcbaGaamOAaaqabaaakiaawI cacaGLPaaadaahaaWcbeqaaiaadsfaaaGccaWHnbWaaeWaaeaacaWH 1bWaaSbaaSqaaiaadQgaaeqaaaGccaGLOaGaayzkaaaaleqaaaaa@45A0@

The time steps are automatically adjusted based on the following conditions (TOL is the user-defined tolerance set on the TSTEP Bulk Data Entry), if:
  • e ˜ > TOL - Reject current step, cutback and redo the current step.
  • TOL > e ˜ > 0.5 * TOL - Accept current step, cutback the next step.
  • 0.5 * TOL > e ˜ > 1/16 * TOL - No changes.
  • 1/16 * TOL > e ˜ - The next time step is doubled.

The MREF continuation line on TSTEP entry can be used to control automatic time stepping, so that the time step h is adjusted according to the LTE of the current step. As shown above, when error is "large" when compared to the tolerance (TOL), h will be reduced by half and the current step is re-calculated. The maximum number of such operations within each step is controlled by the TN1 field. On the other hand, when h is "small" compared to the tolerance (TOL), h is requested to be increased, but only after TN2 contiguous steps with such a request, will h be actually increased.

Damping

This Nonlinear Transient solution sequence currently only provides support for Rayleigh Damping. The global Rayleigh damping parameters can be input using PARAM, ALPHA1 and PARAM, ALPHA2. Subcase-dependent Rayleigh damping parameters can be input using the Alpha and Beta fields on the TSTEP Bulk Data Entry. Material-dependent Rayleigh damping parameters can be input via the RAYL continuation line on the corresponding material entry (MAT1, MAT2, MAT8, MAT9, and MAT9OR Bulk Data Entries).

Baseline Correction

Baseline correction can be applied to correct the acceleration record to avoid displacement shift.

When acceleration is applied as input via SPCD, the integration of displacement results may diverge at the last point of simulation. Baseline correction can be used to correct the acceleration input. The BASELIN Subcase and BASELIN Bulk Data Entry pair can be used to apply baseline correction.

There are two methods by which Baseline correction can be applied. For both cases, the acceleration load is corrected by fitting the load to a polynomial with coefficients C 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ , C 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ , and C 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ .

Acceleration Method (METHOD=ACCE on BASELIN Bulk Data)

The baseline correction to the acceleration load is calculated as:

a b (t)= C 0 + C 1 t+ C 2 t 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGIbaabeaakiaacIcacaWG0bGaaiykaiabg2da9iaadoea daWgaaWcbaGaaGimaaqabaGccqGHRaWkcaWGdbWaaSbaaSqaaiaaig daaeqaaOGaamiDaiabgUcaRiaadoeadaWgaaWcbaGaaGOmaaqabaGc caWG0bWaaWbaaSqabeaacaaIYaaaaaaa@451B@

The corrected acceleration load is:

a c (t)=a(t) a b (t) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGJbaabeaakiaacIcacaWG0bGaaiykaiabg2da9iaadgga caGGOaGaamiDaiaacMcacqGHsislcaWGHbWaaSbaaSqaaiaadkgaae qaaOGaaiikaiaadshacaGGPaaaaa@43CC@

The coefficients of the polynomial are evaluated by minimizing the square of the corrected acceleration load.

S( C 1 , C 2 , C 3 )=min[ i=1 n ( a c ( t i ) ) 2 ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uamaabm aabaGaam4qamaaBaaaleaacaaIXaaabeaakiaacYcacaWGdbWaaSba aSqaaiaaikdaaeqaaOGaaiilaiaadoeadaWgaaWcbaGaaG4maaqaba aakiaawIcacaGLPaaacqGH9aqpciGGTbGaaiyAaiaac6gadaWadaqa amaaqahabaWaaeWaaeaacaWGHbWaaSbaaSqaaiaadogaaeqaaOGaai ikaiaadshadaWgaaWcbaGaamyAaaqabaGccaGGPaaacaGLOaGaayzk aaWaaWbaaSqabeaacaaIYaaaaaqaaiaadMgacqGH9aqpcaaIXaaaba GaamOBaaqdcqGHris5aaGccaGLBbGaayzxaaaaaa@527F@

Where,
C 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ , C 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ , and C 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@
Coefficients of the polynomial, which are unknowns to be calculated from the above minimization functions.
S MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaaaa@36CF@
Function of polynomial coefficients.
a ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyaiaacI cacaWG0bGaaiykaaaa@392E@
Input acceleration load.
a b ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGIbaabeaakiaacIcacaWG0bGaaiykaaaa@3A4B@
Baseline correction of the acceleration.
a c ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGIbaabeaakiaacIcacaWG0bGaaiykaaaa@3A4B@
Corrected acceleration after baseline correction.
t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaaaa@36CF@
Time associated with each acceleration data point in the record.
n MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaaaa@36CF@
Total number of acceleration data points in the record (points on TABLED# entry referenced on TLOAD1 entry).

Velocity Method (METHOD=VELO on BASELIN Bulk Data)

The baseline correction to the acceleration load is calculated as:

a b (t)= C 0 + C 1 ( t T 1 T 2 T 1 )+ C 2 ( t T 1 T 2 T 1 ) 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGIbaabeaakiaacIcacaWG0bGaaiykaiabg2da9iaadoea daWgaaWcbaGaaGimaaqabaGccqGHRaWkcaWGdbWaaSbaaSqaaiaaig daaeqaaOWaaeWaaeaadaWcaaqaaiaadshacqGHsislcaWGubWaaSba aSqaaiaaigdaaeqaaaGcbaGaamivamaaBaaaleaacaaIYaaabeaaki abgkHiTiaadsfadaWgaaWcbaGaaGymaaqabaaaaaGccaGLOaGaayzk aaGaey4kaSIaam4qamaaBaaaleaacaaIYaaabeaakmaabmaabaWaaS aaaeaacaWG0bGaeyOeI0IaamivamaaBaaaleaacaaIXaaabeaaaOqa aiaadsfadaWgaaWcbaGaaGOmaaqabaGccqGHsislcaWGubWaaSbaaS qaaiaaigdaaeqaaaaaaOGaayjkaiaawMcaamaaCaaaleqabaGaaGOm aaaaaaa@56BF@

The corrected acceleration load is:

a c (t)=a(t) a b (t) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGJbaabeaakiaacIcacaWG0bGaaiykaiabg2da9iaadgga caGGOaGaamiDaiaacMcacqGHsislcaWGHbWaaSbaaSqaaiaadkgaae qaaOGaaiikaiaadshacaGGPaaaaa@43CC@

The corrected velocity is:

v c (t)=v( T 1 )+ T 1 T 2 a c (t)dt MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamODamaaBa aaleaacaWGJbaabeaakiaacIcacaWG0bGaaiykaiabg2da9iaadAha caGGOaGaamivamaaBaaaleaacaaIXaaabeaakiaacMcacqGHRaWkda WdXbqaaiaadggadaWgaaWcbaGaam4yaaqabaGccaGGOaGaamiDaiaa cMcacaWGKbGaamiDaaWcbaGaamivamaaBaaameaacaaIXaaabeaaaS qaaiaadsfadaWgaaadbaGaaGOmaaqabaaaniabgUIiYdaaaa@4C96@

The coefficients of the polynomial are evaluated by minimizing the integration of square of the corrected velocity load.

S( C 1 , C 2 , C 3 )=min[ T 1 T 2 ( v c (t) ) 2 dt ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uamaabm aabaGaam4qamaaBaaaleaacaaIXaaabeaakiaacYcacaWGdbWaaSba aSqaaiaaikdaaeqaaOGaaiilaiaadoeadaWgaaWcbaGaaG4maaqaba aakiaawIcacaGLPaaacqGH9aqpciGGTbGaaiyAaiaac6gadaWadaqa amaapehabaWaaeWaaeaacaWG2bWaaSbaaSqaaiaadogaaeqaaOGaai ikaiaadshacaGGPaaacaGLOaGaayzkaaWaaWbaaSqabeaacaaIYaaa aOGaamizaiaadshaaSqaaiaadsfadaWgaaadbaGaaGymaaqabaaale aacaWGubWaaSbaaWqaaiaaikdaaeqaaaqdcqGHRiI8aaGccaGLBbGa ayzxaaaaaa@537A@

Where,
C 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ , C 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ , and C 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@
Coefficients of the polynomial, which are unknowns to be calculated from the above minimization functions.
S MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaaaa@36CF@
Function of polynomial coefficients.
a b ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGIbaabeaakiaacIcacaWG0bGaaiykaaaa@3A4B@
Baseline correction of the acceleration.
a c ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBa aaleaacaWGIbaabeaakiaacIcacaWG0bGaaiykaaaa@3A4B@
Corrected acceleration after baseline correction.
v ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamODaiaacI cacaWG0bGaaiykaaaa@3943@
Velocity load.
v c ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamODamaaBa aaleaacaWGJbaabeaakiaacIcacaWG0bGaaiykaaaa@3A61@
Corrected velocity after baseline correction.
t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaaaa@36CF@
Time associated with each acceleration data point in the record.
T 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ and T 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@
Initial and final total time.
For both continuation subcases or independent subcases, the baseline correction calculation is conducted for each subcase, and separate C 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ , C 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ , and C 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ coefficients are calculated for each subcase for baseline correction. In the case of continuation subcases, the values of T 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ and T 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaaIWaaabeaaaaa@37A5@ used in the calculation are still the total time.

Problem Setup

Input

For Nonlinear Transient Analysis, the loads and boundary conditions are defined in the Bulk Data Entry section of the input deck. The static loads and boundary conditions are defined in the Bulk Data Entry section of the input deck. They should be referenced in the Subcase Information Entry section using a SPC and DLOAD statement in a subcase, respectively.

For Nonlinear One-step Transient Thermal Stress Analysis (OSTTS), if the Nonlinear Transient Analysis subcase has a temperature load defined (via TEMPT), then a TEMP Subcase Information Entry referencing the TEMPT Bulk Data Entry can be directly specified in the nonlinear transient subcase.

Large Displacement Analysis can be activated via PARAM, LGDISP, 1, NLPARM(LGDISP), or NLCTRL(LGDISP) Subcase Information Entries.

To indicate that a nonlinear transient solution is required for any subcase, the combination of NLPARM (or NLCTRL) and TSTEP entries should be present (an alternative definition is to just use TSTEPNL instead of both NLPARM (or NLCTRL) and TSTEP). The NLPARM (or NLCTRL) Subcase Information Entry, in turn, points to the NLPARM (or NLCTRL) Bulk Data Entry, and the TSTEP entry points to the TSTEP Bulk Data Entry. If TSTEPNL is used in the subcase, then it points to the TSTEPNL Bulk Data Entry.

Using both TSTEPNL and NLPARM, or TSTEPNL and NLCTRL in a subcase is not supported (refer to the TSTEPNL Bulk Data Entry for more information). If NLPARM or NLCTRL is present in a Nonlinear Transient subcase, then a TSTEP entry is mandatory. However, if TSTEPNL entry is present instead of NLPARM or NLCTRL, then the TSTEP entry is not compulsory.

The NLOUT Bulk Data Entry and NLOUT Subcase Information Entry can be used to control incremental output. The NLADAPT Bulk Data Entry and Subcase Information Entry can be used to define parameters for time-stepping and convergence criteria.

The DLOAD Subcase Entry can reference corresponding DLOAD or TLOAD1/TLOAD2 Bulk Data Entries. The TLOAD1/TLOAD2 entries define time-dependent loading wherein the EXCITEID field can point to applied loads (DAREA, FORCEx, PLOADx, etc.), enforced motion (SPC/SPCD) or temperature loads (TEMP/TEMPD/TEMPSEC).

Example (NLPARM)

SUBCASE 10
   ANALYSIS = DTRAN
   SPC = 1
   DLOAD = 2
   NLPARM = 99
   TSTEP = 2
   NLOUT = 23
   IC = 12
.
.
BEGIN BULK
PARAM,LGDISP,1
NLPARM  99
TSTEP, 2

Example (NLCTRL)

SUBCASE 10
   ANALYSIS=DTRAN
   SPC = 1
   DLOAD = 2
   NLCTRL = 99
   TSTEP = 2
   NLOUT = 23
   IC = 12
.
.
BEGIN BULK
PARAM,LGDISP,1
NLCTRL  99 
TSTEP, 2

Detailed convergence check information for implicit nonlinear subcases, increments, and Newton-Raphson iterations are printed in the _nl.out file. Additional options for checking Nonlinear Analysis run information and monitoring progress of Nonlinear Implicit Analysis jobs are listed in Runtime Monitoring in the User Guide.

The NLPARM or NLCTRL data provides Nonlinear Analysis parameters. DLOAD inputs the time-dependent loading, and IC provides the initial condition (if absent, the initial velocity is set to zero). NLOUT issues the output controls. The TSTEP entry contains the Nonlinear Transient method parameters, damping, and time-stepping parameters.

The parameters NLDEBUG, CONT2TIE, NLDEBUG, NLP2NLC, and NLDEBUG, RMNLMAT are available to change the model in certain ways to aid in debugging.

Sensors can be defined using the SENSOR Bulk Data Entry and Subcase Information Entry pair to automatically stop a nonlinear analysis subject to specific conditions. See the SENSOR Bulk Data Entry for more information.

Subcase Continuation

Subcase continuation is supported between Nonlinear Transient subcases in the same model. It is also supported from a NLSTAT subcase to a Nonlinear Transient subcase.

With subcase continuation, subcase 2 follows subcase 1. When both are of Nonlinear Transient type, the initial condition (IC entry, if present) in subcase 2 is ignored. Instead, subcase 2 obtains its initial condition from the last time step of subcase 1.

On the other hand, if subcase 1 is of NLSTAT type, IC entry is allowed for subcase 2. In this situation you must also be cautious in defining the loading of subcase 1. If LOAD is used to define loading for subcase 1, its magnitude will be gradually ramped down to zero at the beginning of subcase 2. If DLOAD is used in subcase 1, the load will not continue into subcase 2. Instead, it becomes zero immediately when entering subcase 2. In other words, for NLSTAT(LOAD) + Nonlinear Transient, the load in the 2nd subcase is the combination of DLOAD of Nonlinear Transient and the (ramping down) LOAD of NLSTAT. Whereas, for NLSTAT(DLOAD) + Nonlinear Transient, the load in the 2nd subcase is only from DLOAD of Nonlinear Transient.

One typical use case for subcase continuation is the post-buckling analysis. One can load the structure to a point close to buckling using a NLSTAT Subcase, followed by a Nonlinear Transient subcase (using the Backward Euler method to improve stability) for the rest of buckling process.

Subcase continuation between Small and Large Displacement Nonlinear Transient Analyses (or vice-versa) is not allowed.

Output

The typical output entries (DISPLACEMENT, VELOCITY, and ACCELERATION) can be used to request corresponding output for Nonlinear Transient Analysis. The NLOUT Subcase and Bulk Data Entries can be used to request intermediate results.

Energy output is available by default and additional control is available via NLENRG Bulk/Subcase pair.

If automatic time-stepping is not active, then the typical convergence results are output for each time step similar to NLSTAT Subcases. However, if automatic time-stepping is activated, then, in addition to the convergence results, an extra line is added to each time step.
Figure 1.

time_stepping

The fields show the current time step number (itstp), current time (time), norm of acceleration error (norm(da)), the norm of the reference displacement (u_ref), the acceleration error (err_da) and adjustment request (adjustment), respectively.

The adjustment request field shows the following syntax:
Adjustment Printing (.out file)
Information
Reduce Next
Indicates a reduction of next time step for h
Enlarge Next
Indicates an increase of next time step for h
No Change
Indicates no change in h for current and subsequent step.
Cutback
Indicates reduction of current step. The current step is then redone.
Note:
  1. The same limitations that apply to Nonlinear Static Analysis also apply to Nonlinear Transient Analysis.
  2. The NLADAPT entry is recommended to be used to control the maximum and minimum time-step in Nonlinear Transient Analysis.
  3. Linear and large rotation DMIG (Direct Matrix Input at Grid points) are supported in Nonlinear Transient Analysis.

Plot Convergence During Analysis

The convergence plot can be obtained for all nonlinear iterations across all subcases for nonlinear static and nonlinear transient analysis, using the Altair Compute Console (ACC).

The progress can also be monitored for all the optimization iterations.

For more details, refer to Plot Convergence During Analysis in the Runtime Monitoring section.