CN106771627B - 一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法 - Google Patents
一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法 Download PDFInfo
- Publication number
- CN106771627B CN106771627B CN201611031502.6A CN201611031502A CN106771627B CN 106771627 B CN106771627 B CN 106771627B CN 201611031502 A CN201611031502 A CN 201611031502A CN 106771627 B CN106771627 B CN 106771627B
- Authority
- CN
- China
- Prior art keywords
- equivalent
- phase
- impedance
- frequency
- pole
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/04—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant in circuits having distributed constants, e.g. having very long conductors or involving high frequencies
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法,测量多芯屏蔽动力电缆的共模阻抗与差模阻抗,建立共模阻抗与差模阻抗的等效电路模型,间接推导出多芯屏蔽动力电缆多导体传输线中频变单位长度电参数,分别反演出电阻、电感、电容与电导的极点‑留数形式的有理函数多项式逼近,此时系统为稳定系统,当该稳定系统与无源网络、有源网络连接时,分别处理,实现单位长度频变电参数宽频等效电路中元件的无源化。本发明不仅考虑其结构不对称性对多导体传输线模型中电参数频变特性的影响,而且有效隔离了周围电磁干扰对阻抗测量结果的影响,提高了其测量精度;该方法适用于结构简单的动力电缆,通信线缆以及分裂导线等模型的频变多导体传输线建模。
Description
技术领域
本发明涉及动力电缆宽频建模,具体涉及一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法。
背景技术
根据国际无线电干扰特别委员会CISPR.22标准,传导电磁干扰的频率范围为150kHz~30MHz。在该频率范围内,电机变频驱动系统中,连接变频器输入与输出的动力电缆的尺寸相对于传导电磁干扰最小波长一般属于“电大尺寸”。为了准确建立其模型,应采用多导体传输线分布参数模型。
动力电缆结构通常较复杂,呈现一定的不对称性;除动力线与地线外,往往还有温度控制线等,多根导体按照一定节距扭绞在一起;此外,控制线及动力线周围还有一定覆盖面积的编织屏蔽层。
由于集肤效应与邻近效应的共同影响,使得导体内部电流密度分布不均匀;频率越大,电流密度分布越不均匀,造成电阻与电感参数的频变特性越明显。
而且,电缆绝缘材料介电常数的频变特性导致其单位长度电容与电导也随频率变化。因此,为了建立准确的动力电缆多导体传输线模型,应考虑其单位长度电阻、电感、电容以及电导等分布参数的频变特性。
多导体传输线分布参数的获得主要有解析解法,数值法以及实验测量法等三种。解析解法能够精确处理结构简单的电缆模型,易考虑集肤效应对频变电阻的影响,但是电参数的频变特性很难考虑;数值方法便于计算结构较复杂的电缆,同时易考虑集肤效应与邻近效应,但对于具有编织屏蔽层及绞合缠绕的屏蔽电缆,其建模较复杂,且材料电参数的频变特性往往难以考虑;而实验测量法对传输线频变电参数的获取具有准确、快速、直接等优势,但测量结果也易受阻抗分析仪精度及周围电磁干扰的影响。
发明内容
本发明的目的在于提供一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法,该方法可以同时考虑动力电缆结构不对称性,集肤效应,邻近效应以及介质损耗等因素的影响。且频变电参数的等效电路宏模型是端口无源的,其中元件也是无源的,即严格无源的。
为实现上述目的,本发明采取以下技术方案实现。
一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法,包括以下步骤:第一步:根据共模阻抗与差模阻抗的定义,在电磁屏蔽室中,采用阻抗分析仪分别测量多芯屏蔽动力电缆不同接线方式时的共模阻抗与差模阻抗;
第二步:在第一步基础上,建立共模阻抗与差模阻抗的等效电路模型,间接推导出多芯屏蔽动力电缆多导体传输线中频变单位长度电参数,包括电阻、电感、电容与电导;
第三步:依据第二步得出的多芯屏蔽动力电缆多导体传输线中频变单位长度电参数,分别反演出电阻、电感、电容与电导的极点-留数形式的有理函数多项式逼近,该有理函数多项式包括常数项和比例项,实数极点项以及共轭复数极点对项;由于该有理函数多项式的极点的实数部分小于零,因此该多项式逼近系统为稳定系统;
当该稳定系统与无源网络连接时生成不稳定的新系统时,此时的系统端口是有源的,则进行步骤四;
当该稳定系统与无源网络连接时生成稳定的新系统时,此时的系统端口是无源的,则进行步骤五;
第四步:在第三步的基础上,使用留数摄动法,即保持极点不变,根据端口无源约束条件,修改相应留数值,以误差最小作为目标函数,从而获得端口无源的频变参数模型;
第五步:基于电路综合原理,分别获得单位长度频变电阻、电感、电容与电导的极点-留数形式的有理函数多项式逼近中的常数项和比例项,实数极点项以及共轭复数极点对项对应的等效电路,将所得3个等效电路并联之后得到单位长度频变电参数的宽频等效电路;
当单位长度频变电参数的宽频等效电路中元件有负值时,即是有源的,则进行步骤六;
当单位长度频变电参数的宽频等效电路中元件无负值时,即是无源的,实现了单位长度频变电参数宽频等效电路中元件的无源化;
第六步:以单位长度频变电参数的宽频等效电路中元件值大于0作为约束变量,采用优化算法进行优化,以优化后的多芯屏蔽动力电缆多导体传输线中频变单位长度电参数与测量得出的多芯屏蔽动力电缆多导体传输线中频变单位长度电参数之间的误差最小作为目标函数,从而实现单位长度频变电参数宽频等效电路中元件的无源化。
本发明进一步的改进在于,共模阻抗包括共模等效电阻、等效电感、等效电容以及等效电导,具体通过以下过程得到:
以含有A、B、C三相动力线导体与地线导体并采用编织屏蔽层铠装的动力线缆作为待测电缆,在共模阻抗测量时,在待测线缆首端,将三相动力线与地线短接后和阻抗分析仪一端连接,将编织屏蔽层与阻抗分析仪另一端连接;在待测线缆尾端,将三相动力线与地线短接后通过单刀单掷开关与编织屏蔽层连接;单刀单掷开关闭合时,测量出共模等效电阻与等效电感;单刀单掷开关打开时,测量出共模等效电容与等效电导。
本发明进一步的改进在于,差模阻抗测量时,为了考虑动力线缆结构不对称性对单位长度频变电参数的影响,设置三种不同的接线方式,具体如下:
第一种差模阻抗测量接线方式:在待测线缆首端,将A相与B相动力线短接后与阻抗分析仪相连,C相与地线短接后与阻抗分析仪另一端短接;在待测线缆尾端,将A相与B相动力线短接后通过单刀单掷开关和C相与地线短接后的电路连接;当单刀单掷开关闭合时,测量出差模阻抗等效电阻与等效电感;单刀单掷开关打开时,测量出差模阻抗等效电容与电导;
第二种差模阻抗测量接线方式:在待测线缆首端,将A相与C相动力线短接后与阻抗分析仪相连,B相与地线短接后与阻抗分析仪另一端短接;在待测线缆尾端,将A相与C相动力线短接后通过单刀单掷开关和B相与地线短接后的电路连接,分别测量出差模阻抗等效电阻、电感、电容以及电导;
第三种差模阻抗测量接线方式:在待测线缆首端,将A相与地线短接后与阻抗分析仪相连,B相与C相短接后与阻抗分析仪另一端短接;在待测线缆尾端,将A相与地线短接后通过单刀单掷开关和B相与C相短接后的电路连接,分别测量出差模阻抗等效电阻、电感、电容以及电导。
本发明进一步的改进在于,待测线缆的各相导体之间的互感耦合系数关系为
式中:K12——动力线导体L1与L2之间的互感耦合系数;K13——动力线导体L1与L3之间的互感耦合系数;K23——动力线导体L2与L3之间的互感耦合系数;K1G——动力线导体L1与地线G之间的互感耦合系数;K2G——动力线导体L2与地线G之间的互感耦合系数;K3G——动力线导体L3与地线G之间的互感耦合系数。
本发明进一步的改进在于,分别反演出电阻、电感、电容与电导的极点-留数形式的有理函数多项式逼近通过运用矢量匹配法实现。
本发明进一步的改进在于,电阻、电感、电容与电导的极点-留数形式的有理函数多项式逼近如式(2)
式中:Yrat(s)——Y(s)的有理函数拟合表达式;P——实数极点个数;s——复频域变量;Q——共轭复数极点对个数;ap——实数极点;cp——实数留数;aq、aq *——一对共轭复数极点;cq、cq *——一对共轭复数留数;d——常数项;h——比例项。
本发明进一步的改进在于,优化算法为序列二次规划法与模拟退火算法。
与现有技术相比,本发明具有显著的有益效果为:由于本发明采用阻抗分析仪在电磁屏蔽室中测量多芯屏蔽动力电缆阻抗,不仅考虑其结构不对称性对多导体传输线模型中电参数频变特性的影响,而且有效隔离了周围电磁干扰对阻抗测量结果的影响,提高了其测量精度;此外,采用阻抗分析仪的扫频功能兼顾了动力电缆集肤效应、邻近效应以及介电损耗等因素对电参数频变特性的影响;运用电路综合原理可以精确获得这些频变电参数的等效电路模型,且便于分析单位长度动力电缆分段数对仿真结果精度的影响,适合时域电路仿真分析;同时,频变参数的等效电路宏模型是端口无源的,且其中元件值也均是无源的,保证了与其他无源网络级联时整个电路系统的稳定性;该方法同样适用于结构简单的动力电缆,通信线缆以及分裂导线等模型的频变多导体传输线建模。
进一步的,基于阻抗测量的方法可以计及动力电缆编织屏蔽层及其扭绞节距对电参数频变特性的影响。
进一步的,采用矢量匹配法既可以逼近四个电参数的模值,同时也可以近似电参数的相位,进而反演出四个电参数极点-留数形式的多项式逼近,拟合度大大提高。
附图说明
图1是单根导体频变传输线模型;
图2是实施例多芯屏蔽动力电缆结构;
图3是实施例动力电缆多导体传输线模型;
图4是实施例动力电缆共模阻抗测量接线方式;
图5是实施例动力电缆共模阻抗测量等效电路;
图6是实施例动力电缆差模阻抗测量三种接线方式;其中,(a)接线方式1:L1与L2短接,L3与G短接,(b)接线方式2:L1与L3短接,L2与G短接,(c)接线方式3:L1与G短接,L2与L3短接;
图7是实施例动力电缆差模阻抗测量三种接线方式的等效电路;其中,(a)接线方式1:L1与L2短接,L3与G短接,(b)接线方式2:L1与L3短接,L2与G短接,(c)接线方式3:L1与G短接,L2与L3短接;
图8是频变电参数多项式近似后的常数项与比例项,实数极点以及共轭复数极点对等三项的等效电路;其中,(a)常数项与比例项,(b)实数极点项,(c)共轭复数极点对项;
图9是频变电参数的等效电路。
图中,1为外护套,2为总屏蔽,3为总扎捆带,4为羊毛填充物,5为单元屏蔽层,6为屏蔽衬底,7为单元捆扎带,8为内护套,9为芯线绝缘层,10为动力线导体,11为温度测量线。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法,包括以下步骤:
第一步:根据共模阻抗与差模阻抗的定义,在电磁屏蔽室中,采用阻抗分析仪分别测量多芯屏蔽动力电缆不同接线方式时的共模阻抗与差模阻抗。共模阻抗包括共模等效电阻、等效电感、等效电容以及等效电导;
以含有A、B、C三相动力线导体与地线导体并采用编织屏蔽层铠装的动力线缆作为待测电缆,在共模阻抗测量时,在待测线缆首端,将三相动力线与地线短接后和阻抗分析仪一端连接,将编织屏蔽层与阻抗分析仪另一端连接;在待测线缆尾端,将三相动力线与地线短接后通过单刀单掷开关与编织屏蔽层连接;单刀单掷开关闭合时,测量出共模等效电阻与等效电感;单刀单掷开关打开时,测量出共模等效电容与等效电导。
差模阻抗测量时,为了考虑动力线缆结构不对称性对单位长度频变电参数的影响,设置三种不同的接线方式,具体如下:
第一种差模阻抗测量接线方式:在待测线缆首端,将A相与B相动力线短接后与阻抗分析仪相连,C相与地线短接后与阻抗分析仪另一端短接;在待测线缆尾端,将A相与B相动力线短接后通过单刀单掷开关和C相与地线短接后的电路连接;当单刀单掷开关闭合时,测量出差模阻抗等效电阻与等效电感;单刀单掷开关打开时,测量出差模阻抗等效电容与电导;
第二种差模阻抗测量接线方式:在待测线缆首端,将A相与C相动力线短接后与阻抗分析仪相连,B相与地线短接后与阻抗分析仪另一端短接;在待测线缆尾端,将A相与C相动力线短接后通过单刀单掷开关和B相与地线短接后的电路连接,分别测量出差模阻抗等效电阻、电感、电容以及电导;
第三种差模阻抗测量接线方式:在待测线缆首端,将A相与地线短接后与阻抗分析仪相连,B相与C相短接后与阻抗分析仪另一端短接;在待测线缆尾端,将A相与地线短接后通过单刀单掷开关和B相与C相短接后的电路连接,分别测量出差模阻抗等效电阻、电感、电容以及电导。
待测线缆的各相导体之间的互感耦合系数关系为
式中:K12——动力线导体L1与L2之间的互感耦合系数;K13——动力线导体L1与L3之间的互感耦合系数;K23——动力线导体L2与L3之间的互感耦合系数;K1G——动力线导体L1与地线G之间的互感耦合系数;K2G——动力线导体L2与地线G之间的互感耦合系数;K3G——动力线导体L3与地线G之间的互感耦合系数。
第二步:在第一步基础上,建立共模阻抗与差模阻抗的等效电路模型,间接推导出多芯屏蔽动力电缆多导体传输线中频变单位长度电参数,包括电阻、电感、电容与电导;
第三步:依据第二步得出的多芯屏蔽动力电缆多导体传输线中频变单位长度电参数,通过矢量匹配法分别反演出电阻、电感、电容与电导的极点-留数形式的有理函数多项式逼近,电阻、电感、电容与电导的极点-留数形式的有理函数多项式逼近如式(2)
式中:Yrat(s)——Y(s)的有理函数拟合表达式;P——实数极点个数;s——复频域变量;Q——共轭复数极点对个数;ap——实数极点;cp——实数留数;aq、aq *——一对共轭复数极点;cq、cq *——一对共轭复数留数;d——常数项;h——比例项。
该有理函数多项式包括常数项和比例项,实数极点项以及共轭复数极点对项;由于该有理函数多项式的极点的实数部分小于零,因此该多项式逼近系统为稳定系统;
当该稳定系统与无源网络连接时生成不稳定的新系统时,此时的系统端口是有源的,则进行步骤四;
当该稳定系统与无源网络连接时生成稳定的新系统时,此时的系统端口是无源的,则进行步骤五;
第四步:在第三步的基础上,使用留数摄动法,即保持极点不变,根据端口无源约束条件,修改相应留数值,以误差最小作为目标函数,从而获得端口无源的频变参数模型;
第五步:基于电路综合原理,分别获得单位长度频变电阻、电感、电容与电导的极点-留数形式的有理函数多项式逼近中的常数项和比例项,实数极点项以及共轭复数极点对项对应的等效电路,将所得3个等效电路并联之后得到单位长度频变电参数的宽频等效电路;
当单位长度频变电参数的宽频等效电路中元件有负值时,即是有源的,则进行步骤六;
当单位长度频变电参数的宽频等效电路中元件无负值时,即是无源的,实现了单位长度频变电参数宽频等效电路中元件的无源化;
第六步:以单位长度频变电参数的宽频等效电路中元件值大于0作为约束变量,采用优化算法为序列二次规划法与模拟退火算法进行优化,以优化后的多芯屏蔽动力电缆多导体传输线中频变单位长度电参数与测量得出的多芯屏蔽动力电缆多导体传输线中频变单位长度电参数之间的误差最小作为目标函数,从而实现单位长度频变电参数宽频等效电路中元件的无源化。
下面通过一个实施例进行详细说明。
本发明采用图1中的L型电路来建立考虑动力电缆单位长度电参数频变的单根导体传输线模型。该电路由电阻,电感,电容与电导的并联共同串接而成;电容与电导并联部分的分压作为输出端口。其中,R0(f),L0(f),C0(f),与G0(f)为动力电缆单位长度1m的电阻(Ω/m),电感(H/m),电容(F/m)与电导(S/m);n为1m长传输线模型的子单元个数。
该实施例屏蔽电缆由八根绞线构成,如图2所示。屏蔽电缆包括外护套1,外护套1内设置有总屏蔽2,总屏蔽2与外护套1之间设置有总扎捆带3,总屏蔽2内设置有内护套8,外护套1内中心设置有羊毛填充物4,羊毛填充物4的上方为一对温度测量线11,羊毛填充物4的下方为一对温度测量线11,羊毛填充物4的一侧为两根动力线,另一侧为一根动力线和一根地线,两个温度测量线、三根动力线(三根动力线为三相动力线)以及地线围绕在羊毛填充物的周围。两根温度测量线的结构相同,均包括单元屏蔽层5,单元屏蔽层5内设置有屏蔽衬底6,屏蔽衬底6内部设置有单元捆扎带7;每根动力线包括芯线绝缘层9以及设置在芯线绝缘层内的动力线导体10。
其中,L1,L2与L3为三相动力线,标有字母G的为地线,羊毛填充物4下方和上方的四根导线为两对温度测量线。外护套1由强力挤压无卤素PUR制成;总屏蔽2和温度测量线的单元屏蔽层5由强抗扭力、镀铜锡双层编织层构成,平均线覆盖率为70%,平均光学覆盖率90%;动力线导体10、芯线绝缘层9和温度测量线绝缘层由高机械性能的低电容TPE混合物构成;动力线导体10与温度测量线的导体由高柔性耐弯折多股细绞精细裸铜丝构成;动力线的导体和温度测量线围绕中间加强芯分束以短节距绞合而成;内护套为甲胄式强力挤压成型PUR构成;总捆扎带3与单元捆扎带7由高强度铝复合塑料构成;屏蔽衬底6由塑料薄膜构成。
考虑动力电缆结构的不对称性,动力电缆的多导体传输线模型如图3所示。其中R0表示动力线导体单位长度电阻/Ω·m-1;L0表示动力线导体单位长度电感/H·m-1;C0S与G0S分别为动力线导体与电缆总屏蔽层之间的单位长度电容/F·m-1与单位长度电导/S·m-1。动力电缆总屏蔽与温度控制线编织屏蔽层S双端接地,视作理想导体。根据动力电缆结构不对称性,可知各相导体之间的互感耦合系数关系为
式中:K12——动力线导体L1与L2之间的互感耦合系数;K13——动力线导体L1与L3之间的互感耦合系数;K23——动力线导体L2与L3之间的互感耦合系数;K1G——动力线导体L1与地线G之间的互感耦合系数;K2G——动力线导体L2与地线G之间的互感耦合系数;K3G——动力线导体L3与地线G之间的互感耦合系数。
首先进行动力电缆共模阻抗测量。图4中,在首端,将L1、L2、L3与G短接并与阻抗分析仪一端相接,S与阻抗分析仪另一端相接;在末端,将L1、L2、L3与G短接后通过单刀单掷开关(Single pole single throw,SPST)与S相连。
图5是动力电缆共模阻抗测量的等效电路。其中,RCM,LCM,CCM与GCM分别为其等效电阻、电感、电容以及电导。
共模阻抗测量原理:当SPST闭合时,在被测动力电缆末端将动力线L1、L2、L3与G短接后与S相连,由此可以测量出图4中的等效电阻RCM与等效电感LCM;当SPST断开时,可以测量出等效电容CCM与等效电导GCM。
为了考虑动力电缆结构的不对称性对单位长度电感及其耦合系数以及单位长度电容的影响,根据差模阻抗定义,设置三组差模测量接线方式,分别如图6(a)、6(b)与6(c)所示。图6(a)中,在动力电缆首端,将L1与L2短接,L3与G短接并分别与阻抗分析仪两端相接,在其末端,将L1与L2短接及L3与G短接后的电路通过SPST连接;图6(b)中,在其首端,将L1与L3短接,L2与G短接并分别与阻抗分析仪两端相接,在末端,将L1与L3短接及L2与G短接后的电路通过SPST连接;图6(c)中,在其首端,将L1与G短接,L2与L3短接并分别与阻抗分析仪两端相接,在其末端,将L1与G短接及L2与L3短接后的电路通过SPST连接。
图6(a)、6(b)与6(c)中动力电缆差模阻抗测量三种不同接线方式的等效电路分别如图7(a),7(b)与7(c)所示。其中,RDM1,RDM2,RDM3,LDM1,LDM2,LDM3,CDM1,CDM2,CDM3,以及GDM1,GDM2,GDM3分别为差模阻抗测量时三种不同接线方式下的等效电阻、电感、电容以及电导。
差模阻抗测量原理与共模阻抗测量原理相似:以图6(a)测量接线方式1为例:当SPST闭合时,可以测量出等效电感LDM1,LDM2,LDM3;当SPST断开时,可以测量出等效电容CDM1,CDM2,CDM3与等效电导GDM1,GDM2,GDM3。
根据实施例动力电缆共模阻抗与三种差模阻抗测量接线方式所测参数及所测参数的等效电路可以反推出图3中多导体传输线中各频变电参数。
频变电参数Y(s)的极点-留数形式的有理函数多项式逼近如式(2)
式中:Yrat(s)——Y(s)的有理函数拟合表达式;P——实数极点个数;s——复频域变量;Q——共轭复数极点对个数;ap——实数极点;cp——实数留数;aq、aq *——一对共轭复数极点;cq、cq *——一对共轭复数留数;d——常数项;h——比例项。
式(2)中的未知参数P、Q、ap、cp、aq、aq *、cq、cq *、d及h可采用矢量匹配法来确定。通常,拟合阶数越高,Yrat(s)与Y(s)的拟合度越高,精度越高。此时的系统为稳定系统,但是,当其与无源网络连接时可能生成不稳定的新系统,即此时的系统端口可能是有源的;
因此,本发明使用留数摄动法,即保持极点不变,根据宏模型端口无源约束条件,修改相应留数值,以获得端口无源的频变参数宏模型;
当式(2)中的未知参数确定后,可采用电路综合原理求网络函数Yrat(s)各部分的等效电路。
常数项与比例项的等效电路:式(2)中常数项与比例项可采用图8(a)中的R0与C0并联电路来等效。图8(a)中电路的导纳函数为
由式(2)中的常数项与比例项及式(3)可求得图8(a)中电路参数
R0=1/d,C0=h(4)
实数极点项的等效电路:式(2)中实数极点与实数留数可采用图8(b)中的R1与L1串联电路来等效。图8(b)中电路的导纳函数为
由式(2)中的实数极点与实数留数项及式(5)可知图8(b)中电路参数
共轭复数极点对项的等效电路:式(2)中一对共轭复数极点与留数可采用图8(c)中的R1s、Ls、Cs及R2s的串并联电路来等效,图8(c)中电路的导纳函数为
由式(2)中的一对共轭复数极点与留数项及式(7)可知图8(c)中电路参数
至此,图3中多芯屏蔽动力电缆中的频变电参数的等效电路反演结束,公式(2)中Yrat(s)的等效电路整体结构如图9所示。此时,等效电路中可能有些元件值是负值,即元件是有源的;
最后,先后采用序列二次规划法与模拟退火算法,以优化后的频变曲线与测量得出的频变曲线之间的误差最小作为目标函数,等效电路中元件值作为约束变量,从而实现单位长度频变电参数宽频等效电路中元件的无源化。
Claims (7)
1.一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法,其特征在于,包括以下步骤:第一步:根据共模阻抗与差模阻抗的定义,在电磁屏蔽室中,采用阻抗分析仪分别测量多芯屏蔽动力电缆不同接线方式时的共模阻抗与差模阻抗;
第二步:在第一步基础上,建立共模阻抗与差模阻抗的等效电路模型,间接推导出多芯屏蔽动力电缆多导体传输线中频变单位长度电参数,包括电阻、电感、电容与电导;
第三步:依据第二步得出的多芯屏蔽动力电缆多导体传输线中频变单位长度电参数,分别反演出电阻、电感、电容与电导的极点-留数形式的有理函数多项式逼近,该有理函数多项式包括常数项和比例项,实数极点项以及共轭复数极点对项;由于该有理函数多项式的极点的实数部分小于零,因此该多项式逼近系统为稳定系统;
当该稳定系统与无源网络连接时生成不稳定的新系统时,此时的系统端口是有源的,则进行步骤四;
当该稳定系统与无源网络连接时生成稳定的新系统时,此时的系统端口是无源的,则进行步骤五;
第四步:在第三步的基础上,使用留数摄动法,即保持极点不变,根据端口无源约束条件,修改相应留数值,以误差最小作为目标函数,从而获得端口无源的频变参数模型;
第五步:基于电路综合原理,分别获得单位长度频变电阻、电感、电容与电导的极点-留数形式的有理函数多项式逼近中的常数项和比例项,实数极点项以及共轭复数极点对项对应的等效电路,将所得3个等效电路并联之后得到单位长度频变电参数的宽频等效电路;
当单位长度频变电参数的宽频等效电路中元件有负值时,即是有源的,则进行步骤六;
当单位长度频变电参数的宽频等效电路中元件无负值时,即是无源的,实现了单位长度频变电参数宽频等效电路中元件的无源化;
第六步:以单位长度频变电参数的宽频等效电路中元件值大于0作为约束变量,采用优化算法进行优化,以优化后的多芯屏蔽动力电缆多导体传输线中频变单位长度电参数与测量得出的多芯屏蔽动力电缆多导体传输线中频变单位长度电参数之间的误差最小作为目标函数,从而实现单位长度频变电参数宽频等效电路中元件的无源化。
2.根据权利要求1所述的一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法,其特征在于,共模阻抗包括共模等效电阻、等效电感、等效电容以及等效电导,具体通过以下过程得到:
以含有A、B、C三相动力线导体与地线导体并采用编织屏蔽层铠装的动力线缆作为待测电缆,在共模阻抗测量时,在待测线缆首端,将三相动力线与地线短接后和阻抗分析仪一端连接,将编织屏蔽层与阻抗分析仪另一端连接;在待测线缆尾端,将三相动力线与地线短接后通过单刀单掷开关与编织屏蔽层连接;单刀单掷开关闭合时,测量出共模等效电阻与等效电感;单刀单掷开关打开时,测量出共模等效电容与等效电导。
3.根据权利要求2所述的一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法,其特征在于,差模阻抗测量时,为了考虑动力线缆结构不对称性对单位长度频变电参数的影响,设置三种不同的接线方式,具体如下:
第一种差模阻抗测量接线方式:在待测线缆首端,将A相与B相动力线短接后与阻抗分析仪相连,C相与地线短接后与阻抗分析仪另一端短接;在待测线缆尾端,将A相与B相动力线短接后通过单刀单掷开关和C相与地线短接后的电路连接;当单刀单掷开关闭合时,测量出差模阻抗等效电阻与等效电感;单刀单掷开关打开时,测量出差模阻抗等效电容与等效电导;
第二种差模阻抗测量接线方式:在待测线缆首端,将A相与C相动力线短接后与阻抗分析仪相连,B相与地线短接后与阻抗分析仪另一端短接;在待测线缆尾端,将A相与C相动力线短接后通过单刀单掷开关和B相与地线短接后的电路连接,分别测量出差模阻抗等效电阻、等效电感、等效电容以及等效电导;
第三种差模阻抗测量接线方式:在待测线缆首端,将A相与地线短接后与阻抗分析仪相连,B相与C相短接后与阻抗分析仪另一端短接;在待测线缆尾端,将A相与地线短接后通过单刀单掷开关和B相与C相短接后的电路连接,分别测量出差模阻抗等效电阻、等效电感、等效电容以及等效电导。
4.根据权利要求2所述的一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法,其特征在于,待测线缆的各相导体之间的互感耦合系数关系为
式中:K12——动力线导体L1与L2之间的互感耦合系数;K13——动力线导体L1与L3之间的互感耦合系数;K23——动力线导体L2与L3之间的互感耦合系数;K1G——动力线导体L1与地线G之间的互感耦合系数;K2G——动力线导体L2与地线G之间的互感耦合系数;K3G——动力线导体L3与地线G之间的互感耦合系数。
5.根据权利要求1所述的一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法,其特征在于,分别反演出电阻、电感、电容与电导的极点-留数形式的有理函数多项式逼近通过运用矢量匹配法实现。
6.根据权利要求1所述的一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法,其特征在于,电阻、电感、电容与电导的极点-留数形式的有理函数多项式逼近如式(2)
式中:Yrat(s)——Y(s)的有理函数拟合表达式;P——实数极点个数;s——复频域变量;Q——共轭复数极点对个数;ap——实数极点;cp——实数留数;aq、aq *——一对共轭复数极点;cq、cq *——一对共轭复数留数;d——常数项;h——比例项。
7.根据权利要求1所述的一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法,其特征在于,优化算法为序列二次规划法与模拟退火算法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611031502.6A CN106771627B (zh) | 2016-11-18 | 2016-11-18 | 一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611031502.6A CN106771627B (zh) | 2016-11-18 | 2016-11-18 | 一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106771627A CN106771627A (zh) | 2017-05-31 |
CN106771627B true CN106771627B (zh) | 2019-04-30 |
Family
ID=58970853
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611031502.6A Active CN106771627B (zh) | 2016-11-18 | 2016-11-18 | 一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106771627B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107895078A (zh) * | 2017-11-09 | 2018-04-10 | 西安交通大学 | 强电磁脉冲对长电力输电线路耦合响应的宏模型构建方法 |
CN108763802A (zh) * | 2018-06-04 | 2018-11-06 | 天津农学院 | 计及电力传输线的电力电子换流系统电磁干扰估算方法 |
CN109932567B (zh) * | 2019-03-27 | 2021-02-02 | 国网山东省电力公司青岛供电公司 | 架空输电线路参数的无源测量方法 |
CN109977557B (zh) * | 2019-03-28 | 2023-05-05 | 广东志成冠军集团有限公司 | 一种适用于电声纵振换能器等效电路的建模方法 |
CN110188381B (zh) * | 2019-04-18 | 2023-04-18 | 中国北方车辆研究所 | 一种用于电磁干扰预测的仿真模型的构建方法及系统 |
CN112904252A (zh) * | 2019-12-04 | 2021-06-04 | 清华大学 | 一种分析光学电场传感器的频率响应的方法 |
CN111444592B (zh) * | 2020-02-20 | 2024-02-23 | 重庆大学 | 一种变压器宽频导纳模型及建立方法 |
CN112083237B (zh) * | 2020-07-31 | 2021-06-29 | 西安交通大学 | 一种用于大尺度电气设备宽频特性时域测量方法及系统 |
CN113191036B (zh) * | 2021-03-19 | 2023-03-14 | 西安电子科技大学 | 一种共面波导等效电路结构及其参数提取方法 |
CN112906335B (zh) * | 2021-03-22 | 2022-04-15 | 北京华大九天科技股份有限公司 | 一种用于集成电路系统的无源性校正方法及装置 |
CN117034496B (zh) * | 2023-08-22 | 2024-05-28 | 山东兖矿集团长龙电缆制造有限公司 | 一种基于Grasshopper的电缆参数化建模方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002197136A (ja) * | 2000-10-30 | 2002-07-12 | Hewlett Packard Co <Hp> | 物理回路の電気的振舞いをシミュレートする方法 |
CN102693364A (zh) * | 2012-05-17 | 2012-09-26 | 西安电子科技大学 | 广义通信/控制总线端接模型等效转换方法 |
CN103488839A (zh) * | 2013-09-27 | 2014-01-01 | 中国东方电气集团有限公司 | 一种电动汽车电源线瞬态传导发射的建模方法 |
CN106055837A (zh) * | 2016-06-29 | 2016-10-26 | 中国人民解放军军械工程学院 | 外场激励下有损大地上架空线缆电路模型建立方法和系统 |
CN106096176A (zh) * | 2016-06-23 | 2016-11-09 | 北京航空航天大学 | 屏蔽线缆的软辫线的天线模型的建模方法及装置 |
-
2016
- 2016-11-18 CN CN201611031502.6A patent/CN106771627B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002197136A (ja) * | 2000-10-30 | 2002-07-12 | Hewlett Packard Co <Hp> | 物理回路の電気的振舞いをシミュレートする方法 |
CN102693364A (zh) * | 2012-05-17 | 2012-09-26 | 西安电子科技大学 | 广义通信/控制总线端接模型等效转换方法 |
CN103488839A (zh) * | 2013-09-27 | 2014-01-01 | 中国东方电气集团有限公司 | 一种电动汽车电源线瞬态传导发射的建模方法 |
CN106096176A (zh) * | 2016-06-23 | 2016-11-09 | 北京航空航天大学 | 屏蔽线缆的软辫线的天线模型的建模方法及装置 |
CN106055837A (zh) * | 2016-06-29 | 2016-10-26 | 中国人民解放军军械工程学院 | 外场激励下有损大地上架空线缆电路模型建立方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106771627A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106771627B (zh) | 一种建立多芯屏蔽动力电缆频变多导体传输线模型的方法 | |
CN105243197B (zh) | 基于有限元仿真和等效电路的单芯电力电缆谐振分析方法 | |
Papadopoulos et al. | Narrowband power line communication: Medium voltage cable modeling and laboratory experimental results | |
Pagnetti et al. | An improved method for the calculation of the internal impedances of solid and hollow conductors with the inclusion of proximity effect | |
Gustavsen et al. | A finite-element approach for calculating electrical parameters of umbilical cables | |
Gudmundsdottir | Proximity effect in fast transient simulations of an underground transmission cable | |
Li et al. | Measurement-based modeling and worst-case estimation of crosstalk inside an aircraft cable connector | |
Kruizinga et al. | High frequency modeling of a shielded four-core low voltage underground power cable | |
CN103745054B (zh) | 一种对电缆及电缆束在电磁兼容中的建模及信号串扰分析方法 | |
Wang et al. | Estimation of multi‐conductor powerline cable parameters for the modelling of transfer characteristics | |
CN106777444B (zh) | 一种mvb总线网络传输介质模型设计方法 | |
CN109254189B (zh) | 一种海底电缆耐压试验终端电压计算方法 | |
CN116859182A (zh) | 考虑电缆接头的频域反射谱高压电缆缺陷定位方法及系统 | |
Li et al. | Power cable joint model: based on lumped components and cascaded transmission line approach | |
Guo et al. | Channel model for low voltage three-core power line communication | |
CN110472366B (zh) | 一种载波信号在中压配电网地埋电缆线路传输的解耦方法 | |
Badzagua et al. | Effective computational techniques for EMC analysis of cable harness | |
del-Pino-López et al. | A 3D parametric analysis of three-core armored power cables series impedance | |
Rotgerink et al. | Multi-conductor transmission line modelling of transfer impedance measurement methods | |
Arentsen et al. | MVAC submarine cable, impedance measurements and analysis | |
Buccella et al. | Three-dimensional FEM approach to model twisted wire pair cables | |
Buccella et al. | Prediction of voltage and current propagation in twisted wire pairs (TWPs) by a circuit model | |
Wang et al. | Transfer characteristics modeling of four-conductor cables in power-line communications | |
Hassen et al. | A stranded Unshielded Twisted Pair modeling for online fault location using OMTDR-based diagnosis sensor | |
Ridel et al. | Numerical simulation of aeronautic cable topology and experimental validation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |