CN110132188B - 一种基于多元素x射线特征光谱综合分析的涂渗层厚度计算方法 - Google Patents

一种基于多元素x射线特征光谱综合分析的涂渗层厚度计算方法 Download PDF

Info

Publication number
CN110132188B
CN110132188B CN201910533596.4A CN201910533596A CN110132188B CN 110132188 B CN110132188 B CN 110132188B CN 201910533596 A CN201910533596 A CN 201910533596A CN 110132188 B CN110132188 B CN 110132188B
Authority
CN
China
Prior art keywords
thickness
coating
analyzed
fluorescence intensity
comprehensive analysis
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
Application number
CN201910533596.4A
Other languages
English (en)
Other versions
CN110132188A (zh
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.)
Jiangsu Elite Instruments Co ltd
Original Assignee
Air Force Engineering University of PLA
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 Air Force Engineering University of PLA filed Critical Air Force Engineering University of PLA
Priority to CN201910533596.4A priority Critical patent/CN110132188B/zh
Publication of CN110132188A publication Critical patent/CN110132188A/zh
Application granted granted Critical
Publication of CN110132188B publication Critical patent/CN110132188B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B15/00Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
    • G01B15/02Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons for measuring thickness
    • G01B15/025Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons for measuring thickness by measuring absorption

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明提出一种基于多元素X射线特征光谱综合分析的涂渗层厚度计算方法,基于X射线荧光探测技术中的吸收法,利用了多种X射线特征光谱对涂渗层厚度敏感的元素,对多种元素进行综合分析,建立多元素综合分析下的计算方法,保证了测试结果的准确性。其创新之处在于利用各种元素之间的相互补偿,避免了单一元素可能带来的偶然性误差,提高测试结果的准确性。

Description

