University of Pennsylvania

Department of Electrical Engineering and Computer Science
EEL6482 Electromagnetic Theory
Homework/Project Set No. 2
Problem 7.
A line current I flows along an elliptical ring located in the xoy plane. If we put thumb in
the z direction, current will follow right-hand rule to flow. The elliptical equation is
๐‘ฅ2
๐‘ฆ2
given as ๐‘Ž2 + ๐‘2 = 1 . Let a = 1, b = 0.75, ๐ผ = 107 A. Find magnetic flux density B on
yoz plane defined by [y,z] = meshgrid(-2:0.5:2, -2:0.5:2) and plot using quiver. Discuss
whatever physics you can find through investigation.
(Note: let ๐‘ฅ โ€ฒ = acos(๐œ™ โ€ฒ ) , ๐‘ฆ โ€ฒ = ๐‘ sin(๐œ™ โ€ฒ ) so that ๐๐ฅโ€ฒ = ๐‘‘๐‘ฅ โ€ฒ ๐ฎ๐ฑ + ๐‘‘๐‘ฆโ€ฒ๐ฎ๐ฒ )
Problem 8.
The following figure shows a simple linear actuator. By controlling the DC voltage
source, we can control the movement of the bar. Assuming the mass of the bar is M, the
DC magnetic field is B, the length of the bar is l, and the electrical loss of the system is
described using resistor R. The force from the load is ๐…load.
(1) Before ๐‘ก = 0,, the bar is stationary. There must be a current going through the bar to
get a force ๐…ind = ๐…load . Find the DC voltage ๐‘‰๐ต1 of this case.
(2) At ๐‘ก = 0, apply ๐‘‰๐ต2 > ๐‘‰๐ต1 so that the bar will start to accelerate. Find the the
dynamical equation for transient velocity ๐‘ฃ(๐‘ก)and solve using MatLab function
dsolve.
(3) At ๐‘ก = ๐‘กโ€ฒ, apply ๐‘‰๐ต3 < ๐‘‰๐ต1 (can be negative) so that the bar will start to decelerate.
Find the the dynamical equation for transient velocity ๐‘ฃ(๐‘ก)and solve using MatLab
function dsolve.
(4) When the velocity decelerates to 0, suddenly apply ๐‘‰๐ต1so that the bar will be held at
a new position. Find the total actuation time from the old to new position. Also find
the stroke (distance from old to new position) by integrating the velocity.
(5) If M = 1 kg, B = 0.2 T, l = 10 m, ๐‘‰๐ต2 = 120V, ๐‘‰๐ต3 = โˆ’120๐‘‰, R = 2โ„ฆ, ๐นload =
30 N, and stroke s = 1 m. Find ๐‘‰๐ต1, actuation time ๐‘‡๐‘Ž๐‘๐‘ก and ๐‘กโ€ฒ. Plot velocity ๐‘ฃ(๐‘ก)
and current ๐‘–(๐‘ก)versus time from 0 to ๐‘‡๐‘Ž๐‘๐‘ก .
(6) There is another way to drive the device. Instead of changing the DC voltage source,
we can change the magnetic field B. If this is your research project, discuss how you
are going to do.
Problem 9.
A two legged core with all the dimensions is shown in the following figure. The windings
on the left leg of the core has 600 turns, and the windings on the right has 200 turns. The
coils are wound in the directions shown in the figure. Neglect magnetic leakage.
Assume the relative permeability of the core is 1200 (constant).
(1) Find reluctance of each leg (left, right, top and bottom);
(2) Find the self inductances of these two windings and mutual inductance between them;
(3) Find the flux produced by currents i1 ๏€ฝ 0.5A and i2 ๏€ฝ 1A .
(4) If we do not assume a constant permeability for the core and the material will be M5
steel, find the flux produced by currents i1 ๏€ฝ 0.5A and i2 ๏€ฝ 1A . The BH relationship
of M5 steel (M5DC_B2H.m) was given in class. (Hint: You need to solve the
transcendental equation:
๏†
๏†
๏†
๏†
H(
)lleft ๏€ซ H (
)lright ๏€ซ H (
)ltop ๏€ซ H (
)lbottom ๏€ฝ N1i1 ๏€ซ N 2i2 to find ๏† .
Aleft
Aright
Atop
Abottom
MatLab function fsolve can help to solve transcendental equation. )
Problem 10.
The time harmonic magnetic field inside a source-free perfect magnetic conducting
(PMC) rectangular waveguide is given by
๏ฐ
H ๏€ฝ H 0 sin( y )e ๏€ญ jkz z a x , 0 ๏‚ฃ x ๏‚ฃ a, 0 ๏‚ฃ y ๏‚ฃ b
b
Assume that the media inside the waveguide is free space, and four side walls at x ๏€ฝ 0, a
and y ๏€ฝ 0, b are perfect magnetic conductor. For a section of the waveguide of length d
along z axis from z ๏€ฝ 0 to z ๏€ฝ d , find:
(1) time domain expression of electric fields;
(2) k z ;
(3) the exiting complex power;
(4) the supplied complex power;
(5) the dissipated real power;
(6) time average electric energy;
(7) time average magnetic energy;
(8) verify the conservation of energy equation for this set of fields.
Problem 11.
The phasor form expression for the electric field of a uniform plane wave propagating in
lossless isotropic medium is
E(x) = 60๏ฐ๏›๏€จ๏€ฑ๏€ซj)ay -๏€ (1-j)az]exp(j4x) V/m.
The magnetic field at point (0,0,0) when t = 0 is 2 A/m. The unit of t is second and the
unit of x is meter. The radian frequency is 3๏‚ด 10 8 rad/s. Find:
(1) the time domain expression of the electric field;
(2) the wave impedance ๏จ๏€ of the media;
(3) the relative dielectric constant ๏ฅr and relative permeability ๏ญr;
(4) the time domain expression for magnetic field;
(5) the time-average power density.
Problem 12.
The magnetic field intensity of a plane wave traveling in a lossy earth is given by
H ๏€ฝ (a y ๏€ซ ja z ) H 0 e ๏€ญ๏กx e ๏€ญ j๏ขx
where H 0 ๏€ฝ 1 ฮผA/m. Assuming the lossy earth has a conductivity of 10 ๏€ญ4 S/m, a
dielectric constant of 9, and the frequency of operation is 1 GHz, find inside the earth:
(1) the corresponding time domain electric field;
(2) the average power density;
(3) the phase constant (rad/m);
(4) the phase velocity (m/s);
(5) the wavelength;
(6) the attenuation constant (Np/m);
(7) the skin depth.
Problem 13.
For sea water, assuming ๏ฅ r ๏€ฝ 81, ๏ญ r ๏€ฝ 1, ๏ณ ๏€ฝ 4S/m, and f ๏€ฝ 10 4 GHz, find:
(1) the complex propagation constant;
(2) the phase velocity (m/s);
(3) the wavelength;
(4) the attenuation constant (Np/m);
(5) the skin depth.