19.11.2014 Views

Calculating Iron Losses Taking into Account Effects ... - COMSOL.com

Calculating Iron Losses Taking into Account Effects ... - COMSOL.com

Calculating Iron Losses Taking into Account Effects ... - COMSOL.com

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

2<br />

1.5<br />

s=0MPa<br />

s=20MPa<br />

s=40MPa<br />

B[T]<br />

1<br />

0.5<br />

Figure 3. von mises stress field<br />

0<br />

0 500 1000 1500<br />

H[A/m]<br />

Figure 4. B-H curve as function of stress<br />

2<br />

1.5<br />

B[T]<br />

1<br />

0.5<br />

Figure 4. von mises stress along a radius<br />

The line plot shows <strong>com</strong>pressive stresses inside<br />

the stator and tensile stresses in the housing, as<br />

expected.<br />

This stress state is subsequently used to modify<br />

the magnetic properties of the stator material.<br />

6. Modelling the magnetic properties<br />

For simulation purposes, the magnetization<br />

curve of electrical steel can be expressed as a<br />

function of just two parameters: the saturation<br />

polarization Js and the relative permeability µ r :<br />

2 J<br />

s<br />

π ( µ<br />

r<br />

−1)<br />

µ<br />

0<br />

H<br />

B( H ) = µ<br />

0<br />

H + arctg<br />

π<br />

2 J<br />

s<br />

The effect of the stress can be introduced by a<br />

correction factor:<br />

B H S B H K B<br />

( , ) = ( ,0) (1 + S)<br />

With the factor KS = 0.0035 determined by<br />

fitting to experimental data, stresses being<br />

expressed in MPa.<br />

S<br />

0<br />

0<br />

-20<br />

S[MPa]<br />

-40<br />

-60<br />

0<br />

200<br />

400<br />

H[A/m]<br />

As the emqa equation system needs H as a<br />

function of B, the above function has to be<br />

inversed.<br />

H[A/m]<br />

10000<br />

8000<br />

6000<br />

4000<br />

2000<br />

0<br />

0<br />

-20<br />

S[MPa]<br />

-40<br />

-60<br />

0<br />

This H(B,S) function is implemented in the<br />

Comsol model as a 2D interpolated table-lookup<br />

function.<br />

7. Solving for the magnetic flux density<br />

2D models of electrical machines are well suited<br />

for the perpendicular induction currents, vector<br />

potential (emqa) application mode: the currents<br />

are perpendicular to the modelling plane. This<br />

0.5<br />

B[T]<br />

1<br />

600<br />

1.5<br />

800<br />

1000<br />

2

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!