Neste trabalho é estudada a convergência no método dos elementos finitos (MEF) quando da sua aplicação a dispositivos eletromagnéticos com características não lineares, ou seja, constituídos de materiais ferromagnéticos com curva de magnetização BxH que apresenta o fenômeno da saturação magnética. O método de Newton-Raphson é o método numérico utilizado para a resolução do sistema de equações não linear gerado. São abordadas versões desse método que utilizam fator de relaxação, bem como técnicas para a otimização desse fator, com o intuito de acelerar a sua convergência, mostrando sua efetividade. Aborda-se também a influência da escolha do modelo para a característica B(H) dos materiais ferromagnéticos, tanto em termos de desempenho computacional quanto de precisão. Duas classes de modelos para a curva vxB², bem como para sua derivada dv/dB², são analisados, quais sejam: modelos analíticos e modelos interpoladores. Dentre os modelos analíticos destacam-se os modelos de O'Connor, Marrocco e de Brauer (exponencial). Os modelos interpoladores utilizados foram de Akima e Spline cúbica. Além do tratamento de materiais com características não lineares, também se realizou um breve estudo da aplicação de materiais do tipo ímã permanente e sua formulação pelo MEF. Outro aspecto avaliado foi a influência de diferentes métodos de resolução do sistema linear de equações algébricas na convergência do método de Newton-Raphson, como o ICCG, BiCGStab, decomposição LU e Wavelet-based Multigrid. Os resultados obtidos mostraram a potencialidade desta última técnica em termos de ganho de desempenho global da solução.

This work presents a study on the convergence of non-linear finite analysis applied to the modeling of electromagnetic devices with plane symmetry. Three aspects have been thoroughly investigated which have influence on the convergence performance, namely: the relaxation factor in the Newton-Raphson algorithm, the model of the magnetization curve and the choice of the linear solver. It is shown that convergence acceleration of the Newton-Raphson algorithm can be achieved with the use of a relaxation factor, along with numerical techniques for optimizing this factor, such as Line Search or the Brent method. Several models for the magnetization curve, both analytical (O'Connor, Marrocco, Brauer) and interpolation-based (Akima and Cubic Splines), have been analised and comparisons in terms of their computational efficiency, as well as accuracy, are presented in order to give some guidance to the choice of the model. The influence of several linear solvers in the convergence of the Newton-Raphson algorithm has finally been investigated. The following linear solvers have been compared: LU decomposition , ICCG, BiCGStab and Wavelet-based Multigrid. The results show the superiority of the latter to improve the overall computational performance of the solution.