CN106777817B - 城市轨道交通电能质量仿真分析***的谐波分析建模方法 - Google Patents

城市轨道交通电能质量仿真分析***的谐波分析建模方法 Download PDF

Info

Publication number
CN106777817B
CN106777817B CN201710038345.XA CN201710038345A CN106777817B CN 106777817 B CN106777817 B CN 106777817B CN 201710038345 A CN201710038345 A CN 201710038345A CN 106777817 B CN106777817 B CN 106777817B
Authority
CN
China
Prior art keywords
module
train
harmonic
urban rail
voltage
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
CN201710038345.XA
Other languages
English (en)
Other versions
CN106777817A (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.)
Qingdao University
Qingdao Power Supply Co of State Grid Shandong Electric Power Co Ltd
Original Assignee
Qingdao University
Qingdao Power Supply Co of State Grid Shandong Electric Power Co Ltd
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 Qingdao University, Qingdao Power Supply Co of State Grid Shandong Electric Power Co Ltd filed Critical Qingdao University
Priority to CN201710038345.XA priority Critical patent/CN106777817B/zh
Publication of CN106777817A publication Critical patent/CN106777817A/zh
Application granted granted Critical
Publication of CN106777817B publication Critical patent/CN106777817B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Electric Propulsion And Braking For Vehicles (AREA)
  • Inverter Devices (AREA)
  • Train Traffic Observation, Control, And Security (AREA)

Abstract

本发明提供了一种城市轨道交通电能质量仿真分析***的谐波分析建模方法，包含如下步骤：步骤1，构建城市轨道交通的供电***仿真模块；步骤2，构建城市轨道交通列车运行状况仿真模块；步骤3，构建FFT谐波分析模块，用于对所述城市轨道交通列车运行状况仿真模块中的录波器模块接收到的波形信息中的谐波进行分析比较。本发明可以模拟城市轨道交通列车的三种工况(牵引、惰行、制动)下的运行状态对电网的谐波影响，从而为获得更好的电源质量提供分析基础。

Description