一种基于多元素X射线特征光谱综合分析的涂渗层厚度计算 方法
技术领域
本发明涉及一种涂渗层厚度计算方法,属于无损检测领域。
背景技术
X射线荧光测试技术是目前广泛使用的一种无损检测技术,具有方便、快捷、准确的特点,运用前景十分广阔。但是,就目前来看,这种技术主要运用于有明显界限单一均匀镀层的无损检测,在扩散型渗层厚度的检测上基本上没有相关研究。如果仅使用原有的涂层厚度计算方法对扩散型渗层厚度进行研究,测试得到的结果与实际值就偏差较大,准确性受到极大的限制。
目前,对于扩散型渗层厚度的检测主要依赖于断面扫描的方法,即将试样线切割,然后通过电镜进行扫描观察,根据形貌的不同测出渗层的厚度,但是通过此方法不仅耗时长,破坏材料的完整性,还只能对部分产品进行抽样检查,抽样检查结果偶然性较大。而采用单一渗层元素进行X射线厚度检测,其缺点是精确度较低或需要大量采集标准试件数据。
发明内容
为解决现有技术存在的问题,本发明提出一种基于多元素X射线特征光谱综合分析的涂渗层厚度计算方法,基于X射线荧光探测技术中的吸收法,利用了多种X射线特征光谱对涂渗层厚度敏感的元素,对多种元素进行综合分析,建立多元素综合分析下的计算方法,保证了测试结果的准确性。其创新之处在于利用各种元素之间的相互补偿,避免了单一元素可能带来的偶然性误差,提高测试结果的准确性。
本发明的技术方案为:
所述一种基于多元素X射线特征光谱综合分析的涂渗层厚度计算方法,其特征在于:包括以下步骤:
步骤1:基于以下规则在材料中选择待分析元素:
所选择的元素满足:1)、元素的荧光光子能够穿透涂渗层,到达探测器;2)、元素的荧光强度随厚度的变化不小于设定要求;
步骤2:获取步骤1中选择的n个待分析元素在材料中的浓度,并分别记为Ci,i=1,2…n;并对含量进行校正,消除基体效应,校正后的各个元素浓度记为Ci′;并对校正后的各个元素浓度进行归一化处理,得到归一化后的各个元素浓度为
Figure BDA0002100498370000021
步骤3:对于第i个待分析元素,获取其在不同涂渗层厚度t下的强度比Ri,其中
Figure BDA0002100498370000022
其中Ii为探测器接收到的第i个待分析元素在涂渗层厚度t下的X射线荧光强度,Ii0为涂渗层厚度为0时探测器接收到的第i个待分析元素的X射线荧光强度;
步骤4:根据步骤3得到的不同涂渗层厚度t下的强度比Ri,利用公式y=-kx进行线性拟合,其中y对应lnRi,x对应t,拟合得到第i个待分析元素对应的μiρ=k;
步骤5:对于n个待分析元素,重复步骤3和步骤4,得到每个待分析元素的μiρ值;
步骤6:利用步骤5得到的每个待分析元素的μiρ值,建立多元素综合分析下的合金涂渗层的厚度优化计算模型为:
Figure BDA0002100498370000023
将n个待分析元素对应的Ii、Ii0、Ci″、μiρ代入模型即可计算得到多元素综合分析下的合金涂渗层厚度t。
进一步的优选方案,所述一种基于多元素X射线特征光谱综合分析的涂渗层厚度计算方法,其特征在于:步骤1中判断元素的荧光强度随厚度的变化不小于设定要求的具体过程为:
建立坐标系,横坐标为涂渗层厚度,纵坐标为当前厚度下该元素的荧光强度与无渗层的荧光强度之比的对数;
该元素在无渗层的荧光强度为I0,不同厚度D1,D2…Dn下的荧光强度分别为ID1,ID2…IDn,若(D1,ln(ID1/I0)),(D2,ln(ID2/I0))…(Dn,ln(IDn/I0))这些点与原点(0,0)点连线的斜率变化在设定值以内,则认为该元素的荧光强度随涂渗层厚度的变化不小于设定要求。
有益效果
本发明提供了一种基于吸收法的多元素综合分析涂渗层厚度的优化计算模型,先通过对满足理论计算要求的多种元素进行单独分析,然后再将这些满足要求的元素进行综合分析,得到一个精度高并且能够适用于各类涂渗层厚度计算的计算模型,有效解决了现有的计算方法对某些涂渗层厚度计算不准确的问题,提高了测试精度。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1多元素X射线荧光强度测试图。
具体实施方式
本发明以单一元素特征光谱吸收法计算涂渗层厚度为基础,将多种元素进行综合分析,得到一个全新的基于多种元素综合的涂渗层厚度计算模型。
本发明的原理及具体过程如下:
1)多元素的选取
根据Mosely定律:元素的X射线特征波长倒数的平方根与原子序数成正比。当原子序数较低时,激发其K系电子所需的X能量比较低。其激发出来的的特征X荧光光子能量同样较低,容易被其他原子序数更大的原子吸收,降低其发射能力。材料中的某些元素荧光强度较强,当厚度变化时,探测器接收到的元素的荧光强度随厚度的变化不明显或者基本没有变化;另外一些元素由于元素本身的特征荧光能量较小,探测器基本都接收不到。因此,必须对元素种类进行分析和筛选,选取合适的元素作为建立数学模型的基础。
这些元素的选取两个条件是:a)元素的荧光光子能够穿透涂渗层,到达探测器;b)荧光强度随厚度有一个较为明显的变化。具体的选择方法为:假设无渗层的荧光强度为I0,不同厚度D1,D2…Dn下的荧光强度分别为ID1,ID2…IDn,则(D1,ln(ID1/I0)),(D2,ln(ID2/I0))…(Dn,ln(IDn/I0))这些点与坐标系原点(0,0)点连线的斜率变化在设定值以内,设定值一般取10%。
2)多元素综合计算方法
通过步骤1)分析可以得到,假设有n个元素满足计算要求,这些元素在材料中所占的百分含量分别记为Ci(i=1,2…n)。
1、基体效应的校正
由于基体效应的存在,元素的浓度和分析线的强度之间的关系一般不是线性的。因此需要通过一些方法来消除、减少或校正基体效应,提高测试的精密度和准确度。在这里主要运用常见的数学校正法来对基体效应进行校正。通过校正使测得的X射线荧光强度与计算出来的理论X射线荧光强度相一致,然后将校正得到的浓度值保存下来,为后续多元素综合分析奠定数据和理论基础。将校正后的各个元素浓度记为Ci′(i=1,2…n)。
2、单元素参数拟合分析计算模型
以吸收法计算公式I=I0e-μρt=I0e-μm为理论基础。
其中,t为涂渗层厚度,m为单位截面积内涂渗层材料的质量,I为探测器接收到的基体元素X射线荧光强度,I0为涂渗层厚度为0时探测器接收到的基体元素X射线荧光强度,μ为涂渗层材料的质量吸收系数。
基底元素所产生的X射线荧光穿透涂渗层时会产生强度衰减,并且随着涂渗层厚度的增加,衰减的程度逐渐增加,探测器所能探测到的荧光强度的大小也有所减小。
下面以i(i=1,2…n)元素为例对参数拟合过程进行推导:
首先,将吸收法计算公式做一个变形,将等式两边同时除以Ii0可以得到,
Figure BDA0002100498370000041
下标i表示该符号为i元素的对应符号;再将等式两边取对数,得到
ln(Ii/Ii0)=-μiρt=-μim (1)
Figure BDA0002100498370000051
将Ri记做i元素的强度比。
通过分析得到,该表达式仅仅是一个与厚度t相关的一次表达式,形如y=-kx。其中,系数k对应的是μiρ,y对应的lnRi,x对应的t,所以对该方程求解时只需要对吸收系数μi进行拟合即可得到最后的计算表达式。即通过一元线性回归的方法对lnRi与t拟合得到k=μiρ,然后计算得到质量吸收系数的μi值。
通过以上计算可以得到i元素基于吸收法计算涂渗层厚度的计算模型。
Figure BDA0002100498370000052
同理,可以得到与筛选出来的敏感元素相关的计算模型,
Figure BDA0002100498370000053
3、多元素综合分析计算模型
以单元素的计算模型和校正后的浓度值为基础,对多种元素综合下的优化计算模型进行分析和推导。
首先,将(3)式得到的单元素计算模型进行综合,得到
Figure BDA0002100498370000054
然后,将元素i校正后的浓度Ci′做归一化处理,得到
Figure BDA0002100498370000055
其次,将(4)式开方可以得到
Figure BDA0002100498370000056
再次,用(5)式得到的归一化后各元素的浓度值Ci″对(6)式中各元素对应的质量系数修正,即
Figure BDA0002100498370000057
最后,将(7)式带入(6)式即可得到多元素综合下的涂渗层厚度优化计算模型,即
Figure BDA0002100498370000061
(8)式即作为多元素综合分析下的优化计算模型。
下面详细描述本发明的实施例,所述实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
本实施例以某型合金为例进行测试计算:
1)进行X荧光测试,得到如表1所示为Ti和Cr荧光强度的强度比值,通过对强度比值取对数的值如表2所示,选取Ti和Cr为计算元素。
表1不同厚度涂渗层的某合金Ti和Cr荧光强度的强度比值
Figure BDA0002100498370000062
表2强度比取对数后的值
Figure BDA0002100498370000063
2)根据以上数据计算该合金表面涂渗层计算公式,具体计算步骤如下:
1、基体效应校正后的浓度CTi′=0.14,CCr′=0.23,根据校正以后的浓度值带入(5)式进行计算得到CTi″=0.378,CCr″=0.622。
2、对Ti元素通过一元线性回归的方法拟合得到μTiρ=-0.1682,然后,根据公式(1)得到
ln(ITi/ITi0)=-μTiρt=-0.1682t (9)
对Cr元素通过一元线性回归的方法拟合得到μCrρ=-0.3210,然后根据式(1)得到
ln(ICr/ICr0)=-μCrρt=-0.3210t (10)
3、多元素综合后的计算
将式(8)等价变化以后得到
Figure BDA0002100498370000071
通过上式带入第2步和第3步的计算结果得到该种合金涂渗层的计算方程为
Figure BDA0002100498370000072
方程(12)即为多元素综合下的涂渗层厚度的计算方程。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (2)

