CN114563607A - Pcb罗氏线圈电流传感器基于集总参数的带宽快速计算方法 - Google Patents

Pcb罗氏线圈电流传感器基于集总参数的带宽快速计算方法 Download PDF

Info

Publication number
CN114563607A
CN114563607A CN202210087784.0A CN202210087784A CN114563607A CN 114563607 A CN114563607 A CN 114563607A CN 202210087784 A CN202210087784 A CN 202210087784A CN 114563607 A CN114563607 A CN 114563607A
Authority
CN
China
Prior art keywords
rogowski coil
pcb
turn
point
wire turn
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.)
Pending
Application number
CN202210087784.0A
Other languages
English (en)
Inventor
何怡刚
刘晓宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202210087784.0A priority Critical patent/CN114563607A/zh
Publication of CN114563607A publication Critical patent/CN114563607A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R15/00Details of measuring arrangements of the types provided for in groups G01R17/00 - G01R29/00, G01R33/00 - G01R33/26 or G01R35/00
    • G01R15/14Adaptations providing voltage or current isolation, e.g. for high-voltage or high-current networks
    • G01R15/18Adaptations providing voltage or current isolation, e.g. for high-voltage or high-current networks using inductive devices, e.g. transformers
    • G01R15/181Adaptations providing voltage or current isolation, e.g. for high-voltage or high-current networks using inductive devices, e.g. transformers using coils without a magnetic core, e.g. Rogowski coils
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/0092Arrangements for measuring currents or voltages or for indicating presence or sign thereof measuring current only
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Computing Systems (AREA)
  • Power Engineering (AREA)
  • Measuring Instrument Details And Bridges, And Automatic Balancing Devices (AREA)

Abstract

本发明公开了一种PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法,该方法包括:将PCB罗氏线圈内部的一点设为三维直角坐标系的原点,写出被测导体的参数方程;以PCB罗氏线圈的一个端点为起点,按照罗氏线圈内部电流的流向对PCB罗氏线圈线匝走线的交点和线圈两端点依次进行编号;将线匝点编号以及每个线匝点的坐标储存在矩阵中;分别计算每个线匝走线段与被测导体之间的互感;计算线匝走线段的自感系数、线匝段之间的电容、电阻和集肤深度;得到相邻两点之间的电气关系;将电气关系写成矩阵形式并将矩阵合并为一个方程,计算得到等效互感。本发明通用性强,计算方便,易于编制程序进行计算,节省设计时间;计算过程中没有经验公式,计算精度高。

Description

PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法
技术领域
本发明涉及交流电流测量技术领域,尤其涉及一种PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法。
背景技术
电流传感器在电气安全工程中有着重要作用。罗果夫斯基线圈电流传感器简称罗氏线圈(Rogowski coil)电流传感器,也叫微分电流传感器,以德国物理学家WalterRogowski来命名,已经应用了100多年。被认为可能是测量交流电流和瞬态大电流最好的测量工具。近年来,各种类型的基于印刷电路板(PCB)技术的罗氏线圈电流互感器不仅成为研究热点,由于其不容易变形、灵活的形状和尺寸、无磁饱、易于生产和低成本、以及不需要昂贵的绕线机就能生产出线匝分布均匀的罗氏线圈等优点被广泛应用于大电流、宽频带电流测量场合。罗氏线圈电流传感器的工作原理是基于安培定律和楞次定律。罗氏线圈与被测导体之间的互感决定了传感器的灵敏度的大小,并对传感器的测量精度有重要影响。罗氏线圈与临近的载流导体间的互感大小决定了罗氏线圈的灵敏度的大小以及抗干扰能力,而载流导线的位移对磁通量分布影响很大。
大多数文献中PCB罗氏线圈的计算方法是基于耗时的有限元理论计算方法。而采用电磁仿真软件对具有不规则骨架且密细线匝的PCB罗氏线圈进行仿真比较困难。罗氏线圈具有线匝间寄生电容、线匝自感、线匝间互感、载流导体与每个线匝之间的互感、线匝电阻、端口外接阻抗等诸多电气量。特别是在频率较高或绕线匝数较大时,这些电气量会影响罗氏线圈的测量带宽、灵敏度和测量精度。而常用的罗氏线圈集总参数电路模型仅适用于理想情况下的带宽计算方法。理想情况下的线圈骨架为规则的圆形,绕线线匝在骨架上分布均匀,载流导体为垂直穿过线圈几何中心的无限长直导体。而罗氏线圈的分布参数模型只是适合于上述理想工况且线圈匝数很少的罗氏线圈。然而,在罗氏线圈的实际应用中,罗氏线圈的匝数很多,以及任意形状的载流导体偏离线圈中心等情况。因此,研究一种适用于测量位于任意位置的载流导体时,针对任意骨架形状和任意线匝走线构成的密细线匝走线的PCB罗氏线圈测量带宽的计算方法具有重要意义。
发明内容
本发明要解决的技术问题在于针对现有技术中的缺陷,提供了一种计算方便、计算精度高、计算周期短的任意线匝走线结构的PCB罗氏线圈电流传感器基于集总参数的带宽的快速计算方法。
本发明解决其技术问题所采用的技术方案是:
本发明提供一种PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法,该方法包括以下步骤:
步骤1、将PCB罗氏线圈内部的一点设为三维直角坐标系的原点,若PCB罗氏线圈为几何对称的,将PCB罗氏线圈的几何中心设为三维直角坐标系的原点,PCB罗氏线圈几何对称轴为坐标系的z轴,写出被测导体的参数方程;
步骤2、PCB罗氏线圈的线匝走线为直线段,以PCB罗氏线圈的一个端点为起点,按照罗氏线圈内部电流的流向对PCB罗氏线圈线匝走线的交点和线圈两端点依次进行编号,若线匝走线段为M个,则线匝点编号的总数为M+1;
步骤3、将线匝点编号以及每个线匝点的坐标储存在矩阵EM+1,4中,矩阵EM+1,4的每一行4个元素分别存储一个线匝点的编号i和这个线匝点的直角坐标系的空间坐标(xi,yi,zi);
步骤4、写出被测导体的参数方程,分别计算每个线匝走线段与被测导体之间的互感,第i个线匝走线段与被测导体之间的互感为Mi
步骤5、由于PCB线圈的线匝走线很细,线匝走线的截面形状视为矩形,计算第m个线匝走线段的自感系数;
步骤6、考虑相邻的平行线匝走线段之间的电容,计算线匝段b(i-1)bi和线匝段b(m-1)bm之间的电容;
步骤7、计算电阻r,由于集肤效应,高频会改变走线的电阻,并计算集肤深度δ;
步骤8、b0,b1,b2,b3,…bi-1,bi,…,bm,…,bn,…bN-1,bN是存储在矩阵E中的交点元素序号,b0和bN是PCB罗氏线圈的两个端点,Zi是第i个直走线段的阻抗,Zi=ri+jωLi,Z0是线圈端口并联阻抗,得到相邻两点之间的电气关系;
步骤9、将矩阵E中的交点元素的电气关系写成矩阵形式;
步骤10、将步骤9得到的矩阵合并为一个方程,通过该方程计算得到等效互感MF。
进一步地,本发明的步骤2中的方法具体包括:
针对具有对称相同几何形状骨架的罗氏线圈:1)将其骨架按照其对称结构分为若干结构相同的骨架单元;将骨架单元中几何形状相同的线匝走线分为若干线匝单元;2)分别写出罗氏线圈的第一个骨架单元的第一个线匝单元的线匝交点的编号和坐标,以PCB罗氏线圈的一个端点为起点,按照罗氏线圈内部电流的流向对PCB罗氏线圈的第一个骨架单元的第一个线匝单元的线匝交点依次进行编号;3)根据步骤2)的坐标和骨架单元的对称结构的数学表达式写出第j个骨架单元的第i个线匝单元的线匝交点的坐标表达式;若每个线匝单元的线匝段端点的总数为N1,每个骨架单元的线匝单元数为N2,每个线圈的骨架单元为N3;第j个骨架单元的第i个线匝单元的相关点的坐标表达式为:
Figure BDA0003487759210000041
其中,θ’为x轴与坐标原点和第j个骨架单元几何中心在xOy面上的投影连线间的夹角;下标ij+s代表第j个骨架单元的第i个线匝单元的第s个点,其中ij=(j-1)*N2*N1+i1,s∈(1,2,3,…,N1),i1=N1*(i-1),也是这个点在整个PCB罗氏线圈的点的序列号;PCB罗氏线圈的线匝段的所有交点数量为M+1=N3*N2*N1;将序列号和N3*N2*N1个点的x,y,z坐标储存在一个N3*N2*N1行和4列的矩阵
Figure BDA0003487759210000043
中;线匝段的交点按照PCB罗氏线圈内部电流流向的顺序来计算和编号。
进一步地,本发明的步骤5中计算第m个线匝走线段的自感系数Lm的公式为:
Figure BDA0003487759210000042
其中,lm是第m个线段的长度,k1是线匝矩形截面的宽度,u0表示真空磁导率。
进一步地,本发明的步骤6中计算电容的公式为:
线匝段b(i-1)bi和线匝段b(m-1)bm之间的电容计算表达式为:
Figure BDA0003487759210000051
其中,ε是介电常数,Se和he分别是两个线匝等效面积和两个线匝走线段之间的距离,参数le是两个平行线匝走线段中最短线匝段的长度,w1表示走线的宽度。
进一步地,本发明的步骤7中电阻和集肤深度的计算公式为:
集肤深度δ为:
Figure BDA0003487759210000052
电阻r的计算公式为:
Figure BDA0003487759210000053
其中,ρ是电阻率,l和s分别为导体的长度和导体的截面积;因为PCB走线的截面形状为矩形,当集肤深度δ大于0.5*k1时,电流有效截面面积s近似等于厚度d1和宽度k1的乘积;由于PCB线通常很细,这里面积s近似等于导体的横截面面积减去由集肤效应引起的电流空心面积;如果δ小于0.5*k1并大于0.5*d1,面积s近似等于2δd1;如果δ小于0.5*d1,面积s近似等于k1d1-(d1-2δ)(k1-2δ);第i个线匝走线段的电阻ri通过计算矩阵中相邻两点之间的距离和PCB绕线截面的面积来获得,相邻两点之间的距离即线匝段两个端点的距离。
进一步地,本发明的步骤8中相邻两点之间的电气关系为:
Ui-1-Ui=jωMiI+ZiIi+jωLi,1I1+jωLi,2I2+…+jωLi,i-1Ii-1+…+jωLi,mIm+…+jωLi,NIN
Ii-Ii-1=jωCi-1,1(Ui-1-U1)+jωCi-1,2(Ui-1-U2)+…+jωCi-1,m(Ui-1-Um)+…+jωCi-1,N(Ui-1-UN)
其中,Ui和Ii-1分别是点bi相对于参考电位b0的电压和从点bi-1流到点bi的电流,Li,m和Ci,m分别是流过电流Ii的线匝走线段bi-1bi和流过电流Im的线匝走线段bm-1bm之间互感和电容,I是被测电流,第i个线匝走线段与被测导体之间的互感为Mi
进一步地,本发明的步骤9中的电气关系的矩阵形式为:
Figure BDA0003487759210000061
Figure BDA0003487759210000062
Figure BDA0003487759210000063
Figure BDA0003487759210000064
Figure BDA0003487759210000065
Figure BDA0003487759210000066
其中,B、C、Z表示计算过程的中间矩阵。
进一步地,本发明的步骤10中的方法具体包括:
将步骤9得到的矩阵合并为一个方程:
Figure BDA0003487759210000071
则矩阵元素VN与I之间的关系为:
Figure BDA0003487759210000072
其中,MF为不同频率下罗氏线圈端口输出电压VN与电流为I的载流导体之间的等效互感;令频率为f的电流I幅值为单位1时,可计算得到线圈两端的电压UN和电流I之间的比例关系;频率低于一定阈值时,MF保持不变,这时罗氏线圈工作在理想工作状态;如果没有PCB罗氏线圈的内部电容和电感参数时,等效互感参数MF在任意频率f下均保持不变,即罗氏线圈在任意频率下均工作在理想状态;当频率不断增大时,受线圈内部寄生电容、电感参数和端口并联阻抗参数影响的等效互感MF开始发生改变,这时罗氏线圈不再处于理想工作状态。
本发明产生的有益效果是:本发明的PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法,通用性强,适用于任意走线结构的高密度线匝的PCB罗氏线圈电流传感器的测量带宽。计算方便,不需要借助有限元等其他计算软件;易于编制程序进行计算,节省设计时间;计算过程中没有经验公式,计算精度高。最后,用实验验证了理论研究的正确性。
其能解决:1、现有常见的电磁仿真软件无法对高密度细走线的罗氏线圈进行网格刨分计算的问题;2、现有的理论计算方法仅适用于圆形骨架和方形骨架的稀疏线匝走线的PCB罗氏线圈,不适用于任意走线结构的高密度线匝的PCB罗氏线圈。
本发明提出的计算等效互感MF的方法对探究线圈带宽和测量精度的影响因素有重要意义。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1为本发明实施例高线匝密度的PCB罗氏线圈;
图2为本发明图1中PCB罗氏线圈的简图(a)图1所示的PCB罗氏线圈走线从上侧向下看视角图;(b)带有螺旋线圈的双层PCB小板的水平视角图;
图3为本发明PCB罗氏线圈的等效电路;
图4为本发明PCB罗氏线圈与载流导体的等效互感MF的计算结果和实验结果(10kHz~2MHz)(a)等效互感的幅值(b)等效互感的相位。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本发明实施例选取常见的细线匝走线且结构复杂的高线匝密度的PCB罗氏线圈,利用改进的基于集总参数模型来进行说明。
如图1所示,实施例中PCB罗氏线圈由36个40匝的两层PCB小板和2个PCB底座组成。2个PCB底座将这36个PCB小板固定并连接在一起。为了避免PCB线匝走线直角处的电磁辐射和电磁干扰,通常在线匝的转角处用钝角代替直角,PCB线匝走线的厚度d1和宽度k1分别为1oz和0.254mm。为了验证理论研究的正确性,设计的PCB罗氏线圈的线匝段不仅包括与载流导体平行的线匝走线段,还包括与导体不平行的任意方向的线匝走线段。选取载流直导体的长度为3m,实施例中PCB罗氏线圈在位于载流导体的中点处测量。从坐标原点(线圈中心)到带有螺旋线匝的PCB小板的最短距离为20mm。
步骤1、如图2(a)所示,图1的几何简图。将PCB罗氏线圈的几何中心设为三维直角坐标系的原点,PCB罗氏线圈几何对称轴为坐标系的z轴;
步骤2、为了更好地说明实施例的计算方法,以如图2(a)所示为图1的简图为例来进行描述,实线段和虚线段分别表示位于PCB底座顶部和PCB底座底部上的线匝走线段。如图2(b)所示,点o1,o2,o3,o4,o5,a1,a2,a3,a4,a5,a6,a7和a8位于最外圈的线匝上。h1是两个相邻线匝圈之间的距离。组成整个PCB罗氏线圈的带有螺旋线匝的PCB小板(如图2(b)所示)的数量是N2。而每个带有螺旋线匝的PCB小板的线匝圈数是N3。第j个线匝圈的线匝段交点可由最外圈(第一个线匝圈)线匝上的相应线匝段交点(a1,a2,a3,a4,a5,a6,a7和a8)来计算。例如,带有螺旋线匝的PCB小板的线匝段左下角的交点为:
Figure BDA0003487759210000091
其中,下标(j-1)*8+s代表点的序列号。位于角θ1'位置上的PCB小板的第j个线匝圈的第s个点的表达式可以通过(5-3)来计算,其在整个PCB罗氏线圈上的序列号是(θ1′N2/2π)(5+8N3)+3+8(j-1)+s。PCB罗氏线圈的线匝段的交点总数为N2(5+8N3)。将这N2(5+8N3)个交点的序列号和相应的x,y,z坐标储存在N2(5+8N3)行和4列的矩阵
Figure BDA0003487759210000101
里。矩阵E中PCB罗氏线圈所有线匝段的交点沿电流路径的方向排列(也沿电势递减或递增方向排列)。点Vi(xi,yi,zi)和Vi-1(xi-1,yi-1,zi-1)(1<i)分别是第i个线匝段的两个端点坐标,同时也是矩阵E里的元素。
步骤4、写出被测导体的参数方程,分别计算每个线匝走线段与被测导体之间的互感,第i个线匝走线段与被测导体之间的互感为Mi
步骤5、由于PCB线圈的线匝走线通常很细,线匝走线的截面形状为矩形,第m个线匝走线段的自感系数Lm
Figure BDA0003487759210000102
这里,lm是第m个线段的长度,k1是线匝矩形截面的宽度。
步骤6、为了增加线圈匝数,PCB罗氏线圈的线匝走线通常很细(例如,PCB线匝走线的厚度d1和宽度k1分别为1oz和0.254mm)。线匝走线之间的正对有效面积非常小。所以这里唯一要考虑的电容是相邻的平行线匝走线段之间的电容。线匝段b(i-1)bi和线匝段b(m-1)bm之间的电容计算表达式为:
Figure BDA0003487759210000103
其中,ε是介电常数,Se和he分别是两个线匝等效面积和两个线匝走线段之间的距离。参数le是两个平行线匝走线段中最短线匝段的长度,w1表示走线的宽度。
步骤7、计算电阻r,由于集肤效应,高频会改变走线的电阻。集肤深度δ为:
Figure BDA0003487759210000111
电阻r的计算公式为:
Figure BDA0003487759210000112
其中,ρ是电阻率,l和s分别为导体的长度和导体的截面积。因为PCB走线的截面形状为矩形,当集肤深度δ大于0.5*k1时,电流有效截面面积s近似等于厚度d1和宽度k1的乘积。由于PCB线通常很细,这里面积s近似等于导体的横截面面积减去由集肤效应引起的电流空心面积。如果δ小于0.5*k1并大于0.5*d1,面积s近似等于2δd1。如果δ小于0.5*d1,面积s近似等于k1d1-(d1-2δ)(k1-2δ)。第i个线匝走线段的电阻ri可以通过计算矩阵中相邻两点之间的距离(线匝段两个端点的距离)和PCB绕线截面的面积来获得。
步骤8、如图3所示,黄色的点b0,b1,b2,b3,…bi-1,bi,…,bm,…,bn,…bN-1,bN是存储在矩阵E中的交点元素序号。b0和bN是PCB罗氏线圈的两个端点。Zi(Zi=ri+jωLi)是第i个直走线段的阻抗。Z0是线圈端口并联阻抗,相邻两点之间的电气关系为:
Ui-1-Ui=jωMiI+ZiIi+jωLi,1I1+jωLi,2I2+…+jωLi,i-1Ii-1+…+jωLi,mIm+…+jωLi,NIN
Ii-Ii-1=jωCi-1,1(Ui-1-U1)+jωCi-1,2(Ui-1-U2)+…+jωCi-1,m(Ui-1-Um)+…+jωCi-1,N(Ui-1-UN)
其中,Ui和Ii-1分别是点bi相对于参考电位b0的电压和从点bi-1流到点bi的电流。Li,m和Ci,m分别是流过电流Ii的线匝走线段(线匝走线段bi-1bi)和流过电流Im的线匝走线段(线匝走线段bm-1bm)之间互感和电容。I是被测电流。
步骤9、将矩阵E中的交点元素的电气关系写成矩阵形式
Figure BDA0003487759210000121
Figure BDA0003487759210000122
Figure BDA0003487759210000123
Figure BDA0003487759210000124
Figure BDA0003487759210000125
Figure BDA0003487759210000126
其中,B、C、Z表示计算过程的中间矩阵。
步骤10、将步骤9合并为一个方程:
Figure BDA0003487759210000131
则矩阵元素VN与I之间的关系为:
Figure BDA0003487759210000132
其中,MF为不同频率下罗氏线圈端口输出电压VN与电流为I的载流导体之间的等效互感。令频率为f的电流I幅值为单位1时,线圈两端的电压UN和电流I之间的比例关系可以通过上式计算。频率较低时,MF保持不变,这时罗氏线圈工作在理想工作状态。如果没有PCB罗氏线圈的内部电容和电感参数时,等效互感参数MF在任意频率f下均保持不变。即罗氏线圈在任意频率下均工作在理想状态。当频率不断增大时,受线圈内部寄生电容、电感参数和端口并联阻抗等参数影响的等效互感MF开始发生改变,这时罗氏线圈不再处于理想工作状态。本发明提出的计算等效互感MF的方法对探究线圈带宽和测量精度的影响因素有重要意义。
当电流频率较大时,PCB罗氏线圈内部参数开始影响罗氏线圈的输出电压。如图4所示,参数MF随着频率的变化而变化,特别是在频率较高时。当参数MF的幅值和相位不随频率变化时,PCB罗氏线圈工作在理想状态。当频率在10kHz~2MHz范围内时,罗氏线圈的带宽的实验结果与理论计算结果一致。实验结果验证了理论研究的正确性。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims (8)

1.一种PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法,其特征在于,该方法包括以下步骤:
步骤1、将PCB罗氏线圈内部的一点设为三维直角坐标系的原点,若PCB罗氏线圈为几何对称的,将PCB罗氏线圈的几何中心设为三维直角坐标系的原点,PCB罗氏线圈几何对称轴为坐标系的z轴,写出被测导体的参数方程;
步骤2、PCB罗氏线圈的线匝走线为直线段,以PCB罗氏线圈的一个端点为起点,按照罗氏线圈内部电流的流向对PCB罗氏线圈线匝走线的交点和线圈两端点依次进行编号,若线匝走线段为M个,则线匝点编号的总数为M+1;
步骤3、将线匝点编号以及每个线匝点的坐标储存在矩阵EM+1,4中,矩阵EM+1,4的每一行4个元素分别存储一个线匝点的编号i和这个线匝点的直角坐标系的空间坐标(xi,yi,zi);
步骤4、写出被测导体的参数方程,分别计算每个线匝走线段与被测导体之间的互感,第i个线匝走线段与被测导体之间的互感为Mi
步骤5、由于PCB线圈的线匝走线很细,线匝走线的截面形状视为矩形,计算第m个线匝走线段的自感系数;
步骤6、考虑相邻的平行线匝走线段之间的电容,计算线匝段b(i-1)bi和线匝段b(m-1)bm之间的电容;
步骤7、计算电阻r,由于集肤效应,高频会改变走线的电阻,并计算集肤深度δ;
步骤8、b0,b1,b2,b3,…bi-1,bi,…,bm,…,bn,…bN-1,bN是存储在矩阵E中的交点元素序号,b0和bN是PCB罗氏线圈的两个端点,Zi是第i个直走线段的阻抗,Zi=ri+jωLi,Z0是线圈端口并联阻抗,得到相邻两点之间的电气关系;
步骤9、将矩阵E中的交点元素的电气关系写成矩阵形式;
步骤10、将步骤9得到的矩阵合并为一个方程,通过该方程计算得到等效互感MF。
2.根据权利要求1所述的PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法,其特征在于,所述步骤2中的方法具体包括:
针对具有对称相同几何形状骨架的罗氏线圈:1)将其骨架按照其对称结构分为若干结构相同的骨架单元;将骨架单元中几何形状相同的线匝走线分为若干线匝单元;2)分别写出罗氏线圈的第一个骨架单元的第一个线匝单元的线匝交点的编号和坐标,以PCB罗氏线圈的一个端点为起点,按照罗氏线圈内部电流的流向对PCB罗氏线圈的第一个骨架单元的第一个线匝单元的线匝交点依次进行编号;3)根据步骤2)的坐标和骨架单元的对称结构的数学表达式写出第j个骨架单元的第i个线匝单元的线匝交点的坐标表达式;若每个线匝单元的线匝段端点的总数为N1,每个骨架单元的线匝单元数为N2,每个线圈的骨架单元为N3;第j个骨架单元的第i个线匝单元的相关点的坐标表达式为:
Figure FDA0003487759200000021
其中,θ’为x轴与坐标原点和第j个骨架单元几何中心在xOy面上的投影连线间的夹角;下标ij+s代表第j个骨架单元的第i个线匝单元的第s个点,其中ij=(j-1)*N2*N1+i1,s∈(1,2,3,…,N1),i1=N1*(i-1),也是这个点在整个PCB罗氏线圈的点的序列号;PCB罗氏线圈的线匝段的所有交点数量为M+1=N3*N2*N1;将序列号和N3*N2*N1个点的x,y,z坐标储存在一个N3*N2*N1行和4列的矩阵
Figure FDA0003487759200000031
中;线匝段的交点按照PCB罗氏线圈内部电流流向的顺序来计算和编号。
3.根据权利要求1所述的PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法,其特征在于,所述步骤5中计算第m个线匝走线段的自感系数Lm的公式为:
Figure FDA0003487759200000032
其中,lm是第m个线段的长度,k1是线匝矩形截面的宽度,u0表示真空磁导率。
4.根据权利要求1所述的PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法,其特征在于,所述步骤6中计算电容的公式为:
线匝段b(i-1)bi和线匝段b(m-1)bm之间的电容计算表达式为:
Figure FDA0003487759200000033
其中,ε是介电常数,Se和he分别是两个线匝等效面积和两个线匝走线段之间的距离,参数le是两个平行线匝走线段中最短线匝段的长度,w1表示走线的宽度。
5.根据权利要求1所述的PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法,其特征在于,所述步骤7中电阻和集肤深度的计算公式为:
集肤深度δ为:
Figure FDA0003487759200000034
电阻r的计算公式为:
Figure FDA0003487759200000041
其中,ρ是电阻率,l和s分别为导体的长度和导体的截面积;因为PCB走线的截面形状为矩形,当集肤深度δ大于0.5*k1时,电流有效截面面积s近似等于厚度d1和宽度k1的乘积;由于PCB线通常很细,这里面积s近似等于导体的横截面面积减去由集肤效应引起的电流空心面积;如果δ小于0.5*k1并大于0.5*d1,面积s近似等于2δd1;如果δ小于0.5*d1,面积s近似等于k1d1-(d1-2δ)(k1-2δ);第i个线匝走线段的电阻ri通过计算矩阵中相邻两点之间的距离和PCB绕线截面的面积来获得,相邻两点之间的距离即线匝段两个端点的距离。
6.根据权利要求1所述的PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法,其特征在于,所述步骤8中相邻两点之间的电气关系为:
Ui-1-Ui=jωMiI+ZiIi+jωLi,1I1+jωLi,2I2+…+jωLi,i-1Ii-1+…+jωLi,mIm+…+jωLi,NIN
Ii-Ii-1=jωCi-1,1(Ui-1-U1)+jωCi-1,2(Ui-1-U2)+…+jωCi-1,m(Ui-1-Um)+…+jωCi-1,N(Ui-1-UN)
其中,Ui和Ii-1分别是点bi相对于参考电位b0的电压和从点bi-1流到点bi的电流,Li,m和Ci,m分别是流过电流Ii的线匝走线段bi-1bi和流过电流Im的线匝走线段bm-1bm之间互感和电容,I是被测电流,第i个线匝走线段与被测导体之间的互感为Mi
7.根据权利要求6所述的PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法,其特征在于,所述步骤9中的电气关系的矩阵形式为:
Figure FDA0003487759200000051
Figure FDA0003487759200000052
Figure FDA0003487759200000053
Figure FDA0003487759200000054
Figure FDA0003487759200000055
Figure FDA0003487759200000056
其中,B、C、Z表示计算过程的中间矩阵。
8.根据权利要求7所述的PCB罗氏线圈电流传感器基于集总参数的带宽快速计算方法,其特征在于,所述步骤10中的方法具体包括:
将步骤9得到的矩阵合并为一个方程:
Figure FDA0003487759200000061
则矩阵元素VN与I之间的关系为:
Figure FDA0003487759200000062
其中,MF为不同频率下罗氏线圈端口输出电压VN与电流为I的载流导体之间的等效互感;令频率为f的电流I幅值为单位1时,可计算得到线圈两端的电压UN和电流I之间的比例关系;频率低于一定阈值时,MF保持不变,这时罗氏线圈工作在理想工作状态;如果没有PCB罗氏线圈的内部电容和电感参数时,等效互感参数MF在任意频率f下均保持不变,即罗氏线圈在任意频率下均工作在理想状态;当频率不断增大时,受线圈内部寄生电容、电感参数和端口并联阻抗参数影响的等效互感MF开始发生改变,这时罗氏线圈不再处于理想工作状态。
CN202210087784.0A 2022-01-25 2022-01-25 Pcb罗氏线圈电流传感器基于集总参数的带宽快速计算方法 Pending CN114563607A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210087784.0A CN114563607A (zh) 2022-01-25 2022-01-25 Pcb罗氏线圈电流传感器基于集总参数的带宽快速计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210087784.0A CN114563607A (zh) 2022-01-25 2022-01-25 Pcb罗氏线圈电流传感器基于集总参数的带宽快速计算方法

