单/三相dq解耦控制与特定次谐波抑制

news/2024/12/31 3:13:34/

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=usLdtdisRis
  将电压和电流写成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=usLdtdisRis。需要注意的是,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=UsmLdtdId+IqωLRIdUq=0LdtdIqIdωLRIq  上述的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. LdtdIdLTsIdrefId=Kp(IdrefId)RId=LRLdtdIddtLRKp(IdrefId)dt=Ki(IdrefId)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ωLKp(IdrefId)Ki(IdrefId)dtUd=0IdωLKp(IqrefIq)Ki(IqrefIq)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ωLKp(IL5dIc5d)+Ki(IL5dIc5d)dtU5q=0Ic5dωLKp(IL5qIc5q)+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=0Ic5dω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=U5dsin(5ωt)+U5dcos(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

http://www.ppmy.cn/news/1229649.html

相关文章

2023亚太杯数学建模思路 - 案例:FPTree-频繁模式树算法

文章目录 算法介绍FP树表示法构建FP树实现代码 建模资料 ## 赛题思路 (赛题出来以后第一时间在CSDN分享) https://blog.csdn.net/dc_sinor?typeblog 算法介绍 FP-Tree算法全称是FrequentPattern Tree算法,就是频繁模式树算法&#xff0c…

Linux下快速确定目标服务器支持哪些协议和密码套件

实现原理是利用TLS协议的特点和握手过程来进行测试和解析响应来确定目标服务器支持哪些TLS协议和密码套件。 在TLS握手过程中,客户端和服务器会协商并使用相同的TLS协议版本和密码套件来进行通信。通过发送特定的握手请求并分析响应,可以确定目标服务器…

指令系统、流水线

指令系统 分类 寻址方式 设计 能够改变控制流的指令:分支、跳转、过程调用、过程返回 操作码设计 MIPS 流水线 MIPS流水线 改进后 取指(IF) 译码(ID) 执行(EX) 存储器访问 寄存器-寄存器 A…

【51单片机系列】C51基础

本文内容是关于C51语言的基础内容的,包括C51的数据类型、变量、运算符、函数以及reg52.h文件中的内容,有些与C中相同的内容没有记录在此,比如常量、某些变量、表达式、程序结构、数组等没有涉及。 文章目录 C51的数据类型1. C51中的基本数据类…

arcgis中投影文件(.prj)和地理转换文件(.gtf)存储路径

1、投影文件(自定义的.prj)的存储路径 C:\Users\14635\AppData\Roaming\ESRI\Desktop10.5\ArcMap\Coordinate Systems 2、地理转换文件(.gtf)--自定义 C:\Users\14635\AppData\Roaming\ESRI\Desktop10.5\ArcToolbox\CustomTransfo…

【Zabbix监控二】之zabbix自定义监控内容案例(自动发现、自动注册)

一、自定义监控内容 案例:自定义监控客户端服务器登录的人数 需求:限制登录人数不超过3个人,超过5个人就发出报警 1、在客户端创建自定义key 明确需要执行的linux命令 创建zabbix监控项配置文件,用于自定义Key #在zabbix的配…

类型体系与基本数据类型(第五节)

目录 前言 一、列表 1.1 Batch模板 如何定位子矩阵中的元素 1.2 Array模板 1.3 重复与Duplicate模板 总结 前言 一个深度学习框架的初步实现为例,讨论如何在一个相对较大的项目中深入应用元编程,为系统优化提供更多的可能。 以下内容结合书中原…

mysql5.7安装

下载地址: https://downloads.mysql.com/archives/community/ 引用: https://blog.csdn.net/antma/article/details/79224279