WO2015035758A1 - 一种基于多回路法的场路瞬态-瞬态耦合仿真方法 - Google Patents

Abstract

一种基于多回路法的场路瞬态-瞬态耦合仿真方法，包括以下步骤：1）、依照电机的有限元模型，根据多回路理论，建立电机的多回路模型，如下（1）~ （6）：2）通过调用有限元模型，在每一个系统仿真步长中，利用电机电磁场有限元模型求解每一个绕组、包括阻尼回路的自感、互感参数，然后将该参数返回到电机多回路模型中，对电机及外电路模型进行迭代求解，实现在系统仿真中耦合有限元多回路模型的联合仿真。本发明提供了一种有效适用于多个有限元模型场合、提升收敛性能的基于多回路法的场路瞬态-瞬态耦合仿真方法。

Description

一种基于多回路法的场路瞬态 -瞬态耦合仿真方法 技术领域

本发明涉及一种电机的场路耦合方法。

背景技术

现代功率电力电子系统中， 电机是主要的发电与用电设备之一。 为了在设 计、 分析、 优化等各个阶段， 对系统进行精确模拟， 必须在系统仿真过程中建 立详细可靠的电机模型。 由于电机的结构、 材料会对系统的谐波、 负载、 瞬变 等特性产生影响， 因此需要建立电机的电磁场模型耦合到系统仿真器中进行联 合仿真。

目前通用的场路耦合算法有两种， 一种是直接耦合算法， 这种算法将电路 模型耦合到电机有限元方程中进行求解， 因此只能实现一个电机有限元模型与 系统电路的耦合仿真。 另一种是间接耦合， 通过将有限元模型集成到系统电路 中实现瞬态-瞬态联合仿真。 目前这种方法通用的算法是采用戴维南法和诺顿 法。 戴维南法将电机的每个绕组等效成一个恒压源与阻抗模型， 通过时步算法 实时调用电压与阻抗参数实现联合仿真。 诺顿法则将电机等效成恒流源与电导 参数模型， 并进行实时迭代计算。

诺顿法与戴维南法其物理原理非常简单， 也是电路分析中最常见的等效技 术，因此在联合仿真中得到广泛应用。但是由于引进了激励源（电压源 /电流源）， 当系统中包含多个有限元模型时， 系统电路求解时容易造成由于该激励源变化 引起的数值计算不收敛问题， 尤其是当有限元模型端子之间存在相互连接的情 况时， 这中情况尤为明显。 而这样的应用场合又是目前系统中普遍存在的情况， 例如在三级无刷交流发电系统中， 励磁机的输出电压直接供给给主发电机的励 磁绕组。 因此目前现有技术的瞬态-瞬态联合仿真主要局限于单个电机有限元模 型与系统电路的联合仿真。

发明内容

为了克服已有的电机场路耦合方法的不能适用于多个有限元模型场合、收敛性能较差的不足，本发 明提供了一种有效适用于多个有限元模型场合、提升收敛性能的基于多回路法的场路瞬态 -瞬态耦合仿 真方法。 本发明解决其技术问题所采用的技术方案是：

一种基于多回路法的场路瞬态-瞬态耦合仿真方法， 包括以下步骤：

1 )、 依照电机的有限元模型， 根据多回路理论， 建立电机的多回路模型， 如下 ( 1 ) ~ (6)：

t - ) + ¾₊₁ ( _t - ) + R_{Rt t} + P _L (2)

άω + (6)

dt

其中， S为摩擦系数， /为转动惯量， 为阻尼绕组电流， 为电机绕组电流， M为绕组互感， 为转子位置， ¾为阻尼条电阻， r_em为电磁转矩， 7 为负载 转矩， W为电机转速， 为阻尼绕组磁链， W为电机绕组磁链；

2)通过调用有限元模型， 在每一个系统仿真步长中， 利用电机电磁场有限元模 型求解每一个绕组、 包括阻尼回路的自感、 互感参数， 然后将该参数返回到电 机多回路模型中， 对电机及外电路模型进行迭代求解， 实现在系统仿真中耦合 有限元多回路模型的联合仿真。

进一步， 所述步骤 2 ) 中， 采用有限元模型进行电感参数的求解如下： 2.1 )首先根据上一时刻各绕组激励， 求解电机内地电磁场分布和电磁转矩， 求解完成后， 保持定转子铁芯各三角单元中的相对磁导率 值不变；

2.2) 然后依次在各绕组、 导条中分别通入单位电流， 同时保持其余绕组、 导条为零激励， 此时各回路中的磁链即为各回路的互感参数， 将该电感矩阵参 数和电磁转矩提取后返回到系统仿真器中建立多回路电机模型， 并将该模型与 外部电路、 机械运动方程耦合进行求解。

