CN106017218A - 一种抗风干扰的飞行器攻角指令补偿方法及装置 - Google Patents

一种抗风干扰的飞行器攻角指令补偿方法及装置 Download PDF

Info

Publication number
CN106017218A
CN106017218A CN201610331421.1A CN201610331421A CN106017218A CN 106017218 A CN106017218 A CN 106017218A CN 201610331421 A CN201610331421 A CN 201610331421A CN 106017218 A CN106017218 A CN 106017218A
Authority
CN
China
Prior art keywords
aircraft
attack
angle
instruction
guidance
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.)
Granted
Application number
CN201610331421.1A
Other languages
English (en)
Other versions
CN106017218B (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.)
Beijing Aerospace Automatic Control Research Institute
Original Assignee
Beijing Aerospace Automatic Control Research Institute
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 Beijing Aerospace Automatic Control Research Institute filed Critical Beijing Aerospace Automatic Control Research Institute
Priority to CN201610331421.1A priority Critical patent/CN106017218B/zh
Publication of CN106017218A publication Critical patent/CN106017218A/zh
Application granted granted Critical
Publication of CN106017218B publication Critical patent/CN106017218B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F41WEAPONS
    • F41GWEAPON SIGHTS; AIMING
    • F41G3/00Aiming or laying means
    • F41G3/22Aiming or laying means for vehicle-borne armament, e.g. on aircraft
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
    • G05D1/0816Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
    • G05D1/0825Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using mathematical models

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • General Engineering & Computer Science (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明提供一种抗风干扰的飞行器攻角指令补偿方法及装置,涉及飞行器控制应用技术领域,用于解决由于风干扰等因素影响,导致制导精度降低的问题。该方法包括获取飞行器的加速度平滑值;根据飞行器制导产生的纵向过载力和侧向过载力,获取基本攻角指令;根据飞行器制导产生的纵向过载力、侧向过载力及加速度平滑值,获得攻角补偿指令;根据所述基本攻角指令及攻角补偿指令计算实际的制导攻角指令。上述方案,能在有风条件下通过对飞行器制导输出的攻角指令进行补偿,克服了风对飞行器制导的干扰问题,从而达到了提高制导精度的目的。

Description

一种抗风干扰的飞行器攻角指令补偿方法及装置
技术领域
本发明涉及飞行器控制应用技术领域,特别是指一种抗风干扰的飞行器攻角指令补偿方法及装置。
背景技术
目前,对于导弹、滑翔飞行器等飞行器,如果制导采用攻角控制方案,在飞行器实际飞行过程中,会由于风干扰等因素影响,导致制导精度降低,不能满足目标打击、定点着陆等任务要求。为解决由于风干扰等因素影响,导致制导精度降低的问题,近来研究人员提出利用过载控制方案解决该问题,虽然该方案具有抗风干扰的能力,但一些特定飞行器或特定工作条件对攻角参数会有严格要求的情况,如采用吸气式发动机的飞行器发动机工作时要求攻角控制在特定范围内,又如飞行器在飞行末段对姿态有严格约束等问题,过载控制方案就不在适用。
故急需一种抗风干扰的飞行器控制方法及装置,能具有抗风干扰的能力,从而提高飞行器制导的精度。
发明内容
为了解决由于风干扰等因素影响,导致制导精度降低的问题,本发明提供一种抗风干扰的飞行器攻角指令补偿方法及装置,能在有风条件下通过对飞行器制导输出的攻角指令进行补偿,克服了风对飞行器制导的干扰问题,从而达到了提高制导精度的目的。
本发明提供的一种抗风干扰的飞行器攻角指令补偿方法,包括如下步骤:
获取飞行器的加速度平滑值;
根据飞行器制导产生的纵向过载力和侧向过载力,获取基本攻角指令;
根据飞行器制导产生的纵向过载力、侧向过载力及加速度平滑值,获得攻角补偿指令;
根据所述基本攻角指令及攻角补偿指令计算实际的制导攻角指令。
其中,所述获取飞行器的加速度平滑值,包括:
测量n个周期所述飞行器的加速度值
对所述n个周期所述飞行器的加速度值按照下式进行平滑处理,获得所述飞行器的加速度平滑值
W ‾ y 1 = 1 n Σ i = 1 n W y 1 i
其中,表示第i个周期所述飞行器的加速度值。
其中,所述根据飞行器制导产生的纵向过载力、侧向过载力及加速度平滑值,获得攻角补偿指令Dαcx的方法为:
Δα c x = F y c x 2 + F z c x 2 - W ‾ y 1 · m q · S · b
其中,为飞行器制导产生的纵向和侧向过载力,为飞行器的加速度平滑值,m为飞行器质量,q为飞行器动压,S为飞行器参考面积,b为飞行器气动系数。
其中,所述飞行器气动系数根据所述飞行器的气动特性拟合得到。
其中,所述根据所述基本攻角指令及攻角补偿指令计算实际的制导攻角指令,包括:将所述基本攻角指令和所述攻角补偿指令之和作为所述实际的制导攻角指令。
本发明还提供了一种抗风干扰的飞行器攻角指令补偿装置,包括加速度平滑模块、基本攻角指令获取模块、攻角补偿指令获取模块和制导攻角指令获取模块;
所述加速度平滑模块,用于获取飞行器的加速度平滑值,并将所述飞行器的加速度平滑值发送给所述攻角补偿指令获取模块;
所述基本攻角指令获取模块,用于根据飞行器制导产生的纵向过载力和侧向过载力获取基本攻角指令,并将所述基本攻角指令发送给所述制导攻角指令获取模块;
所述攻角补偿指令获取模块,用于根据飞行器制导产生的纵向过载力、侧向过载力及所述加速度平滑模块提供的加速度平滑值,获得攻角补偿指令并发送给所述制导攻角指令获取模块;
所述制导攻角指令获取模块,用于根据收到的所述基本攻角指令及攻角补偿指令,计算实际的制导攻角指令。
其中,所述加速度平滑模块包括:
惯性器件,用于测量n个周期所述飞行器的加速度值并将测量的每个周期的飞行器的加速度值输出给平滑单元;
平滑单元,用于将收到的n个周期所述飞行器的加速度值按照式进行平滑处理,获得所述飞行器的加速度平滑值并将所述加速度平滑值发送给所述攻角补偿指令获取模块;其中,表示第i个周期所述飞行器的加速度值。
其中,所述攻角补偿指令获取模块根据下式获得攻角补偿指令Dαcx
Δα c x = F y c x 2 + F z c x 2 - W ‾ y 1 · m q · S · b
其中,为飞行器制导产生的纵向和侧向过载力,为飞行器的加速度平滑值,m为飞行器质量,q为飞行器动压,S为飞行器参考面积,b为飞行器气动系数。
本发明的上述技术方案的有益效果如下:
上述方案中,基于飞行器上惯性器件测量得到的加速度信息,结合制导产生的纵向过载力和侧向过载力,能够很方便的得到由风干扰引起的需补偿的攻角指令。通过该方案可有效消除风干扰对制导的影响,提高制导精度,可直接用于采用攻角控制的各类飞行器制导方案设计中。
与采用真实攻角估计方法相比,该方案不会受到攻角估计误差的影响,因为攻角估计精度除风影响外,还会受到方向力和大气密度误差的影响。
附图说明
图1为本发明实施例提供的一种抗风干扰的飞行器攻角指令补偿方法流程示意图;
图2为获取飞行器的加速度平滑值方法的流程示意图;
图3为本发明提供的一种抗风干扰的飞行器攻角指令补偿装置结构示意图;
图4为加速度平滑模块的结构示意图。
[附图标记说明]
31、加速度平滑模块;
32、基本攻角指令获取模块;
33、攻角补偿指令获取模块;
34、制导攻角指令获取模块;
311、惯性器件;
312、平滑单元。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
图1为本发明实施例提供的一种抗风干扰的飞行器攻角指令补偿方法流程示意图,包括如下步骤S101-S104:
步骤S101:获取飞行器的加速度平滑值;
步骤S102:根据飞行器制导产生的纵向过载力和侧向过载力,获取基本攻角指令;
步骤S103:根据飞行器制导产生的纵向过载力、侧向过载力及加速度平滑值,获得攻角补偿指令;
步骤S104:根据基本攻角指令及攻角补偿指令计算实际的制导攻角指令。
上述方法,可以应用于飞行器在有风条件下通过对制导输出的攻角指令进行补偿来提高制导精度,通过该方法,能够根据飞行器上惯性器件测得的加速度信息间接得到由风干扰引起的攻角,从而基于该攻角对制导指令进行补偿,达到提高制导精度目的。
在一个实施例中,如图2所示,步骤S101可包括如下步骤S201-S202:
步骤S201:测量n个周期飞行器的加速度值
此步骤中,可以通过飞行器上装有的惯性器件测量得到飞行器的加速度信息,实现方式较简单。另外n为平滑周期,n可根据需要选取合适的值。
步骤S202:对n个周期飞行器的加速度值按照下式进行平滑处理,获得飞行器的加速度平滑值
W ‾ y 1 = 1 n Σ i = 1 n W y 1 i - - - ( 1 )
在公式(1)中,表示第i个周期飞行器的加速度值。
在一个实施例中,步骤S103中采用以下公式(2)获得攻角补偿指令Dαcx
Δα c x = F y c x 2 + F z c x 2 - W ‾ y 1 · m q · S · b - - - ( 2 )
在公式(2)中,为飞行器制导产生的纵向和侧向过载力,为飞行器的加速度平滑值,m为飞行器质量,q为飞行器动压,S为飞行器参考面积,b为飞行器气动系数。
其中,飞行器气动系数b可以根据飞行器的气动特性拟合得到。
在一个实施例中,步骤S104中将获得的基本攻角指令和攻角补偿指令之和作为实际的制导攻角指令,即:
αcx=αcx0+Dαcx (3)
在公式(3)中,αcx为实际的制导攻角指令,αcx0为基本攻角指令,Δαcx为攻角补偿指令。
对应于上述实施例提供的一种抗风干扰的飞行器攻角指令补偿方法,本发明实施例还提供一种抗风干扰的飞行器攻角指令补偿装置,如图3所示,包括:加速度平滑模块31、基本攻角指令获取模块32、攻角补偿指令获取模块33和制导攻角指令获取模块34;
加速度平滑模块31,用于获取飞行器的加速度平滑值,并将飞行器的加速度平滑值发送给攻角补偿指令获取模块33;
基本攻角指令获取模块32,用于根据飞行器制导产生的纵向过载力和侧向过载力获取基本攻角指令,并将基本攻角指令发送给制导攻角指令获取模块34;
攻角补偿指令获取模块33,用于根据飞行器制导产生的纵向过载力、侧向过载力及加速度平滑模块提供的加速度平滑值,获得攻角补偿指令并发送给制导攻角指令获取模块34;
制导攻角指令获取模块34,用于根据收到的基本攻角指令及攻角补偿指令,计算实际的制导攻角指令。
在一个实施例中,如图4所示,加速度平滑模块31包括:
惯性器件311,用于测量n个周期飞行器的加速度值并将测量的每个周期的飞行器的加速度值输出给平滑单元312;
平滑单元312,用于将收到的n个周期飞行器的加速度值按照式进行平滑处理,获得飞行器的加速度平滑值并发送给攻角补偿指令获取模块33。其中,表示第i个周期飞行器的加速度值。
上述实施例提供的装置中,攻角补偿指令获取模块33可根据下式获得攻角补偿指令Δαcx
Δα c x = F y c x 2 + F z c x 2 - W ‾ y 1 · m q · S · b
其中,为飞行器制导产生的纵向和侧向过载力,为飞行器的加速度平滑值,m为飞行器质量,q为飞行器动压,S为飞行器参考面积,b为飞行器气动系数。
上述方案中,基于飞行器上惯性器件测量得到的加速度信息,结合制导产生的纵向过载力和侧向过载力,能够很方便的得到由风干扰引起的需补偿的攻角指令。通过该方案可有效消除风干扰对制导的影响,提高制导精度,可直接用于采用攻角控制的各类飞行器制导方案设计中。
与采用真实攻角估计方法相比,该方案不会受到攻角估计误差的影响,因为攻角估计精度除风影响外,还会受到方向力和大气密度误差的影响。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (8)

1.一种抗风干扰的飞行器攻角指令补偿方法,其特征在于,包括如下步骤:
获取飞行器的加速度平滑值;
根据飞行器制导产生的纵向过载力和侧向过载力,获取基本攻角指令;
根据飞行器制导产生的纵向过载力、侧向过载力及加速度平滑值,获得攻角补偿指令;
根据所述基本攻角指令及攻角补偿指令计算实际的制导攻角指令。
2.如权利要求1所述的抗风干扰的飞行器攻角指令补偿方法,其特征在于,所述获取飞行器的加速度平滑值,包括:
测量n个周期所述飞行器的加速度值
对所述n个周期所述飞行器的加速度值按照下式进行平滑处理,获得所述飞行器的加速度平滑值
W ‾ y 1 = 1 n Σ i = 1 n W y 1 i
其中,表示第i个周期所述飞行器的加速度值。
3.如权利要求1或2所述的抗风干扰的飞行器攻角指令补偿方法,其特征在于,所述根据飞行器制导产生的纵向过载力、侧向过载力及加速度平滑值,获得攻角补偿指令Δαcx的方法为:
Δα c x = F y c x 2 + F z c x 2 - W ‾ y 1 · m q · S · b
其中,为飞行器制导产生的纵向和侧向过载力,为飞行器的加速度平滑值,m为飞行器质量,q为飞行器动压,S为飞行器参考面积,b为飞行器气动系数。
4.如权利要求3所述的抗风干扰的飞行器攻角指令补偿方法,其特征在于,所述飞行器气动系数根据所述飞行器的气动特性拟合得到。
5.如权利要求3所述的抗风干扰的飞行器攻角指令补偿方法,其特征在于,所述根据所述基本攻角指令及攻角补偿指令计算实际的制导攻角指令,包括:将所述基本攻角指令和所述攻角补偿指令之和作为所述实际的制导攻角指令。
6.一种抗风干扰的飞行器攻角指令补偿装置,其特征在于,包括加速度平滑模块、基本攻角指令获取模块、攻角补偿指令获取模块和制导攻角指令获取模块;
所述加速度平滑模块,用于获取飞行器的加速度平滑值,并将所述飞行器的加速度平滑值发送给所述攻角补偿指令获取模块;
所述基本攻角指令获取模块,用于根据飞行器制导产生的纵向过载力和侧向过载力获取基本攻角指令,并将所述基本攻角指令发送给所述制导攻角指令获取模块;
所述攻角补偿指令获取模块,用于根据飞行器制导产生的纵向过载力、侧向过载力及所述加速度平滑模块提供的加速度平滑值,获得攻角补偿指令并发送给所述制导攻角指令获取模块;
所述制导攻角指令获取模块,用于根据收到的所述基本攻角指令及攻角补偿指令,计算实际的制导攻角指令。
7.如权利要求6所述的抗风干扰的飞行器攻角指令补偿装置,其特征在于,所述加速度平滑模块包括:
惯性器件,用于测量n个周期所述飞行器的加速度值并将测量的每个周期的飞行器的加速度值输出给平滑单元;
平滑单元,用于将收到的n个周期所述飞行器的加速度值按照式进行平滑处理,获得所述飞行器的加速度平滑值并将所述加速度平滑值发送给所述攻角补偿指令获取模块;其中,表示第i个周期所述飞行器的加速度值。
8.如权利要求6或7所述的抗风干扰的飞行器攻角指令补偿装置,其特征在于,所述攻角补偿指令获取模块根据下式获得攻角补偿指令Δαcx
Δα c x = F y c x 2 + F z c x 2 - W ‾ y 1 · m q · S · b
其中,为飞行器制导产生的纵向和侧向过载力,为飞行器的加速度平滑值,m为飞行器质量,q为飞行器动压,S为飞行器参考面积,b为飞行器气动系数。
CN201610331421.1A 2016-05-18 2016-05-18 一种抗风干扰的飞行器攻角指令补偿方法及装置 Expired - Fee Related CN106017218B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610331421.1A CN106017218B (zh) 2016-05-18 2016-05-18 一种抗风干扰的飞行器攻角指令补偿方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610331421.1A CN106017218B (zh) 2016-05-18 2016-05-18 一种抗风干扰的飞行器攻角指令补偿方法及装置

Publications (2)

Publication Number Publication Date
CN106017218A true CN106017218A (zh) 2016-10-12
CN106017218B CN106017218B (zh) 2017-03-01

Family

ID=57097535

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610331421.1A Expired - Fee Related CN106017218B (zh) 2016-05-18 2016-05-18 一种抗风干扰的飞行器攻角指令补偿方法及装置

Country Status (1)

Country Link
CN (1) CN106017218B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109976378A (zh) * 2019-03-13 2019-07-05 南京航空航天大学 风扰下无人机栖落机动的轨迹控制方法
CN112278325A (zh) * 2020-09-21 2021-01-29 北京控制工程研究所 一种基于法向过载的总攻角控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2625809A1 (fr) * 1988-01-12 1989-07-13 Aerospatiale Dispositif pour annuler les forces aerodynamiques sur un corps en rotation, et application a l'equilibrage dynamique
CN102749851A (zh) * 2012-07-24 2012-10-24 北京航空航天大学 一种挠性高超声速飞行器的精细抗干扰跟踪控制器
CN103984237A (zh) * 2014-06-04 2014-08-13 西北工业大学 基于运动状态综合识别的轴对称飞行器三通道自适应控制***设计方法
CN105182984A (zh) * 2015-07-16 2015-12-23 孙明玮 飞行器俯仰姿态的线性自抗扰控制器设计与参数整定
CN105278545A (zh) * 2015-11-04 2016-01-27 北京航空航天大学 适用于高超声速机动飞行的自抗扰轨迹线性化控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2625809A1 (fr) * 1988-01-12 1989-07-13 Aerospatiale Dispositif pour annuler les forces aerodynamiques sur un corps en rotation, et application a l'equilibrage dynamique
CN102749851A (zh) * 2012-07-24 2012-10-24 北京航空航天大学 一种挠性高超声速飞行器的精细抗干扰跟踪控制器
CN103984237A (zh) * 2014-06-04 2014-08-13 西北工业大学 基于运动状态综合识别的轴对称飞行器三通道自适应控制***设计方法
CN105182984A (zh) * 2015-07-16 2015-12-23 孙明玮 飞行器俯仰姿态的线性自抗扰控制器设计与参数整定
CN105278545A (zh) * 2015-11-04 2016-01-27 北京航空航天大学 适用于高超声速机动飞行的自抗扰轨迹线性化控制方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109976378A (zh) * 2019-03-13 2019-07-05 南京航空航天大学 风扰下无人机栖落机动的轨迹控制方法
CN112278325A (zh) * 2020-09-21 2021-01-29 北京控制工程研究所 一种基于法向过载的总攻角控制方法
CN112278325B (zh) * 2020-09-21 2022-05-24 北京控制工程研究所 一种基于法向过载的总攻角控制方法

Also Published As

Publication number Publication date
CN106017218B (zh) 2017-03-01

Similar Documents

Publication Publication Date Title
CN102680201B (zh) 基于视频测量的抖振风洞试验方法
CN110017929B (zh) 基于子结构灵敏度分析的船撞桥荷载与损伤同步识别方法
CN112683446B (zh) 一种飞机实时重心位置估计方法
CN105136423B (zh) 考虑摩擦力的自由振动动导数试验的数据分析方法
CN105628051A (zh) 一种嵌入式大气测量装置性能评估方法
CN106017218A (zh) 一种抗风干扰的飞行器攻角指令补偿方法及装置
CN110532621A (zh) 一种飞行器气动参数在线辨识方法
CN102901613B (zh) 一种再入飞行器压力中心确定方法
CN103034125A (zh) 一种飞船返回舱气动与喷流控制力矩参数辨识方法
Cavaliere et al. Method to account for estimator-induced previewed information losses-application to synthesis of lidar-based gust load alleviation functions
Goodman et al. Investigation at low speeds of the effect of aspect ratio and sweep on rolling stability derivatives of untapered wings
CN104787361B (zh) 升力式飞行器再入制导的阻力加速度变化率的确定方法
CN105509689B (zh) 一种用于无人机机载武器发射的三轴校准方法
Nielsen et al. Lift and Pitching-Moment Interference Between a Pointed Cylindrical Body and Triangular Wings of Various Aspect Ratios at Mach Numbers of 1.50 and 2.02
Mayer et al. Flight Measurements with the Douglas D-558-II (buaero No. 37974) Research Airplane: Measurements of the Buffet Boundary and Peak Airplane Normal-force Coefficients at Mach Numbers up to 0.90
Fu et al. Helicopter system identification in the frequency domain
Siu et al. Flight test results of an angle of attack and angle of sideslip calibration method using Output-Error optimization
Colgren et al. Flight Test validation of sideslip estimation using inertial accelerations
Haering Jr Airdata calibration techniques for measuring atmospheric wind profiles
Jandali et al. A potential flow theory for airfoil spoilers
Mineck Wind-tunnel investigation of aerodynamic efficiency of three planar elliptical wings with curvature of quarter-chord line
Carman et al. Lift and Drag Coefficients for the Bell X-1 Airplane (8-Percent-Thick Wing) in Power-Off Transonic Flight
Shukla et al. Modeling & Identification of High Performance Aircraft in Frequency Domain
WIESEMAN Methodology for using steady experimental aerodynamic data to improve steady and unsteady aerodynamic analysis(M. S. Thesis)
Landers et al. Hinge moment coefficient prediction for nose-mounted canard controls at supersonic speeds

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

Granted publication date: 20170301

Termination date: 20200518