CN111668826B - 用于纯电动汽车直流供电系统的稳定性分析方法及系统 - Google Patents

用于纯电动汽车直流供电系统的稳定性分析方法及系统 Download PDF

Info

Publication number
CN111668826B
CN111668826B CN202010438023.6A CN202010438023A CN111668826B CN 111668826 B CN111668826 B CN 111668826B CN 202010438023 A CN202010438023 A CN 202010438023A CN 111668826 B CN111668826 B CN 111668826B
Authority
CN
China
Prior art keywords
electric vehicle
pure electric
direct
power supply
supply system
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.)
Active
Application number
CN202010438023.6A
Other languages
English (en)
Other versions
CN111668826A (zh
Inventor
林程
孙建侠
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202010438023.6A priority Critical patent/CN111668826B/zh
Publication of CN111668826A publication Critical patent/CN111668826A/zh
Application granted granted Critical
Publication of CN111668826B publication Critical patent/CN111668826B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • H02J1/00—Circuit arrangements for dc mains or dc distribution networks
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
    • B60L1/00—Supplying electric power to auxiliary equipment of vehicles
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
    • B60L50/00—Electric propulsion with power supplied within the vehicle
    • B60L50/50—Electric propulsion with power supplied within the vehicle using propulsion power supplied by batteries or fuel cells
    • B60L50/51—Electric propulsion with power supplied within the vehicle using propulsion power supplied by batteries or fuel cells characterised by AC-motors
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
    • B60L50/00—Electric propulsion with power supplied within the vehicle
    • B60L50/50—Electric propulsion with power supplied within the vehicle using propulsion power supplied by batteries or fuel cells
    • B60L50/60—Electric propulsion with power supplied within the vehicle using propulsion power supplied by batteries or fuel cells using power supplied by batteries
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B60—VEHICLES IN GENERAL
    • B60R—VEHICLES, VEHICLE FITTINGS, OR VEHICLE PARTS, NOT OTHERWISE PROVIDED FOR
    • B60R16/00—Electric or fluid circuits specially adapted for vehicles and not otherwise provided for; Arrangement of elements of electric or fluid circuits specially adapted for vehicles and not otherwise provided for
    • B60R16/02—Electric or fluid circuits specially adapted for vehicles and not otherwise provided for; Arrangement of elements of electric or fluid circuits specially adapted for vehicles and not otherwise provided for electric constitutive elements
    • B60R16/03—Electric or fluid circuits specially adapted for vehicles and not otherwise provided for; Arrangement of elements of electric or fluid circuits specially adapted for vehicles and not otherwise provided for electric constitutive elements for supply of electrical power to vehicle subsystems or for
    • B60R16/033—Electric or fluid circuits specially adapted for vehicles and not otherwise provided for; Arrangement of elements of electric or fluid circuits specially adapted for vehicles and not otherwise provided for electric constitutive elements for supply of electrical power to vehicle subsystems or for characterised by the use of electrical cells or batteries
    • 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
    • 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
    • H02J7/00—Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries
    • H02J7/0063—Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries with circuits adapted for supplying loads from the battery
    • 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
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • 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
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/60—Other road transportation technologies with climate change mitigation effect
    • Y02T10/70—Energy storage systems for electromobility, e.g. batteries

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Mechanical Engineering (AREA)
  • Transportation (AREA)
  • Physics & Mathematics (AREA)
  • Sustainable Energy (AREA)
  • Sustainable Development (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Electric Propulsion And Braking For Vehicles (AREA)

Abstract

本发明提供用于纯电动汽车直流供电系统的稳定性分析方法及系统，所述方法包括步骤建立纯电动汽车直流供电系统的数学模型及直流母线电压传递函数；基于Nyquist定理，提出根据闭环传递函数相角条件判定系统稳定性的方法；根据所述判定系统稳定性的方法，对所述纯电动汽车直流母线电压传递函数进行分析，确定其稳定条件。本发明所得的稳定性条件对系统中的线路阻抗、滤波器阻抗等参数做出了限制，为纯电动汽车直流供电系统的设计提供了依据。

Description

用于纯电动汽车直流供电系统的稳定性分析方法及系统
技术领域
本发明涉及稳定性分析领域，具体涉及用于纯电动汽车直流供电系统的稳定性分析方法及系统。
背景技术
随着汽车性能的提升，用电形式的不断变化，系统越来越复杂，系统稳定性问题上也面临着巨大的挑战。纯电动汽车(EV)的一些负载，如主驱电机、低压负载等，其前端的变换器采用了严格的闭环控制，可以视为恒功率负载(CPL)；其他负载，如加热器等，可以视为恒电压负载(CVL)。EV的电气系统中多个变换器之间级联或并联，相互耦合影响严重。线路上分布参数阻抗的影响，各变换器滤波回路的相互关联，使得整个EV直流供电系统形成了一个复杂的高阶系统。并且，CPL存在较强的负阻抗问题，降低了系统的稳定性，严重时可能会导致整个直流供电系统发生谐振或失稳。系统谐振或失稳时将直接导致系统停止工作，并损坏电路器件，造成严重的影响。
目前，对直流系统的稳定性分析多集中于单电源单负载的级联系统稳定性分析，缺乏对于多负载并联系统的稳定性分析，导致对纯电动汽车直流供电系统的稳定性分析不足，因此，本发明提出一种用于纯电动汽车直流供电系统的稳定性分析方法及系统，用于对纯电动汽车直流供电系统的稳定性进行分析。
发明内容
为解决以上技术问题，本发明提供了一种用于纯电动汽车直流供电系统的稳定性分析方法及系统，通过建立系统的闭环传递函数模型，提出一种根据闭环传递函数相角条件判定系统稳定性的方法，并将这种稳定性判定方法应用于纯电动汽车直流供电系统，得到了系统在多负载并联时的稳定条件。所得的稳定性条件对系统中的线路阻抗、滤波器阻抗等参数做出了限制，为纯电动汽车直流供电系统的设计提供了依据。
本发明的第一方面提供一种用于纯电动汽车直流供电系统的稳定性分析方法，包括步骤：
建立纯电动汽车直流供电系统的数学模型及纯电动汽车直流母线电压传递函数；
基于Nyquist定理，提出根据闭环传递函数相角条件判定系统稳定性的方法；
根据所述判定系统稳定性的方法，对所述纯电动汽车直流母线电压传递函数进行分析，确定其稳定条件。
优选的是，所述纯电动汽车直流供电系统的数学模型为：将电源侧阻抗Zs等效为电阻Rs及电感Ls的组合，将负载前端滤波器等效为电阻RLi、电感Li及电容Ci的组合，将CPL等效为负电阻-Ri，将CVL等效为电阻Rr，纯电动汽车直流母线电压的传递函数Gc(s)如下式所示：
Figure BDA0002503024760000021
其中，Udc为母线电压，Us为电源电压，Y∑为母线上并联的总负载导纳，Fi(s)如下式所示：
Figure BDA0002503024760000022
且Mi(s)＝1/n+Fi(s)，即
Figure BDA0002503024760000023
其中
Figure BDA0002503024760000024
且Rr(s)表示为
Figure BDA0002503024760000031
在上述任一方案中优选的是，基于Nyquist定理，提出根据闭环传递函数相角条件判定系统稳定性的方法包括步骤：
建立系统的复数域模型，用系统闭环传递函数G1c(s)定义系统的输入输出关系；
设定所述系统闭环传递函数G1c(s)满足前提条件：A、所有的零点均位于s平面的左半平面，B、在s平面的虚轴上没有极点，C、传递函数的分母阶数比分子阶数大一阶；
根据Nyquist定理进而确定系统稳定条件D、
Figure BDA0002503024760000032
其中
Figure BDA0002503024760000033
表示传递函数G1c(s)的相角，b0表示传递函数G1c(s)的系数。
在上述任一方案中优选的是，根据所述判定系统稳定性的方法，对所述纯电动汽车直流母线电压传递函数进行分析，确定其稳定条件包括步骤：
根据纯电动汽车直流母线电压的传递函数Gc(s)需满足条件A、所有的零点均位于s平面的左半平面，确定纯电动汽车直流供电系统的稳定条件1；
根据纯电动汽车直流母线电压的传递函数Gc(s)需满足条件B、在s平面的虚轴上没有极点，确定纯电动汽车直流供电系统的稳定条件2；
根据纯电动汽车直流母线电压的传递函数Gc(s)需满足条件C、传递函数的分母阶数比分子阶数大一阶，确定纯电动汽车直流供电系统的稳定条件3；
进一步地，根据纯电动汽车直流母线电压的传递函数Gc(s)需满足条件D、
Figure BDA0002503024760000034
确定纯电动汽车直流供电系统的稳定条件4；
综合所述稳定条件1、稳定条件2、稳定条件3和稳定条件4，确定纯电动汽车直流供电系统的综合稳定条件，在满足所述综合稳定条件时，系统一定是稳定的。
根据所述纯电动汽车直流供电系统的综合稳定条件对所述纯电动汽车直流供电系统的稳定性进行分析。
在上述任一方案中优选的是，根据纯电动汽车直流供电系统的传递函数Gc(s)需满足条件A、所有的零点均位于s平面的左半平面，确定纯电动汽车直流供电系统的稳定条件1为：
Figure BDA0002503024760000041
在上述任一方案中优选的是，根据纯电动汽车直流供电系统的传递函数Gc(s)需满足条件B、在s平面的虚轴上没有极点，确定纯电动汽车直流供电系统的稳定条件2为：
Figure BDA0002503024760000042
在上述任一方案中优选的是，根据纯电动汽车直流供电系统的传递函数Gc(s)需满足条件C、传递函数的分母阶数比分子阶数大一阶，确定纯电动汽车直流供电系统的稳定条件3为：无限制。
在上述任一方案中优选的是，根据纯电动汽车直流供电系统的传递函数Gc(s)需满足条件D、
Figure BDA0002503024760000043
确定纯电动汽车直流供电系统的稳定条件4为：
Figure BDA0002503024760000044
在上述任一方案中优选的是，综合所述稳定条件1、稳定条件2、稳定条件3和稳定条件4，确定纯电动汽车直流供电系统的综合稳定条件为：
Figure BDA0002503024760000045
本发明的第二方面提供一种用于纯电动汽车直流供电系统的稳定性分析系统，其用于运行所述的用于纯电动汽车直流供电系统的稳定性分析方法，对所述纯电动汽车直流供电系统的稳定性进行分析。
本发明通过建立系统的闭环传递函数模型，提出一种根据闭环传递函数相角条件判定系统稳定性的方法，并将这种稳定性判定方法应用于纯电动汽车直流供电系统，得到了纯电动汽车直流供电系统在多负载并联时的稳定条件。所得的稳定性条件对系统中的线路阻抗、滤波器阻抗等参数做出了限制，为纯电动汽车直流供电系统的设计提供了依据，本发明提出的稳定性分析方法同样适用于其他具有分布式结构的直流系统。
附图说明
图1为按照本发明的用于纯电动汽车直流供电系统的稳定性分析方法的流程示意图。
图2为按照本发明的用于纯电动汽车直流供电系统的稳定性分析方法的纯电动汽车直流供电系统示意图。
图3为按照本发明的用于纯电动汽车直流供电系统的稳定性分析方法的纯电动汽车直流供电系统等效电路图。
图4为按照本发明的用于纯电动汽车直流供电系统的稳定性分析方法的修正后的Nyquist轨迹及对应的G1c(s)曲线图。
图5为按照本发明的用于纯电动汽车直流供电系统的稳定性分析方法的系统稳定与失稳时的G1c(s)曲线示意图。
图6为按照本发明的用于纯电动汽车直流供电系统的稳定性分析方法的N(s)、D(s)的相频示意图。
具体实施方式
为了更好地理解本发明，下面结合具体实施例对本发明作详细说明。
实施例1
如图1所示，一种用于纯电动汽车直流供电系统的稳定性分析方法，包括步骤：
S1、建立纯电动汽车直流供电系统的数学模型及纯电动汽车直流母线电压传递函数；
S2、基于Nyquist定理，提出根据闭环传递函数相角条件判定系统稳定性的方法；
S3、根据所述判定系统稳定性的方法，对所述纯电动汽车直流母线电压传递函数进行分析，确定其稳定条件。
图2所示为纯电动汽车直流供电系统示意图。纯电动汽车(EV)的直流供电系统是一个复杂的高阶系统，存在多个变换器之间级联或并联的情况，各变换器的滤波回路相互的关联，相互耦合影响严重。EV的一些典型负载，如主驱电机、DC/DC负载等，其前端的变换器采用了严格的闭环控制，可以视为恒功率负载(CPL)；其他负载，如加热器等阻性负载，可以视为恒电压负载(CVL)。由于实际的纯电动汽车直流供电系统是一个复杂的非线性系统，因此为了进行稳定性分析，需要对系统进行简化，图3为纯电动汽车直流供电系统的等效电路图。考虑母线并联的负载为CPL和CVL，在步骤S1中，建立的所述纯电动汽车直流供电系统的数学模型为：将电源侧阻抗Zs等效为电阻Rs及电感Ls的组合，将负载前端滤波器等效为电阻RLi、电感Li及电容Ci的组合，将CPL等效为负电阻-Ri，将CVL等效为电阻Rr，纯电动汽车直流母线电压的传递函数Gc(s)如下式所示：
Figure BDA0002503024760000061
其中，Udc为母线电压，Us为电源电压，Y∑为母线上并联的总负载导纳，Fi(s)如下式所示：
Figure BDA0002503024760000062
且Mi(s)＝1/n+Fi(s)，即
Figure BDA0002503024760000071
其中
Figure BDA0002503024760000072
且Rr(s)表示为
Figure BDA0002503024760000073
n为母线上并联的CPL的个数，i＝1，2，…，n。
步骤S2中，基于Nyquist定理，根据闭环传递函数相角条件判定系统稳定性的方法包括步骤：
S21、建立系统的复数域模型，用系统闭环传递函数G1c(s)定义系统的输入输出关系，系统闭环传递函数G1c(s)的一般表达式为：
Figure BDA0002503024760000074
其中z1，z2，…，zm表示G1c(s)的m个零点，λ1，λ2，…，λn表示G1c(s)的n个极点。
S22、设定所述系统闭环传递函数G1c(s)满足前提条件：A、所有的零点均位于s平面的左半平面，B、在s平面的虚轴上没有极点，C、传递函数的分母阶数比分子阶数大一阶。设定条件C是因为：纯电动汽车的直流供电系统在并联CPL及CVL负载时，系统传递函数的分母比分子高一阶。
S23、根据Nyquist定理，进而确定系统稳定条件D、
Figure BDA0002503024760000075
其中
Figure BDA0002503024760000076
表示传递函数G1c(s)的相角，b0表示传递函数G1c(s)的系数。
步骤S23中，(1)将s平面的Nyquist围线分为三部分，分别为低频部分：s＝0，对应G1c(s)的起点K，由于s＝0不是G1c(s)的零点或极点，因此K为一个非0且有限的实数，所以G1c(s)曲线的起点K位于实轴。中频部分：s＝jω，由于虚轴上没有零点或极点，G1c(s)的幅值不为0且有限。高频部分：s＝∞，对应G1c(s)的终点M，由于G1c(s)的分母阶数比分子阶数大，因此M为0，终点M位于原点。(2)应用Nyquist围线映射定理的前提是Nyquist轨迹不包含任何零极点，为应用Nyquist围线映射定理，对s平面的Nyquist围线轨迹进行修正，用一个半径ρ趋近于∞的半圆代替s平面上的无限大半径圆，经过修正的Nyquist围线轨迹半径足够大且趋于∞，且不包含零点和极点。图4为修正后的Nyquist轨迹及对应的G1c(s)曲线图，以G1c(s)具有正系数b0为例，此时经过修正后的Nyquist围线映射到G1c(s)平面为一个半径趋近于0的半圆，且位于G1c(s)平面的右半平面。s平面的点A,B,C,D对应G1c(s)平面的点A*,B*,C*及D*。经过修正的轨迹半径足够大且趋于∞，且不包含零极点，可以应用围线映射定理。(3)对于半径ρ的半圆轨迹，复变量s可以表示为s＝ρejθ，θ从90°变化到-90°。因n－m＝1，且假设G1c(s)的系数b0为正时，则有
Figure BDA0002503024760000081
用s＝ρejθ代替上式中的s＝∞，则G1c(s)可以表示为Glc(ρejθ)＝e-jθ/ρ。此时，G1c(s)的幅值为正，且随着ρ趋于∞，G1c(s)的幅值趋于0，且G1c(s)的相角从-90°变化到90°。因此，s右平面上趋于无限大半径圆轨迹映射到G1c(s)右平面一个半径趋近于0的圆。假设G1c(s)的系数b0为负时，则有Glc(ρejθ)＝-e-jθ/ρ。此时，G1c(s)的幅值为负，且随着ρ趋于∞，G1c(s)的幅值趋于0，且G1c(s)的相角从-90°变化到90°。此时，s右平面上趋于无限大半径圆轨迹映射到G1c(s)左平面一个半径趋近于0的圆。(4)根据Nyquist围线映射定理，G1c(s)曲线围绕原点的次数等于Nyquist围线轨迹包含的零极点个数差，根据系统稳定性要求，Nyquist曲线不应包含极点，又因为s右半面不包含零点，因此，G1c(s)曲线应围绕原点的次数应为0，因此G1c(s)曲线与实轴的所有交点应均位于正实轴或负实轴，否则，如果G1c(s)曲线与实轴的交点位于不同的实半轴，那么G1c(s)曲线一定会包围原点，系统一定失稳。图5所示为系统稳定与失稳时的G1c(s)曲线示意图，图5(a)所示为G1c(s)曲线与实轴的交点均位于同样的实半轴，系统稳定，图5(b)所述为G1c(s)曲线与实轴的交点均位于不同的实半轴，包围原点，系统失稳。(5)当传递函数的表达式确定后，G1c(s)的系数b0的正负也可以确定。若b0为正，则G1c(s)的相角
Figure BDA0002503024760000093
应始终不等于±180°以保证G1c(s)曲线与实轴的交点不位于负实轴。若b0为负，则G1c(s)的相角
Figure BDA0002503024760000094
应始终不等于0°以保证G1c(s)曲线与实轴的交点不位于正实轴。因此，对于满足条件A、B、C的闭环传递函数G1c(s)，其稳定条件D可以表述为：
Figure BDA0002503024760000091
同时传递函数的表达式确定后，G1c(s)曲线起点K的位置也可以确定。若G1c(s)的系数b0为正，曲线起点K的位置应位于实轴的正半轴；若G1c(s)的系数b0为负，曲线起点K的位置应位于实轴的负半轴。否则，G1c(s)曲线将包围原点，导致系统失稳。
步骤3中，根据所述判定系统稳定性的方法，对所述纯电动汽车直流母线电压传递函数进行分析，确定其稳定条件包括步骤：
S31、根据纯电动汽车直流母线电压的传递函数Gc(s)需满足条件A所有的零点均位于s平面的左半平面，确定纯电动汽车直流供电系统的稳定条件1；
S32、根据纯电动汽车直流母线电压的传递函数Gc(s)需满足条件B在s平面的虚轴上没有极点，确定纯电动汽车直流供电系统的稳定条件2；
S33、根据纯电动汽车直流母线电压的传递函数Gc(s)需满足条件C传递函数的分母阶数比分子阶数大一阶，确定纯电动汽车直流供电系统的稳定条件3；
S34、根据纯电动汽车直流母线电压的传递函数Gc(s)需满足相角条件
Figure BDA0002503024760000092
确定纯电动汽车直流供电系统的稳定条件4；
S35、综合所述稳定条件1、稳定条件2、稳定条件3和稳定条件4，确定纯电动汽车直流供电系统的综合稳定条件，在满足所述综合稳定条件时，系统一定是稳定的。
S36、根据所述纯电动汽车直流供电系统的综合稳定条件对所述纯电动汽车直流供电系统的稳定性进行分析。
步骤S31中，根据纯电动汽车直流母线电压的传递函数Gc(s)需满足条件A、所有的零点均位于s平面的左半平面，确定纯电动汽车直流供电系统的稳定条件1为：
Figure BDA0002503024760000101
根据Gc(s)的表达式可以确定，Gc(s)的零点即为Fi(s)的极点，Gc(s)的零点位于s平面的左半平面，则Fi(s)的极点位于s平面的左半平面，因此需要
Figure BDA0002503024760000102
步骤S32中，根据纯电动汽车直流母线电压的传递函数Gc(s)需满足条件B、在s平面的虚轴上没有极点，确定纯电动汽车直流供电系统的稳定条件2为：
Figure BDA0002503024760000103
在
Figure BDA0002503024760000104
的条件限制下，Mi(s)的分母多项式系数均为正。如果可以使Mi(s)的分子多项式系数均为正，则∑Mi(s)的分母、分子多项式系数均为正且不缺项，这是因为所有项的系数均为正时，所求的和也分别为正，不存在抵消的情况。同理，由于Rr(s)自身具有正系数，因此∑Mi(s)+Rr(s)的分母、分子多项式系数均为正，此时，Gc(s)可以表示为
Figure BDA0002503024760000105
其中，α0…αm及β0…βm+1均大于0。由于Gc(s)的分母多项式系数均为正且不缺项，因此，Gc(s)的分母不包含s、Ts2+1因子，Gc(s)也就不存在虚轴上的极点，因此，要求Mi(s)的分子多项式系数a、b、c均为正，在满足
Figure BDA0002503024760000111
时，a、b、c一定为正，此时，Gc(s)不存在虚轴上的极点。
步骤S33中，根据纯电动汽车直流母线电压的传递函数Gc(s)需满足条件C、传递函数的分母阶数比分子阶数大一阶，确定纯电动汽车直流供电系统的稳定条件3为：无限制。根据Gc(s)的表达式，Gc(s)曲线的高频部分为
Figure BDA0002503024760000112
可得Gc(s)的分母阶数比分子阶数大一阶，因此稳定条件3为无限制。
步骤S34中，根据纯电动汽车直流母线电压的传递函数Gc(s)需满足条件D、
Figure BDA0002503024760000113
确定纯电动汽车直流供电系统的稳定条件4为：
Figure BDA0002503024760000114
根据Gc(s)的表达式，Gc(s)的系数b0为正。为满足系统稳定性要求，Gc(s)的相角
Figure BDA0002503024760000115
应始终不等于±180°。Gc(s)的分母多项式由Rr(s)和∑Mi(s)两部分组成，其中，Rr(s)的相角范围为[0°，90°]。接下来分析∑Mi(s)的相角范围。为简化分析，将Mi(s)表示为
Figure BDA0002503024760000116
在低频段，Mi(s)＝c/E，相角为0°；在高频段，Mi(s)＝a/A，相角为0°。Mi(s)的分子N(s)、分母D(s)均为二阶环节，相角范围在[0°，180°]。N(s)、D(s)的转折频率分别为
Figure BDA0002503024760000117
可知，必有ωc1<ωc2，因此，在ωc1处有∠N(s)>∠D(s)。如图6所示为N(s)、D(s)的相频示意图，在中频段内，若可以使N(s)、D(s)的相角始终不相等，则必有∠N(s)>∠D(s)，使得Mi(s)的相角范围为(0°，180°)。进而在整个频域内，Mi(s)的相角范围为[0°，180°)。根据N(s)、D(s)的表达式，在其相角相等时有
Figure BDA0002503024760000118
可以求解得
Figure BDA0002503024760000121
ω无解时N(s)、D(s)的相角始终不相等。由上式可知，此时一定有bE>Bc，因此，在Ab<Ba时，ω无解。代入公式可以求得
Figure BDA0002503024760000122
化简后可得
Figure BDA0002503024760000123
此时，可以使得Mi(s)在整个频域内的相角范围为[0°，180°)。根据矢量和的原理，∑Mi(s)+Rr(s)的相角范围为[0°，180°)，因此，Gc(s)的相角范围为(-180°，0°]，满足不等于±180°的要求，因此稳定条件4为
Figure BDA0002503024760000124
步骤S35中，综合所述稳定条件1、稳定条件2、稳定条件3和稳定条件4，确定纯电动汽车直流供电系统的综合稳定条件为：
Figure BDA0002503024760000125
步骤S36中，根据所述纯电动汽车直流供电系统的综合稳定条件对所述纯电动汽车直流供电系统的稳定性进行分析。稳定条件1保证了各个CPL前端电容电压的稳定性，在此基础上，稳定性条件2、3、4保证了直流母线电压的稳定性，因此，根据所述综合稳定条件，纯电动汽车直流供电系统各项取值在满足相应的条件时可以保证系统全局稳定性，为实际的纯电动汽车直流供电系统设计提供依据。
实施例2
一种适用于纯电动汽车直流供电系统稳定性分析的系统，其用于运行所述的适用于纯电动汽车直流供电系统稳定性分析的方法，对所述纯电动汽车直流供电系统的稳定性进行分析。应当可以理解，实现所述适用于纯电动汽车直流供电系统稳定性分析的方法的全部或部分步骤，可以通过计算机程序来指示相关的硬件完成，所述的计算机程序可存储于非易失性的计算机可读存储介质中，该计算机程序在执行时，可实施如上述方法的全部或者部分步骤。如在所述非易失性的计算机可读存储介质中存储有所述综合判断条件，当输入某纯电动汽车直流供电系统的结构参数后，系统自动对所述纯电动汽车直流供电系统进行数学模型建立，根据所述综合稳定条件，判断某一个或几个部件的取值范围以使所述纯电动汽车直流供电系统保持稳定。
需要说明的是，以上实施例仅用于说明本发明的技术方案，而非对其限制；尽管前述实施例对本发明进行了详细的说明，本领域的技术人员应该理解：其可以对前述实施例记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换，而这些替换，并不使相应技术方案的本质脱离本发明技术方案的范围。

Claims (2)

1.用于纯电动汽车直流供电系统的稳定性分析方法，其特征在于：包括步骤：
建立纯电动汽车直流供电系统的数学模型及纯电动汽车直流母线电压传递函数；
基于Nyquist定理，提出根据闭环传递函数相角条件判定系统稳定性的方法；
根据所述判定系统稳定性的方法，对所述纯电动汽车直流母线电压传递函数进行分析，确定其稳定条件；
所述纯电动汽车直流供电系统的数学模型为：将电源侧阻抗Zs等效为电阻Rs及电感Ls的组合，将负载前端滤波器等效为电阻RLi、电感Li及电容Ci的组合，将恒功率负载等效为负电阻-Ri，将恒电压负载等效为电阻Rr，纯电动汽车直流母线电压的传递函数Gc(s)如下式所示：
Figure FDA0003347227410000011
其中，Udc为母线电压，Us为电源电压，Y∑为母线上并联的总负载导纳，Fi(s)如下式所示：
Figure FDA0003347227410000012
且Mi(s)＝1/n+Fi(s)，即
Figure FDA0003347227410000013
其中
Figure FDA0003347227410000014
且Rr(s)表示为
Figure FDA0003347227410000021
所述n为母线上并联的CPL的个数；
基于Nyquist定理，提出根据闭环传递函数相角条件判定系统稳定性的方法包括步骤：
建立系统的复数域模型，用系统闭环传递函数G1c(s)定义系统的输入输出关系，其中
Figure FDA0003347227410000022
其中z1，z2，…，zm表示G1c(s)的m个零点，λ1，λ2，…，λn表示G1c(s)的n个极点；
设定所述系统闭环传递函数G1c(s)满足前提条件：A、所有的零点均位于s平面的左半平面，B、在s平面的虚轴上没有极点，C、传递函数的分母阶数比分子阶数大一阶；
根据Nyquist定理进而确定系统稳定条件D、
Figure FDA0003347227410000023
其中
Figure FDA0003347227410000024
表示G1c(s)的相角，b0表示G1c(s)的系数；
根据所述判定系统稳定性的方法，对所述纯电动汽车直流母线电压传递函数进行分析，确定其稳定条件包括步骤：
根据纯电动汽车直流供电系统的传递函数Gc(s)需满足条件A所有的零点均位于s平面的左半平面，确定纯电动汽车直流供电系统的稳定条件1，所述稳定条件1为：
Figure FDA0003347227410000025
根据纯电动汽车直流供电系统的传递函数Gc(s)需满足条件B在s平面的虚轴上没有极点，确定纯电动汽车直流供电系统的稳定条件2，所述稳定条件2为：
Figure FDA0003347227410000026
根据纯电动汽车直流供电系统的传递函数Gc(s)需满足条件C传递函数的分母阶数比分子阶数大一阶，确定纯电动汽车直流供电系统的稳定条件3，所述稳定条件3为：无限制；
根据纯电动汽车直流供电系统的传递函数Gc(s)需满足条件
Figure FDA0003347227410000031
确定纯电动汽车直流供电系统的稳定条件4，所述稳定条件4为：
Figure FDA0003347227410000032
综合所述稳定条件1、稳定条件2、稳定条件3和稳定条件4，确定纯电动汽车直流供电系统的综合稳定条件，所述综合稳定条件为：
Figure FDA0003347227410000033
根据所述纯电动汽车直流供电系统的综合稳定条件对所述纯电动汽车直流供电系统的稳定性进行分析。
2.用于纯电动汽车直流供电系统的稳定性分析系统，其特征在于：其用于运行如权利要求1所述的用于纯电动汽车直流供电系统的稳定性分析方法，对所述纯电动汽车直流供电系统的稳定性进行分析。
CN202010438023.6A 2020-05-21 2020-05-21 用于纯电动汽车直流供电系统的稳定性分析方法及系统 Active CN111668826B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010438023.6A CN111668826B (zh) 2020-05-21 2020-05-21 用于纯电动汽车直流供电系统的稳定性分析方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010438023.6A CN111668826B (zh) 2020-05-21 2020-05-21 用于纯电动汽车直流供电系统的稳定性分析方法及系统

Publications (2)

Publication Number Publication Date
CN111668826A CN111668826A (zh) 2020-09-15
CN111668826B true CN111668826B (zh) 2022-03-15

Family

ID=72384321

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010438023.6A Active CN111668826B (zh) 2020-05-21 2020-05-21 用于纯电动汽车直流供电系统的稳定性分析方法及系统

Country Status (1)

Country Link
CN (1) CN111668826B (zh)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10090764B1 (en) * 2018-01-10 2018-10-02 National Technology & Engineering Solutions Of Sandia, Llc Method to provide meta-stable operation of a DC microgrid comprising a pulsed load
CN109459613A (zh) * 2018-09-19 2019-03-12 中国南方电网有限责任公司 一种换流变压器铁心饱和型谐波不稳定分析方法
CN110011295B (zh) * 2019-03-06 2023-11-03 中国电力科学研究院有限公司 一种分布式直流电网的稳定性分析方法及系统
CN110729729B (zh) * 2019-11-14 2021-01-05 山东大学 弱电网连接下电压源型换流器电压前馈环节作用分析方法

Also Published As

Publication number Publication date
CN111668826A (zh) 2020-09-15

Similar Documents

Publication Publication Date Title
CN108649780B (zh) 一种考虑弱电网下逆变器稳定性的lcl滤波器参数优化方法
US11901830B2 (en) Method for controlling the input voltage frequency of a DC-DC converter
CN107294225B (zh) 副边并联补偿的ipt系统负载与互感参数识别方法
CN112054694B (zh) 基于最小电流应力的双向变换器优化控制方法及装置
CN108322050A (zh) 一种适用于谐振网络的拓扑优化与元件参数优化方法
CN106169879A (zh) 修正注入零序分量的vienna整流器调制方法、控制器及系统
CN111668826B (zh) 用于纯电动汽车直流供电系统的稳定性分析方法及系统
CN110266044B (zh) 一种基于储能变流器的微电网并网控制系统及方法
CN107623341B (zh) 向无源网络供电的vsc逆变站数学模型及内模控制器
CN116799928B (zh) 一种s-s型无线供电系统的电容参数补偿方法
WO2023236624A1 (zh) 一种并联apf的控制方法及装置
CN115622377B (zh) 一种基于分数阶元件的阻抗匹配方法及系统
CN109165474B (zh) 一种基于神经网络的逆变器重复控制设计方法
CN111654016B (zh) 一种适用于分布式直流电网稳定性分析的方法及系统
CN110134004A (zh) 一种基于电力弹簧电路结构的pi控制参数整定方法
CN115940675A (zh) 一种整流器无模型预测电流控制方法及系统
CN109831114A (zh) 一种llc类双向有源桥逆变器谐振参数设计方法
CN116388214A (zh) 构网型变流器控制参数整定方法及装置
CN110165677B (zh) 一种单相交流电力弹簧的电压平稳控制方法
CN114597975A (zh) 并联系统的功率分配方法及控制设备
CN112054696A (zh) 基于最小回流功率的多电平变换器优化控制方法及装置
CN112332689A (zh) 基于几何变换的中点箝位型三电平变换器的调制方法
CN112421605A (zh) 基于无源积分的直流微电网改进下垂控制方法
CN112615378A (zh) 配电网高频谐振频移方法、装置以及计算机可读存储介质
Stanisławski et al. The steady-state error issue in Laguerre-based modeling of the GL fractional difference

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