How to extract Stiffness matrix(s) for Non-linear FE Dynamic Implicit case in ABAQUS

As discussed in the previous posts, it does not seem to be possible, as it only outputs for linear perturbation analysis which can include initial stresses/stiffness effects, but that is not the same as obtaining the FE matrices from a full nonlinear analysis using a NR iterative solution

I don’t think it’s possible (at least without some kind of subroutine). That’s because *Matrix generate keyword will only give you the stiffness matrix from the end of the step.

Well it might be sort of possible if you have run your general nonlinear analysis up to the point to where you want to extract your FE matrices, and then input that last step/iteration as an initial condition to a linear perturbation analysis where you can output them (so one should get the correct stiffness matrix which was generated in the last increment of the NR iterations prior to the linear perturbation, because that will be used in the pert. analysis, not sure about other vectors).

Try it out on a simple model