一种分布参数电路放电电流的频域计算方法及系统技术方案

本发明专利技术公开了一种分布参数电路放电电流的频域计算方法及系统，考虑到放电过程的物理特性，使用电阻、电感和开关来模拟放电电路，放电的发生等效为放电电路的开关闭合。开关闭合前分布参数电路的波过程可以使用经典传输线理论中的链参数矩阵在频域中计算；放电发生时，即开关闭合后，由于分布参数电路中分布电感和分布电容的能量不能瞬间释放，因此其中的能量同样影响放电的波过程。考虑到分布电感和分布电容在放电前一瞬间储存的能量，推导了频域中放电后分布参数电路的波过程。本发明专利技术可用于分析快前沿、高幅值脉冲下，线圈、传输线等结构发生绝缘击穿发电时的波过程，频域方法较时域方法占用更少时间资源、计算效率更高。计算效率更高。计算效率更高。

【技术实现步骤摘要】
一种分布参数电路放电电流的频域计算方法及系统


[0001]本专利技术涉及分布参数电路的波过程分析，特别涉及一种分布参数电路放电电流的频域计算方法及系统。

技术介绍

[0002]电气、电子设备面临的电磁威胁日益严重，快前沿、高幅值的电磁脉冲的产生方式越来越方便和便携，其可以通过架空线耦合进入线圈等设备。因此线圈、架空线等在该脉冲下发生绝缘击穿放电的概率也在增加。当脉冲的波长和电路的尺寸相当时，传输线理论是分析分布参数电路的有效工具。
[0003]由于绝缘击穿放电的瞬间相当于放电通路的导通，可能导致电路中产生新的过冲，进而存在威胁设备可靠性的可能。因此对于分布参数电路放电的波过程分析是十分重要的。文献《通信与信号电缆分布参数及放电特性的研究》研究了信号电缆放电时的放电波形，以及电缆长度、分布电感、分布电容对于放电波形的影响。然而，该文献采取的计算模型较简单，即使用集中电感和集中电容来代替电缆的分布电感和分布电容，仅仅将计算聚焦于放电回路的波形计算，忽略了放电前后分布参数电路的波过程。同时，该文献仅仅考虑了等效的集总参数电路在直流电源下、在振荡状态下的放电电流，而实际上的线圈、传输线等并不处于振荡状态下，且电源为快脉冲，因此该文献的方法对于分析快脉冲下分布参数电路放电的波过程是不完备的。

技术实现思路

