Next Article in Journal
Enhancing Cross-Lingual Sarcasm Detection by a Prompt Learning Framework with Data Augmentation and Contrastive Learning
Next Article in Special Issue
Long-Range Wireless Power Transfer for Moving Wireless IoT Devices
Previous Article in Journal
Security and Trust in the 6G Era: Risks and Mitigations
Previous Article in Special Issue
ClosedLoop Control of Pulse-Density-Modulated Wireless Power Transfer with Fast MEPT
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Wireless Power Transfer System Model Reduction with Split Frequency Matching

1
Electric Power Research Institute, Yunnan Power Grid Co., Ltd., Kunming 650217, China
2
School of Electrical Engineering, Xi’an Jiaotong University, Xi’an 710049, China
*
Author to whom correspondence should be addressed.
Electronics 2024, 13(11), 2160; https://doi.org/10.3390/electronics13112160
Submission received: 18 April 2024 / Revised: 19 May 2024 / Accepted: 30 May 2024 / Published: 1 June 2024
(This article belongs to the Special Issue Wireless Power Transfer Technology and Its Applications)

Abstract

:
Reduced-order dynamic models of wireless power transfer (WPT) systems are desired to simplify the analysis and design of power control, phase synchronization, and maximum efficiency tracking. The reduced-order dynamic phasor model is a good choice because of its straightforward physical meaning and concise mathematical formula. However, the model relies on the assumption of loose coupling and loses accuracy when the coupling becomes stronger. In this paper, a model reduction method with split frequency matching is proposed to improve model accuracy under relatively strong coupling conditions, which is suitable for most short-distance WPT applications, such as wireless electrical vehicle charging. Split frequency matching is achieved through a pair of conjugate equivalent mutual inductances, which are derived from the asymmetry characteristics of the full-order dynamic phasor model in the positive and negative frequency domains. The proposed model retains the advantages of the existing model while significantly improving the accuracy under strong coupling conditions. Its characteristics are verified by comparing the experimental results and model predictions under both large step changes and small-signal perturbations.

1. Introduction

Wireless power transfer (WPT) systems based on magnetic resonance coupling play an important role in applications such as portable device charging and electrical vehicle (EV) charging [1,2,3]. The power regulation [4], online efficiency maximization [5], and dual-side phase synchronization [6] of a WPT system require closed-loop control and dynamic modeling. Conventionally, WPT systems are modeled using the generalized state-state averaging (GSSA) method [7,8,9], the extended describing function method (EDF) [10], or the dynamic phasor method [11,12]. However, as discussed by many previous studies [13,14,15], one of the inconveniences of these conventional modeling methods is the order doubling problem, because each ac variable in the system is represented by a complex-valued Fourier coefficient, whose magnitude and phase, or real and imaginary parts, are taken as two independent variables in the model. For example, the order of the model of a resonant tank with four energy-storage elements (two inductors and two capacitors for series–series topology) is eight. This problem leads to complicated model formulae, high computational cost, and difficulties in data-driven identification [16,17].
Classically, a high-order WPT system model can be simplified to a balanced reduced-order model with states eliminated based on Hankel singular values [18]. However, this numerical method loses the physical meaning of the original model and can only be performed on the real-valued linear state-space equations.
With the consideration of the periodical switching operation in WPT systems, a sampled data model is presented in [15]. The model takes the sampled currents and voltages at the end of switching periods as discrete state variables and calculates the piecewise state transition matrices with switching time constraints to obtain a difference equation that fully represents the system dynamics. The z-domain transfer functions are then derived from small-signal perturbation analysis. The sampled-data model avoids the order doubling problem and exhibits high accuracy in a wide frequency band. The main limitation of the model is the computational burden of determining the switching time, especially in closed-loop situations.
From the circuit point of view, an averaged model derived from coupled modes [13] avoids order doubling by describing each resonator with only two real-valued stable variables. The model is time-continuous and autonomous. It can be easily combined with classical control models to analyze closed-loop characteristics. However, the coupled-mode description is accurate only for loosely coupled resonators. Using the Laplace phasor transform [12], the dynamic phasor model of a pair of coupled resonators can be simplified to a pair of coupled equivalent inductances [14], and therefore, a complex-valued second-order equation that takes resonant current phasors as state variables is sufficient to describe a fourth-order resonant tank. This model is very simple and elegant in mathematics, but like the model based on coupled modes, the accuracy decreases as the coupling becomes stronger.
Currently, most high-power WPT applications, such as wireless EV charging [19], focus on efficient short-distance power transfer with relatively high coupling coefficients (for example, close to or even greater than 0.1). Under these conditions, the resonant tank exhibits a frequency splitting phenomenon and the existing reduced-order model [14] have large errors. In order to improve model accuracy under strong coupling conditions, this paper proposes a new dynamic phasor-based model reduction technique with split frequency [20,21] matching. The matching is achieved through a pair of conjugate equivalent mutual inductances, which can completely retain the low-frequency characteristics determined by the split beat frequencies. Compared with the existing reduced-order dynamic phasor model, the proposed model has similar simple second-order complex-valued formulae and a low computational cost, while the accuracy under strong coupling conditions is significantly improved.
The rest of the paper is organized as follows: Section 2 reviews the frequency splitting characteristics of a pair of coupled resonators. Section 3 provides the frequency domain analysis of the original dynamic phasor model. The proposed reduced-order model is derived in Section 4 and verified in Section 5. The conclusion is drawn in Section 6.

2. Split Frequencies of Coupled Resonators

The frequency splitting phenomenon of coupled resonators is manifested as the gain peaks of the network appear at the so-called split frequencies, which are different from the resonant frequencies. The reason is that the eigenvalues of the coupled resonators drift with the increase in the coupling coefficient [20,21]. To show this more clearly, let us consider the coupled resonators shown in Figure 1, where L1, C1 and L2, C2 are the primary- and secondary-side inductance and capacitance, respectively; M is the mutual inductance; u1 and u2 are the ac excitation voltages; uC1 and uC2 are the resonant voltages; and iL1 and iL2 are the resonant currents.
The instantaneous descriptor state-space model of the coupled resonators is given by
E d d t i L 1 t i L 2 t u C 1 t u C 2 t = A i L 1 t i L 2 t u C 1 t u C 2 t + B u 1 t u 2 t
where
E = L 1 M 0 0 M L 2 0 0 0 0 C 1 0 0 0 0 C 2 A = 0 0 1 0 0 0 0 1 1 0 0 0 0 1 0 0 B = 1 0 0 1 0 0 0 0
In (1), iL1 and iL2 are taken as the output because they are directly coupled with the sources. The transfer functions from u1 and u2 to iL1 and iL2 deduced from (1) are expressed as
i L 1 s u 1 s i L 1 s u 2 s i L 2 s u 1 s i L 2 s u 2 s = L 2 s s 2 + ω 2 2 M s 3 M s 3 L 1 s s 2 + ω 1 2 L 1 L 2 1 k 2 s 2 + ω o 2 s 2 + ω e 2
where k is the coupling coefficient:
k = M L 1 L 2
ω1 and ω2 are the resonant frequencies:
ω 1 = 1 L 1 C 1 ,     ω 2 = 1 L 2 C 2
ωo and ωe are the odd- and even-mode split frequencies:
ω o , e = ω 1 2 + ω 2 2 ω 1 2 ω 2 2 2 + 4 k 2 ω 1 2 ω 2 2 2 1 k 2
According to (3), the zeros of the coupled resonators locate at ±jω1, ±jω2, and 0, and the poles locate at ±jωo and ±jωe. Consequently, the gain valleys appear at ω1 and ω2, while the gain peaks appear at ωo and ωe, as shown by the bode diagrams in Figure 2.
The split frequencies ωo and ωe are functions of the resonant frequencies ω1 and ω2 and the coupling coefficient k. When k increases, ωo and ωe gradually move away from ω1 and ω2, as shown in Figure 3.

3. Frequency Domain Analysis of the Dynamic Phasor Model

The state-space model in Section II takes the instantaneous voltages and currents at state variables. The model is intuitive in nature, but it does not reveal the magnitude and phase relationships between the ac excitation voltages and the resonant currents, which are of interest in dynamic control. This gap can be filled by the generalized state-space model or the dynamic phasor model when considering only the fundamental harmonic. In this section, we review the dynamic phasor model and provide the analysis of the frequency domain as the foundation of model reduction.
The dynamic phasor, X, of a signal, x, is a complex-valued number given by
X t = 2 x 1 t = 2 T s t T s t x τ e j ω s τ d τ
where ωs is the switching frequency and Ts the corresponding period. Usually, ωs is set to the source frequency so that the dynamic phasor of the source is constant. The magnitude of X(t) corresponds to the amplitude of x, and the angle of X(t) is the phase of x. The derivative of X(t) with respect to time can be derived from (7) and expressed as
d X t d t = 2 d x d t 1 t 2 j ω s x 1 t
Let IL1(t), IL2(t), UC1(t), UC2(t), U1(t), and U2(t) represent the dynamic phasors of iL1(t), iL2(t), uC1(t), uC2(t), u1(t), and u2(t), respectively. By substituting the dynamic phasors and their derivatives into (1), it yields
E d d t I L 1 t I L 2 t U C 1 t U C 2 t = A j ω E I L 1 t I L 2 t U C 1 t U C 2 t + B U 1 t U 2 t
Consequently, the transfer functions from U1 and U2 to IL1 and IL2 are
I L 1 s U 1 s I L 1 s U 2 s I L 2 s U 1 s I L 2 s U 2 s = L 2 s s 2 + ω 2 2 M s 3 M s 3 L 1 s s 2 + ω 1 2 L 1 L 2 1 k 2 s 2 + ω o 2 s 2 + ω e 2
where
s = s + j ω s
which means that the transfer functions between the dynamic phasors are obtained by shifting the transfer functions between the instantaneous signals in the frequency domain by −ωs [22]. To illustrate this, Figure 4 compares the bode diagrams of (10) and (3). The bode diagrams are plotted along both positive and negative frequencies because (10) is complex-valued and is asymmetric in the frequency domain.
Furthermore, Figure 5 compares the zeros and poles of (3) and (10) on the complex plane. The zeros of (10) locate at ±jω1 − jωs, −jωs ± jω2, and −jωs, and the poles locate at −jωs ± jωo and −jωs ± jωe, which are the translations of the zeros and poles of (3) along the imaginary axis, as indicated by (11). Considering that ω1, ω2, ωo, ωe, and ωs are very close to each other, the dominant zeros of (10) are −jωs + jω1 and −jωs + jω2, and the dominant poles are −jωs + jωo and −jωs + jωe. On the complex plane, these dominant zeros and poles are within the circle whose radius is equal to the Nyquist frequency, i.e., ωs/2, as shown in Figure 5. In a physical sense, the dominant zeros correspond to the beat frequencies, i.e.,
Δ ω 1 = ω s ω 1 ,     Δ ω 2 = ω s ω 2
and the dominant poles correspond to the split beat frequencies, i.e.,
Δ ω o = ω s ω o ,     Δ ω e = ω s ω e

4. Order Reduction with Split Frequency Matching

The dynamic phasor model focuses on the signal magnitude and phase, which are low-frequency characteristics as compared with the switching frequency, ωs. From this point of view, it is desired to find a reduced-order model that only contains dominant zeros and poles, i.e., −jωs+jω1, −jωs+jω2, −jωs+jωo, and −jωs+jωe. As there are two dominant poles and two dominant zeros, a reasonable formula for the reduced-order model is
E r d d t I L 1 t I L 2 t = A r I L 1 t I L 2 t + B r U 1 t U 2 t
where Er, Ar, and Br are 2-by-2 coefficient matrices. For simplicity, one can let Br be
B r = 1 0 0 1
and then determine Er and Ar by matching the steady-state gains and dominant zeros and poles of (14) with those of (9).
According to (9), the steady-state values of IL1, IL2, U1, and U2 satisfy
I ¯ L 1 I ¯ L 2 = A j ω s E 1 B U ¯ 1 U ¯ 2
On the other hand, according to (14):
I ¯ L 1 I ¯ L 2 = A r 1 B r   U ¯ 1 U ¯ 2
To match the results of (16) and (17), let
A r 1 B r = A j ω s E 1 B
Then, it yields
A r = j Δ ω 1 L ω 1 j ω s M j ω s M j Δ ω 2 L ω 2
where Lω1 and Lω2 are the equivalent inductances [14] given by
L ω 1 = ω s + ω 1 ω s L 1 ,     L ω 2 = ω s + ω 2 ω s L 2
Let e11, e12, e21, and e22 be the undermined elements of Er, i.e.,
E r = e 11 e 12 e 21 e 22
With (21), the reduced transfer functions are deduced from (14) and expressed as
I L 1 s U 1 s I L 1 s U 2 s I L 2 s U 1 s I L 2 s U 2 s = e 22 s + j Δ ω 2 L ω 2 e 12 s j ω s M e 21 s j ω s M e 11 s + j Δ ω 1 L ω 1 Δ s
where Δ(s) is the characteristic polynomial:
Δ s = e 11 s + j Δ ω 1 L ω 1 e 22 s + j Δ ω 2 L ω 2 e 12 s + j ω s M e 21 s + j ω s M
To set the dominant zeros of (22) at −jωs + jω1 and −jωs + jω2, it is solved that
e 11 = L ω 1 ,     e 22 = L ω 2
With (24), to set the dominant poles of (22) at −jωs + jωo and −jωs + jωe, it is solved that
e 12 = M ω ,     e 21 = M ω *
where Mω and Mω* are a pair of conjugate equivalent mutual inductances, which can be calculated using (26) and (27), derived from to the eigen equations of (9) and (22).
M ω + M ω * = 1 k 2 · ω s + ω 1 ω s · ω s + ω 2 ω s ·     Δ ω 1 + Δ ω 2 ω s + k 2 1 k 2 · ω s + ω o ω s · ω s + ω e ω s · Δ ω o + Δ ω e ω s M
M ω M ω * = 1 k 2 · ω s + ω 1 ω s · ω s + ω 2 ω s + k 2 1 k 2 · ω s + ω o ω s · ω s + ω e ω s M 2
Figure 6 shows the real and imaginary parts of Mω as functions of coupling coefficient k with ω1 = ωs ±5% and ω2 = ωs ±5%. It is found that Re[Mω] and Im[Mω] are almost linear with k in a wide range, like the mutual inductance M. The proportional factors of Re[Mω] and Im[Mω] to k depend on the relative error of ω1 and ω2 to ωs but do not have large dispersions. Therefore, one can calculate Mω and Mω* using (28) and (29) instead of (26) and (27) for simplicity
lim k 0 M ω + M ω * = 2 + ω 1 + ω 2 2 ω s + ω 1 ω 2 2 ω 1 + ω 2 ω s ω 1 2 ω 2 2 2 ω 1 + ω 2 ω s 3 M 3 M
lim k 0 M ω M ω * = 1 + ω 1 + ω 2 2 ω s + ω 1 ω 2 2 ω 1 + ω 2 ω s + ω 1 ω 2 2 ω s 2 M 2 2.75 M 2
and obtain
M ω 3 2 + 1 2 j M ,     M ω * 3 2 1 2 j M
With (24) and (25), Er is rewritten as:
E r = L ω 1 M ω M ω * L ω 2
and the reduced transfer function is rewritten as
I L 1 s U 1 s I L 1 s U 2 s I L 2 s U 1 s I L 2 s U 2 s = L ω 2 s + j Δ ω 2 L ω 2 M ω s j ω s M M ω * s j ω s M L ω 1 s + j Δ ω 1 L ω 1 L ω 1 L ω 2 M ω M ω * s + j Δ ω o s + j Δ ω e
where
Δ s = L ω 1 L ω 2 M ω M ω * s + j Δ ω o s + j Δ ω e
Figure 7 compares the bode diagrams of transfer function (32) derived from the reduced-order model to the bode diagrams of transfer function (10) derived from the original dynamic phasor model (10) within the frequency range from −0.1ωs to 0.1ωs. It is shown that the poles exactly locate at the split beat frequencies, and the two groups of transfer functions match with each other well, even when ω2, ω2, and ωs are not identical. This means that the low-frequency characteristics of the original dynamic phasor model (9) can be precisely represented by the reduced-order model (14).

5. Reduced Wpt System Model

The coupled resonant network is essentially a linear system, and therefore, it is analyzed in the frequency domain in the above sections for order reduction. However, in a dc-to-dc WPT system, as shown in Figure 8, the most concerned relationships, e.g., the transfer functions between magnitudes used for power control, and the transfer functions between phases used for synchronization, can only be derived from a locally linearized small-signal model because the magnitudes and phases are nonlinear functions of dynamic phasors.
In Figure 8, the switching functions of the two converters are denoted by s1 and s2. For full-bridge converters, s1 and s2 may equal to 1, 0, and −1, depending on the switching control (for an active converter) or resonant current direction (for a diode rectifier). Consequently, the ac excitation voltages u1 and u2 are represented by the products of the switch functions and the dc side voltages v1 and v2 as
u 1 t = s 1 t v 1 t u 2 t = s 2 t v 2 t
and the dc side currents i1 and i2 are represented by the products of the switch functions and the resonant currents iL1 and iL2 as
i 1 t = s 1 t i L 1 t i 2 t = s 2 t i L 2 t
With the fundamental harmonic approximation, (33) and (34) can be rewritten as
U 1 t = S 1 t V 1 t U 2 t = S 2 t V 2 t  
and
I 1 t = R e S 1 * t I L 1 t I 2 t = R e S 2 * t I L 2 t
where S1 and S2 are the dynamic phasors of s1 and s2, respectively, the superscript asterisk means to take the conjugate, and V1, V2, I1, and I2 are the moving averages of v1, v2, i1, and i2, respectively.
According to (35), the Jacobian matrix of R e U 1 , R e U 2 , I m U 1 , I m U 2 T with respect to V 1 , S 1 , S 1 , V 2 , S 2 , S 2 T is
J B = 1 0 0 0 0 0 1 0 0 1 0 0 0 0 0 1 U 1 r V 1 U 1 r S 1 U 1 i U 1 i V 1 U 1 i S 1 U 1 r 0 0 U 2 r V 2 U 2 r S 2 U 2 i U 2 i V 2 U 2 i S 2 U 2 r
According to (36), the Jacobian matrix of I 1 , I L 1 , I L 1 , I 2 , I L 2 , I L 2 T with respect to R e I L 1 , R e I L 2 , I m I L 1 , I m I L 2 T is
J C = S 1 r S 1 i I L 1 r I L 1 I L 1 i I L 1 I L 1 i I L 1 2 I L 1 r I L 1 2 0 0 S 2 r S 2 i I L 2 r I L 2 I L 2 i I L 2 I L 2 i I L 2 2 I L 2 r I L 2 2 1 0 0 0 0 0 1 0 0 1 0 0 0 0 0 1
and the Jacobian matrix of I 1 , I L 1 , I L 1 , I 2 , I L 2 , I L 2 T with respect to V 1 , S 1 , S 1 , V 2 , S 2 , S 2 T is
J D = 0 I 1 S 1 I m S 1 * I L 1 0 0 0 0 0 0 0 0 0 I 2 S 2 I m S 2 * I L 2 0 0 0 0 0 0
Meanwhile, the reduced-order model (14) is rewritten in an expanded real-valued form as
f E r d d t R e I L 1 t R e I L 1 t I m I L 2 t I m I L 2 t = f A r R r R e I L 1 t R e I L 1 t I m I L 2 t I m I L 2 t + f B r R e U 1 t R e U 2 t I m U 1 t I m U 2 t
where Rr is the equivalent series resistance (ESR) matrix introduced by R1 and R2:
R r = R 1 0 0 R 2
and f maps a complex matrix Z to an expanded real-valued matrix as
f Z = R e Z I m Z I m Z R e Z
By substituting (37)–(39) into (40), the small-signal model of the WPT system is obtained:
f E r d d t R e I ^ L 1 t R e I ^ L 1 t I m I ^ L 2 t I m I ^ L 2 t = f A r R r R e I ^ L 1 t R e I ^ L 1 t I m I ^ L 2 t I m I ^ L 2 t + f B r J B V ^ 1 t S ^ 1 t S ^ 1 t V ^ 2 t S ^ 2 t S ^ 2 t I ^ 1 t I ^ L 1 t I ^ L 1 t I ^ 2 t I ^ L 2 t I ^ L 2 t = J C R e I ^ L 1 t R e I ^ L 1 t I m I ^ L 2 t I m I ^ L 2 t + J D V ^ 1 t S ^ 1 t S ^ 1 t V ^ 2 t S ^ 2 t S ^ 2 t
Equation (43) is a 4th-order real-valued small-signal model with 6 input signals and 6 output signals, including both magnitudes and phases of the control signals, dc voltages, and resonant currents.
Based on (43), the 6-by-6 small-signal transfer function matrix from the input vector V 1 ,   S 1 , S 1 , V 2 , S 2 , S 2 T to the output vector I 1 ,   I L 1 , I L 1 , I 2 , I L 2 , I L 2 T , which may contain all the concerned information of the WPT system, is derived from (43) and expressed as
I ^ 1 s V ^ 1 s I ^ 1 s S ^ 1 s I ^ 1 s S ^ 1 s I ^ L 1 s V ^ 1 s I ^ L 1 s S ^ 1 s I ^ L 1 s S ^ 1 s I ^ L 1 s V ^ 1 s I ^ L 1 s S ^ 1 s I ^ L 1 s S ^ 1 s I ^ 1 s V ^ 2 s I ^ 1 s S ^ 2 s I ^ 1 s S ^ 2 s I ^ L 1 s V ^ 2 s I ^ L 1 s S ^ 2 s I ^ L 1 s S ^ 2 s I ^ L 1 s V ^ 2 s I ^ L 1 s S ^ 2 s I ^ L 1 s S ^ 2 s I ^ 2 s V ^ 1 s I ^ 2 s S ^ 1 s I ^ 2 s S ^ 1 s I ^ L 2 s V ^ 1 s I ^ L 2 s S ^ 1 s I ^ L 2 s S ^ 1 s I ^ L 2 s V ^ 1 s I ^ L 2 s S ^ 1 s I ^ L 2 s S ^ 1 s I ^ 2 s V ^ 2 s I ^ 2 s S ^ 2 s I ^ 1 s S ^ 2 s I ^ L 2 s V ^ 2 s I ^ L 2 s S ^ 2 s I ^ L 2 s S ^ 2 s I ^ L 2 s V ^ 2 s I ^ L 2 s S ^ 2 s I ^ L 2 s S ^ 2 s = J C s f E r f A r R r 1 f B r J B + J D

