The difference between the stiffness matrix at different harmonic indices is the constraint equations created to enforce the harmonic index shape. These constraint equations are created during solution. A possible way to create a model for a harmonic index greater than 0 is to issue the command cycopt,hindex,N in a commands object. Then the system matrices will be only for harmonic index N.
You are navigating away from the AIS Discovery experience