where \(\mathbf{M}\) and \(\mathbf{C}\) is the mass matrix and damping matrix, respectively. The \(\mathbf{Q}\) and \(\mathbf{F}\) are the internal force vector, and external force vector, respectively. Note that, the mass density is consistent in the structural analysis, thus the mass matrix keep constants regardless of the deformation in non-linearity.
In the Newmark-\(\beta\) method, the displacement, velocity, and acceleration at the each time incremental \(\triangle t\) are
where \(\gamma\) and \(\beta\) are time solver parameters. Given the specific values, the numerical algorithm becomes linear acceleration method, or the trapezoid rule.
substituting equations \(\eqref{eq:ch5_time_solver_imp2}\) and \(\eqref{eq:ch5_time_solver_imp3}\) into equation \(\eqref{eq:ch5_time_solver_imp1}\), the following equation can be acquired
when we use linear stiffness matrix \(\mathbf{K}_{L}\) for a linear problem, the equation above becomes \(\mathbf{Q}(t+\triangle t)=\mathbf{K}_{L}\mathbf{U}(t+\triangle t)\). Substituting this term into the equation (), we have
In the analysis practice, the acceleration and velocity boundary conditions are imposed. Then the displacement of the following equation can be derived from equation \(\eqref{eq:ch5_time_solver_imp1}\).
where \(u_{is}(t+\triangle t)\) is the nodal displacement at time \(t+\triangle t\), \(\dot{u}{}_{is}(t+\triangle t)\) is the nodal velocity, \(\ddot{u}{}_{is}(t+\triangle t)\) is the nodal acceleration, i is the degree of freedom per node, s is the node number.
The mass and damping terms are treated as follows
The lumped mass matrix is used at most of cases in this program.
The damping matrix is treated using Rayleigh algorithm \(\mathbf{C}=R_{m}\mathbf{M}+R_{k}\mathbf{K}_{L}\).
where the displacement at the time \(t+\triangle t\) and \(t-\triangle t\) can be expressed by the Taylor's expansion at time t with the second order truncation.
For the linear problem, we also have condition \(\mathbf{Q}(t)=\mathbf{K}_{L}\mathbf{U}(t)\) for equation. Finally, the displacement at \(t+\triangle t\) is: