Vol.

17, 1, 219248 (2001)

Revista Internacional de
Metodos Numericos para
C alculo y Dise no en Ingeniera
El Metodo de los Elementos Finitos en problemas
electromagneticos: planteamiento y aplicaciones
Josep Sarrate y Ramon Clariso
Universidad Politecnica de Catalu na
ETS de Ingenieros de Caminos, Canales y Puertos
Departamento de Matematica Aplicada III
Jordi Girona 13, Edicio C2, Campus Norte
08034 Barcelona, Espa na
Tel.: 34-93-401 69 11, Fax: 34-93-401 18 25
e-mail: [email protected]
Resumen
En el presente trabajo se analizan y comparan diferentes alternativas para resolver mediante el Metodo de los
Elementos Finitos (MEF) problemas electromagneticos. Asimismo, se repasan diversos conceptos fsicos a n
de justicar los posteriores desarrollos numericos. Las formulaciones para la resolucion de dichos problemas
se clasican en dos grandes grupos: 1. formulaciones basadas en potenciales escalares y/o vectoriales y
2. formulaciones basadas en la resolucion de las ecuaciones de Maxwell. Puesto que la mayora de aplicaciones
industriales requieren un calculo preciso de las fuerzas electromagneticas, en este trabajo tambien se analizan
diferentes metodos para calcularlas. Finalmente, se presentan diversos ejemplos numericos que ilustran
algunos de los metodos presentados. Estos ejemplos incluyen desde problemas academicos hasta aplicaciones
industriales.
Palabras clave: Palabras clave:
Ecuaci on de Maxwell, electromagnetismo, fuerzas magneticas, Metodo de los Elementos
Finitos.
FINITE ELEMENT METHOD ON ELECTROMAGNETICS: AN OVERVIEW AND
APPLICATIONS
Summary
This work reviews and compares dierent methods of solving electromagnetic problems by means of the
Finite Element Method (FEM). In order to justify the numerical developments, a brief review on the physics
of the problem is given. Numerical techniques are classied into two main categories: 1. methods that
use scalar and vectorial potentials and 2. methods that solve directly the Maxwell equations. Moreover,
in engineering applications, the design of electromagnetic devices requires accurate calculation of magnetic
forces. Therefore, dierent methods for computing the magnetic forces are analyzed in this paper. Finally,
several magnetostatic examples, that range from academic problems to industrial applications, are included
in order to illustrate some of the described methods.
Keywords: Keywords:
Maxwell equations, electromagnetism, magnetic forces, nite element method.
c Universitat Polit`ecnica de Catalunya (Espa na). ISSN: 02131315 Recibido: Julio 2000
220 J. Sarrate y R. Clariso
INTRODUCCI

ON
En la actualidad, la inferencia de los fen omenos electromagneticos en los procesos in-
dustriales es muy importante. En particular, la inducci on electromagnetica se utiliza para
medir, calentar, modicar, manipular y formar una amplia gama de materiales. Es m as, en
las ultimas decadas, la demanda de nuevas aplicaciones m as potentes, versatiles y economicas
(como por ejemplo, la levitaci on de trenes mediante superconductores, el dise no de aviones
capaces de minimizar el scattering de una onda incidente, el desarrollo de antenas de alta po-
tencia y precisi on, la utilizaci on de l aseres en sistemas de telecomunicaciones opticas, entre
otras) ha inducido un cambio radical en el dise no de diferentes aparatos y dispositivos, intro-
duciendo nuevos materiales de mayor conductividad electrica o capaces de soportar mayores
tensiones estructurales. En este sentido, cabe resaltar la introducci on de superconductores
en la fabricaci on de motores
1
o imanes para un reactor de fusi on
2
.
La necesidad de calcular rapidamente y con precisi on estas modicaciones, que gen-
eralmente tambien introducen cambios importantes en la geometra de los dispositivos,
ha creado la necesidad de desarrollar nuevos metodos de calculo mucho m as veloces,
versatiles y precisos. Las formulaciones para la resolucion de problemas en electro-
magnetismo computacional se pueden clasicar en dos grandes grupos, seg un el tipo de
ecuaci on que resuelven: 1. las basadas en potenciales escalares y/o vectoriales
3,4,5,6,7,8,9
y
2. las basadas en la resoluci on directa de las ecuaciones de Maxwell
6,10,11,12,13
. Las formula-
ciones basadas en potenciales son ampliamente utilizadas en problemas estaticos, arm onicos
y especialmente en problemas de corrientes de Foucault. Cabe destacar su relevancia en la
resolucion de problemas en electrotecnia. Adem as, los potenciales vectoriales presentan la
gran ventaja de vericar autom aticamente las ecuaciones de la divergencia (leyes de Gauss).
Sin embargo, presentan una perdida notable de precisi on en el calculo de las intensidades de
campo electrico y magnetico debido a la propagaci on de errores en la diferenciaci on numerica.
Ademas, su utilizaci on requiere la denici on de una condici on de contraste (gauge) adecuada
que depende de la formulaci on utilizada. Las formulaciones basadas en la resoluci on directa
de las ecuaciones de Maxwell se utilizan tanto en problemas est aticos como en transitorios.
En este sentido, se debe mencionar que este tipo de formulaciones ocupa un lugar desta-
cado en fen omenos de propagaci on de ondas como el dise no de antenas y dispositivos de
microondas, entre otros. Estas formulaciones permiten hallar directamente las intensidades
de campo electrico y magnetico sin necesidad de utilizar derivaci on numerica, obteniendo,
generalmente, resultados m as exactos.
Uno de los aspectos fundamentales en el dise no y optimizaci on de equipos electro-
magneticos, como sensores, motores o electroimanes, es el calculo de la fuerza magnetica.
En este sentido, el acoplamiento entre problemas mecanicos y electromagneticos es una de
las areas de mayor interes industrial en el electromagnetismo computacional. En efecto,
la mayora de aparatos electricos contienen una parte cuyo movimiento es debido y esta
gobernado por las fuerzas magneticas. Por consiguiente, el resultado del an alisis de este
problema acoplado permite al ingeniero calcular los desplazamientos, visualizar las defor-
maciones, evaluar las vibraciones inducidas y ensayar diferentes soluciones a n de mejorar
el rendimiento del equipo.
En general, la simulaci on y evaluaci on de las fuerzas magneticas se basa en el calculo
previo del campo magnetico mediante el MEF u otra tecnica numerica. Por consiguiente, en
este trabajo se revisan y analizan diversas formulaciones para calcular las fuerzas magneticas
en funci on de los valores nodales de los campos
14,15,16,17,18,19,20,21,22,23,24,25,26
. Asimismo,
se presenta una clasicacion de los metodos para evaluar dichas fuerzas a partir de las
formulaciones anteriores. Se demuestra que estos metodos son teoricamente equivalentes.
Sin embargo, tambien se constata que su precision depende fuertemente de los valores de la
permeabilidad relativa del medio y de los valores normales y tangenciales a la supercie de
los campos previamente calculados.
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 221
La estructura del artculo es la siguiente. En la segunda seccion se plantean detallada-
mente el dominio, las ecuaciones y las condiciones de contorno que rigen el problema electro-
magnetico a resolver. En la tercera secci on se presenta una clasicacion de las formulaciones
mas utilizadas para la resoluci on de problemas de electromagnetismo. Seguidamente, en la
cuarta seccion se analizan y clasican los diferentes metodos de evaluar las fuerzas electro-
magneticas a partir del calculo previo de los campos electromagneticos mediante el MEF.
En la quinta secci on se utilizan algunas de las formulaciones presentadas anteriormente en
la resolucion de diferentes problemas pr acticos que van desde ejemplos academicos hasta
aplicaciones industriales. Finalmente, se presentan las conclusiones de este trabajo.
PLANTEAMIENTO DEL PROBLEMA
Considerese un dominio cualquiera compuesto por dos medios disjuntos,
1
y
2
, cuya
intercara es
12
(Figura 1). Sea su contorno. De acuerdo con los dos tipos de condiciones
de contorno que se imponen (por simplicidad y sin perdida de generalidad se supondr an
homogeneas) este tambien se divide en dos partes,
1
y
2
, tales que

1

2
=

1

2
=
Figura 1. Esquema del dominio y contornos para un problema general de electromagnetismo
Los campos electromagneticos existentes en el dominio estan gobernados por el siguien-
te sistema de ecuaciones diferenciales en derivadas parciales junto con las correspondientes
condiciones de contorno sobre la intercara y sobre las fronteras exteriores
5,27,6
EEEEEEE EEEEE EE +
B BBBBBB BBBBB BB
t
= KKKKKKK KKKKK KK
eqv
en
HHHHHHH HHHHH HH
D DDDDDD DDDDD DD
t
JJJJJJJ JJJJJ JJ = JJJJJJJ JJJJJ JJ
imp
en
DDDDDDD DDDDD DD =
imp
en
BBBBBBB BBBBB BB = 0 en
(1)
222 J. Sarrate y R. Clariso
nnnnnnn nnnnn nn EEEEEEE EEEEE EE = 0000000 00000 00 sobre
1
nnnnnnn nnnnn nn BBBBBBB BBBBB BB = 0 sobre
1
nnnnnnn nnnnn nn HHHHHHH HHHHH HH = 0000000 00000 00 sobre
2
nnnnnnn nnnnn nn DDDDDDD DDDDD DD = 0 sobre
2
nnnnnnn nnnnn nn
1
EEEEEEE EEEEE EE
1
+nnnnnnn nnnnn nn
2
EEEEEEE EEEEE EE
2
= 0000000 00000 00 sobre
12
nnnnnnn nnnnn nn
1
BBBBBBB BBBBB BB
1
+nnnnnnn nnnnn nn
2
BBBBBBB BBBBB BB
2
= 0 sobre
12
nnnnnnn nnnnn nn
1
DDDDDDD DDDDD DD
1
+nnnnnnn nnnnn nn
2
DDDDDDD DDDDD DD
2
=
s
sobre
12
nnnnnnn nnnnn nn
1
HHHHHHH HHHHH HH
1
+nnnnnnn nnnnn nn
2
HHHHHHH HHHHH HH
2
= JJJJJJJ JJJJJ JJ
s
sobre
12
donde EEEEEEE EEEEE EE y HHHHHHH HHHHH HH son las intensidades de campo electrico y magnetico respectivamente, DDDDDDD DDDDD DD y
BBBBBBB BBBBB BB son las densidades de ujo electrico y magnetico respectivamente, JJJJJJJ JJJJJ JJ es la densidad de
corriente,
imp
es una fuente de densidad de carga electrica impuesta, JJJJJJJ JJJJJ JJ
imp
es una fuente
impuesta de densidad de corriente electrica, KKKKKKK KKKKK KK
eqv
es una densidad de corriente magnetica
equivalente debida a truncamiento de dominios o discontinuidades de medios,
s
es una
densidad de carga supercial y JJJJJJJ JJJJJ JJ
s
es una corriente electrica supercial en la intercara. Se
debe resaltar que los valores impuestos
imp
, JJJJJJJ JJJJJ JJ
imp
y KKKKKKK KKKKK KK
eqv
deben cumplir ciertas condiciones
de compatibilidad
11
.
Las cuatro primeras ecuaciones de (1) constituyen las ecuaciones de Maxwell: ley de
Faraday, ley de MaxwellAmpere, ley de Gauss para el campo electrico y ley de Gauss para
el campo magnetico.
El sistema (1) queda completamente denido al especicar las ecuaciones constitutivas
de los materiales que conguran los medios
1
y
2
DDDDDDD DDDDD DD = (EEEEEEE EEEEE EE)EEEEEEE EEEEE EE
BBBBBBB BBBBB BB = (HHHHHHH HHHHH HH)HHHHHHH HHHHH HH
JJJJJJJ JJJJJ JJ = (EEEEEEE EEEEE EE)EEEEEEE EEEEE EE
(2)
donde es la permitividad, es la permeabilidad y es la conductividad. Frecuentemente,
la segunda expresi on de (2) se escribe como HHHHHHH HHHHH HH = BBBBBBB BBBBB BB, siendo = 1/ la reluctividad del
medio. En general, estos par ametros constitutivos pueden depender de la posici on, de la
orientaci on, del tiempo y de la intensidad del campo correspondiente. Si dichos par ametros
no dependen de la intensidad de campo, se obtiene un comportamiento lineal del medio
(que conducira en el problema discreto a la resoluci on de sistemas de ecuaciones lineales).
Si adem as el medio es isotropo y homogeneo, dichos par ametros son constantes. Por el
contrario, si dependen de la intensidad de campo se obtiene un comportamiento no lineal
(que en el problema discreto implicar a la resoluci on de sistemas no lineales de ecuaciones).
Esta relacion no lineal generalmente se conoce experimentalmente o mediante la teora
microscopica de la materia. Por ejemplo, para materiales ferromagneticos la permeabilidad
depende fuertemente de la intensidad de campo y de su historia (curvas de histeresis) como
muestra la Figura 2. En el modelado de estos materiales, el valor del campo coercitivo
desempe na un papel fundamental. Por ultimo, debe resaltarse que puede existir otro tipo
de no linealidad, generada por ejemplo mediante fuentes de carga no lineales.
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 223
Figura 2. Curva de histeresis para un material ferromagnetico
CLASIFICACI

ON DE LAS FORMULACIONES PARA LA RESOLUCI

ON DE
PROBLEMAS ELECTROMAGN

ETICOS
Las formulaciones para la resoluci on de problemas de campo en electromagnetismo
computacional se pueden clasicar seg un el tipo de ecuaci on diferencial en derivadas parciales
que resuelven. En general, se pueden dividir en dos grandes grupos: las formulaciones que
utilizan potenciales escalares y/o vectoriales y las formulaciones basadas en la resoluci on de
las ecuaciones de Maxwell u otro sistema equivalente.
Basadas en potenciales escalares y/o vectoriales
Las formulaciones basadas en la utilizacion de potenciales son ampliamente utilizadas
en el calculo numerico de campos electromagneticos. Basicamente se utilizan en problemas
estaticos o arm onicos y especialmente en problemas de corrientes de Foucault
4,3
. Esto
se debe a dos motivos: por una parte, en problemas estaticos y arm onicos, los terminos
temporales de las ecuaciones de Maxwell (1) pueden tratarse de forma algebraica; y por otra
parte, las ecuaciones en terminos de potenciales vectoriales no incluyen de forma explcita las
ecuaciones en divergencia (leyes de Gauss) puesto que estas se verican automaticamente.
Un ejemplo de esta aproximacion para problemas de corrientes de Foucault en dominios
acotados puede hallarse en la referencia
3
, para dominios no acotados ver la referencia
8
. En
ambos trabajos y con el objetivo de asegurar la unicidad de la soluci on se presentan varias
formulaciones, en las que se introducen diferentes tipos de potencial dependiendo del n umero
de materiales que se consideren. La nomenclatura utilizada para distinguir las diferentes
formulaciones se basa en los siguientes criterios:
1. los potenciales vectores se denotan mediante letras may usculas y en negrita, (AAAAAAA AAAAA AA), mientras
que para los potenciales escalares se utilizan smbolos romanos o griegos (A, , . . .);
2. los potenciales se ordenan seg un el medio al que est an asignados, separando cada medio
por un gui on ();
3. el primer grupo de potenciales se referira siempre al medio conductor, mientras que el
resto de potenciales estaran asignados a las regiones no conductoras.
La primera formulaci on se denomina AAAAAAA AAAAA AA, V AAAAAAA AAAAA AA, donde en el medio conductor se utiliza
un potencial vector AAAAAAA AAAAA AA y un potencial escalar V , mientras que en la regi on no conductora
solo se emplea el potencial vector. En este caso, los potenciales se denen como
BBBBBBB BBBBB BB = AAAAAAA AAAAA AA (3)
EEEEEEE EEEEE EE =
A AAAAAA AAAAA AA
t
V (4)
224 J. Sarrate y R. Clariso
Este tipo de formulaciones se aplica generalmente a problemas compuestos por dos
medios: un material conductor
1
, donde aparecen las corrientes de Foucault, rodeado
por otro medio no conductor
2
, que puede contener fuentes de corriente pero que esta libre
de corrientes de Foucault (Figura 3). La intercara entre ambos medios se denota por
12
.
El contorno exterior de
2
se divide en dos partes seg un el tipo de condicion de contorno
que se impone:
B
donde se prescribe componente normal de la densidad de ujo magnetico
nula (BBBBBBB BBBBB BB nnnnnnn nnnnn nn = 0) y
H
donde se prescribe componente tangencial de la intensidad de campo
magnetico nula (HHHHHHH HHHHH HH nnnnnnn nnnnn nn = 0). Es importante resaltar la relaci on existente entre el presente
enunciado y el problema general planteado en la segunda secci on (Figura 1). Puesto que
2
es un medio libre de corrientes de Foucault, EEEEEEE EEEEE EE no aparece en las condiciones de contorno
sobre la frontera exterior. Por consiguiente, los contornos
1
y
2
de la Figura 1 se reducen
a los contornos
B
y
H
respectivamente y ambos problemas coinciden.
Figura 3. Dominio y condiciones de contorno para un problema de corrientes de Foucault
compuesto por dos medios
Reescribiendo el sistema (1) en terminos de los potenciales (3) y (4) se obtiene la siguiente
formulacion, ver detalles en la referencia
3
AAAAAAA AAAAA AA AAAAAAA AAAAA AA +
A AAAAAA AAAAA AA
t
+V = 0 000000 00000 00 en
1

A AAAAAA AAAAA AA
t
V

= 0 en
1
AAAAAAA AAAAA AA AAAAAAA AAAAA AA = JJJJJJJ JJJJJ JJ
imp
en
2
nnnnnnn nnnnn nn AAAAAAA AAAAA AA = 0 000000 00000 00 sobre
B
AAAAAAA AAAAA AA = 0 sobre
B
AAAAAAA AAAAA AA nnnnnnn nnnnn nn = 0 000000 00000 00 sobre
H
n AAAAAAA AAAAA AA = 0 sobre
H
AAAAAAA AAAAA AA
1
= AAAAAAA AAAAA AA
2
sobre
12

1
AAAAAAA AAAAA AA
1
nnnnnnn nnnnn nn
1
+
2
AAAAAAA AAAAA AA
2
nnnnnnn nnnnn nn
2
= 0 000000 00000 00 sobre
12

1
AAAAAAA AAAAA AA
1

2
AAAAAAA AAAAA AA
2
= 0 sobre
12
n

A
t
V

= 0 sobre
12
(5)
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 225
Esta formulaci on es una de las mas utilizadas, puesto que para medios lineales, isotropos y
homogeneos, es constante. En consecuencia, en la primera y tercera ecuacion de (5), el
termino: AAAAAAA AAAAAAA AAAAAAA AAAAA AA se reduce a: AAAAAAA AAAAA AA que es facilmente implementable en un
codigo de elementos nitos. Esta simplicaci on se utiliza en la resoluci on de alguno de los
problemas presentados en los ejemplos numericos. Es mas, en problemas magnetostaticos
bidimensionales las ecuaciones de Maxwell (1) se reducen a
6,5
HHHHHHH HHHHH HH = JJJJJJJ JJJJJ JJ
imp
(6)
BBBBBBB BBBBB BB = 0 (7)
con JJJJJJJ JJJJJ JJ
imp
solo en la direcci on del eje z. Ademas, no existe potencial electrico y el potencial
vector magnetico es AAAAAAA AAAAA AA = (0, 0, A
z
) constante respecto al tiempo. Entonces, de (6) se deduce
que AAAAAAA AAAAA AA verica la ecuaci on de Poisson: ((0, 0, A
z
)) = (0, 0, J
imp
z
). Por consiguiente,
despues de imponer la condicion de contraste (gauge de Coulomb) AAAAAAA AAAAA AA = 0, necesaria para
asegurar unicidad de la solucion, el sistema de ecuaciones anterior se reduce a
6
A
z
= J
imp
z
en (8)
junto con las condiciones de contorno
A
z
= 0 sobre
B
A
z
n
= 0 sobre
H
A
z1
= A
z2
sobre
12
1

1
A
z1
n
=
1

2
A
z2
n
sobre
12
(9)
Es importante resaltar que una vez se ha hallado el potencial vector magnetico, la densidad
de ujo magnetico se obtiene mediante diferenciaci on numerica de acuerdo con (3). Ademas,
en problemas no lineales bidimensionales, = (BBBBBBB BBBBB BB), y por consiguiente la ecuaci on (8) es
A
z
= (BBBBBBB BBBBB BB)J
imp
z
.
Notese que las expresiones (8) y (9) son v alidas para el problema bidimensional plano.
En el caso de analizar un problema bidimensional axisimetrico, las ecuaciones anteriores se
pueden escribir como (ver detalles en la referencia
6
).

= J
imp

en (10)
junto con las condiciones de contorno
A

= 0 sobre
B
A

n
= 0 sobre
H
A
1
= A
2
sobre
12
1

1
A
1
n
=
1

2
A
2
n
sobre
12
(11)
La segunda formulaci on se denomina AAAAAAA AAAAA AA, V , . La diferencia esencial respecto a la
formulacion anterior es la utilizacion en la regi on no conductora de potenciales escalares
226 J. Sarrate y R. Clariso
magneticos en lugar del potencial vector magnetico, ver detalles en las referencias
7,3
. La
idea basica es descomponer la zona no conductora,
2
, en dos regiones: una que contiene
las fuentes de corriente,

2
, y la otra compuesta por el resto del medio no conductor,

2
,
separadas por una intercara,

(Figura 4). Esta formulaci on se utilizara en el cuarto

ejemplo numerico.
En la regi on no conductora donde existen fuentes de corriente se dene el potencial
escalar reducido como
HHHHHHH HHHHH HH = HHHHHHH HHHHH HH
imp
(12)
donde HHHHHHH HHHHH HH
imp
es la intensidad de campo magnetico debido a las fuentes de corriente,
HHHHHHH HHHHH HH
imp
= JJJJJJJ JJJJJ JJ
imp
y es por tanto conocido. En el resto de la zona no conductora se dene el
potencial escalar total magnetico como
HHHHHHH HHHHH HH = (13)
Figura 4. Dominio y condiciones de contorno para un problema de corrientes de Foucault
compuesto por dos medios donde la zona no conductora se ha dividido en

2
y

2
La descripcion del medio no conductor mediante unicamente el potencial escalar reducido
conduce en algunas ocasiones a errores de redondeo importantes. Puesto que HHHHHHH HHHHH HH = HHHHHHH HHHHH HH
imp

y el termino es del mismo orden que HHHHHHH HHHHH HH

imp
(ref.
7
), entonces los problemas de cancelaci on
pueden conducir a valores de HHHHHHH HHHHH HH muy inexactos. Ademas, para materiales con valores altos
de , y de acuerdo con BBBBBBB BBBBB BB = HHHHHHH HHHHH HH, los errores de cancelacion pueden verse dr asticamente
amplicados, conduciendo a valores de la densidad de ujo magnetico completamente falsos.
El sistema de ecuaciones equivalente a (1) que se debe resolver en funci on de los nuevos
potenciales es
7,3
AAAAAAA AAAAA AA AAAAAAA AAAAA AA +
A AAAAAA AAAAA AA
t
+V = 0 000000 00000 00 en
1

A
t
V

= 0 en
1
() = 0 en

2
((HHHHHHH HHHHH HH
imp
)) = 0 en

2
n () = 0 sobre
B
(14)
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 227
= 0 sobre
H
nnnnnnn nnnnn nn

2
+

2
(HHHHHHH HHHHH HH
imp
)) = 0 sobre

( + (HHHHHHH HHHHH HH
imp
)) nnnnnnn nnnnn nn

= 0 000000 00000 00 sobre

n
1
AAAAAAA AAAAA AA n
2

2
= 0 sobre
12

1
AAAAAAA AAAAA AA n
1
n
2
= 0 000000 00000 00 sobre
12
n AAAAAAA AAAAA AA = 0 sobre
12
n

A
t
V

= 0 sobre
12
Es importante resaltar que, en general, debido a la utilizaci on de y en lugar de AAAAAAA AAAAA AA,
en la regi on no conductora se ha reducido el n umero de inc ognitas de tres a una (n otese
que, por ejemplo, en problemas bidimensionales planos con las componentes del campo en
el plano de trabajo, el potencial vectorial se reduce a su componente perpendicular, y por
consiguiente el ahorro anteriormente mencionado no es totalmente cierto). En la resoluci on
del problema discreto, esto se traduce en un ahorro considerable del coste computacional.
Una simplicacion de esta formulaci on consiste en utilizar unicamente el potencial escalar
total magnetico en la zona no conductora. Esto da lugar a una formulaci on denominada
AAAAAAA AAAAA AA, V (ref.
3
).
Cuando estas formulaciones se utilizan sobre dominios con regiones conductoras no
simplemente conexas se generan soluciones esp ureas. Esto es debido a la imposibilidad de
demostrar la unicidad de la soluci on bajo dichas condiciones. Por este motivo se introduce
una tercera formulaci on denominada: AAAAAAA AAAAA AA, V AAAAAAA AAAAAAA, .

Esta es mas general y se puede entender
como una combinacion de las formulaciones AAAAAAA AAAAA AA, V AAAAAAA AAAAA AA y AAAAAAA AAAAA AA, V .
En la Figura 5 se presenta un esquema del dominio y sus contornos. El dominio total
se subdivide en tres regiones: la primera es una regi on conductora no simplemente conexa,

1
; la segunda es una region no conductora
3
, tal que su uni on con
1
genera una regi on
simplemente conexa; y la tercera,
2
, esta formada por el resto de la regi on no conductora.
Las intercaras entre las tres regiones se denotan por:
12
,
13
y
23
, mientras que los
contornos exteriores de las regiones no conductoras se subdividen en dos partes a semejanza
de la divisi on realizada en la primera formulaci on:
B2
,
H2
,
B3
y
H3
. De acuerdo con
la exposici on presentada en la referencia
3
, el potencial vector magnetico se utiliza en
3
,
mientras que en
2
se utiliza el potencial escalar total magnetico. El sistema de ecuaciones
que conduce a la unicidad de los potenciales es
3
Figura 5. Dominio y condiciones de contorno para un problema de corrientes de Foucault
compuesto por tres medios
228 J. Sarrate y R. Clariso
AAAAAAA AAAAA AA AAAAAAA AAAAA AA+
A
t
+ V = 0000000 00000 00 en
1

A
t
V

= 0 en
1
AAAAAAA AAAAA AA AAAAAAA AAAAA AA = 0000000 00000 00 en
3
= 0 en
2
nnnnnnn nnnnn nn AAAAAAA AAAAA AA = 0000000 00000 00 sobre
B3
AAAAAAA AAAAA AA = 0 sobre
B3
AAAAAAA AAAAA AA n = 0000000 00000 00 sobre
H3
nnnnnnn nnnnn nn AAAAAAA AAAAA AA = 0 sobre
H3
nnnnnnn nnnnn nn = 0 sobre
B2
= 0 sobre
H2
AAAAAAA AAAAA AA
1
= AAAAAAA AAAAA AA
3
sobre
13

1
AAAAAAA AAAAA AA
1
nnnnnnn nnnnn nn
1
+
3
AAAAAAA AAAAA AA
3
nnnnnnn nnnnn nn
3
= 0000000 00000 00 sobre
13

1
AAAAAAA AAAAA AA
1

3
AAAAAAA AAAAA AA
3
= 0 sobre
13
nnnnnnn nnnnn nn

A
t
V

= 0 sobre
12
y
13
nnnnnnn nnnnn nn
13
AAAAAAA AAAAA AAn
2

2
= 0 sobre
12
y
23
AAAAAAA AAAAA AA nnnnnnn nnnnn nn
13
nnnnnnn nnnnn nn
2
= 0000000 00000 00 sobre
12
y
23
n AAAAAAA AAAAA AA = 0 sobre
12
y
23
(15)
Existen varias alternativas de resolver mediante el MEF los sistemas de ecuaciones (5),
(14) o (15). Una de las m as utilizadas es la formulaci on de Galerkin, en la que coinciden las
funciones de peso y las funciones de forma. Puesto que las expresiones (5) y (14) pueden
ser vistas como un caso particular de la formulaci on (15), a continuaci on se presentan las
ecuaciones correspondientes a la discretizacion mediante el MEF de la tercera formulaci on
(15)
3
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 229

1
+
3
(NNNNNNN NNNNN NN
i
AAAAAAA AAAAA AA + NNNNNNN NNNNN NN
i
AAAAAAA AAAAA AA + NNNNNNN NNNNN NN
i

A AAAAAA AAAAA AA
t
+NNNNNNN NNNNN NN
i
V )d

12
+
23
NNNNNNN NNNNN NN
i
( nnnnnnn nnnnn nn
2
)d = 0

1
(NNNNNNN NNNNN NN
i

A
t
+ NNNNNNN NNNNN NN
i
V )d = 0

2
(NNNNNNN NNNNN NN
i
)d

12
+
23
NNNNNNN NNNNN NN
i
(AAAAAAA AAAAA AA nnnnnnn nnnnn nn
13
)d = 0
(16)
donde i es el ndice asociado a los nodos de la malla y N
i
son las funciones de forma.
Las formulaciones basadas en potenciales permiten modelizar adecuadamente materiales
no homogeneos as como discontinuidades materiales. Adem as, los potenciales vectoriales
presentan la gran ventaja de vericar autom aticamente la ecuacion de la divergencia. Sin
embargo, presentan una perdida notable de precisi on en el calculo de las intensidades
de campo electrico y magnetico debido a la diferenciaci on numerica. Ademas, es preciso
imponer para cada formulaci on el gauge adecuado (e.g. AAAAAAA AAAAA AA = 0 para casos cuasiest aticos;
AAAAAAA AAAAA AA+
V
t
= 0 para casos acoplados o de propagaci on; ...).
Basadas en la resolucion de las ecuaciones de Maxwell
Desde el punto de vista computacional, la resoluci on numerica de las ecuaciones de
Maxwell (1) puede presentar soluciones esp ureas. La caracterstica basica de estas solu-
ciones es que no verican las leyes de Gauss para el campo electrico y magnetico, ecuaciones
que contienen el operador divergencia en (1). En general, se ha considerado que la causa de
estas discrepancias estaba en el metodo numerico utilizado. Sin embargo, la referencia
11
de-
muestra que la aparici on de estas soluciones cticias radica en un planteamiento err oneo del
problema: tradicionalmente, las leyes de Gauss han sido vistas como ecuaciones auxiliares.
Esto se debe a dos motivos. Por una parte, el sistema de ecuaciones (1) parece sobredetermi-
nado: para el caso tridimensional el sistema presenta 8 ecuaciones y 6 inc ognitas, mientras
que para el caso bidimensional el sistema presenta 6 ecuaciones y 4 incognitas. Por otra
parte, si se aplica el operador divergencia a la ley de Faraday y a la ley de Maxwell-Ampere,
se obtiene que las leyes de Gauss para el campo electrico y magnetico se verican siem-
pre si estas se verican inicialmente
12
. Asimismo, en la referencia
11
se demuestra que, si
se incluyen las leyes de Gauss, entonces las ecuaciones de Maxwell son equivalentes a un
sistema elptico bien determinado en el espacio. Siendo necesario, por consiguiente, tenerlas
en cuenta.
En la ultima decada, la utilizacion de este tipo de formulaciones en la resoluci on de prob-
lemas electromagneticos ha experimentado un gran auge. Su aplicaci on va desde el calculo
de campos estaticos hasta problemas transitorios, pasando por problemas arm onicos. En
general, existen tres formulaciones para resolver (1) sin utilizar potenciales (por consigu-
iente, sin tener que diferenciar numericamente para obtener EEEEEEE EEEEE EE y BBBBBBB BBBBB BB). Con el objetivo de
introducir las caractersticas basicas de cada una de ellas, considerese un dominio , aco-
tado, simplemente conexo y convexo (Figura 6), donde su contorno se divide en dos partes,

1
y
2
, tales que

1

2
=

1

2
=
230 J. Sarrate y R. Clariso
Figura 6. Dominio y condiciones de contorno para un problema de corrientes de Foucault
compuesto por un unico medio
En cada una de ellas, se imponen las siguientes condiciones de contorno (como en la segunda
seccion, por simplicidad y sin perdida de generalidad estas se supondr an homogeneas)
nnnnnnn nnnnn nn EEEEEEE EEEEE EE = 0000000 00000 00 sobre
1
nnnnnnn nnnnn nn (HHHHHHH HHHHH HH) = 0 sobre
1
nnnnnnn nnnnn nn HHHHHHH HHHHH HH = 0000000 00000 00 sobre
2
nnnnnnn nnnnn nn (EEEEEEE EEEEE EE) = 0 sobre
2
(17)
La primera formulaci on se basa en la resolucion de ecuaciones de primer orden.

Esta
se basa en la resolucion directa del sistema (1) junto con las condiciones de contorno
(17). B asicamente, existen dos razones por las que esta formulacion es, en ocasiones,
desaconsejable. La primera es la dicultad que tienen los metodos numericos convencionales
en tratar operadores no autoadjuntos. La segunda es que las intensidades de campo electrico
y magnetico aparecen acopladas en (1), aumentando considerablemente la complejidad de
los calculos.
La segunda formulaci on se basa en la resolucion de ecuaciones de segundo orden.

Esta
se deriva del sistema (1) aplicando el operador rotacional, obteniendo de esta forma un
sistema de ecuaciones de segundo orden que desde el punto de vista numerico es mas f acil
de abordar. En la referencia
11
se demuestra que la resolucion de (1) es equivalente a resolver
(E) +

t

(E)
t
+ E

= K
eqv

J
imp
t
en
(E) =
imp
en
n E = 0000000 00000 00 sobre
1
n (E) = 0 sobre
2
n (E) = n K
eqv
sobre
2
(H) +

t
+

(H)
t
= J
imp

t
+

K
eqv
en
(H) = 0 en
n H = 0000000 00000 00 sobre
2
n (H) = 0 sobre
1
n (H) = n J
imp
sobre
1
(18)
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 231
Como puede observarse en (18), las intensidades de campo electrico y magnetico aparecen
desacopladas simplicando considerablemente los calculos. Sin embargo, de acuerdo con la
discusi on presentada en la referencia
11
, el sistema (18) puede presentar soluciones esp ureas
si no se consideran todas sus ecuaciones. Esto es debido a que generalmente las ecuaciones
en divergencia, segunda y septima ecuaciones de (18), as como las que se derivan de ellas,
quinta y decima ecuaciones de (18), eran consideradas ecuaciones complementarias y no se
tenan en cuenta en la resoluci on del problema. Bajo estas aproximaciones, el sistema(18)
admite mas soluciones que el sistema original (1).
La tercera formulacion se basa en la resolucion de la ecuaci on de Helmholtz. En la
referencia
11
se demuestra que el sistema (18) es equivalente a
EEEEEEE EEEEE EE +

t

(EEEEEEE EEEEE EE)

t
+EEEEEEE EEEEE EE

= KKKKKKK KKKKK KK
eqv

J JJJJJJ JJJJJ JJ
imp
t

imp
en
(EEEEEEE EEEEE EE) =
imp
sobre
1
nnnnnnn nnnnn nn E = 0000000 00000 00 sobre
1
nnnnnnn nnnnn nn (EEEEEEE EEEEE EE) = 0 sobre
2
nnnnnnn nnnnn nn (EEEEEEE EEEEE EE) = nnnnnnn nnnnn nn KKKKKKK KKKKK KK
eqv
sobre
2
HHHHHHH HHHHH HH +

t
+

(HHHHHHH HHHHH HH)

t
= J
imp

t
+

KKKKKKK KKKKK KK
eqv
en
(H) = 0 sobre
2
n H = 0000000 00000 00 sobre
2
nnnnnnn nnnnn nn (HHHHHHH HHHHH HH) = 0 sobre
1
nnnnnnn nnnnn nn (HHHHHHH HHHHH HH) = nnnnnnn nnnnn nn JJJJJJJ JJJJJ JJ
imp
sobre
1
(19)
La ventajas del sistema (19) son que: 1. la ecuaci on de la divergencia no debe vericarse
en todo el dominio, sino que s olo debe cumplirse en una porcion del dominio y 2. como en
(18) las ecuaciones para EEEEEEE EEEEE EE y HHHHHHH HHHHH HH aparecen desacopladas.
Merece la pena indicar que los metodos basados en la resoluci on directa de las ecua-
ciones de Maxwell se emplean basicamente para aplicaciones donde es importante analizar
el acoplamiento electrico y magnetico, mientras que los metodos basados en potenciales
se utilizan primordialmente cuando es suciente simular la din amica mediante corrientes
inducidas.
Es importante resaltar que la resoluci on directa de las ecuaciones de Maxwell permite
obtener directamente los vectores de intensidad electrica y magnetica sin necesidad de
utilizar derivaci on numerica, obteniendose por lo tanto resultados m as exactos. Ademas, se
obtienen campos continuos a traves de los elementos. Tambien se debe resaltar que mediante
este tipo de formulaciones no ha sido necesario denir ning un gauge.
Finalmente, debe mencionarse que de forma similar a los problemas de uidos o s olidos
incompresibles, tambien se han desarrollado formulaciones mixtas para problemas electro-
magneticos
13,12
. A diferencia del trabajo presentado en la referencia
11
, estas formulaciones
requieren la utilizaci on de elementos de diferente orden para las inc ognitas.
FUERZAS MAGN

ETICAS
Conceptos preliminares
Desde el punto de vista mec anico, el ingeniero est a interesado en el calculo de dos tipos
de fuerzas. La primera es la fuerza local (tambien denominada en algunos ambitos densidad
232 J. Sarrate y R. Clariso
de fuerza). El c alculo de las fuerzas locales es b asico en el estudio y an alisis de sistemas
formados por cuerpos deformables. En estos casos, el objetivo es conocer el comportamiento
local del material. El segundo tipo de fuerzas es la fuerza total.

Esta es fundamental en
el an alisis del movimiento de un s olido rgido. En este caso, el objetivo es conocer el
comportamiento global del cuerpo. Pero incluso en este ultimo supuesto, el conocimiento
de la distribuci on de fuerzas locales puede proporcionar una informaci on muy importante
en el desarrollo y mejora del dise no nal.
Desde el punto de vista conceptual, ambos tipos de fuerzas se hallan relacionados medi-
ante la siguiente expresi on
FFFFFFF FFFFF FF
total
=

fffffff fffff ff
local
d
Es decir, la fuerza total se puede calcular como la suma de las fuerzas locales que act uan
sobre el cuerpo. Cuando las fuerzas locales act uan sobre la supercie de un medio, la forma
usual de realizar dicha suma es desacoplando la fuerza local en su componente normal y
tangencial
FFFFFFF FFFFF FF =

fffffff fffff ff
t
ttttttt ttttt tt +fffffff fffff ff
n
nnnnnnn nnnnn nn d
donde ttttttt ttttt tt y nnnnnnn nnnnn nn denotan los vectores normalizados tangencial y normal a la supercie.
Con el objetivo de presentar los metodos mas ampliamente utilizados en el calculo
de fuerzas electromagneticas, a continuacion se presentan tres formulaciones diferentes de
calcularla
27,20,22,5
.
En la primera formulaci on, la fuerza total electromagnetica FFFFFFF FFFFF FF
em
se expresa en funci on
de las cargas en movimiento, la intensidad de campo electrico EEEEEEE EEEEE EE y la densidad de ujo
magnetico BBBBBBB BBBBB BB. En particular, si s olo existe una carga puntual en movimiento q, dicha relaci on
esta determinada por la ley de Lorentz
FFFFFFF FFFFF FF
em
= FFFFFFF FFFFF FF
e
+FFFFFFF FFFFF FF
m
= qEEEEEEE EEEEE EE + q(vvvvvvv vvvvv vv BBBBBBB BBBBB BB) (20)
siendo FFFFFFF FFFFF FF
e
la fuerza electrica, FFFFFFF FFFFF FF
m
la fuerza magnetica y vvvvvvv vvvvv vv la velocidad de la partcula.
En conductores, donde en lugar de partculas existen densidades de corrientes, la expresi on
(20) se debe modicar adecuadamente. Entonces, para una intensidad de corriente I (Amp),
una densidad supercial de corriente kkkkkkk kkkkk kk (Amp/m) y una densidad vol umica de corriente JJJJJJJ JJJJJ JJ
(Amp/m
2
), la ecuaci on (20) se escribe respectivamente como
FFFFFFF FFFFF FF
m
=

C
Iddddddd ddddd ddL BBBBBBB BBBBB BB (21)
FFFFFFF FFFFF FF
m
=

kkkkkkk kkkkk kk BBBBBBB BBBBB BBd (22)

FFFFFFF FFFFF FF
m
=

JJJJJJJ JJJJJ JJ BBBBBBB BBBBB BBd (23)

En algunas aplicaciones es posible hallar expresiones m as simplicadas que las ecuaciones
(21), (22) y (23). As por ejemplo, considerando que: 1. para campos electricos estacionarios
la ley de Maxwell-Ampere en el vaco es BBBBBBB BBBBB BB =
0
JJJJJJJ JJJJJ JJ y 2. la identidad BBBBBBB BBBBB BB BBBBBBB BBBBB BB =
(
1
2
BBBBBBB BBBBB BB
2
) (BBBBBBB BBBBB BB )BBBBBBB BBBBB BB, el integrando de la ecuaci on (23) puede escribirse como
JJJJJJJ JJJJJ JJ BBBBBBB BBBBB BB = (
BBBBBBB BBBBB BB
2
2
0
) +
1

0
(BBBBBBB BBBBB BB )BBBBBBB BBBBB BB (24)
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 233
donde la cantidad p
m
=
BBBBBBB BBBBB BB
2
2
0
es la presi on magnetica. Integrando se obtiene
FFFFFFF FFFFF FF
m
=

JJJJJJJ JJJJJ JJ BBBBBBB BBBBB BBd =

(
BBBBBBB BBBBB BB
2
2
0
) +
1

0
(BBBBBBB BBBBB BB )BBBBBBB BBBBB BB d =

1
2
0
BBBBBBB BBBBB BB
2
nnnnnnn nnnnn nn +
1

0
BBBBBBB BBBBB BB(BBBBBBB BBBBB BB nnnnnnn nnnnn nn) d (25)
La importancia de la ecuaci on (25) radica en que ahora es posible calcular directamente la
fuerza total a partir del campo BBBBBBB BBBBB BB calculado previamente mediante el MEF.
Si en el an alisis, ademas de medios conductores se deben considerar materiales
magneticos, entonces la expresi on (23) es
FFFFFFF FFFFF FF
m
=

JJJJJJJ JJJJJ JJ BBBBBBB BBBBB BBd

1
2
HHHHHHH HHHHH HH
2
d (26)
que a semejanza de (25) se puede expresar en terminos de HHHHHHH HHHHH HH, calculado anteriormente
mediante el MEF.
En la segunda formulaci on, la fuerza total se expresa en funci on de la energa magnetica
U (ref.
27,22
). En particular, es posible expresar la energia magnetica como
U =

B
0
HHHHHHH HHHHH HHdB

d =

1
2
HHHHHHH HHHHH HH
2
d (27)
Entonces, la fuerza magnetica total en la direccion p es
FFFFFFF FFFFF FF
m,p
=

U
p

(28)
donde es una restriccion de ujo constante ( =

BBBBBBB BBBBB BB d = cte).

La tercera formulacion se basa en reemplazar los materiales magneticos por corrientes
magneticas o cargas magneticas equivalentes (Equivalent Magnetizing Currents, Equiva-
lent Magnetizing Charges, EMC) y es ampliamente utilizada en problemas con campos
estacionarios
20,22,16
. En magnetost atica, la densidad de ujo magnetico puede expresarse
como
BBBBBBB BBBBB BB =
0

r
HHHHHHH HHHHH HH =
0
(HHHHHHH HHHHH HH +MMMMMMM MMMMM MM)
donde MMMMMMM MMMMM MM es la densidad de magnetizaci on y
r
es la permeabilidad relativa del medio.
Entonces, la ley de Maxwell-Ampere, HHHHHHH HHHHH HH = JJJJJJJ JJJJJ JJ
imp
, puede escribirse como

0
BBBBBBB BBBBB BB = JJJJJJJ JJJJJ JJ
imp
+MMMMMMM MMMMM MM (29)
En la expresi on anterior, se dene la densidad de corriente magnetica equivalente como
JJJJJJJ JJJJJ JJ
M
= MMMMMMM MMMMM MM (30)
La importancia de esta nueva magnitud radica en que permite plantear el problema (29)
en el vaco: ahora existe una nueva intensidad de campo magnetico

H

H, relacionada con la
densidad de ujo magnetico BBBBBBB BBBBB BB, mediante la permeabilidad del vaco
0
. Es decir, la ley
constitutiva del material es: BBBBBBB BBBBB BB =
0

H

H y la expresi on (29) puede escribirse como

H = JJJJJJJ JJJJJ JJ
imp
+JJJJJJJ JJJJJ JJ
M
En consecuencia, la fuerza resultante puede calcularse mediante la ley de Biot-Savart
234 J. Sarrate y R. Clariso
F
m
=

(JJJJJJJ JJJJJ JJ
imp
+JJJJJJJ JJJJJ JJ
M
) BBBBBBB BBBBB BBd
En el caso de cargas magneticas, el calculo de la fuerza magnetica se realiza mediante
un desarrollo semejante al anterior deniendo la densidad de carga equivalente
Q
M
=
0
M (31)
Puesto que
0
M = (1
1

r
)B, tomando divergencia se obtiene
Q
M
=
0
M = (1
1

r
) B+(
1

r
) B
Entonces, considerando la ley de Gauss para el campo magnetico, B = 0, se obtiene que:
Q
M
= (
1

r
) BBBBBBB BBBBB BB. Finalmente, la fuerza magnetica se calcula mediante una analoga con la
ley de Coulomb para el campo electrico (FFFFFFF FFFFF FF
e
= qEEEEEEE EEEEE EE)
FFFFFFF FFFFF FF
m
= Q
M
HHHHHHH HHHHH HH = (
1

r
) BBBBBBB BBBBB BBHHHHHHH HHHHH HH
Metodos: clasicacion y descripcion
El c alculo numerico de las fuerzas magneticas generalmente se basa en el calculo previo
del campo magnetico mediante el MEF. La mayora de metodos numericos para calcu-
larlas se basan en las tres formulaciones presentadas en el subapartado anterior
23,20,17,24
.
Basicamente, existen cinco metodos de calcular las fuerzas magneticas.
El primero se denomina metodo del tensor de tensiones de Maxwell (FFFFFFF FFFFF FF
Ttm
m
). Este metodo
se basa en integrar la expresi on (25) sobre una supercie arbitraria que englobe el dominio
de interes . Notese que debido a la deducci on de (25), dicha supercie debe estar situada
en el vaco. El c alculo de dicha integral se realiza separando su componente tangencial:
f
t
= H
t
B
n
, de su componente normal: f
n
=
1
2
(
1

0
B
2
n

0
H
2
t
). Debido a su simplicidad y
similitud con el fenomeno fsico, es uno de los metodos mas utilizados. Sin embargo, es muy
susceptible a errores numericos y depende considerablemente de la supercie de integraci on
escogida
26
.
El segundo es el metodo del trabajo virtual de Coulomb (FFFFFFF FFFFF FF
Tvc
m
). Este metodo se basa en
la expresi on de la fuerza en funci on de la energa magnetica y puede ser implementado de
diversas maneras. La mas sencilla consiste en aproximar la expresi on (28) mediante diferen-
cias nitas. Otra posibilidad es la utilizaci on de las matrices jacobianas de la transformaci on
isoparametrica del MEF en el calculo de la integral
14,15
que aparece en (28). Por ultimo,
tambien es posible calcularla descomponiendo la fuerza en sus componentes tangencial y
normal, de forma semejante a la realizada en el metodo del tensor de tensiones de Maxwell.
En este caso
20
, se obtiene: fffffff fffff ff
t
= 0, fffffff fffff ff
n
=
1
2
[(
1

)BBBBBBB BBBBB BB
2
n
(
0
)HHHHHHH HHHHH HH
2
t
].
El tercer metodo es el denominado metodo de las corrientes de magnetizaci on (FFFFFFF FFFFF FF
Jm
m
).

Este se basa en substituir un medio magnetico por un medio no magnetico donde existe
una distribuci on de corrientes magneticas. La fuerza magnetica total se calcula a partir de
las componentes tangenciales y normales de la fuerza local como
20
: fffffff fffff ff
t
= (1
r
)BBBBBBB BBBBB BB
n
HHHHHHH HHHHH HH
t
y
fffffff fffff ff
n
=

0
2
(
2
r
1)HHHHHHH HHHHH HH
2
t
. El c alculo de la fuerza local mediante este metodo depende fuertemente
de la calidad de la malla y en general sus resultados no son muy precisos. Paralelamente,
en la referencia
19
se presenta una versi on diferente que permite mejorar los resultados.
El cuarto metodo es el metodo de las cargas magneticas (FFFFFFF FFFFF FF
Qm
m
). Consiste en substituir
un medio magnetico por un medio no magnetico donde existe una distribuci on de cargas
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 235
magneticas
20
. A semejanza del metodo anterior, la fuerza magnetica total se calcula a
partir de fffffff fffff ff
t
= (1
1

r
)BBBBBBB BBBBB BB
n
HHHHHHH HHHHH HH
t
y fffffff fffff ff
n
=
1
2
0
(1
1

2
r
)BBBBBBB BBBBB BB
2
n
. La precisi on de la fuerza total depende
esencialmente del calculo previo de BBBBBBB BBBBB BB
n
. Por lo tanto, en la implementaci on del metodo debe
realizarse especial enfasis en su calculo
20,22
.
Por ultimo, el quinto metodo es el metodo de la densidad de fuerza supercial (FFFFFFF FFFFF FF
Dfs
m
).

Este se deriva del metodo del tensor de tensiones de Maxwell, generando los mismos resul-
tados que el metodo del trabajo virtual de Coulomb
18,22
. En este caso, la densidad de fuerza
solo depende de la componente normal: fffffff fffff ff
t
= 0, fffffff fffff ff
n
=
1
2
[(
1

)BBBBBBB BBBBB BB
2
n
(
0
)HHHHHHH HHHHH HH
2
t
].
En la Tabla I se resumen las componentes normal y tangencial de la densidad de fuerza
para cada uno de los cinco metodos anteriores.
Metodos fffffff fffff ff
t
fffffff fffff ff
n
FFFFFFF FFFFF FF
Ttm
m
HHHHHHH HHHHH HH
t
BBBBBBB BBBBB BB
n
1
2
(
1

0
BBBBBBB BBBBB BB
2
n

0
HHHHHHH HHHHH HH
2
t
)
FFFFFFF FFFFF FF
Tvc
m
0
1
2
[(
1

)BBBBBBB BBBBB BB
2
n
(
0
)HHHHHHH HHHHH HH
2
t
]
FFFFFFF FFFFF FF
Jm
m
(1
r
)BBBBBBB BBBBB BB
n
HHHHHHH HHHHH HH
t

0
2
(
2
r
1)HHHHHHH HHHHH HH
2
t
FFFFFFF FFFFF FF
Qm
m
(1
1

r
)BBBBBBB BBBBB BB
n
HHHHHHH HHHHH HH
t
1
2
0
(1
1

2
r
)BBBBBBB BBBBB BB
2
n
FFFFFFF FFFFF FF
Dfs
m
0
1
2
[(
1

)BBBBBBB BBBBB BB
2
n
(
0
)HHHHHHH HHHHH HH
2
t
]
Tabla I. Densidades de fuerza superciales (fuerza local) para los diferentes metodos
Es importante resaltar que te oricamente todos estos metodos son equivalentes
20,22
, puesto
que se verica
FFFFFFF FFFFF FF
Jm
m
= FFFFFFF FFFFF FF
Ttm
m

r
FFFFFFF FFFFF FF

FFFFFFF FFFFF FF
Qm
m
= FFFFFFF FFFFF FF
Ttm
m

FFFFFFF FFFFF FF

FFFFFFF FFFFF FF
Dfs
m
= FFFFFFF FFFFF FF
Ttm
m
FFFFFFF FFFFF FF

FFFFFFF FFFFF FF
Dfs
m
= FFFFFFF FFFFF FF
Tvc
m
donde el termino FFFFFFF FFFFF FF

S
(HHHHHHH HHHHH HH
t
BBBBBBB BBBBB BB
n
)ttttttt ttttt tt +
1
2

BBBBBBB BBBBB BB
2
n
HHHHHHH HHHHH HH
2
t

nnnnnnn nnnnn nn es la integral sobre una supercie

interior arbitraria, y por consiguiente debera ser nula
22
. Sin embargo, en las aplicaciones
numericas la precisi on de los metodos es diferente porque estos se basan en una aproximaci on
previa, mediante el MEF, de BBBBBBB BBBBB BB
n
y HHHHHHH HHHHH HH
t
. Por consiguiente, los cinco metodos anteriores
proporcionan valores diferentes de la fuerza magnetica
21,25
.
EJEMPLOS NUM

ERICOS
En este apartado se presentan cuatro ejemplos de c alculos magnetostaticos basados en
algunas de las formulaciones expuestas anteriormente. En todos ellos, se han utilizado
diferentes formulaciones basadas en el c alculo de potenciales. Asimismo, todos los c alculos
se han realizado mediante el c odigo de elementos nitos orientado al objeto CASTEM2000.
En el primer y tercer ejemplo solo se ha estudiado un medio lineal. Por el contrario, en
el segundo y cuarto ejemplo se ha considerado tanto un material con una permeabilidad
constante diferente de
0
, como un material no lineal (hierro) cuya ley constitutiva se
presenta en la Figura 7.
236 J. Sarrate y R. Clariso
0 1 2 3 4 5 6
x 10
4
0
0.5
1
1.5
2
2.5
H (Amp/m)

B

(
T
)

Figura 7. Ley constitutiva no lineal para el hierro
En el primer ejemplo se analiza un dipolo magnetico. Este problema, extremadamente
sencillo, permitir a vericar la calidad de los resultados obtenidos. En el segundo ejemplo
se analiza un conductor de secci on cuadrada por el que circula una corriente prescrita
5
. Se
considera que el medio puede ser lineal o no lineal. Ambos supuestos se estudian utilizando
dos mallas diferentes a n de constatar la mejora de los resultados. Asimismo, se calcula la
densidad de fuerza magnetica en el interior del conductor mediante (23). El tercer ejemplo
consiste en una aplicaci on industrial. En el se analiza un transformador de dos bobinas
5
.
En el cuarto ejemplo se comparan los resultados obtenidos mediante el an alisis axisimetrico
de una bobina con los generados mediante un estudio tridimensional.
Ejemplo 1
Como se ha comentado anteriormente, el primer ejemplo consiste en el analisis de un
dipolo magnetico: dos conductores paralelos e innitamente largos que transportan una
corriente de 10
4
Amp. En el presente estudio se analizan dos casos. En el primero se supone
que las corrientes tienen sentidos opuestos, mientras que en el segundo se supone que tienen
el mismo sentido. El an alisis bidimensional de este problema se ha realizado de acuerdo con
la aproximaci on potencial denida por las ecuaciones (8) y (9). La Figura 8a presenta la
geometra del problema as como las condiciones de contorno. El dominio se ha discretizado
mediante una malla estructurada formada por 2542 elementos cuadrangulares y 2645 nodos
(Figura 8b). Puesto que este ejemplo posee solucion analtica, se analiza la exactitud de
los resultados numericos. En particular, se estudia el error del potencial magnetico y de la
densidad de ujo magnetico a lo largo de la lnea de union de los dos polos.
X
Y
a) b)
Figura 8. a) Dominio y condiciones de contorno para un dipolo magnetico; b) Malla de
c alculo
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 237
En la primera columna de la Figura 9 se presentan los resultados obtenidos considerando
que las direcciones de corriente en los conductores son opuestas. Como es sabido, existe
una interacci on destructiva entre los dos polos. En la Figura 9a se muestran las lneas
equipotenciales de A
z
. Como era de esperar, se obtienen lneas equipotenciales concentricas
que parten de los polos y se anulan a lo largo del eje x = 0, 5. En las Figuras 9b y 9c se
muestran para y = 0, 5 los valores numericos y analticos correspondientes al potencial A
z
y al m odulo de la densidad de ujo magnetico B.
-8.77E-03
-7.42E-03
-6.07E-03
-4.72E-03
-3.37E-03
-2.02E-03
-6.75E-04
6.75E-04
2.02E-03
3.37E-03
4.72E-03
6.07E-03
7.42E-03
8.77E-03
4.19E-04
1.26E-03
2.10E-03
2.93E-03
3.77E-03
4.61E-03
5.45E-03
6.29E-03
7.12E-03
7.96E-03
8.80E-03
9.64E-03
1.05E-02
1.13E-02
a) lneas equipotenciales para corrientes
opuestas
d) lneas equipotenciales para corrientes
paralelas
0.35 0.50 0.65
0.006
0.000
0.006
X (m)

A
z

(
T

m
)

0.35 0.50 0.65
0.007
0.009
0.011
X (m)

A
z

(
T

m
)

b) valor del potencial A
z
en y = 0, 5 e) valor del potencial A
z
en y = 0, 5
0.35 0.50 0.65
0.000
0.065
0.130
X (m)

B

(
T
)

0.35 0.50 0.65
0.000
0.063
0.130
X (m)

B

(
T
)

c) valor de B en y = 0, 5 f ) valor de B en y = 0, 5
Figura 9. Valores numericos (lnea continua) y valores analticos (lnea discontinua) obteni-
dos para el dipolo magnetico
238 J. Sarrate y R. Clariso
En la segunda columna de la Figura 9 se presentan los resultados obtenidos cuando
la corriente en ambos conductores tiene el mismo sentido. Como en el caso anterior los
resultados obtenidos coinciden con los esperados: lneas equipotenciales concentricas que
parten de los polos con interacci on constructiva. Por consiguiente, el ujo uye alrededor
de los polos (Figura 9d). En las Figuras 9e y 9f se muestran para y = 0, 5 los valores
numericos y analticos correspondientes al potencial A
z
y al m odulo de la densidad de ujo
magnetico B.
En la Figura 10 se analiza el error de la soluci on numerica respecto a la solucion analtica
para el caso de corrientes opuestas (primera columna) y para el caso de corrientes paralelas
(segunda columna). En las Figuras 10a y 10d se presenta el error absoluto del potencial A
z
para ambos casos. En las Figuras 10b y 10e se presenta el error absoluto del m odulo de la
densidad de ujo magnetico B para ambos casos.
0.35 0.50 0.65
0.0000
0.0001
0.0002
X (m)

E
r
r
o
r

A
b
s
o
l
u
t
o

A
z

(
T

m
)

0.35 0.50 0.65
0.0000
0.0001
0.0002
X (m)

E
r
r
o
r

A
b
s
o
l
u
t
o

A
z

(
T

m
)

a) error absoluto de A
z
para corrientes op-
uestas
d) error absoluto de A
z
para corrientes
paralelas
0.35 0.50 0.65
0.000
0.0005
0.0010
0.0015
X (m)

E
r
r
o
r

A
b
s
o
l
u
t
o

B

(
T
)

0.35 0.50 0.65
0.000
0.0005
0.0010
0.0015
X (m)

E
r
r
o
r

A
b
s
o
l
u
t
o

B

(
T
)

b) error absoluto de B para corrientes op-
uestas
e) error absoluto de B para corrientes par-
alelas
0.35 0.50 0.65
0.0
3.0
6.0
X (m)

E
r
r
o
r

R
e
l
a
t
i
v
o

B

(
%
)

0.35 0.50 0.65
0.0
3.0
6.0
X (m)

E
r
r
o
r

R
e
l
a
t
i
v
o

A
z

(
%
)

c) error relativo de B para corrientes op-
uestas
f ) error relativo de A
z
para corrientes par-
alelas
Figura 10. Comparaci on de los resultados obtenidos para el dipolo magnetico
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 239
Para el caso de corrientes opuestas, el potencial se anula para x = 0, 5 y se obtiene un error
relativo del potencial indeterminado. Consecuentemente en la Figura 10c s olo se muestra
el error relativo de B.

Este resulta ser siempre inferior al 6 %. Para el caso de corrientes
paralelas, la densidad de ujo B se anula en el punto medio (0,5; 0,5) y consecuentemente su
error relativo es indeterminado. Por lo tanto, la Figura 10f s olo muestra el error relativo del
potencial A
z
. Como puede observarse, se calcula el potencial A
z
con una precision inferior
al 1,5 %.
Debe resaltarse que los resultados numericos reproducen cuantitativamente los valores
analticos. Asimismo, en ambos casos estos resultados conservan la simetra del problema.
Debe resaltarse que los resultados numericos reproducen cuantitativamente los valores
analticos. Asimismo, en ambos casos estos resultados conservan la simetra del problema.
Ejemplo 2
En el segundo ejemplo se analiza un conductor de secci on cuadrada por el que circula
una densidad de corriente JJJJJJJ JJJJJ JJ
imp
= (0, 0, J
imp
z
) = (0, 0, 10
4
) Amp/m
2
. El an alisis se realiza
tanto para un conductor lineal (
r
= 10) como para un conductor no lineal. N otese que la
no linealidad del material est a en su permeabilidad (su ley constitutiva es la presentada en
la Figura 7). Como en el ejemplo anterior, mediante las ecuaciones (8) y (9) se ha realizado
el an alisis bidimensional de este problema. Debido a la simetra del problema s olo se ha
estudiado un cuarto del problema. En la Figura 11 se presenta la geometra del modelo y
las condiciones de contorno. El problema se ha discretizado utilizando dos mallas diferentes
a n de observar la mejora obtenida en la soluci on.
Figura 11. Dominio y condiciones de contorno para un conductor de seccion cuadrada
En la Figura 12 se presentan las lneas equipotenciales obtenidas utilizando ambas mallas
tanto para un conductor lineal como no lineal. En ambos casos, al aumentar la discretizaci on
del dominio se captura mejor el gradiente del potencial. En particular, para el problema
lineal los resultados concuerdan con los presentados en la referencia
5
y la mayor parte del
gradiente del potencial magnetico se concentra en el interior del conductor.
Los vectores de densidad de fuerza en el interior del conductor, calculados de acuerdo
con (24) y utilizando la malla na, se presentan superpuestos a la lneas de equipotencial
magnetico en la Figura 13. Como puede observarse, el vector de densidad de fuerza es
ortogonal a las lneas de equipotencial magnetico y su sentido esta determinado por el signo
del potencial: puesto que en este ejemplo J
imp
z
> 0, el vector densidad de fuerza esta
orientado hacia el centro del conductor.
240 J. Sarrate y R. Clariso
4.49E-08
1.35E-07
2.25E-07
3.14E-07
4.04E-07
4.94E-07
5.84E-07
6.74E-07
7.64E-07
8.54E-07
9.43E-07
1.03E-06
1.12E-06
1.21E-06
4.47E-08
1.34E-07
2.24E-07
3.13E-07
4.03E-07
4.92E-07
5.82E-07
6.71E-07
7.61E-07
8.50E-07
9.40E-07
1.03E-06
1.12E-06
1.21E-06
a) material lineal con malla grosera b) material lineal con malla na
1.33E-05
4.00E-05
6.67E-05
9.34E-05
1.20E-04
1.47E-04
1.73E-04
2.00E-04
2.27E-04
2.53E-04
2.80E-04
3.07E-04
3.33E-04
3.60E-04
1.33E-05
3.98E-05
6.63E-05
9.28E-05
1.19E-04
1.46E-04
1.72E-04
1.99E-04
2.25E-04
2.52E-04
2.78E-04
3.05E-04
3.31E-04
3.58E-04
c) material no lineal con malla grosera d) material no lineal con malla na
Figura 12. Lineas equipotenciales correspondientes al conductor de seccion cuadrada del
segundo ejemplo
2.96E-07
3.67E-07
4.38E-07
5.09E-07
5.79E-07
6.50E-07
7.21E-07
7.92E-07
8.63E-07
9.34E-07
1.00E-06
1.08E-06
1.15E-06
1.22E-06
Figura 13. Densidad de fuerza y lneas equipotenciales en el conductor de seccion cuadrada
del segundo ejemplo
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 241
Ejemplo 3
El tercer ejemplo consiste en una aplicacion industrial
5
. En particular se analiza el
ujo magnetico generado en un transformador de doble bobina de 3 4 10 cm. Puesto
que la longitud del transformador es mayor que el resto de medidas, es posible estudiar su
comportamiento mediante un analisis bidimensional. La Figura 14 presenta una secci on del
mismo. El an alisis de este problema se ha realizado mediante las ecuaciones (8) y (9).
Figura 14. Seccion del transformador correspondiente al tercer ejemplo
En este ejemplo solo se considera un acero de permeabilidad relativa
r
= 100. Por la
bobina primaria circula una densidad de corriente de 100 Amp/cm
2
en una seccion de 0,6
cm
2
. Para que el ujo magnetico este orientado en el sentido de las manecillas del reloj
dentro del acero, la corriente electrica debe entrar en el papel (direccion de Z negativa)
en el lado derecho del transformador y salir del papel (direcci on de Z positiva) en el lado
izquierdo del mismo. La densidad de corriente en el lado derecho de la bobina primaria es
de 100 Amp/cm
2
, mientras que en el lado izquierdo es de 100 Amp/cm
2
. Puesto que el
problema es simetrico respecto el eje de abcisas, solo se analiza la mitad del problema. En
la Figura 15 se presentan las condiciones de contorno.
Figura 15. Dominio computacional y condiciones de contorno para el transformador presen-
tado en el tercer ejemplo
El problema se ha discretizado utilizando una malla estructurada formada por 3000
elementos cuadrangulares y 3131 nodos. En la Figura 16 se presenta la densidad de ujo
magnetico. Estos resultados concuerdan con los presentados por la referencia
5
.
242 J. Sarrate y R. Clariso
-3.12E-04
-2.87E-04
-2.61E-04
-2.36E-04
-2.10E-04
-1.85E-04
-1.59E-04
-1.34E-04
-1.08E-04
-8.26E-05
-5.71E-05
-3.16E-05
-6.11E-06
1.94E-05
Figura 16. Lneas de densidad de ujo magnetico en el interior del transformador
Ejemplo 4
En el cuarto ejemplo se analiza el campo magnetico dentro y alrededor de un cilindro
de hierro con una bobina en su interior por la que circula una densidad de corriente de
800 Amp/mm
2
. Se considera tanto el problema de un material lineal (
r
= 2900) como el
de material no lineal (Figura 7). La Figura 17 presenta un esquema del modelo.
Figura 17. Modelo geometrico correspondiente a la bobina presentada en el cuarto ejemplo
Puesto que el problema es axisimetrico, en este ejemplo se comparan los resultados
obtenidos mediante un c alculo bidimensional axisimetrico en el plano X Z (coincidiendo
la coordenada radial con el eje X) con los obtenidos mediante un c alculo tridimensional.
El an alisis bidimensional axisimetrico de este problema se ha realizado de acuerdo con la
aproximaci on potencial denida por las ecuaciones (10) y (11). Las condiciones de contorno
para el problema bidimensional axisimetrico se presentan en la Figura 18. Se ha impuesto
potencial nulo en las fronteras exteriores (lejos del conductor), condiciones de Neumann
homogeneas en la lnea de simetra y potencial nulo en el eje de simetra, puesto que en el
debe aparecer una lnea de ujo y esta ha de ser paralela a las lneas de equipotencial (como
consecuencia directa de imponer simetra de revoluci on). El problema se ha discretizado
utilizando una malla estructurada compuesta por 1216 cuadril ateros y 1287 nodos.
En la Figura 19a se presentan las lneas de equipotencial para un material lineal y en la
Figura 19b para un material no lineal. Puesto que la unica diferencia entre ambos casos es
el tipo de hierro, la diferencia entre ambas guras se centra en el hierro y no en la regi on
interna de aire o en la bobina.
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 243
Figura 18. Condiciones de contorno para el analisis axisimetrico de la bobina
5.11E06
3.58E05
6.64E05
9.71E05
1.28E04
1.58E04
1.89E04
2.20E04
2.50E04
2.81E04
3.12E04
3.42E04
3.73E04
4.04E04
4.34E04
4.65E04
4.96E04
5.26E04
5.57E04
5.88E04
6.18E04
6.49E04
a)
4.96E06
3.47E05
6.45E05
9.43E05
1.24E04
1.54E04
1.84E04
2.13E04
2.43E04
2.73E04
3.03E04
3.32E04
3.62E04
3.92E04
4.22E04
4.52E04
4.81E04
5.11E04
5.41E04
5.71E04
6.00E04
6.30E04
b)
Figura 19. An alisis axisimetrico de la bobina. Lneas equipotenciales en la zona hierroaire
bobina: a) material lineal; b) material no lineal
244 J. Sarrate y R. Clariso
Figura 20. Denicion de los dominios y condiciones de contorno para el analisis tridimen-
sional de la bobina
7.75E+03
7.38E+03
7.01E+03
6.65E+03
6.28E+03
5.92E+03
5.55E+03
5.18E+03
4.82E+03
4.45E+03
4.09E+03
3.72E+03
3.35E+03
2.99E+03
2.62E+03
2.26E+03
1.89E+03
1.52E+03
1.16E+03
7.93E+02
4.27E+02
61.
a)
7.74E+03
7.37E+03
7.01E+03
6.64E+03
6.28E+03
5.91E+03
5.54E+03
5.18E+03
4.81E+03
4.45E+03
4.08E+03
3.72E+03
3.35E+03
2.99E+03
2.62E+03
2.25E+03
1.89E+03
1.52E+03
1.16E+03
7.92E+02
4.26E+02
61.
b)
Figura 21. An alisis tridimensional de la bobina. Lneas equipotenciales en la zona

2
:
a) material lineal; b) material no lineal
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 245
Los calculos tridimensionales se han realizado en una regi on formada por un corte de 20

utilizando la formulaci on - presentada en las ecuaciones(14) de la tercera secci on. Puesto

que se trata de un problema magnetost atico, no existe la zona
1
. La regi on

2
contiene
la bobina y parte de la zona de aire interior. La regi on

2
esta formada por el resto del
dominio (Figura 20).
La malla est a formada por 2699 elementos y 5452 nodos. La regi on donde se dene el
potencial (

2
) esta formada por 270 elementos y 616 nodos. En la Figura 21a y 21b se
presenta la distribuci on del potencial en el dominio

2
para el modelo lineal y no lineal
respectivamente. Debido a que esta zona no incluye el hierro ambos resultados son muy
similares.
En la Figura 22a y 22b se presenta el potencial sobre el hierro. Por consiguiente, en
esta zona ambas distribuciones dieren substancialmente.
0.0
1.2
2.5
3.7
5.0
6.2
7.4
8.7
9.9
11.
12.
14.
15.
16.
17.
19.
20.
21.
22.
24.
25.
26.
a)
0.0
1.2
2.5
3.7
5.0
6.2
7.4
8.7
9.9
11.
12.
14.
15.
16.
17.
19.
20.
21.
22.
24.
25.
26.
b)
Figura 22. An alisis tridimensional de la bobina. Lneas equipotenciales en el hierro:
a) material lineal; b) material no lineal
246 J. Sarrate y R. Clariso
En la Figura 23a se muestra, para el an alisis bidimensional axisimetrico con material
lineal, la densidad de ujo magnetico representado sobre las lneas de equipotencial A

. En la
Figura 23b se muestra, para el an alisis tridimensional con material lineal, la densidad de ujo
magnetico representado sobre las lneas de equipotencial y , en la seccion y = 0. Como
puede observarse, los valores y la distribuci on de los potenciales son totalmente distintos,
puesto que su denici on es diferente (ecuaci on (3) para el potencial AAAAAAA AAAAA AA y ecuaciones (12)y (13)
para los potenciales y , respectivamente). Sin embargo, el comportamiento de la densidad
de ujo magnetico coincide satisfactoriamente tanto cualitativa como cuantitativamente.
a)
b)
Figura 23. Densidad de ujo magnetico y lneas equipotenciales en la zona hierroaire
bobina: a) c alculo 2D lineal utilizando potencial A

; b) c alculo 3D lineal
utilizando potenciales y
CONCLUSIONES
En este trabajo se ha presentado una clasicaci on de las formulaciones m as relevantes
para la resoluci on de problemas en electromagnetismo computacional.

Estas se han clasi-
cado en dos grandes grupos: 1. las formulaciones que utilizan potenciales y 2. las formu-
laciones que resuelven directamente las ecuaciones de Maxwell. Asimismo, se ha realizado
una presentacion de las caractersticas basicas de cada una de ellas, discutiendo sus ventajas
y limitaciones. Tambien se han presentado tres formulaciones diferentes para el c alculo de
fuerzas electromagneticas. A la vista de estos resultados, se han analizado y clasicado los
metodos mas utilizados en el calculo de dichas fuerzas. Finalmente se han presentado cuatro
El Metodo de los Elementos Finitos en problemas electromagneticos: planteamiento y aplicaciones 247
ejemplos numericos que incluyen problemas lineales y no lineales, tanto en dos como en tres
dimensiones. En su resoluci on mediante el MEF se han utilizado dos de las formulaciones
presentadas anteriormente que se basan en diferentes potenciales magneticos. Los resultados
obtenidos concuerdan tanto con la soluci on analtica como con las soluciones halladas en
estudios anteriores.
AGRADECIMIENTOS
Este trabajo ha sido parcialmente nanciado por el proyecto TAP98-0421 de la CICYT.
REFERENCIAS
1 J. Pallares, X. Granados, R. Bosch y X. Obradors, 1 kW bearingless superconducting motor
with axial eld exitation, 4th. European Conference on Applied Superconductivity, Sitges, Spain,
(1999).
2 H. Hashizume, T. Sugiura, K. Miya y S. Toda, Numerical analysis of electromagnetic phenomena
in superconductors, IEEE Transactions on Magnetics, Vol. 25, pp. 31453159, (1989).
3 O. Br o y K. Preis, On the use of the magnetic vector potential in the nite element analysis
of three-dimensional eddy currents, IEEE Transactions on Magnetics, Vol. 25, pp. 31453159,
(1989).
4 J. Donea, S. Giuliani y A. Philippe, Finite elements in the solution of electromagnetic induction
Problems, International Journal for Numerical Methods in Engineering, Vol. 8, pp. 359367,
(1974).
5 S. Hoole y P. Hoole, A modern short course in engineering electromagnetics, Oxford University
Press, New York, (1996).
6 J. Jin, The nite element method in electromagnetics, John Wiley and Sons, New York, (1993).
7 J. Simkin y C. W. Trowbridge, On the use of the total scalar potential in the numerical solution of
eld problems in electromagnetics, International Journal for Numerical Methods in Engineering,
Vol. 14, pp. 423440, (1979).
8 L. Hui L. y S. Saigal, Mapped innite elements for 3-D vector potencial problems, International
Journal for Numerical Methods in Engineering, Vol. 37, pp. 343356, (1994).
9 P.P. Silvester y R.L. Ferrari, Finite elements for electrical engineers, Cambridge University
Press, Cambridge, (1983).
10 C.W. Steele, Numerical Computation of Electric and Magnetic Fields, Van Nostrand Reinhold
Cop., New York, (1987).
11 B.N. Jiang, J. Wu y L.A. Povinelli, The origin of spurious solutions in computational electro-
magnetics, Journal of Computational Physics, Vol. 125, pp. 104123, (1996).
12 F. Assous, P. Degond, E. Heintze, P. A. Raviart y J. Seger, On a nite-element method for
solving the three-dimensional Maxwell equations, Journal of Computational Physics,Vol. 109,
pp. 222237, (1993).
13 P. Monk, A comparison of three mixed methods for the timedependent Maxwells equations,
SIAM Journal on Scientic and Statistical Computing, Vol. 13, pp. 10971122, (1992).
248 J. Sarrate y R. Clariso
14 J.L. Coulomb, A methodology for the determination of global electromechanical quantities from
a nite element analysis and its application to the evaluation of magnetic forces, torques and
stiness, IEEE Transactions on Magnetics, Vol. 19, pp. 25142519, (1983).
15 J.L. Coulomb y G. Meunier, Finite element implementation of virtual work principle for mag-
netic or electric force and torque computation, IEEE Transactions on Magnetics, Vol. 20,
pp. 18941896, (1984).
16 T. Kabashima, A. Kawahara y T. Goto, Force calculation using magnetizing currents, IEEE
Transactions on Magnetics, Vol. 24, pp. 451454, (1988).
17 W. M uller, Comparison of dierent methods of force calculation, IEEE Transactions on Mag-
netics, Vol. 26, pp. 10581061, (1990).
18 G. Henneberger, P.K. Sattler, W. Hadrys y D. Shen, Procedure for the numerical computation
of mechanical vibrations in electrical machines, IEEE Transactions on Magnetics, Vol. 28, pp.
13511354, (1992).
19 G. Henneberger, P.K. Sattler y D. Shen, Nature of the equivalent magnetizing current for the
force calculation, IEEE Transactions on Magnetics, Vol. 28, pp. 10681071, (1992).
20 N. Sadowski, Y. Lef`evre, M. Lajoie-Mazenc y J.P.A. Bastos, Sur le calcul des forces
magnetiques, J. Physique III , Vol. 2, pp. 859870, (1992).
21 M. Bartsch y T. Weiland, 2D and 3D calculation of forces, IEEE Transactions on Magnetics,
Vol. 30, pp. 34673470, (1994).
22 Z. Ren, Comparison of dierent force calculation methods in 3D nite element modelling, IEEE
Transactions on Magnetics, Vol. 30, pp. 34713474, (1994).
23 M. Kasper y J. Franz, Highly accurate computation of eld quantities and forces by supercon-
vergence in nite elements, IEEE Transactions on Magnetics, Vol. 31, pp. 14241427, (1995).
24 K. Komeza, A. Pelikant, J. Tegopoulos y S. Wiak, Comparative computation of forces and
torques of electromagnetic devices by means of dierent formulae, IEEE Transactions on Mag-
netics, Vol. 30, pp. 34753477, (1994).
25 R. Mertens, U. Pahner, K. Hameyer y R. Belmans, Force calculation based on a local solution
of Laplaces equation, IEEE Transactions on Magnetics, Vol. 33, pp. 12161218, (1997).
26 S. Salon, S. Bhatia y D. Burow, Some aspects of torque calculations in electrical machines,
IEEE Transactions on Magnetics, Vol. 33, pp. 20182021, (1997).
27 J.R. Reitz, F.J. Milford y R.W. Christy, Fundamentos de la teora electromagnetica, Addison-
Wesley, Wilmington, (1986).