CN113720904B - 一种霍尔磁梯度张量结构设计方法和霍尔磁梯度张量结构 - Google Patents
一种霍尔磁梯度张量结构设计方法和霍尔磁梯度张量结构 Download PDFInfo
- Publication number
- CN113720904B CN113720904B CN202111048273.XA CN202111048273A CN113720904B CN 113720904 B CN113720904 B CN 113720904B CN 202111048273 A CN202111048273 A CN 202111048273A CN 113720904 B CN113720904 B CN 113720904B
- Authority
- CN
- China
- Prior art keywords
- axis
- magnetic field
- hall
- prewitt
- gradient tensor
- 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
- 230000005291 magnetic effect Effects 0.000 title claims abstract description 180
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000013461 design Methods 0.000 title claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 38
- 238000012545 processing Methods 0.000 claims abstract description 19
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 238000001514 detection method Methods 0.000 claims description 15
- 238000005339 levitation Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000005856 abnormality Effects 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 abstract description 21
- 238000005259 measurement Methods 0.000 description 38
- 238000004088 simulation Methods 0.000 description 6
- 230000002159 abnormal effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 230000007797 corrosion Effects 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000003822 epoxy resin Substances 0.000 description 2
- 230000005358 geomagnetic field Effects 0.000 description 2
- 229920000647 polyepoxide Polymers 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/72—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables
- G01N27/82—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws
- G01N27/83—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws by investigating stray magnetic fields
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/022—Measuring gradient
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/06—Measuring direction or magnitude of magnetic fields or magnetic flux using galvano-magnetic devices
- G01R33/07—Hall effect devices
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Electrochemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本申请涉及一种霍尔磁梯度张量结构设计方法和霍尔磁梯度张量结构。所述方法包括:在三维空间中,建立笛卡尔坐标系,笛卡尔坐标系以正方形的中心为原点,x轴的正方向朝北,y轴的正方向朝东,z轴的正方向朝下;将图像处理算法中的Prewitt梯度算子进行实例化,设计将十字形结构与正方形结构1:1结合的Prewitt结构;根据三轴霍尔传感器测量的磁场强度在x轴、y轴、z轴方向的分量和基线距离,采用Prewitt梯度算子进行计算,得到原点处的霍尔磁梯度张量矩阵。本设计方法的Prewitt结构可降低磁梯度张量系统对多传感器性能指标一致性的要求,适用于零点偏移大、灵敏度漂移范围广的磁梯度张量系统。
Description
技术领域
本申请涉及磁探测技术领域,特别是涉及一种霍尔磁梯度张量结构设计方法和霍尔磁梯度张量结构。
背景技术
与轮轨交通中单纯承载车辆的“轨道”不同,常导高速磁浮轨道实际上是一种同步直线电机的初级,更是悬浮系统、牵引系统和运控系统的重要电气组成部件。但是受车轨磕碰及环境腐蚀的影响,表面环氧树脂层破损的定子铁心容易进一步腐蚀、短路,导致牵引磁场异常,影响车辆运行的舒适性、稳定性及安全性。
行波漏磁场背景下的长定子铁心异常检测与地磁场背景下的矿藏勘探、潜艇侦查等磁异常检测极为相似,都存在严重的共模干扰。采用磁梯度张量测量的方式,可以抑制强背景噪声及测量系统本体姿态变化导致的共模干扰,提高对异常磁场的辨识能力。然而,与地磁场nT、uT级的异常检测不同,行波漏磁场的检测存在量级大(mT级)、安装空间小(mm级)的问题,传统超导量子干涉仪和磁通门传感器便无法使用,因此,需采用测量范围广(10-7~10T),集成尺寸小(1mm×1mm)的霍尔传感器构建小体积、mT量级的磁梯度张量测量系统。但霍尔传感器的精度更低,制造工艺也相对较差,它存在较大的不等位电压偏差、温漂、时漂等,导致存在较大的零点偏移与三轴灵敏度漂移。当使用多个传感器进行磁场梯度张量测量时,不可避免地存在多传感器性能参数不一致的情况,进而引入更大误差。目前,针对上述问题多采用针对单传感器的校正方法与针对多传感器测量系统的对准方法,这两种方法不但对校正环境与设备有较高要求,其校正策略与算法也较为复杂,同时也只是针对现有的十字形、正方形等测量结构,未从多传感空间结构设计的角度来减小系统测量误差。
现有的磁场梯度张量测量结构还不能应用于零点偏移大、灵敏度漂移范围广的霍尔磁梯度张量系统。
发明内容
基于此,有必要针对上述技术问题,提供一种霍尔磁梯度张量结构设计方法和霍尔磁梯度张量结构。采用该方法设计的霍尔磁梯度张量结构可以降低磁梯度张量系统对多传感器性能指标一致性的要求,同时也为传统磁梯度张量系统的结构设计提供了一种新的思路。
一种霍尔磁梯度张量结构设计方法,所述结构适用于高速磁浮轨道定子铁芯行波漏磁场的异常检测;所述方法包括:
在三维空间中,建立笛卡尔坐标系,所述笛卡尔坐标系以正方形的中心为原点,x轴的正方向朝北,y轴的正方向朝东,z轴的正方向朝下;
将图像处理算法中的Prewitt梯度算子进行实例化,设计将十字形结构与正方形结构1:1结合的Prewitt结构;所述Prewitt结构是在同一个平面内,十字形结构中相对两个三轴霍尔传感器连线的中点与笛卡尔坐标系的原点重合,十字形结构的四个三轴霍尔传感器在x轴或y轴上,正方形结构的相邻两个三轴霍尔传感器连线的中点与十字形的其中一个三轴霍尔传感器重合;正方形结构的基线距离和十字形结构的基线距离均等于的Prewitt结构的基线距离;
根据三轴霍尔传感器测量的磁场强度在x轴、y轴、z轴方向的分量和Prewitt结构的基线距离,采用Prewitt梯度算子进行计算,得到原点处的霍尔磁梯度张量矩阵。
一种霍尔磁梯度张量结构,所述霍尔磁梯度张量结构是上述霍尔磁梯度张量结构设计方法设计的;所述结构包括:8个三轴霍尔传感器、安装支架以及数据处理装置。
所述安装支架是边长与Prewitt结构的基线距离相等的正方形框架。
所述三轴霍尔传感器分别安装在正方形框架的四个顶点和正方形四条边的中点上。
所述数据处理装置用于根据采集的8个所述三轴霍尔传感器测量的x轴、y轴、z轴方向的磁场分量和结构的基线距离,采用Prewitt算子进行计算,得到原点处的霍尔磁梯度张量矩阵。
上述一种霍尔磁梯度张量结构设计方法和霍尔磁梯度张量结构,所述方法包括:在三维空间中,建立笛卡尔坐标系,笛卡尔坐标系以正方形的中心为原点,x轴的正方向朝北,y轴的正方向朝东,z轴的正方向朝下;将图像处理算法中的Prewitt梯度算子进行实例化,设计将十字形结构与正方形结构1:1结合的Prewitt结构;所述Prewitt结构是在同一个平面内,十字形结构中相对两个三轴霍尔传感器连线的中点与笛卡尔坐标系的原点重合,十字形结构的四个三轴霍尔传感器在x轴或y轴上,正方形结构的相邻两个三轴霍尔传感器连线的中点与十字形的其中一个三轴霍尔传感器重合;正方形结构的基线距离和十字形结构的基线距离均等于的Prewitt结构的基线距离;根据三轴霍尔传感器测量的磁场强度在x轴、y轴、z轴方向的分量和Prewitt结构的基线距离,采用Prewitt梯度算子进行计算,得到原点处的霍尔磁梯度张量矩阵。设计方法中采用价格低廉、体积小的霍尔传感器替代传统超导量子干涉仪、磁通门传感器构建了mT量级的磁梯度张量结构,可应用于零点偏移大、灵敏度漂移范围广的霍尔磁梯度张量系统;本方法设计的Prewitt结构可以降低磁梯度张量系统对多传感器性能指标一致性的要求,同时也为传统磁梯度张量系统的结构设计提供了一种新的思路。
附图说明
图1为一个实施例中搭载式磁场检测方法示意图;
图2为一个实施例中霍尔磁梯度张量结构设计方法的流程示意图;
图3为一个实施例中Prewitt结构的空间结构示意图;
图4为另一个实施例中零点偏移不一致对Ax的影响;
图5为一个实施例中灵敏度不一致对Ax的影响;
图6为另一个实施例中灵敏度漂移范围与Ax相对测量误差的关系。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
磁场测量的物理量主要有总磁场B的大小或模量B、模量的空间变化率、磁场的三个分量(Bx、By、Bz)和三个分量的空间变化率(Bxx、Bxy、Bxz、Byx、Byy、Byz、Bzx、Bzy、Bzz),这9个元素构成了二阶梯度张量。
常导高速磁浮轨道实际上是一种同步直线电机的初级,更是悬浮系统、牵引系统和运控系统的重要电气组成部件。但是受车轨磕碰及环境腐蚀的影响,表面环氧树脂层破损的定子铁芯容易进一步腐蚀、短路,导致牵引磁场异常,影响车辆运行的舒适性、稳定性及安全性。如图1所示,为实现动态、快速、高效的长定子状态监测,采用搭载式的磁敏阵列采集长定子产生的与空气闭合的行波漏磁场。该磁场幅值约30mT,易受悬浮磁场、列车姿态变化变化及外部噪声的干扰。本方法设计的霍尔磁梯度张量结构可作为测量行波漏磁场的磁敏阵列,用于高速磁浮轨道行波漏磁场的异常检测。
在一个实施例中,如图2所示,提供了霍尔磁梯度张量结构设计方法,所述霍尔磁梯度张量结构适用于高速磁浮轨道定子铁芯行波漏磁场的异常检测;霍尔磁梯度张量结构还适用于基于磁场的铁磁构件损伤检测、电机定子铁心短路检测等领域。该方法包括以下步骤:
步骤100:在三维空间中,建立笛卡尔坐标系。
笛卡尔坐标系以正方形的中心为原点,x轴的正方向朝北,y轴的正方向朝东,z轴的正方向朝下。
步骤102:将图像处理算法中的Prewitt梯度算子进行实例化,设计将十字形结构与正方形结构1:1结合的Prewitt结构。
Prewitt结构是在同一个平面内,十字形结构中相对两个三轴霍尔传感器连线的中点与笛卡尔坐标系的原点重合,十字形结构的四个三轴霍尔传感器在x轴或y轴上,正方形结构的相邻两个三轴霍尔传感器连线的中点与十字形的其中一个三轴霍尔传感器重合;正方形结构的基线距离和十字形结构的基线距离均等于的Prewitt结构的基线距离。
梯度计算在图像处理的边缘检测领域有着广泛的应用,常见的梯度检测算子如:Prewitt梯度算子,Prewitt梯度算子的两个卷积核为:
其中:Gx为水平方向的卷积核,Gy为垂直方向的卷积核。
从卷积核Gx和Gy结构上看,Prewitt对应的空间结构是十字形与正方形结构的1:1的结合。这种复合空间结构很好地保留了十字形与正方形的正交特性,在计算G各元素时能够引入更多的传感器参与,平均测量误差,同时,也并未大大提升硬件电路设计及元器件布局难度。Prewitt对应的空间结构如图3所示。
步骤104:根据三轴霍尔传感器测量的磁场强度在x轴、y轴、z轴方向的分量和Prewitt结构的基线距离,采用Prewitt梯度算子进行计算,得到原点处的霍尔磁梯度张量矩阵。
上述一种霍尔磁梯度张量结构设计方法中,所述方法包括:在三维空间中,建立笛卡尔坐标系,笛卡尔坐标系以正方形的中心为原点,x轴的正方向朝北,y轴的正方向朝东,z轴的正方向朝下;将图像处理算法中的Prewitt梯度算子进行实例化,设计将十字形结构与正方形结构1:1结合的Prewitt结构;所述Prewitt结构是在同一个平面内,十字形结构中相对两个三轴霍尔传感器连线的中点与笛卡尔坐标系的原点重合,十字形结构的四个三轴霍尔传感器在x轴或y轴上,正方形结构的相邻两个三轴霍尔传感器连线的中点与十字形的其中一个三轴霍尔传感器重合;正方形结构的基线距离和十字形结构的基线距离均等于的Prewitt结构的基线距离;根据三轴霍尔传感器测量的磁场强度在x轴、y轴、z轴方向的分量和Prewitt结构的基线距离,采用Prewitt梯度算子进行计算,得到原点处的霍尔磁梯度张量矩阵。设计方法中采用价格低廉、体积小的霍尔传感器替代传统超导量子干涉仪、磁通门传感器构建了mT量级的磁梯度张量结构,可应用于零点偏移大、灵敏度漂移范围广的霍尔磁梯度张量系统;本方法设计的Prewitt结构可以降低磁梯度张量系统对多传感器性能指标一致性的要求,同时也为传统磁梯度张量系统的结构设计提供了一种新的思路。
在其中一个实施例中,步骤104包括:将Prewitt结构的三轴霍尔传感器,按照从位于右上角的三轴霍尔传感器为起点顺时针方向进行编号;将三轴霍尔传感器测量的磁场强度在x轴、y轴、z轴方向的分量按照三轴霍尔传感器的编号顺序按照预定的规则分别表示为x轴、y轴、z轴方向的磁场强度分量矩阵;磁场强度分量矩阵是第二行第二列为0的3×3矩阵;预定的规则为:第1行第3列是编号为1的三轴霍尔传感器测量的磁场强度,按照顺时针方向,以此类推,第1行第2类是编号为8的三轴霍尔传感器测量的磁场强度;根据x轴、y轴、z轴方向的磁场强度分量矩阵和Prewitt结构的基线距离,采用Prewitt梯度算子进行计算,得到原点处的霍尔磁梯度张量矩阵。
在其中一个实施例中,步骤104中原点处的霍尔磁梯度张量矩阵的表达式为:
霍尔磁梯度张量矩阵中各个元素的计算公式为:
其中:G表示霍尔磁梯度张量矩阵;Bx、By、Bz分别表示原点处磁场强度在x轴、y轴、z轴方向的分量;d表示Prewitt结构的基线距离;Bxi、Byi、Bzi分别表示第i个三轴霍尔传感器测量x轴、y轴、z轴方向的磁场分量,其中i为大于等于1且小于等于8的整数;Bxx、Bxy、Bxz、Byy、Byz分别表示x方向磁场沿x方向的变化率、x方向磁场沿y方向的变化率、x方向磁场沿z方向的变化率,y方向磁场沿y方向的变化率及y方向磁场沿z方向的变化率。
应该理解的是,虽然图2的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图2中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在一个方法验证性仿真实施例中,传统磁异常检测方法常采用磁梯度张量解析信号Ax、Ay、Az极大值定位法,因此以Ax为例,计算十字形、正方形、Prewitt测量结构在多个传感器零点偏移与灵敏度不一致的情况下磁梯度张量解析信号测量误差的分布。其中,Ax的计算公式为:
(1)误差分析
理想的三轴霍尔传感器在无磁环境中应满足Bx=By=Bz=0,重点考虑传感器零点偏移与灵敏度不一致,以O-xyz为标准参考坐标系,则实测量B'=[B'x B'y B'z]T可表示为:
其中:o表示传感器零点偏移矩阵,k表示灵敏度矩阵,简记:B'=k·B+o。认为一批传感器的零点偏移o服从正态分布灵敏度k服从正态分布/>则实测量B':
考虑磁场在x,y,z三个方向的空间变化率,磁梯度张量矩阵G(magnetic gradienttensor)如式(1)所示,无源空间中磁场的散度与旋度均为0,因此G只有5个独立元素。
常见的十字形与正方形磁梯度张量测量结构,磁梯度张量矩阵G各元素B′ij同样服从正态分布:
矩阵各元素测量误差e(B′ij):
其中:n=1对应十字形结构,n=2对应正方形结构。
根据式(3)、(7),可知Ax服从χ2分布,以上侧α分位数c表示Ax的测量误差,取α=0.0001。
(2)仿真结果
仿真参数:均匀磁场Bx=By=Bz=30mT,自然噪声±0.5mT,霍尔传感器零点偏移o±1mT,灵敏度k:1±0.2。单传感器三轴灵敏度约束条件:Ky<4.8%,Kz<18.2%,单组仿真次数1000。
仿真结果如图4、图5所示,图中灵敏度一致是指:十字形结构在多个传感器不存在零点偏移及灵敏度漂移条件下的误差分布,并将其他情况相对灵敏度一致曲线的误差变化率定义为相对误差e,统计结果见表1。
表1不同测量结构误差
结果表明:当多传感器存在1mT范围内的零点偏移不一致时,Prewitt结构的测量误差小于参数一致情况,性能明显优于传统十字形、正方形结构,即Prewitt测量结构可完全消除传感器零点偏移带来的测量误差,这一点得益于Prewitt的冗余结构。当多传感器存在20%范围内的灵敏度不一致时,Prewitt结构的测量误差虽大于参数一致情况,但相对传统结构,其测量误差c可从340.5降低至168.4,约为理想状态的58%,测量结构引起的误差明显降低。如图6所示,Prewitt磁梯度张量测量结构理论上可允许系统内各传感器存在13%以内不同程度的灵敏度漂移,该指标覆盖了约95%的传感器。通过对测得信号进一步的滤波处理,可完全消除由传感器小概率参数漂移过大带来的影响,避免对复杂高成本的传感器标定筛选环节。
在一个实施例中,如图3所示,提供了一种霍尔磁梯度张量结构,该张量结构是采用上述霍尔磁梯度张量结构设计方法设计的;该结构包括:8个三轴霍尔传感器、安装支架以及数据处理装置。
安装支架是边长与Prewitt结构的基线距离相等的正方形框架。
三轴霍尔传感器分别安装在正方形框架的四个顶点和正方形四条边的中点上。
数据处理装置用于根据采集的8个三轴霍尔传感器测量的x轴、y轴、z轴方向的磁场分量和结构的基线距离,采用Prewitt算子进行计算,得到原点处的霍尔磁梯度张量矩阵。
霍尔磁梯度张量结构采用价格低廉、体积小的霍尔传感器替代传统超导量子干涉仪、磁通门传感器构建了mT量级的磁梯度张量系统;从结构设计的角度出发,将图像处理算法中的Prewitt梯度算子进行实例化,设计了将十字形与正方形结合的Prewitt结构;在高速磁浮轨道行波漏磁场异常检测的背景下,综合考虑所选霍尔传感器的各项性能指标,分别从理论和仿真的角度证明了Prewitt结构可以降低磁梯度张量系统对多传感器性能指标一致性的要求;最后,通过直线电机行波漏磁场测量实验,证明了Prewitt结构的测量误差更小,对传感器参数一致性的要求也更低。
在其中一个实施例中,数据处理装置,还用于将Prewitt结构的三轴霍尔传感器,按照从位于右上角的三轴霍尔传感器为起点顺时针方向进行编号;将三轴霍尔传感器测量的磁场强度在x轴、y轴、z轴方向的分量按照三轴霍尔传感器的编号顺序按照预定的规则分别表示为x轴、y轴、z轴方向的磁场强度分量矩阵;磁场强度分量矩阵是第二行第二列为0的3×3矩阵;预定的规则为:第1行第3列是编号为1的三轴霍尔传感器测量的磁场强度,按照顺时针方向,以此类推,第1行第2类是编号为8的三轴霍尔传感器测量的磁场强度;根据x轴、y轴、z轴方向的磁场强度分量矩阵和Prewitt结构的基线距离,采用Prewitt梯度算子进行计算,得到原点处的霍尔磁梯度张量矩阵。
在其中一个实施例中,数据处理装置中原点处的霍尔磁梯度张量矩阵的表达式如式(1)所示;霍尔磁梯度张量矩阵中各个元素计算公式如式(2)所示。
在一个传感器参数不一致对结构性能影响的验证性实施例中,为验证Prewitt结构相对十字形、正方形结构对传感器参数不一致不敏感的优势,采用标准磁源对霍尔传感器的零点偏移与灵敏度进行标定,筛选出各项参数相对集中的A组与参数相对分散的B组,分组构建Prewitt测量结构,进行对比实验。共标定传感器32只,A、B两组传感器参数分别见表1、2。A、B两组各8只,根据式(1)可实现磁梯度张量各元素的测量,同时可获得基线距离相同的十字形、正方形结构对应测量结果,与Prewitt进行对比。
表1 A组传感器参数
表2 B组传感器参数
利用变频器激励无转子的直线电机初级模拟高速磁浮轨道长定子产生的行波漏磁场,通过垂向移动模组控制磁梯度测量阵列与定子面的检测距离,驱动水平移动模组对直线电机定子面进行低速扫描,通过对交变磁场的进一步采样,获得搭载同步测量情况下行波漏磁场对应的磁梯度张量解析信号Ax,A、B两组不同测量结构测得Ax,对测量结果进行分析,得到参数不一致对Ax测量误差及信噪比的影响程度如表3所示。
表3参数不一致对Ax测量误差及信噪比的影响
十字形 | 正方形 | Prewitt结构 | |
测量误差变化率(%) | 63.2 | 60.1 | 16.9 |
信噪比变化率(%) | -7.0 | -5.4 | -1.7 |
实验表明:在磁梯度张量系统各传感器参数分散的情况下,Prewitt结构的测量误差与信噪比明显优于十字形与正方形结构,且与参数集中情况下接近。这说明Prewitt结构对霍尔传感器零点偏移大,灵敏度易漂移的特点不敏感,传感器可直接使用,不必进行筛选与误差校正。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (5)
1.一种霍尔磁梯度张量结构设计方法,所述结构适用于高速磁浮轨道定子铁芯行波漏磁场的异常检测;其特征在于,所述方法包括:
在三维空间中,建立笛卡尔坐标系,所述笛卡尔坐标系以正方形的中心为原点,x轴的正方向朝北,y轴的正方向朝东,z轴的正方向朝下;
将图像处理算法中的Prewitt梯度算子进行实例化,设计将十字形结构与正方形结构1:1结合的Prewitt结构;所述Prewitt结构是在同一个平面内,十字形结构中相对两个三轴霍尔传感器连线的中点与笛卡尔坐标系的原点重合,十字形结构的四个三轴霍尔传感器在x轴或y轴上,正方形结构的相邻两个三轴霍尔传感器连线的中点与十字形的其中一个三轴霍尔传感器重合;正方形结构的基线距离和十字形结构的基线距离均等于的Prewitt结构的基线距离;
根据三轴霍尔传感器测量的磁场强度在x轴、y轴、z轴方向的分量和Prewitt结构的基线距离,采用Prewitt梯度算子进行计算,得到原点处的霍尔磁梯度张量矩阵;
其中,根据三轴霍尔传感器测量的磁场强度在x轴、y轴、z轴方向的分量和Prewitt结构的基线距离,采用Prewitt梯度算子进行计算,得到原点处的霍尔磁梯度张量矩阵,包括:
将Prewitt结构的三轴霍尔传感器,按照从位于右上角的三轴霍尔传感器为起点顺时针方向进行编号;
将三轴霍尔传感器测量的磁场强度在x轴、y轴、z轴方向的分量按照三轴霍尔传感器的编号顺序按照预定的规则分别表示为x轴、y轴、z轴方向的磁场强度分量矩阵;所述磁场强度分量矩阵是第二行第二列为0的3×3矩阵;所述预定的规则为:第1行第3列是编号为1的三轴霍尔传感器测量的磁场强度,按照顺时针方向,以此类推,第1行第2类是编号为8的三轴霍尔传感器测量的磁场强度;
根据x轴、y轴、z轴方向的磁场强度分量矩阵和Prewitt结构的基线距离,采用Prewitt梯度算子进行计算,得到原点处的霍尔磁梯度张量矩阵;原点处的霍尔磁梯度张量矩阵的表达式为:
所述霍尔磁梯度张量矩阵中各个元素的计算公式为:
其中:G表示霍尔磁梯度张量矩阵;Bx、By、Bz分别表示原点处磁场强度在x轴、y轴、z轴方向的分量;d表示Prewitt结构的基线距离;Bxi、Byi、Bzi分别表示第i个三轴霍尔传感器测量x轴、y轴、z轴方向的磁场分量,其中i为大于等于1且小于等于8的整数;Bxx、Bxy、Bxz、Byy、Byz分别表示x方向磁场沿x方向的变化率、x方向磁场沿y方向的变化率、x方向磁场沿z方向的变化率,y方向磁场沿y方向的变化率及y方向磁场沿z方向的变化率。
2.根据权利要求1所述的方法,其特征在于,将图像处理算法中的Prewitt梯度算子进行实例化,设计将十字形结构与正方形结构结合的Prewitt结构。
3.一种霍尔磁梯度张量结构,其特征在于,所述霍尔磁梯度张量结构是采用权利要求1所述霍尔磁梯度张量结构设计方法设计的;所述结构包括:8个三轴霍尔传感器、安装支架以及数据处理装置;
所述安装支架是边长与Prewitt结构的基线距离相等的正方形框架;
所述三轴霍尔传感器分别安装在正方形框架的四个顶点和正方形四条边的中点上;
所述数据处理装置用于根据采集的8个所述三轴霍尔传感器测量的x轴、y轴、z轴方向的磁场分量和结构的基线距离,采用Prewitt算子进行计算,得到原点处的霍尔磁梯度张量矩阵。
4.根据权利要求3所述的结构,其特征在于,所述数据处理装置,还用于将Prewitt结构的三轴霍尔传感器,按照从位于右上角的三轴霍尔传感器为起点顺时针方向进行编号;将三轴霍尔传感器测量的磁场强度在x轴、y轴、z轴方向的分量按照三轴霍尔传感器的编号顺序按照预定的规则分别表示为x轴、y轴、z轴方向的磁场强度分量矩阵;所述磁场强度分量矩阵是第二行第二列为0的3×3矩阵;所述预定的规则为:第1行第3列是编号为1的三轴霍尔传感器测量的磁场强度,按照顺时针方向,以此类推,第1行第2类是编号为8的三轴霍尔传感器测量的磁场强度;根据x轴、y轴、z轴方向的磁场强度分量矩阵和Prewitt结构的基线距离,采用Prewitt梯度算子进行计算,得到原点处的霍尔磁梯度张量矩阵。
5.根据权利要求4所述的结构,其特征在于,所述数据处理装置中原点处的霍尔磁梯度张量矩阵的表达式为:
所述霍尔磁梯度张量矩阵中各个元素的计算公式为:
其中:G表示霍尔磁梯度张量矩阵;Bx、By、Bz分别表示原点处磁场强度在x轴、y轴、z轴方向的分量;d表示Prewitt结构的基线距离;Bxi、Byi、Bzi分别表示第i个三轴霍尔传感器测量x轴、y轴、z轴方向的磁场分量,其中i为大于等于1且小于等于8的整数;Bxx、Bxy、Bxz、Byy、Byz分别表示x方向磁场沿x方向的变化率、x方向磁场沿y方向的变化率、x方向磁场沿z方向的变化率,y方向磁场沿y方向的变化率及y方向磁场沿z方向的变化率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111048273.XA CN113720904B (zh) | 2021-09-08 | 2021-09-08 | 一种霍尔磁梯度张量结构设计方法和霍尔磁梯度张量结构 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111048273.XA CN113720904B (zh) | 2021-09-08 | 2021-09-08 | 一种霍尔磁梯度张量结构设计方法和霍尔磁梯度张量结构 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113720904A CN113720904A (zh) | 2021-11-30 |
CN113720904B true CN113720904B (zh) | 2024-05-10 |
Family
ID=78682425
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111048273.XA Active CN113720904B (zh) | 2021-09-08 | 2021-09-08 | 一种霍尔磁梯度张量结构设计方法和霍尔磁梯度张量结构 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113720904B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114609559B (zh) * | 2022-02-22 | 2023-06-02 | 电子科技大学 | 一种三轴霍尔角度传感器 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7038458B1 (en) * | 2004-10-12 | 2006-05-02 | The United States Of America As Represented By The Secretary Of The Navy | Magnetic anomaly homing system and method using rotationally invariant scalar contractions of magnetic gradient tensors |
CN102928885A (zh) * | 2012-11-22 | 2013-02-13 | 吉林大学 | 球形反馈三分量磁通门磁梯度全张量探头 |
KR20150102403A (ko) * | 2014-02-28 | 2015-09-07 | 주식회사 우진산전 | 마그네틱 바와 rfid 태그를 이용한 자기부상 열차용 열차 위치 검출장치 |
CN107817457A (zh) * | 2017-10-13 | 2018-03-20 | 北京工业大学 | 一种地磁梯度张量测量阵列的设计方法 |
CN108333551A (zh) * | 2018-02-14 | 2018-07-27 | 中国科学院电子学研究所 | 一种磁力计的校正方法 |
CN109270475A (zh) * | 2018-09-27 | 2019-01-25 | 中国人民解放军国防科技大学 | 高速磁浮长定子牵引行波磁场检测系统 |
CN109515479A (zh) * | 2018-11-16 | 2019-03-26 | 中国人民解放军国防科技大学 | 基于色标传感器的高速磁浮列车相对位置定位方法 |
CN110095738A (zh) * | 2019-04-15 | 2019-08-06 | 哈尔滨工业大学 | 磁梯度仪参数设计方法 |
CN112649766A (zh) * | 2020-11-27 | 2021-04-13 | 北京自动化控制设备研究所 | 基于全张量磁测系统实现合作磁目标识别的方法及系统 |
CN113267817A (zh) * | 2021-05-17 | 2021-08-17 | 国网福建省电力有限公司莆田供电公司 | 一种基于磁梯度张量的水下磁性物质的定位方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110007350B (zh) * | 2019-04-15 | 2020-08-25 | 哈尔滨工业大学 | 一种磁探测方法盲区的分析方法 |
-
2021
- 2021-09-08 CN CN202111048273.XA patent/CN113720904B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7038458B1 (en) * | 2004-10-12 | 2006-05-02 | The United States Of America As Represented By The Secretary Of The Navy | Magnetic anomaly homing system and method using rotationally invariant scalar contractions of magnetic gradient tensors |
CN102928885A (zh) * | 2012-11-22 | 2013-02-13 | 吉林大学 | 球形反馈三分量磁通门磁梯度全张量探头 |
KR20150102403A (ko) * | 2014-02-28 | 2015-09-07 | 주식회사 우진산전 | 마그네틱 바와 rfid 태그를 이용한 자기부상 열차용 열차 위치 검출장치 |
CN107817457A (zh) * | 2017-10-13 | 2018-03-20 | 北京工业大学 | 一种地磁梯度张量测量阵列的设计方法 |
CN108333551A (zh) * | 2018-02-14 | 2018-07-27 | 中国科学院电子学研究所 | 一种磁力计的校正方法 |
CN109270475A (zh) * | 2018-09-27 | 2019-01-25 | 中国人民解放军国防科技大学 | 高速磁浮长定子牵引行波磁场检测系统 |
CN109515479A (zh) * | 2018-11-16 | 2019-03-26 | 中国人民解放军国防科技大学 | 基于色标传感器的高速磁浮列车相对位置定位方法 |
CN110095738A (zh) * | 2019-04-15 | 2019-08-06 | 哈尔滨工业大学 | 磁梯度仪参数设计方法 |
CN112649766A (zh) * | 2020-11-27 | 2021-04-13 | 北京自动化控制设备研究所 | 基于全张量磁测系统实现合作磁目标识别的方法及系统 |
CN113267817A (zh) * | 2021-05-17 | 2021-08-17 | 国网福建省电力有限公司莆田供电公司 | 一种基于磁梯度张量的水下磁性物质的定位方法 |
Non-Patent Citations (2)
Title |
---|
十字形磁梯度张量系统的误差校正;迟铖;吕俊伟;黄婧丽;;光学精密工程;20170715(07);全文 * |
基于Prewitt 算子的自适应弱小目标检测;许四祥 等;红外技术;20190531;189-193 * |
Also Published As
Publication number | Publication date |
---|---|
CN113720904A (zh) | 2021-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106405658B (zh) | 一种基于矢量磁梯度计的运动式磁性目标定位方法 | |
US20210161420A1 (en) | Magnetic field measuring apparatus, magnetic field measuring method, and recording medium storing magnetic field measuring program | |
CN103080755B (zh) | 电流传感器 | |
CN105891755B (zh) | 航空器挂载式磁通门磁梯度张量仪的校正方法 | |
DE102017110197A1 (de) | Außeraxiale Magnetfeld-Winkelsensoren | |
DE102014116844A1 (de) | Außeraxiale Magnetfeld-Winkelsensoren | |
JP7487156B2 (ja) | 磁場計測装置 | |
CN110146839A (zh) | 一种移动平台磁梯度张量系统校正方法 | |
JP7222814B2 (ja) | 測定装置および測定方法 | |
CN113720904B (zh) | 一种霍尔磁梯度张量结构设计方法和霍尔磁梯度张量结构 | |
CN112050799B (zh) | 一种基于磁梯度张量缩并量之比的测距定位方法 | |
CN115542225B (zh) | 一种提升悬丝扭秤装置磁性测量精度的校正方法 | |
CN112050800A (zh) | 一种基于日字型三轴磁强计对称配置平面阵列的磁梯度张量定位方法 | |
US20220057470A1 (en) | Magnetic field calibration device and method of calibrating magnetism measurement device using the same | |
CN113341195A (zh) | 一种基于隧道磁阻元件的电流测量方法、测量装置及设备 | |
CN110596619B (zh) | 一种全张量磁梯度测量组件及其优化方法 | |
CN114444299B (zh) | 一种基于距离加权多极展开方法的磁场重构方法 | |
Gheorghe et al. | Calibration of skew redundant sensor configurations in the presence of field alignment errors | |
Xu et al. | Factor analysis on integrated error model of magnetic gradient tensor | |
CN113608152B (zh) | 磁传感器 | |
US20240133714A1 (en) | Device and method for determining an encoder magnet rotation angle | |
CN115526121B (zh) | 基于单轴tmr传感信号深度学习的电流计算方法及系统 | |
US20230176147A1 (en) | Magnetic sensor, electric control device, correction method, and magnetic sensor manufacturing method | |
US11598827B2 (en) | Method of designing and manufacturing magnetic sensor | |
Yuan et al. | Small current measurement based on cross-type magnetic sensor array |
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 |