[0004]本专利技术的目的在于提供一种分布参数电路放电电流的频域计算方法及系统，以解决现有技术存在的问题，本专利技术可用于分析快前沿、高幅值脉冲下，线圈、传输线等结构发生绝缘击穿发电时的波过程，本专利技术采用的频域方法较时域方法占用更少时间资源、计算效率更高。
[0005]为达到上述目的，本专利技术采用如下技术方案：一种分布参数电路放电电流的频域计算方法，包括：步骤一：对分布参数电路击穿放电发生前的波过程进行分析，得到放电发生前的波形；步骤二：拟合分布参数电路击穿放电发生时刻的电压、电流随距离的变化曲线；步骤三：结合步骤二所得变化曲线对分布参数电路击穿放电发生后的波过程进行分析，得到放电发生后的波形；步骤四：合并放电发生前的波形和放电发生后的波形，得到完整的分布参数电路放电的波过程。
[0006]进一步地，所述步骤一具体为：设放电发生的时间为t
d
，脉冲源的波形为，则放电发生前频域中用于激励的脉冲源为时刻0至时刻t
d
的脉冲源波形，即
，其中，表示击穿时刻前的激励源，t表示时间，表示信号的傅里叶变换，为阶跃函数；此时频域传输线方程为：其中，j为虚数单位，ω为角频率，x为位置变量，分别为电压相量、电流相量，Z,Y分别为传输线的阻抗矩阵、导纳矩阵；则此时分布参数电路首末端电压和电流关系即为：其中，为分布参数电路的传播常数，R、L
s
、G、C、j、ω分别为电阻、电感、电导、电容、虚数单位、角频率，L为分布参数电路首末端的距离，为分布参数电路的特征阻抗；分别为位置L、位置0处的电压相量，分别为位置L、位置0处的电流相量；在放电发生前的脉冲源激励下，结合分布参数电路首末端的边界条件，即能够求解击穿放电发生前的分布参数电路上的电流、电压，获得放电发生前的波形。
[0007]进一步地，所述步骤二中采用以幂函数为底的基函数拟合电流、电压在t
d
时刻随距离x的变化曲线，如下式所示：其中，y
v
、y
i
、、b
i
、i、n分别为电压变量、电流变量、拟合电压随x变化用到的系数、拟合电流随x变化用到的系数、角标、拟合阶数，通过链参数矩阵在长度为L的分布参数电路取N个样本，即得：
其中，Δx表示采样间隔，m为采样的样本量，k1表示第k1个样本编号；由此得到距离自变量对应的样本数据为，因变量电压对应的样本数据为，因变量电流对应的样本数据为，将样本数据代入以幂函数为底的基函数中，得到如下等式：的基函数中，得到如下等式：拟合电压随x变化用到的系数、拟合电流随x变化用到的系数b
i
即为：
令，，分别表示t
d
时刻电压随x的变化规律、t
d
时刻电流随x的变化规律，进而得分布参数电路击穿放电发生时刻的电压、电流随距离的变化曲线。
[0008]进一步地，所述步骤三具体为：设放电发生的时间为t
d
，脉冲源的波形为，则放电发生后频域中用于激励的脉冲源为时刻t
d
之后的脉冲源波形，即，表示击穿时刻后的激励源，其中，表示信号的傅里叶变换，为阶跃函数；此时频域传输线方程为：其中，L
S
为电感矩阵，C为电容矩阵；则首末端电压和电流关系为：则首末端电压和电流关系为：其中：其中：其中：、分别表示等效电压源、等效电流源，的自变量为
距离x，分别表示分布电压源、分布电流源，其表达式为：，分别表示在t
d
时刻电流随x的变化关系、电压随x的变化关系；在放电发生后的脉冲源激励下，结合分布参数电路首末端的边界条件，即能够求解击穿放电发生后的分布参数电路上的电流、电压，获得放电发生后的波形。
[0009]一种分布参数电路放电电流的频域计算系统，包括：放电发生前的波形获取模块：用于对分布参数电路击穿放电发生前的波过程进行分析，得到放电发生前的波形；变化曲线获取模块：用于拟合分布参数电路击穿放电发生时刻的电压、电流随距离的变化曲线；放电发生后的波形获取模块：用于结合所得变化曲线对分布参数电路击穿放电发生后的波过程进行分析，得到放电发生后的波形；分布参数电路放电的波过程获取模块：用于合并放电发生前的波形和放电发生后的波形，得到完整的分布参数电路放电的波过程。
[0010]进一步地，对分布参数电路击穿放电发生前的波过程进行分析，得到放电发生前的波形，具体为：设放电发生的时间为t
d
，脉冲源的波形为，则放电发生前频域中用于激励的脉冲源为时刻0至时刻t
d
的脉冲源波形，即，其中，表示击穿时刻前的激励源，t表示时间，表示信号的傅里叶变换，为阶跃函数；此时频域传输线方程为：其中，j为虚数单位，ω为角频率，x为位置变量，分别为电压相量、电流相量，Z,Y分别为传输线的阻抗矩阵、导纳矩阵；则此时分布参数电路首末端电压和电流关系即为：其中，为分布参数电路的传播常数，R、L
s
、
G、C、j、ω分别为电阻、电感、电导、电容、虚数单位、角频率，L为分布参数电路首末端的距离，为分布参数电路的特征阻抗；分别为位置L、位置0处的电压相量，分别为位置L、位置0处的电流相量；在放电发生前的脉冲源激励下，结合分布参数电路首末端的边界条件，即能够求解击穿放电发生前的分布参数电路上的电流、电压，获得放电发生前的波形。
[0011]进一步地，所述拟合分布参数电路击穿放电发生时刻的电压、电流随距离的变化曲线时，采用以幂函数为底的基函数拟合电流、电压在t
d
时刻随距离x的变化曲线，如下式所示：其中，y
v
、y
i
、、b...

【技术保护点】

【技术特征摘要】
1.一种分布参数电路放电电流的频域计算方法，其特征在于，包括：步骤一：对分布参数电路击穿放电发生前的波过程进行分析，得到放电发生前的波形；步骤二：拟合分布参数电路击穿放电发生时刻的电压、电流随距离的变化曲线；步骤三：结合步骤二所得变化曲线对分布参数电路击穿放电发生后的波过程进行分析，得到放电发生后的波形；步骤四：合并放电发生前的波形和放电发生后的波形，得到完整的分布参数电路放电的波过程。2. 根据权利要求1所述的一种分布参数电路放电电流的频域计算方法，其特征在于，所述步骤一具体为：设放电发生的时间为t
d
，脉冲源的波形为，则放电发生前频域中用于激励的脉冲源为时刻0至时刻t
d
的脉冲源波形，即，其中，表示击穿时刻前的激励源，t表示时间，表示信号的傅里叶变换，为阶跃函数；此时频域传输线方程为：其中，j为虚数单位，ω为角频率，x为位置变量，分别为电压相量、电流相量，Z,Y分别为传输线的阻抗矩阵、导纳矩阵；则此时分布参数电路首末端电压和电流关系即为：其中，为分布参数电路的传播常数，R、L
s
、G、C、j、ω分别为电阻、电感、电导、电容、虚数单位、角频率，L为分布参数电路首末端的距离，为分布参数电路的特征阻抗；分别为位置L、位置0处的电压相量，分别为位置L、位置0处的电流相量；
在放电发生前的脉冲源激励下，结合分布参数电路首末端的边界条件，即能够求解击穿放电发生前的分布参数电路上的电流、电压，获得放电发生前的波形。3.根据权利要求1所述的一种分布参数电路放电电流的频域计算方法，其特征在于，所述步骤二中采用以幂函数为底的基函数拟合电流、电压在t
d
时刻随距离x的变化曲线，如下式所示：其中，y
v
、y
i
、、b
i
、i、n分别为电压变量、电流变量、拟合电压随x变化用到的系数、拟合电流随x变化用到的系数、角标、拟合阶数，通过链参数矩阵在长度为L的分布参数电路取N个样本，即得：个样本，即得：个样本，即得：其中，Δx表示采样间隔，m为采样的样本量，k1表示第k1个样本编号；由此得到距离自变量对应的样本数据为，因变量电压对应的样本数据为，因变量电流对应的样本数据为，将样本数据代入以幂函数为底的基函数中，得到如下等式：
拟合电压随x变化用到的系数、拟合电流随x变化用到的系数b
i
即为：即为：令，，分别表示t
d
时刻电压随x的变化规律、t
d
时刻电流随x的变化规律，进而得分布参数电路击穿放电发生时刻的电压、电流随距离的变化曲线。4. 根据权利要求3所述的一种分布参数电路放电电流的频域计算方法，其特征在于，所述步骤三具体为：设放电发生的时间为t
d
，脉冲源的波形为，则放电发生后频域中用于激励的脉冲源为时刻t
d
之后的脉冲源波形，即，表示击穿时刻后的激励源，其中，表示信号的傅里叶变换，为阶跃函数；此时频域传输线方程为：
其中，L
S
为电感矩阵，C为电容矩阵；则首末端电压和电流关系为：则首末端电压和电流关系为：其中：其中：其中：、分别表示等效电压源、等效电流源，的自变量为距离x，分别表示分布电压源、分布电流源，其表达式为：，分别表示在t
d
时刻电流随x的变化关系、电压随x的变化关系；在放电发生后的脉冲源激励下，结合分布参数电路首末端的边界条件，即能够求解击穿放电发生后的分布参数电路上的电流、电压，获得放电发生后的波形。5.一种分布参数电路放电电流的频域计算系统，其特征在于，包括：放电发生前的波形获取模块：用于对分布参数电路击穿放电发生前的波过程进行分析，得到放电发...

【专利技术属性】
技术研发人员：曾卫东，王宗扬，谢彦召，田爽，高海东，肖勇，于在松，
申请(专利权)人：西安交通大学华能集团技术创新中心有限公司，
类型：发明
国别省市：