Dynamic Analysis Method
Dynamic Analysis Method
In this section, the dynamic problem analysis method with a direct time integration method applied is described. As presented below, with this development code, it is possible to perform time history response analysis by an implicit or explicit method.
Formulation of the implicit method
Focusing on dynamic problems, the direct time integration method was applied to solve the following equation motion indicated:
where and the mass matrix, is the damping matrix, and is the internal stress vector, and is the external force vector. This software does not consider the changes in mass; thus, the mass matrix is non-linear and constant regardless of deformation.
The displacement within time increment , and the change in speed and acceleration are approximated with the Newmark- method, as expressed in Eq. and Eq.:
where and are the parameters of the Newmark- method.
if and have the following values, it coincides with the linear acceleration method or the trapezoidal rule.
As it is already known, when and are substituted into the following values, it will match the linear acceleration method, or trapezoid rule.
- , (Linear acceleration method)
- , (Trapezoid rule)
If Eq. and Eq. are substituted into Eq., the following equation is obtaind:
is linear stiffness matrix for linear problem; thus, . If this equation is substituted into the equation above, the following equation is obtained:
In the portion, where the acceleration is specified as a geometric boundary condition, the displacement of the following equation is obtained from Eq..
Similarly, if the speed is specified, the displacement of the following equation is obtained from Eq.:
- is the nodal displacement at ,
- is the nodal speed at time
- is the nodal acceleration at time
- is the nodal degree-of-freedom (DOF)
- is the Node number.
Furthermore, the mass and damping terms were handled as follows:
(1) Handling of mass term
The mass matrices are handled as concentrated mass matrices.
(2) Handling of damping terms
The damping terms are handled as Rayleigh damping expressed in Eq..
where and are parameters of Rayleigh damping.
Formulation of Explicit Method
The explicit method is based on the motion equation at time expressed in the following equation:
where the displacement at time and that at time are expressed by the Taylor expansion at time . If it is expanded up to the secondary term with , it becomes
From the difference and sum of Eq. and Eq., the following equation is obtained:
If Eq. and Eq. are substituted into Eq., the following equation is obtained:
For linear problems, specifically, , the above equation becomes
In this case, if mass matrix is set as a concentrated mass matrix, and the damping matrix as the proportional damping matrix , Eq. eliminates the requirement of solving operations for simultaneous equations.
Therefore, from Eq., can be determined by the following equation: