A dynamic free energy hysteresis model in magnetostrictive actuators is presented. It is the free energy hysteresis model coupled to an ordinary different equation in an unusual way. According to its special structure, numerical implementation method of the dynamic model is provided. The resistor parameter in the dynamic model changes according to different frequency ranges. This makes numerical implementation results reasonable in the discussed operating frequency range. The validity of the dynamic free energy model is illustrated by comparison with experimental data.