CN105799953B - 一种高精度抗辐射微小卫星三轴磁力矩器及其工作方法 - Google Patents

一种高精度抗辐射微小卫星三轴磁力矩器及其工作方法 Download PDF

Info

Publication number
CN105799953B
CN105799953B CN201610256614.5A CN201610256614A CN105799953B CN 105799953 B CN105799953 B CN 105799953B CN 201610256614 A CN201610256614 A CN 201610256614A CN 105799953 B CN105799953 B CN 105799953B
Authority
CN
China
Prior art keywords
axis
magnetic
processor
bar magnet
magnetic torquer
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
CN201610256614.5A
Other languages
English (en)
Other versions
CN105799953A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201610256614.5A priority Critical patent/CN105799953B/zh
Publication of CN105799953A publication Critical patent/CN105799953A/zh
Application granted granted Critical
Publication of CN105799953B publication Critical patent/CN105799953B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/32Guiding or controlling apparatus, e.g. for attitude control using earth's magnetic field
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/36Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors

Abstract

本发明涉及一种高精度抗辐射微小卫星三轴磁力矩器,包括两个有金属芯的磁棒、一个空心磁力矩器、处理器和磁强计,其中,两个金属芯的磁棒分别形成X轴磁棒和Y轴磁棒,并在其上缠有线圈,空心磁力矩器形成Z轴空心线圈,将X轴磁棒、Y轴磁棒、Z轴空心线圈、处理器和磁强计集中在一块PCB板子上,并形成模块化,三个执行机构由PWM信号独立驱动。本发明还涉及一种高精度抗辐射微小卫星三轴磁力矩器的工作方法。本发明的优点体现在:增加磁力矩器的可靠性;提高磁力矩器的控制精度;具有电流监控与温度测量的功能,磁棒与紧固件一体设计,方便装配,增加了强度。

Description

一种高精度抗辐射微小卫星三轴磁力矩器及其工作方法
技术领域
本发明涉及卫星姿态控制技术领域,具体涉及一种高精度抗辐射微小卫星三轴磁力矩器及其工作方法。
背景技术
磁力矩器是一种卫星姿态控制的执行机构,广泛应用于各种小卫星中,尤其是对于在磁场较强的地球低轨道运行的微小卫星。其工作原理主要是依靠控制线圈中电流的大小和方向来产生相应的磁矩并与地磁场相互作用,从而产生可以控制卫星姿态的磁力矩,进而达到控制卫星姿态的目的。
现有的技术缺点如下:
1.输出精度差;
2.数据处理能力差;
3.缺乏自动温度保护功能;
4.抗辐射性能不佳。
发明内容
本发明以应用于微小卫星的磁力矩器为研究对象,设计了一种新型的磁力矩器,针对目前基于商用器件的磁力矩器普遍存在的处理能力弱、控制精度低、易受空间辐射影响等缺点,采用新型抗辐射处理器,增强处理能力,提高响应速度与控制精度,采用磁强计进行磁场闭环控制,能够有效提高磁力矩器的控制精度。
为实现上述目的,本发明公开了如下技术方案:
一种高精度抗辐射微小卫星三轴磁力矩器,包括两个有金属芯的磁棒、一个空心磁力矩器、处理器和磁强计,其中,两个金属芯的磁棒分别形成X轴磁棒和Y轴磁棒,并在其上缠有线圈,空心磁力矩器形成Z轴空心线圈,将X轴磁棒、Y轴磁棒、Z轴空心线圈、处理器和磁强计集中在一块PCB板子上,并形成模块化,三个执行机构由PWM信号独立驱动。
进一步的,所述Z轴空心线圈设置在PCB板中间,X轴磁棒设置在Z轴空心线圈下侧,Y轴磁棒设置在Z轴空心线圈左侧,处理器和磁强计设置在Z轴空心线圈内侧。
进一步的,所述PCB板子的尺寸为90mm*90mm。
进一步的,所述X轴磁棒和Y轴磁棒选取磁芯最大长度70mm,其中磁芯两端各留出5mm的长度用于固定,磁芯的尺寸确定为Ф4mm×70mm。
进一步的,X轴磁棒和Y轴磁棒上的线圈为具有抗辐照能力的聚酰亚胺包铜线,绝缘封装固定薄膜采用聚酰亚胺薄膜。
进一步的,所述Z轴空心线圈为正方形绕线框,绕线框的内边框长为53mm,外边框长为70mm,线框可绕线的宽度为8.5mm,满足理论计算的约束要求。
本发明还公开了一种高精度抗辐射微小卫星三轴磁力矩器的工作方法,所述磁力矩器采用双磁强计进行磁场反馈,第一磁强计采集磁力矩器工作时三轴线圈附近的磁场强度,第二磁强计安装在卫星其他位置,在磁力矩器关机时采集地磁场强度,两个磁强计采集的数据反馈至处理器,对输出磁矩进行校正。
进一步的,磁力矩器的电路包括处理器,磁强计,时钟,驱动电路,温度传感器,电流监测芯片,处理器是电路中心,时钟、驱动电路直接接入处理器,磁强计通过I2C总线与处理器俩姐,温度传感器、电流检测芯片通过SPI总线与处理器连接,由处理器产生PWM信号,经过驱动电路驱动三轴产生磁力矩,磁力矩器通过I2C总线接收星载计算机传送的控制指令,并将温度与电流信息传送给星载计算机。
进一步的,所述处理器为采用Microsemi生产的SmartFusion2SoC FPGA;所述驱动电路采用TI公司的低压H桥集成芯片drv8837,同时,采用BB公司的宽温仪表用放大器INA337,对流经线圈上的电流信号进行采集,从而反馈给处理器,这样使得磁力矩器板具有检测自身电流的功能;所述温度传感器采用国家半导体公司的低压数字温度传感器LM70,温度传感器实时监测磁力矩器工作温度,并将数据传回处理器进行处理,一旦出现温度过高的情况则关闭磁力矩器,提高磁力矩器的可靠性。
本发明公开的一种高精度抗辐射微小卫星三轴磁力矩器及其工作方法,具有以下有益效果:
1.解决目前磁力矩器控制精度低的问题,采用新型抗辐射处理器,增加磁力矩器的可靠性;
2.采用双磁强计进行输出磁矩的闭环控制,提高磁力矩器的控制精度。
3.具有电流监控与温度测量的功能,在温度过高情况下自动关闭磁力矩器,提高磁力矩器的可靠性。
4.磁棒与紧固件一体设计,方便装配,增加了强度。
附图说明
图1是三轴磁力矩器示意图;
图2是空芯磁力矩器绕线框示意图;
图3是磁力矩器磁场反馈电路;
图4是磁力矩器电路设计整体方案图;
图5是磁力矩器驱动电路;
图6是磁力矩器测温电路;
其中:
1-Y轴磁棒 2-Z轴空心线圈
3-磁强计 4-X轴磁棒
5-处理器
具体实施方式
下面结合实施例并参照附图对本发明作进一步描述。
请参见图1。一种高精度抗辐射微小卫星三轴磁力矩器,包括两个有金属芯的磁棒、一个空心磁力矩器、处理器5和磁强计3,其中,两个金属芯的磁棒分别形成X轴磁棒4和Y轴磁棒1,并在其上缠有线圈,空心磁力矩器形成Z轴空心线圈2,将X轴磁棒4、Y轴磁棒1、Z轴空心线圈2、处理器5和磁强计3集中在一块PCB板子上,并形成模块化,三个执行机构由PWM信号独立驱动。
作为具体实施例,所述Z轴空心线圈2设置在PCB板中间,X轴磁棒4设置在Z轴空心线圈2下侧,Y轴磁棒1设置在Z轴空心线圈2左侧,处理器5和磁强计3设置在Z轴空心线圈2内侧。
作为具体实施例,所述PCB板子的尺寸为90mm*90mm。
作为具体实施例,所述X轴磁棒4和Y轴磁棒1选取磁芯最大长度70mm,其中磁芯两端各留出5mm的长度用于固定,磁芯的尺寸确定为Ф4mm×70mm。
作为具体实施例,考虑到磁力矩器工作的空间环境条件的恶劣性,X轴磁棒4和Y轴磁棒1上的线圈为具有抗辐照能力的聚酰亚胺包铜线,绝缘封装固定薄膜采用聚酰亚胺薄膜。
见图2。作为具体实施例,所述Z轴空心线圈2为正方形绕线框,绕线框的内边框长为53mm,外边框长为70mm,平均边长近似60mm,线框可绕线的宽度为8.5mm,满足理论计算的约束要求。
本磁力矩器是一个三轴磁力矩器,针对微小卫星尺寸、质量和功耗都较小的特点,本磁力矩器的设计方案由两个有金属芯的磁棒和一个空芯磁力矩器组成,将三轴磁力矩器集中在一块90mm×90mm的PCB板子上,节省空间,并形成模块化,三个执行机构由PWM信号独立驱动,如图1所示。该磁力矩器体积很小,可在三个方向提供相当大的磁矩,用于立方星的姿态控制。
见图3。本发明公开了一种高精度抗辐射微小卫星三轴磁力矩器的工作方法,所述磁力矩器采用双磁强计进行磁场反馈,第一磁强计采集磁力矩器工作时三轴线圈附近的磁场强度,第二磁强计安装在卫星其他位置,在磁力矩器关机时采集地磁场强度,两个磁强计采集的数据反馈至处理器,对输出磁矩进行校正。
见图4。作为具体实施例,磁力矩器的电路包括处理器,磁强计,时钟,驱动电路,温度传感器,电流监测芯片,处理器是电路中心,时钟、驱动电路直接接入处理器,磁强计通过I2C总线与处理器俩姐,温度传感器、电流检测芯片通过SPI总线与处理器连接,由处理器产生PWM信号,经过驱动电路驱动三轴产生磁力矩,磁力矩器通过I2C总线接收星载计算机传送的控制指令,并将温度与电流信息传送给星载计算机。
作为具体实施例,所述处理器为采用Microsemi生产的SmartFusion2SoC FPGA;所述驱动电路采用TI公司的低压H桥集成芯片drv8837,同时,采用BB公司的宽温仪表用放大器INA337,对流经线圈上的电流信号进行采集,从而反馈给处理器,这样使得磁力矩器板具有检测自身电流的功能;所述温度传感器采用国家半导体公司的低压数字温度传感器LM70,温度传感器实时监测磁力矩器工作温度,并将数据传回处理器进行处理,一旦出现温度过高的情况则关闭磁力矩器,提高磁力矩器的可靠性。
此方案有以下一些特点:
(1)SmartFusion2处理器处理能力强,且抗辐射性能良好,可以提高磁力矩器在空间工作的可靠性;
(2)采用双磁强计结构,可以测量磁力矩器工作时的磁场强度与磁力矩器关机时卫星当地磁场,从而进行磁力矩器反馈控制,提高控制精度
(3)磁力矩器板根据指令自主产生PWM信号,避免了模拟信号在板间传递;
(4)三轴磁力矩独立驱动。
其中处理器采用Microsemi生产的SmartFusion2SoC FPGA,集成了基于Flash的FPGA构件、一个166MHz ARM CortexTM-M3处理器、先进的安全处理加速器、DSP模块、SRAM、eNVM以及SPI,I2C,MMUART等高性能通信接口。SmartFusion2有丰富的时钟资源,并且最大支持200M外部时钟,相对于目前磁力矩器处理器普遍采用的6M或8M时钟,其控制精度与处理速度将提高数十倍。SmartFusion2对单粒子翻转(SEU)免疫,因此具有很高的可靠性,适合作为航天器处理器。SmartFusion2功耗极低,M2S050系列工作时仅为10mW,待机状态低于1mW。由于其强大的计算能力、高可靠性以及低功耗的特点,SmartFusion2被广泛应用于通信、国防、航空、航天以及医疗等领域。
而驱动电路采用TI公司的低压H桥集成芯片drv8837,该芯片专门用于驱动电机、步进电机或其他封装的绕组,最大可输出1.8A的驱动电流,适合1.8V~11V的电机运行电源电压范围。同时,采用BB公司的宽温仪表用放大器INA337,对流经线圈上的电流信号进行采集,从而反馈给处理器,这样使得磁力矩器板具有检测自身电流的功能。见图5,磁力矩器驱动电路。
见图6磁力矩器测温电路。温度测量采用国家半导体公司的低压数字温度传感器LM70,该芯片使用SPI总线,10位分辨率,精度±2℃,测量范围-55℃~125℃。温度传感器实时监测磁力矩器工作温度,并将数据传回处理器进行处理,一旦出现温度过高的情况则关闭磁力矩器,提高磁力矩器的可靠性。
本发明的关键点在于基于此类架构的三轴磁力矩器抗辐射方法,基于温度反馈的自动保护方法,双磁强计磁场反馈控制方法设计。
由于该磁力矩器X,Y,Z三轴与控制电路在同一个平面内,因此使用该磁力矩器可以有效降低磁力矩器的高度,节省微小卫星的空间,使磁力矩器结构更紧凑;可以控制输出的磁矩大小,并且检测各轴的电流、测量磁力矩器周围的温度,在温度过高时自动关机,从而保障卫星在轨正常运行;抗辐射性能良好,适合空间工作;加入磁场反馈,控制精度提高。
以上所述仅是本发明的优选实施方式,而非对其限制;应当指出,尽管参照上述各实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,其依然可以对上述各实施例所记载的技术方案进行修改,或对其中部分或者全部技术特征进行等同替换;而这些修改和替换,并不使相应的技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (1)

