In [1]:
    
%pylab inline
%precision 4
from scipy import constants as c   # we like to use some constants
    
    
The core shown in Figure P1-4 below:
and is made of a steel whose magnetization curve is shown in Figure P1-9 below:
Repeat Problem 1-7, but this time do not assume a constant value of $\mu_r$.
In [2]:
    
N1 = 600
N2 = 200
i1 = 0.5           # A
i2 = 1.0           # A
    
The two coils on this core are wound so that their magnetomotive forces are additive, so the total magnetomotive force on this core is $$\mathcal{F}_\text{TOT} = N_1 i_1 + N_2 I_2$$
In [3]:
    
F_tot = N1 * i1 + N2 * i2 
print('F_tot = {:.1f} At'.format(F_tot))
    
    
Therefore, the magnetizing intensity $H$ is: $$ H = \frac{\mathcal{F_\text{TOT}}}{l_c}$$
In [4]:
    
lc = 4 * (0.075 + 0.5 + 0.075)   # [m] core length on all 4 sides.
H = F_tot / lc
print('H = {:.1f} At/m'.format(H))
    
    
From the magnetization curve:
In [5]:
    
B = 0.17 # [T]
    
and the total flux in the core is: $$\phi_\text{TOT} = BA$$
In [6]:
    
A = 0.15**2  # [m²] 
phi_tot = B*A
print('ϕ = {:.3f} mWb'.format(phi_tot*1000))
    
    
The relative permeability of the core can be found from the reluctance as follows:
$$\mathcal{R}_\text{TOT} = \frac{\mathcal{F}_\text{TOT}}{\mathcal{\phi}_\text{TOT}} = \frac{l_c}{\mu_0 \mu_r A}$$Solving for $\mu_r$ yields:
$$\mu_r = \frac{\mathcal{\phi}_\text{TOT}l_c}{\mathcal{F}_\text{TOT}\mu_0 A}$$
In [7]:
    
mu_r = phi_tot * lc / (F_tot * c.mu_0 * A)
print('μ = {:.1f}'.format(mu_r))
    
    
The assumption that $\mu_r = 1200$ is not very good here. Actually, it is not very good in general.