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
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