1.一种基于多元素X射线特征光谱综合分析的涂渗层厚度计算方法,其特征在于:包括以下步骤:
步骤1:基于以下规则在材料中选择待分析元素:
所选择的元素满足:1)、元素的荧光光子能够穿透涂渗层,到达探测器;2)、元素的荧光强度随厚度的变化不小于设定要求;
步骤2:获取步骤1中选择的n个待分析元素在材料中的浓度,并分别记为Ci,i=1,2…n;并对含量进行校正,消除基体效应,校正后的各个元素浓度记为Ci′;并对校正后的各个元素浓度进行归一化处理,得到归一化后的各个元素浓度为
Figure FDA0002668772020000011
步骤3:对于第i个待分析元素,获取其在不同涂渗层厚度t下的强度比Ri,其中
Figure FDA0002668772020000012
其中Ii为探测器接收到的第i个待分析元素在某一涂渗层厚度下的X射线荧光强度,Ii0为涂渗层厚度为0时探测器接收到的第i个待分析元素的X射线荧光强度;
步骤4:根据步骤3得到的不同涂渗层厚度t下的强度比Ri,利用公式y=-kx进行线性拟合,其中y对应lnRi,x对应t,拟合得到第i个待分析元素对应的μiρ=k,其中μi表示第i个待分析元素的质量吸收系数;
步骤5:对于n个待分析元素,重复步骤3和步骤4,得到每个待分析元素的μiρ值;
步骤6:利用步骤5得到的每个待分析元素的μiρ值,建立多元素综合分析下的涂渗层的厚度优化计算模型为:
Figure FDA0002668772020000013
将n个待分析元素对应的Ii、Ii0、Ci″、μiρ代入模型即可计算得到多元素综合分析下的涂渗层厚度T。
2.根据权利要求1所述一种基于多元素X射线特征光谱综合分析的涂渗层厚度计算方法,其特征在于:步骤1中判断元素的荧光强度随厚度的变化不小于设定要求的具体过程为:
建立坐标系,横坐标为涂渗层厚度,纵坐标为当前厚度下该元素的荧光强度与无渗层的荧光强度之比的对数;
该元素在无渗层的荧光强度为I0,不同厚度D1,D2…Dn下的荧光强度分别为ID1,ID2…IDn,若(D1,ln(ID1/I0)),(D2,ln(ID2/I0))…(Dn,ln(IDn/I0))这些点与原点(0,0)点连线的斜率变化在设定值以内,则认为该元素的荧光强度随涂渗层厚度的变化不小于设定要求。
CN201910533596.4A 2019-06-19 2019-06-19 一种基于多元素x射线特征光谱综合分析的涂渗层厚度计算方法 Active CN110132188B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910533596.4A CN110132188B (zh) 2019-06-19 2019-06-19 一种基于多元素x射线特征光谱综合分析的涂渗层厚度计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910533596.4A CN110132188B (zh) 2019-06-19 2019-06-19 一种基于多元素x射线特征光谱综合分析的涂渗层厚度计算方法

