How to create stiffness matrix for a model both with load and without load?

I am solving a FDU unit, where the requirement is to extract stiffness matrix in bearing locations which can be imported to Kisssoft or Sabr for bearing evaluation.