A kinetics model describing macroscopic behavior of martensitic variants rearrangement in ferromagnetic shape memory alloys is proposed. The model is based on a tensor description of thermodynamic continuum mechanics taking into account magnetomechanical coupling on the rearrangement process. A generalized thermodynamic driving force which rearranges the martensitic variants is derived. The thermodynamic evolution equations of the kinetics are obtained by means of the dissipation, which is induced by the twin boundary motion during rearrangement process, reaching a maximum. The ability of the model to characterize the macroscopic stress-strain behavior and magnetic field-induced strain of a ferromagnetic shape memory alloy rod under both magnetic field and stress field is demonstrated. The application of the theoretical kinetics model shows the good agreement with experimental observation.