城市轨道交通电能质量仿真分析***的谐波分析建模方法
技术领域
本发明涉及一种建模方法，尤其是涉及一种城市轨道交通电能质量仿真分析***的谐波分析建模方法。
背景技术
随着经济发展和城市现代化进程的加快，世界各国都面临同样的问题，即城市人口迅猛增长、地域不断扩大，原有的城市地面交通无法满足市民日益增长的出行需求。城市轨道交通正是现代化城市中理想的、能满足市民出行需求的一种公共交通工具。城市轨道交通供电***是为轨道交通运营提供动力能源的一个***，在众多的城市电网用户中是一个重要的用户。它从城市电网接受电源，通过变配电，将合适的电源提供给城市轨道交通各个子***。随着城市轨道交通规模不断扩大，地铁作为城市电网的一级负荷，该***中含有大量的非线性负荷，其投运给电网带来的影响更不容忽视，尤其是对***公共连接点处产生的谐波影响，城市轨道交通接入对电网的谐波影响越来越引起关注。目前，国内外对城轨供电***电能质量的研究主要集中在直流电源整流***所产生的影响，城市轨道交通在运行过程中应分为牵引、惰行、制动三个过程，但是从列车实际投运后，在牵引、惰行、制动三种具体工况下对网侧谐波的分析和研究却很少，这对地铁***谐波的精准化治理带来了不便。
本发明在城市轨道交通电能质量仿真分析***中主要提出不同运行阶段城市轨道交通对电网的谐波影响建模方法。
发明内容
针对现有技术的上述缺陷，本发明提供了一种城市轨道交通电能质量仿真分析***的谐波分析建模方法，其特征在于包含如下步骤：
步骤1，构建城市轨道交通的供电***仿真模块，该供电***仿真模块包括高压母线模块A和高压母线模块B，主变压器模块，所述主变压器模块包括两个110KV/35KV变压器模块，分别用于接收来自高压母线模块A和高压母线模块B的高压电，所述两个110KV/35KV变压器模块中的母线采用单母线分段的接线方式，所述主变压器模块将所述高压母线模块A和高压母线模块B的110KV三相交流电经过降压处理后变成35KV三相交流电送入牵引变电站模块，然后经过牵引变电站模块的降压、整流步骤，得到模拟驱动列车的直流电；
步骤2，构建城市轨道交通列车运行状况仿真模块，该仿真模块包括：直流电源模块，该直流电源模块即为供电模块中的牵引变电站模块，用于向列车运行状况仿真模块输送直流电；IGBT逆变器模块，用于模拟列车逆变器，所述IGBT逆变器模块分别连接至四个并联的感应电动机模块，所述感应电动机模块用于模拟向终端输出机械能；还包括录波器模块，该录波器模块用于接收IGBT逆变器模块的输出电压信息、感应电动机模块输出的电流信息、转速信息、电磁转矩信息；还包括矢量控制模块，所述矢量控制模块的一端与所述IGBT逆变器模块连接，通过向IGBT逆变器模块发送开、断控制脉冲，对列车的运行状态起到控制作用；还包括逻辑控制模块，所述逻辑控制模块的一端与矢量控制模块连接，另一端连接至感应电动机模块的转速信息输出端；所述逻辑控制模块接收所述感应电动机模块输出的转速信息以获得列车速度的信息，通过获得的列车速度信息，确定所述IGBT逆变器模块的开关断脉冲的发送时间；
步骤3，构建FFT谐波分析模块，用于对所述录波器接收到的信息的波形的谐波进行分析比较。
进一步的，所述牵引变电站模块在整流时，采用24脉波整流***进行整流。
进一步的，在步骤2中，所述城市轨道交通列车运行状况仿真模块进行如下操作：从0s开始,通过矢量控制模块,将IGBT逆变器模块打开,控制牵引电机模块转速达到3500转/分钟,负载转矩为300牛*米；列车牵引到70公里/小时；所述录波器模块记录IGBT逆变器模块的电压信号、感应电动机模块的电流信息、转速信号和电磁转矩信号的波形特性信息；
当列车速度达到70公里/小时，通过矢量控制模块,将IGBT逆变器模块关闭,控制牵引电机模块转速下降到2500转/分钟，录波器模块记录IGBT逆变器模块的电压信号、感应电动机模块的电流信息、转速信号和电磁转矩信号的波形特性信息；
当牵引电机模块转速下降到2500转/分钟,即当列车速度达到50公里/小时，逻辑控制模块控制IGBT逆变器模块打开，负载转矩为-300牛*米，录波器模块记录IGBT逆变器模块的电压信号、感应电动机模块的电流信息、转速信号和电磁转矩信号的波形特性信息。
进一步的，所述FFT谐波分析模块，用于对谐波进行分析比较，从而得到列车在各个运行工况下由地铁供电***向***公共连接点处的谐波侵扰情况，包括：分析***公共连接点处的谐波电流值以及谐波电压含有率；分析电压波形畸变率和电流波形畸变率；分析各次谐波对应的电压含有率和谐波电流值。
本发明专利可以对城市轨道交通电能质量进行仿真分析，从而对电流的调整获得直观的分析。
附图说明
图1是本发明的城轨供电***模块仿真图；
图2是本发明的城轨列车三种运行状况进行仿真的模块图；
图3是本发明的城轨列车接入城轨供电***后对三种工况下的运行状态进行跟踪后的仿真结果示意图。
具体实施方式
本发明的城市轨道交通电能质量仿真分析***的谐波分析建模方法包括对城轨供电***的供电模型进行建模以及对列车的运行状况进行建模以及对测试结果进行建模。
图1是城轨供电***的供电***仿真模块，包括高压母线模块A和高压母线模块B，主变压器模块，主变压器模块包括两个110KV/35KV变压器模块，分别用于接收来自高压母线模块A和高压母线模块B的高压电源。主变压器模块采用从不同高压母线模块引入供电线路的双回路供电的方式，两个110KV/35KV变压器模块中的母线则采用单母线分段的接线方式，这样提高了供电的可靠性。主变压器模块将城轨供电***的供电模块的110KV三相交流电经过降压处理后变成35KV三相交流电送入牵引变电站模块，然后经过牵引变电站模块的降压、整流等步骤，得到能够驱动列车的直流电。
在整流时，采用了24脉波整流***，即一个周期内输出的脉波数为24，其目的是减小直流电压的纹波系数并降低低次谐波的出现，从而有效降低了谐波源对城市电网的影响。
在列车运营的过程中，共分为牵引、惰行和制动三种运行工况。牵引时，列车通过位于车顶的受电弓将来自架空线路的且合乎城市轨道牵引电压的直流电压输送给牵引逆变器，从而得到用于驱动交流电机的三相交流电。列车内部的四个异步电动机则通过并联运行，将来自逆变部分的交流电转化成供机车运行的机械能，并通过位于编组列车的转向架为其四个轴提供转向力，进而带动整个列车组迅速地完成牵引过程。
列车完成牵引后，通过开关控制关闭列车的牵引逆变器，列车进入惰行状态。由于逆变器关断，列车的电机未施加电源，没有能量的供给，故整个惰行过程没有功率的消耗，列车只凭惯性的作用继续前进。
制动时，逆变器打开，电机转矩反向，进入发电机状态，将产生的三相交流电通过逆变器转变为直流电能。
基于以上原理，图2构建了本发明地铁列车三种运行状况的仿真模块。
该仿真模块包括直流电源模块，该直流电源模块即为仿真供电模块中的牵引变电站模块，其向列车仿真模块输送直流电，IGBT逆变器模块用于模拟列车逆变器，所述IGBT逆变器模块分别连接至四个并联的感应电动机模块，所述感应电动机模块用于模拟向终端输出机械能。还包括录波器模块，该录波器模块用于接收IGBT逆变器模块的输出电压信息、感应电动机模块输出的电流信息、转速信息、电磁转矩信息。
还包括矢量控制模块，感应电机模块作为列车驱动力的源头，且电机本质上为非线性、多变量、强耦合、参数时变、大干扰的复杂对象，通过矢量控制模块能够较好地实现电机的有效控制。该矢量控制模块的输出端一端与IGBT逆变器模块连接，一端与逻辑控制模块连接。矢量控制模块通过向IGBT逆变器模块发送开、断控制脉冲，对列车的运行状态起到控制作用。逻辑控制***的一端与IGBT逆变器模块连接另一端连接至感应电动机模块的转速信息输出端。该转速信息输出端同时还连接到录波器模块。
逻辑控制模块接收所述感应电动机模块输出的转速信息以获得列车速度的信息，通过获得的列车速度信息，决定逆变器模块的开关断脉冲的发送时间，具体步骤如下：
第一阶段为牵引过程,从0s开始,通过矢量控制模块,将IGBT逆变器模块打开,控制牵引电机转速达到3500转/分钟,负载转矩为300牛*米；列车牵引到70公里/小时；此时录波器模块记录IGBT逆变器模块的电压信号、感应电动机模块的电流信息、转速信号和电磁转矩信号的波形特性信息；
第二阶段为惰行过程，当列车速度达到70公里/小时，通过矢量控制模块,将IGBT逆变器模块关闭,控制牵引电机模块转速下降到2500转/分钟。此时录波器模块记录IGBT逆变器模块的电压信号、感应电动机模块的电流信息、转速信号和电磁转矩信号的波形特性信息；
第三阶段为制动过程，当牵引电机模块转速下降到2500转/分钟,即当列车速度达到50公里/小时，逻辑控制模块控制IGBT逆变器模块打开，负载转矩为-300牛*米。此时录波器模块记录IGBT逆变器模块的电压信号、感应电动机模块的电流信息、转速信号和电磁转矩信号的波形特性信息；
这就完成了对列车在三种工况(牵引、惰行、制动)下的运行状态的跟踪，录波器记录的仿真结果如图3所示。
从图3中电机转矩(electromagnetic torque Te)及转速(N)可以看出，随着时间(time)的前进，电机的启停一次先后经过了牵引过程、惰行过程以及制动过程。列车首先进入的是牵引工况，在异步电动机提供的牵引力的作用下，带动转子逐步升速直至最高转速，此时，列车也相应达到了行驶的最高时速，约60～80公里/小时；然后列车的牵引逆变器关闭，进入惰行工况，列车凭借惯性前行，只受到来自空气和地面的阻力的作用，速度逐渐下降，当转速下降到设定值的时候，重新开启列车的逆变器；最后，列车进入制动工况，此时电机作发电机运行，在制动过程中，速度急速降落，直到列车车速降为零。
还包括FFT谐波分析模块，用于对谐波进行分析比较，可以得到列车在各个运行工况下由地铁供电***向***公共连接点处的谐波侵扰情况。
包括分析***公共连接点处的谐波电流值以及谐波电压含有率：可以发现***公共连接点处的谐波电流值以及谐波电压含有率都是在第23次、第25次谐波处数值最大，在第47次、第49次谐波处数值较大，该结果与24脉波整流的方式相吻合，故可以得出在24脉波整流方式下，***公共连接点处的特征次谐波值较大，而非特征次谐波处的谐波值较小的结论。
分析电压波形畸变率和电流波形畸变率：从三种运行工况的角度分析，牵引工况下电压波形畸变率和电流波形畸变率相对较大，制动过程的对应值次之，惰行过程所产生谐波畸变率最小。
分析各次谐波对应的电压含有率和谐波电流值：牵引工况下各次谐波对应的电压含有率和谐波电流值最大，制动过程的对应值次之，惰行过程数值最小。
尽管已经结合实施例对本发明进行了详细地描述，但是本领域技术人员应当理解地是，在不背离本发明精神和实质下的各种修正、形变都是允许的，它们都落入本发明权利要求的保护范围之中。

Claims (3)

1.一种城市轨道交通电能质量仿真分析***的谐波分析建模方法，其特征在于包含如下步骤：
步骤1，构建城市轨道交通的供电***仿真模块，该供电***仿真模块包括高压母线模块A和高压母线模块B，主变压器模块，所述主变压器模块包括两个110KV/35KV变压器模块，分别用于接收来自高压母线模块A和高压母线模块B的高压电，所述两个110KV/35KV变压器模块中的母线采用单母线分段的接线方式，所述主变压器模块将所述高压母线模块A和高压母线模块B的110KV三相交流电经过降压处理后变成35KV三相交流电送入牵引变电站模块，然后经过牵引变电站模块的降压、整流步骤，得到模拟驱动列车的直流电；
步骤2，构建城市轨道交通列车运行状况仿真模块，该仿真模块包括：直流电源模块，该直流电源模块即为供电模块中的牵引变电站模块，用于向列车运行状况仿真模块输送直流电；IGBT逆变器模块，用于模拟列车逆变器，所述IGBT逆变器模块分别连接至四个并联的感应电动机模块，所述感应电动机模块用于模拟向终端输出机械能；还包括录波器模块，该录波器模块用于接收IGBT逆变器模块的输出电压信息、感应电动机模块输出的电流信息、转速信息、电磁转矩信息；还包括矢量控制模块，所述矢量控制模块的一端与所述IGBT逆变器模块连接，通过向IGBT逆变器模块发送开、断控制脉冲，对列车的运行状态起到控制作用；还包括逻辑控制模块，所述逻辑控制模块的一端与矢量控制模块连接，另一端连接至感应电动机模块的转速信息输出端；所述逻辑控制模块接收所述感应电动机模块输出的转速信息以获得列车速度的信息，通过获得的列车速度信息，确定所述IGBT逆变器模块的开关断脉冲的发送时间；其中，所述城市轨道交通列车运行状况仿真模块进行如下操作：从0s开始,通过矢量控制模块,将IGBT逆变器模块打开,控制牵引电机模块转速达到3500转/分钟,负载转矩为300牛*米；列车牵引到70公里/小时；所述录波器模块记录IGBT逆变器模块的电压信号、感应电动机模块的电流信息、转速信号和电磁转矩信号的波形特性信息；
当列车速度达到70公里/小时，通过矢量控制模块,将IGBT逆变器模块关闭,控制牵引电机模块转速下降到2500转/分钟，录波器模块记录IGBT逆变器模块的电压信号、感应电动机模块的电流信息、转速信号和电磁转矩信号的波形特性信息；
当牵引电机模块转速下降到2500转/分钟,即当列车速度达到50公里/小时，逻辑控制模块控制IGBT逆变器模块打开，负载转矩为-300牛*米，录波器模块记录IGBT逆变器模块的电压信号、感应电动机模块的电流信息、转速信号和电磁转矩信号的波形特性信息；
步骤3，构建FFT谐波分析模块，用于对所述录波器接收到的信息的波形的谐波进行分析比较。
2.根据权利要求1所述城市轨道交通电能质量仿真分析***的谐波分析建模方法，其特征在于，所述牵引变电站模块在整流时，采用24脉波整流***进行整流。
3.根据权利要求2所述的城市轨道交通电能质量仿真分析***的谐波分析建模方法，其特征在于，所述FFT谐波分析模块，用于对谐波进行分析比较，从而得到列车在各个运行工况下由地铁供电***向***公共连接点处的谐波侵扰情况，包括：分析***公共连接点处的谐波电流值以及谐波电压含有率；分析电压波形畸变率和电流波形畸变率；分析各次谐波对应的电压含有率和谐波电流值。
CN201710038345.XA 2017-01-19 2017-01-19 城市轨道交通电能质量仿真分析***的谐波分析建模方法 Expired - Fee Related CN106777817B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710038345.XA CN106777817B (zh) 2017-01-19 2017-01-19 城市轨道交通电能质量仿真分析***的谐波分析建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710038345.XA CN106777817B (zh) 2017-01-19 2017-01-19 城市轨道交通电能质量仿真分析***的谐波分析建模方法

Publications (2)

Publication Number Publication Date
CN106777817A CN106777817A (zh) 2017-05-31
CN106777817B true CN106777817B (zh) 2020-04-07

Family

ID=58943411

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710038345.XA Expired - Fee Related CN106777817B (zh) 2017-01-19 2017-01-19 城市轨道交通电能质量仿真分析***的谐波分析建模方法

Country Status (1)

Country Link
CN (1) CN106777817B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019200514A1 (zh) * 2018-04-16 2019-10-24 北京千驷驭电气有限公司 城市轨道交通牵引供电动模仿真测试***
CN108279350B (zh) * 2018-04-16 2020-03-17 北京千驷驭电气有限公司 城市轨道交通牵引供电动模仿真测试***
CN113486541B (zh) * 2021-09-08 2021-11-23 中铁电气化勘测设计研究院有限公司 基于列车安全距离和电压的城轨直流供电***仿真方法
CN114325021A (zh) * 2021-12-29 2022-04-12 国网上海市电力公司 一种考虑高铁谐波特性的电费计算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103793563A (zh) * 2014-01-24 2014-05-14 广西电网公司电力科学研究院 一种风电场等效模拟模块与交直机车及其供电***
CN104281736A (zh) * 2014-08-01 2015-01-14 国家电网公司 一种城市轨道交通牵引供配电***供电电缆的选型方法
CN104346525A (zh) * 2014-09-26 2015-02-11 广东电网有限责任公司电力科学研究院 城轨交通供电***牵引变电所多列车次累计功率计算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10134883A1 (de) * 2001-07-18 2003-01-30 Abb Research Ltd Verfahren und Vorrichtung zur drehzahlstellbaren leistungselektronischen Regelung einer getriebelosen Windkraftanlage

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103793563A (zh) * 2014-01-24 2014-05-14 广西电网公司电力科学研究院 一种风电场等效模拟模块与交直机车及其供电***
CN104281736A (zh) * 2014-08-01 2015-01-14 国家电网公司 一种城市轨道交通牵引供配电***供电电缆的选型方法
CN104346525A (zh) * 2014-09-26 2015-02-11 广东电网有限责任公司电力科学研究院 城轨交通供电***牵引变电所多列车次累计功率计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于城市轨道交通运行图的供配电***谐波预评估仿真;朱明星 等;《城市轨道交通研究》;20151231(第6期);56-61 *

