320 likes | 481 Vues
第 3 ç« ç”µåŠ›ç”µå电路与系统的建模与仿真. 3.1 概述. 电路与系统建模与仿真的概念 电路的计算机仿真需è¦è§£å†³çš„基本问题是: ①建立电路方程和仿真模型; ②求解电路方程的算法。 除æ¤ä¹‹å¤–,å¯è§†åŒ–电路录入ã€ä»¿çœŸç»“果的分æžä¸Žå¤„ç†ä»¥åŠæ³¢å½¢åˆ†æžç‰é—®é¢˜ä¹Ÿæ˜¯è®¡ç®—机仿真必须很好解决的关键问题。. 常常用的电路系统仿真方法有如下五ç§ï¼š å°ä¿¡å·åˆ†æž 离散时域仿真方法 ç‰æ•ˆç”µè·¯æ³• Laplace å˜æ¢æ³• 周期时间åºåˆ—分æž. 1. å°ä¿¡å·åˆ†æž.
E N D
第3章 电力电子电路与系统的建模与仿真 3.1 概述 电路与系统建模与仿真的概念 电路的计算机仿真需要解决的基本问题是: ①建立电路方程和仿真模型; ②求解电路方程的算法。 除此之外,可视化电路录入、仿真结果的分析与处理以及波形分析等问题也是计算机仿真必须很好解决的关键问题。
常常用的电路系统仿真方法有如下五种: 小信号分析 离散时域仿真方法 等效电路法 Laplace变换法 周期时间序列分析
1. 小信号分析 小信号分析是要建立某一个工作点附近的近似小信号线性模型,该工作点由变换器参数、输入电压和负载决定。这样,就可以采用线性系统分析方法对开关变换器的调节控制回路进行设计。但是当输入电压或者负载发生变化时,由于工作点将发生变化,这时该分析方法将受到限制。当变换器的自然频率远小于开关频率的一半时,该方法带来的误差是可以接受的。 开关电源的小信号分析是分析电源变换器动态性能的有力工具,也是系统动态设计的依据。其显著优点在于物理概念清楚,可用伯德(Bode)图设计校正环节,因此该方法在电力电子系统的分析、仿真与设计中得到普遍重视。
2. 离散时域仿真方法 电力电子系统是一种强非线性动态系统,要准确地找到其解析解是相当困难的。为此,1979年美国弗吉尼亚电力电子技术中心(VPEC)的李泽元教授首先提出了开关DC-DC变换器的离散时域仿真方法。80年代后期,清华大学蔡宣三教授对该方法进行了深入的研究。离散时域仿真方法是研究拓扑变化及元件参数变化对系统瞬态特性影响的有力工具。离散时域仿真方法的基本思路是:利用状态空间法列出非线性系统的分段线性方程,找出状态转移规律,并得出非线性差分方程,用计算机进行求解。运用该方法可以较精确地对开关电源等电力电子系统进行数字分析。
3. 等效电路法 在开关变换器中,开关元件的作用是使某一支路以一定的占空比接通或断开,所以这些元件的电压和电流平均值常常与电路中另外某条支路的电流或电压的平均值有关。因此,这些元件可近似用一个与占空比有关的受控源来代替。 等效电路法就是应用一个载波周期内平均值的概念,把开关变换器变为一个含有受控源的线性电路,然后用求解线性电路的方法对开关变换器进行稳态和小信号分析。 例如图所示的三端开关器件电路可用图中所示的受控源模型代替。 三端开关电路及其受控源模型
4. Laplace变换法 Laplace变换法首先写出网络的频域方程式,得出所求电压或电流的频域表达式,再将其进行Laplace反变换,从而得出网络解的时域表达式。 使用Laplace变换法对开关变换器进行仿真时,要写出系统的频域方程式是非常困难的。考虑到开关器件的特殊性,通常用时域频域混合方程式来描述电源变换器。 处理开关器件时,引入开关函数 如果以上图所示的理想开关元件模拟上页图中的开关电路,由图可以写出如下的表达式:
那么图中的三端开关器件电路可以写为如下的时域表达式那么图中的三端开关器件电路可以写为如下的时域表达式 式中[x, y]定义为: 其中阶段I指晶体管T导通,二极管D不导通,阶段II指晶体管T不导通,二极管D导通,阶段III指晶体管T和二极管D都不导通。
例1-1 如图所示为BUCK变换器,试用Laplace变换法写出网络的频域方程式,以便得出网络解的时域表达式。 BUCK变换器
解:首先写出时域额域混合方程式为 在这五个方程式中,第一个和第二个方程是考虑晶体管T为理想开关时的时域方程,第三个和第五个方程是频域方程,第四个方程的时域和频域方程形式是相似的。为了区别,方程中的频域量包含s。
求解时可以总体求解,也可以分阶段求解。分阶段求解可得出三个阶段的解分别是:求解时可以总体求解,也可以分阶段求解。分阶段求解可得出三个阶段的解分别是: 阶段I
5. 周期时间序列分析 周期时间序列分析仿真方法是一个稳态仿真方法。由该方法可以获得变换器在一个开关周期的稳态响应波形,由此可以对其进行稳态特性研究及谐波分析。 开关变换器由于开关的接通和断开使得在不同的时刻变换器中各处的电流电压关系完全不同,不能够用一个解析表达式来描述。但是,在每一个时间点处,电流电压的关系却有着确定的关系,因此,如果我们将电流和电压量用一个周期时间序列来表示,则可以得出一个固定的矩阵方程式来描述变换器的特性。 周期时间序列分析仿真方法将随时间变化的有源电路部分单独进行处理,而非时变线性无源网络则单独处理,很容易得出它的矩阵方程式,其方法如下:
对于一个M端非时变线性无源网络,正弦稳态方程式的一般形式如下:对于一个M端非时变线性无源网络,正弦稳态方程式的一般形式如下: (1) 式中i,m为节点编号;k为谐波次数;Uik为第i个节点电位的第k次谐波电压相量,k=0时为直流分量;Zimk为第i个节点与第m个节点之间的阻抗;Imk为流入第m个节点的电流。 同时将Ui表示为傅氏级数,即 (2) (3) 其中:
则可以得到: (4) 可以看出,上式包括直流分量和与开关时变有源电路接口处电流作用分量,同时完成了频域到时域的转换。 在进行计算机仿真分析时,计算机能够处理的都是离散量,因此,还必须对上式进行离散化处理,以满足计算机仿真的需要。为此,在每一个基波周期中对端口电流进行P次采样,以 ( )表示采样值,则上式中的直流分量时域离散值可表示为:
式(4)中的第二项表示与开关时变有源电路接口处电流,该分量进行时域离散后为:式(4)中的第二项表示与开关时变有源电路接口处电流,该分量进行时域离散后为: 式中,P为一个周期的离散点数;N为谐波次数;M为电压或电流的个数;uia表示第i个电压量的第a个离散值;imb表示第m个电流量的第b个离散值。
这样uia可写为: 如果记: 则 上式写为矩阵形式为: 时变线性无源网络稳态时的周期时间序列模型
如果按上式直接进行计算,其数据存储量是十分惊人的。对一个M端非时变线性无源网络,一个周期的离散点数为P,数据存储量为M2P2。如果按上式直接进行计算,其数据存储量是十分惊人的。对一个M端非时变线性无源网络,一个周期的离散点数为P,数据存储量为M2P2。 事实上,该矩阵具有一个很好的性质,它的对角元素是相等的,可写为如下形式: 这样,对一个M端非时变线性无源网络,一个周期的离散点数为P,数据存储量变为M2(2P一1),其数据存储量大大减少。
3.2 通用电路建模与仿真方法 1. 网络分析的基本概念和基本定律 任何集中参数网络都服从三条基本定律:基尔霍夫电压定律(KVL)、基尔霍夫电流定律(KCL)和元件定律(支路特性)。KVL和KCL是网络中各支路电压、电流的约束。元件定律通常指元件的伏安特性。根据三条基本定律就可以建立求解网络的方程式。 网络拓扑分析从理论上已经十分成熟,电路的计算机仿真程序就是在此基础上建立的。为了便于计算机的存储和加工,通常将网络拓扑图用矩阵来描述,常用的有关联矩阵、回路矩阵和割集矩阵。
2 改进节点法 节点分析法在处理独立电压源和受控电压源时,由于支路的导纳为无穷大,不能直接写出支路导纳矩阵Yb,所以现在通常采用所谓“改进节点法”或混合节点法。改进节点分析法在节点分析法的基础上,保留原来的未知量节点电压,增加独立电压源和受控电压源支路的电流作为电路的另一组未知量。利用改进节点分析法的一个方便之处是,欲求得流过某一条支路的电流,只要在该支路中串联一个零值独立电压源就行。即零值独立电压源具有电流表的功能。该方法成为目前电路的计算机辅助分析中应用最为广泛的一种,PSpice就是采用改进节点分析法。 改进节点分析法最终列出的方程式具有形式:WX=K,其中W为系数矩阵,K为右端向量,X为未知向量。
把支路电流方程分块成 式中下标代表电路器件:Y为导纳元件;J为电流源;E为电压源;VC为电压控制电流源VCCS;CC为电流控制电流源CCCS;EC为电压源支路电流控制的电流源IECCS;VV为电压控制电压源VCVS;CV为电流控制电压源CCVS;EV为电压源支路电流控制的电压源IECVS。
按基耳霍夫电流定律列出各节点的电流平衡方程按基耳霍夫电流定律列出各节点的电流平衡方程 方程组的未知变量除了五个节点电位v1~v5外,尚有iE1,iE2,共7个未知变量,因此需根据电压源支路特性再补上两个支路电压方程:
合并成一个方程组,并将各支路电阻用电导表示写成矩阵形式:合并成一个方程组,并将各支路电阻用电导表示写成矩阵形式: 记做: vn为电路的各节点(除参考节点外)电位矢量,iE是电压源支路电流矢量,I为独立电流源的电流矢量,E是独立电压源的电压矢量。系数矩阵A可以写成分块形式
Yn称为节点导纳矩阵。为了便于计算机自动形成A和B,下面详细讨论一下各电路元件的参数对A和B的贡献。Yn称为节点导纳矩阵。为了便于计算机自动形成A和B,下面详细讨论一下各电路元件的参数对A和B的贡献。
(1)电阻R 假设电阻R1接在节点1和3之间 对节点1:其电流 对节点3:其电流 因此,它对a11,a13,a31,a33四个元素有贡献,其贡献值为: a11=a13= G1,a31=a33=-G1,即: 对于接在任意两个节点e和f之间的电导G=1/R,对矩阵A的贡献 当一个节点上有几个电阻相连接时,矩阵Yn中与该节点相关的元素值为各电阻分别对该元素贡献值的代数和。可见,Yn中各对角元素的值为与该对应节点相连的各电阻值倒数(电导)之和,而非对角元素值为对应两节点间所接电阻值倒数的负值。
(2) 独立电压源 参照例子,有电压方程v2-v5=E2,同样可以推得接在任意节点e,f间的电压源将对子阵D,D1以及B有贡献。
(3)独立电流源 按上述相同的方法可以推得独立电流源仅对B中与它相连的节点所对应的行的元素有贡献,见下图。
(4)电压控制电流源VCCS 例子中的VCCS对节点e和f的电流贡献分别为G(vj-vk)和- G(vj-vk),所以它对Yn中相应元素的贡献如下所示。
(5)电压控制电压源VCVS 如下图所示,若Eef是受控源,受(vj-vk)的控制,即 与独立电压源相仿,需增加一个未知变量iE。因此矩阵应相应地也应增加一个方程,即上式。同时在与节点e和f相应的节点方程中也应增加iE的贡献。因此,上述电压控制电压源对A和B的贡献见下。
(6)电流控制电流源CCCS 设支路e,f是一个受控电流源。它的电流受支路j,k的电流ijk控制,控制系数为β。对于一般的阻抗支路,ijk并不包含在式(3-25)的未知变量中,为方便起见且不是一般性,可以在支路j,k中插入一个电势为零值的独立电压源E,这样并不影响原电路的特性和解。但是增加一个零值电压源E后,ijk作为流经E的电流就成为式(3-25)的一个未知变量,当然方程组也相应地增加了一个方程,控制支路也相应地由支路j,k变为支路m,k了。CCCS对A的贡献也可以方便地求得。
(7)电流控制电压源CCVS 与CCCS相仿,对于CCVS也可以在控制支路中插入一个零值电势的电压源E,由于受控电压源的存在,式(3-25)的未知变量中也应包括它的电流,所以这时ief及ijk都包含在未知变量中。因为ve-vf=γijk,很容易推得CCVS对A的贡献。