【技术实现步骤摘要】
本专利技术涉及电力仿真，尤其涉及一种传输线电磁暂态仿真计算方法及装置。
技术介绍
1、时域有限差分法(fdtd)被广泛应用于解决各种电磁场问题，可以通过求解电报方程获得有损多导体传输线瞬态响应。fdtd方法具有灵活性与简单性，在建模时不需要对传输线进行解耦，可以直接处理非均匀传输线。此外，它也被成功地拓展到对具有频变参数的传输线方面的问题的分析。
2、尽管fdtd具有很多优点，该方法却受到了克朗条件(courant–friedrichs–lewycondition，即cfl条件，确保数值格式收敛稳定)的限制。在克朗条件的限制下，时间段必须足够小才能保证所要求的精度。但是一些电气暂态过程相对缓慢且持续时间较长，比如，超高压输电线路操作暂态上升时间可达几百微秒。采用经典的fdtd对其进行仿真会受限于cfl条件而十分耗时。基于交变隐式差分方向方法的时域有限差分法(adi-fdtd)克服了cfl条件的限制，但是这一方法存在严重的数值色散和分裂误差。而采用精确的时步积分法求解fdtd模型可以在大的时步下获得稳定而精确的结果，但是该方法只适用于线性电路，在求解非线性问题时，首先要对其进行线性化处理，并重复计算雅各比矩阵，因此，对于大规模的非线性电路，该方法的效率会大打折扣。
3、基于谱延迟预估校正的时域有限差分法(sdc-fdtd)利用谱延迟预估校正方法(spectral deferred correction，sdc)来提高fdtd在大时步线性电路的稳定性和精度。当输电线路与外电路耦合时，模型在数学形式上由常微分方
技术实现思路
1、本专利技术所要解决的技术问题是：提供一种传输线电磁暂态仿真计算方法及装置，能够提升计算效率和稳定性。
2、为了解决上述技术问题，本专利技术采用的技术方案为：
3、一种传输线电磁暂态仿真计算方法，包括步骤：
4、确定多耦合的传输线电报方程，并根据预设yee氏网格使用中心差分法对所述多耦合的传输线电报方程进行离散，得到传输线的半离散化模型；
5、基于所述传输线的半离散化模型将多电阻网络与所述多耦合的传输线相连，得到耦合模型的常微分方程组，并根据基尔霍夫定律确定电路系统的代数方程；
6、根据所述耦合模型的常微分方程组和所述电路系统的代数方程构成微分代数方程，并对所述微分代数方程进行离散，得到离散后的微分代数方程；
7、对所述离散后的微分代数方程进行谱延迟预估校正，得到校正后的方程，并使用jfnk算法结合krylov子空间法对所述校正后的方程进行计算，得到电磁暂态的解。
8、进一步地，所述确定多耦合的传输线电报方程包括：
9、
10、
11、v(x,t)＝[v1(x,t),v2(x,t),……,vn(x,t)]；
12、i(x,t)＝[i1(x,t),i2(x,t),……,in(x,t)]；
13、式中，t表示时间，v(x,t)表示传输线上每一节点的电压，i(x,t)表示传输线上每一节点的电流，vn(x,t)表示传输线上节点n的电压，in(x,t)表示传输线上节点n的电流，x表示传输线上的位置，r(x)表示传输线上位置x处的电阻，l(x)表示示传输线上位置x处的电感，g(x)表示传输线上位置x处的导纳，c(x)表示传输线上位置x处的电容。
14、进一步地，所述传输线的半离散化模型具体为：
15、
16、
17、
18、
19、式中，in表示预设yee氏网格中离散节点n处的电流，ln表示在预设yee氏网格的离散节点n处，用差分方式来近似表示的电感，in+1表示预设yee氏网格中离散节点n+1处的电流，vn+1表示预设yee氏网格中离散节点n+1处的电压，vn表示预设yee氏网格中离散节点n处的电压，rn表示在预设yee氏网格的离散节点n处，用差分方式来近似表示的电阻，cn表示在预设yee氏网格的离散节点n处，用差分方式来近似表示的电容，gn表示在预设yee氏网格的离散节点n处，用差分方式来近似表示的导纳，vm表示预设yee氏网格中离散节点m处的电压，cm表示在预设yee氏网格的离散节点m处，用差分方式来近似表示的电容，im-2表示预设yee氏网格中离散节点m-2处的电流，im-1表示预设yee氏网格中离散节点m-1处的电流，im表示预设yee氏网格中离散节点m处的电流，cm表示在预设yee氏网格的离散节点m处，用差分方式来近似表示的电容，gm表示在预设yee氏网格的离散节点m处，用差分方式来近似表示的导纳。
20、进一步地，所述基于所述传输线的半离散化模型将多电阻网络与所述多耦合的传输线相连，得到耦合模型的常微分方程组包括：
21、将所述传输线的半离散化模型转换为矩阵形式的常微分方程组；
22、确定电阻的v-i特性；
23、考虑一含有多个电阻和电压源的电阻网络，作为多电阻网络；
24、基于所述矩阵形式的常微分方程组、所述电阻的v-i特性将多电阻网络与所述多耦合的传输线相连，得到耦合模型的常微分方程组。
25、进一步地，所述将所述传输线的半离散化模型转换为矩阵形式的常微分方程组包括：
26、
27、x'＝[v1,i1,v2,i2,…,vm-1,im-1,vm]t；
28、
29、
30、式中，x′表示第一列向量，a表示系数矩阵，f表示第二列向量；
31、所述确定电阻的v-i特性包括：
32、vrs-vre-f(ir)＝0；
33、式中，vrs表示负载电阻始端电位，vre表示负载电阻末端的电位，f(ir)表示电阻的电压与电流关系；
34、所述基于所述矩阵形式的常微分方程组、所述电阻的v-i特性将多电阻网络与所述多耦合的传输线相连，得到耦合模型的常微分方程组包括：
35、
36、x'＝[v1,i1,v2,i2,…,vm-1,im-1,vm]；
37、
38、式中，b表示外部施加影响的矩阵，z表示向量，表示第一段传输线始端电压，表示第一段传输线末端电压，表示第一段传输线始端电流，表示第一段传输线末端电流，表示多电阻网络各个负载以及电源的v-i特性方程各个参数组成的向量。
39、进一步地，所述根据基尔霍夫定律确定电路系统的代数方程包括：
40、根据所述电阻的v-i特性和所述多电阻网络基于基尔霍夫定律得到电本文档来自技高网...
【技术保护点】
1.一种传输线电磁暂态仿真计算方法，其特征在于，包括步骤：
2.根据权利要求1所述的一种传输线电磁暂态仿真计算方法，其特征在于，所述确定多耦合的传输线电报方程包括：
3.根据权利要求2所述的一种传输线电磁暂态仿真计算方法，其特征在于，所述传输线的半离散化模型具体为：
4.根据权利要求3所述的一种传输线电磁暂态仿真计算方法，其特征在于，所述基于所述传输线的半离散化模型将多电阻网络与所述多耦合的传输线相连，得到耦合模型的常微分方程组包括：
5.根据权利要求4所述的一种传输线电磁暂态仿真计算方法，其特征在于，所述将所述传输线的半离散化模型转换为矩阵形式的常微分方程组包括：
6.根据权利要求4所述的一种传输线电磁暂态仿真计算方法，其特征在于，所述根据基尔霍夫定律确定电路系统的代数方程包括：
7.根据权利要求5所述的一种传输线电磁暂态仿真计算方法，其特征在于，所述根据所述耦合模型的常微分方程组和所述电路系统的代数方程构成微分代数方程包括：
8.根据权利要求1所述的一种传输线电磁暂态仿真计算方法，其特征在于，
9.根据权利要求7所述的一种传输线电磁暂态仿真计算方法，其特征在于，所述离散后的微分代数方程具体为：
10.一种传输线电磁暂态仿真计算装置，其特征在于，包括：
...【技术特征摘要】
1.一种传输线电磁暂态仿真计算方法，其特征在于，包括步骤：
2.根据权利要求1所述的一种传输线电磁暂态仿真计算方法，其特征在于，所述确定多耦合的传输线电报方程包括：
3.根据权利要求2所述的一种传输线电磁暂态仿真计算方法，其特征在于，所述传输线的半离散化模型具体为：
4.根据权利要求3所述的一种传输线电磁暂态仿真计算方法，其特征在于，所述基于所述传输线的半离散化模型将多电阻网络与所述多耦合的传输线相连，得到耦合模型的常微分方程组包括：
5.根据权利要求4所述的一种传输线电磁暂态仿真计算方法，其特征在于，所述将所述传输线的半离散化模型转换为矩阵形式的常微分方程组包括：...
【专利技术属性】
技术研发人员：庄池杰，黄友聪，孙永森，傅智为，郑钟楠，张莹，林梓圻，吴晓杰，张世炼，彭利强，许军，
申请(专利权)人：国网福建省电力有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。