- #1
SunnyBoyNY
- 63
- 0
(this is not a hw)
Assume you have a magnet of dimensions x_m, h_m, remanent flux density Br, and coercive field density Hc. The magnet is placed in a magnetic "C" structure (perfect iron) such that it is connected on one side but there is an airgap on the other side.
xxxxxxxx
xx... xx
xx...gg
xx...mm
xx...mm
xxxxxxxx
Stack length is 1 m for simplicity.
I know how to calculate gap flux density as a function of airgap length. I am struggling, however, with using the Maxwell Stress tensor to calculate the force between the magnet and the structure through the airgap.
This is what I tried:
[itex]
Bm = Br / (1+Br*g(h_m*u0*Hc))
[/itex]
[itex]
\nabla\cdot\sigma_{xyz}= \frac{1}{\mu_{0}} \begin{pmatrix}
\frac{\partial 0.5B_{x}^{2}}{\partial x} & \frac{\partial B_{x}B_{y}}{\partial x} & \frac{\partial B_{x}B_{z}}{\partial x}\\
\frac{\partial B_{y}B_{x}}{\partial y} & \frac{\partial 0.5B_{y}^{2}}{\partial y} & \frac{\partial B_{y}B_{z}}{\partial y}\\
\frac{\partial B_{z}B_{x}}{\partial z} & \frac{\partial B_{z}B_{y}}{\partial z}& \frac{\partial 0.5B_{z}^{2}}{\partial z} \\
\end{pmatrix}
[/itex]
care only about one direction, which simplifies the equation to 1/2u0 * dBm^2/dx.
Now integrate over volume:
[itex]
F = \int_{V} \mathbf{f} \mathrm{d} V = \int_{V} (\nabla\cdot\sigma)\mathrm{d} V = \oint_{S} (\mathbf{f}\cdot\mathbf{n})\mathrm{d} A
[/itex]
But the numbers consistently come out wrong - with respect to a FEA simulation. Do I use the MST incorrectly?
Thank you.
Assume you have a magnet of dimensions x_m, h_m, remanent flux density Br, and coercive field density Hc. The magnet is placed in a magnetic "C" structure (perfect iron) such that it is connected on one side but there is an airgap on the other side.
xxxxxxxx
xx... xx
xx...gg
xx...mm
xx...mm
xxxxxxxx
Stack length is 1 m for simplicity.
I know how to calculate gap flux density as a function of airgap length. I am struggling, however, with using the Maxwell Stress tensor to calculate the force between the magnet and the structure through the airgap.
This is what I tried:
[itex]
Bm = Br / (1+Br*g(h_m*u0*Hc))
[/itex]
[itex]
\nabla\cdot\sigma_{xyz}= \frac{1}{\mu_{0}} \begin{pmatrix}
\frac{\partial 0.5B_{x}^{2}}{\partial x} & \frac{\partial B_{x}B_{y}}{\partial x} & \frac{\partial B_{x}B_{z}}{\partial x}\\
\frac{\partial B_{y}B_{x}}{\partial y} & \frac{\partial 0.5B_{y}^{2}}{\partial y} & \frac{\partial B_{y}B_{z}}{\partial y}\\
\frac{\partial B_{z}B_{x}}{\partial z} & \frac{\partial B_{z}B_{y}}{\partial z}& \frac{\partial 0.5B_{z}^{2}}{\partial z} \\
\end{pmatrix}
[/itex]
care only about one direction, which simplifies the equation to 1/2u0 * dBm^2/dx.
Now integrate over volume:
[itex]
F = \int_{V} \mathbf{f} \mathrm{d} V = \int_{V} (\nabla\cdot\sigma)\mathrm{d} V = \oint_{S} (\mathbf{f}\cdot\mathbf{n})\mathrm{d} A
[/itex]
But the numbers consistently come out wrong - with respect to a FEA simulation. Do I use the MST incorrectly?
Thank you.