3-1

Soluzioni per i problemi del Capitolo 3

1. Campi Magnetici e Prodotti Cruzati

P3.1: Trova AxB per il seguente:
a. A = 2ax– 3ay+ 4az,B= 5ay- 1az
b. A=aρ + 2aφ + 4az,B= 2aρ + 6az
c. A= 2ar+ 5aθ + 1aφ,B=ar+ 3aφ

(a)
a x a ya z
A× B = 2 −3 4 = ( 3 − 20 ) a x − ( − 2 ) a y + ( 10 ) a z = − 17ax + 2u n+y10az
0 5 −1
(b)
a a un z
ρ φ

A× B = 1 2 4 = 12a − ( 6 − 8 ) aφ + ( 0 − 4 ) a z = 12a + 2a − 4az

ρ ρ φ

2 0 6
(c)
a run un θ φ

A× B = 2 5 1 = 15ar − ( 6 − 1) aθ + ( 0 − 5 ) aφ = 15ar − 5a − 5a
θ φ

1 0 3

P3.2: Se un parallelogramma ha un lato corto a, un lato lungo b e un angolo interno θ (il

l'area del parallelogramma è data da
area = assenzio
θ.
Determina come utilizzeresti il prodotto vettoriale di una coppia di vettori per trovare l'area di un
parallelogramma definito dai punti O(0,0,0), P(6,0,0), Q(8,12,0) e R(2,12,0).
(Assumere le dimensioni in metri)

area = assenzio
θ
= A × B = A B seno
θ
A = 6ax,B=2ax+ 12ay
AxB= 72az,
Area = 72 m2

Fig. P3.2
 3-2

P3.3: Dati i vertici di un triangolo P(1,2,0), Q(2,5,0) e R(0,4,7), trova (a) l'interno
angoli, (b) un vettore unitario normale alla superficie che contiene il triangolo e (c) l'area di
il triangolo.

(a) PQ = 1a x + 3un; y PQ = 3.16

PR = − 1ax + 2a y + 7a z ; PR7.348
=
PQ × PR = 21ax − 7a y + 5az
PQ × PR = 22,69
A × B = A B seno
θ ab
22,69
peccato
θP = ;
( 3.16 ) ( 7.348 )
θ P = 78o
QP = − 1a x − 3a y ; QP = 3.16 Fig. P3.3

QR = − 2a x − 1a y + 7az ; QR = 7.348
QP × QR = 22,69 θ Q = 78o ; θ R = 180o−θ−θ P Q = 24o
PQ × PR
(b) a n = = 0.93ax − 0.31ay + 0.22az
PQ × PR
1
(c) area = PQ × PR = 11,4 metri2
2

2. Legge di Biot-Savart
P3.4: Un segmento di conduttore sull'asse z si estende da z = 0 a z = h. Se questo segmento
conduce corrente I nel +azdirection, findH(0,y,0). Compare your answer to that of
Esempio 3.2.

Utilizziamo l'equazione (3.7) e cambiamo i limiti:

h
Ioρ a z
H=
φ

4π ρ 2
z2 + ρ 2
0

Iy ( − a x ) h
=
4π y 2h2 + y2
− Ihun x
H=
4π y h2 + y 2
Nota che se la linea di corrente è semi- Fig. P3.4
infinito (va da z = 0 a z = ∞), noi
avere
 3-3

ax
− Io
H=
4π y

P3.5: Una linea di lunghezza infinita con una corrente di 2.0 A nel +axla direzione esiste aty= -3.0 m,z
= 4,0 m. Una seconda linea di lunghezza infinita con una corrente di 3,0 A nel +azla direzione esiste a tx=
0,y= 3.0 m. Trova H(0, 0, 0).

Questa situazione è mostrata nella figura P3.5a.

Ho=H1+H2
Riferendosi alla figura,
R = 3ay- 4az, R = 5,
aR= 0.6ay-0,8az
aφ =axxaR= 0.80ay+ 0,60az
Io1
H1 = a
2π ρ
φ

2A ( 0,80ay + 0.60az ) Fig. P3.5a

=
2π 5m
= 51ay + 38az mA
m
Io 2
H2 = a
2π ρ
φ

3Aa x
= = 159ax mA
2π 3m m
Ho= 159ax+ 51ay+38azmA/m

O con 2 cifre significative

Ho= 160ax+ 51ay+38azmA/m

Fig. P3.5b

P3.6: Un circuito conduttore a forma di triangolo equilatero con lato di 8,0 cm è centrato
nel piano x-y. Porta una corrente di 20,0 mA in senso orario quando vista dal +a z
direzione. TrovaH(0, 0, 16cm).

La situazione è illustrata nella Figura P3.6a. Il disegno nella Figura P3.6b è utilizzato per trovare il
Intercetta sull'asse x per il triangolo. Con semplice trigonometria abbiamo:
b
= tan30o ;e dato che a = 4 cm troviamo b = 2,31 cm.
a
 3-4

Ora per un segmento adattiamo l'Equazione (3.7):

+a
a
Io z
H1 =
φ

4π ρ z2 + ρ 2
−a

conρaρ = RaR,
R = − 2.31ax16a
+ z , R = 16,17 cm, a φ
= − u n×aay R = − 0.99ax0,143a
− z

H1 =
( 20x10 A ) ( − 0.99a 0.143a
−3
− ) x z 2 ( 4)
4π ( 16,17x10 m ) −2
42 + 16.172
H1= -4,7ax - 0,68az mA/m
Ora, per simmetria, il totaleH contiene solo la azcomponent:
Htot= -2,0 mA/maz.

Fig. P3.6a Fig. P3.6b

P3.7: Un circuito conduttivo quadrato di lato 10,0 cm è centrato nel piano x-y. Trasporta
Corrente di 10.0 mA in senso orario se vista dal +azdirezione. Trova H(0, 0, 10cm).

Troviamo H per una sezione del quadrato

adottando l'equazione (3.7):
+a
a
Io z
H=
φ

4π ρ z2 + ρ 2
−a

Iaun
H=
φ

a
2π ρ 2 + ρ 2
Conρaρ = RaR, abbiamo
R= -5ax+ 10az, |R| = 11.18x10-2m
unR= -0.447ax+ 0.894az,
aφ = -ayxaR= -0.894ax- 0,447az Fig. P3.7
 3-5

H1 =
( 10x10 ) ( 5) ( − 0.894a − 0,447a ) = − 5.2a
3
−
x z
− 2.6az mA
2π ( 11,18x10 ) 5 + 11.18 − 2 2 2
x
m

Ora, per simmetria, il totaleH contiene solo il azcomponent:

HTOT= -10.4azmA/m

P3.8: Un circuito conduttivo nel piano x-y è delimitato da ρ= 2.0 cm, ρ= 6.0 cm, φ= 0°
eφ= 90°. 1.0 A di corrente scorre nel circuito, andando verso laaφ direzione su ρ= 2.0
cm braccio. Determina l'origine.

Dall'ispezione della figura, vediamo

che solo le porzioni dell'arco di
il ciclo contribuisce a H.
Da un esempio di anello abbiamo:
2π
Ia2a z
H= 3 ∫ dφ
4π ( h 2a+2 ) 2
0

Per il ρ= un segmento del ciclo:

/2
Ia2a z π Io
3 ∫
Ha = dφ = az
4π a 0 8a
Ib2a z − Io
0
Fig. P3.8
Atρ= b: H b = 3 ∫
dφ = az ;
4π b π / 2 8b
I1 1 1 1 1 A
HTOT =
Quindi − az = − a z = 4.2az
8a b 8 0,02 0,06 m

P3.9: MATLAB: Quanto vicino devi essere al centro di una lunghezza finita di corrente-
portare la linea prima che appaia infinita in lunghezza? Considera Hf(0,a, 0) è il campo per il
linea finita di lunghezza 2h centrata sull'asse z, e che Hio(0,a, 0) è il campo per un
linea di corrente di lunghezza infinita sull'asse z. In entrambi i casi considera corrente I nella direzione +az
direzione. Disegna Hf/Hiovsh/a.

Adattando l'equazione (3.7), per la linea di lunghezza finita abbiamo:

a
Io h
Hf =
φ

2π ρ h2 + ρ 2
ρ =a

Per la lunghezza infinita della linea:

a
Io
H io=
φ

2π ρ
 3-6

Il rapporto che desideriamo tracciare è:

Hf
=
h
=
( ha )

Hi
( ha
2
h2 + a2 ) +1

La routine MATLAB segue.

% M-File: MLP0309

%Considera il campo per una linea finita di lunghezza 2h

%orientato sull'asse z con corrente I in direzione +z.

Il campo si trova a una distanza a da

%la corrente sull'asse y (punto (0,a,0)).

Vogliamo confrontare questo campo con quello di un infinito

%lunghezza linea di corrente.

%Traccia Hf/Hi contro h/a. Ci aspettiamo che man mano che h/a cresce grande,
%il linea volontà apparire di più infinito to un osservazione
%puntare a (0,a,0).

hova=0.01:.01:100;
HfovHi=hova./sqrt(1+(hova).^2);

semilogx(hova,HfovHi)
etichetta asse x('h/a')
etichetta asse y('Hf/Hi')
griglia attiva

Fig. P3.9a Fig. P3.9b

 3-7

P3.10: MATLAB: Per l'anello di corrente descritto in MATLAB 3.2, trova

following points (a) (0, 0, 1m), (b) (0, 2m, 0), and (c) (1m, 1m, 0).

%M-File: MLP0310

Trova l'intensità del campo magnetico in qualsiasi punto di osservazione

%risultante da un anello di raggio a e corrente I,
%nella direzione aphi centrata nel piano x-y.
df=1; %incremento in gradi
a=1; %raggio dell'anello in m
I=1; %corrente in A
Ro=input('vector location of observation point: ');

per j=1:df:360;
Fr=j*pi/180;
Rs=[a*cos(Fr) a*sin(Fr) 0];
as=unitvector(Rs);
dL=a*df*(pi/180)*cross([0 0 1],as);
Rso=Ro-Rs;
aso=unitvector(Rso);
dH=I*cross(dL,aso)/(4*pi*(magvector(Rso))^2);
dHx(j)=dH(1);
dHy(j)=dH(2);
dHz(j)=dH(3);
fine
H=[sum(dHx) sum(dHy) sum(dHz)]

Ora per eseguire il programma:

MLP0310
posizione vettoriale del punto di osservazione: [0 0 1]

H = -0.0000 -0.0000 0.1768

MLP0310
posizione vettoriale del punto di osservazione: [0 2 0]

H= 0 0 -0.0431

MLP0310
posizione vettoriale del punto di osservazione: [1 1 0]

H= 0 0 -0.1907
 3-8

MLP0310
posizione vettoriale del punto di osservazione: [0 1 1]

H = 0.0000 0.0910 0.0768

>>

Quindi vediamo:
(a)H= 0,18azA/m
(b)H= -0.043azA/m
(c)H= -0.19azA/m
(extra)H= 9.1ay+ 7.7azmA/m

P3.11: Un solenoide ha 200 avvolgimenti, è lungo 10.0 cm e ha un raggio di 1.0 cm. Assumendo
1.0 A di corrente, determina l'intensità del campo magnetico al centro esatto del solenoide.
Come si confronta con la tua soluzione se fai l'assunzione che 10 cm >> 1
cm?

Eq. (3.10):

NI h− z z
H= + az
2h 2
( h − z) + a 2 z2 + a2

200 ( 1A) 0,1 − 0,05 0,05 A

H= + a z = 1961az
2 ( 0,1m ) ( 0,1 − 0,05)
2
+ 0,012 0,0520,01
+ 2 m

OrH= 1960 A/maz

La soluzione approssimativa, assumendo 10cm >> 1cm, è

NI 200 ( 1A) A
H= az = a z = 2000az
h 0,1m m

P3.12: MATLAB: Per il solenoide del problema precedente, traccia la magnitudine della
campo rispetto alla posizione lungo l'asse del solenoide. Includere l'asse 2 cm oltre ciascuno
fine del solenoide.

% File M: MLP0312
%
% Traccia H vs lunghezza attraverso il centro di un solenoide
%
 3-9

clc
chiaro

% inizializzare le variabili
N=200; %numero di giri
h=0,10; %altezza del solenoide
a=0.01; %raggio del solenoide
I=1; %corrente
dz=0.001; %cambiamento passo in z

z=-0,02:dz:h+0,02;
zcm=z.*100;
A1=(h-z)./sqrt((h-z).^2+a^2);
A2=z./radice(z.^2+a^2);
Atot=A1+A2;
H=N*I.*Atot/(2*h);

% generare tracciato
traccia(zcm,H)
etichetta asse x('z(cm)')
etichetta dell'asse y('H (A/m)')
griglia attivata

Fig. P3.12

P3.13: Un nastro di corrente largo 4,0 cm è centrato attorno all'asse y nel piano x-y e
ha una densità di corrente superficiale K = 2πayA/m. Determina l'intensità del campo magnetico a
the point (a) P(0, 0, 2cm), (b) Q(2cm, 2cm, 2cm).

(a) Because of the symmetry (Figure

P3.13a), possiamo usare un'equazione modificata (3.14):
Ky 1 d
H= ax
abbronzato −

π a
2π 1 2 A
= a x = 1,57ax
abbronzare −

π 2 m
(b) Facendo riferimento alla Figura P3.13b;
Io
dH = a ;
2π ρ
φ

Fig. P3.13a
Ra r = ρ un ;dove R = ( d − x ) a x + aa z
ρ
 3-10

2
quindi
ρ = ( d − x) + a 2, e
aR =
( d − x) a x + a a z. Quindi
2
( d − x) + a 2
a una + x( x − d ) a z
a = a y × un R=
φ
2
( d − x) + a 2
K y ( a a x + xa z − da zdx )
H=
2π ∫ ( d − x) + a
2 2
Fig. P3.13b
Questo è separato in 3 integrali, ognuno
uno risolto tramite integrazione numerica,
risultando in:
H=1.1083ax+0.3032az- 1.1083az; or
H= 1.1ax- 0,80azA/m

3. Legge del Circuito di Ampère

P3.14: Una coppia di fogli di corrente di estensione infinita esiste a z = -2,0 m e a z = +2,0 m.
il foglio superiore ha una densità di corrente uniforme K = 3.0ayA/m
e quello in basso ha K= -3.0
ayA/m. Trova Hat (a) (0,0,4m), (b) (0,0,0) e (c) (0,0,-4m).

1
ApplichiamoH = K × aN ,
2
(a) H = 1( 3a y ) × a z + 1( -3ay ) × a z = 0
2 2
(b)
1 1
H = ( 3ay×) − az + ( -3ay ) × un z
2 2
A
= − 3ax
m
(c)H= 0 Fig. P3.14

P3.15: Una lamina di corrente di estensione infinita con K = 6.0ayA/m

esiste a z=0. Un conduttore
un anello di raggio 1,0 m, nel piano y-z centrato a z= 2,0 m, ha un campo magnetico nullo
intensità misurata al suo centro. Determina l'entità della corrente nel circuito e
mostra la sua direzione con un disegno.

Hfino a=HS+HL
 3-11

1 A
HS = K × a N = 1( 6ay ) × a z = 3ax
2 2 m
Per il ciclo, usiamo l'equazione (3.10):
Io
H= az
2a
where here
Io − Io
HL = ( − ax ) = ax
2a 2
(il segno è scelto oppostoH)S).
Quindi, I/2 = 3 e I = 6A.
Fig. P3.15

