CN108416093B - 一种通用的星箭分离时刻卫星飞行状态确定系统及方法 - Google Patents
一种通用的星箭分离时刻卫星飞行状态确定系统及方法 Download PDFInfo
- Publication number
- CN108416093B CN108416093B CN201810089905.9A CN201810089905A CN108416093B CN 108416093 B CN108416093 B CN 108416093B CN 201810089905 A CN201810089905 A CN 201810089905A CN 108416093 B CN108416093 B CN 108416093B
- Authority
- CN
- China
- Prior art keywords
- satellite
- rocket
- variation
- separation
- trajectory
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种通用的星箭分离时刻卫星飞行状态确定系统及方法,主要解决在卫星发射试验任务中,实时数据处理系统在根据火箭的实时测量数据进行指挥决策信息处理时,计算确定星箭分离时刻卫星飞行状态,以进一步进行卫星初轨预报,提高卫星初轨预报的可信度和精度,为各级指挥员提供正确的指挥决策支持信息。包括:1.提出了卫星安装方式对星箭分离时刻卫星飞行状态的影响模型。2.提出了星箭分离方式对星箭分离时刻卫星飞行状态的影响模型。3.确定了星箭分离时刻卫星相对于火箭的发射系弹道变化量的计算方法。4.确定了星箭分离时刻卫星飞行状态的计算方法。
Description
技术领域
本发明属于航天测控领域火箭和航天器的数据处理分析技术,特别是卫星与火箭分离时刻的数据处理分析技术,具体是一种通用的星箭分离时刻卫星飞行状态确定方法。
背景技术
在卫星发射任务中,火箭将卫星推送到预定轨道时,执行星箭分离预定指令。因此卫星的初始轨道参数(以下简称为:卫星初轨根数)是衡量卫星发射任务是否成功的关键指标。
卫星初轨根数由星箭分离时刻卫星飞行状态计算得到。初轨计算算法已经成熟。因此,准确确定卫星初轨的关键是确定星箭分离时刻卫星的运动状态。
受任务特点、政治、经济和时间等因素的影响,发射任务可能是一箭一星、一箭双星、一箭三星或者其它形式。星箭分离方式可能为分离弹簧箭体轴正向弹射、分离弹簧箭体轴斜向弹射、火箭自旋甩分或多种方式的组合分离。此前,针对不同形式的星箭分离方式,实时数据处理系统有单独的卫星分离状态处理模型和函数,导致软件越来越庞大,维护细节较多,容易人为引入缺陷。在高密度发射任务期间,时间紧任务重,软件维护时间短,避免软件维护引入软件缺陷显得尤其困难。
开展通用的星箭分离时刻卫星飞行状态确定方法研究具有较大的现实意义。主要表现在:
1)能够为实时数据处理系统提供统一的星箭分离时刻卫星飞行状态确定方法,适应不同形式的发射任务,从而精简实时数据处理系统,提高软件可靠性。
2)统一的数据输入接口能实现实时数据处理系统的可配置化,通过配置卫星个数、卫星安装参数、星箭分离参数,就能确定星箭分离时刻卫星飞行状态,进而进行卫星初轨预报,降低了软件维护难度,缩短软件维护周期,满足高密度发射任务需求。
发明内容
本发明针对现有实时数据处理系统现状,为了适应发射任务中的不同星箭分离方式,提高实时数据处理系统可靠性,实现软件可配置化,降低软件维护难度,缩短维护周期,提出了一种通用的星箭分离时刻卫星飞行状态确定系统及方法。本发明的技术方案如下:
一种通用的星箭分离时刻卫星飞行状态确定系统,其特征在于,包括位置变化量计算模块、速度变化量计算模块、弹道变化量计算模块及弹道确定模块,其中所述位置变化量计算模块,用于根据根据卫星安装方式,计算星箭分离时刻卫星相对于火箭的箭体系位置变化量;速度变化量计算模块,用于根据星箭分离方式,计算星箭分离时刻卫星相对于火箭的箭体系速度变化量;弹道变化量计算模块,用于将箭体系位置变化量和速度变化量通过箭体系到发射系的坐标转换,计算得到星箭分离时刻卫星相对于火箭的发射系弹道变化量;弹道确定模块,用于将发射系弹道变化量通过弹道合成,计算得到星箭分离时刻卫星在发射系下的弹道,完成星箭分离时刻卫星飞行状态的确定。(位置变化量计算模块和速度变化量计算模块的输出结果是弹道变化量计算模块的输入,弹道变化量计算模块的输出是弹道确定模块的输入。括号内的关系内容在这段文字中已表达清楚)
进一步的,所述位置变化量计算模块中根据卫星安装方式计算星箭分离时刻卫星相对于火箭的箭体系位置变化量具体为:
卫星安装方式对星箭分离时刻卫星飞行状态的影响模型,即卫星相对于火箭的箭体系位置变化量BDX计算模型如下:
其中,Install_Ay为在箭体横截面即YOZ面上箭体Y轴正向与卫星安装点与箭体中心的连线的夹角,逆时针为正;Install_Ax为在箭体纵截面即XOY面上箭体X轴正向与卫星轴向的夹角,逆时针为正;Install_x为卫星质点中心在箭体轴上的坐标;Install_r为卫星质点中心与箭体轴的距离;Mi[a](i=1,2,3)为两个坐标系的方向余弦矩阵,表示第i轴平行,另外两个坐标轴的夹角相同,均为α。
进一步的,
进一步的,所述速度变化量计算模块根据星箭分离方式计算星箭分离时刻卫星相对于火箭的箭体系速度变化量具体为:
星箭分离方式对星箭分离时刻卫星飞行状态的影响模型,即卫星在星箭分离时刻获得的箭体系速度变化量BDV计算模型如下:
其中,FenLi_DVad为卫星在分离时刻从分离弹簧的分离动作获得的卫星轴向速度增量,沿轴向正向方向为正;FenLi_dA为分离时刻箭体轴向与卫星轴向的夹角,逆时针方向为正;FenLi_DVrd为卫星在分离时刻从火箭的自旋获得的径向速度增量,从火箭底部看逆时针方向为正。
进一步的,所述弹道变化量计算模块计算星箭分离时刻卫星相对于火箭的发射系弹道变化量具体为:
进一步的,所述弹道确定模块得到星箭分离时刻卫星在发射系下的弹道,完成星箭分离时刻卫星飞行状态的确定,具体内容为:
弹道结构定义为
tagOrbit=[T X V]
其中,T表示弹道时间,X表示弹道位置矢量,V表示弹道速度矢量;
星箭分离时刻,卫星在发射系下的弹道Orbits如下计算得到
其中,txf为星箭分离时间;Orbittxf、GJXtxf、GJVtxf分别为火箭在txf时刻的发射系弹道、位置和速度;GXXtxf、GXVtxf分别为卫星在txf时刻的发射系弹道位置和速度。
一种基于所述系统的通用的星箭分离时刻卫星飞行状态确定方法,其包括以下步骤:
101、首先,根据卫星安装方式,计算星箭分离时刻卫星相对于火箭的箭体系位置变化量;
102、其次,根据星箭分离方式,计算星箭分离时刻卫星相对于火箭的箭体系速度变化量;
103、将箭体系位置变化量和速度变化量,通过箭体系到发射系的坐标转换,得到星箭分离时刻卫星相对于火箭的发射系弹道变化量;
104、最后,将火箭的发射系弹道变化量通过弹道合成,计算得到星箭分离时刻卫星在发射系下的弹道,完成星箭分离时刻卫星飞行状态的确定。
本发明的优点及有益效果如下:
本发明法对不同卫星安装方式和星箭分离方式进行了统一建模,能够适应不同星箭分离方式的卫星发射任务,从而精简实时数据处理系统结构,提高实时数据处理系统的可维护性和可靠性;为实现软件可配置化提供支持,避免引入人为软件缺陷,降低维护难度,缩短维护周期,满足高密度发射任务需求。
本方法已应用于本中心的实时任务数据处理软件中,成功执行了五次卫星发射任务。
在任务实施过程中,基于本方法成功完成了北斗IGSO卫星(一箭一星发射方式)、北斗MEO卫星(一箭双星发射方式)、遥感卫星(一箭三星发射方式)的初轨预报,与卫星测控中心综合卫星测量数据得到的卫星初轨根数基本一致,方法的正确性和有效性得到验证,为各级指挥员提供了正确的指挥决策支持信息。
附图说明
图1是本发明提供优选实施例实时数据处理系统应用本方法的处理流程示意图;
表1是某次卫星发射中实时数据处理系统应用本方法基于理论数据计算的卫星初轨根数精度;
表2是某次卫星发射中实时数据处理系统应用本方法基于实测数据计算的卫星初轨根数精度。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、详细地描述。所描述的实施例仅仅是本发明的一部分实施例。
本发明解决上述技术问题的技术方案是:
根据上述研究背景,本发明的目的是提出一种通用的星箭分离时刻卫星飞行状态确定方法。本方法包括:
1、卫星安装方式对星箭分离时刻卫星飞行状态的影响模型
卫星安装方式对星箭分离时刻卫星飞行状态的影响模型,即卫星相对于火箭的箭体系位置变化量BDX计算模型为:
其中,Install_Ay为在箭体横截面(YOZ面)上箭体Y轴正向与卫星安装点与箭体中心的连线的夹角,逆时针为正;Install_Ax为在箭体纵截面(XOY面)上箭体X轴正向与卫星轴向的夹角,逆时针为正;Install_x为卫星质点中心在箭体轴上的坐标;Install_r为卫星质点中心与箭体轴的径向距离;Mi[a](i=1,2,3)为两个坐标系的方向余弦矩阵,表示第i轴平行,另外两个坐标轴的夹角相同,均为α,即
Install_Ay、Install_Ax、Install_x、Install_r的具体数值在卫星发射方案设计阶段已经确定。
2、星箭分离方式对星箭分离时刻卫星飞行状态的影响模型
星箭分离方式对星箭分离时刻卫星飞行状态的影响模型,即卫星在星箭分离时刻获得的箭体系速度变化量BDV计算模型如下:
其中,FenLi_DVad为卫星在分离时刻从分离弹簧的分离动作获得的卫星轴向速度增量,沿轴向正向方向为正;FenLi_dA为分离时刻箭体轴向与卫星轴向的夹角,逆时针方向为正;FenLi_DVrd为卫星在分离时刻从火箭的自旋获得的径向速度增量,从火箭底部看逆时针方向为正。FenLi_DVad、FenLi_dA具体数值在卫星发射方案设计阶段已经确定;FenLi_DVrd根据实时遥测参数计算得到。
3、星箭分离时刻卫星相对于火箭的发射系弹道变化量的计算方法
通过箭体系到发射系的坐标转换模型,计算星箭分离时刻卫星相对于火箭的发射系弹道变化量的计算方法如下:
其中,GDX、GDV分别为卫星在分离时刻的发射系弹道位置变化量和弹道速度变化量;M_B2G为分离时刻箭体系到发射系的转换矩阵;γ、ψ、分别为分离时刻箭体相对于发射系的滚动角、偏航角、俯仰角。4、星箭分离时刻卫星飞行状态的计算方法
弹道结构定义为
tagOrbit=[T X V]
其中,T表示弹道时间,X表示弹道位置矢量,V表示弹道速度矢量。
星箭分离时刻,卫星在发射系下的弹道Orbits如下计算得到。
其中,txf为星箭分离时间;Orbittxf、GJXtxf、GJVtxf分别为火箭在txf时刻的发射系弹道、位置和速度;GXXtxf、GXVtxf分别为卫星在txf时刻的发射系弹道位置和速度。
5、方法步骤
下面以伪代码形式,描述实时数据处理系统利用本方法计算分离时刻卫星飞行状态及卫星初轨的具体步骤。
实时数据处理系统中利用本方法的流程如图1所示。
如图1所示,实时数据处理系统利用本方法在星箭分离后计算每颗卫星初轨的流程。首先利用本方法确定星箭分离时刻每颗卫星飞行状态,然后计算该卫星初轨根数。为了将卫星分离时刻卫星飞行状态确定函数做得通用,每颗卫星的安装参数和分离参数单独配置,以适应搭载不同个数卫星和不同分离方式的发射任务及卫星初轨预报任务。
如表1所示,实时数据处理系统应用本方法在某次卫星发射任务准备阶段基于理论数据计算的卫星初轨根数精度。本方法根据理论数据计算的卫星初轨根数与精确弹道设计的卫星初轨根数基本一致。以精确弹道设计的卫星初轨根数为真值,本方法计算的卫星初轨根数的预报精度在计算误差范围之内。
表1
卫星初轨根数 | 绝对误差 | 相对误差 |
半长轴 | 11.353m | 1.627×10-6 |
偏心率 | 0.000003 | 0.00328 |
倾角 | 0.000618<sup>°</sup> | 1.765×10-5 |
升交点经度 | 0.000966<sup>°</sup> | 2.907×10-6 |
如表2所示,实时数据处理系统应用本方法在某次卫星发射任务实施阶段基于实测数据计算的卫星初轨根数精度。在卫星发射过程中,利用实测数据应用本方法实时计算预报的卫星初轨根数与卫星测控中心计算的初轨根数基本一致。卫星测控中心计算的卫星初轨综合了星箭分离时刻火箭测量数据和分离后的卫星测量数据。以卫星测控中心计算结果为卫星初轨真值,本方法计算的卫星初轨根数的预报精度在计算误差范围之内。
表2
卫星初轨根数 | 绝对误差 | 相对误差 |
半长轴 | 554.881m | 7.95×10<sup>-5</sup> |
偏心率 | 0.0000734 | 0.05258 |
倾角 | 0.00049<sup>°</sup> | 1.399×10<sup>-5</sup> |
升交点经度 | 0.00048<sup>°</sup> | 1.445×10<sup>-6</sup> |
以上这些实施例应理解为仅用于说明本发明而不用于限制本发明的保护范围。在阅读了本发明的记载的内容之后,技术人员可以对本发明作各种改动或修改,这些等效变化和修饰同样落入本发明权利要求所限定的范围。
Claims (7)
1.一种通用的星箭分离时刻卫星飞行状态确定系统,其特征在于,包括位置变化量计算模块、速度变化量计算模块、弹道变化量计算模块及弹道确定模块,其中所述位置变化量计算模块,用于根据卫星安装方式,计算星箭分离时刻卫星相对于火箭的箭体系位置变化量;速度变化量计算模块,用于根据星箭分离方式,计算星箭分离时刻卫星相对于火箭的箭体系速度变化量;弹道变化量计算模块,用于将箭体系位置变化量和速度变化量通过箭体系到发射系的坐标转换,计算得到星箭分离时刻卫星相对于火箭的发射系弹道变化量;弹道确定模块,用于将发射系弹道变化量通过弹道合成,计算得到星箭分离时刻卫星在发射系下的弹道,完成星箭分离时刻卫星飞行状态的确定,位置变化量计算模块和速度变化量计算模块的输出结果是弹道变化量计算模块的输入,弹道变化量计算模块的输出是弹道确定模块的输入。
2.根据权利要求1所述的通用的星箭分离时刻卫星飞行状态确定系统,其特征在于,所述位置变化量计算模块中根据卫星安装方式计算星箭分离时刻卫星相对于火箭的箭体系位置变化量具体为:
卫星安装方式对星箭分离时刻卫星飞行状态的影响模型,即卫星相对于火箭的箭体系位置变化量BDX计算模型如下:
其中,Install_Ay为在箭体横截面即YOZ面上箭体Y轴正向与卫星安装点与箭体中心的连线的夹角,逆时针为正;Install_Ax为在箭体纵截面即XOY面上箭体X轴正向与卫星轴向的夹角,逆时针为正;Install_x为卫星质点中心在箭体轴上的坐标;Install_r为卫星质点中心与箭体轴的距离;Mi[α](i=1,2,3)为两个坐标系的方向余弦矩阵,表示第i轴平行,另外两个坐标轴的夹角相同,均为α。
7.一种基于权利要求1-6之一所述系统的通用的星箭分离时刻卫星飞行状态确定方法,其特征在于,包括以下步骤:
101、首先,根据卫星安装方式,计算星箭分离时刻卫星相对于火箭的箭体系位置变化量;
102、其次,根据星箭分离方式,计算星箭分离时刻卫星相对于火箭的箭体系速度变化量;
103、将箭体系位置变化量和速度变化量,通过箭体系到发射系的坐标转换,得到星箭分离时刻卫星相对于火箭的发射系弹道变化量;
104、最后,将火箭的发射系弹道变化量通过弹道合成,计算得到星箭分离时刻卫星在发射系下的弹道,完成星箭分离时刻卫星飞行状态的确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810089905.9A CN108416093B (zh) | 2018-01-30 | 2018-01-30 | 一种通用的星箭分离时刻卫星飞行状态确定系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810089905.9A CN108416093B (zh) | 2018-01-30 | 2018-01-30 | 一种通用的星箭分离时刻卫星飞行状态确定系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108416093A CN108416093A (zh) | 2018-08-17 |
CN108416093B true CN108416093B (zh) | 2021-10-12 |
Family
ID=63127461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810089905.9A Expired - Fee Related CN108416093B (zh) | 2018-01-30 | 2018-01-30 | 一种通用的星箭分离时刻卫星飞行状态确定系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108416093B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115358098B (zh) * | 2022-10-20 | 2023-05-12 | 北京宇航系统工程研究所 | 一种远场安全性分析方法、系统、电子设备及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009186355A (ja) * | 2008-02-07 | 2009-08-20 | Ihi Aerospace Co Ltd | 航行体の姿勢検出方法及び装置 |
CN105588568A (zh) * | 2015-12-11 | 2016-05-18 | 中国人民解放军63791部队 | 火箭发射无动力飞行段弹道外推与滤值方法 |
CN105631229A (zh) * | 2016-01-19 | 2016-06-01 | 中国人民解放军63796部队 | 一种运载器滑行段测量盲区运行轨迹填补方法 |
CN106570316A (zh) * | 2016-10-20 | 2017-04-19 | 北京空间飞行器总体设计部 | 基于推进剂预算的低轨椭圆轨道卫星成功入轨判定方法 |
CN107607947A (zh) * | 2017-08-22 | 2018-01-19 | 西安电子科技大学 | 基于卡尔曼滤波的星载雷达成像参数在线估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9091552B2 (en) * | 2011-10-25 | 2015-07-28 | The Boeing Company | Combined location and attitude determination system and methods |
-
2018
- 2018-01-30 CN CN201810089905.9A patent/CN108416093B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009186355A (ja) * | 2008-02-07 | 2009-08-20 | Ihi Aerospace Co Ltd | 航行体の姿勢検出方法及び装置 |
CN105588568A (zh) * | 2015-12-11 | 2016-05-18 | 中国人民解放军63791部队 | 火箭发射无动力飞行段弹道外推与滤值方法 |
CN105631229A (zh) * | 2016-01-19 | 2016-06-01 | 中国人民解放军63796部队 | 一种运载器滑行段测量盲区运行轨迹填补方法 |
CN106570316A (zh) * | 2016-10-20 | 2017-04-19 | 北京空间飞行器总体设计部 | 基于推进剂预算的低轨椭圆轨道卫星成功入轨判定方法 |
CN107607947A (zh) * | 2017-08-22 | 2018-01-19 | 西安电子科技大学 | 基于卡尔曼滤波的星载雷达成像参数在线估计方法 |
Non-Patent Citations (3)
Title |
---|
GPS-Based Attitude Determination for a Spinning Rocket;Karen Chiang等;《IEEE Transactions on Aerospace and Electronic Systems》;20141031;第50卷(第4期);第2654-2663页 * |
一种可实现星箭分离前后测轨数据联合定轨的初轨确定方法;茅永兴等;《宇航学报》;20141231;第35卷(第12期);第1359-1366页 * |
不同星箭分离方式下整星冲击环境特征分析;杨艳静等;《装备环境工程》;20170831;第14卷(第8期);第70-78页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108416093A (zh) | 2018-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107544467B (zh) | 整星条件下的双星编队控制闭环测试系统及方法 | |
CN107765226B (zh) | 一种sar卫星雷达回波模拟方法、系统和介质 | |
CN111301715B (zh) | 基于霍曼变轨的同轨道特定相位分布的星座布局与轨道调整方法、装置及计算机存储介质 | |
CN108562295B (zh) | 一种基于同步卫星二体模型的三站时差定轨方法 | |
CN109781059B (zh) | 星载点波束天线对地指向精度评估系统 | |
CN104236546A (zh) | 一种卫星星光折射导航误差确定与补偿方法 | |
CN104898642A (zh) | 一种用于航天器姿态控制算法的集成测试仿真系统 | |
CN107300700B (zh) | 敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法 | |
CN110816896B (zh) | 一种卫星星上简易轨道外推方法 | |
CN104751012A (zh) | 沿飞行弹道的扰动引力快速逼近方法 | |
CN103645489A (zh) | 一种航天器gnss单天线定姿方法 | |
CN116382124B (zh) | 一种运载火箭的姿态控制仿真方法和系统 | |
CN114936471A (zh) | 一种基于并行计算的航天器碰撞预警分层快速筛选方法 | |
CN108416093B (zh) | 一种通用的星箭分离时刻卫星飞行状态确定系统及方法 | |
CN105487405A (zh) | 低低跟踪重力测量卫星半物理仿真系统 | |
CN115560757B (zh) | 随机姿态误差条件下基于神经网络的无人机直接定位校正方法 | |
CN114781275B (zh) | 基于人工智能的航天器轨道拦截的燃料控制方法、装置及介质 | |
CN111319801A (zh) | 一种适用于火星探测的中途修正策略制定及实施方法 | |
CN114935277B (zh) | 一种滑翔增程制导炮弹理想弹道的在线规划方法 | |
CN115902980A (zh) | 可在线更新的飞行目标天基监视与预测方法及装置 | |
CN113108779B (zh) | 一种基于1553b总线三冗余的自主箭测系统 | |
CN111272336B (zh) | 一种基于gnss观测实现大型低轨航天器质心位移估算的方法 | |
CN115268390A (zh) | 卫星高精度跟踪指向控制地面仿真系统ahp效能评估方法 | |
CN112257172B (zh) | 一种基于离散点计算的卫星天线指向角正确性的验证方法 | |
CN114154355A (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: 20211012 Termination date: 20220130 |