For steady state dynamics calculations the steady state response of the system to a harmonic excitation is again assumed to be a linear combination of the lowest eigenmodes. To this end, the eigenvalues and eigenmodes must have been calculated, either in the same run, or in a previous run. At the end of a frequency calculation this data, including the stiffness and mass matrix, is stored in binary form in a .eig file, provided the STORAGE=YES option is activated on the *FREQUENCY or *HEAT TRANSFER,FREQUENCY card. This file is read at the beginning of file steadystate.c.

In file steadystate.c the response is calculated in an explicit way, for details the reader is referred to [19]. Modal damping is allowed in the form of Rayleigh damping. Within file steadystate the following routines are used:

- tempload, to calculate the instantaneous loading
- rhs, to determine the external force vector of the system
- results, to calculate all displacements, stresses and/or any other variables selected by the user

Notice that if nonzero boundary conditions are prescribed (base loading, e.g for earthquake calculations) the stiffness matrix of the system is used to calculate the steady state response to these nonzero conditions. It serves as particular solution in the modal dynamic solution procedure.