A free energy hysteresis model for giant magnetostrictive actuators was briefly presented, and the free energy hysteresis model was analysed. Four different numerical calculation methods are advanced. Computational precision and computational efficiency of the results obtained by using the four different numerical calculation methods were compared. These results also were compared to the calculation data and the experimental data in the published papers. It states that the method in which the approximately Gauss-Legendre integral method is used to integral discretization and matrix representation is used to implement integral kernel function not only increases computational efficiency but also guarantee computational precision. It is the optimum method for a free energy hysteresis model in giant magnetostrictive actuators.