本发明专利技术公开了一种获取圆柱形导体电磁暂态过程中时域电流响应的方法和系统，以解决圆柱形导体电磁暂态过程中现有的计算方法难以实用化的问题。所述方法包括：利用电压信号的阶跃函数级数逼近方法，将施加在导体两端的电压信号按照阶跃函数级数进行时域上的分解，将其表示成一系列阶跃电压的叠加形式；根据导体对单位阶跃电压信号的时域电流响应，再将分解后的每个阶跃电压所对应的时域电流响应进行组合；从而获得圆柱形导体对所施加电压信号的时域电流响应。本发明专利技术不仅对施加在圆柱形导体或等效圆柱形导体两端的任意形式的电压信号都能进行处理，而且对电压信号进行了离散化处理，可用计算机进行数值运算，提高了计算精度和速度。
【技术实现步骤摘要】
本专利技术涉及导体的电磁场暂态计算领域，更具体地涉及一种获取圆柱形导体电磁暂态过程中时域电流响应的方法和系统。
技术介绍
目前在圆柱形导体电磁场暂态过程的计算中，一般所采用的方法主要可分为时域计算和频域计算两种。以上方法的核心都是求解贝塞尔方程。对于截面为非圆形的柱形导体，也可以将非圆形截面等效为等面积的圆形后再进行计算。圆柱形导体的时域计算方法由H. S. Carslaw与J. C. Jaeger为解决导体热传导问题最先提出，并由Edward J. Tuohy等人引入电学领域，但该方法仅限定为导体表面电动势恒定时电流密度的扩散过程，并未建立导体两端电压与电流之间的关系。圆柱形导体的频域计算方法由R. J.Hill提出，该方法利用圆柱形导体电流信号频率稳定时导体的电阻与电感值不变这一特性，将贝塞尔方程转换至频域，将电路整体进行频域计算，最后将结果采用傅立叶反变换法转至时域，获得其时域解。该方法计算较为准确，但由于在实际工程利用时，大量元件无法转换至频域表示，如整流器、电动机等，所以该方法难以在实际工程计算中实现。综上所述，目前的时域计算方法只能计算导体表面电动势恒定时电流密度的扩散过程，但在实际应用中因为导体表面的电动势通常是不确定的，所以限制了该方法的应用范围；而频域计算方法在实际的工程计算中，由于大量元器件无法转换至频域表示而不能进行计算，所以能求解的问题也相当有限。而随着高压直流输电及轨道交通领域的快速发展，对导体电磁场暂态过程的计算在实用性上提出了更高的要求，特别是在电力系统的电磁暂态仿真、继电保护算法研究与定值整定等方面有着重要的意义。总之，需要本领域技术人员迫切解决的一个技术问题是如何能够实现圆柱形导体电磁暂态过程计算方法的实用化。
技术实现思路
本专利技术所要解决的技术问题是如何能够实现圆柱形导体电磁暂态过程计算方法在实际工程应用中的实用化。为了解决上述问题，本专利技术公开了一种获取圆柱形导体在电磁暂态过程中时域电流响应的方法，包括对施加在圆柱形导体或等效圆柱形导体两端的电压信号在采样测量前首先进行滤波处理及转换处理；对电压信号进行采样测量，并用每一采样时刻得到的电压值减去上一采样时刻的电压值，得到当前时刻的电压变化值，所述的电压变化值也是对电压信号按照阶跃函数级数逼近方法进行分解后的某一个分量电压值；3按照前述方法得到的每个分量电压值都对应一个时域电流响应，该时域电流响应为单位阶跃电压信号所产生的电流响应与该分量电压值的乘积；将所有分量电压值对应的时域电流响应进行时域内的组合，就能得到在电磁暂态过程中任意电压信号下圆柱形导体或等效圆柱形导体中的时域电流响应。所述的电压信号在采样前的滤波处理及转换处理是指高次谐波滤除及模数转换处理。所述方法还包括，获取目标圆柱形导体的详细参数，所述参数包括电导率、相对磁导率、半径和导体长度；若目标导体为非圆柱形，则须为可等效为圆柱形导体的导体；等效圆柱形导体是指横截面为非圆形的均勻柱形导体，所述横截面为非圆形的均勻柱形导体可等效为等面积的圆形横截面的均勻柱形导体进行计算，并按照等效后的圆柱形导体获取其详细参数。所述方法还包括，利用圆柱形导体在阶跃电压下的集肤效应和电流密度扩散公式，获得单位阶跃电压信号施加在单位长度圆柱形导体上的时域电流响应。本专利技术还公开了一种获取圆柱形导体在电磁暂态过程中时域电流响应的系统，所述系统包括电压信号处理模块，用于对施加在圆柱形导体或等效圆柱形导体两端的电压信号在采样测量前进行滤波处理及转换；所述电压信号处理模块包括高次谐波滤除装置及模数转换器件。电压信号采集模块，用于采集和记录不同采样时刻的电压值；所述电压信号采集模块包括电压信号的采集设备及记录采样数值的存储设备。阶跃电压级数模块，用于计算和记录每一采样时刻相对于上一采样时刻的电压变化值，所述的电压变化值是指电压信号按照阶跃函数级数逼近方法进行分解后的某一个分量电压值；所述阶跃电压级数模块包括能计算相邻采样时刻的电压变化值的计算设备及记录分量电压值的存储设备。时域电流组合模块，用于计算各个分量电压值对应的时域电流响应，并将得到的时域电流响应进行组合；所述时域电流组合模块包括能计算和组合每个分量电压值对应的时域电流响应的处理设备。与现有技术相比，本专利技术具有以下优点首先，本专利技术采用阶跃函数级数逼近电压信号的方法，将问题转化成求解单位阶跃电压信号施加在圆柱形导体两端时导体中的时域电流响应，而将电压信号表示成阶跃函数级数的形式后，级数中的每一项阶跃电压所对应的时域电流响应都能进行求解，然后将所有这些时域电流响应进行时域内的组合，就能得到圆柱形导体电磁暂态过程的时域电流响应。该方法将复杂的电磁暂态计算过程分成了几个易于实现的步骤，提高了使用上的灵活性。其次，由于任意形式的电压信号都能利用阶跃函数级数逼近方法进行近似，近似的精度依赖于采样时刻的间隔，采样越密集，精度越高。因此，施加在圆柱形导体上的任意形式的电压信号都能利用本方法进行计算，这在直流馈电线路、地铁线路等的电磁场暂态计算中具有较强的实用性。再次，电压信号利用阶跃函数级数进行逼近后相当于对电压信号进行了离散化处4理，后续的计算都能利用计算机进行数值运算，大大提高了计算的速度和准确度。附图说明图1为本专利技术实施例所述获取圆柱形导体电磁暂态过程的时域电流响应的方法流程图；图2为本专利技术实施例中对圆柱形导体施加的单位阶跃电压信号及其产生的电流响应的示意图；图3为本专利技术实施例中对电压信号采用阶跃函数级数逼近方法进行表示的示意图；图4为本专利技术实施例中对电压采用阶跃函数级数表示后所对应的圆柱形导体电流响应的示意图；图5为本专利技术实施例所述的系统所包含的单元模块结构示意图。具体实施例方式为使本专利技术的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本专利技术作进一步详细的说明。本专利技术可以用多种方法实现，可以是具有专一功能的专用设备，也可以和其它功能(如MATLAB、RTDS)共用仿真计算平台。参照图1，本专利技术实施例中所述的获取圆柱形导体电磁暂态过程的时域电流响应的方法流程图。步骤101，首先获取目标圆柱形导体的详细参数；所述目标圆柱形导体的详细参数，包括电导率σ、相对磁导率μ、半径R、导体长度L等；若目标导体为非圆柱形，则须为可等效为圆柱形导体的导体，如钢轨等，将非圆柱形导体截面等效为等面积的圆形，再获得上述参数，以进行计算。步骤102，获得单位阶跃电压信号施加在单位长度圆柱形导体上时所产生的时域电流响应；参照图2，单位阶跃电压信号施加在单位圆柱形导体上时，计算相应的时域电流响应的具体方法如下在0时刻，电压信号发生突变，这时电压信号的频率可被认为是无穷大。根据集肤效应的定义，可认为电流基本集中于圆柱形导体的表层，在圆柱形导体内部分布无限少。在计算过程中可认为圆柱形导体内部电流近似等于0 ；圆导体表面电流密度Js可以由欧姆定律Js = Es/ο =U/ο计算获得。其他时亥lj，由于导体内部电流密度J(r，t)可以被视为符合Maxwell方程且表面电流密度不变的扩散过程。这一过程，Edward J. Tuohy等人引入R. J. Hill提出的贝塞尔方程级数解法，其电流密度J (r，t)关于时间t与扩散半径r的解为J(r,t) = Js\本文档来自技高网...
【技术保护点】
【技术特征摘要】
【专利技术属性】
技术研发人员：和敬涵，喻乐，胡俊，徐金凤，薄志谦，
申请(专利权)人：北京交通大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。