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
relative permeability
magnetic core
frequency
circuit
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
Original Assignee
Harbin Institute of Technology
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 filed Critical Harbin Institute of Technology
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

  • Measurement Of Resistance Or Impedance (AREA)
  • Measuring Magnetic Variables (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) 集成电感器的软磁铁氧体磁芯宽频相对磁导率的测量方法
US7479864B2 (en) Total fluid conductivity sensor system and method
Kulkarni et al. Design and analysis of shielded vertically stacked ring resonator as complex permittivity sensor for petroleum oils
de Lima Ávila et al. Improving the performance of an RF resonant cavity water-cut meter using an impedance matching network
CN102004121B (zh) 同时测量陶瓷收缩率和介电常数的装置及方法
Radonić et al. Measurement of complex permeability using short coaxial line reflection method
CN105929346B (zh) 一种非接触式矢量网络高温测试薄膜磁导率的装置及其测量方法
Liu et al. A pair of parallel differential magnetic‐field probes with high measurement accuracy and high electric‐field suppression ratio
CN113009254B (zh) 一种大功率高线性度的电流注入探头
WO2013142388A1 (en) A device for direct microwave measurement of permeability as a function of high dc voltage
Sivaraman et al. Broad band PCB probes for near field measurements
Shenhui et al. An alternative method to determine the initial permeability of ferrite core using network analyzer
Wang et al. A resonant magnetic probe capable of simultaneously measuring two orthogonal magnetic-field components
Cuellar et al. High frequency current probes for common-mode impedance measurements of power converters under operating conditions
Blaz et al. Modeling and characterization of frequency and temperature variation of complex permeability of ferrite LTCC material
Huang et al. Using a single toroidal sample to determine the intrinsic complex permeability and permittivity of Mn–Zn ferrites
Hu et al. Assessment of electromagnetic interference suppression performance of ferrite core loaded power cord
Wang et al. A high-sensitivity composite probe capable of simultaneously measuring electric-and magnetic-field components
Cheng et al. Multifrequency measurements of dielectric properties using a transmission-type overmoded cylindrical cavity
Op't Land et al. Printed circuit board permittivity measurement using waveguide and resonator rings
CN202101949U (zh) 测量低温共烧结陶瓷收缩率和介电常数的装置
Zhang et al. Ultra broadband planar transmission line transformers with ferrite based bandwidth extension
Hao et al. Design of a Sensor Based on CSRR-Derived Structures for Characterizing Permittivity and Permeability Simultaneously
Ahlawat et al. Improved reflection technique for the permittivity measurement of ferrites using slotted coaxial samples
Zhu et al. MATLAB-based multi-parameter optimization of bulk current injection probe

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171117

Termination date: 20180520