P3.16: Data il campo H= 3y2axtrova la corrente che passa attraverso un quadrato nel piano x-y
piano che ha un angolo all'origine e l'angolo opposto a (2, 2, 0).

Riferendosi alla Figura P3.6, valutiamo la circolazione di H attorno al percorso quadrato.

b c d a

Ñ
∫ HgdL = Io inc= ∫ +∫+∫+∫
a b c d
b

∫ 3 ( 0 )a gdxun
2
x x = 0
a
c

∫ 3y ungdya
b
2
x y = 0

d 0

∫ 3 ( 2 )a xgdxun x = 12dx
2

c
∫ = − 24 2
a

∫
d
= 0
Fig. P3.16
Quindi abbiamo ioenc= 24 A. Il negativo
Il segno indica che la corrente sta andando verso il
-azdirezione.

P3.17: Data un filo solido con raggio di 3.0 mm centrato sull'asse z con una distribuzione uniforme
2,0 ampere di corrente nel +azdirection, plot the magnetic field intensityHversus radial
distanza dall'asse z nell'intervallo 0 ≤ρ≤ 9 mm.

La figura P3.17 mostra la situazione insieme ai percorsi amperiani. Abbiamo:

∫ HgdL = Io inc,dove H = H a e d L = ρ dφ un; Q 2π ρ Hφ = Io enc
Ñ φ φ φ
 3-12

Questo sarà vero per ogni percorso amperiano.

ρ 2π
Io Io Iρ 2
AP1: Io inc= ∫ J gdS , J =
π un
a2 z,Ioenc =
π a 2 ∫0
ρ dρ ∫
0
dφ =
a2
Io
H=
Quindi: ρ a per ρ ≤ a
2 a2
φ
π

Io
AP2: IoincIo, H = un per ρ ≥ a
2π ρ
φ

Fig. P3.17a Fig. P3.17b

% MLP0317
% generare grafico per problema ACL

a=3e-3; %raggio del filo solido (m)

I=2; %corrente (A)
N=30; %numero di punti dati da tracciare
rmax=9e-3; %raggio massimo per il grafico (m)
dr=rmax/N;

per i=1:round(a/dr)
r(i)=i*dr;
H(i)=(I/(2*pi*a^2))*r(i);
fine
per i=round(a/dr)+1:N
r(i)=i*dr;
H(i)=I/(2*pi*r(i));
fine
traccia(r, H)
xlabel('rho(m)')
etichette_y('H (A/m)')
griglia attivata
 3-13

P3.18: Data un filo solido con raggio di 2,0 cm centrato sull'asse z con una densità di corrente J=
3ρA/cm2az(per ρ in cm) traccia l'intensità del campo magnetico H in funzione della distanza radiale da
l'asse z nell'intervallo 0 ≤ρ≤ 8 cm.

Lasceremo a = 2 cm.
Ñ
∫ HgdL = Io inc,dove H = H un e d φ φ
L = ρ dφ a ; Q 2π ρ Hφ = Io enc
φ

∫ J gd S = ∫ 3ρ a
3
AP1 (ρ<a): Io inc= ρ
gρ d ρ dφ a z = 2π ρ

