Este trabalho aborda a modelagem numérica de dispositivos eletromagnéticos incluindo o fenômeno de histerese em três dimensões através da formulação da permeabilidade diferencial. esta formulação possibilita o acoplamento com circuitos elétricos e o tratamento de materiais anisotrópicos e histeréticos. o modelo de histerese utilizado é o vetorial de jiles-atherton. são discutidos aspectos sobre a implementação bem como melhorias para a convergência do método. propõe-se uma mudança na discretização temporal substituindo o método de euler por crank-nicolson. um benchmark proposto pela ics (international compumag society) é utilizado para validação do método, além da comparação com outras formulações já implementadas. o software é incorporado ao sistema feecad, desenvolvido pelo grucad (grupo de concepção e análise de dispositivos eletromagnéticos) em colaboração com outros grupos de pesquisa, que já conta com diversos módulos implementados utilizando diferentes formulações.