CN110429839A - 一种三相电压型pwm整流器的分数阶建模方法 - Google Patents
一种三相电压型pwm整流器的分数阶建模方法 Download PDFInfo
- Publication number
- CN110429839A CN110429839A CN201910783755.6A CN201910783755A CN110429839A CN 110429839 A CN110429839 A CN 110429839A CN 201910783755 A CN201910783755 A CN 201910783755A CN 110429839 A CN110429839 A CN 110429839A
- Authority
- CN
- China
- Prior art keywords
- phase
- fractional order
- coordinate system
- pwm rectifier
- type pwm
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000004088 simulation Methods 0.000 claims abstract description 63
- 230000003068 static effect Effects 0.000 claims abstract description 51
- 230000001360 synchronised effect Effects 0.000 claims abstract description 35
- 239000003990 capacitor Substances 0.000 claims abstract description 27
- 238000013178 mathematical model Methods 0.000 claims description 65
- 230000009466 transformation Effects 0.000 claims description 29
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000007935 neutral effect Effects 0.000 claims description 8
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 6
- 230000001939 inductive effect Effects 0.000 claims description 6
- 230000001131 transforming effect Effects 0.000 claims description 6
- POXUQBFHDHCZAD-MHTLYPKNSA-N (2r)-2-[(4s)-2,2-dimethyl-1,3-dioxolan-4-yl]-3,4-dihydroxy-2h-furan-5-one Chemical compound O1C(C)(C)OC[C@H]1[C@@H]1C(O)=C(O)C(=O)O1 POXUQBFHDHCZAD-MHTLYPKNSA-N 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 238000005094 computer simulation Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 16
- 230000004044 response Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003989 dielectric material Substances 0.000 description 1
- 238000010292 electrical insulation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M7/00—Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
- H02M7/02—Conversion of ac power input into dc power output without possibility of reversal
- H02M7/04—Conversion of ac power input into dc power output without possibility of reversal by static converters
- H02M7/12—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
- H02M7/21—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal
- H02M7/217—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only
- H02M7/219—Conversion of ac power input into dc power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only in a bridge configuration
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Rectifiers (AREA)
Abstract
本发明涉及电力电子系统建模技术领域，具体涉及一种三相电压型PWM整流器的分数阶建模方法。本发明考虑电感和电容的分数阶特性，提出了在三相静止（a，b，c）坐标系、两相静止（D，Q）坐标系以及同步旋转（d，q）坐标系中建立三相电压型PWM整流器分数阶模型的方法及其对应的在三相静止（a，b，c）坐标系、两相静止（D，Q）坐标系以及同步旋转（d，q）坐标系中的三相电压型PWM整流器分数阶模型，可以方便的用于三相电压型PWM整流器及其分数阶控制器的设计与分析，改善运行性能与控制效果。其中在三相静止（a，b，c）坐标系中的三相电压型PWM整流器分数阶模型，可以实现三相电压型PWM整流器的高精度波形仿真。
Description
技术领域
本发明涉及电力电子系统建模技术领域，具体涉及一种三相电压型PWM整流器的分数阶建模方法。
背景技术
分数阶微积分是整数阶微积分的微积分阶次从整数扩展到非整数(分数、无理数、复数)后得到的。自然界本质上是分数维的，采用分数阶微积分可以建立系统准确或更准确的分数阶数学模型。
文献Westerlund S,Ekstam L.Capacitor theory[J].IEEE TRANSACTIONS ONDIELECTRICS AND ELECTRICAL INSULATION和文献Westerlund S.Dead Matter HasMemory[M].Kalmar,Sweden:Causal Consulting的研究表明：电容和电感本质上都是分数阶的，整数阶的电感和电容在实际工程中并不存在，以往用于描述电感和电容的电特性的整数阶模型是不够准确的，甚至可能是错误的，须建立相应的分数阶模型来描述它们的电特性。
开关型电力电子变换器在电能转换、存储与传输中发挥着越来越重要的作用，电感和电容是开关型电力电子变换器的关键元器件。传统的开关型电力电子变换器建模都是把电感和电容当作整数阶元器件，从而建立系统的整数阶数学模型。近年来，逐渐有学者在考虑电感和电容分数阶特性的基础上，研究开关型电力电子变换器的分数阶建模问题。文献王发强,马西奎.电感电流连续模式下Boost变换器的分数阶建模与仿真分析[J]、王发强,马西奎.基于分数阶微积分的电感电流断续模式下Boost变换器的建模与分析[J]、Faqiang Wang,Xikui Ma.Transfer function modeling and analysis of the open-loop Buck converter using the fractional calculus[J]、Faqiang Wang,XikuiMa.Fractional order Buck–Boost converter in CCM:modelling,analysis andsimulations[J]、郑征,马方军,韦延方.单相PWM整流器分数阶建模与仿真分析[J]、韦延方,郑征,王新良等.基于分数阶微积分的PWM整流器建模方法的研究结果表明，在开关型电力电子变换器建模时应考虑电感和电容的分数阶特性，采用分数阶模型能更真实地描述变换器的实际动力学行为。其中，文献郑征,马方军,韦延方.单相PWM整流器分数阶建模与仿真分析[J]、韦延方,郑征,王新良等.基于分数阶微积分的PWM整流器建模方法考虑电感和电容的分数阶特性，初步研究了PWM整流器的分数阶建模方法，但仅讨论了单相电压型PWM整流器在静止坐标系下的建模方法，没有解决三相电压型PWM整流器(Three PhaseVoltage Source PWM Rectifier，简称三相VSR)的分数阶建模问题。
发明内容
为了解决传统三相电压型PWM整流器建模方法没有考虑电感、电容的分数阶特性，存在模型不够准确的问题，本发明提供了一种三相电压型PWM整流器的分数阶建模方法，具体技术方案如下：
一种三相电压型PWM整流器的分数阶建模方法，包括以下步骤：
步骤1：建立三相电压型PWM整流器的电感和电容的分数阶数学模型：
采用Caputo型分数阶微积分定义，建立三相电压型PWM整流器的电感和电容的通用分数阶数学模型，具体如下：
其中，ul为电感电压，il为电感电流，uC为电容电压，iC为电容电流，L、C分别为电感值和电容值，α、β分别为电感阶数和电容阶数，并且满足0＜α,β＜2；
步骤2：建立三相静止坐标系中三相电压型PWM整流器的交流侧分数阶数学模型：
采用基尔霍夫电压、电流定律，并考虑电感和电容的分数阶特性，建立三相电压型PWM整流器的a、b、c三相回路方程，具体如下：
其中，ua、ub、uc分别为电网的三相电压；ia、ib、ic为三相输入相电流，由于三相电压型PWM整流器的直流侧负极点N一般不与交流侧中性点0相连，故无论三相电网电压是否平衡，均有ia+ib+ic＝0；La、Lb、Lc为三相进线滤波电抗器的电感，αa、αb、αc为三相进线滤波电抗器的电感阶数，且0＜αa,αb,αc＜2；Ra、Rb、Rc为包括进线滤波电抗器等效电阻和功率开关管损耗等效电阻在内的三相进线电阻；Udc为直流侧电压；vN0为直流侧负极点N与三相电网中性点0之间的电压；sa、sb、sb为各相桥臂的单极性二值逻辑开关函数，且定义上桥臂元件导通时值为1，下桥臂元件导通时值为0；
步骤3：将ia+ib+ic＝0代入式(3)，整理后vN0可表示为：
步骤4：将式(4)代入式(3)，得到用开关函数描述的三相静止坐标系中三相电压型PWM整流器的交流侧分数阶数学模型：
步骤5：建立三相静止坐标系中三相电压型PWM整流器的直流侧分数阶数学模型：
对三相电压型PWM整流器的直流侧电容正极节点处应用基尔霍夫电流定律，得到在三相静止坐标系中的三相电压型PWM整流器直流侧分数阶数学模型，具体如下：
其中，C和β分别为直流侧滤波电容及其阶数，且0＜β＜2；idc为直流侧电流，idc＝iasa+ibsb+icsc；iL为负载电流，iL＝(Udc-eL)/RL，RL为负载电阻，eL为负载等效直流电动势；
步骤6：建立三相静止坐标系中三相电压型PWM整流器的分数阶数学模型：
联立式(5)和式(6)，得到三相电压型PWM整流器三相静止坐标系分数阶数学模型表达式为：
步骤7：对式(7)进行整理，可以得到三相电压型PWM整流器的三相静止坐标系分数阶数学模型表达式为：
其中：L＝diag[La,Lb,Lc]；R＝diag[Ra,Rb,Rc]；
步骤8：在零初始条件下，对式(7)进行Laplace变化，得三相电压型PWM整流器的三相静止坐标系分数阶数学模型的s域表达式为：
步骤9：在MATLAB/Simulink仿真平台，采用改进的Oustaloup分数阶微积分滤波器算法模拟分数阶电感和分数阶电容，建立三相电压型PWM整流器在三相静止坐标系中的分数阶数字仿真电路模型；
步骤10：运行仿真，当误差分析满足要求时，记录仿真结果
优选地，还包括以下步骤：
步骤11：建立两相静止坐标系中三相电压型PWM整流器的分数阶数学模型：
忽略各相进线电抗器的电阻、电感及其阶数间的微小差异，即Ra＝Rb＝Rc＝R，La＝Lb＝Lc＝L，αa＝αb＝αc＝α；
采用坐标变换关系对式(8)所示的三相静止坐标系中的三相电压型PWM整流器分数阶数学模型进行变换，可以得到三相电压型PWM整流器的两相静止坐标系分数阶数学模型，该步骤的坐标变换关系为：
iabc＝C2s/3siDQ，uabc＝C2s/3suDQ，sabc＝C2s/3ssDQ；
三相电压型PWM整流器的两相静止坐标系分数阶数学模型为：
式中：
其中，为两相静止坐标系到三相静止坐标系的变换矩阵，C3s/2s为三相静止坐标系到两相静止坐标系的变换矩阵；uD、uQ分别为三相电网电压的D、Q分量；iD、iQ分别为输入电流的D、Q分量；sD、sQ分别为开关函数的D、Q分量；
步骤12：对式(10)进行Laplace变换，得三相电压型PWM整流器的两相静止坐标系分数阶数学模型的s域表达式，具体如下：
步骤13：在MATLAB/Simulink仿真平台，采用改进的Oustaloup分数阶微积分滤波器算法模拟分数阶电感和分数阶电容，建立三相电压型PWM整流器在两相静止坐标系中的分数阶数字仿真电路模型；
步骤14：运行仿真，当误差分析满足要求时，记录仿真结果。
优选地，还包括以下步骤：
步骤15：建立同步旋转坐标系中三相电压型PWM整流器的分数阶数学模型：
采用如下坐标变换关系：
iDQ＝C2r/2sidq，uDQ＝C2r/2sudq，sDQ＝C2r/2ssdq；
将式(10)所示的三相电压型PWM整流器的两相静止坐标系分数阶数学模型变换到同步旋转坐标系：
式中：
其中，为同步旋转坐标系到两相静止坐标系的变换矩阵，C2s/2r为两相静止坐标系到同步旋转坐标系的变换矩阵；ud、uq分别为三相电网电压的d、q分量；id、iq分别为输入电流的d、q分量；sd、sq分别为开关函数的d、q分量；θ为d轴与D轴之间的夹角，θ＝ωt+θ0，θ0为t＝0时刻的初始相位角，ω为电网电压角频率；
基于Caputo型分数阶导数定义，当0＜α＜2时，有：
将式(13)代入式(12)，整理后得到三相电压型PWM整流器的同步旋转坐标系分数阶数学模型：
步骤16：对式(14)进行Laplace变换，得到三相电压型PWM整流器的同步旋转坐标系分数阶数学模型的s域表达式为：
步骤17：在MATLAB/Simulink仿真平台，采用改进的Oustaloup分数阶微积分滤波器算法模拟分数阶电感和分数阶电容，建立三相电压型PWM整流器在同步旋转坐标系中的分数阶数字仿真电路模型；
步骤18：运行仿真，当误差分析满足要求时，记录仿真结果。
本发明的有益效果为：
1、本发明考虑电感和电容的分数阶特性，提出了在三相静止(a，b，c)坐标系中建立三相电压型PWM整流器分数阶模型的方法，同时提供了在三相静止(a，b，c)坐标系中的三相电压型PWM整流器分数阶模型，该模型可以实现三相电压型PWM整流器的高精度波形仿真。
2、本发明进一步提出了在两相静止(D，Q)坐标系以及同步旋转(d，q)坐标系中建立三相电压型PWM整流器分数阶模型的方法；同时提供了在两相静止(D，Q)坐标系以及同步旋转(d，q)坐标系中的三相电压型PWM整流器分数阶模型。
3、本发明提供的在三相静止(a，b，c)坐标系、两相静止(D，Q)坐标系以及同步旋转(d，q)坐标系中的三相电压型PWM整流器分数阶模型，考虑了电感和电容的分数阶特性，是三相电压型PWM整流器模型的通用表达式，可以替代三相电压型PWM整流器的传统整数阶模型，也可以方便的用于三相电压型PWM整流器及其分数阶控制器的设计与分析，改善运行性能与控制效果。
附图说明
图1为具有分数阶特性的三相电压型PWM整流器主电路结构图；
图2为三相静止坐标系中三相电压型PWM整流器的分数阶模型结构图；
图3为两相静止坐标系中三相电压型PWM整流器的分数阶模型结构图；
图4为同步旋转坐标系中三相电压型PWM整流器的分数阶模型结构图；
图5为三相电压型PWM整流器三相静止坐标系分数阶数字仿真模型；
图6为三相电压型PWM整流器两相静止坐标系分数阶数字仿真模型；
图7为三相电压型PWM整流器同步旋转坐标系分数阶数字仿真模型；
图8三相电压型PWM整流器三相静止坐标系分数阶模型的直流输出电压仿真波形图；
图9为三相电压型PWM整流器三相静止坐标系分数阶模型的a相电感电流仿真波形图；
图10为三相电压型PWM整流器两相静止坐标系分数阶模型的直流输出电压仿真波形图；
图11为三相电压型PWM整流器两相静止坐标系分数阶模型的D轴电感电流仿真波形图；
图12为三相电压型PWM整流器两相静止坐标系分数阶模型的Q轴电感电流仿真波形图；
图13为三相电压型PWM整流器同步旋转坐标系分数阶模型的直流输出电压仿真波形图；
图14为三相电压型PWM整流器同步旋转坐标系分数阶模型的d轴电感电流仿真波形图；
图15为三相电压型PWM整流器同步旋转坐标系分数阶模型的q轴电感电流仿真波形图。
具体实施方式
为了更好的理解本发明，下面结合附图和具体实施例对本发明作进一步说明：
图1是本发明所述的三相电压型PWM整流器主电路拓扑结构图。图中，ua、ub、uc为三相电网的相电压；va、vb、vc为三相电压型PWM整流器的交流侧输入相电压；ia、ib、ic为三相输入相电流由于三相电压型PWM整流器的直流侧负极点N一般不与交流侧中性点0相连，故无论三相电网电压是否平衡，均有ia+ib+ic＝0；La、Lb、Lc为三相进线滤波电抗器的电感，其阶数分别为αa、αb、αc，且0＜αa,αb,αc＜2；Rla、Rlb、Rlc为三相进线滤波电抗器等效电阻；Tl～T6为功率开关元件；Dl～D6为续流二极管；idc为直流侧电流；Udc为直流侧电压；C为直流侧滤波电容，其阶数为β，且0＜β＜2；iL为负载电流；RL为负载电阻；eL为负载等效直流电动势；0为三相电网中性点，P为直流侧正极点，N为直流侧负极点。
本发明提供的一种三相电压型PWM整流器的分数阶建模方法，包括以下步骤：
步骤1：建立三相电压型PWM整流器的电感和电容的分数阶数学模型：
采用Caputo型分数阶微积分定义，建立三相电压型PWM整流器的电感和电容的通用分数阶数学模型，具体如下：
其中，ul为电感电压，il为电感电流，uC为电容电压，iC为电容电流，L、C分别为电感值和电容值，α、β分别为电感阶数和电容阶数，并且满足0＜α,β＜2；
步骤2：建立三相静止(a，b，c)坐标系中三相电压型PWM整流器的交流侧分数阶数学模型：
采用基尔霍夫电压、电流定律，并考虑电感和电容的分数阶特性，建立三相电压型PWM整流器的a、b、c三相回路方程，具体如下：
其中，ua、ub、uc分别为电网的三相电压；ia、ib、ic为三相输入相电流，由于三相电压型PWM整流器的直流侧负极点N一般不与交流侧中性点0相连，故无论三相电网电压是否平衡，均有ia+ib+ic＝0；La、Lb、Lc为三相进线滤波电抗器的电感，αa、αb、αc为三相进线滤波电抗器的电感阶数，且0＜αa,αb,αc＜2；Ra、Rb、Rc为包括进线滤波电抗器等效电阻和功率开关管损耗等效电阻在内的三相进线电阻；Udc为直流侧电压；vN0为直流侧负极点N与三相电网中性点0之间的电压；sa、sb、sb为各相桥臂的单极性二值逻辑开关函数，且定义上桥臂元件导通时值为1，下桥臂元件导通时值为0；
步骤3：将ia+ib+ic＝0代入式(3)，整理后vN0可表示为：
步骤4：将式(4)代入式(3)，得到用开关函数描述的三相静止坐标系中三相电压型PWM整流器的交流侧分数阶数学模型：
步骤5：建立三相静止坐标系中三相电压型PWM整流器的直流侧分数阶数学模型：
对三相电压型PWM整流器的直流侧电容正极节点处应用基尔霍夫电流定律，得到在三相静止坐标系中的三相电压型PWM整流器直流侧分数阶数学模型，具体如下：
其中，C和β分别为直流侧滤波电容及其阶数，且0＜β＜2；idc为直流侧电流，idc＝iasa+ibsb+icsc；iL为负载电流，iL＝(Udc-eL)/RL，RL为负载电阻，eL为负载等效直流电动势；
步骤6：建立三相静止坐标系中三相电压型PWM整流器的分数阶数学模型：
联立式(5)和式(6)，得到三相电压型PWM整流器三相静止坐标系分数阶数学模型表达式为：
步骤7：对式(7)进行整理，可以得到三相电压型PWM整流器的三相静止坐标系分数阶数学模型表达式为：
其中：L＝diag[La,Lb,Lc]；R＝diag[Ra,Rb,Rc]；
步骤8：在零初始条件下，对式(7)进行Laplace变化，得三相电压型PWM整流器的三相静止坐标系分数阶数学模型的s域表达式为：
三相电压型PWM整流器在三相静止坐标系中的分数阶模型结构如图2所示。
两相静止(D，Q)坐标系中三相电压型PWM整流器分数阶建模与模型：
步骤9：建立两相静止(D，Q)坐标系中三相电压型PWM整流器的分数阶数学模型：
忽略各相进线电抗器的电阻、电感及其阶数间的微小差异，即Ra＝Rb＝Rc＝R，La＝Lb＝Lc＝L，αa＝αb＝αc＝α；
采用坐标变换关系对式(8)所示的三相静止坐标系中的三相电压型PWM整流器分数阶数学模型进行变换，可以得到三相电压型PWM整流器的两相静止坐标系分数阶数学模型，该步骤的坐标变换关系为：
iabc＝C2s/3siDQ，uabc＝C2s/3suDQ，sabc＝C2s/3ssDQ；
三相电压型PWM整流器的两相静止坐标系分数阶数学模型为：
式中：
其中，为两相静止坐标系到三相静止坐标系的变换矩阵，C3s/2s为三相静止坐标系到两相静止坐标系的变换矩阵；uD、uQ分别为三相电网电压的D、Q分量；iD、iQ分别为输入电流的D、Q分量；sD、sQ分别为开关函数的D、Q分量；
步骤10：对式(10)进行Laplace变换，得三相电压型PWM整流器的两相静止坐标系分数阶数学模型的s域表达式，具体如下：
三相电压型PWM整流器在两相静止坐标系中的分数阶模型结构如图3所示。
同步旋转(d，q)坐标系中三相电压型PWM整流器分数阶建模与模型：
步骤11：建立同步旋转坐标系中三相电压型PWM整流器的分数阶数学模型：
采用如下坐标变换关系：
iDQ＝C2r/2sidq，uDQ＝C2r/2sudq，sDQ＝C2r/2ssdq；
将式(10)所示的三相电压型PWM整流器的两相静止坐标系分数阶数学模型变换到同步旋转坐标系：
式中：
其中，为同步旋转坐标系到两相静止坐标系的变换矩阵，C2s/2r为两相静止坐标系到同步旋转坐标系的变换矩阵；ud、uq分别为三相电网电压的d、q分量；id、iq分别为输入电流的d、q分量；sd、sq分别为开关函数的d、q分量；θ为d轴与D轴之间的夹角，θ＝ωt+θ0，θ0为t＝0时刻的初始相位角，ω为电网电压角频率；
基于Caputo型分数阶导数定义，当0＜α＜2时，有：
将式(13)代入式(12)，整理后得到三相电压型PWM整流器的同步旋转坐标系分数阶数学模型：
步骤12：对式(14)进行Laplace变换，得到三相电压型PWM整流器的同步旋转坐标系分数阶数学模型的s域表达式为：
三相电压型PWM整流器在同步旋转坐标系中的分数阶模型结构如图4所示。
仿真与分析：
步骤13：在MATLAB/Simulink仿真平台，采用改进的Oustaloup分数阶微积分滤波器算法模拟分数阶电感和分数阶电容，分别建立三相电压型PWM整流器在三相静止、两相静止和同步旋转坐标系中的分数阶数字仿真电路模型；
步骤14：运行仿真，当误差分析满足要求时，记录仿真结果。
对于图1所示的三相电压型PWM整流器，令基本电路参数为：三相电网的相电压有效值为220V，且三相电压平衡；电网电压频率50Hz；三相进线滤波电抗器参数相同，电感L＝4mH，等效电阻R＝0.2Ω；直流侧电容C＝4700μH；负载电阻RL＝5Ω，负载等效直流电动势eL＝470V。用SPWM波对三相VSR进行开环控制，开关频率fs＝5kHz。
为了反映分数阶电感和分数阶电容对三相电压型PWM整流器运行特性的影响，设置分数阶电感阶数和分数阶电容阶数为四组不同的值进行仿真、对比，参数如表1所示。
表1分数阶电感阶数和分数阶电容阶数取值表
根据图2，建立三相静止坐标系中三相电压型PWM整流器的分数阶数字仿真模型如图5所示。图5中，“Frac L”和“Frac C”分别为在三相静止坐标系中的分数阶电感和分数阶电容仿真模块，“SPWM”表示正弦脉宽调制波形生成模块。
根据图3，建立两相静止坐标系中三相电压型PWM整流器的分数阶数字仿真模型如图6所示。图6中，“Frac LDQ”和“Frac CDQ”分别为在两相静止坐标系中的分数阶电感和分数阶电容仿真模块，“Clark”模块是实现三相静止坐标系到两相静止坐标系之间坐标变换的仿真模块。
根据图4，建立三相电压型PWM整流器同步旋转坐标系分数阶模型的分数阶数字仿真模型如图7所示。图7中，“Frac Ldq”和“Frac Cdq”分别为在同步旋转坐标系中的分数阶电感和分数阶电容仿真模块，“Park”模块是实现三相静止坐标系到同步旋转坐标系之间坐标变换的仿真模块。
运行仿真，当误差分析满足要求时，记录仿真结果。
其中，图8为三相电压型PWM整流器三相静止坐标系分数阶模型的直流输出电压波形仿真波形，图9为三相电压型PWM整流器三相静止坐标系分数阶模型的a相电感电流波形仿真波形。由图8-图9可知，当α＝0.8、β＝0.8时，系统响应相对较快；当α＝1.0，β＝1.0时，系统响应相对较慢；当α＝0.8、β＝1.0时，直流电压的稳定值相对较大；当α＝1.0、β＝0.8时，直流电压的稳定值相对较小。
图10为三相电压型PWM整流器两相静止坐标系分数阶模型的直流输出电压仿真波形图，图11为三相电压型PWM整流器两相静止坐标系分数阶模型的D轴电感电流仿真波形图，图12为三相电压型PWM整流器两相静止坐标系分数阶模型的Q轴电感电流仿真波形图。对比图8和图10可知，三相电压型PWM整流器的三相静止坐标系分数阶模型与两相静止坐标分数阶模型的直流输出特性基本相同。
图13为三相电压型PWM整流器同步旋转坐标系分数阶模型的直流输出电压仿真波形图，图13为三相电压型PWM整流器同步旋转坐标系分数阶模型的直流输出电压仿真波形图，图15为三相电压型PWM整流器同步旋转坐标系分数阶模型的q轴电感电流仿真波形图。对比图8、图10和图13可知，三相电压型PWM整流器的三相静止坐标系分数阶模型、两相静止坐标分数阶模型和同步旋转坐标系模型的直流输出特性基本相同。
仿真结果验证了本文所提的三相电压型PWM整流器分数阶建模方法与所建模型的正确性。
本发明不局限于以上所述的具体实施方式，以上所述仅为本发明的较佳实施案例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (3)
1.一种三相电压型PWM整流器的分数阶建模方法，其特征在于：包括以下步骤：
步骤1：建立三相电压型PWM整流器的电感和电容的分数阶数学模型：
采用Caputo型分数阶微积分定义，建立三相电压型PWM整流器的电感和电容的通用分数阶数学模型，具体如下：
其中，ul为电感电压，il为电感电流，uC为电容电压，iC为电容电流，L、C分别为电感值和电容值，α、β分别为电感阶数和电容阶数，并且满足0＜α,β＜2；
步骤2：建立三相静止坐标系中三相电压型PWM整流器的交流侧分数阶数学模型：
采用基尔霍夫电压、电流定律，并考虑电感和电容的分数阶特性，建立三相电压型PWM整流器的a、b、c三相回路方程，具体如下：
其中，ua、ub、uc分别为电网的三相电压；ia、ib、ic为三相输入相电流，由于三相电压型PWM整流器的直流侧负极点N一般不与交流侧中性点0相连，故无论三相电网电压是否平衡，均有ia+ib+ic＝0；La、Lb、Lc为三相进线滤波电抗器的电感，αa、αb、αc为三相进线滤波电抗器的电感阶数，且0＜αa,αb,αc＜2；Ra、Rb、Rc为包括进线滤波电抗器等效电阻和功率开关管损耗等效电阻在内的三相进线电阻；Udc为直流侧电压；vN0为直流侧负极点N与三相电网中性点0之间的电压；sa、sb、sb为各相桥臂的单极性二值逻辑开关函数，且定义上桥臂元件导通时值为1，下桥臂元件导通时值为0；
步骤3：将ia+ib+ic＝0代入式(3)，整理后vN0可表示为：
步骤4：将式(4)代入式(3)，得到用开关函数描述的三相静止坐标系中三相电压型PWM整流器的交流侧分数阶数学模型：
步骤5：建立三相静止坐标系中三相电压型PWM整流器的直流侧分数阶数学模型：
对三相电压型PWM整流器的直流侧电容正极节点处应用基尔霍夫电流定律，得到在三相静止坐标系中的三相电压型PWM整流器直流侧分数阶数学模型，具体如下：
其中，C和β分别为直流侧滤波电容及其阶数，且0＜β＜2；idc为直流侧电流，idc＝iasa+ibsb+icsc；iL为负载电流，iL＝(Udc-eL)/RL，RL为负载电阻，eL为负载等效直流电动势；
步骤6：建立三相静止坐标系中三相电压型PWM整流器的分数阶数学模型：
联立式(5)和式(6)，得到三相电压型PWM整流器三相静止坐标系分数阶数学模型表达式为：
步骤7：对式(7)进行整理，可以得到三相电压型PWM整流器的三相静止坐标系分数阶数学模型表达式为：
其中：L＝diag[La,Lb,Lc]；R＝diag[Ra,Rb,Rc]；
步骤8：在零初始条件下，对式(7)进行Laplace变化，得三相电压型PWM整流器的三相静止坐标系分数阶数学模型的s域表达式为：
步骤9：在MATLAB/Simulink仿真平台，采用改进的Oustaloup分数阶微积分滤波器算法模拟分数阶电感和分数阶电容，建立三相电压型PWM整流器在三相静止坐标系中的分数阶数字仿真电路模型；
步骤10：运行仿真，当误差分析满足要求时，记录仿真结果。
2.根据权利要求1所述的一种三相电压型PWM整流器的分数阶建模方法，其特征在于：还包括以下步骤：
步骤11：建立两相静止坐标系中三相电压型PWM整流器的分数阶数学模型：
忽略各相进线电抗器的电阻、电感及其阶数间的微小差异，即Ra＝Rb＝Rc＝R，La＝Lb＝Lc＝L，αa＝αb＝αc＝α；
采用坐标变换关系对式(8)所示的三相静止坐标系中的三相电压型PWM整流器分数阶数学模型进行变换，可以得到三相电压型PWM整流器的两相静止坐标系分数阶数学模型，该步骤的坐标变换关系为：
iabc＝C2s/3siDQ，uabc＝C2s/3suDQ，sabc＝C2s/3ssDQ；
三相电压型PWM整流器的两相静止坐标系分数阶数学模型为：
式中：
其中，为两相静止坐标系到三相静止坐标系的变换矩阵，C3s/2s为三相静止坐标系到两相静止坐标系的变换矩阵；uD、uQ分别为三相电网电压的D、Q分量；iD、iQ分别为输入电流的D、Q分量；sD、sQ分别为开关函数的D、Q分量；
步骤12：对式(10)进行Laplace变换，得三相电压型PWM整流器的两相静止坐标系分数阶数学模型的s域表达式，具体如下：
步骤13：在MATLAB/Simulink仿真平台，采用改进的Oustaloup分数阶微积分滤波器算法模拟分数阶电感和分数阶电容，建立三相电压型PWM整流器在两相静止坐标系中的分数阶数字仿真电路模型；
步骤14：运行仿真，当误差分析满足要求时，记录仿真结果。
3.根据权利要求1或权利要求2所述的一种三相电压型PWM整流器的分数阶建模方法，其特征在于：还包括以下步骤：
步骤15：建立同步旋转坐标系中三相电压型PWM整流器的分数阶数学模型：
采用如下坐标变换关系：
iDQ＝C2r/2sidq，uDQ＝C2r/2sudq，sDQ＝C2r/2ssdq；
将式(10)所示的三相电压型PWM整流器的两相静止坐标系分数阶数学模型变换到同步旋转坐标系：
式中：
其中，为同步旋转坐标系到两相静止坐标系的变换矩阵，C2s/2r为两相静止坐标系到同步旋转坐标系的变换矩阵；ud、uq分别为三相电网电压的d、q分量；id、iq分别为输入电流的d、q分量；sd、sq分别为开关函数的d、q分量；θ为d轴与D轴之间的夹角，θ＝ωt+θ0，θ0为t＝0时刻的初始相位角，ω为电网电压角频率；
基于Caputo型分数阶导数定义，当0＜α＜2时，有：
将式(13)代入式(12)，整理后得到三相电压型PWM整流器的同步旋转坐标系分数阶数学模型：
步骤16：对式(14)进行Laplace变换，得到三相电压型PWM整流器的同步旋转坐标系分数阶数学模型的s域表达式为：
步骤17：在MATLAB/Simulink仿真平台，采用改进的Oustaloup分数阶微积分滤波器算法模拟分数阶电感和分数阶电容，建立三相电压型PWM整流器在同步旋转坐标系中的分数阶数字仿真电路模型；
步骤18：运行仿真，当误差分析满足要求时，记录仿真结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910783755.6A CN110429839B (zh) | 2019-08-23 | 2019-08-23 | 一种三相电压型pwm整流器的分数阶建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910783755.6A CN110429839B (zh) | 2019-08-23 | 2019-08-23 | 一种三相电压型pwm整流器的分数阶建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110429839A true CN110429839A (zh) | 2019-11-08 |
| CN110429839B CN110429839B (zh) | 2021-07-09 |
Family
ID=68415697
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910783755.6A Active CN110429839B (zh) | 2019-08-23 | 2019-08-23 | 一种三相电压型pwm整流器的分数阶建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110429839B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111478565A (zh) * | 2020-05-23 | 2020-07-31 | 西安科技大学 | Vienna整流器的高次谐波抑制控制器的设计方法 |
| CN112630497A (zh) * | 2020-12-03 | 2021-04-09 | 苏州英威腾电力电子有限公司 | 一种pwm整流器的自检方法、装置及系统 |
| CN116702575A (zh) * | 2023-08-08 | 2023-09-05 | 石家庄铁道大学 | 一种柔性互联装置分数阶建模与优化控制方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105608266A (zh) * | 2015-12-10 | 2016-05-25 | 河南理工大学 | 基于分数阶微积分的pwm整流器建模方法 |
-
2019
- 2019-08-23 CN CN201910783755.6A patent/CN110429839B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105608266A (zh) * | 2015-12-10 | 2016-05-25 | 河南理工大学 | 基于分数阶微积分的pwm整流器建模方法 |
Non-Patent Citations (1)
| Title |
|---|
| 姜春燕 等: "基于分数阶的航空用三相电压型整流器设计", 《航空计算技术》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111478565A (zh) * | 2020-05-23 | 2020-07-31 | 西安科技大学 | Vienna整流器的高次谐波抑制控制器的设计方法 |
| CN111478565B (zh) * | 2020-05-23 | 2021-10-26 | 西安科技大学 | Vienna整流器的高次谐波抑制控制器的设计方法 |
| CN112630497A (zh) * | 2020-12-03 | 2021-04-09 | 苏州英威腾电力电子有限公司 | 一种pwm整流器的自检方法、装置及系统 |
| CN112630497B (zh) * | 2020-12-03 | 2024-03-26 | 苏州英威腾电力电子有限公司 | 一种pwm整流器的自检方法、装置及系统 |
| CN116702575A (zh) * | 2023-08-08 | 2023-09-05 | 石家庄铁道大学 | 一种柔性互联装置分数阶建模与优化控制方法 |
| CN116702575B (zh) * | 2023-08-08 | 2023-10-20 | 石家庄铁道大学 | 一种柔性互联装置分数阶建模与优化控制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110429839B (zh) | 2021-07-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110429839B (zh) | 一种三相电压型pwm整流器的分数阶建模方法 | |
| Liu et al. | Modeling, analysis, and parameters design of LC-filter-integrated quasi-Z-source indirect matrix converter | |
| CN106208737B (zh) | 基于三次谐波注入矩阵变换器的模型预测电流控制方法 | |
| CN107732959B (zh) | 用于分布式光储并网系统的非线性微分平滑前馈控制方法 | |
| Palanisamy et al. | Maximum Boost Control for 7-level z-source cascaded h-bridge inverter | |
| Ellabban et al. | Model predictive control of a grid connected quasi-Z-source inverter | |
| Garcia et al. | Modeling and control design of the six-phase interleaved double dual boost converter | |
| Shieh et al. | Modelling and design of a reversible three-phase switching mode rectifier | |
| Cho et al. | Three phase sine wave voltage source inverter using the soft switched resonant poles | |
| CN113659905B (zh) | 基于时变扰动补偿的三电平发电系统模型预测控制方法 | |
| CN108964119B (zh) | 一种t型三电平并网逆变器共模电压预估方法 | |
| CN104362870B (zh) | 整流器的控制方法、装置及整流器 | |
| Balikci et al. | A three‐phase four‐wire static synchronous compensator with reduced number of switches for unbalanced loads | |
| CN110391726B (zh) | 单向三相星接可控整流器输入电流过零畸变的抑制方法 | |
| CN107846154B (zh) | Z源t型三电平逆变器的无源控制方法 | |
| Pan et al. | Fractional-order sliding mode control strategy for quasi-Z source photovoltaic grid-connected inverter | |
| CN103078348B (zh) | 一种无信号互联线并联的三相逆变器及其控制方法 | |
| CN110676860B (zh) | 一种基于扩展瞬时有功理论的快速预测不平衡控制方法 | |
| CN106160421B (zh) | 电流源型变流器直流电感的新型设计方法 | |
| Siwek et al. | Improvement of the torque control dynamics of the PMSM drive using the FOC-controlled simple boost QZSDMC converter | |
| Darwish et al. | Operation and control design of new Three-Phase inverters with reduced number of switches | |
| KR102288216B1 (ko) | 토크 예측 제어를 통한 pmsm 구동 제어 장치 및 방법 | |
| Kim et al. | DC link voltage control of reduced switch VSI-PWM rectifier/inverter system | |
| CN110867864A (zh) | 一种有源三次谐波注入矩阵变换器离网运行控制方法 | |
| Talavat et al. | Direct predictive control of asynchronous machine torque using matrix converter |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |