1. 单相整流器dq坐标系下建模
单相整流器的拓扑如图所示,可知
u
a
b
=
u
s
−
L
d
i
s
d
t
−
R
i
s
{u_{ab}} = {u_{s}} - L\frac{{d{i_s}}}{{dt}} - R{i_s}
uab=us−Ldtdis−Ris。
将电压和电流写成dq的形式。
{
u
s
=
U
s
m
sin
(
ω
t
)
i
s
=
I
d
sin
(
ω
t
)
+
I
q
cos
(
ω
t
)
u
a
b
=
U
d
sin
(
ω
t
)
+
U
q
cos
(
ω
t
)
\left\{\begin{array}{l} u_{\mathrm{s}}=U_{\mathrm{sm}} \sin (\omega \mathrm{t}) \\ i_{\mathrm{s}}=I_{\mathrm{d}} \sin (\omega \mathrm{t})+I_{\mathrm{q}} \cos (\omega \mathrm{t}) \\ u_{\mathrm{ab}}=U_{\mathrm{d}} \sin (\omega \mathrm{t})+U_{\mathrm{q}} \cos (\omega \mathrm{t}) \end{array}\right.
⎩
⎨
⎧us=Usmsin(ωt)is=Idsin(ωt)+Iqcos(ωt)uab=Udsin(ωt)+Uqcos(ωt)
将dq表达式代入
u
a
b
=
u
s
−
L
d
i
s
d
t
−
R
i
s
{u_{ab}} = {u_{s}} - L\frac{{d{i_s}}}{{dt}} - R{i_s}
uab=us−Ldtdis−Ris。需要注意的是,Idsin(wt)的求导包含两部分,一部分是对Id求偏导,另一部分是对sinwt求偏导。dq表达式中的其他项也是如此。
最终得到表达式。
{
U
d
=
U
s
m
−
L
d
I
d
d
t
+
I
q
ω
L
−
R
I
d
U
q
=
0
−
L
d
I
q
d
t
−
I
d
ω
L
−
R
I
q
\left\{ \begin{array}{l} {U_d} = {U_{sm}} - L\frac{{d{I_d}}}{{dt}} + {I_q}\omega L - R{I_d}\\ {U_q} = 0 - L\frac{{d{I_q}}}{{dt}} - {I_d}\omega L - R{I_q} \end{array} \right.
{Ud=Usm−LdtdId+IqωL−RIdUq=0−LdtdIq−IdωL−RIq 上述的Ud和Uq是实际值,为了求得目标的Ud和Uq,以d轴为例,需要将PI控制器替代
L
d
I
d
d
t
+
R
I
d
L\frac{{d{I_d}}}{{dt}} +R{I_d}
LdtdId+RId,能够得到以下表达式:
{
L
d
I
d
d
t
≈
L
I
d
r
e
f
−
I
d
T
s
=
K
p
(
I
d
r
e
f
−
I
d
)
R
I
d
=
R
L
∫
L
d
I
d
d
t
d
t
≈
R
L
∫
K
p
(
I
d
r
e
f
−
I
d
)
d
t
=
K
i
∫
(
I
d
r
e
f
−
I
d
)
d
t
K
i
=
R
L
K
p
\left\{ \begin{array}{l} L\frac{{d{I_d}}}{{dt}} \approx L\frac{{{I_{dref}} - {I_d}}}{{{T_s}}} = {K_p}({I_{dref}} - {I_d})\\ R{I_d} = \frac{R}{L}\int L{\frac{{d{I_d}}}{{dt}}} dt \approx \frac{R}{L}\int {{K_p}({I_{dref}} - {I_d})} dt = {K_i}\int {({I_{dref}} - {I_d})} dt\\ {K_i} = \frac{R}{L}{K_p} \end{array} \right.
⎩
⎨
⎧LdtdId≈LTsIdref−Id=Kp(Idref−Id)RId=LR∫LdtdIddt≈LR∫Kp(Idref−Id)dt=Ki∫(Idref−Id)dtKi=LRKp
同理,
L
d
I
q
d
t
+
R
I
q
L\frac{{d{I_q}}}{{dt}} +R{I_q}
LdtdIq+RIq也能得到类似的表达式。简单来看就是用PI控制器代替目标的电感电压,由于整流器本质上控制的是端口电压,而端口电压又是和电网电压,电感电压形成矢量三角形的,所以控制电感电压本质上就是控制电感电压。
将PI的表达式代入Ud和Uq表达式中,得到:
{
U
d
∗
=
U
s
m
+
I
q
ω
L
−
K
p
(
I
d
r
e
f
−
I
d
)
−
K
i
∫
(
I
d
r
e
f
−
I
d
)
d
t
U
d
∗
=
0
−
I
d
ω
L
−
K
p
(
I
q
r
e
f
−
I
q
)
−
K
i
∫
(
I
q
r
e
f
−
I
q
)
d
t
\left\{\begin{array}{l} U_{d}^{*}=U_{s m}+I_{q} \omega L-K_{p}\left(I_{d r e f}-I_{d}\right)-K_{i} \int\left(I_{d r e f}-I_{d}\right) d t \\ U_{d}^{*}=0-I_{d} \omega L-K_{p}\left(I_{q r e f}-I_{q}\right)-K_{i} \int\left(I_{q r e f}-I_{q}\right) d t \end{array}\right.
{Ud∗=Usm+IqωL−Kp(Idref−Id)−Ki∫(Idref−Id)dtUd∗=0−IdωL−Kp(Iqref−Iq)−Ki∫(Iqref−Iq)dt
控制框图如图所示。
2. 单相APFdq坐标系下特定次谐波抑制
单相APF的拓扑如图所示,变换器的目标是补偿电流ic等于iL中谐波电流相反数,由于is=iL+ic,因此is谐波电流为0,由于APF的容量有限,往往仅仅需要补偿低次谐波,就能够达到不错的补偿效果,所以特定次谐波补偿的需求应运而生。
参考上面的dq解耦控制,可以得到5次谐波抑制的控制表达式,需要注意的是,补偿谐波电流的参考值是负的谐波负载电流,所以
I
c
5
d
I_{c5d}
Ic5d的参考值是
−
I
L
5
d
-I_{L5d}
−IL5d。
{
U
5
d
∗
=
U
s
m
+
I
c
5
q
ω
L
−
K
p
(
−
I
L
5
d
−
I
c
5
d
)
+
K
i
∫
(
−
I
L
5
d
−
I
c
5
d
)
d
t
U
5
q
∗
=
0
−
I
c
5
d
ω
L
−
K
p
(
−
I
L
5
q
−
I
c
5
q
)
+
K
i
∫
(
I
L
5
q
+
I
c
5
q
)
d
t
\left\{\begin{array}{l} U_{5d}^{*}=U_{s m}+I_{c5q} \omega L-K_{p}\left(-I_{L5d}-I_{c5d}\right)+K_{i} \int\left(-I_{L5d}-I_{c5d}\right) d t \\ U_{5q}^{*}=0-I_{c5d} \omega L-K_{p}\left(-I_{L5q}-I_{c5q}\right)+K_{i} \int\left(I_{L5q}+I_{c5q}\right) d t \end{array}\right.
{U5d∗=Usm+Ic5qωL−Kp(−IL5d−Ic5d)+Ki∫(−IL5d−Ic5d)dtU5q∗=0−Ic5dωL−Kp(−IL5q−Ic5q)+Ki∫(IL5q+Ic5q)dt
重新整理表达式,得:
{
U
5
d
∗
=
U
s
m
+
I
c
5
q
ω
L
+
K
p
(
I
L
5
d
+
I
c
5
d
)
+
K
i
∫
(
I
L
5
d
+
I
c
5
d
)
d
t
U
5
q
∗
=
0
−
I
c
5
d
ω
L
+
K
p
(
I
L
5
q
+
I
c
5
q
)
+
K
i
∫
(
I
L
5
q
+
I
c
5
q
)
d
t
\left\{\begin{array}{l} U_{5d}^{*}=U_{s m}+I_{c5q} \omega L+K_{p}\left(I_{L5d}+I_{c5d}\right)+K_{i} \int\left(I_{L5d}+I_{c5d}\right) d t \\ U_{5q}^{*}=0-I_{c5d} \omega L+K_{p}\left(I_{L5q}+I_{c5q}\right)+K_{i} \int\left(I_{L5q}+I_{c5q}\right) d t \end{array}\right.
{U5d∗=Usm+Ic5qωL+Kp(IL5d+Ic5d)+Ki∫(IL5d+Ic5d)dtU5q∗=0−Ic5dωL+Kp(IL5q+Ic5q)+Ki∫(IL5q+Ic5q)dt
由于
U
s
m
U_{s m}
Usm,
I
5
d
ω
L
I_{5d} \omega L
I5dωL,
I
5
q
ω
L
I_{5q} \omega L
I5qωL稳态时为常数,并且可以由后面的PI控制器生成,所以实际应用中可以忽略。得到以下表达式:
{
U
5
d
∗
=
K
p
(
I
L
5
d
+
I
c
5
d
)
+
K
i
∫
(
I
L
5
d
+
I
c
5
d
)
d
t
U
5
q
∗
=
K
p
(
I
L
5
q
+
I
c
5
q
)
+
K
i
∫
(
I
L
5
q
+
I
c
5
q
)
d
t
\left\{\begin{array}{l} U_{5d}^{*}=K_{p}\left(I_{L5d}+I_{c5d}\right)+K_{i} \int\left(I_{L5d}+I_{c5d}\right) d t \\ U_{5q}^{*}=K_{p}\left(I_{L5q}+I_{c5q}\right)+K_{i} \int\left(I_{L5q}+I_{c5q}\right) d t \end{array}\right.
{U5d∗=Kp(IL5d+Ic5d)+Ki∫(IL5d+Ic5d)dtU5q∗=Kp(IL5q+Ic5q)+Ki∫(IL5q+Ic5q)dt
再经过dq反变换,得到:
U
5
∗
=
U
5
d
∗
sin
(
5
ω
t
)
+
U
5
d
∗
cos
(
5
ω
t
)
U_5^* = U_{5d}^*\sin (5\omega t) + U_{5d}^*\cos (5\omega t)
U5∗=U5d∗sin(5ωt)+U5d∗cos(5ωt),将
U
5
∗
U_5^*
U5∗加入到基波电流控制的调制波中,就能实现谐波补偿的功能。
3. 单相整流器dq坐标系下特定次谐波抑制
单相整流器由于存在直流电压二次纹波,经过调制,从而使电网电流产生3次谐波电流,为了抑制电网电流的3次谐波,可以参考APF的谐波补偿思想,假设负载5次谐波电流为零,则补偿的五次谐波电流参考值为0,从而可得到以下表达式:
{
U
5
d
∗
=
K
p
(
0
+
I
c
5
d
)
+
K
i
∫
(
0
+
I
c
5
d
)
d
t
U
5
q
∗
=
K
p
(
0
+
I
c
5
q
)
+
K
i
∫
(
0
+
I
c
5
q
)
d
t
\left\{\begin{array}{l} U_{5d}^{*}=K_{p}\left(0+I_{c5d}\right)+K_{i} \int\left(0+I_{c5d}\right) d t \\ U_{5q}^{*}=K_{p}\left(0+I_{c5q}\right)+K_{i} \int\left(0+I_{c5q}\right) d t \end{array}\right.
{U5d∗=Kp(0+Ic5d)+Ki∫(0+Ic5d)dtU5q∗=Kp(0+Ic5q)+Ki∫(0+Ic5q)dt
控制框图如图所示,将5次谐波抑制的调制波加入到基波控制器的调制波中,实现了5次谐波电流抑制的功能。
4. 三相系统特定次谐波抑制
和单相系统的唯一区别就是abc-dq,和dq-abc矩阵的差异,建模思想完全通用。
5. 参考文献
[1] 三相并联型有源电力滤波器电流重复控制_耿攀
[2] 并联型有源电力滤波器按容量比例分频段补偿并联控制策略_张树全
[3] 多同步旋转坐标系下指定次谐波电流控制_张树全
6. 反馈与建议
- 邮箱:upqcyjp@q63.com