In [1]:
%pylab inline
%precision %.4g
from scipy import constants as c # we like to use some constants
A ferromagnetic core is shown in Figure P1-3:
The relative permeability of the core is:
In [2]:
mu_r = 1500
mu = mu_r * c.mu_0
The dimentions are as shown in the diagram, and the depth of the core is 5 cm. The air gaps on the left and right sides of the core are 0.050 and 0.070 cm, respectively. Because of fringing effects, the effective area of the air gaps is 5 percent larger than their physical size.
If there are 300 turns in the coil wrapped around the center leg of the core and if the current in the coil is 1.0 A,
This core can be divided up into five regions. Let:
If we assume that the mean path length of the flux is in the center of each leg of the core, and if we ignore spreading at the corners of the core, then the path lengths are:
In [3]:
l1 = 1.11 # m
l2 = 0.0007 # m
l3 = 1.11 # m
l4 = 0.0005 # m
l5 = 0.37 # m
The areas can be calculated as:
In [4]:
A1 = 0.07 * 0.05 # m^2
A2 = 0.07 * 0.05 * 1.05 # m^2 5 percent larger than the physical size
A3 = 0.07 * 0.05 # m^2
A4 = 0.07 * 0.05 * 1.05 # m^2 5 percent larger than the physical size
A5 = 0.07 * 0.05 # m^2
And the reluctances are hence:
In [5]:
R1 = l1 / (mu * A1) # At /Wb = At/Vs
R2 = l2 / (c.mu_0 * A2) # At /Wb = At/Vs
R3 = l3 / (mu * A3) # At /Wb = At/Vs
R4 = l4 / (c.mu_0 * A4) # At /Wb = At/Vs
R5 = l5 / (mu * A5) # At /Wb = At/Vs
print('R1 = {:>5.1f} kAt/Wb'.format(R1/1000) )
print('R2 = {:>5.1f} kAt/Wb'.format(R2/1000) )
print('R3 = {:>5.1f} kAt/Wb'.format(R3/1000) )
print('R4 = {:>5.1f} kAt/Wb'.format(R4/1000) )
print('R5 = {:>5.1f} kAt/Wb'.format(R5/1000) )
Then the total reluctance of the core is $\mathcal{R}_\text{TOT} = \mathcal{R}_5 + \frac{(\mathcal{R}_1 + \mathcal{R}_2)(\mathcal{R}_3 + \mathcal{R}_4)}{\mathcal{R}_1 + \mathcal{R}_2 + \mathcal{R}_3 + \mathcal{R}_4}$.
In [6]:
Rtot = R5 + ((R1 + R2) * (R3 + R4)) / (R1 + R2 + R3 + R4)
print('Rtot = {:.1f} kAt/Wb'.format(Rtot/1000))
and the magnetomotive force is $\mathcal{F} = \mathcal{N} \mathcal{I}$:
In [7]:
N = 300 # t given in description
I = 1.0 # A given in description
F = N * I
The total flux in the core $\phi_\text{TOT}$ is equal to the flux in the center leg $\phi_\text{center} = \frac{\mathcal{F}}{\mathcal{R}_\text{TOT}}$ , which is:
In [8]:
phi_cent = F / Rtot
print('phi_center = {:.3f} mWb'.format(phi_cent*1000))
The fluxes in the left and right legs can be found by the "flux divider rule", which is analogous to the current divider rule.
Thus the flux in the left leg $\phi_\text{left} = \frac{(\mathcal{R}_3 + \mathcal{R}_4)}{\mathcal{R}_1 + \mathcal{R}_2 + \mathcal{R}_3 + \mathcal{R}_4}\phi_\text{TOT}$ is:
In [9]:
phi_left = (R3 + R4) / (R1 + R2 + R3 + R4) * phi_cent
print('phi_left = {:.3f} mWb'.format(phi_left*1000))
The flux in the right leg $\phi_\text{right} = \frac{(\mathcal{R}_1 + \mathcal{R}_2)}{\mathcal{R}_1 + \mathcal{R}_2 + \mathcal{R}_3 + \mathcal{R}_4}\phi_\text{TOT}$ is:
In [10]:
phi_right = (R1 + R2) / (R1 + R2 + R3 + R4) * phi_cent
print('phi_right = {:.3f} mWb'.format(phi_right*1000))
The flux density $B = \frac{\phi}{A}$ in the left air gap is:
In [11]:
B2 = phi_left / A2
print('B2 = {:.3f} T'.format(B2))
The flux density $B = \frac{\phi}{A}$ in the right air gap is:
In [12]:
B4 = phi_right / A4
print('B4 = {:.3f} T'.format(B4))