e H = ρ 2a per
φ
ρ ≤ un
a3
3
AP2 (ρ>a):Ienc= 2πa, quindi H= a per ρ ≥ a
ρ
φ

La routine di tracciamento di MATLAB è la seguente:

% MLP0318
% genera grafico per problema ACL

a=2; %raggio del filo solido (cm)

N=40; %numero di punti dati da tracciare
rmax=8; %raggio massimo per il grafico (cm)
dr=rmax/N;

per i=1:round(a/dr)
r(i)=i*dr;
H(i)=r(i)^2;
fine
per i=round(a/dr)+1:N
r(i)=i*dr;
H(i)=a^3/r(i);
fine
traccia(r,H)
etichetta dell'asse x('rho(cm)')
etichetta dell'asse y('H (A/cm)')Fig. P3.18
griglia accesa

P3.19: Un guscio cilindrico metallico infinitesimamente sottile di raggio 4,0 cm è centrato su

l'asse z e trasporta una corrente distribuita uniformemente di 10,0 mA nell'+azdirezione. (a)
Determina il valore della densità di corrente superficiale sulla superficie conduttrice e (b) traccia H
come funzione della distanza radiale dall'asse z nell'intervallo 0 ≤ρ≤ 12 cm.
 3-14

Io 10mA mA; mA
(a) K s = = = 39.8 quindiK = 40az
2π a2 π ( 0,04 m) m m
(b) per ρ < a, H = 0. Per ρ > a abbiamo:
Io
H= a
2π ρ
φ

La routine MATLAB per generare il grafico è la seguente:

% MLP0319
% generare grafico per il problema ACL
a=4; %raggio del filo solido (cm)
N=120; %numero di punti dati da tracciare
I=10e-3; %corrente (A)
rmax=12; %raggio massimo del grafico(cm)
dr=rmax/N;

per i=1:tondo(a/dr)
r(i)=i*dr;
H(i)=0;
fine
per i=round(a/dr)+1:N
r(i)=i*dr;
H(i)=100*I/(2*pi*r(i));
fine
traccia(r,H)
etichetta asse x('rho(cm)')
etichetta asse y('H (A/m)')
griglia attiva

Fig. P3.19a Fig. P3.19b

 3-15

P3.20: Un tubo cilindrico con uno spessore della parete di 1,0 cm e un raggio interno di 4,0 cm è
centrato sull'asse z e ha una corrente di 3,0 ampere distribuita uniformemente in +az
direzione. Traccia l'intensità del campo magnetico H rispetto alla distanza radiale dall'asse z su
l'intervallo 0 ≤ ρ ≤ 10 cm.

Per ogni percorso amperiano:

Ñ
∫ HgdL = Io ,dove H = H a e d
insegna φ φ
L = ρ dφ a ; Q 2π ρ Hφ = Io enc
φ

Ora, per ρ < a, Ienc= 0 quindi H = 0.

Io

)a ze
Fora <ρ< b, Io enc= ∫ J gdS, dove J =
π ( b − a2
2
S =ρ d ρ d φ a z

Io
ρ 2π
ρ 2 − a2 Io ( ρ 2 − un)2
π ( b − a 2 ) ∫a ∫
Io inc= ρ dρ dφ = Io 2 2 , H = a
b −a 2π ( b 2 − a 2 ) ρ
φ
2
0

Fig. P3.20a Fig. P3.20b

% MLP0320
% generare grafico per il problema ACL

a=4; %raggio interno del tubo (cm)

b=5; %raggio esterno del tubo (cm)
N=120; %numero di punti dati da tracciare
I=3; %corrente (A)
rmax=10; %raggio massimo per il grafico (cm)
dr=rmax/N;
aoverdr=a/dr
boverdr=b/dr

per i=1:round(a/dr)
r(i)=i*dr;
H(i)=0;
 3-16

end

per i=round(a/dr)+1:round(b/dr)
r(i)=i*dr;
num(i)=I*(r(i)^2-a^2);
den(i)=2*pi*(b^2-a^2)*r(i);
H(i)=100*num(i)/den(i);
fine

per i=round(b/dr)+1:N
r(i)=i*dr;
H(i)=100*I/(2*pi*r(i));
end

traccia(r,H)
etichetta asse x('rho(cm)')
etichetta asse y('H (A/m)')
griglia accesa

P3.21: Una linea di lunghezza infinita trasporta corrente I nella direzione +azdirezione lungo l'asse z, e questo
è circondato da una copertura cilindrica di lunghezza infinita (centrata attorno all'asse z) di raggio a
trasportando la corrente di ritorno Iin nella -azdirezione come corrente di superficie. Trova espressioni per
l'intensità del campo magnetico ovunque. Se la corrente è di 1.0 A e il raggio è di 2.0 cm,
tracciare la magnitudine di H rispetto alla distanza radiale dall'asse z da 0,1 cm a 4 cm.

Io
∫ HgdL = Io enc; per 0<ρ < a H =
2π ρ
a e per
φ
ρ > a H = 0.

Fig. P3.21a Fig. P3.21b

La routine MATLAB utilizzata per generare

La figura P3.21b è la seguente:
% MLP0321
 3-17

% generare grafico per il problema ACL

clc
clear

a=2; %raggio interno del cilindro (cm)

N=80; %numero di punti dati da tracciare
I=1; %corrente (A)
rmax=4; %raggio massimo per il grafico (cm)
dr=rmax/N;

per i=1:40
r(i)=.1+(i-1)*dr;
H(i)=100*I/(2*pi*r(i));
fine

per i=40:N
r(i)=i*dr;
H(i)=0;
fine

traccia(r,H)
etichetta asse x('rho(cm)')
etichettaY('H (A/m)')
griglia attiva

P3.22: Considera una coppia di gusci cilindrici collineari centrati sull'asse z. L'interno
Il guscio ha un raggio e trasporta una corrente superficiale totale di I ampere nella direzione +azdirezione mentre il
il guscio esterno di raggio b trasporta la corrente di ritorno I nell'azdirezione. Trova espressioni
per l'intensità del campo magnetico ovunque. Se a = 2 cm, b = 4 cm e I = 4 A, traccia il
magnitudo di H rispetto alla distanza radiale dall'asse z da 0 a 8 cm.

∫ HgdL = Io enc; per 0<ρ < a H = 0;

Io
per a<ρ < b H = un;
2π ρ
φ

e per ρ > b H = 0.
La routine MATLAB utilizzata per generare la Figura P3.22b è la seguente:
% MLP0322
% generare grafico per il problema ACL

a=2; %raggio interno del coax (cm)

b=4; %raggio esterno del coassiale (cm)
 3-18

N=160; %numero di punti dati da tracciare

Fig. P3.22a Fig. P3.22b

I=4; %corrente (A)
rmax=8; %raggio massimo per il grafico (cm)
dr=rmax/N;
aoverdr=a/dr
boverdr=b/dr

per i=1:round(a/dr)
r(i)=i*dr;
H(i)=0;
fine

per i=round(a/dr)+1:round(b/dr)
r(i)=i*dr;
H(i)=100*I/(2*pi*r(i));
fine

per i=round(b/dr)+1:N
r(i)=i*dr;
H(i)=0;
fine

traccia(r,H)
etichetta x('rho(cm)')
etichetta dell'asse y('H (A/m)')
griglia accesa

P3.23: Considera il toroide nella Figura 3.55 che è avvolto strettamente con N spire di
filo conduttivo. Per un percorso amperiano con raggio inferiore a a, nessuna corrente è racchiusa
e quindi il campo è zero. Allo stesso modo, per raggi maggiori di c, la corrente netta
 3-19

incluso è zero e di nuovo il campo è zero. Usa la Legge Circuitale di Ampère per trovare un
espressione per il campo magnetico a raggio b, il centro del toro.

∫ HgdL = Io enc;
All'interno del toroide, H=Hφ unφ, quindi

∫ Hgd L = ∫ H a φ φ
daρ bH
b φ
ρ =
.φ φ
= 2π
Poi, ioencper il percorso amperiano è: Iinc=NI.
NI
H= a .
2π b
φ

4. Curl e la Forma Puntuale della Legge Circuitale di Ampère

P3.24: Trova × A per i seguenti campi:
a.A=3xy2/zax
b.A=ρsin2φaρ –ρ2zcosφaφ
c.A=r2sinθar+r/cosφaθ

∂ Ax ∂ Ax − 3xy2 6xy
(a) × ( 3xy2z ) un
x = ay − az = ay − az
∂z ∂y z 2
z
∂ Aφ ∂ Aρ 1 ∂ ∂ Aρ
(b) × ( ρ peccato
φ a −ρ2 2zcosφa )=− a + a + (ρ A)−
φ
az
∂z ∂z ρ ∂ρ ∂φ
ρ φ ρ φ

3ρ 2zcosφ 1
= ρ 2cosφa + 0a − a z − ρ 2sinφ cosaφz
ρ ρ
ρ φ

= ρ 2cosφa − [ 3ρ zcos φ + 2senφ cos φ ] a z

ρ

(c)
r
θ ar +
× r 2peccato a
cos φ
θ

−1 ∂ Aθ 1 1 ∂ Ar 1 ∂ ( rAθ ) ∂ Ar
= ar + a + − a
rsin θ ∂φ rsin θ ∂φ r ∂r ∂θ
θ φ

− r ∂ ( cos φ )
1 −
1 ∂ r2 1 ∂
=
rsin θ ∂φ
un+r
r ∂ rcos φ
a −
r ∂θ
φ ( r sinθ ) a
2
φ

− peccato
φ 2
= a r + − rcos θ a
peccato
θ cos2 φ cos φ
φ
 3-20

P3.25: Trova Jat (3m, 60°, 4m) per H= (z/sinφ)aρ – (ρ2/cosφ)azA/m.

1 ∂ Az ∂ Aρ ∂ Az − 1 ∂ Aρ
× H= a + − a + un z
ρ ∂φ ∂ z ∂ρ ρ ∂φ
ρ φ

-ρsen 1φ 2ρ zcos φ
= + a + a + az
cos φ 2
peccato
φ cos φ ρ peccato
φ 2
ρ φ

Ora trova J valutando × H al punto dato:

A
J = − 10a 13a
+ 0,89a+z
m2
ρ φ

P3.26: Supponi H=y2ax+x2ayA/m.

[Link] Ñ
∫ HgdL attorno al percorso A → B → C → D → A, dove A(2m,0,0),
B(2m,4m,0), C(0,4m,0) e D(0,0,0).
[Link] questo Ñ
∫ HgdL dalla superficie S (2m*4m = 8m2).
[Link] × H al punto centrale.
d. Commenta i tuoi risultati per (b) e (c).

B C D A

(a) Riferendosi alla figura, valutiamo ∫ H gdL = Io enc= ∫

A
+∫+∫+∫
B C D
B 4

∫ = ∫ x dy
2
= 4 ( 4 ) = 16C
A 0 x= 2
C 0

∫
B
= ∫y
2
2
dx = − 32C
y= 4
D 0

∫
C
= ∫ x dy
4
2
= 0
x= 0
A 2

∫
D
= ∫ y dy
0
2
= 0
y= 0

Fig. P3.26
∫
Quindi abbiamo H gdL = − 16C

(b) dividendo per S = 8m2abbiamo -2 C/m2

(c) Valutare il rotore di H:
∂ Ay ∂ Ax
× H= − a z = ( 2x − 2 y ) a z, e al punto centrale (x = 1 e y = 2) abbiamo
∂x ∂y
C
× Hcentro = − 2 az
m2
(d) In questo caso particolare, × H = Ñ
∫ Hgd L S anche se S è di dimensioni apprezzabili.
 3-21

P3.27: Per l'esempio del cavo coassiale 3.8, abbiamo trovato:

Ioρ
per ρ ≤ a H= un,
2π a 2
φ

Io
per a<ρ ≤ b H= a ,
2π ρ
φ

2
Io c 2−ρ
per b<ρ ≤ c H= a ,
2π ( c 2 − b 2 ) ρ
φ

e per c< ρ, H =0.

a. Valuta il rotore in tutte e 4 le regioni.
b. Calcola la densità di corrente nelle regioni conduttive dividendo la corrente per
l'area. Questi risultati sono gli stessi di quelli che hai trovato in (a)?

1 ∂ Ioρ 2 Io
(a) × H= a =
2 z
a zper ρ ≤ a
ρ ∂ρ 2π un π a2
1 ∂ Io
× H= a z = 0 per un
< ρ < b
ρ ∂ρ 2π
1 ∂ I (c2−ρ 2
) − Io
× H=
ρ ∂ρ
× H = 0 per ρ > c
2π ( c − b 2 2
)
az =
π ( c2 − b2
)z a per b<ρ ≤ c

Io
(b) per ρ ≤ a,S=πa2, J = Io ) a z =
S
(
π a2
az

− Io
per b<ρ ≤ c,S=πc2 ( − b 2,) J =
π ( c2 − b2 )
az

Comment: × H = J è confermato.

P3.28: Supponiamo che tu abbia il campo H = r cosθaφ A/m. Ora considera il cono specificato
byθ=π/4, con un'altezza a come mostrato nella Figura 3.56. La parte superiore circolare del cono ha un
raggio.
a. Valuta il lato destro del teorema di Stokes attraverso il dS= dSa.θ superficie.
b. Valuta il lato sinistro del teorema di Stokes integrando attorno al circuito.

1 ∂ 1 ∂ ( rH φ )
(a) × H=
rsin θ ∂θ
( peccato )
θ H un φ r −
r ∂r
a θ

arderivata
∂ 1 ∂ ( cos θ
2
θ ) 2
− peccato
∂θ
( rsin θ (
cos θ ) = rcos 2
θ );
θ − peccato
rsin θ
2

∂θ
( peccato
θ H a)
φ r =
peccato
θ
ar
 3-22

Fig. P3.28
aθ derivata
− 1 ∂ ( r cosθ ) un = − 2 c o sθ un
2

r ∂r
θ θ

Quindi,
× H=
( cos θ
θ ) 2
− peccato 2

ar − 2 c o s θ a
peccato
θ
θ

Ora dobbiamo integrare questo su aθ surface:

( cos θ 2
θ ) 2
− peccato
a r − 2 c o sθ a grsin θ drd φ ( − a
∫ × HgdS = ∫ )
sin θ
θ θ

2a 2π
= ∫ 2rsin θ cos θ drd φ = 2 ( seno
θ cos θ ) θ =π ∫ rdr ∫ d φ = 2π a 2
4 0 0
2π

(b) Ñ
∫ Hgd L = ∫ r r= a 2
cos θ a gannuncio
φ
φ a = a 2 cosθ
φ
θ =
π ∫ dφ = 2π a 2
4 0

Chiaramente in questo caso la circolazione è l'approccio più semplice.

5. Densità di Flusso Magnetico

P3.29: Una linea di lunghezza infinita con una corrente di 3,0 A nel +ayla direzione si trova sull'asse y. Trova
the magnetic flux density at P(7.0m,0,0) in (a) Teslas, (b) Wb/m2, e (c) Gauss.

Io 3A mA
H= a = ( − a z ) = − 68 a z
2π ρ 2π ( 7m) m
φ

H A Wb Wb
B = µ o H = 4π x10 − 7
− 68 az = − 86x10 9 2 a z = − 86x10 9Ta z
− −

m m HA m

B = ( − ottantasei
10.000G
9 6
volte dieci=T−a860x10 Ga z
) T
− −
z
 3-23

P3.30: Supponiamo una lastra infinita di corrente con K = 12axA/m giace nel piano x-y
piano atz= 0. Trova B per qualsiasi punto sopra il foglio. Trova il flusso magnetico che passa
attraverso un 2m2area nel piano x-z per z > 0.

1
H= K × aN ;
2

B=
µo
K × aN =
( 4π x10 H7 m −
) 12ax
A
× a
6Wb
z = − 7.54x10 ay
−

2 2 m m2
Questo è valido in qualsiasi punto sopra il foglio.

Wb
Ora, φ = ∫ BgdS = BgS = -7,54x10 -6

m2
a y2m2 ( ( − un) ) = 15µWb
y

P3.31: Esiste un cavo coassiale di lunghezza infinita lungo l'asse z, con un guscio interno di
raggio che trasporta corrente I nell' +azdirezione e guscio esterno di raggio che trasporta il
restituire corrente. Trova il flusso magnetico che attraversa un'area di lunghezza lungo il z-
asse delimitato dal raggio tra a e b.

Io µ Io
H= a , B= o a ,
2π ρ 2π ρ
φ φ

µ oIo un µ Io
Per un <ρ< b, φ = ∫ BgdS = ∫ gd ρ dza = o ln ( ρ
φ
) bah
2π ρ 2π
φ

µ oIh b
φ = ln Wb
2π a

6. Forze Magnetiche
P3.32: Una carica di 1.0 nC con velocità di 100. m/sec in direzione y entra in una regione dove
L'intensità del campo elettrico è 100 V/[Link] la densità del flusso magnetico è 5,0 Wb/m2ax.
Determina il vettore forza che agisce sulla carica.
m Wb Wb
F = q ( E + u × B ) ; u × B = 100ay ×5 = − 500 az
s m ax
2
sm
9 V Wb Vs mN
F = 10 C100
−
az + − 500 az = − 400nNaz
m sm Wb VC

P3.33: Una carica di 10 nC con velocità 100 m/sec nella direzione z entra in una regione dove
l'intensità del campo elettrico è 800 V/mAxe la densità di flusso magnetico 12,0 Wb/m2ay.
Determina il vettore della forza che agisce sulla carica.
 3-24

V m Wb
F = q ( E + u × B ) = 10x10 9C800a x
−
+ 100az × 12 = − 4 µ Na x
m s m 2a y

P3.34: Una particella carica di 10 nC ha una velocità v = 3,0ax+ 4.0ay+ 5.0azm/sec così com'è
entra in un campo magnetico B = 1000. Tay(ricorda che un tesla T = Wb/m)2). Calcola la forza
vettore sulla carica.

−
( x + 4ay + 5az
F = q ( u × B ) = 10x10 9C3un ) ms × 1000
Wb
m2
ay

Il prodotto vettoriale:
a xa y az
3 4 5 = − 5000ax3000a
+ z

0 1000 0
Valutando troviamo: F = -50ax+ 30azµN

P3.35: Quale campo elettrico è necessario affinché la velocità della particella carica nel
il problema precedente rimane costante?

dv
F = ma = m = 0 (velocità costante)
dt
F = q ( E + u × B ) = 0;

E = − u × B = - (3a x + 4ay + 5az ) ms × 1000

Wb
m 2
V
a y = − 3000az
m
V
+ 5000ax
m
kV
E = 5ax − 3az
m

P3.36: Un elettrone (con massa a riposo Me= 9.11x10-31kg e carica q = -1,6 x 10-19C) ha un
velocità di 1,0 km/sec mentre entra in un campo magnetico di 1,0 nT. Il campo è orientato normale a
la velocità dell'elettrone. Determina la grandezza dell'accelerazione sull'elettrone
causato dal suo incontro con il campo magnetico.

F = ma = q ( u × B ) ;
q
a= ( u × B) ;
m
 3-25

un=
quB ( − 1.6x10C1000
=
− 19
) s) (
m 10 9 Wb
(
m 2 = 175x103 m
−
)
m ( 9.11x10 31kg ) − s2

P3.37: Supponiamo di avere una corrente superficiale K = [Link]/m lungo il piano z = 0. Circa un
un metro o giù di lì sopra questo piano, una particella caricata di 5,0 nC si sta muovendo insieme alla velocità v =
-[Link]/s. Determina il vettore forza su questa particella.

1 1 A
H= K × un N= 20ax × a z = − 10ay
2 2 m
B = µ o H = − 10 ( 4π x10 7 ) un
y
−Wb = − 40π x10 7a y Wb −

m2 m2
m Wb
F = qu × b = ( 5x10 9C ) − 10ax × − 40π x10 7a y 2 = 0.63pNaz
− −

s m

P3.38: A un metro circa sopra la corrente superficiale del problema precedente c'è un infinito
linea di lunghezza che conduce 1.0 A di corrente nell –a xdirezione. Determina la forza per unità
lunghezza che agisce su questa linea di corrente.

0
Benvinuto
F12 = ∫ Io dL× B
2 2 1 = Io∫2dxa x×
L
− 40 π x10 7a y
−

m2
F12 = Io(2 − L) ( − 40π x10 7a z ) ;−

F12 µ N
= 40π x10 7un z= 12.6 a z
−

L m

P3.39: Ricorda che la forza gravitazionale su una massa m F =è mg, dove, sulla terra
superficie, g = 9.8 m/s (-azUna linea di corrente di 2.0 A con una massa di 100 g per metro di lunghezza è
2

orizzontale rispetto alla superficie terrestre e è diretto da ovest a est. Quale magnitudo e
Quale direzione della densità di flusso magnetico uniforme sarebbe necessaria per levitare questo filo?

m Ns 2kg
Fg = mg; per 1 mF= g ( 100g9.8
) ( − az ) = − 0.98Naz
s2 kg m 1000g
1m
F = ∫
0
IdL× B = − Fg = + 0.98Naz

Per ispezione, B = Bo(-ax)

 3-26

∫ Idya y × Bo ( − a x ) = ILBo a z
o

Wb
= ( 2A1m
) ( )B o (
= 2Bo N )
m2
La conversione dell'unità per arrivare a
Newtons è come segue:
Am Wb Vs W J Nm
= N
m 2Wb VA Ws J

Quindi abbiamo Bo0.490

Wb/m2, e
B= 0,490 Wb/m2(-ax) (diretto a nord) Fig. P3.39

P3.40: Supponiamo di avere una coppia di linee parallele ciascuna con una massa per unità di lunghezza di 0,10
kg/m. Una linea si trova a terra e conduce 200. A nel +axdirezione, e l'altro
uno, 1.0 cm sopra il primo (e parallelo), ha una corrente sufficiente per levitare. Determina il
corrente e la sua direzione per la linea 2.

7
F12 µ oIo1 2Io 4π x10 H m200A −

Qui utilizzeremo = a y = Io 2a z= 4x10 3Io 2a z −

L 2π y 2π 0,01m
F mg Ns2 N
L
=
L
(
= ( 0,10 kg m )9,8 m s2
) kgm
= 0,98
m
Quindi risolvendo per I2:
0,98
Io 2= = 245Ain nel -axdirezione.
4x10 3 −

P3.41: Nella Figura 3.57, una corrente di 2.0 A è mostrata lungo l'asse z con la corrente in
il +azdirezione. Un circuito di corrente esiste nel piano x-y (z = 0) che ha 4 fili (etichettati
1 fino a 4) e porta 1,0 mA come mostrato. Trova la forza su ogni braccio e la forza totale
agire sul loop del campo della linea 2.0 A.

µ oIo1
F12 = ∫ Io dL× B ;
2 2
2π ρ
a 1 B1 = φ

A → B:dL2 = 5dφ a , a × a = 0 φ φ φ

C → D:un × a = 0 φ φ

3
µ oIo1 µ Io Io 3
B → C: F12 = Io∫2d ρ a × a = o 1 2ln az
2π ρ 2π 5
ρ φ

5
 3-27

Fig. P3.41
7
4π x10 H
−
m 3
F12 = ( ) ( Aln
2A10− 3
) un=z − 204pNaz
2π 5
Quindi per B a C:F12= -0,20 nNaz
Allo stesso modo, da D a A:F12 = +0,20 nNaz

P3.42: MATLAB: Modificare MATLAB 3.4 per trovare la forza differenziale che agisce da ciascuno
segmento differenziale individuale sul ciclo. Traccia questa forza rispetto alla posizione phi del
segmento.

%MLP0342
%modifica ML0304 per trovare dF che agisce dal campo
% di ciascun segmento di corrente; grafico vs phi

chiaro
clc
I=1; %corrente in A
a=1; %raggio del ciclo, in m
mu=pi*(4e-7); %permeabilità del vuoto
az=[0 0 1]; %vettore unitario nella direzione z
DL1=a*2*(pi/180)*[0 1 0]; %Assume 2 degree increments
%DL1 è il vettore dell'elemento di prova
%F è l'angolo phi in radianti
%xi,yi è la posizione dell'i-esimo elemento nel ciclo
%Ai & ai = vettore e vettore unitario dall'origine
% a xi,yi
%DLi è il vettore del i-esimo elemento
%Ri1 & ri1 = vettore e vettore unitario dall'i-esimo
% punto da testare

per i=1:179
phi(i)=i*2;
 3-28

F=2*i*pi/180;
xi=a*cos(F);
yi=a*sin(F);
Ai=[xi yi 0];
ai=unitario(Ai);
DLi=(pi*a/90)*incrocia(az,ai);
Ri1=[a-xi -yi 0];
ri1=unitario(Ri1);
num=mu*I*incrocia(DLi,ri1);
den=4*pi*(magvector(Ri1)^2);
B=num/den;
dFvect=I*cross(DL1,B);
dF(i)=dFvect(1);
fine
traccia(phi,dF)
etichetta asse x('angolo in gradi')
etichetta dell'asse y('la forza differenziale, N')

Fig. P3.42

P3.43: MATLAB: Considera un anello conduttore circolare di raggio 4,0 cm nel piano y-z
centrato in (0,6cm,0). Il circuito conduce una corrente di 1,0 mA in senso orario se visto da
+asse x. Una linea di lunghezza infinita sull'asse z conduce 10. Una corrente nell'asse +azdirezione.
Trova la forza netta sul loop.

La seguente routine MATLAB mostra la forza in funzione della posizione radiale attorno a
il ciclo. Nota che mentre c'è una forza netta nella direzione -y, le forze nella direzione z-
direzione annulla.
 3-29

% MLP0343
% trovare la forza totale e il momento su un anello di
% corrente accanto a una linea di corrente

% variabili
% I1,I2 corrente nella linea e nel circuito (A)
% yo centro del ciclo sull'asse y (m)
% uo permeabilità dello spazio libero (H/m)
% N numero di segmenti sul loop
% un raggio di loop (m)
% dalpha elemento di loop differenziale
% dL lunghezza della sezione differenziale
% DL differenza sezione vettore
% B1 I1's vettore di flusso magnetico (Wb/m^2)
% Rv vettore dal centro del ciclo
% al segmento di differenza
% ar vettore unitario per Rv
% y,z la posizione del segmento di differenza

clc
chiaro

% inizializza variabili
I1=10;
I2=1e-3;
yo=.06;
uo=pi*4e-7;
N=180;
a=0.04;
dalpha=360/N;
dL=a*dalpha*pi/180;
ax=[1 0 0];

% eseguire calcoli
per i=1:N
dalpha=360/N; Fig. P3.43
alpha=(i-1)*dalpha;
phi(i)=alpha;
z=a*sin(alpha*pi/180);
y=yo+a*cos(alpha*pi/180);
B1=-(uo*I1/(2*pi*y))*[1 0 0];
Rv=[0 y-.06 z];
 3-30

ar=unitvector(Rv);
aL=cross(ar,ax);
DL=dL*aL;
dF=cross(I2*DL,B1);
dFx(i)=dF(1);
dFy(i)=dF(2);
dFz(i)=dF(3);
fine
traccia(phi,dFy,phi,dFz,'--k')
legenda('dFy','dFz')
Fnet=somma(dFy)

Eseguendo il programma otteniamo:

Fnet = -4.2932e-009
>>
SoFrete= -4.3 nNay

P3.44: MATLAB: Un circuito quadrato con una corrente di 1,0 A e lato di 4,0 cm è centrato sul piano x-y
piano. Assumere un filo con diametro di 1 mm e stimare il vettore di forza su un braccio risultante
dal campo delle altre 3 braccia.

% MLP0344 V2
%
% Anello quadrato di corrente è centrato nel piano x-y. Visto
% dall'asse +z, lascia che la corrente fluisca in senso orario. Vogliamo
% trovare la forza sul braccio a x = +2 cm risultante da
il
% corrente nelle braccia a y = -2 cm, x = -2 cm e y = +2 cm.

% Wentworth, 03/12/03

% Variabili
% lunghezza del lato (m)
% b raggio del filo (m)
% I corrente nel ciclo (A)
% uo permeabilità dello spazio libero (H/m)
% N numero di segmenti per ogni braccio
% xi,yi posizione del segmento del braccio di prova (a x = +2 cm)
% xj,yj posizione del segmento del braccio sorgente (a y = -2 cm)
% xk,yk posizione del segmento dell'arm e sorgente (a x = -2 cm)
% xL,yL posizione del segmento del braccio sorgente (a y = +2 cm)
 3-31

% dLi vettore segmento di test differenziale

% dLj, dLk, dLL vettori di differenza su fonti
% Rji vettore dal punto sorgente j al punto di test i
% aji vettore unitario di Rji

clc;clear;

a=0.04;
b=.0005;
I=1;
uo=pi*4e-7;
N=80;

per i=1:N
xi=(a/2)+b;
yi=-(a/2)+(i-0.5)*a/N;
ypos(i)=yi;
dLi=(a/N)*[0 -1 0];
per j=1:N
xj=-(a/2)+(j-0.5)*a/N;
yj=-a/2-b;
dLj=(a/N)*[-1 0 0];
Rji=[xi-xj yi-yj 0];
aji=unitario(Rji);
num=I*CROSS(dLj,aji);
den=4*pi*(magvector(Rji))^2;
H=num/den;
dHj(j)=H(3);
fine
per k=1:N
yk=(-a/2)+(k-0.5)*a/N;
xk = -(a/2) - b;
dLk=(a/N)*[0 1 0];
Rki=[xi-xk yi-yk 0];
aki=unitario(Rki);
num=I*CROSS(dLk,aki);
den=4*pi*(magvector(Rki))^2;
H=num/den;
dHk(k)=H(3);
fine
per L=1:N
xL = (-a/2) + (L - 0.5) * a / N;
yL=(a/2)+b;
 3-32

dLL=(a/N)*[1 0 0];
RLi=[xi-xL yi-yL 0];
aLi=unitvector(RLi);
num=I*CROSS(dLL,aLi);
den=4*pi*(magvector(RLi))^2;
H=num/den;
dHL(L)=H(3);
fine
H=somma(dHj)+somma(dHk)+somma(dHL);
B=uo*H*[0 0 1];
F=I*INCROCIATE(dLi,B);
dF(i)=F(1);
fine
Ftot=sum(dF)
traccia(ypos,dF)

Eseguendo il programma:

Ftot =

7.4448e-007
Fig. P3.44
Quindi Ftot= 740 nN

P3.45: Un foglio di corrente K= 100axA/m esiste a z = 2.0 cm. Un anello di 2.0 cm di diametro
centrato nel piano x-y a z = 0 conduce una corrente di 1.0 mA in +a φ direzione. Trova il
torsione su questo anello.

τ = m × B;
m = Èa N = ( 10 Un
3 −
(
) π ( 0,01m) 2
)a z = 314x10 9Am2a z
−

1 1 A A
B = µ o H; H = K × aN = 100ax × ( − a z ) = 50ay ;
2 2 m m
Wb
B = 50 µ o a y ;
m2
τ = m × B = -20pNma x

P3.46: 10 spire di filo isolato in una bobina di 4.0 cm di diametro sono centrate nel piano x-y.
Ogni filo della bobina conduce 2,0 A di corrente nellaφ direzione. (a) Qual è la
 3-33

momento dipolare magnetico di questa bobina? Ora supponiamo che questa bobina sia in un campo magnetico uniforme
fieldB= 6.0ax+ 3.0ay+ 6.0azWb/m2, (b) qual è la coppia sulla bobina?

) ( π ( 0,02m) ) = 25,1mAm a
2
(a) m = NISa x = ( 10 ) (2A 2
z

Wb
(b) τ = m × B = 25.1mAaz × ( 6a x + 3ay + 6a z ) = 0,151ay0.075a
− xNm
m2
τ = − 75ax + 151aymNm

P3.47: Un loop quadrato conduttivo di lato 2,0 cm è libero di ruotare attorno a un lato che è
fisso sull'asse z. C'è una corrente di 1,0 A nel circuito, che scorre in -azdirezione su
campo B uniforme esiste in modo tale che quando il circuito è posizionato a φ= 90°, non
la coppia agisce sull'anello, e quando l'anello è posizionato a φ = 180° una coppia massima di
8,0µN-mazsi verifica. Determina la densità del flusso magnetico.

Aφ= 90° = m × B = 0, Bis in direzione di m,

m = Èa N = 0,0004Am2a x. Anche, dato τche
e quindi B = ±Boax.
Atφ= 180°, m = Èa Nun = 0.0004Am2 y, e τ = m × B = 8x10 6a zNm. −

Quindi, B = -BoaxandmBo= 8x10-6, quindi

8x10 6 −
mWb mWb
Bo = = 20 , e B = − 20 2 a x .
0,0004 m 2
m

Materiali Magnetici
P3.48: Un filo di nichel solido con un diametro di 2,0 mm conduce uniformemente 1,0 ampere di corrente.
Determina la grandezza della densità di flusso magnetico B come funzione della distanza radiale
dal centro del filo. Traccia a un raggio di 2 mm.

Io 1A kA
J= az = 2un z= 31.8az 2
πa π ( 1x10 m ) m
2 3 −

∫ HgdL = Io enc= ∫ J gdS

Io Io
2π ρ Hφ = 2∫
ρ d ρ dφ = 2 ρ 2
π un a
Ioρ µ µ Ioρ
per ρ ≤ a H = a ; B= r o2 a
2π a 2
2π a
φ φ

Io µ Io
per ρ ≥ a H = a ; B= o a
2π ρ 2π ρ
φ φ
 3-34

% MLP0348
% generare grafico per il problema ACL

a=2e-3; %raggio del filo solido (m)

I=1; %corrente (A)
N=30; %numero di punti dati da tracciare
rmax=4e-3; %raggio massimo per il grafico (m)
dr=rmax/N;
uo=pi*4e-7;
ur=600;

per i=1:round(a/dr)
r(i)=i*dr;
B(i)=(ur*uo*I/(2*pi*a^2))*r(i);
fine
per i=round(a/dr)+1:N
r(i)=i*dr;
B(i)=uo*I/(2*pi*r(i));
Fine
rmm=r*1000;
trama(rmm,B)
etichetta asse x('rho(cm)')
Etichetta asse Y ('B (Wb/m^2)')
griglia attivata

Fig. P3.48

8. Condizioni al contorno
P3.49: Un'interfaccia piana separa due mezzi magnetici. Il campo magnetico nel mezzo 1
(conµr1) forma un angolo α1con una normale all'interfaccia. (a) Trova un'equazione per α2,
l'angolo il campo in media 2 (che ha µr2) fa con una normale all'interfaccia, in termini
diα1e le permeabilità relative nei due mezzi. (b) Supponiamo che il mezzo 1 sia nichel e
il mezzo 2 è aria, e il campo magnetico nel nichel forma un angolo di 89° con una normale
alla superficie. Trova α2.

H1 = H1Na N + H1T aT;H2TaT = H1TaT;B1Na N = µ r1 µ oH1Na N

µ r1
B2 Na N = B1Na N = µ r1 µ o H1Nun N= µ r 2 µ o H 2 Na N ; H2N = H1N
µ r2
 3-35

H 2T H1T
abbronzare
α2= = ;
H2 N µ r1
H1N
µ r2
µ H1T µ
abbronzare
α 2 = r2 = r 2 abbronzare
α1
µ r1 H1N µ r1

1 µ r2
α 2 = abbronzato abbronzare
−
α1
µ r1
1 1
α 2 tan89o = 5.5o
= abbronzato −

600 Fig. P3.49

P3.50: MATLAB: Supponiamo che il piano z = 0 separi due mezzi magnetici, e che non
La corrente superficiale esiste all'interfaccia. Costruisci un programma che richieda all'utente di inserire µr1
(forz< 0),µr2(forz> 0), e uno dei campi, eitherH1orH2. Il programma è per
calcolare l'ignoto H. Verifica il programma utilizzando l'Esempio 3.11.

% M-File: MLP0350
%
% Data H1 al confine tra una coppia di
% materiali senza corrente superficiale al confine,
% calcolare H2.
%
clc
chiaro
% inserire variabili
inserisci le quantità di vettori tra parentesi,
disp('ad esempio: [1 2 3]')
ur1=input('permeabilità relativa nel materiale 1: ');
ur2=input('permeabilità relativa nel materiale 2: ');
a12=input('vettore unitario da mtrl 1 a mtrl 2: ');
F=input('materiale in cui il campo è noto (1 o 2): ');
Ha=input('vettore di intensità del campo magnetico noto: ');

se F==1
ura=ur1;
urb=ur2;
a=a12;
altrimenti
ura=ur2;
urb=ur1;
 3-36

a=-a12;
fine

% eseguire calcoli
Hna=dot(Ha,a)*a;
Hta=Ha-Hna;
Htb=Hta;
Bna=ura*Hna; %ignora uo poiché verrà fattorizzato
Bnb=Bna;
Hnb=Bnb/urb;
Il campo magnetico nell'altro mezzo è:
Hb=Htb+Hnb

Esegui ora il programma (per l'Esempio 3.11):

inserisci le quantità dei vettori tra parentesi
[1 2 3]
relative permeability in material 1: 6000
relative permeability in material 2: 3000
vettore unitario da mtrl 1 a mtrl 2: [0 0 1]
material where field is known (1 or 2): 1
vettore di intensità del campo magnetico noto: [6 2 3]

ans =
Il campo magnetico nell'altro mezzo è:
Hb = 6 2 6

Per un secondo test, eseguire il programma per il problema P3.52(a).

inserire le quantità dei vettori tra parentesi

[1 2 3]
permeabilità relativa nel materiale 1: 4
permeabilità relativa nel materiale 2: 1
vettore unitario da mtrl 1 a mtrl 2: [0 0 -1]
materiale dove il campo è conosciuto (1 o 2): 1
vettore di intensità del campo magnetico noto: [3 0 4]

ans =
Il campo magnetico nell'altro mezzo è:
Hb = 3 0 16
 3-37

P3.51: Il piano y = 0 separa due mezzi magnetici. Il mezzo 1 (y < 0) ha µr1= 3.0 e
media 2 (y> 0) ha µr2= 9.0. Un foglio currentK= (1/µo)axA/m esiste all'interfaccia,
eB1= 4.0ay+ 6.0azWb/m2. (a) TrovaB2. (b) Quali angoli faB1eB2fai con un
normale rispetto alla superficie?

(a) B N1 = 4a y
(b) B N2 = B N1 = 4ay
(c) BT1 = 6az
BT 1 2
(d) HT1 = = un z
µ r1 µ o µ o
3
(e) HT2 = a z(vedi sotto)
µo
(f) BT2 = µ r µ o H T2 = 27az
Wb Fig. P3.51
(g) B 2 = 4a y + 27az 2
m
Ora per il passaggio (e):
1
a 21 × ( H1 − H 2 ) = K ; − a y × ( H T 1a z − HT2a z ) = un x
µo
1 1 1 2 3
− ( H T1 − H T2 ) a x = a x;HT2 − HT 1 = ;HT2 = + =
µo µo µo µo µo

1 BT 1 1 BT2
Angles: α 1 = abbronzare = 56o ; α 2 = abbronzarsi = 82o
− −

BN1 BN2

P3.52: Sopra il piano x-y (z> 0), esiste un materiale magnetico con µ r1= 4.0 e un
campoH1= 3.0ax+ 4.0azA/m. Sotto il piano (z< 0) c'è spazio libero. (a) Trova H 2,

assumendo che il confine sia privo di corrente superficiale. Quale angolo fa H2fai con una
normale alla superficie? (b) Trova H2, assumendo che il confine abbia una corrente superficiale K = 5.0
unxA/m.

(a)
B N2
(1) H N1 = 4a z, (2) B N1 = 16 µ o un
z, (3) B N2 = B N1 = 16µo a z,
H N2 = = 16az
µo
A
(5) HT1 = 3a x(6) HT2 = HT 1 = 3a x, (7) H 2 = 3ax + 16az
m
 3-38

1 H T2
α 2 = abbronzare = 10.6o
−

H N2
(b)
Ora il passo (6) diventa a 21 × ( H1 − H 2 ) = K , sherea21=az.

Lasciamo H
che2 = Aa x + Ba y16a
+ z, poi a 21 × ( 3a x + 4az − Aa x − Ba y − 16az ) = 5ax
Risolvi per A e B:
ax a y az
0 0 1 = Ba x + ( A− 3) a y = 5a x; so A = 3 and B = 5
3 − A − B − 12

A
H 2 = 3ax + 5a y16a
Finalmente, + z .
m

P3.53: Il piano x-z separa il materiale magnetico con µr1= 2.0 (per y < 0) da magnetico
materiale conµr2= 4.0 (se y> 0). Nel mezzo 1, c'è un campo H1= 2.0ax+ 4.0ay+ 6.0az
A/m. Trova H2supponendo che il confine abbia una corrente superficiale K = 2.0ax- 2.0azA/m.

8µ o
(1)H1N = 4a y(2)B1N = 8µ oun y(3)B 2 N= B1N = 8µ oun y(4)H2 N = a y = 2ay
4µ o
(5)H1T = 2a x + 6un z,(6)a21× ( H1 − H 2 ) = K, lasciatHT2 = H xa x + H z un, z
× (y ( 2 − H x ) a x + ( 6 − H z ) un
− una
quindi z ) = 2un−x 2un
,z
ax ay az
0 −1 0 = ( H z − 6 ) a x + ( 2 − H x ) a z = 2un−x 2un z, soHz = 8,Hx = 4
2 − H x0 6 − Hz
A
(7)H2Ta=x 4 + 8a z , H 2 = 4a x + 2a y + 8az
m

P3.54: Una linea di lunghezza infinita di corrente 2πA nella direzione +azla
direzione esiste sull'asse z.
Questo è circondato dall'aria per ρ≤ 50 cm, a quel punto il mezzo magnetico ha µr2= 9,0
per ρ > 50 cm. Se il campo in mezzo 2 a ρ = 1,0 m è H = 5,0 aφ A/m, trova il foglio
il vettore densità di corrente a ρ = 50 cm, se esistente.

Metodo 1:
Io1
Dalla sola linea di corrente avremmo H1 =
2π ( 1) una= 1a .
φ
φ
 3-39

Ora, dato cheHTOT = 5a = H1 + H 2, poiH2 φ

= 4a φ è il contributo dal foglio
attuale.
Io 2 Io 2 8π A
H2 =
2π ( 1) a φ
= 4a , quindi Io2= 8π A, poi K =
φ
az =
2π una 2π ( 0,5)
= 8 az
m
A
K= 8 az
m

Metodo 2:
Da me1al confine abbiamo
2π a
H1 = = 2a , ma 5a a ρ = 1,0m corrisponde a 10a a ρ =0.5m poiché H varia
φ

2π ( 0,5 )
φ φ φ

come 1/ρ. Quindi

− a × ( 2a − 10a
ρ φ φ ) = K, − a × ( − 8a
ρ φ ) = 8a z = K
A
K= 8 una z
m

9. Induttanza ed Energia Magnetica

P3.55: Considera una lunga coppia di fili paralleli dritti, ciascuno di raggio a, con i loro centri
separati da una distanza. Supponendo che d>>a, trova l'induttanza per unità di lunghezza per questo
coppia di fili.

Ogni filo può essere risolto usando il

approccio energetico
L µ µ
= 2 =
h fili 8π 4π
I campi non sono confini a un
volume, quindi dobbiamo utilizzare il flusso
collegamentoapproccio a trova
induttanza al di fuori dei fili.
Io Io
H1 = a , H2 = a ,
2π ρ 2π ( d −ρ )
φ φ

φ1 = ∫ BgdS , φ TOT = 2∫ BgdS

Fig. P3.55
d h
µ Io 2µ Io − ad ρ µ Io d − a
φ TOT = 2 ∫
2π ∫un ρ ∫ dz =
ungd ρ dza = ln h
2π ρ π a
φ φ

φ TOT L µ 1 d− a
L=
Finalmente, con , arriviamo a = + ln .
Io hTOT π 4 un
 3-40

P3.56: Nel problema P3.23, il compito era trovare il campo al centro (raggio b) di un N-giri.
toroide. Se il raggio del toroide è grande rispetto al diametro della bobina (cioè, se
b>>c-a), quindi il campo è approssimativamente costante da raggio a raggio c. (a) Ottieni
un'espressione per l'induttanza del toroide. (b) Trova L se ci sono 600 spire attorno a un
Nucleo di ferro 99,8% con a = 8,0 cm e c = 9,0 cm.

NI
Da P3.23 abbiamo trovato:H = a al raggio per il toro.
2π b
φ

µ NI
φ = ∫ Bgd S = ∫ 2π b
ga
a dS
φ φ

µ NI 2
( b − a.)
2
e con area π ( b − un) abbiamo φ =
2b
µ N 2 Io 2
λ = Nφ = (b − a )
2b
λ µ N2 2
L= = (b− a)
Io 2b
2
( 5000 ) ( 4π )( 2
7
x10 600
−
) (
.085 − .080 )
Inserendo i nostri numeri abbiamo: L = = 0,33H
2 ( 0,085 )

P3.57: MATLAB: Considera che un filo solido di raggio a = 1,0 mm sia piegato in un cerchio.
anello di raggio 10 cm. Trascurando l'induttanza interna del filo, scrivi un programma per trovare
l'induttanza per questo circuito.

% M-File: MLP0357
%
% Induttanza all'interno di un circuito conduttivo
%
% Questo modifica ML0302 per calcolare l'induttanza di
% un circuito conduttivo. Lo fa calcolando il
% campo magnetico in punti discreti lungo un settore di torta,
% calcola quindi il flusso attraverso ciascuna parte del
% fetta. Poi moltiplica per il numero di fette
% nell' 'torta'.
%
% Wentworth, 19/01/03
%
% Variabili:
% I corrente(A) in direzione +phi sul anello
 3-41

% a raggio dell'anello (m)

% b raggio del filo (m)
% Ndeg numero di incrementi per phi
% f angolo di phi in radianti
% df cambiamento differenziale in phi
% dL vettore di lunghezza differenziale sull'anello
% dLmag magnitudine di dL
% dLuv vettore unitario nella direzione di dL
[xL,yL,0] posizione del punto sorgente
% Ntest numero di punti di test
% Rsuv vettore unitario dall'origine al punto sorgente
% R vettore dalla sorgente al punto di prova
% Ruv vettore unitario per R
% Rmag magnitude of R
% dH contributo differenziale a H
% dHmag magnitudine di dH
% raggio distanza radiale dall'origine
% Hz campo magnetico totale al punto di test
% Bz densità di flusso magnetico totale al punto di prova
% flusso flusso attraverso ciascun segmento differenziale

clc %pulisce la finestra dei comandi

cancella %cancella le variabili

% Inizializza Variabili
a=0.1;
b=1e-3;
I=1;
Ndeg=180;
Ntest=60;
uo=pi*4e-7;
df=360/Ndeg;
dLmag=(df*pi/180)*a;
dr=(a-b)/Ntest;

% Calcola il flusso attraverso ciascun segmento del cuneo di torta

per j=1:Ntest
x=(j-0.5)*dr;
per i=(df/2):df:360
f=i*pi/180;
xL=a*cos(f);
yL=a*sin(f);
Rsuv=[xL yL 0]/a;
 3-42

dLuv=cross([0 0 1],Rsuv);
dL=dLmag*dLuv;
R=[x-xL -yL 0];
Rmag=magvector(R);
Ruv=R/Rmag;
dH=I*croce(dL,Ruv)/(4*pi*Rmag^2);
dHmag(i)=magvector(dH);
fine
Hz(j)=sum(dHmag);
Bz(j)=uo*Hz(j);
dSz(j)=x*df*(pi/180)*dr;
flux(j)=Bz(j)*dSz(j);
fine

fluxwedge = somma(flux);
Inductance=Ndeg*fluxwedge

Ora esegui il programma:

Inductance = 5.5410e-007
o
L = 550 nH

P3.58: Trova l'induttanza mutua tra un filo infinito e un filo rettangolare

con le dimensioni mostrate nella Figura 3.58.

Io1 µ Io
H1 = un, B1 = o 1a
2π ρ 2π ρ
φ φ

ρ + b
µ oIo1uno d ρ µ oIo1b ρ o + a
φ 12 = ∫ B1 gdS 2 =
2π ρ∫o ρ ∫ dz =
0
2π
ln
ρo
φ 12 µ ob ρ + a
M 12 = = ln o
Io1 2π ρo

P3.59: Considera una coppia di anelli conduttori concentrici, centrati nella stessa pianura, con
radiiaandb. Determina l'induttanza mutua tra questi anelli se b >> a.

In questo caso si guiderà il loop b-radius con corrente I. Qui, al centro del b-radius
ciclo che abbiamo dall'Eqn. (3.10):
 3-43

µ oIo
B1 = az
2b

Poi, φ 12 = ∫
µ Io 2
B1 gdS 2 = o 1∫ ρ d ρ dφ = π µ o 1
ρ Io 2 a
= π µoa Io
2
1
2b 2b2 ρ =0
2b

M 12 =
Quindi
φ 12 π µ
= o un 2
Io1 2b

P3.60: Un filo di nichel solido con un diametro di 4,0 cm, centrato sull'asse z, conduce corrente con
una densità J= 4ρA/cm2az(dove ρ è in cm). Trova l'induttanza interna per unità di lunghezza
per il filo con questa distribuzione di corrente.

ρ 2π
8π
Ñ
∫ HgdL = Io enc, 2π ρ H= 4∫ ρ 2d ρ φ ∫ dφ =
3
ρ 3

0 0
3
8π ρ 4 4µ 2
H = = ρ 2, B = ρ
6π ρ 3 3
φ φ

a 2π h
1 4µ ρ2 4 ρ 2 1 16 8
Wm = ∫ a g a ρ d ρ dφ dz = µ
∫ρ dρ
5
∫ dφ ∫ dz = µ a π h
6

2 3 3 2 9 27
φ φ

0 0 0

1 2
Wm = LI ;
2
8π a 3
Ora risolviamo per I: Io= ∫ J gdS =
3
, e poi risolvi per L:

1 µh
L=
12 π

7
L 1 µ 1 4π x10 −
µ H
= L= = 600 = 20 .
h 12 π 12 π m

10. Circuiti Magnetici

P3.61: Riferendosi alla Figura 3.48(a), supponiamo che 2.0 Ampere fluiscano attraverso 80 spire di un toroide.
che ha un'area della sezione trasversale del nucleo di 2,0 cm2e
un raggio medio di 80 cm. Il nucleo è
{"text":{"a":"Quanto flusso magnetico esiste nel toroide?","b":"Quanta energia"}}
è immagazzinato nel campo magnetico contenuto dal toroide?

µ NI ( 5000 ) ( 4π 7
x10 80) ( 2) (
−
)
(a) φ = BA= A;B = = 0.200 Wb , φ = 40µWb
2π ρ o 2π ( 0,8 ) m2
 3-44

1 1B2
(b) Wm = BH ∫ dv = 2π ρoA = 3,2 mJ
2 2 µ

P3.62: Nella Figura 3.59, un nucleo toroidale di diametro 2,0 cm con µr110.000 è avvolto con
uno strato spesso 1,0 cm di µr2= 3000. Il toroide ha un raggio medio di 1,0 m. Per 20 A di
corrente condotta attraverso 50 anelli di filo, trova l'intensità del campo magnetico in ogni materiale
del toroide.

Questo toroidale può essere modellato con il circuito della Figura P3.62.

Fig. P3.62

Vm= NI = (50)(20A)=1000 A-giri

l 2π ρ o 2π ρo
R= ; R1 = = 1.592x106 ; R= = 1,768x106 ;
µ r2 µ oπ ( b − a )
2
µA µ r1 µ oπ a 2 2 2

Vm 6 628x10 6 −
B A
φ1 = = 628x10 = BA;B =
−
2
= 2;H1 = = 159
R1 π ( 0,01) µ r1 µ o m
Vm 566x10 6 −

φ2 = = 566x10 6 = BA;B =
−
= 0.6001
R2 π ( ( 0,02) 2
− ( 0.01)
2
)
B A
H2 = = 159
µ r2 µ o m

P3.63: Riferendosi alla Figura 3.49(a), il nucleo di 2,0 cm di diametro è caratterizzato da

curva di magnetizzazione della Figura 3.60. Il toroide ha un raggio medio di 60 cm. Per 10 A
 3-45

di corrente attraversata da 100 spire di filo, trova l'intensità del campo magnetico in 1.0
mm intervallo.

Riferendosi al modello della Figura 3.49(b),

Vm= NI = (100)(10A)=1000 A-giri; A =πa2= 314x10-6
lG
RG= = 2,53x106 ;
µ oA
Approccio:
(1) Assumere B: B = 0,4
(2) Leggi H dal grafico: H = 200
2π ρ o 12x10 3H
(3) Valuta R C = = :Rc=6x106
µA B
(4)RTOT=Rc+ 2,53x106:RTOT=8,53x106
(5)φ= Vm/RTOT:φ=117x10-6
(6) B =φ/A: B = 0.373

Inserisci questo valore di B nel passo (1) e ripeti. Dopo 4 o 5 iterazioni, la routine
converges to a solution B = 0.382 Wb/m2e H = 190 kA/m. Quindi, nel gap, abbiamo
B kA
H= = 300 .
µo m

P3.64: Riferendosi alla Figura 3.52, supponiamo che l'area della sezione trasversale della barra sia di 3.0 cm2,
mentre quello del nucleo dell'elettromagnete è di 2,0 cm2Inoltre, la barra ha una permittività relativa
of 3000, while that of the magnetic core is 10,000. The dimensions forhandware 12.
cm e 16. cm, rispettivamente. Se la massa della barra è 20. kg, quanta corrente deve essere
guidato attraverso 24 anelli per sostenere il bar contro la gravità?

Seguirà un approccio simile all'Esempio 3.18, ma ci rendiamo conto che Ac≠ Ab.

2
l 2(.12) + 0,16 100cm
R=
c = = 159x103
µ Ac 10.000 µ o(2cm2) m
2
l 0.16 100 cm
R=
b = = 141x103
µ Ab 3.000 µ o(3cm )
2
m
Vm NI NI
Vm = NI φ = = = µ oHAc; H=
RTOTRTOT RTOT µ oAc
2 2
NI NI 1
F = µ oH A c = µ o
2
Ac = = mg
RTOTALE
µ oAc RTOT µ oAc
 3-46

RTOT
I=
Risolvere per I: mgµ oUn
c = 2.8A
N

P3.65: Considera uno spazio d'aria di 1.0 mm nella Figura 3.49(a). Il raggio medio del toroide e la sezione trasversale

l'area sezionale è di 50 cm e 2,0 cm2, rispettivamente. Se il nucleo magnetico ha aµr= 6000,

e 4,0 A vengono fatti passare attraverso 30 spire, determina la forza che tira il gap chiuso.

6(
2π ρ o l A) ( svolte)
F = µ oH2ARTOT= + = 6,06x10
µ r µ oA µ oA Wb

φ =
Vm
=
( 30 ) ( 4 ) = 19,79x10 6 = BA −

RTOT6.06x106
Nello spazio, B= 19.79x10-6/A = 98,97x10-3=µoH,H= 78,76x103A/m
2
1
F = µ oH A = ( 4π x10 78,76
) ( x 10 )
2
2 − 7 3
( 2) = 1.6
100
F = 1.6N