Calculation of eddy currents in moving structures by a sliding mesh finite element method

Description

Calculation of eddy currents in moving structures by a sliding mesh-nite element method A. Bua , Y.Maday , F. Rapetti ASCI - UPR 9029 CNRS, Paris Sud University, Building 506, 91403 Orsay cedex, France LGEP - UMR 8507 CNRS, Supelec University Plateau de Moulon, 91192 Gif-Sur-Yvette cedex, France Abstract { In this paper we present some theoreti- cal and numerical results concerning the simulation of eddy currents in non-stationary structures. Both 2D and 3D models are considered. The approximation is based on the sliding-mesh mortar method combined with node elements in 2D and edge elements in 3D. An implicit Euler scheme is used to discretize in time. Key words { eddy currents, nite element approxima- tion on non-matching grids, magnetic eld as primary variable, moving structures. I. Introduction The subject of our research activity is the analysis and development of simulation tools to eectively predict the induced current distribution in non-stationary geometries with sliding interfaces. We then study the ux of eddy currents in a domain composed of two rigid conductive parts, a rotating one (rotor) in sliding contact with a xed one (stator). Our system is a sample of induction motor.

• edge elements


• both


• dimensional eddy


• constraint


• lagrangian variables


• current density


• conforming domain


• method

Sujets

• Constraint
• Current density
• Method

Informations

• Publié par: pefav
• Nombre de lectures: 10
• Langue: English

Extrait

t > 0t = 0
rt
x
r
x θr
x xs srotor rotor
r-t
stator stator
Bu

,
Yada
in
t
form
y
conuration

Euler
,
b
F
v
Rap
curl
etti
has

follo

v
ASCI
t
-
,
UPR
reference
9029
r
CNRS
H
P
n
aris
tions
Sud
sparse
Univ
in
ersit
ted
y
to
,
set
Building
is
506,
1
91403
the
Orsa
ving
y
with
cedex
=
F
as
rance
;

t
LGEP
n
-
or
UMR
are
8507
sub
CNRS
is
Sup
The
elec
e
Univ
w
ersit
ergence
y
the
Plateau
and
de
continuous
Moulon
least
91192
mo
Gifurv

ette
in
cedex
n
F
b
rance
1
A
(
bstr
ving
act
time
{
t
In
w
this
.
pap
x
er
problem
w
i
e

presen
0
t
;
some

theoreti
H
cal
(
and
t
n
distortion
umerical
adopted
results
e
concerning
a
the
prop
sim
the
ulation
y
of
ositiv
eddy
pap
curren
tion
ts
con
in
Section
nontationary
on
structures
b
Both
umerical
2
w
D
m
and
stated
3
VI
D
conclusions
mo

dels
hitz
are
o
considered

The

appro
2
ximation
The
is
sub
based
W
on
@
the
ard
slidingesh

mortar
angle
metho

d
erator
com
of
bined
r
with
b
no
domain
de
the
elemen
is
ts
problem
in
rotor
2
=
D
k
and
op
edge
r
elemen
x
ts
r
in
The
3
mag
D
read
.
1
An
J
implicit
(1)
Euler
n
sc

heme
0)
is
0
used
conditions
to
H
discretize
)
in
;
time
t
Key
;
wor
2
ds
:
{
lo
eddy
partial
curren
b
ts
these
ite
exp
elemen
with
t
[7]
appro
on
xima
W
tion
an
on
in
nonatc
linear
hing
haracterized
grids
ving
magnetic
and
ld
deite
as
of
primary
is
v
in
ariable
I
mo
t
ving
uous
structures
form
I
I
Intr
concen
oduction
2
The
tion
sub
con
ject
and
of
are
our
section
researc
presen
h
discrete
activit
A
y
result
is
Section
the
,
analysis
dev
and
considerations
dev
I
elopmen
oblem
t
R
of
a
sim
ounded
ulation
t
to
op
ols
1
to
e
ectiv
part
ely
its
predict

,
the

induced
1
curren

t
axisym
distribution
strictly
in
2
nontationary
call
geometries
namely
with
1
sliding
in
in
to
terfaces
.
W
t
e
the
then
time
study
t
the
!
x
rotation
of
h
eddy

curren
angle
ts
)
in
its
a
ee
domain
d
comp
metho
osed
the
of
form
t
conuration
w
=
o
en
rigid
for
conductiv
at
e
0,
parts
has
a
angle
rotating
t
one
go
otor
the
in
to
sliding
r
con
particular
tact
t
with
t
a
(0))
ed
(0)
one
(
tator
t
Our
curren
system
with
is
ld
a
v
sample
t
of
curl
induction
H
motor

W
)
e

w
(
ork
)
in

spaceime
@
v
;
ariables
x
and
H
w

e
(3)
consider
o
b
osed
oth
r
the
(
2
;
D
n
transv
(
erse
)
electric
(4)
E
H
form
t
ulation
)
and
curl
the
x
3
^
D
[8])
complete
the
system
cal
of
and
the
reeshing
eddy
to
curren
e
t
All
problem
remedies
These
rather
for
ensiv
m
and
ulations
limita
are
ee
obtained
for
from
review
Maxw
the
ell
ject
equations
e
where
ose
the
implicit
displacemen
discretization
t
time
curren
al
ts
system
are
c
neglected
b
with
ha
resp
a
ect
symmetric
to
p
the
e
conducting
matrix
ones
outline
The
the
magnetic
er
ld
the
H
wing
is
Sec
c
I
hosen
w
as
presen
primary
the
v
tin
ariable
problem
The
transmission
problem
in
equations
I
are
I
solv
e
ed
trate
in
the
their
D
o
ula
wn
and
frames
oth
(i.e.
v
w
results
e
n
w
tests
ork
presen
in
In
Lagrangian
IV
v
e
ariables
t
This
3D
c
for
hoice
ulation
al
con
lo
ergence
ws
is
us
in
to
V
a
ally
v
Section
oid
is
the
oted
presence
some
of
and
con
I
v
The
ection
pr
terms
Let
whic

h
3
giv
e
e
Lipsc
rise
b
to
t
instabilities
in
in
w
the
directions
n
en
umerical

solution

when
b
the
its
sp
ving
eed
and
of
2
the
complemen
mo
in
ving

part
=
is
n
high

In
.
this
set
con
1
text
a
w
metric
e
domain
prop
included
ose

a
.
new
e
sliding
the
mesh
terface
metho
=
d

based
,
on
the
the
w
mor
normal
tar

elemen
at
t
Let
metho
(
d
)
tec
e
hnique
rotation
[1],
at
[4-3]
t
whic
r
h
:
handles
1
in

an
the
optimal
op
w
whic
a
turns
y
domain
n
1
a
the
sense

that
t
will
and
b
t
y
in
elemen
erse
meshnite
Fig
in
structures
A
sliding
curren
mo
of
t
This
ts
a
a
ersal
eddy
of
e
Calculation
clear
so
Fig
on
is
non
transv
matc
section
hing
the
grids
In
at
mo
the
geometries
in
ulation
terface
system
b
at
et
t
w
0
een
tak
the
as
static
conuration
and
the
mo
solution
ving
time
part
>
No
the
constrain
domain
t
rotated
needs
an
to

b
!
e
and
imp
e
osed
bac
b
to
et
reference
w
thanks
een
the
the
erator
mesh
t
size
In
h
x
and
(
the
)
time
r
step
(

r
t
and
at
r
the
=
statorotor
t
in
x
terface
(
In
)).
space
eddy
the
t
2
equations
D
the
problem
netic
form
H
ulation
primary
is
ariable
dis
@
cretized
H
b
+
y

means
curl
of
=
Lagrangian
(
ite
1
elemen
i
ts
;
and
i
the
I
3
H
D
x
one
t
b
^
y
@
applying
=
the
;
second

family
I
of
(2)
st
(
order
;
edge
=
elemen
0
ts

ee
f
[6]
g
for
with
the
w
deition
transmission
If
imp
the
on
appro
I
ximation
t
do
1
es
r
not
x
allo
t
w
^
for
=
nonatc
2
hing
x
grids
t
either
^
a
;
o
r
c
curl
ktep
1
tec
r
hnique
x
ee
t
[5]),
^
or
=
an
H
o
(
v
;
erlapping
)
one
n
ee
(5)I