1. Introduction
Over the last decade, interest in wireless power transfer has grown. The main contributing factor to this was increased demand for wireless charging of electric vehicles and portable electronics [
1,
2,
3]. The most frequently used method for wireless power transfer is inductive wireless power transfer, or IPT. It has several advantages over other wireless power transfer methods, including safety, a wide power transfer range, scalability, and high efficiency at short and medium ranges. The power range of an IPT is usually between a few watts and several kilowatts, depending on the application.
The most important parts of an IPT system are the transmitter and receiver coils, which are separated by an air gap that can range from a few millimeters to a few decimeters. The transmitter coil generates an alternating magnetic field that induces a voltage in the receiver coil and supplies power to the load. In a wireless charging application, the load is usually a battery. The most common type of coil used in IPT is a planar spiral coil with a circular, quadrature, or square footprint. It can be used in different applications, from mobile phones to electric car chargers. The dimensions of the planar spiral coil are usually determined by the application constraint, such as the available area, the air gap between the coils, and the thickness of the wire. The inductance of the coil is determined by the maximum transmitted power. The number of turns is determined by the inductance of the coil. The coil can be wound with high-frequency litz wire or printed on a circuit board. In high-frequency, high-power applications, the litz wire is necessary to reduce the AC resistance of the coil, and thus reduce the thermal losses of the IPT system.
To better shape the radiated magnetic field, and to prevent organic and metallic foreign matter from interacting with the magnetic field, ferrite material is placed under the transmitter and on the receiver coil. Ferrite material ranges in shape from whole ferrite plates or blocks to ferrite rods. Ferrite plates are generally used in coils with smaller areas, while ferrite rods are used in larger coil areas. Ferrite rods have poorer shielding properties than ferrite plates but are less expensive and lower in weight than large ferrite plates.
The shape and direction of the magnetic field radiated by the transmitter coil is determined by the shape of the coil, the number of turns, the diameter of the wire, and other physical parameters [
4,
5,
6,
7,
8,
9]. The strength of the transmitter coil field is determined by the physical parameters of the coil and the primary excitation current. The calculation of the magnetic field and the optimization of the coil structure are usually preformed with an EM field simulation software, such as Ansys Maxwell or COMSOL Multiphysics. The coil parameters can be optimized numerically and then used in the fabrication process of the IPT coils. The numerical calculations are demanding on time and computer resources. The practical measurements of magnetic flux density are not very common due to their inconvenience. In practical tests, the magnetic field is usually measured at specific, discrete points [
10,
11,
12], to verify that the EM radiation emitted by the IPT system complies with safety guidelines and regulations [
13,
14].
This study aimed at presenting a measurement platform for automated magnetic flux density measurement in the space around the transmitter coil. The measurements were based on a search coil magnetometer mounted onto a 3D positioning mechanism. In the proposed platform, the positioning mechanism moves the search coil around the transmitter coil, and measures the field radiated from the transmitter coil. The system scans the magnetic field around the coil and records the data for visualization. The custom developed app communicates with the 3D positioning mechanism and the search coil measurement circuit. The app is used to determine the scan area, as well as for receiving, recalculating, and logging data. Data are visualized using a 3D magnetic flux density diagram, similar to the graphical representation of magnetic flux density in the EM simulation software. The measurement data can be used to test and develop new IPT transmitter coil topologies and to test the fabricated IPT transmitter coils. This could be especially useful if the fabricated coil were to have a more complex structure, which could not be easily modelled in the EM simulation software. The automated measurement of the magnetic field could also help to test the integration of the IPT transmitter coil in systems, such as mobile devices and laptops. An additional advantage of the measurement platform is that it could serve as a replacement for time-consuming computer simulations that require expensive programs and computer hardware.
The magnetic field measurement system was based on the search coil magnetometer, which works on the principle of electromagnetic induction. The measurement system enables flexibility by choosing sense coils with different radii and number of turns, depending on the required resolution and magnetic field strength. In the case of transmitter coils with a smaller radius, the resolution of the magnetic field measurement should be larger. Likewise, in the case of larger transmitter coils, the resolution of the magnetic field measurement should be smaller. The search coil radius also impacts the minimum vertical z distance between the transmitter coil and the measurement search coil. The 3D search coil enables the measurement of all three magnetic flux density components in the point in space. On the other hand, if a sensor were composed of three sensors based on the hall effect, additional corrections to the measurements would be required in order to compensate for the position offset of the hall sensors. The measurement system was also designed with the replicability of the 3D search coil in mind. New 3D search coils are easy to fabricate; thus, measurements can be adapted to specific requirements. On the other hand, the development and fabrication of MEMS sensors require special manufacturing equipment and processes.
The paper is organized as follows.
Section 2 described the magnetic flux measurement platform and was divided into three main parts. In the first part, we described the transmitter part of the IPT system, which drives the transmitter coil representing the device under test. In the second part of the section, we described the 3D positioning mechanism with mounted search coil and transmitter coil. The third part focused on the measurement circuit used to measure the voltage induced in the 3D search coil. In
Section 3, we described the theory behind the search coil magnetometer.
Section 4 contains the results of measuring the magnetic flux density of two different coil topologies. The measured results were compared with results obtained from the EM field simulation software. In
Section 5, a discussion of the measured and simulated results was presented. Lastly,
Section 6 served as a conclusion of the measurement circuit presented in this paper.
2. Measurement Platform Description
A measurement platform was fabricated to measure the magnetic flux density of the transmitter coil under test in the 3D space. The proposed magnetic flux density measurement platform consisted of three interconnected systems:
A transmitter coil driving circuit;
A search coil with the positioning mechanism;
A search coil voltage measurement circuit.
The measurement platform was designed to work as follows: first, the transmitter coil driving circuit is connected to the transmitter coil and drives the coil with a sinusoidal input current. Then, the magnetic field generated by the transmitter coil is then picked up by a 3D search coil, which can be positioned in x, y, and z directions in the scan area around the coil using the 3D positioning mechanism. The 3D search coil is connected to a 3-channel high impedance voltage measurement circuit, which measures the amplitude of the induced voltage and transmits the value to the PC. The PC logs the data and controls the position of the search coil. The voltage data are converted into a *.csv file. The transmitter coil driving circuit, the transmitter coils, and the search coil were designed with replaceability in mind in case different coil or converter topologies were used.
The fabricated platform is presented in
Figure 1. The main parts of the system included a high-frequency inverter, the transmitter coil and 3D search coil, and the voltage measurement system. The parts were mounted on the 3D positioning mechanism. The platform was connected to the PC.
2.1. Transmitter Coil Driving Circuit
A classic full-bridge high-frequency inverter was used to drive the transmission coil. The inverter converted the input DC voltage to a high-frequency AC voltage, which was used to drive the transmitter coil. To drive the transmitter coil with a sinusoidal current, the coil must be connected to a resonant compensation network. The main purpose of the compensation network was to compensate for the reactive power of the system’s IPT system. A compensation network usually consists of passive components, such as capacitors and inductors. In the system shown in
Figure 1, the series (S) compensation scheme was used on the transmitter side. The transmitter coil was connected in series to a compensation capacitor with value
CT. The compensation capacitor was selected based on the resonant frequency of the IPT system. The resonant frequency of the system (inductor
LT, capacitor
CT, and coil resistance
RT) was adjusted to the high-frequency inverter switching frequency by designing the capacitor
CT as follows:
where
fs is the switching frequency of the high-frequency inverter. The circuit of the inverter, compensation circuit, and transmitter coil are presented in
Figure 1. The inverter consisted of a DC voltage source with voltage
Uin and four MOSFET transistors in full-bridge topology. The inverter was connected to the compensation capacitor
CT and transmitter coil with inductance
LT. The DC resistance of the coil and wires was presented with resistor
RT.
The current and voltage at the input of the resonant IPT circuit are presented in
Figure 2a,b.
Figure 2a presents the theoretical waveforms. Each leg of the full-bridge inverter generated a square voltage with a 50% duty cycle. To control the resonator current
i1, a phase shift was introduced between the half-bridge legs. The phase shift angle could be between 0 and π. In the figure, input voltage
VA −
VB is marked with a red line, while the main harmonic component that drives the current through the circuit is marked with the dashed red line
u1. The resonant circuit served as a bandpass filter. The current through the resonant circuit
i1 is marked with the blue line.
Figure 2b presents the waveforms, measured on a real system. The red line represents the voltage at the input of the resonant circuit, and the blue line represents the sinusoidal current at the input of the resonant circuit.
The first harmonic component of the output voltage of the inverter was calculated from inverter input voltage
Uin and the phase-shift angle between the half-bridge legs of the inverter [
2].
The series compensated transmitter resonator circuit can be described with the voltage equation:
where
U1 is the amplitude of the first harmonic component of the inverter output voltage.
I1 is the amplitude of the resonator current.
ULT is the amplitude of the inductor voltage, and
UCT is the amplitude of the capacitor voltage. The voltage
UM was induced due to the currents in the 3D search coil. Because the search coil was connected to the high impedance measurement circuit, the current through each of its windings was very small. Therefore, the voltage
UM induced by the search coil was almost zero and could be neglected. The
UM stayed at zero regardless of the position of the search coil in the space.
By replacing the capacitor voltage and the inductor voltage with primary current and impedances, Equation (3) now takes form as:
where ω is the angular frequency of the high-frequency inverter. The Equation (4) was derived for the harmonically excited transmitter resonant circuit. The input impedance of the circuit can be expressed from (4) as:
In the case of resonance, the impedance of the capacitor and transmitter coil are the same and can be described with the condition:
Under the condition of resonance, the resistance of the transmitter circuit is defined only by the resistance of the circuit:
If the transmitter coil does not transfer power to the receiver coil, the input current is only limited by the resistance of the primary side. This includes wire resistance on the resistance of the transistors, resistance of the compensating capacitor, and, finally, the resistance of the transmitter coil. In the resonance, the current through the unloaded transmitter coil is relatively high due to the low resonator resistance.
2.2. Fabricated Search Coils and Positioning Mechanism
The mechanical part of the magnetic flux density measurement platform consisted of a 3D positioning mechanism and a 3D search coil magnetometer. The 3D positioning mechanism was fabricated from an open-source 3D printer, Velleman K8200. The 3D printer was stripped of its printing nozzle and heated bed and accompanying thermistors. The metallic positioning platforms were replaced with platforms fabricated from a 1 cm thick acrylic plate. This was a necessity because the magnetic field from the IPT coil could interact with metal objects in the nearby space, which could alter the magnetic flux density and would present additional losses.
The stepper motors of the 3D printer mechanisms were controlled with a driver board based on Arduino UNO. The driver board was connected to the computer with a USB cable. The 3D printer driver communicated with the PC using G-code commands sent via a virtual serial COM port. The commands, important for the interaction between the PC and 3D positioning platforms, were homing and movement to an absolute or relative position in the 3D space.
The 3D search coil was mounted on the top platform of the 3D positioning mechanism. The base of the 3D coil was fabricated from plastic using the 3D printer. Each of the single-axis 3D coils had the same dimensions and the number of turns. The radius of each coil was 5.5 mm and number of turns N was 30. The 3D search coil and the transmitter coil positioned in the 3D space are presented in
Figure 3. The bottom platform with the transmitter coil moves in both
x and
y directions, and the top platform with the 3D search coil moves in the
z-direction.
The parameters of the fabricated search coil used by the magnetometer are presented in
Table 1. The self-inductance of the coils was not required as the system did not calculate mutual inductance. Only the number of turns and surface of each coil were required to calculate magnetic flux density from the induced voltage. The diameter of the search coil was 11 mm. This resulted in a small surface area of the coil itself, allowing more accurate measurement of the magnetic flux density at a point in space. The number of turns was chosen so that even a small magnetic flux density could induce a sufficiently large voltage on the search coil.
2.3. Fabricated Search Coils and Positioning Mechanism
If the 3D search coil was placed in the varying magnetic field, the voltage was induced across the terminals. This voltage was then measured and converted to a digital value with the analogue–digital converter of the microprocessor. The induced voltage was sinusoidal, due to sinusoidal magnetic flux density. The measurement circuit converted the open circuit AC voltage into a DC voltage, which was then sampled using an AD converter. The value of the DC voltage was then used to calculate the amplitude of the input voltage, and, later, the magnetic flux density, based on the coil’s parameters. Each of the three single-axis search coils required its own measurement circuit. A block scheme of the measurement system is presented in
Figure 4.
The measurement circuit consisted of three parts: an inverting amplifier, a precision rectifier, and an averaging circuit with operational amplifiers. The inverting amplifier, together with an input buffer circuit, is presented in
Figure 5. The single-axis search coil was denoted with inductance
Lc,i, where
i can be
x,
y, or
z. One terminal of the coil was connected to the ground and the other terminal was connected to the input of the measurement circuit.
The buffer amplifier served as a high impedance load. The resistance of the coil could be neglected. Because the input resistance was so large, the input current was so small that it could also be neglected. Therefore, the measurement search coil did not induce a voltage in the transmitter coil circuit. The second part of the input amplifier circuit was the inverting amplifier. The gain of the amplifying part of the circuit is denoted with
Gi, and is defined with the values of the resistors
R1,i, and
R2,i. The gain
Gi can be expressed with:
The measurement circuit had variable gain, to enable the measurements of different voltage ranges and, therefore, fields of different strengths. In the case of this paper, the gain was set to
Gi = 1 with two 10 kΩ resistors in all three measurement circuits. The voltage induced on the search coils was sufficiently high and did not require additional amplification. In the case that coils with fewer turns would be used instead, the unity gain could be appropriately set a higher value. The input and the output voltages are presented in
Figure 6a,b. In
Figure 6a, the red line is the induced voltage, and the blue line is the amplified voltage of the first stage. The output of the first stage was connected to the input of the second stage. In
Figure 6b, the red signal is the input voltage induced in the search coil, and the blue line is the voltage at the output of the inverting amplifier.
The second part of the measurement circuit was a precision full wave-rectifier, which served to rectify the amplified input signal. It consisted of two amplifiers, fast signal diodes, and resistors. The full-wave rectifier is presented in
Figure 7. The operational amplifier compensated for the voltage drop across the diode. The precision rectifier consisted of two stages: a precision half-wave rectifier and an addition circuit.
The gain of the precision half-wave rectifier with a gain was defined by
R1,i and R
2,i. The second stage was the addition circuit, which added the original input signal to the output of the first stage. The gain of the addition circuit was determined with
R3,i, R5,i, and
R4,i. In the case of the measurement circuit, the gain of the measurement rectifier was 1. The values of the resistors required for unity gain rectification can be calculated using:
The input and the output voltages are presented in
Figure 8a,b. In
Figure 8a, the red line represents the voltage at the output of the inverting amplifier and the blue line represents the rectified voltage. The output of the second stage was connected to the input of the third stage. In
Figure 8b, the red signal is the output voltage of the inverting amplifier and the blue line is the rectified voltage. Both signals measured on the amplifier circuit exhibited oscillatory behavior. The cause of the damped oscillations was the self-resonance of the measurement coil. The transmitter coil was driven using square voltage, which also had spikes due to the self-resonance of the transmitter coil. Those voltage spikes caused the impulse response in the measurement coils, which resulted in the oscillations present in measured voltages.
The third and last part of the measurement circuit was the averaging circuit. The purpose of this circuit was to convert the rectified AC voltage to DC voltage. The output of the circuit was the average value of the input voltage. The averaging circuit could also be replaced with a peak detector circuit. The circuit of the third and last stage is presented in
Figure 9. It consisted of an input buffer amplifier, a low-pass RC filter, and an output buffer amplifier. The input voltage was rectified voltage
uabs,i, and the output voltage was the DC representation of the rectified AC voltage
UDC,i.
The input and output signals of the averaging circuit are presented in
Figure 10a,b. In
Figure 10a, the input rectified sinusoidal signal is marked with the red line and the output average value is marked with the blue line. In
Figure 10b, the rectified voltage signal is marked with the red color and the DC average value is marked with the blue color. The oscillatory behavior of the induced and rectified AC voltage was filtered using the RC low-pass filter. As a result, the self-resonance of the measurement search coil did not impact the induced voltage measurement.
The circuit in
Figure 9 averages the input waveform. The input of the circuit was the absolute sine waveform. Therefore, the average value
UDC,i of the voltage can be calculated using:
The
UDC,i was measured using the AD converter in the microcontroller. The peak value of the voltage
Uo,i can be calculated using:
The averaging window of the circuit (
Figure 9) was defined using the filter resistor and capacitor. The time constant of the circuit was defined with
R1,i and
C1,i:
The accurate average value of the input voltage could be measured after a duration of five time constants. The value of the voltage was measured using the AD converter on the microprocessor. Because the measurement circuit was used in discrete measurement intervals, the microprocessor could sample the voltage multiple times and then calculate the average value:
2.4. Bandwidth Evaluation of the First Stage of Measurement System
When measuring magnetic flux density, the bandwidth of the measuring circuit must also be considered, especially the bandwidth of the operational amplifier (OA) system that presented in
Figure 5. The input circuit included a buffer amplifier for impedance conditioning and an inverting amplifier. The frequency characteristics can be calculated as presented in [
15].
The open-loop gain
A(
jω) of the operational amplifier used in the measurement circuit was defined as:
where
Ao is the open-loop gain of the buffer amplifier at the frequency ω = 0 rad/s, and ω
t is the cross-over frequency when the gain of the open-loop is 0 dB. When the resistances
R1,i and
R2,i are included in the OA system, the closed-loop gain can be calculated for the non-inverting OA system
and for the inverting system (14):
After a short manipulation with (14) and (15) the gain for the non-inverting amplifier system was obtained:
The buffer amplifier with unity gain had the characteristics of the low-pass filter with the cutoff frequency at ω
t, and after manipulation with (14) and (16), it was obtained:
where
Gi =
R2,i/
R1,i and ω
b,i = ω
t/(1 +
R2,i/
R1,i). The gain
Gi represents the closed-loop gain of the inverting amplifier at the frequency of zero. The transfer function of the first stage of the input measurement scheme could then be defined with the equation:
A bode plot of the equations is presented in
Figure 11. The blue line presents the characteristics of the amplifier gain, the red line is the characteristics of the buffer amplifier, the yellow line is the characteristics of the inverting amplifier with a gain of 1, and the purple line is the characteristics of the transfer function of the first stage of the input measurement scheme. The cut-off frequency of the buffer amplifier of the measurement circuit was equal to cross-over frequency (for chosen OA
ft = 16.28 MHz) ω
t = 1.02 × 10
8, and the cut-off frequency of the second OA system (with resistors
R2,i and
R1,i), ω
b,i = ω
t/2 = 5.11 × 10
7 rad/s.
From (19) and
Figure 11, it follows that the pass-band was defined when the transfer function
Gi(
jω) drop for −3 dB, which was equal to ω
b,i and, according to the chosen measurement system frequency
fs = 87 kHz (ω
s = 5.47 × 10
5 rad/s), the selected OA system satisfied the measurement requirements. Other OAs in the measurement chain had the same cross-over frequency.
Due to component tolerances, the frequency response of an inverting amplifier can be different from the ideal one. Thus, a sensitivity analysis was performed to predict the margin of error in the closed-loop gain due to the tolerances of the resistors in the circuit. The sensitivity function used in the analysis was defined as:
where
x denotes the value of the component, or, in this case, closed-loop gain
Gi at frequency
ωb,i and
y denotes the circuit parameter in interest, which is the closed-loop gain
G2,i in absolute form. The absolute value of (19) is:
The gain
GA,i has sensitivities on two parameters,
Gi and ω
b,i. The sensitivities were evaluated as follows:
The low-frequency gain of the inverting amplifier was dependent on resistances
R2,i and
R1,i. If the resistors were chosen with ±1% tolerance, the gain
Gi had ±2% tolerance. From sensitivity Equation (22), the variation in closed-loop gain caused the following tolerance:
The sensitivity analysis for Equation (23) could be performed at two different frequencies of ω. The first frequency was at ω = 0 and the second was at ω = ω
b,i:
The variation of the closed-loop gain can be expressed as:
As stated before, because of the variation in the feedback loop resistance, the closed-loop gain was also variable. From the sensitivity analyses (24), (27) and (28) it was concluded that the maximum change in the closed-loop gain was in the range of ±2%.
2.5. Measurement Application
The computer application for measuring magnetic flux density was developed using the Microsoft Visual Studio using C#. The user interface of the application is presented in
Figure 12. The interface included the measurement circuit control panel and the transmitter coil position grid. The main purpose of the application was the interaction between the 3D positioning mechanism and the search coil magnetometer. The graphic positioning grid interface was used for positioning the bottom platform with the transmitter coil in the
x and
y-axes and setting the vertical
z distance between the 3D search coil and the transmitter coil.
To scan the magnetic field of the transmitter coil, the boundaries of the scan area must be defined using the start and the endpoints. The user can also determine the scanning step size. The smallest movement step in the
x,
y, or
z direction was set at 1 mm. With the scan area defined, the measurement of magnetic flux density can start. The search coil moves incrementally within the designated area. At each point in the space, the application sends the command to the search coil magnetometer to perform the voltage measurement. Data are then transferred from the search coil magnetometer to the application. The magnetic flux density is then calculated and saved in the Random Access Memory (RAM). If the measurement data is transferred correctly, the coil can be moved into the next position. At the end of the measurement, data from the RAM can be exported as a .csv file for further analysis and visualization. The measurement procedure is described using the flowchart in
Figure 13.
The app also included a 3D visualization function. The saved data, with extension .csv, can be imported and visualized in the 3D space. Each of the three components and the absolute value of the magnetic flux density can be represented with a magnetic flux density diagram.
3. Search Coil Magnetometer
The main part of the search coil magnetometer was a 3D search coil without a ferromagnetic core. The purpose of the search coil was magnetic field measurement in the 3D space, i.e., in the
x,
y, and
z directions. The 3D search coil consisted of three single axis search coils. The single axis search coils ere wound around the sphere and rotated 90° to each other. Each of the three coils was used to measure magnetic flux density along one specific axis. Because of the position and angle between the coils, the coupling coefficient between the coils was zero. This meant that the
x component of magnetic flux density induced voltage only in the
x-axis sense coil, i.e., it induced no voltage in the other two coils. The same was true for the
y and
z components of the magnetic flux density, which induced voltage only in the
y and
z single axis search coils, respectively. The modeled and fabricated 3D search coils are presented in
Figure 14.
A wireless power transmitter coil, measured using the search coil magnetometer, usually generates magnetic flux in all three directions. Each component of an electric field induces a voltage in one of the search coils. The induced voltage is proportional to the magnetic flux density.
In this paper, the 3D search coil was used in connection with the transmitter side of a typical IPT system. The primary inverter side of the IPT system was used to drive the transmitter coil with a sinusoidal current. The excited transmitter coil generated a magnetic field. The shape of the magnetic field was defined by the geometric properties of the coil. The strength of the field was determined by the primary transmitter coil current. To measure the generated magnetic field, the 3D search coil was moved in the space around the transmitter coil. The voltages induced on the 3D search coils were measured in each point. The connection between magnetic flux and the induced voltage was described using Faraday’s induction law,
where
ex, ey and
ez are the voltage induced across the
x,
y, and
z axis search coil, and
ϕx,
ϕy, and
ϕz are the magnetic flux in the
x,
y, and
z-axes directions. The connection between the magnetic flux and magnetic flux density can be described with the equation
where
Nx,
Ax,
Ny,
Ay, and
Nz,
Az are the number of turns and cross-sectional area of the
x,
y, and
z axis search coil, and
bi is the varying magnetic flux density in the
x,
y, and
z-axes. The cross-sectional area of the search coil magnetometer is presented in
Figure 15. The cross-sectional area of all three search coils was the same.
By manipulating Equation (30), the magnetic flux density can be calculated by integrating the search coil voltage:
The magnetic field was generated by the sinusoidal excitation current. Therefore, the induced voltage was also sinusoidal, and could be described using:
where
Êx,
Êy, and
Êz are the amplitude of the voltage induced in the
x,
y, and
z search coil, and ω is the angular frequency of the induced voltage. The integral of induced voltage could be calculated with the angular frequency of the inverter. Furthermore, the integral of the induced voltages (32) could be inserted in (31). The equation for calculating the magnetic flux density takes the form:
At each point in the
x,
y, and
z space, due to the change in magnetic flux density, the value of the induced voltage also changed. The amplitude of the magnetic flux density in the specific point in space can be calculated using:
The Equation (34) could then be used to calculate the magnetic flux density in three axes. To calculate the actual value of magnetic flux density in a specific point in the 3D space, the absolute value must be calculated from the measurements in all three axes:
The voltage induced in the search coil was measured with the previously described high impedance measurement circuit. The circuit transformed the AC voltage into DC voltage, which could be measured using the microcontroller. The measured DC voltage was used to calculate the amplitude of the magnetic flux density:
Because the amplitude of the induced voltage, measured during testing, was high enough, no additional amplification was required. The resistors R2,x, and R1,x were set to the fixed value of 10 kΩ. Therefore, the gain of the measurement circuit was 1.
3.1. Bandwidth and Self-Resonance of the Search Coil Magnetometer
The search coil of the magnetometer can be modeled using the AC voltage source, self-inductance
Lc,i, parasitic capacitance
Cc,i, which includes search coil capacitance and parasitic capacitance of the first amplifier stage, and DC resistance
Rc,i, where
i =
x,
y,
z. The circuit is presented in
Figure 16. Induced voltage can be measured across the parasitic capacitance
Cc,i using high impedance amplifier.
The transfer function of the search coil can be defined as:
The resonant frequency of the coil can be described with:
In the case of a 3D flux magnetometer, all coils have the same parameters. The self-inductance
Lc,i was 17.6 μH, the parasitic capacitance
Cc,i was 17.8 pF and the resistance of the wire was 1.4 Ω. Therefore, the self-resonant frequency of the single-axis search coil was 9 MHz (56.5 × 10
6 rad/s). The bode plot of transfer function (37) is presented in
Figure 17. The blue line represents the real transfer function with resonant frequency and the red line represents the ideal transfer function of the search coil. Under the resonant frequency, the ideal search coil behaved similarly to the nonideal search coil. The magnetic field measurements performed with the proposed system were under the self-resonant frequency of the coil. Therefore, the bandwidth of the magnetic field measurement circuit was limited by the precision rectifier and the RC filter.
3.2. Voltage Conversion Error of the Search Coil Magnetometer
As described in [
16], the error of the magnetic field to voltage conversion can be expressed with noise equivalent magnetic induction (
NEMI).
NEMI is defined as the output voltage noise of the search coil magnetometer in relation to the transfer function of the search coil magnetometer, expressed in
.
NEMI can be expressed with the equation:
where
vi,noise is the variance of the noise, due to the noise at the input of the first amplifier stage and the thermal noise:
where
vA,i is the variance of the noise at the input of the amplifier,
k is Boltzmann constant and
T is the temperature in Kelvin. Using the transfer function (37),
NEMI can be expressed with:
Because the search coil magnetometer measures varying magnetic fields with frequency under the resonant frequency of the coil, the Equation (40) can be simplified, by excluding the resonant frequency due to the search coil capacitance and self-inductance:
Form (42) it wasconcluded that NEMI increased with temperature and resistance of the search coil and decreased with the number of turns, and in a larger coil cross-section. It also decreased with the increment of the frequency of the measured magnetic field.
4. Measurement Results
To measure the magnetic flux density of the transmitter coil under the test in the 3D space, a measurement of the magnetic flux density was set up using the proposed system. The measurements were performed on two different coils: an ordinary, square planar coil on a square ferrite pad, and a directional DD coil on a square ferrite pad. Both coils are presented in
Figure 18.
Figure 18a presents the square planar coil with 10 turns, and
Figure 18b presents the DD coil, which consisted of two rectangular coils with 9 turns. Both coils had nearly the same inductance.
The parameters of the used system are presented in
Table 2, including the inductance of the square coil and the inductance of the DD coil, the transmitter coil excitation currents, input voltage of the inverter, and switching frequency of the phase-shifted PWM. The coils were designed to have similar inductance and footprints. Both coils were shielded with a 100 mm × 10 mm × 6 mm ferrite plate.
The measurements were performed in the
x-y planes at different z distances. The search coil was moved on the
x-y plane along with the user-defined path. The scan path of the search coil is presented in
Figure 19 with red lines. The arrows indicate the direction of the movement. The search coil followed the path in discrete incremental movements. At the end of each incremental movement, the magnetic field was measured in the 3D space.
The extract of the experimental measurements is presented in
Table 3 and
Table 4 for the square planar coil and DD coil, respectively. The measurement was performed on the
x-y plane at discrete distances. Three induced voltages were measured at each point on the plane.
Êx was the voltage induced in the
x-axis search coil,
Êy was the voltage induced in the
y-axis search coil, and
Êz was the voltage induced in the
z-axis search coil. The induced voltage was proportional to the magnetic flux density. The relation between induced voltage and magnetic flux density is described with Equation (36).
The comparisons between the simulated and measured results are presented in
Figure 20 and
Figure 21, respectively. The diagrams present the absolute magnetic flux density in three different
x-y planes, located at 15 mm, 20 mm, and 25 mm. The figures on the left side present the measured results and the figures on the right side present the simulated results. The results are presented using a color scale. The blue color defines the weakest magnetic flux density, and the red color defines the strongest magnetic field.
As described in the previous Section, the C# application was used to define the scan area and z-axis distance between the transfer coil and the 3D search coil. In both cases, the coils had the same, 10 cm × 10 cm footprint. The scan area was then defined as a 10 cm × 10 cm area around the centre of the coils. The measurements were performed at three z planes, at 15 mm, 20 mm, and 25 mm. The scan step on each x-y plane was set to the smallest value of 1 mm.
Figure 20 presents the measurement results of the spiral square coil, which radiated a nondirectional field. The current through the transmitter coil was set to 4 A.
Figure 20a,b presents the measured and simulated magnetic flux density in the
x-y plane at z = 15 mm. The biggest difference between the measured results and the simulation was at the center of the square coil. The measured magnetic field density was higher compared to the simulated value.
Figure 20c,d presents the difference between the measured and simulated magnetic flux density at z = 20 mm. Similar to the previous case, the measured magnetic flux density at the center of the coil was higher than the simulated results. Finally,
Figure 18 presents the measured and simulated magnetic flux density at z = 25 mm. Another notable difference in magnetic flux density was observed in the strength of the magnetic field in the center of the square coil.
The maximum value of magnetic flux density was above the wire at a distance of 15 mm. At the larger z distance, the magnetic flux density decreased, as can be seen in the measured and simulated figures. At all three distances, the measured magnetic flux density in the center of the coil was higher than the simulated values. When simulated, the maximum value of the magnetic flux density was above the conducting wire.
Figure 21 presents the measurement results of the DD coil, which radiated a directional field. The direction of the magnetic field was set along the
x-axis. The current through the transmitter coil was set to 4 A.
Figure 21a,b presents the measured and simulated magnetic flux density in the
x-y plane at
z = 15 mm. The biggest difference between the simulated and measured values of magnetic flux density was in the shape of the magnetic field. This difference was mainly due to the difference in the shape of the wound and 3D-modeled coils.
Figure 21c,d presents the difference between the measured and simulated magnetic flux density at
z = 20 mm. There was also some difference between the measured and simulated magnetic flux density. Finally,
Figure 21e,f presents the measured and simulated magnetic flux density at
z = 25 mm.
The maximum value of magnetic flux density was above the wire in the center of the coil at a distance of 15 mm. The maximum measured value was 1.55 mT. Compared to the magnetic flux density of the planar spiral coil, the magnetic flux density of DD coil was larger. In both cases, the coils were excited by a current with the same value. At greater distances, the magnetic flux density decreased, as can be seen in the measured and simulated figure data. Similar to the square planar coil, the measured magnetic flux density in the z-axis was more uniform than the simulated magnetic field, especially at greater z distances.
5. Discussion
The measured magnetic field could be used for wireless power transfer coil optimization. It could be used to visualize magnetic flux density in x-y planes at different z-axis distances. The measured magnetic field incuded x, y and, z-axis components. Each component could be displayed graphically. On the other hand, the simulation software was often only able to be used to present the absolute value of magnetic flux density. In the measurement results, only the absolute value of magnetic flux density was presented for the sake of comparison. The simulated magnetic flux density differed from the measured magnetic flux density. The largest differences were observed in the z-axis. The measured magnetic flux density in the z-axis was more uniform than the simulated value.
In IPT transmitter design, the most important direction of the field is the z-direction. This is because the receiver coil is usually oriented in the direction of the z-axis and in parallel to the transmitter coil; therefore, only the z-axis component of the magnetic field induces voltage in the receiver coil. If the receiver coil is not in parallel with transmitter coil and is rotated to it, the x and y components of the magnetic field also induce a voltage in the receiver coil and transfer power.
When designing wireless power transfer pads with better misalignment tolerance, the most important parameter is the magnetic flux density distribution in the z-axis. The transfer coils with the best magnetic flux density distributions are coils with a uniform directional magnetic field, such as DD coils.
The measurements were also performed with transmitter coils with a 12 mm thick ferrite plate. There was little to no difference observed between the measurement results with 6 mm thick ferrite. From a certain thickness on, ferrite did not significantly impact the magnetic field shape and strength.
A shortcoming of this system was that it was only able to evaluate performance and a field generated by a single transmitter coil. In this paper’s experimental configuration, it could not be used to measure the magnetic field between the transmitter and receiver coils during wireless power transfer. This was because the receiver coil could not be mounted on the current system. The measurement platform could, however, be expanded in the future, in order to allow the measurement of the magnetic flux during wireless power transfer.
6. Conclusions
This paper described the fabrication of a magnetic flux density measurement platform, designed for measurement and evaluation of different IPT transfer coil topologies. The measurement platform was designed based on a 3D search coil magnetometer and a 3D positioning mechanism. Due to Faraday’s Law of Induction, the varying magnetic field of the transmitter coil induced a voltage in the 3D search coil. The high impedance measurement circuit was used to measure and convert the induced AC voltage to DC voltage, which could then be measured using the AD converter and microcontroller. An application was also presented for interaction between a PC and the measurement platform. The measurement platform was used to evaluate the magnetic flux density of transmitter coils. The magnetic flux densities of two different transmitter coil structures were presented for comparison. The first was a planar square coil and the second was a DD coil. Additionally, the measurement results were compared to the results from EM simulation. The measurement platform could be used to evaluate and optimize wireless power transfer coils at the IPT design stage, together with EM simulation software. The platform would also be suitable for experimentation and development of new coil concepts. The platform was designed with flexibility in mind; for instance, different magnetometers could be used for magnetic field measurement. The search coil magnetometer could be replaced by a magnetometer, based on MEMS technology. In this paper, the search coil magnetometer enabled additional flexibility, which may not have been possible using different magnetic field measurement techniques.