6. Experimental Verification

The accuracy of the proposed reduced-order model is verified in this section by the low-power experimental system shown in Figure 9. The system consisted of a pair of coupled resonators made of AWG 46 Litz wires, high-Q ceramic capacitors, and two gallium-nitride full-bridge converters controlled by Zynq development boards. To focus on the dynamics of the WPT stage, the converters were directly fed by bidirectional dc supplies so that the characteristics predicted by the reduced-order state-space model (14) and the transfer functions in (44) could be directly measured and verified. The nominal operating parameters of the system are listed in Table 1, and the steady-state operating waveforms are shown in Figure 10. The transferred power was 100 W and the dc-to-dc efficiency was 93.5%. Thanks to the high-efficiency GaN devices and zero-voltage switching (ZVS) operation [23], the losses on the converters were insignificant and most of the system power loss came from the ESR of the coils. For model verification, the representative dynamic characteristics of the system, including both large-signal step responses under ON–OFF control, and the small-signal bode diagrams at steady states, were measured and compared to the model predictions under different test conditions, as listed in Table 2.

6.1. Large-Signal Characteristics

A simple control method for low-power WPT applications is the ON–OFF keying method introduced in [24]. The average power is controlled by periodically operating the primary-side converter between ON and OFF periods, e.g., for a full-bridge converter, switching s1 between 1 and −1 during ON periods and keeping s1 = 0 during OFF periods.
Figure 11 shows the waveforms of secondary-side resonant current iL2 at the beginning of an ON period. The magnitude of iL2, i.e., |IL2|, oscillated quickly when the coupling was strong and oscillated slowly when the coupling was weak. The reason is that the oscillations depend on the split beat frequencies Δωo and Δωe, and the split beat frequencies depend on the coupling coefficient. When the system was slightly detuned, the oscillation mode became more complicated because the beat frequencies Δω1 and Δω2 were no longer zero. Under all tested operation conditions, the proposed reduced-order model (14) and the original dynamic phasor model matched the experimental results (envelopes) well. In contrast, the existing reduced-order model [14] was accurate, only with weak coupling, and had accumulated errors when the coupling was relatively strong.