Publications (2)

Publication Number Publication Date
CN110132188A CN110132188A (zh) 2019-08-16
CN110132188B true CN110132188B (zh) 2020-11-10

Family

ID=67578119

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910533596.4A Active CN110132188B (zh) 2019-06-19 2019-06-19 一种基于多元素x射线特征光谱综合分析的涂渗层厚度计算方法

Country Status (1)

Country Link
CN (1) CN110132188B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109556541A (zh) * 2019-01-14 2019-04-02 汪诚 一种基于x射线的金属表面渗层厚度无损检测装置及方法
CN110530912B (zh) * 2019-09-12 2022-01-04 岛津企业管理(中国)有限公司 一种含镀层贵金属成分的x射线荧光光谱分析方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102187208A (zh) * 2009-09-07 2011-09-14 株式会社理学 荧光x射线分析方法
JP2011185881A (ja) * 2010-03-11 2011-09-22 Hitachi Maxell Ltd 層厚測定装置
CN104880475A (zh) * 2014-02-28 2015-09-02 赫尔穆特费希尔有限责任公司电子及测量技术研究所 一种用于借助于x射线荧光进行测量对象的测量的方法
CN108226201A (zh) * 2016-12-22 2018-06-29 马尔文帕纳科公司 用xrf分析层状样品

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8811459D0 (en) * 1988-05-13 1988-06-15 Dmc Boyle Ltd Method & apparatus for measuring thickness of coating on substrate
BR9916368A (pt) * 1998-12-21 2001-10-09 Coring Inc Análise de emissão fluorescente de raio x para determinar concentração de material
JP2003050115A (ja) * 2001-08-07 2003-02-21 Seiko Instruments Inc X線膜厚計
US6507634B1 (en) * 2001-09-19 2003-01-14 Therma-Wave, Inc. System and method for X-ray reflectometry measurement of low density films
DE602005013435D1 (de) * 2005-09-26 2009-04-30 Jfe Steel Corp Verfahren zur messung der filmdicke eines oberflächenoxidfilms eines mit zink galvanisierten stahlblechs
US8513603B1 (en) * 2010-05-12 2013-08-20 West Virginia University In-situ determination of thin film and multilayer structure and chemical composition using x-ray fluorescence induced by grazing incidence electron beams during thin film growth
MX361123B (es) * 2012-04-25 2018-11-28 Nippon Steel & Sumitomo Metal Corp Método dispositivo para determinar el espesor de la fase de aleación de hierro - zinc de una lámina de acero galvanizada por inmersión en caliente.
GB201214448D0 (en) * 2012-08-14 2012-09-26 Element Six Ltd Electrochemical deposition and x-ray fluorescence spectroscopy
WO2014160637A1 (en) * 2013-03-26 2014-10-02 Momentive Performance Materials Inc. Coating composition and method for determining the uniformity and thickness of a no-rinse silane pretreatment
JP6172054B2 (ja) * 2014-05-29 2017-08-02 株式会社島津製作所 蛍光x線分析装置
US9784699B2 (en) * 2015-03-03 2017-10-10 Panalytical B.V. Quantitative X-ray analysis—matrix thickness correction
CN105203054B (zh) * 2015-10-28 2018-01-12 武汉理工大学 基于γ射线的柴油机缸套活塞之间油膜测厚模拟测试平台
JP6641898B2 (ja) * 2015-11-04 2020-02-05 日本製鉄株式会社 放射線厚み測定装置及びその校正方法
CN108508052B (zh) * 2018-06-11 2023-10-20 西北核技术研究所 基于参考元素的x射线荧光薄层质量厚度测量系统及方法
CN109556541A (zh) * 2019-01-14 2019-04-02 汪诚 一种基于x射线的金属表面渗层厚度无损检测装置及方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102187208A (zh) * 2009-09-07 2011-09-14 株式会社理学 荧光x射线分析方法
JP2011185881A (ja) * 2010-03-11 2011-09-22 Hitachi Maxell Ltd 層厚測定装置
CN104880475A (zh) * 2014-02-28 2015-09-02 赫尔穆特费希尔有限责任公司电子及测量技术研究所 一种用于借助于x射线荧光进行测量对象的测量的方法
CN108226201A (zh) * 2016-12-22 2018-06-29 马尔文帕纳科公司 用xrf分析层状样品

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于X荧光涂层厚度测量的试验研究;黄明 等;《核电子学与探测技术》;20140720;第34卷(第7期);第866-868、873页 *

