Abstract: A novel method for calculating the magnetic stiffness matrix was proposed for the numerical analysis of the magneto-elastic stability of complicated current-carrying structures aim- ing for application in the magneto-elastic behavior of the tokamak system. A code based on the proposed method was developed and applied to the numerical analysis of two typical current- carrying structures. The good consistency of the numerical and analytical results validated the proposed method and the related numerical code.