6.2. Small-Signal Characteristics

Synchronizing the secondary-side converter of a WPT system with its primary-side converter is crucial for stable power transfer. For an active secondary-side converter, the synchronization requires a phase-locked loop that adjusts the phase of the converter drive signal according to the phase of the measured resonant current [25]. The stability of the synchronization depends on the closed-loop phase dynamics, which involve the small-signal transfer function from the drive signal phase angle, i.e., S 2 , to the resonant current phase angle, i.e., I L 2 .
In the experiment, to obtain the small-signal bode diagram of the transfer function I L 2 s / S 2 s , a low-frequency (10~100 kHz) phase perturbation signal of 7.2° magnitude (corresponding to 20 ns) was injected into the steady-state S 2 using a variable time-delay unit, and the responded I L 2 was measured by a high-bandwidth current transformer (CST4835), a zero-crossing detector (LMH7220), and a digital phase detector on the Zynq board. Then, the measured time-domain data of S 2 and I L 2 were processed with discrete Fourier transform to obtain I L 2 s / S 2 s at a particular perturbation frequency. Figure 12 shows the small-signal bode diagrams at the steady states of different operating conditions. When the coupling was strong, the gain curve had two peaks that corresponded to the split beat frequencies Δωo and −Δωe, as shown in Figure 12a,b. When the coupling was weak, the two gain peaks merged into one, as shown in Figure 12c. Under all tested operation conditions, both the proposed reduced-order model and the original dynamic phasor model accurately predicted the small-signal characteristics. However, the existing reduced-order model failed to predict the two gain peaks under the tuned and strong coupling condition, and had a bit peak frequency error under the detuned condition.

7. Conclusions

This paper improved the accuracy of the reduced-order dynamic phasor model of WPT systems with strongly coupled resonators, where the split frequencies are dominant in the dynamic characteristics. Compared with the existing reduced-order model, the major change is a pair of conjugate equivalent mutual inductance derived from split beat frequency matching. Due to the equivalent mutual inductances, the asymmetry characteristics of the full-order dynamic phasor model in positive and negative frequency domains are completely preserved. The proposed model maintains a simple formula and low computational cost, while significantly improving accuracy under strong coupling conditions, as revealed by the large-signal and small-signal experimental tests. The model is useful for the dynamic analysis and control design of WPT systems with either weak or strong coupling, such as short-distance wireless EV charging.

Author Contributions

Conceptualization, K.W. and H.L.; methodology, H.L. and Q.W.; hardware design and software, K.W. and J.P.; validation, K.W. and Q.W.; writing—original draft preparation, Q.W.; writing—review and editing, H.L.; supervision, H.L.; project administration, K.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by China Southern Power Grid under grant YNKJXM20222524.

Data Availability Statement

Some or all of the data during this study are available from the corresponding author upon request.

Conflicts of Interest

Authors K.W. and J.P. were also employed by Yunnan Grid Electric Power Research Institute, Yunnan Power Grid Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial of financial relationships that could be construed as a potential conflict of interest.