Also Published As

Publication number Publication date
CN110132188A (zh) 2019-08-16

Similar Documents

Publication Publication Date Title
CN110312928B (zh) 荧光x射线分析方法以及荧光x射线分析装置
Tardio et al. Comparative study of the native oxide on 316L stainless steel by XPS and ToF-SIMS
Ernst et al. Signal‐to‐noise ratios in forensic glass analysis by micro X‐ray fluorescence spectrometry
CN110132188B (zh) 一种基于多元素x射线特征光谱综合分析的涂渗层厚度计算方法
CA2393885C (en) Device for measurement of metal sheet thickness and clad layer thickness and use of the device
Malzer et al. A fingerprint model for inhomogeneous ink–paper layer systems measured with micro‐x‐ray fluorescence analysis
WO2015056305A1 (ja) 蛍光x線分析方法及び蛍光x線分析装置
Bonetto et al. Optimization of parameters in electron probe microanalysis
Queralt et al. Thickness measurement of semiconductor thin films by energy dispersive X-ray fluorescence benchtop instrumentation: Application to GaN epilayers grown by molecular beam epitaxy
Castle et al. Report on the 34th IUVSTA workshop ‘XPS: from spectra to results—towards an expert system’
Barreiros et al. Quality assurance of X-ray spectrometry for chemical analysis
JP3965173B2 (ja) 蛍光x線分析装置およびそれに用いるプログラム
JPH0688792A (ja) 汚染元素分析方法
JP2006313132A (ja) 試料分析方法およびx線分析装置
JP4834613B2 (ja) 蛍光x線分析装置およびその方法
US6845147B2 (en) Scatter spectra method for x-ray fluorescent analysis with optical components
Ito et al. Comparison of the Japanese industrial standard and α‐correction methods for X‐ray fluorescence analysis of steels
JP2013205080A (ja) 蛍光x線分析装置
Jonnard et al. Inter‐laboratory comparison of a WDS–EDS quantitative X‐ray microanalysis of a metallic glass
KR101411595B1 (ko) 불용성 전극 내의 이리듐 코팅량 측정방법
CN109563606B (zh) 镀锌钢板的合金化度和/或镀覆附着量测定方法
JP4513425B2 (ja) 亜鉛系めっき鋼板のプレス成形性評価方法
CN113748333A (zh) 荧光x射线分析装置
CN111551579A (zh) 一种利用空白校正确定x射线背景强度的方法
Hajati et al. Three-dimensional X-ray photoelectron tomography on the nanoscale: limits of data processing by principal component analysis

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231126

Address after: No. 168 Chenggong Road, Yushan Town, Kunshan City, Suzhou City, Jiangsu Province, 215300

Patentee after: JIANGSU ELITE INSTRUMENTS Co.,Ltd.

Address before: 710051, No. 1, Changle East Road, Baqiao District, Shaanxi, Xi'an

Patentee before: AIR FORCE ENGINEERING UNIVERSITY OF THE CHINESE PLA