Also Published As

Publication number Publication date
CN106777817A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106777817B (zh) 城市轨道交通电能质量仿真分析***的谐波分析建模方法
CN102841251B (zh) 一种微网中电动汽车充电站谐波检测法
CN206531936U (zh) 电动汽车电机驱动***测试设备
CN112693327B (zh) 一种降低非工作损耗的新能源永磁整车控制子***、方法及车辆
CN106385215A (zh) 一种交流电机变开关频率pwm转矩脉动控制方法
CN101488728A (zh) 无电解电容器的电梯驱动逆变器***
CN111505963A (zh) 一种磁悬浮列车地面供电***的半实物仿真***及方法
CN113253703B (zh) 一种高速磁浮列车模拟测试***及测试方法
CN101244699A (zh) 城市轨道交通车辆运行控制实验***
CN116742660B (zh) 一种考虑高铁再生制动的全周期工况负序分析方法及***
CN109725222A (zh) 一种牵引箱整机智能负载测试设备及其方法
CN201177914Y (zh) 一种轨道交通车辆运行控制实验装置
CN104428158A (zh) 交流电车的控制装置
Chang et al. Modeling and simulation of DC transit system with VSI-fed induction motor driven train using PSB/MATLAB
CN101997474B (zh) 一种用于汽车的电动机控制器
CN113486541B (zh) 基于列车安全距离和电压的城轨直流供电***仿真方法
CN103856134A (zh) 用于电动汽车的异步电机控制方法及其***
CN113110109B (zh) 一种高速磁浮列车电机模拟测试装置及***
CN111707935B (zh) 电机模拟***
Kovudhikulrungsri et al. Real time simulator for railway traction and auxiliary power unit control applications
CN206583983U (zh) 一种再生能量逆变回馈的小功率验证平台
CN207472582U (zh) 一种列车滚动试验台的能量回收装置
Wang et al. Hybrid modeling and applications of virtual metro systems
Lu et al. Simulation software for CRH2 and CRH3 traction driver systems based on SIMULINK and VC
CN116560261B (zh) 一种电流源型电机模拟器及其控制方法和***

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
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: 20200407

Termination date: 20210119