1.一种高精度抗辐射微小卫星三轴磁力矩器,其特征在于,包括两个有金属芯的磁棒、一个空心磁力矩器、处理器和磁强计,其中,两个金属芯的磁棒分别形成X轴磁棒和Y轴磁棒,并在其上缠有线圈,空心磁力矩器形成Z轴空心线圈,将X轴磁棒、Y轴磁棒、Z轴空心线圈、处理器和磁强计集中在一块PCB板子上,并形成模块化,三个执行机构由PWM信号独立驱动,其中:
所述Z轴空心线圈设置在PCB板中间,X轴磁棒设置在Z轴空心线圈下侧,Y轴磁棒设置在Z轴空心线圈左侧,处理器和磁强计设置在Z轴空心线圈内侧;
所述PCB板子的尺寸为90mm*90mm;
所述X轴磁棒和Y轴磁棒选取磁芯最大长度70mm,其中磁芯两端各留出5mm的长度用于固定,磁芯的尺寸确定为Ф4mm×70mm;
X轴磁棒和Y轴磁棒上的线圈为具有抗辐照能力的聚酰亚胺包铜线,绝缘封装固定薄膜采用聚酰亚胺薄膜;
所述Z轴空心线圈为正方形绕线框,绕线框的内边框长为53mm,外边框长为70mm,线框可绕线的宽度为8.5mm,满足理论计算的约束要求;
磁力矩器采用双磁强计进行磁场反馈,第一磁强计采集磁力矩器工作时三轴线圈附近的磁场强度,第二磁强计安装在卫星其他位置,在磁力矩器关机时采集地磁场强度,两个磁强计采集的数据反馈至处理器,对输出磁矩进行校正,其中:
磁力矩器的电路包括处理器,磁强计,时钟,驱动电路,温度传感器,电流监测芯片,处理器是电路中心,时钟、驱动电路直接接入处理器,磁强计通过I2C总线与处理器连接,温度传感器、电流检测芯片通过SPI总线与处理器连接,由处理器产生PWM信号,经过驱动电路驱动三轴产生磁力矩,磁力矩器通过I2C总线接收星载计算机传送的控制指令,并将温度与电流信息传送给星载计算机;
所述处理器为采用Microsemi生产的SmartFusion2SoC FPGA;所述驱动电路采用TI公司的低压H桥集成芯片drv8837,同时,采用BB公司的宽温仪表用放大器INA337,对流经线圈上的电流信号进行采集,从而反馈给处理器,这样使得磁力矩器板具有检测自身电流的功能;所述温度传感器采用美国国家半导体公司的低压数字温度传感器LM70,温度传感器实时监测磁力矩器工作温度,并将数据传回处理器进行处理,一旦出现温度过高的情况则关闭磁力矩器,提高磁力矩器的可靠性。
CN201610256614.5A 2016-04-22 2016-04-22 一种高精度抗辐射微小卫星三轴磁力矩器及其工作方法 Expired - Fee Related CN105799953B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610256614.5A CN105799953B (zh) 2016-04-22 2016-04-22 一种高精度抗辐射微小卫星三轴磁力矩器及其工作方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610256614.5A CN105799953B (zh) 2016-04-22 2016-04-22 一种高精度抗辐射微小卫星三轴磁力矩器及其工作方法

