专利名称：三重四象限整流器仿真方法
技术领域：
本发明涉及仿真技术领域，更具体地说涉及半实物仿真领域。
背景技术：
随着计算机软、硬件的飞速发展，计算机仿真得到了越来越广泛的应用。半实物仿 真同其他类型的仿真方法相比具有更高真实度的可能性，是仿真技术中置信度最高的一种 仿真方法。交流传动是我国铁道牵引动力发展的一个重要方向，为了给设计开发电力机车 交流传动控制装置提供现代化的计算机仿真测试手段，常使用dSPACE实时仿真系统平台 来研究开发交流传动半实物仿真系统。交流传动系统中的整流器在使用过程中，应在直流侧提供平直的直流电流和直流 电压，且使交流供电网保持很好的功率因数。四象限整流器在PWM信号控制下能够很好地 满足这些要求，并且具有较快的动态响应。在交流传动系统中，为了使变压器一次侧电流的等效干扰电流最小，一般采用多 个四象限整流器并联的方式向异步牵引电机系统供电。每个整流器的控制调制信号的三角 载波互相错开一定的角度，可以使得四象限整流器的输入电流的高次谐波互相错开，在变 压器一次侧电流的谐波总量中部分抵消，所以使用三重四象限整流器。参见图1，四象限整流器系统主要由变压器次边、整流器、中间电路和PWM控制器 组成。三重四象限整流器系统由三个四象限整流器并联而成。为了建立系统级仿真数学模 型，可不考虑GTO和二极管的换流过程，把整流元件看成理想开关，忽略其中的保护电路， 分别用一个电感和一个电阻元件代替变压器的次边漏感和电阻。交流传动半实物仿真试验系统用来开发和验证四象限整流器控制算法和逆变器 控制算法，发现可能的问题并确定最佳参数，节约开发成本，缩短研发周期，提升公司的核 心竞争力。参见图2，现有技术中，利用MATLAB/SIMULINK软件对系统进行仿真。三重四象限 整流器的仿真模型由SIMULINK的SimPowerSystems工具箱中的模块直接构成。选取Power Electronics元件中的Universal Bridge，通过设置作为2桥臂的整流器，取三个整流器模 块，加上中间电容电感，构成整个四象限整流器子系统。但是，使用MATLAB工具箱中模块搭建的三重四象限整流器仿真模型，实现步骤繁 琐，占用系统资源多，实时仿真时采样时间要求比较长，实时性不高。而且使用工具箱中模 块搭建的三重四象限整流器仿真模型，在dSPACE中进行实时仿真时，难以解算，影响计算 结果的精度，甚至出现错误，准确性不高。
发明内容
有鉴于此，本发明实施例提供一种三重四象限整流器系统仿真方法，以提高仿真 时的实时性。为了实现上述目的，本发明提供以下方案
一种三重四象限整流器仿真方法，其特征在于，包括获取PWM控制信号与三重四象限整流器系统主电路拓扑结构的对应关系；建立三重四象限整流器系统的状态数学表达；获得所述三重四象限整流器系统的状态数学表达在控制信号下的状态描述方
程；对所述状态描述方程进行状态方程判断，获得状态方程矩阵及三重四象限整流器 系统的数学模型。优选的有，上述方法中，所述获取PWM控制信号与三重四象限整流器系统主电路 拓扑结构的对应关系具体为获取各种PWM控制信号下的整流元件的导通状态；获取与所述整流元件的导通状态对应的主电路拓扑结构；建立PWM控制信号与主电路拓扑结构的对应关系。优选的有，上述方法中，所述建立三重四象限整流器系统的状态数学表达包括定义A (n) = Tnl-Tn2，表示功率元件Tnl和Τη2开关情况；定义B (η) = Τη3_Τη4， 表示功率元件Τη3和Τη4开关情况；
1，Tnl开通，Τη2关断 0 ,Tnl关断，Τη2关断 Β(η) = < -l，Tnl关断，Τη2开通其中，Α(η):
1，Τη3开通，Τη4关断 0，Τη3关断，Τη4关断 -1，Τη3关断，Τη4开通直流电压和整流器的交流侧电流的接入关系为当Usn = +Ud 时，用 A(η)，B (η)表示为A (η) = 1，Β(η) =_1，或六(11) = Ο,Β(η) = 0 且满足 iNn > = 0，或 A (η) = 1，Β(η) =0 且满足 iNn >= 0，或 A(n) = 0，B(n) = -1 且满足 iNn >= 0 ；当Usn = 0 时，用 A(n)，B(n)表示为=A(η) = _1，Β(η) = -1，或 Α(η) = 1，Β(η)= 1，或 Α(η) =-1，Β(η) = 0 且满足 iNn > 0，或 A(η) = Ο,Β(η) = 1 且满足 iNn > 0，或 A(η) =1，Β(η) = 0 且满足 iNn < 0，或 Α(η) = 0，Β(η) = -1 且满足 iNn < 0 ；当Usn = -Ud 时，用 Α(η)，B(η)表示为Α(η) = -1，B(η) = 1，或 Α(η) = 0，B(η) =0 且满足 iNn < 0，或 Α(η) = -1，Β(η) = 0 且满足 iNn <= 0，或 Α(η) = Ο,Β(η) = 1 且 满足 iNn < = 0 ；
构建三重四象限系统状态描述方程 其中，输入变量为，变压器次边电压UNn (第η个整流器)、中间直流侧负载电流iL、 变压器次边电流iNn、直流电压Ud、滤波电流i2，滤波电感电压u2、第η个整流器决定的电感 值LL(n)和第η个整流器决定的电容值Cd (η)；其中，η为第η个整流器系统，η = 1,2,3；
_2]当 Usn = +Ud 时，fa㈨= 1/Q/、
\ LLiri) = 0当Usn = O时，㈨吣根据各整流器的接通情况获得Tnl，Tn2, Τη3, Τη4的值，计算出Α(η)，B(η)的值， 再结合电感电流的正负，判断出Ud的接入情况，根据Ud的接入情况选择各整流器的LL (η)、 Cd(η)的值，代入所述状态描述方程，获得三重四象限整流器系统的状态数学表达。优选的有，上述方法中，所述状态方程判断包括对所述状态描述方程进行控制信号及电流的判断，得到控制信号和电流的情况；对所述控制信号及电流的情况进行逻辑判断及选择，得到所述状态描述方程的矩 阵和判断信号。优选的有，上述方法中，所述状态方程判断中采用adms算法对所述的状态描述方 程进行解算。本发明实施例在实时仿真中，根据不同PWM控制信号下的功率元件的不同导通情 况，结合电流接入情况，得到变压器交流侧的接入情况，再结合构建好的状态数学表达得到 此种情况下的三重四象限整流器系统的状态描述方程。本发明实时仿真方法将三重四象限 整流器的系统状态空间模型建立问题，简化为根据开关元件导通情况进行相关参数选择的 问题，保证了准确性，降低了数学建模的复杂性，实现步骤简洁，占用系统资源较少，实时仿 真时采样时间较短，具有很高的实时性。对状态方程采用adms算法进行解算，能解决零状 态时的收敛性，快速准确的解算出微分方程的值，保证了解算的速度和精度。
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。
6
图1为三重四象限整流器系统主电路拓扑结构；图2为MATLAB工具箱中模块搭建的三重四象限整流器仿真模型；图3为本发明中步骤流程图；图4为本发明步骤1的流程图；图5为本发明步骤4的模块示意图；图6为本发明步骤4的流程图；图7为本发明步骤401的模块示意图；图8为本发明步骤402的模块示意图；图9为本发明采用的公式解微分方程模块示意图；图10为本发明采用adms算法解微分方程模块一示意图；图11为本发明采用adms算法解微分方程模块另一示意图。
具体实施例方式下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。本发明实施例一。参见图3，本发明三重四象限整流器系统仿真方法包括步骤1，获取PWM控制信号与三重四象限整流器系统主电路拓扑结构的对应关系；步骤2，建立三重四象限整流器系统的状态数学表达；步骤3，获得所述三重四象限整流器系统的状态数学表达在控制信号下的状态描 述方程；步骤4，对所述状态描述方程进行状态方程判断，获得状态方程矩阵及三重四象限 整流器系统的数学模型。每一重整流器存在三种情况，三重整流器系统就存在33 = 27种情况。也就是说， 根据Ud分别接入Usl，Us2和Us3，即Usn = Ud，0，-Ud，其中η = 1,2, 3，可以得到33 = 27种情 况，因此需要分别列出这27种情况下的不同主电路拓扑结构，即进行步骤1。参见图4，步 骤1又进一步包括步骤101，获取各种PWM控制信号下的整流元件的导通状态，即Tnl，Τη2, Τη3, Τη4 的导通状态。步骤102，获取与所述整流元件的导通状态对应的主电路拓扑结构；步骤103，建立PWM控制信号与主电路拓扑结构的对应关系。步骤2建立三重四象限整流器系统的状态数学表达，包括首先，进行定义；η代表 第η个整流器系统，其中η = 1，2，3 ；值1表示GTO功率元件开通，值0表示GTO功率元件 关断。令A (n) = Tnl-Tn2，表示GTO功率元件Tnl和Τη2开关情况，一个桥臂的两个GTO 不能同时开通，则Α(η)的可能取值为1，0和-1，含义为‘1，Tnl开通，Tn2关断 A(n) = | 0，Tnl 关断，Τη2关断 -1，Tnl关断，Τη2开通令B (η) = Τη3-Τη4，表示GTO功率元件Τη3和Τη4开关情况。一个桥臂的两个GTO 不能同时开通，Β(η)的可能取值为1，0和-1，含义为根据各GTO功率元件的导通情况与直流电压接入各整流器之间的关系和上述的 定义，得到直流电压和整流器的交流侧的接入关系的表达1)当 Usn = +Ud 时，用 A(η)，B(η)表示即为Α(η) = 1，Β(η) =-1，或A (η) = 0，B (η) = 0 且满足 iNn >=0，或Α(η) = 1，Β(η) = 0 且满足 iNn >= 0，或A (η) = 0，B (η) = _1 且满足 iNn > = 0 ;2)当 Usn = O 时，用 A(η)，B(η)表示即为Α(η) =-1，Β(η) =-1，或Α(η) = 1，Β(η) = 1，或A (η) = -1，Β(η) =0 且满足 iNn > 0，或A (η) = 0，B (η) = 1 且满足 iNn > 0，或A (η) = 1，B (η) =0 且满足 iNn < 0，或A (η) = 0，Β(η) = _1 且满足 iNn < 0，3)当 Usn = -Ud 时，用 A(η)，B(η)表示即为Α(η) =-1，Β(η) = 1，或A (η) = 0，B (η) =0 且满足 iNn < 0，或A (η) = -1，Β(η) =0 且满足 iNn <=0，或A (η) = 0，B (η) = 1 且满足 iNn <=0，LL(η)、Cd(η)为第η个整流器决定的电感及电容值，这两个量的情况根据对应的 整流器的导通情况决定。每个整流器的电感及电容值存在三种情况当Usn = +Ud 时，
‘1，Τη3开通，Τη4关断 B(n) = J 0，Tn3关断，Τη4关断 -1，Τη3关断，Τη4开通当 Usn = -Ud 时，当 Usn = 0 时，
构建状态描述方程，为了建立整流器系统的数学模型，选择状态描述方程的状态 变量为变压器次边电压UnJ第η个整流器)、中间直流侧负载电流iL、变压器次边电流iNn、 直流电压Ud、滤波电流i2，滤波电感电压u2、第η个整流器决定的电感值LL(η)和第η个整 流器决定的电容值Cd (η)；其中，η为第η个整流器系统，η = 1,2,3；根据不同的电路拓扑，得出不同情况下的系统状态方程描述为 以上为建立的三重四象限整流器系统的状态数学表达，该数学模型的输出为三重 四象限整流器系统的状态描述方程。在不同的PWM控制信号下，根据各整流器的接通可得Tnl，Tn2, Τη3, Τη4的值，由 此可计算出A(η)，B(η)的值，再结合电感电流的正负，判断出Ud的接入情况，即Usn的值， 即可选择各整流器的LL (n)、Cd (η)的值，将PWM控制信号下各整流器的LL (η)、Cd (η)的值， 带入所述状态描述方程，获得该PWM控制信号下三重四象限整流器系统的状态描述方程， 即完成了三重四象限整流器系统的数学模型的构建。本发明实施例一在实时仿真过程中，将三重四象限整流器的系统状态空间模型建 立问题，简化为根据开关元件导通情况进行相关参数选择的问题，保证了状态数学描述的 准确性和真实性，降低了数学建模的复杂性，实现步骤简洁，占用系统资源较少，实时仿真 时采样时间较短，具有很高的实时性。本发明实施例二。参见图5，步骤4在判断过程中，输入PWM控制信号以及电流，输出状态描述方程的 矩阵A及状态信号S,即对应的27种情况中的某一种。参见图6，步骤4包括
步骤401，对得到的状态描述方程进行控制信号及电流的判断，得到控制信号情况 和电流的情况。参见图7，输入为PWM控制信号下的功率元件的接通情况和电流情况，输出 为控制信号情况和电流的情况，具体的，得到的控制信号情况为步骤2中所定义的Al，Bi， A2，B2，A3，B3的值，得到的电流情况为电流大于、等于或小于零。步骤402，参见图8，对得到的控制信号及电流的情况进行逻辑判断及选择，得到 状态描述方程的矩阵A和状态信号S ；对得到的状态信号S进行选择，确定是27种情况中 的第几种，得到选择结果。其中，输入为步骤401得到的控制信号和电流的情况，进行逻辑 判断后输出状态信号S，然后对输出状态信号S进行选择，判断为27种情况中的第几种。其 中，参见图9，本发明实施例二在逻辑判断中，接收步骤401输出的控制信号及电流的情况， 进行计算后得到此PWM控制信号下的三重四象限整流器系统的状态描述方程，对得到的状 态描述方程进行逻辑判断后，再进行与或的运算，来确定所述情况的编号数字，如2，5，6分 别对应第2种，第5种，第6种的状态情况，输出这个状态信号S，即这个数字S。对得到的 状态信号S进行选择，判定为27种情况中的第几种，即得到了三重四象限的数学模型。本发明实施例在步骤402中，参见图10，在对得到的状态方程进行解算时，采用公 式
UU =(Ι-—55ΑΤ\χη+-(55Βυ +1-59/η+37/η_1 -9/ _2))在步骤402中，对得到的状态方程采用adms算法进行解算，隐式解算公式x +1 =X11 +^(55/ +1 -59/ +37/^ -9/ _2)该算法中存在代数环，不能用于实时仿真，因此进行进一步变换。根据公式fn+1 = Axn+1+BUn+1，将该公式代入上述公式中，得到
_]U ~γΑ55— 59/ +37/— -9/ _2)进行进一步化简，得到x +1 =(/-^ 55Ay' (xn + ^(55BUn+l 一 59/ + 37/ ., - 9/ _2))。参见图11，在用上述公式进行逻辑判断时，输入为PWM控制信号下的直流电压Ud，
状态描述方程B，Al至A27计算(/ + ^rSSJr1得到的A，状态信号S，输出为状态方程Xi+1。
24本发明实施例二在对状态描述方程的状态方程判断中，考虑到零状态时的收敛 性，采用adms算法进行解算，能解决零状态时的收敛性，快速准确的解算出微分方程的值， 保证了解算的速度和精度，保证了解算的准确性。以上所述仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在 本发明的精神和原则之内所作的任何修改、等同替换、改进等，均包含在本发明的保护范围 内。
权利要求
一种三重四象限整流器仿真方法，其特征在于，包括获取PWM控制信号与三重四象限整流器系统主电路拓扑结构的对应关系；建立三重四象限整流器系统的状态数学表达；获得所述三重四象限整流器系统的数学模型在控制信号下的状态描述方程；对所述状态描述方程进行状态方程判断，获得状态方程矩阵及三重四象限整流器系统的数学模型。
2.根据权利要求1所述的方法，其特征在于，所述获取PWM控制信号与三重四象限整流 器系统主电路拓扑结构的对应关系具体为获取各种PWM控制信号下的整流元件的导通状态； 获取与所述整流元件的导通状态对应的主电路拓扑结构； 建立PWM控制信号与主电路拓扑结构的对应关系。
3.根据权利要求1所述的方法，其特征在于，所述建立三重四象限整流器系统的状态 数学表达包括定义A (η) 二！^^！^，表示功率元件他丨和！“！^开关情况；定义^!!) =Τη3-Τη4，表示 功率元件Τη3和Τη4开关情况；其中，=`1，Tnl幵通，Τη2关断 0，Tnl关断，Τη2关断 Β(η) = -l，Tnl关断，Τη2开通`1，Τη3幵通，Τη4关断 0，Τη3关断，Τη4关断 -1，Τη3关断，Τη4开通直流电压和整流器的交流侧电流的接入关系为 当Usn = +Ud时，用Α(η)，Β(η)表示为A(η) = 1，Β(η) = -1，或 Α(η) = Ο,Β(η) = 0 且满足 iNn >= 0，或 A(n) = l，B(n)= 0 且满足 iNn >= 0，或 A(n) = 0, B(n) = -1 且满足 iNn >= 0 ；当 Usn = 0 时，用 A(n)，B(η)表示为Α(η) = -1，B(η) = -1，或 Α(η) = 1，B(η) = 1， 或 Α(η) = -1，B (η) = 0 且满足 iNn > 0，或 Α(η) = 0, B (η) = 1 且满足 iNn > `0，或 Α(η) =1，Β(η) = 0 且满足 iNn < 0，或 Α(η) = 0，Β(η) = -1 且满足 iNn < 0 ；当 Usn = -Ud 时，用 Α(η)，B(η)表示为Α(η) = -1，B(η) = 1，或 Α(η) = 0，B(η) = 0 且满足 iNn < 0，或 A(η) =-1，Β(η) = 0 且满足 iNn < = 0，或 A(η) = Ο,Β(η) = 1 且满足 iNn < = 0 ；构建三重四象限系统状态描述方程 其中，输入变量为，变压器次边电压UNn(第n个整流器)、中间直流侧负载电流iL、变 压器次边电流iNn、直流电压Ud、滤波电流i2，滤波电感电压u2、第n个整流器决定的电感值 LL(n)和第n个整流器决定的电容值Cd(n)；其中，n为第n个整流器系统，n = 1，2，3 ；i{n) = 0. ？根据各整流器的接通情况获得Tnl，Tn2,Tn3,Tn4的值，计算出A (n)，B (n)的值，再结合 电感电流的正负，判断出Ud的接入情况，根据Ud的接入情况选择各整流器的LL(n)、Cd(n) 的值，代入所述状态描述方程，获得三重四象限整流器系统的数学模型。
4.根据权利要求1所述的方法，其特征在于，所述状态方程判断包括对所述状态描述方程进行控制信号及电流的判断，得到控制信号和电流的情况； 对所述控制信号及电流的情况进行逻辑判断及选择，得到所述状态描述方程的矩阵和 判断信号。
5.根据权利要求1所述的方法，其特征在于，所述状态方程采用adms法对所述的状态 方程进行解算。
全文摘要
本发明公开了一种三重四象限整流器仿真方法，所述方法包括获取PWM控制信号与三重四象限整流器系统主电路拓扑结构的对应关系；建立三重四象限整流器系统的状态数学描述；获得所述三重四象限整流器系统的数学模型在控制信号下的状态描述方程；对所述状态描述方程进行状态方程判断，获得状态方程矩阵及三重四象限整流器的数学模型。本发明实时仿真方法将三重四象限整流器的系统状态空间模型建立问题，简化为根据开关元件导通情况进行相关参数选择的问题，保证了准确性，降低了数学建模的复杂性，实现步骤简洁，占用系统资源较少，实时仿真时采样时间较短，具有很高的实时性。
文档编号G06F9/455GK101895210SQ200910171929
公开日2010年11月24日 申请日期2009年9月22日 优先权日2009年9月22日
发明者刘可安, 应婷, 王坚 申请人:株洲南车时代电气股份有限公司