CN104849680B - 集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法 - Google Patents

集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法 Download PDF

Info

Publication number
CN104849680B
CN104849680B CN201510259882.8A CN201510259882A CN104849680B CN 104849680 B CN104849680 B CN 104849680B CN 201510259882 A CN201510259882 A CN 201510259882A CN 104849680 B CN104849680 B CN 104849680B
Authority
CN
China
Prior art keywords
integrated inductor
soft ferrite
magnetic permeability
broadband
board
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.)
Expired - Fee Related
Application number
CN201510259882.8A
Other languages
English (en)
Other versions
CN104849680A (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.)
Harbin Institute of Technology Shenzhen
Original Assignee
Harbin Institute of Technology Shenzhen
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 Harbin Institute of Technology Shenzhen filed Critical Harbin Institute of Technology Shenzhen
Priority to CN201510259882.8A priority Critical patent/CN104849680B/zh
Publication of CN104849680A publication Critical patent/CN104849680A/zh
Application granted granted Critical
Publication of CN104849680B publication Critical patent/CN104849680B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measuring Magnetic Variables (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法,属于电磁计量领域。为了解决现有宽频相对磁导率测量方法操作复杂和稳定性差的问题。所述集成电感器包括软磁铁氧体磁芯和板型印制电路;本方法中分别建立了集成电感器的等效电路模型和空载时的等效电路模型,再结合对应的阻抗特性,计算集成电感器等效电感的频率特性;利用几何尺寸测量的方法获取模型中的参数,进一步得到软磁铁氧体磁芯宽频相对磁导率。本发明适用于集成电感器的软磁铁氧体宽频相对磁导率的测量。

Description

集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法
技术领域
本发明属于电磁计量领域。
背景技术
随着半导体器件的发展以及电路拓扑的改进,电力电子变换器的工作频率得以不断提高,然而电磁干扰(EMI)的存在限制了工作效率的进一步提高。在电力电子设备工作的频率范围内,其电磁干扰以传导干扰(0.15~30MHz)为主,抑制传导干扰最有效的措施是加装EMI滤波器。为了在滤波器设计阶段能够准确地预测其插入损耗,为了设计更为有效的电感,需要掌握软磁铁氧体磁芯的高频特性,其中一个重要的电磁参数就是磁芯的相对磁导率。虽然有很多磁性材料可以应用于高频环境下,例如锰锌铁氧体、镍锌铁氧、超微晶软磁磁合金等,但是由于相对磁导率的非线性和高频时电感的寄生效应,仍然很难预测EMI滤波器的插入损耗。通常情况下,生产厂家给出的材料相对磁导率都是低频状态下的初始磁导率,最高频率也不超过3MHz。因此在高频情况下,初始磁导率将不再具有参考价值,必须采用测量的方法重新确定该频率下的相对磁导率,以便更好地掌握材料的非线性磁特性。
目前比较成熟的相对磁导率测量方法按测量原理分主要有两种,分别是网络参数法、谐振腔法。其中,网络参数法是借助对材料散射参数及复反射系数的测量来计算材料的电磁参数,这样传播常数的确定与样品的厚度密切相关,当样品的厚度是测试频率对应的半个波导波长的整数倍时,该方法很不稳定,测量的结果不够精确;谐振腔法是通过测量品质因数及谐振频率来计算材料的电磁参数,对装置的要求较高,波导器件需要较多,装置的微扰影响和负载引起的失谐都会引发实验的误差。因此,提出一种宽频相对磁导率的测量方法,利用集成电感器进行测量,简化实验操作、稳定且具有高度可重复性,使得测量结果具有很高的稳定性和准确性。
发明内容
本发明的目的是为了解决现有宽频相对磁导率测量方法操作复杂和稳定性差的问题,本发明提供一种集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法。
本发明的集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法,
所述集成电感器包括软磁铁氧体磁芯和板型印制电路,所述方法包括如下步骤:
步骤一:根据集成电感器的结构,建立集成电感器的高频等效电路模型一;
步骤二:利用网络分析仪测量集成电感器的阻抗特性,包括幅频特性和相频特性,并得到集成电感器的阻抗特性的曲线;
步骤三:测量网络分析仪空载时的阻抗特性,得到空载时的阻抗特性曲线;
步骤四:根据步骤三得到的空载时的阻抗特性曲线,进行建模,获得空载时的等效电路模型二;
步骤五:将步骤一建立的等效电路模型一与步骤四获得的等效电路模型二相并联,得到等效电路模型三,并获得等效电路模型三的端口阻抗,再结合集成电感器的阻抗特性的曲线反算出待测集成电感器的电感随频率的变化关系;
步骤六:测量待测集成电感器的软磁铁氧体磁芯的有效磁路长度,计算有效磁路面积;
步骤七:根据步骤五获得的电感随频率的变化关系与步骤六获得的有效磁路面积,获得待测集成电感器的软磁铁氧体磁芯宽频相对磁导率。
所述方法还包括:
步骤八:更换集成电感器的板型印制电路,采用步骤一至步骤七进行测量,获得更换板型印制电路后的软磁铁氧体磁芯宽频相对磁导率,利用更换板型印制电路后获得的宽频相对磁导率和步骤七获得的宽频相对磁导率求平均,平均后的宽频相对磁导率为待测集成电感器的宽频相对磁导率;
所述更换集成电感器的板型印制电路包括:更换同一型号的板型印制电路、更换不同绕组匝数的板型印制电路或不同线宽的板型印制电路。
所述方法是基于射频同轴连接器和射频同轴转接线实现,
所述步骤二包括:集成电感器通过射频同轴连接器和射频同轴转接线与网络分析仪相连接,利用网络分析仪测量集成电感器的阻抗特性,包括幅频特性和相频特性,并得到集成电感器的阻抗特性的曲线;
所述步骤三包括:将集成电感器与射频同轴连接器断开,即射频同轴转接线与射频同轴连接器相连,再连接到网络分析仪上,测量射频同轴转接线与射频同轴连接器的阻抗特性,即得到空载时的阻抗特性曲线。
所述步骤七中:获得待测集成电感器的软磁铁氧体磁芯宽频相对磁导率的方法为:
所述宽频相对磁导率
其中,N为板型印制电路的绕组匝数,le为有效磁路长度,Ae为有效磁路面积,L为被测集成电感器的电感。
所述步骤七还包括:对获得的宽频相对磁导率μe进行优化,
优化后的宽频相对磁导率
其中,lg为集成电感器的气隙。
本发明的有益效果在于,建立集成电感器的等效电路模型,利用网络分析仪测量其阻抗特性,同时利用网络分析仪测量射频同轴转接线和射频同轴连接器的阻抗特性,进而建立其等效电路模型。将两个等效电路模型并联,并且与阻抗特性相结合,得到集成电感器的电感L随频率的变化关系。再利用几何测量方法得到磁芯的物理参数,进一步得到磁芯宽频相对磁导率随频率的变化关系。该方法对实验条件和设备要求低,不受外界环境的影响,操作简单,模型中的参数采用几何测量的方法得到,简单易行。采用板型印制电路集成电感器,稳定可靠,对铜导体匝数及位置控制方便,并且具有高度可重复性,进而避免随机性。该平台可以应用于不同形状铁氧体磁芯宽频相对磁导率的测量。
附图说明
图1为具体实施方式一中所述的测量平台的原理示意图。
图2为具体实施方式一中所述的等效电路模型一的电路示意图。
图3为具体实施方式一中集成电感器的增益频率曲线示意图。
图4为具体实施方式一中集成电感器的相角频率曲线示意图。
图5为具体实施方式一中射频同轴转接线与射频同轴连接器幅值频率曲线示意图。
图6为具体实施方式一中射频同轴转接线与射频同轴连接器相角频率曲线示意图。
图7为具体实施方式一中等效电路模型二的电路示意图。
图8为具体实施方式一中所述集成电感器的软磁铁氧体磁芯的原理示意图。
图9为图8的A向剖视图。
图10为具体实施方式一中步骤七获得的有效的宽频相对磁导率的曲线示意图。
具体实施方式
具体实施方式一:结合图1说明本实施方式,本实施方式所述的集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法,所述测量方法基于测量平台实现,所述测量平台包括、射频同轴连接器2、射频同轴转接线3和网络分析仪4;
待测集成电感器1包括软磁铁氧体磁芯和板型印制电路(PCB),
板型印制电路为单面喷镀绕组铜导体,中间镂空,套置在铁氧体磁芯中,集成电感器1的输出端与射频同轴连接器2的一个连接端连接,射频同轴连接器2的另一个连接端经射频同轴转接线3与网络分析仪4连接;
所述测量方法包括:
步骤一、根据待测集成电感器1的结构,建立集成电感器1的高频等效电路模型一:
本实施方式中,首先对集成电感器的结构进行分析,集成电感器结构是一个分布参数结构。铜导体绕组的损耗Rcopper分为直流损耗和交流损耗,其中绕组的直流损耗是由绕线的直流电阻造成的,直流电阻表示为RDC
式中,l为绕组导体长度,γ为绕组导体电导率,S为绕组导体横截面面积。
交流损耗主要是由高频电流的趋肤效应引起的,即导体传送高频交流电流时,电流会集中于导体表面传送,因而呈现出较大的电阻。趋肤深度是指交变电流沿导体表面向导体中心衰减,当衰减到表面电流强度的1/e时所达到的径向深度,其表达式为:
式中,ω为角频率,μ为导体磁导率。由于趋肤效应引起电阻增加,此时应用导体的交流电阻来计算绕组的损耗:
其中,K为趋表系数,不仅与交流电频率有关,而且与材料性质、导线的形状有关。
集成电感器是一个单面PCB结构,故不需要考虑层间分布电容。单面铜绕组各匝之间的分布电容以及绕组与磁芯之间的分布电容可以用一个集总参数Cs来表示,Cs连接在该绕组两端,即:
式中,ε0为真空介电常数,εr为基板的相对介电常数,N为绕组匝数,Ai为每匝绕组导体之间的正对面积,d为绕组匝间距。
磁芯损耗包括磁滞损耗、涡流损耗和剩余损耗,剩余损耗相对于另外两者很小,可以忽略不计。根据铁氧体磁芯的磁滞回线可以发现磁滞损耗Ph与磁感应强度B幅值的平方成正比,与频率f成正比,则有效体积Ve的磁芯产生的磁损可以表示为磁滞损耗和涡流损耗之和,即:Pcore=(af+bf2)VeB2 (5)
式中a、b分别为对曲线线性拟合得到的系数。结合法拉第电磁感应定律,即可得到磁芯损耗的等效电阻:
式中,Ae为磁芯有效磁路面积,le为磁芯有效磁路长度。
这样,集成电感器的高频等效电路模型一建立出来,如图2所示,Rcore为磁芯损耗的等效电阻,Cp为分布电容的集总参数;
步骤二、利用所述测量平台,测量集成电感器1的阻抗特性,包括幅频特性和相频特性,在网络分析仪4上得到集成电感器的阻抗特性的曲线:
首先对网络分析仪(4)进行校准,校准后开始测量,横坐标为频率,纵坐标设为阻抗Z,得到集成电感器的阻抗特性曲线,所述集成电感器的阻抗特性曲线包括图3和图4;
步骤三、将测量平台中集成电感器与射频同轴连接器断开,即射频同轴转接线3与射频同轴连接器2相连,并连接到网络分析仪4上,测量射频同轴转接线3与射频同轴连接器2的阻抗特性,即得到空载时的阻抗特性曲线,所述空载时的阻抗特性曲线包括图5和图6;
步骤四、利用步骤三得到的空载时的阻抗特性曲线,对射频同轴转接线3与射频同轴连接器2进行建模,得到射频同轴转接线3与射频同轴连接器2的等效电路模型二:
基于空载时的阻抗特性曲线,发现射频同轴转接线(3)与射频同轴连接器(2)的阻抗特性呈容性,建立其等效电路模型二,如图7所示,图中的参数R0、C0通过空载时的阻抗特性曲线计算得到,即:
步骤五、将等效电路模型一与等效电路模型二相并联,得到等效电路模型三,计算等效电路模型三的端口阻抗。结合集成电感器的阻抗特性的曲线反算出集成电感器的电感L随频率的变化关系:将得到的两个等效电路模型相并联,得到等效电路模型3,计算等效电路模型3的端口阻抗,即:
Z3=R0+C0//ZAB (8)
式中,ZAB=Cp//(Rcopper+Rcore//L)。这样,结合阻抗集成电感器的阻抗特性的曲线,利用Z3的表达式,反算出集成电感器电感L随频率f的变化关系;
步骤六、利用游标卡尺测量软磁铁氧体磁芯的有效磁路长度le,计算有效磁路面积Ae:
应用几何测量方法,如图8和图9所示,利用游标卡尺测量软磁铁氧体磁芯的有效磁路长度le,并计算有效磁路面积Ae
步骤七、将步骤五与步骤六得到的结果相结合,利用公式计算得到铁氧体磁芯磁导率随频率的变化关系:
由:
反算出有效的宽频相对磁导率μe
有效的宽频相对磁导率的曲线如图10所示。
具体实施方式二:本实施方式是对具体实施方式一所述的集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法的进一步限定,
所述方法还包括:
步骤八:更换集成电感器1的板型印制电路,采用步骤一至步骤七进行测量,获得更换板型印制电路后的软磁铁氧体磁芯宽频相对磁导率,利用更换板型印制电路后获得的宽频相对磁导率和步骤七获得的宽频相对磁导率求平均,平均后的宽频相对磁导率为待测集成电感器(1)的宽频相对磁导率;
所述更换集成电感器1的板型印制电路包括:更换同一型号的板型印制电路、更换不同绕组匝数的板型印制电路或不同线宽的板型印制电路。
步骤八既是对具体实施方式一获得的宽频相对磁导率的一种优化,也可以作为实验验证。
具体实施方式三:本实施方式是对具体实施方式二所述的集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法的进一步限定,
当考虑集成电感器的气隙,宽频相对磁导率μe需进一步优化,优化后的宽频相对磁导率位:
式中,lg为集成电感器的气隙。

Claims (5)

1.一种集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法,所述集成电感器(1)包括软磁铁氧体磁芯和板型印制电路,其特征在于,所述方法包括如下步骤:
步骤一:根据集成电感器(1)的结构,建立集成电感器(1)的高频等效电路模型一;
步骤二:利用网络分析仪(4)测量集成电感器(1)的阻抗特性,包括幅频特性和相频特性,并得到集成电感器的阻抗特性的曲线;
步骤三:测量网络分析仪(4)空载时的阻抗特性,得到空载时的阻抗特性曲线;
步骤四:根据步骤三得到的空载时的阻抗特性曲线,进行建模,获得空载时的等效电路模型二;
步骤五:将步骤一建立的等效电路模型一与步骤四获得的等效电路模型二相并联,得到等效电路模型三,并获得等效电路模型三的端口阻抗,再结合集成电感器的阻抗特性的曲线反算出集成电感器(1)的电感随频率的变化关系;
步骤六:测量集成电感器(1)的软磁铁氧体磁芯的有效磁路长度,计算有效磁路面积;
步骤七:根据步骤五获得的电感随频率的变化关系与步骤六获得的有效磁路面积,获得集成电感器(1)的软磁铁氧体磁芯宽频相对磁导率。
2.根据权利要求1所述的集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法,其特征在于,所述方法还包括:
步骤八:更换集成电感器(1)的板型印制电路,采用步骤一至步骤七进行测量,获得更换板型印制电路后的软磁铁氧体磁芯宽频相对磁导率,利用更换板型印制电路后获得的宽频相对磁导率和步骤七获得的宽频相对磁导率求平均,平均后的宽频相对磁导率为集成电感器(1)的宽频相对磁导率;
所述更换集成电感器(1)的板型印制电路包括:更换同一型号的板型印制电路、更换不同绕组匝数的板型印制电路或不同线宽的板型印制电路。
3.根据权利要求1或2所述的集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法,其特征在于,所述方法是基于射频同轴连接器(2)和射频同轴转接线(3)实现,
所述步骤二包括:集成电感器(1)通过射频同轴连接器(2)和射频同轴转接线(3)与网络分析仪(4)相连接,利用网络分析仪(4)测量集成电感器(1)的阻抗特性,包括幅频特性和相频特性,并得到集成电感器的阻抗特性的曲线;
所述步骤三包括:将集成电感器(1)与射频同轴连接器(2)断开,即射频同轴转接线(3)与射频同轴连接器(2)相连,再连接到网络分析仪(4)上,测量射频同轴转接线(3)与射频同轴连接器(2)的阻抗特性,即得到空载时的阻抗特性曲线。
4.根据权利要求3所述的集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法,其特征在于,所述步骤七中:获得集成电感器(1)的软磁铁氧体磁芯宽频相对磁导率的方法为:
所述宽频相对磁导率
其中,N为板型印制电路的绕组匝数,le为有效磁路长度,Ae为有效磁路面积,L为集成电感器(1)的电感。
5.根据权利要求4所述的集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法,其特征在于,所述步骤七还包括:对获得的宽频相对磁导率μe进行优化,
优化后的宽频相对磁导率
其中,lg为集成电感器(1)的气隙。
CN201510259882.8A 2015-05-20 2015-05-20 集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法 Expired - Fee Related CN104849680B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510259882.8A CN104849680B (zh) 2015-05-20 2015-05-20 集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510259882.8A CN104849680B (zh) 2015-05-20 2015-05-20 集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法

Publications (2)

Publication Number Publication Date
CN104849680A CN104849680A (zh) 2015-08-19
CN104849680B true CN104849680B (zh) 2017-11-17

Family

ID=53849464

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510259882.8A Expired - Fee Related CN104849680B (zh) 2015-05-20 2015-05-20 集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法

Country Status (1)

Country Link
CN (1) CN104849680B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107797083A (zh) * 2017-10-24 2018-03-13 上海交通大学 电磁超材料等效磁导率的测量方法
CN109188324A (zh) * 2018-09-06 2019-01-11 横店集团东磁股份有限公司 一种铁氧体片磁导率的无损全检测方法
CN110954818B (zh) * 2019-12-13 2022-08-09 哈尔滨工业大学 一种中间继电器线圈软故障检测方法
CN113203893B (zh) * 2021-04-29 2022-07-01 福州大学 一种电感器及松耦合变压器的线圈绕组交流电阻提取方法
KR102537967B1 (ko) * 2023-03-02 2023-05-26 김형주 페라이트 코어의 상대 투자율 추출 방법 및 이를 이용한 전자기 시뮬레이션 방법
CN116720469B (zh) * 2023-06-13 2023-12-22 哈尔滨工业大学 基于宽频宏模型的emi滤波器选型方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN203465407U (zh) * 2013-10-12 2014-03-05 北京航空航天大学 基于示波器和信号发生器的软磁材料磁导率测量实验装置
CN104535945A (zh) * 2014-12-01 2015-04-22 国家电网公司 一种饱和电抗器铁心材料的磁特性检测装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8159240B2 (en) * 2009-01-09 2012-04-17 Tdk Corporation Bulk current injection (BCI) probe with multiple, symmetrically spaced feeds

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN203465407U (zh) * 2013-10-12 2014-03-05 北京航空航天大学 基于示波器和信号发生器的软磁材料磁导率测量实验装置
CN104535945A (zh) * 2014-12-01 2015-04-22 国家电网公司 一种饱和电抗器铁心材料的磁特性检测装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种测量软磁铁氧体磁芯磁导率的方法;张近民等;《电测与仪表》;20090631;第46卷(第522期);26-29,67 *

Also Published As

Publication number Publication date
CN104849680A (zh) 2015-08-19

Similar Documents

Publication Publication Date Title
CN104849680B (zh) 集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法
Chou et al. Magnetic near-field probes with high-pass and notch filters for electric field suppression
CN106021811B (zh) 一种磁性材料宽频复数磁导率测定方法
CN106872925B (zh) 一种基于传输线的近场探头空间分辨率驻波校准方法
CN106707210B (zh) 一种基于传输线的近场探头空间分辨率的行波校准方法
CN108982971B (zh) 一种基于矩形腔微扰法测量非磁性材料复介电常数的方法
CN103744039A (zh) 一种测量磁性软磁薄膜复数磁导率的方法
CN112684253A (zh) 一种非接触式负载阻抗测试系统及其工作方法
Jin et al. Improved “Root-Omega” method for transmission-line-based material property extraction for multilayer PCBs
Dosoudil „Determination of permeability from impedance measurement using vector network analyzer,”
Lawson et al. High $ Q $ coil measurement for inductive power transfer
He et al. A wideband tangential electric field probe and a new calibration kit for near-field measurements
CN108872710A (zh) 一种用于测量介电常数的微型双层磁耦合微波传感器
CN105929346B (zh) 一种非接触式矢量网络高温测试薄膜磁导率的装置及其测量方法
US20150035522A1 (en) Device For Direct Microwave Measurement of Permeability as a Function of High DC Voltage
Huang et al. Using a single toroidal sample to determine the intrinsic complex permeability and permittivity of Mn–Zn ferrites
KR101393423B1 (ko) 주파수 종속 복소 유전 상수를 사용하여 주파수 종속 전송선로 파라미터 결정 방법
Illia et al. Finite element method based electromagnetic modeling of three-phase EMI filters
Blaz et al. Modeling and characterization of frequency and temperature variation of complex permeability of ferrite LTCC material
Wu et al. The match game: Numerical de-embedding of field simulation and parameter extraction of circuit models for electromagnetic structure using calibration techniques
JP4743208B2 (ja) 電子部品の電気特性測定方法
Op't Land et al. Printed circuit board permittivity measurement using waveguide and resonator rings
JP2009174951A (ja) 誘電正接評価方法
Hu et al. Assessment of electromagnetic interference suppression performance of ferrite core loaded power cord
Hu et al. Evaluation of ferrite core EMI suppression under realistic working conditions

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171117

Termination date: 20180520

CF01 Termination of patent right due to non-payment of annual fee