Publications (2)

Publication Number Publication Date
CN105799953A CN105799953A (zh) 2016-07-27
CN105799953B true CN105799953B (zh) 2019-06-04

Family

ID=56457437

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610256614.5A Expired - Fee Related CN105799953B (zh) 2016-04-22 2016-04-22 一种高精度抗辐射微小卫星三轴磁力矩器及其工作方法

Country Status (1)

Country Link
CN (1) CN105799953B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106335656B (zh) * 2016-11-16 2019-11-15 哈尔滨工业大学 一种基于pcb的分布式磁力矩器装置
CN107458625A (zh) * 2017-07-28 2017-12-12 中国科学院长春光学精密机械与物理研究所 应用于微纳卫星的高密度模块化磁力矩器及设计方法
CN109050982A (zh) * 2018-08-30 2018-12-21 西北工业大学 磁力矩器极性在轨检测方法
CN110515023A (zh) * 2019-10-11 2019-11-29 深圳航天东方红海特卫星有限公司 一种微小卫星三轴磁强计校正方法
CN110920936B (zh) * 2019-11-19 2023-03-17 南京璟宏电子有限责任公司 一种应用于微纳卫星的高密度集成化磁力矩器及设计方法
CN112977890B (zh) * 2021-04-06 2024-01-16 湖南揽月机电科技有限公司 一种智能姿控组件中的线圈磁力矩器

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2870609A1 (fr) * 2004-05-18 2005-11-25 Alcatel Sa Dispositif electromagnetique de generation d'une force et d'un couple en vue du positionnement d'un corps
CN1983098A (zh) * 2005-12-14 2007-06-20 上海微小卫星工程中心 主动磁控为主的微小卫星姿态控制方法及系统
US7376496B1 (en) * 2005-12-13 2008-05-20 Lockheed Martin Corporation Spacecraft magnetic momentum control system
CN101554926A (zh) * 2009-05-20 2009-10-14 上海微小卫星工程中心 航天器的姿态控制系统及方法
CN103818566A (zh) * 2014-03-18 2014-05-28 西北工业大学 一种三轴磁力矩器的模块化制作方法
CN104176275A (zh) * 2014-07-16 2014-12-03 北京控制工程研究所 一种使用动量轮与磁力矩器联合的速率阻尼方法
CN104638874A (zh) * 2015-01-05 2015-05-20 南京航空航天大学 单片式三轴磁力矩器

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2870609A1 (fr) * 2004-05-18 2005-11-25 Alcatel Sa Dispositif electromagnetique de generation d'une force et d'un couple en vue du positionnement d'un corps
US7376496B1 (en) * 2005-12-13 2008-05-20 Lockheed Martin Corporation Spacecraft magnetic momentum control system
CN1983098A (zh) * 2005-12-14 2007-06-20 上海微小卫星工程中心 主动磁控为主的微小卫星姿态控制方法及系统
CN101554926A (zh) * 2009-05-20 2009-10-14 上海微小卫星工程中心 航天器的姿态控制系统及方法
CN103818566A (zh) * 2014-03-18 2014-05-28 西北工业大学 一种三轴磁力矩器的模块化制作方法
CN104176275A (zh) * 2014-07-16 2014-12-03 北京控制工程研究所 一种使用动量轮与磁力矩器联合的速率阻尼方法
CN104638874A (zh) * 2015-01-05 2015-05-20 南京航空航天大学 单片式三轴磁力矩器

Also Published As

Publication number Publication date
CN105799953A (zh) 2016-07-27

Similar Documents

Publication Publication Date Title
CN105799953B (zh) 一种高精度抗辐射微小卫星三轴磁力矩器及其工作方法
CN101934863B (zh) 基于磁力矩器和飞轮的卫星姿态全方位控制方法
CN103808321B (zh) 一种基于光源冷备份的三轴一体光纤陀螺惯测装置及补偿和安装方法
Ali et al. Innovative power management, attitude determination and control tile for CubeSat standard NanoSatellites
CN109471433A (zh) 一种基于偏振罗盘的航向与姿态参考系统
CN105799951B (zh) 机电一体化微磁力矩器及磁矩测量方法
CN103337704A (zh) Ku波段卫星通信便携站天线自动对星装置和方法
CN106828981A (zh) 斜飞大惯量耦合卫星的常值干扰力矩补偿方法和系统
CN104638874A (zh) 单片式三轴磁力矩器
RU163618U1 (ru) Электромагнитное устройство ориентации космического аппарата
CN213677194U (zh) 一种航磁探头收放装置及航磁探测装置
US3390848A (en) Magnetic torquing of spin axis stabilization
Daffalla et al. Hardware selection for attitude determination and control subsystem of 1U cube satellite
CN110568387B (zh) 一种基于磁梯度张量的航天器磁矩测试方法
CN205405265U (zh) 一种航拍飞行器的电路
CN204368416U (zh) 一种微型四旋翼飞行器
RU2692741C1 (ru) Устройство контроля параметров углового движения космического аппарата по данным бортовых измерений состояния геомагнитного поля
Cui et al. Analysis of driving torque generated by superconducting motor based on the Meissner effect
CN216209807U (zh) 一种磁棒磁矩的快速测试系统
CN205405148U (zh) 一种具有电压检测功能的航拍飞行器电路
CN213693475U (zh) 一种无人机电源的降压控制电路
Li et al. Control system for a light and small pan-tilt based on multirotor UAV for aerial remote sensing applications
CN211236740U (zh) 一种太阳能无人机的飞行及能源系统控制器
Meili et al. Design of Three Axis Magnetorquer for Microsatellites
CN209928032U (zh) 适用于无人机的磁通门磁力仪

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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: 20190604

Termination date: 20200422

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