本发明的有益效果主要表现在： 1、 实现了系统仿真中集成多个电机有限元 模型的瞬态-瞬态联合仿真技术， 解决了系统仿真中多有限元模型耦合仿真的技 术难点； 2、 采用多回路联合仿真技术使得联合仿真的收敛性能更好； 3、 采用 同步迭代技术， 实现瞬态模型与系统电路之间同步迭代， 提高了仿真精度和收 敛性能。

附图说明

图 1 是通过调用有限元模型， 实时迭代求解电机绕组与阻尼回路的电感参 数的流程图。

图 2是多回路场路耦合仿真流程图。

图 3 是采用交流电源时， 仿真结果与实验结果的电流波形比较示意图， 其 中， （a) 为采用多回路场路耦合仿真获得的空载启动过程绕组电流响应， （b ) 为实测空载启动电流响应， （c ) 为实测电流与仿真结果比较。

图 4是采用变频电源时， 仿真结果与实验结果的电流波形比较示意图， （a) 为采用多回路场路耦合仿真获得的绕组稳态电流响应波形， （b ) 为实测稳态电 流响应波形， （c ) 为实测电流与仿真结果比较。 图 5是三级无刷交流发电系统的示意图。

图 6 是仿真得到的主发输出电压、 与励磁机励磁电流响应波形的示意图， 其中， （a) 表示励磁机的励磁电流响应波形， （b) 表示空载建压过程中的电压 响应波形， （c) 表示突加负载时的电压瞬态响应波形。

具体实施方式

下面结合附图对本发明作进一步描述。

参照图 1~图6， 一种基于多回路法的场路瞬态-瞬态耦合仿真方法， 包括以 下步骤：

1)、 依照电机的有限元模型， 根据多回路理论， 建立电机的多回路模型， 如下

(1) ~ (6)： ^uw_k = ^Rw w_k

⁰ = ^RB_k ( -

人 +∑ ^ (4)

1 _r.，<3|

=÷['·] (5)

2^{L J} 3Θ ['ΐ+['·] δθ

j^ _{+ Bm = T} ― _Tl. (6)

dt ^{em L}

其中， S为摩擦系数， /为转动惯量， 为阻尼绕组电流， 为电机绕组电流， M为绕组互感， 为转子位置， ¾为阻尼条电阻， r_em为电磁转矩， 7 为负载 转矩， w为电机转速， 为阻尼绕组磁链， W为电机绕组磁链； 2)通过调用有限元模型， 在每一个系统仿真步长中， 利用电机电磁场有限元模 型求解每一个绕组、 包括阻尼回路的自感、 互感参数， 然后将该参数返回到电 机多回路模型中， 对电机及外电路模型进行迭代求解， 实现在系统仿真中耦合 有限元多回路模型的联合仿; 进一步， 所述步骤 2) 中， 采用有限元模型进行电感参数的求解如下：

2.1 )首先根据上一时刻各绕组激励， 求解电机内地电磁场分布和电磁转矩， 求解完成后， 保持定转子铁芯各三角单元中的相对磁导率 值不变；

2.2) 然后依次在各绕组、 导条中分别通入单位电流， 同时保持其余绕组、 导条为零激励， 此时各回路中的磁链即为各回路的互感参数， 将该电感矩阵参 数和电磁转矩提取后返回到系统仿真器中建立多回路电机模型， 并将该模型与 外部电路、 机械运动方程耦合进行求解。

本发明中， 根据多回路模型的定义， 电机电感参数列表的精确求解是建模 的关键， 目前通常的方法是通过离线计算得到电感参数关于位置的数据表， 然 后在电机的时域多回路模型中， 实时查询电感参数， 迭代求解电机模型。 这种 方法由于忽略电机的饱和因素， 因此难以对电机重载工况下的稳态、 瞬态性能 进行精确模拟。

为了弥补多回路模型难以考虑电机饱和特性的不足， 本方案通过调用有限 元模型， 实时迭代求解电机绕组与阻尼回路的电感参数， 从而实现外电路模型 与电机有限元模型的联合仿真。 在这种方法里， 在每一个系统仿真步长中， 利 用电机电磁场有限元模型求解每一个绕组、 包括阻尼回路的自感、 互感参数， 然后将该参数返回到电机多回路模型中， 对电机及外电路模型进行迭代求解， 其求解流程如图 1所示。

在采用有限元模型进行电感参数的求解时， 首先根据上一时刻各绕组激励， 求解电机内地电磁场分布和电磁转矩， 求解完成后， 保持定转子铁芯各三角单 元中的相对磁导率/值不变， 然后依次在各绕组、 导条中分别通入单位电流， 同 时保持其余绕组、 导条为零激励， 此时各回路中的磁链即为各回路的互感参数， 将该电感矩阵参数和电磁转矩提取后返回到系统仿真器中建立多回路电机模 型， 并将该模型与外部电路、 机械运动方程等耦合进行求解。 为了保证计算精 度， 在系统仿真器迭代求解过程中， 每一次迭代都将调用有限元模型进行参数 求解， 直到当前步长迭代完成后， 将当前绕组的激励和转子位置作为下一个步 长有限元计算的初始值。 经过如此循环迭代， 这样就实现了在系统仿真中耦合 有限元多回路模型的联合仿真。

由于在每一个步长内， 有限元模型均参与了系统仿真器的迭代计算， 因此 在该方法中， 有限元模型与系统仿真器之间实际上实现了同步仿真， 解决了通 常联合仿真中有限元模型落后系统仿真一个步长的问题。

采用多回路法建立电机有限元的联合仿真模型， 由于在模型中仅仅提供电 感、 电阻参数， 而不涉及有限元模型中计算出来的绕组电压和电流参数， 因此 在系统模型中进行联合仿真时具备非常好的收敛性能， 可以支持多个多回路模 型绕组之间的串联、 并联分析， 从而解决了系统仿真器中耦合多个有限元模型 时容易造成的收敛性问题。 另一方面， 由于绕组电感参数反映了电机的饱和、 空间谐波等因素， 因此可以精确对系统瞬态、 稳态性能进行仿真评估。

本方案采用系统仿真软件 Portunus和电磁场有限元仿真软件 EasiMotor建立 了实现了电机的多回路瞬态联合仿真， 如图 2所示。

为了验证该算法， 我们采用某感应电机对其进行了实验验证， 图 3、 图 4分 别为采用交流电源与变频电源时， 仿真结果与实验结果的电流波形比较， 显然 采用多回路算法的仿真结果与实验结果基本吻合。

为了验证该算法在多有限元模型耦合仿真中的应用， 我们对某三级无刷交 流发电系统进行了仿真验证， 系统组成如下图 5所示， 该系统由永磁机、 励磁 机、 主发、 旋转整流桥以及调压器组成， 在仿真模型中， 永磁机、 励磁机、 主 发采用有限元模型进行仿真。 图 6所示为仿真得到的主发输出电压、 与励磁机励磁电流响应波形， 系统首 先进行空载建压过程， 在 0.1s 时突加 75% 额定功率。 仿真结果的励磁电流与 突加负载时电压跌落等瞬态数据与实测结果基本一致。

Claims

权利要求书

1)、 依照电机的有限元模型， 根据多回路理论， 建立电机的多回路模型， 如下 (1) ~ (6)：

lw_k = ^Rw_kiw_k + ΡΨ ( _t - ― ₁ ) + ¾₊₁ ( _t - ) + R_{Rt t} + p .

= ^Lw_t w_t +∑ M _WtWj_w_ +∑ M _Wt人 +∑ _wtM

T =½ [ +

ji^_{+Bo) = T} __T (6)

dt ^{em L}

其中， 下标 表示第 个电机绕组或者阻尼条， 下标 W„、 L„和 分别表示电 机绕组、 阻尼绕组、 永磁体激励第《个回路， p表示求导； S为摩擦系数， /为 转动惯量， ^为阻尼绕组电流， ½为电机绕组电流， 示绕组漏感， Μ为绕 组互感， 为转子位置， ¾为阻尼条电阻， 为端环电阻， R_w为绕组电阻， r_em为电磁转矩， 7 为负载转矩， _Mw为电机绕组端电压， w为电机转速， /_L%阻 尼绕组磁链， ^为电机绕组磁链；

2)通过调用有限元模型， 在每一个系统仿真步长中， 利用电机电磁场有限元模 型求解每一个绕组、 包括阻尼回路的自感、 互感参数， 然后将该参数返回到电 机多回路模型中， 对电机及外电路模型进行迭代求解， 实现在系统仿真中耦合 有限元多回路模型的联合仿真。

2、如权利要求 1所述的基于多回路法的场路瞬态-瞬态耦合仿真方法，其特征在 于： 所述步骤 2) 中， 采用有限元模型进行电感参数的求解如下：

2.1 )首先根据上一时刻各绕组激励， 求解电机内地电磁场分布和电磁转矩， 求解完成后， 保持定转子铁芯各三角单元中的相对磁导率 值不变；

2.2)然后依次在各绕组、 导条中分别通入单位电流， 同时保持其余绕组、 导 条为零激励， 此时各回路中的磁链即为各回路的互感参数， 将该电感矩阵参数 和电磁转矩提取后返回到系统仿真器中建立多回路电机模型， 并将该模型与外 部电路、 机械运动方程耦合进行求解。