References

  1. Buja, G.; Choi, S.Y.; Rim, C.T. Modern Advances in Wireless Power Transfer Systems for Roadway Powered Electric Vehicles. IEEE Trans. Ind. Electron. 2016, 63, 6533–6545. [Google Scholar] [CrossRef]
  2. Hui, S.Y. Planar Wireless Charging Technology for Portable Electronic Products and Qi. Proc. IEEE 2013, 101, 1290–1301. [Google Scholar] [CrossRef]
  3. Zhang, Z.; Pang, H.; Georgiadis, A.; Cecati, C. Wireless power transfer—An overview. IEEE Trans. Ind. Electron. 2019, 66, 1044–1058. [Google Scholar] [CrossRef]
  4. Zhong, W.; Hui, S.R. Charging Time Control of Wireless Power Transfer Systems Without Using Mutual Coupling Information and Wireless Communication System. IEEE Trans. Ind. Electron. 2017, 64, 228–235. [Google Scholar] [CrossRef]
  5. Huang, Z.; Wong, S.-C.; Tse, C.K. Control design for optimizing efficiency in inductive power transfer systems. IEEE Trans. Power Electron. 2018, 33, 4523–4534. [Google Scholar] [CrossRef]
  6. Jiang, Y.; Wang, L.; Wang, Y.; Wu, M.; Zeng, Z.; Liu, Y.; Sun, J. Phase-locked loop combined with chained trigger mode used for impedance matching in wireless high power transfer. IEEE Trans. Power Electron. 2020, 35, 4272–4285. [Google Scholar] [CrossRef]
  7. Hu, A.P. Modeling a contactless power supply using GSSA method. In Proceedings of the 2009 IEEE International Conference on Industrial Technology, Churchill, VIC, Australia, 10–13 February 2009; pp. 1–6. [Google Scholar]
  8. Hao, H.; Covic, G.A.; Boys, J.T. An approximate dynamic model of LCL-T-based inductive power transfer power supplies. IEEE Trans. Power Electron. 2014, 29, 5554–5567. [Google Scholar] [CrossRef]
  9. Sanders, S.R.; Noworolski, J.M.; Liu, X.Z.; Verghese, G.C. Generalized averaging method for power conversion circuits. IEEE Trans. Power Electron. 1991, 6, 251–259. [Google Scholar] [CrossRef]
  10. Yang, E.X.; Lee, F.C.; Jovanovic, M.M. Small-signal modeling of series and parallel resonant converters. In Proceedings of the IEEE Applied Power Electronics Conference and Exposition (APEC), Boston, MA, USA, 23–27 February 1992; pp. 785–792. [Google Scholar]
  11. Rim, C.T.; Cho, G.H. Phasor transformation and its application to the DC/AC analyses of frequency phase-controlled series resonant converters (SRC). IEEE Trans. Power Electron. 1990, 5, 201–211. [Google Scholar] [CrossRef]
  12. Lee, S.; Choi, B.; Rim, C.T. Dynamics characterization of the inductive power transfer system for online electric vehicles by Laplace phasor transform. IEEE Trans. Power Electron. 2013, 28, 5902–5909. [Google Scholar] [CrossRef]
  13. Li, H.; Wang, K.; Huang, L.; Chen, W.; Yang, X. Dynamic Modeling Based on Coupled Modes for Wireless Power Transfer Systems. IEEE Trans. Power Electron. 2015, 30, 6245–6253. [Google Scholar] [CrossRef]
  14. Li, H.; Fang, J.; Tang, Y. Dynamic phasor-based reduced-order models of wireless power transfer systems. IEEE Trans. Power Electron. 2019, 34, 11361–11370. [Google Scholar] [CrossRef]
  15. Tang, J.; Dong, S.; Cui, C.; Zhang, Q. Sampled-Data Modeling for Wireless Power Transfer Systems. IEEE Trans. Power Electron. 2020, 35, 3173–3182. [Google Scholar] [CrossRef]
  16. Chen, F.; Padilla, A.; Young, P.C.; Garnier, H. Data-Driven Modeling of Wireless Power Transfer Systems With Slowly Time-Varying Parameters. IEEE Trans. Power Electron. 2020, 35, 12442–12456. [Google Scholar] [CrossRef]
  17. Chen, F.; Young, P.C.; Garnier, H.; Deng, Q.; Kazimierczuk, M.K. Data-Driven Modeling of Wireless Power Transfer Systems With Multiple Transmitters. IEEE Trans. Power Electron. 2020, 35, 11363–11379. [Google Scholar] [CrossRef]
  18. Gawronski, W.; Juang, J.-N. Model reduction in limited time and frequency intervals. Int. J. Syst. Sci. 1990, 21, 349–376. [Google Scholar] [CrossRef]
  19. SAE J2954_202010; Wireless Power Transfer for Light-Duty Plug-In and Electric Vehicles and Alignment Methodology. SAE International: Warrendale, PA, USA, 2016.
  20. Niu, W.Q.; Chu, J.X.; Gu, W.; Shen, A.D. Exact Analysis of Frequency Splitting Phenomena of Contactless Power Transfer Systems. IEEE Trans. Circuits Syst. I-Regul. Pap. 2013, 60, 1670–1677. (In English) [Google Scholar] [CrossRef]
  21. Haus, H.A.; Huang, W. Coupled-mode theory. Proc. IEEE 1991, 79, 1505–1518. [Google Scholar] [CrossRef]
  22. Tan, T.; Chen, K.; Lin, Q.; Jiang, Y.; Yuan, L.; Zhao, Z. Impedance Shaping Control Strategy for Wireless Power Transfer System Based on Dynamic Small-Signal Analysis. IEEE Trans. Circuits Syst. I 2021, 68, 1354–1365. [Google Scholar] [CrossRef]
  23. Li, H.; Wang, K.; Fang, J.; Tang, Y. Pulse density modulated ZVS full-bridge converters for wireless power transfer systems. IEEE Trans. Power Electron. 2019, 34, 369–377. [Google Scholar] [CrossRef]
  24. Zhong, W.; Hui, S.Y.R. Maximum energy efficiency operation of series-series resonant wireless power transfer systems using on-off keying modulation. IEEE Trans. Power Electron. 2018, 33, 3595–3603. [Google Scholar] [CrossRef]
  25. Jia, S.; Chen, C.; Liu, P.; Duan, S. A digital phase synchronization method for bidirectional inductive power transfer. IEEE Trans. Ind. Electron. 2020, 67, 6450–6460. [Google Scholar] [CrossRef]
Figure 1. Magnetically coupled series resonators.
Figure 1. Magnetically coupled series resonators.
Electronics 13 02160 g001
Figure 2. Bode diagrams of the coupled resonators.
Figure 2. Bode diagrams of the coupled resonators.
Electronics 13 02160 g002
Figure 3. Split frequencies as a function of coupling coefficient under three different cases: ω1 = ω2, ω1 < ω2, and ω2 < ω1.
Figure 3. Split frequencies as a function of coupling coefficient under three different cases: ω1 = ω2, ω1 < ω2, and ω2 < ω1.
Electronics 13 02160 g003
Figure 4. Bode diagrams of instantaneous model and dynamic phasor model when k = 0.1.
Figure 4. Bode diagrams of instantaneous model and dynamic phasor model when k = 0.1.
Electronics 13 02160 g004
Figure 5. Zeros and poles of instantaneous model (left) and dynamic phasor model (right).
Figure 5. Zeros and poles of instantaneous model (left) and dynamic phasor model (right).
Electronics 13 02160 g005
Figure 6. The real and imaginary parts of the equivalent mutual inductance, Mω, as a function of coupling coefficient k.
Figure 6. The real and imaginary parts of the equivalent mutual inductance, Mω, as a function of coupling coefficient k.
Electronics 13 02160 g006
Figure 7. Bode diagrams of original and reduced dynamic phasor models when k = 0.1, ω1 = 0.98ωs, and ω2 = 1.02ωs.
Figure 7. Bode diagrams of original and reduced dynamic phasor models when k = 0.1, ω1 = 0.98ωs, and ω2 = 1.02ωs.
Electronics 13 02160 g007
Figure 8. Dual-voltage-sourced WPT system based on magnetically coupled series resonators.
Figure 8. Dual-voltage-sourced WPT system based on magnetically coupled series resonators.
Electronics 13 02160 g008
Figure 9. Experimental WPT system.
Figure 9. Experimental WPT system.
Electronics 13 02160 g009
Figure 10. Steady-state operating waveforms with nominal operating parameters.
Figure 10. Steady-state operating waveforms with nominal operating parameters.
Electronics 13 02160 g010
Figure 11. Time-domain responses of secondary-side resonant current when turning on the primary-side converter when the system is (a) tuned with strong coupling, (b) detuned with strong coupling, and (c) with weak coupling.
Figure 11. Time-domain responses of secondary-side resonant current when turning on the primary-side converter when the system is (a) tuned with strong coupling, (b) detuned with strong coupling, and (c) with weak coupling.
Electronics 13 02160 g011
Figure 12. Bode diagram of I L 2 s / S 2 s when the system is (a) tuned with strong coupling, (b) detuned with strong coupling, and (c) with weak coupling.
Figure 12. Bode diagram of I L 2 s / S 2 s when the system is (a) tuned with strong coupling, (b) detuned with strong coupling, and (c) with weak coupling.
Electronics 13 02160 g012aElectronics 13 02160 g012b
Table 1. Nominal operating parameters.
Table 1. Nominal operating parameters.
SymbolQuantityValue
L1,2Resonant inductances63.3 μH
C1,2Resonant capacitances400 pF
fr1,2Resonant frequencies1 MHz
fsSwitching frequency1 MHz
R1,2ESRs~1 Ω
Q1,2Quality factors400
kCoupling coefficient0.1
DPower transfer distance15 cm
V1,2DC supply voltages40 V
Table 2. Test conditions.
Table 2. Test conditions.
Operating ConditionCoupling CoefficientSwitching Frequency
Tuned resonance w/strong couplingk = 0.1fs = 1 MHz
Detuned resonance w/strong couplingk = 0.1fs = 1.02 MHz
Tuned resonance w/week couplingk = 0.03fs = 1 MHz
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, K.; Wu, Q.; Peng, J.; Li, H. Wireless Power Transfer System Model Reduction with Split Frequency Matching. Electronics 2024, 13, 2160. https://doi.org/10.3390/electronics13112160

AMA Style

Wang K, Wu Q, Peng J, Li H. Wireless Power Transfer System Model Reduction with Split Frequency Matching. Electronics. 2024; 13(11):2160. https://doi.org/10.3390/electronics13112160

Chicago/Turabian Style

Wang, Ke, Qingyu Wu, Jing Peng, and Hongchang Li. 2024. "Wireless Power Transfer System Model Reduction with Split Frequency Matching" Electronics 13, no. 11: 2160. https://doi.org/10.3390/electronics13112160

APA Style

Wang, K., Wu, Q., Peng, J., & Li, H. (2024). Wireless Power Transfer System Model Reduction with Split Frequency Matching. Electronics, 13(11), 2160. https://doi.org/10.3390/electronics13112160

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop