In these attached papers, the authors said that it can be used frozen permeability (FP) to determine PM flux linkage, Ld and Lq corresponding to rotor position, Fig. 11 and 12 of the first paper, and Fig. 2 of the second paper. But in my knowledge, the FP usually used to separate average PM and reluctance parameters.

Does anyone can suggest the method to determine above parameters in the variation with rotor position, not the average ones like some existed questions in this website? And how to do these in FEA softwares like JMAG, Ansys Maxwell, etc?

Thank you.