I have a .inp file which roughly looks like:

Step1:

*** Apply some loads ***

Step2:

*Matrix Generate, STIFFNESS
*Matrix Output, STIFFNESS, FORMAT=COORDINATE

Step3:

*** Apply more loads ***

I have a couple questions which are unanswered by the documentation I have read.

In the context of my .inp file above, would the generated stiffness matrix in step 2 be valid for the state of model after the loads in step 1 have been applied?

Also, I want to make sure that I am interpreting the global stiffness matrix correctly. Let's say the model is 3D and there are N nodes in my mesh. The stiffness matrix is always 3N x 3N (I have confirmed this experimentally). The degrees of freedom it embodies are the displacements (in terms of the global coordinate system) of each node in the X, Y, and Z directions right? It does explicitly embody rotational degrees of freedom right?

Abaqus ​​​​​​​