CN112836369A - 一种基于降阶模型的无线电能传输控制系统设计分析方法 - Google Patents
一种基于降阶模型的无线电能传输控制系统设计分析方法 Download PDFInfo
- Publication number
- CN112836369A CN112836369A CN202110131799.8A CN202110131799A CN112836369A CN 112836369 A CN112836369 A CN 112836369A CN 202110131799 A CN202110131799 A CN 202110131799A CN 112836369 A CN112836369 A CN 112836369A
- Authority
- CN
- China
- Prior art keywords
- equivalent
- reduced
- power transmission
- wireless power
- resonant
- 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
- 230000005540 biological transmission Effects 0.000 title claims abstract description 35
- 238000004458 analytical method Methods 0.000 title claims abstract description 27
- 239000003990 capacitor Substances 0.000 claims description 27
- 230000001131 transforming effect Effects 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000001052 transient effect Effects 0.000 description 5
- 230000007547 defect Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J50/00—Circuit arrangements or systems for wireless supply or distribution of electric power
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Computer Networks & Wireless Communication (AREA)
- Power Engineering (AREA)
- Transmitters (AREA)
Abstract
本发明涉及一种基于降阶模型的无线电能传输控制系统设计分析方法,该方法包括如下步骤:1)根据无线电能传输系统拓扑选取交流状态变量降阶量;2)对交流状态变量降阶量进行变换获取谐振元件的等效近似模型,合并谐振元件,获得降阶等效电路;3)根据降阶等效电路获得降阶状态空间模型;4)根据降阶状态空间模型计算闭环系统传递函数,完成控制参数设计和稳定性分析。与现有技术相比,本发明提出方法能够精确、直观的描述系统动态行为,有效降低了系统模型阶数,对于控制系统设计分析具有重要意义。
Description
技术领域
本发明涉及无线电能传输技术领域,尤其是涉及一种基于降阶模型的无线电能传输控制系统设计分析方法。
背景技术
无线电能传输系统的闭环控制是该技术应用推广的重要前提,然而,由于无线电能传输系统中储能元件较多,其电压/电流都为交流信号,需要分解为实部和虚部作为系统的状态变量,因此系统模型的阶数通常较高,使得系统的动态行为分析变得非常复杂,给系统控制设计及其稳定性分析带来了一系列问题。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于降阶模型的无线电能传输控制系统设计分析方法。
本发明的目的可以通过以下技术方案来实现:
一种基于降阶模型的无线电能传输控制系统设计分析方法,该方法包括如下步骤:
1)根据无线电能传输系统拓扑选取交流状态变量降阶量;
2)对交流状态变量降阶量进行变换获取谐振元件的等效近似模型,合并谐振元件,获得降阶等效电路;
3)根据降阶等效电路获得降阶状态空间模型;
4)根据降阶状态空间模型计算闭环系统传递函数,完成控制参数设计和稳定性分析。
优选地,步骤2)对交流状态变量降阶量进行变换获取谐振元件的等效近似模型具体为:
21)将交流状态变量降阶量进行傅里叶展开,获取基波傅里叶级数;
22)对交流状态变量降阶量基波傅里叶级数进行局部线性逼近,将带幅值扰动的交流状态变量降阶量波形近似为幅值线性变化的正弦波形;
23)基于近似的正弦波形确定谐振元件的等效近似模型。
优选地,步骤1)中若无线电能传输系统拓扑中含有LCL谐振结构或串联谐振结构,则选取谐振电容电压为降阶量。
优选地,步骤2)中谐振电容电压近似波形方程为:
优选地,步骤3)中谐振电容的等效近似模型确定为等效电感和等效电阻串联形式。
优选地,谐振电容的等效近似模型中等效电感大小记作L′,等效电阻大小记作R′:
其中,C为谐振电容的大小,j为虚部。
优选地,步骤1)中若无线电能传输系统拓扑中含有并联谐振结构,则选取谐振电感电流为降阶量。
优选地,步骤2)中谐振电感电流近似波形方程为:
优选地,步骤3)中谐振电感的等效近似模型确定为等效电容和等效电阻串联形式。
优选地,谐振电感的等效近似模型中等效电容大小记作C′,等效电阻大小记作R″:
其中,L为谐振电感的大小,j为虚部。
与现有技术相比,本发明具有如下优点:
本发明依据局部线性逼近原理实现系统冗余状态变量的近似,通过谐振回路的状态空间合并实现模型的降阶,简化了系统等效模型,便于动态分析和系统控制的实现,对于目前的高阶无线电能传输系统控制具有较高的应用价值和实际意义;
本发明提出等效简化电路模型,模型参数具有明确物理含义,易于理解、应用,可广泛适用于目前常用的各种谐振拓扑,不受参数条件约束;
本发明设计分析方法可精确描述实际全阶模型的动态行为,在时域、频域都具有较高的精度,对于实际系统控制的设计和分析都具有重要意义。
附图说明
图1为本发明一种基于降阶模型的无线电能传输控制系统设计分析方法的流程框图;
图2为局部线性逼近基本原理;
图3为四种谐振回路拓扑的等效降阶电路模型;
其中,(3a)是串联谐振电路的等效降阶电路模型,(3b)是并联谐振电路的等效降阶电路模型,(3c)是LCL型谐振电路的等效降阶电路模型,(3d)是LCC型谐振电路的等效降阶电路模型;
图4为LCL-S型无线电能传输系统实施例;
其中,(4a)是LCL-S型无线电能传输系统拓扑结构图,(4b)是降阶等效电路图;
图5为本发明实例模型与瞬态电路模型、全阶模型对比图;
其中,(5a)是降阶模型与瞬态、全阶模型的串联谐振时域电流波形对比,(5b)是降阶模型与瞬态、全阶模型的LCL谐振电流波形对比;
图6为无线电能传输闭环控制系统等效框图;
图7为系统电流增益的频域特性对比图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。注意,以下的实施方式的说明只是实质上的例示,本发明并不意在对其适用物或其用途进行限定,且本发明并不限定于以下的实施方式。
实施例
如图1所示,一种基于降阶模型的无线电能传输控制系统设计分析方法,该方法包括如下步骤:
1)根据无线电能传输系统拓扑选取交流状态变量降阶量;
2)对交流状态变量降阶量进行变换获取谐振元件的等效近似模型,合并谐振元件,获得降阶等效电路;
3)根据降阶等效电路获得降阶状态空间模型;
4)根据降阶状态空间模型计算闭环系统传递函数,完成控制参数设计和稳定性分析。
步骤2)对交流状态变量降阶量进行变换获取谐振元件的等效近似模型具体为:
21)将交流状态变量降阶量进行傅里叶展开,获取基波傅里叶级数;
22)对交流状态变量降阶量基波傅里叶级数进行局部线性逼近,将带幅值扰动的交流状态变量降阶量波形近似为幅值线性变化的正弦波形;
23)基于近似的正弦波形确定谐振元件的等效近似模型。
步骤1)中若无线电能传输系统拓扑中含有LCL谐振结构或串联谐振结构,则选取谐振电容电压为降阶量。
步骤2)中谐振电容电压近似波形方程为:
步骤3)中谐振电容的等效近似模型确定为等效电感和等效电阻串联形式。
谐振电容的等效近似模型中等效电感大小记作L′,等效电阻大小记作R′:
其中,C为谐振电容的大小,j为虚部。
步骤1)中若无线电能传输系统拓扑中含有并联谐振结构,则选取谐振电感电流为降阶量。
步骤2)中谐振电感电流近似波形方程为:
步骤3)中谐振电感的等效近似模型确定为等效电容和等效电阻串联形式。
谐振电感的等效近似模型中等效电容大小记作C′,等效电阻大小记作R″:
其中,L为谐振电感的大小,j为虚部。
本发明对交流状态变量降阶量基波傅里叶级数进行局部线性逼近过程中,对于各次谐波分别计算,其中基波可近似为谐振频率的正弦波形,幅值波动为低频扰动,如附图2所示,下面以谐振电容电压为例进行说明:
谐振电容的电压基波进行局部线性逼近,在较小的时间步长内,可将带幅值扰动的电压近似为幅值线性变化的正弦波形;
-对上述电压方程求一阶微分,可获得
其中,由于扰动频率远低于谐振频率,在短时间内可视为幅值变化几乎不变,在频域下则有
由此可进一步求得电压方程的二阶微分和其化简形式
结合谐振电流的一次、二次微分方程:
其中,C为谐振电容的容值,iC为流过谐振电容的电流值。
将上述方程联立可求得电容的电压与电流微分形式的关系,此时电容的阻抗可对应等效为阻感元件;
根据所述近似模型将谐振元件进一步化简,可获得降阶等效电路,LC串联谐振、并联谐振和LCL/LCC谐振简化等效电路如附图3所示,有效减少了系统状态变量;
如附图4所示,本实施例提供一种基于降阶模型的无线电能传输系统控制方法,选取原边LCL-副边串联型拓扑结构,此时系统中状态变量为6个(5个交流量,1个直流量),系统动态模型的阶数应为11阶,通过对谐振回路中的状态变量的局部线性近似重构和变量组合,建立了系统的降阶等效电路模型,进一步计算控制系统传递函数,分析系统稳定性。
其建立的具体过程如下:
步骤1对于LCL谐振和串联谐振,分别选取其中的谐振电容电压为降阶量,将状态变量傅里叶展开后,由于系统运行在谐振状态下,基波信号为主导量,此时系统的高频扰动可近似为高次谐波分量,系统的低频扰动等效为幅值波动,
步骤2选取基波进行局部线性近似,将则有:
步骤3对谐振电容电压方程求一阶微分
由于在时间变化较小时,信号幅值变化不明显,则有
由此可进一步求得电压方程的二阶微分和其化简形式
根据电容的阻抗特性,则有电压和电流的关系,对于LCL型和LC串联型分别有:
结合电压方程的微分形式,可分别求出电容的电压电流关系:
即电容可等效为电感和电阻串联形式,为更为直观的表述,可通过等效简化电路实现:
对于LC串联,可将电感和电容合并为阻感元件,如附图3a所示,其等效阻抗为:
此时状态变量仅为谐振电流,降阶模型和瞬态模型、全阶模型的电路波形对比如附图5a所示,模型精度较高;
对于LCL型谐振,系统中存在三个电感连接成T型电路,可通过反向T型变换,等效为两个存在互感的电感,如图3c所示,其等效阻抗为
此时状态变量为电感谐振电流,降阶模型和瞬态模型、全阶模型的电路波形对比如附图5b所示,模型精度较高;
步骤4此时系统的等效电路如附图4b所示,状态变量为4个(3个交流量,1个直流量),系统的差分方程为
由此建立状态空间方程:
步骤5根据建立的降阶状态空间模型,可获得系统增益
G(s)=C(sI-A)-1B
闭环控制系统等效框图如附图6所示,可进一步计算系统传递函数,降阶模型和全阶模型的频域特性如附图7,在谐振频率以下具有较高精度,充分证明了本发明提出方法的有效性和精确性。
上述实施方式仅为例举,不表示对本发明范围的限定。这些实施方式还能以其它各种方式来实施,且能在不脱离本发明技术思想的范围内作各种省略、置换、变更。
Claims (10)
1.一种基于降阶模型的无线电能传输控制系统设计分析方法,其特征在于,该方法包括如下步骤:
1)根据无线电能传输系统拓扑选取交流状态变量降阶量;
2)对交流状态变量降阶量进行变换获取谐振元件的等效近似模型,合并谐振元件,获得降阶等效电路;
3)根据降阶等效电路获得降阶状态空间模型;
4)根据降阶状态空间模型计算闭环系统传递函数,完成控制参数设计和稳定性分析。
2.根据权利要求1所述的一种基于降阶模型的无线电能传输控制系统设计分析方法,其特征在于,步骤2)对交流状态变量降阶量进行变换获取谐振元件的等效近似模型具体为:
21)将交流状态变量降阶量进行傅里叶展开,获取基波傅里叶级数;
22)对交流状态变量降阶量基波傅里叶级数进行局部线性逼近,将带幅值扰动的交流状态变量降阶量波形近似为幅值线性变化的正弦波形;
23)基于近似的正弦波形确定谐振元件的等效近似模型。
3.根据权利要求2所述的一种基于降阶模型的无线电能传输控制系统设计分析方法,其特征在于,步骤1)中若无线电能传输系统拓扑中含有LCL谐振结构或串联谐振结构,则选取谐振电容电压为降阶量。
5.根据权利要求4所述的一种基于降阶模型的无线电能传输控制系统设计分析方法,其特征在于,步骤3)中谐振电容的等效近似模型确定为等效电感和等效电阻串联形式。
7.根据权利要求2所述的一种基于降阶模型的无线电能传输控制系统设计分析方法,其特征在于,步骤1)中若无线电能传输系统拓扑中含有并联谐振结构,则选取谐振电感电流为降阶量。
9.根据权利要求8所述的一种基于降阶模型的无线电能传输控制系统设计分析方法,其特征在于,步骤3)中谐振电感的等效近似模型确定为等效电容和等效电阻串联形式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110131799.8A CN112836369B (zh) | 2021-01-30 | 2021-01-30 | 一种基于降阶模型的无线电能传输控制系统设计分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110131799.8A CN112836369B (zh) | 2021-01-30 | 2021-01-30 | 一种基于降阶模型的无线电能传输控制系统设计分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112836369A true CN112836369A (zh) | 2021-05-25 |
CN112836369B CN112836369B (zh) | 2022-08-19 |
Family
ID=75932496
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110131799.8A Active CN112836369B (zh) | 2021-01-30 | 2021-01-30 | 一种基于降阶模型的无线电能传输控制系统设计分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112836369B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117767288A (zh) * | 2023-12-20 | 2024-03-26 | 安徽农业大学 | 一种基于等效状态变量的vsg并网系统降阶建模方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102467593A (zh) * | 2010-11-09 | 2012-05-23 | 复旦大学 | 非线性电路时域模型降阶方法及装置 |
CN109459934A (zh) * | 2019-01-15 | 2019-03-12 | 华北电力大学 | 一种基于pid控制器整定降阶自抗扰控制器参数的方法 |
CN110266045A (zh) * | 2019-06-30 | 2019-09-20 | 东北电力大学 | 光伏发电系统lcl-t模块的降阶方法 |
CN110262236A (zh) * | 2019-06-20 | 2019-09-20 | 合肥工业大学 | 一种电力电子接口并网系统模型降阶的降阶变量选取方法 |
CN110994666A (zh) * | 2019-10-21 | 2020-04-10 | 合肥工业大学 | 在谐振频率处将阻抗模型降阶为rlc电路的降阶方法 |
-
2021
- 2021-01-30 CN CN202110131799.8A patent/CN112836369B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102467593A (zh) * | 2010-11-09 | 2012-05-23 | 复旦大学 | 非线性电路时域模型降阶方法及装置 |
CN109459934A (zh) * | 2019-01-15 | 2019-03-12 | 华北电力大学 | 一种基于pid控制器整定降阶自抗扰控制器参数的方法 |
CN110262236A (zh) * | 2019-06-20 | 2019-09-20 | 合肥工业大学 | 一种电力电子接口并网系统模型降阶的降阶变量选取方法 |
CN110266045A (zh) * | 2019-06-30 | 2019-09-20 | 东北电力大学 | 光伏发电系统lcl-t模块的降阶方法 |
CN110994666A (zh) * | 2019-10-21 | 2020-04-10 | 合肥工业大学 | 在谐振频率处将阻抗模型降阶为rlc电路的降阶方法 |
Non-Patent Citations (6)
Title |
---|
HAO FENG 等: ""Reduced-Order Modeling and Design of Single-Stage"", 《IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY》 * |
HAO FENG 等: ""Reduced-Order Modeling and Design of Single-Stage"", 《IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY》, 4 March 2020 (2020-03-04) * |
HONGCHANG LI 等: ""Dynamic Phasor-Based Reduced-Order Models"", 《IEEE TRANSAC TIONS ON POWER ELECTRONICS》 * |
HONGCHANG LI 等: ""Dynamic Phasor-Based Reduced-Order Models"", 《IEEE TRANSAC TIONS ON POWER ELECTRONICS》, 11 November 2019 (2019-11-11), pages 11365 * |
曾庆典: ""多路输入高增益DC/DC变换器的研究"", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》 * |
曾庆典: ""多路输入高增益DC/DC变换器的研究"", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》, 15 March 2020 (2020-03-15), pages 46 - 47 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117767288A (zh) * | 2023-12-20 | 2024-03-26 | 安徽农业大学 | 一种基于等效状态变量的vsg并网系统降阶建模方法 |
CN117767288B (zh) * | 2023-12-20 | 2024-06-04 | 安徽农业大学 | 一种基于等效状态变量的vsg并网系统降阶建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112836369B (zh) | 2022-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Input differential-mode EMI of CRM boost PFC converter | |
CN105006839A (zh) | 并网发电系统的弱电网网源荷阻抗模型分析方法 | |
CN110137971B (zh) | 一种三相交流电力弹簧的电压平稳控制方法 | |
CN106357134B (zh) | 一种双向ac-dc-dc单相变换器及其控制方法 | |
CN110768299A (zh) | 负荷虚拟同步机的序阻抗建模与稳定性分析方法 | |
Rong et al. | Output feedback control of single-phase UPQC based on a novel model | |
CN203870168U (zh) | 有源电力滤波器中电抗器电感量的测量装置 | |
CN110138241B (zh) | 一种基于雅可比理论的单相电压源的稳定性分析方法 | |
CN112836369B (zh) | 一种基于降阶模型的无线电能传输控制系统设计分析方法 | |
Kolar et al. | Interference between electric traction supply network and distribution power network-resonance phenomenon | |
Rasekh et al. | Artificial neural network aided loss maps for inductors and transformers | |
CN103795060A (zh) | 有源电力滤波器的输出滤波电路、滤波方法及设计方法 | |
Pantic et al. | Inductive power transfer by means of multiple frequencies in the magnetic link | |
CN104037769B (zh) | 一种单相papf输出lcl滤波器的参数设计方法 | |
CN106130351A (zh) | 一种电动汽车直流充电器输出电压波纹抑制系统及方法 | |
Jayakrishna et al. | Fuzzy logic control based three phase shunt active filter for voltage regulation and harmonic reduction | |
CN103941103B (zh) | 有源电力滤波器中电抗器电感量的测量装置及方法 | |
CN110501919B (zh) | 模块化多电平换流器数字物理混合仿真接口设计方法 | |
CN109492324B (zh) | 基于载波的d类放大器双积分滑模控制器设计方法及电路 | |
Bagchi et al. | Small-signal phasor modeling of an underwater IPT system in constant current distribution | |
CN115236404B (zh) | 一种并网逆变器端口阻抗自测量方法 | |
CN105608285A (zh) | 一种基于电路分析的三相整流设备静态负荷建模方法 | |
CN111259571B (zh) | 一种动车组网侧单相脉冲整流器的稳定性分析方法 | |
CN210578287U (zh) | 一种分数阶d类半波整流器 | |
Almakhles et al. | Single-bit modulator based controller for capacitive power transfer system |
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 |