Publications (1)

Publication Number Publication Date
CN114563607A true CN114563607A (zh) 2022-05-31

Family

ID=81714648

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210087784.0A Pending CN114563607A (zh) 2022-01-25 2022-01-25 Pcb罗氏线圈电流传感器基于集总参数的带宽快速计算方法

Country Status (1)

Country Link
CN (1) CN114563607A (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101916647A (zh) * 2010-08-23 2010-12-15 山东电力研究院 多频段测量数字式线性交流电流互感器及其测量方法
WO2014072760A1 (en) * 2012-11-06 2014-05-15 Freescale Semiconductor, Inc. High bandwidth current sensor and method therefor
CN106019072A (zh) * 2016-05-13 2016-10-12 国网四川省电力公司电力科学研究院 罗氏线圈集中参数的测量方法
CN106093547A (zh) * 2016-08-04 2016-11-09 贵州电网有限责任公司电力科学研究院 基于辅助线圈修正的罗氏线圈电流互感器测量装置及方法
WO2017197269A1 (en) * 2016-05-13 2017-11-16 The University Of North Carolina At Charlotte Wideband contactless magnetoresistive-rogowski current sensing
CN109100558A (zh) * 2018-06-06 2018-12-28 中国电力科学研究院有限公司 一种罗氏线圈及电流测量装置
CN109782062A (zh) * 2017-11-10 2019-05-21 许继集团有限公司 一种罗氏线圈型电子式电流互感器谐波计量特性优化方法
CN113219231A (zh) * 2020-02-05 2021-08-06 弗兰克公司 具有可调尺寸罗戈夫斯基线圈的非接触式电压测量

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101916647A (zh) * 2010-08-23 2010-12-15 山东电力研究院 多频段测量数字式线性交流电流互感器及其测量方法
WO2014072760A1 (en) * 2012-11-06 2014-05-15 Freescale Semiconductor, Inc. High bandwidth current sensor and method therefor
CN106019072A (zh) * 2016-05-13 2016-10-12 国网四川省电力公司电力科学研究院 罗氏线圈集中参数的测量方法
WO2017197269A1 (en) * 2016-05-13 2017-11-16 The University Of North Carolina At Charlotte Wideband contactless magnetoresistive-rogowski current sensing
CN106093547A (zh) * 2016-08-04 2016-11-09 贵州电网有限责任公司电力科学研究院 基于辅助线圈修正的罗氏线圈电流互感器测量装置及方法
CN109782062A (zh) * 2017-11-10 2019-05-21 许继集团有限公司 一种罗氏线圈型电子式电流互感器谐波计量特性优化方法
CN109100558A (zh) * 2018-06-06 2018-12-28 中国电力科学研究院有限公司 一种罗氏线圈及电流测量装置
CN113219231A (zh) * 2020-02-05 2021-08-06 弗兰克公司 具有可调尺寸罗戈夫斯基线圈的非接触式电压测量

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘晓宇: "《用于电流测量的罗氏线圈理论与应用研究》", 《万方数据-学位论文》 *

Similar Documents

Publication Publication Date Title
US20230084942A1 (en) Offset current sensor structure
US20140218018A1 (en) Apparatus and method for in situ current measurement in a conductor
US9983273B2 (en) Magnetic field detecting sensor and magnetic field detecting apparatus using the same
ES2327899T3 (es) Bobina de rogowski de alta precision.
CN105408756B (zh) 多分量磁场传感器
US20180372812A1 (en) Equilibrium-type magnetic field detection device
CN105143903B (zh) 磁场传感器系统
KR100417095B1 (ko) 기판의 배선검사장치 및 배선검사방법
US10704930B2 (en) Electromagnetic induction type displacement detection apparatus and measuring instrument using the same
US10739384B2 (en) PCB Rogowski coil
CN114563607A (zh) Pcb罗氏线圈电流传感器基于集总参数的带宽快速计算方法
CN112857194B (zh) 一种基于电涡流效应的平面二维位移传感器
US20160238645A1 (en) Device for measuring an electric field in a conducting medium
CN114563752A (zh) 一种pcb罗氏线圈电流传感器灵敏度的快速计算方法
JP2016161304A (ja) 抵抗器
Parry Helmholtz coils and coil design
US11162772B2 (en) Magnetic sensor device
CN110233019A (zh) 多层pcb结构三维磁场线圈
AU2015340864B2 (en) Sensor for measuring current in a conductor
CN114638193A (zh) 柔性罗氏线圈电流传感器基于集总参数的带宽快速计算方法
CN112557736B (zh) 一种导体面电流密度分布的测量系统及测量方法
CN104349595A (zh) 一种抗磁场干扰两块pcb开口罗氏线圈设计方法与实现
JP2015057582A (ja) 電流センサ
CN114814702B (zh) 双环骨架罗氏线圈电流传感器灵敏度的快速计算方法
JPH04283676A (ja) 多チャネルsquid磁束計

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20220531

RJ01 Rejection of invention patent application after publication