CN108803646B - 一种变增益章动阻尼器的实现方法 - Google Patents
一种变增益章动阻尼器的实现方法 Download PDFInfo
- Publication number
- CN108803646B CN108803646B CN201810889031.5A CN201810889031A CN108803646B CN 108803646 B CN108803646 B CN 108803646B CN 201810889031 A CN201810889031 A CN 201810889031A CN 108803646 B CN108803646 B CN 108803646B
- Authority
- CN
- China
- Prior art keywords
- nutation
- rotating speed
- damper
- speed
- frequency
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000013139 quantization Methods 0.000 claims abstract description 35
- 230000000694 effects Effects 0.000 claims abstract description 23
- 230000008569 process Effects 0.000 claims abstract description 10
- 238000013461 design Methods 0.000 claims abstract description 5
- 230000005284 excitation Effects 0.000 claims abstract description 5
- 238000002474 experimental method Methods 0.000 claims description 8
- 230000004044 response Effects 0.000 claims description 5
- 238000012546 transfer Methods 0.000 claims description 3
- 238000011426 transformation method Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 2
- 238000001228 spectrum Methods 0.000 claims description 2
- 230000009471 action Effects 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 230000001629 suppression Effects 0.000 description 4
- 238000013016 damping Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Vibration Prevention Devices (AREA)
Abstract
一种变增益章动阻尼器的实现方法,涉及变速倾侧动量轮运动体控制领域。为了解决变速倾侧动量轮由于转子陀螺效应在外加激励作用下产生章动造成系统失稳的问题。技术要点:得到章动频率与动量轮转速的函数关系;根据变速倾侧动量轮运行时的实测转速,确定测速过程中引入的随机噪声变化范围,通过在转速时变和转速稳定条件下章动阻尼器效果的对比,确定章动阻尼器转速量化步长并得到章动阻尼器中心频率;根据章动阻尼器转速量化步长和章动谐振峰幅值大小确定章动阻尼器的参数完成章动阻尼器的设计;将设计的章动阻尼器进行离散化处理以对变速倾侧动量轮进行章动抑制。本发明考虑不同因素的影响,设定转速量化步长,得到效果更优的变增益章动阻尼器。
Description
技术领域
本发明涉及变速倾侧动量轮运动体控制领域,具体涉及一种变增益章动阻尼器的实现方法。
背景技术
变速倾侧动量轮作为一种适用于微小航天器姿态敏感和控制的新型集成化装置,借鉴变速双框架控制力矩陀螺和动力调谐陀螺的思想,通过控制一个扁平转子体角动量幅值和指向变化,能够同时实现三轴姿控力矩输出和两轴姿态角速率测量功能。变速倾侧动量轮的应用将很大程度地提高微小航天器姿控系统的集成度和效率,减小系统的质量、体积、功耗和研制成本。
转子受到陀螺效应的影响,存在进动和章动两种涡动模态,其中转子章动是变速倾侧动量轮失稳的主要形式。由于转子在运行中需要实时改变转速,因此变速倾侧动量轮对应的章动频率也是时变的。为了抑制章动造成的系统失稳,需要设计依赖转速变化的变增益章动阻尼器。在已有专利(专利名称为:基于变参数章动阻尼的变速倾侧动量轮进动控制方法,专利号为:201711498624.0)中提出了一种变增益章动阻尼器的设计方法。但在实际数字控制系统中,一方面控制器是离散的,另一方面由于实际系统转子转速存在小范围波动,使用参数实时随转速变化的阻尼器往往不能达到预期的章动抑制效果。基于此,有必要提出适用于数字控制系统且适应转子转速变化和存在一定程度波动的变增益章动阻尼器的计算机实现方法。
发明内容
针对变速倾侧动量轮由于转子陀螺效应在外加激励作用下产生章动造成系统失稳的问题,本发明提出了一种针对转速测量存在噪声情况下的变增益章动阻尼器的实现方法。
本发明为解决上述技术问题采取的技术方案是:
一种变增益章动阻尼器的实现方法,所述方法具体为:
步骤一、通过动量轮倾侧回路开环扫频实验,确定变速倾侧动量轮在不同转速下的章动频率点,对转速-章动频率进行拟合,得到章动频率与动量轮转速的函数关系;
步骤二、根据变速倾侧动量轮运行时的实测转速,确定测速过程中引入的随机噪声变化范围,进而通过在转速时变和转速稳定条件下章动阻尼器效果的对比,确定章动阻尼器转速量化步长并得到章动阻尼器中心频率,章动阻尼器的转速在所述量化步长范围内不受噪声影响而改变;
步骤三、根据章动阻尼器转速量化步长和章动谐振峰幅值大小确定章动阻尼器的参数:串联陷波器的数量、中心陷波深度和陷波带宽,完成章动阻尼器的设计;将设计的章动阻尼器进行离散化处理以对变速倾侧动量轮进行章动抑制。
步骤一、通过倾侧回路开环实验,对得到的数据进行线性拟合,得到系统章动频率与动量轮转速的关系,具体为:
(a)令变速倾侧动量轮以转速nzi稳定运行,在倾侧控制回路开环的条件下,给倾侧轴x轴方向分别施加幅值为Ai,频率为fi的正弦电压,同时保持与x轴正交的y轴方向施加电压为0,测量正交于激励的y轴方向的开环稳态响应,记录在转速nzi下的幅值比ki,其中,i=1,2,…,用于区别不同转速;下角标z表示z轴;
(b)得到转速nzi下的幅值特性的谐振峰对应的频率为fri,此频率即确定为nzi转速下的动量轮章动频率;
(c)改变设定的运行转速nz,重复上述步骤(a)、(b),得到尽可能多的章动频率-转速数据,对所得数据进行线性拟合,得到章动频率随转速变化的函数关系fn=f(nz)。
在步骤二中,确定章动阻尼器转速量化步长并得到章动阻尼器中心频率的过程为:
(a)令变速倾侧动量轮在旋转轴z轴方向闭环运行,记录测得的变速倾侧动量轮(飞轮)转速,通过对转速数据的分析,得到噪声幅值上限Anoise;
(b)在变速倾侧动量轮转速稳定的情况下,为使章动阻尼器效果最优取转速量化步长nth1≥Anoise;
(c)在变速倾侧动量轮转速存在动态过程的情况下,应取章动阻尼器的转速量化步长nth2小于转速可控制步长;所述转速可控制步长是指变速倾侧动量轮能够实现的最小转速变化量;
(d)综合变速倾侧动量轮转速稳定和转速存在动态过程的情况,得到阻尼器的最终转速量化步长nth;nth的取值介于[nth2,nth1]之间;
在步骤三中,设计变增益章动阻尼器的参数,章动阻尼器由一个或多个串联的二阶陷波器构成,再对连续的章动阻尼器进行离散化处理;
(a)一个二阶陷波器的传递函数如式(1)所示
其中,s为拉普拉斯算子,a/b表示陷波器的中心陷波深度,陷波深度应大于所述章动谐振峰幅值;b表示陷波器的陷波宽度(陷波带宽),陷波宽度应大于章动谐振峰的频谱宽度;fr表示章动阻尼器中心频率;
(b)将设计的章动阻尼器进行离散化,常用的离散化方式为双线性变换法,可利用下式实现
式中,T表示离散控制系统的采样时间,z表示离散变换算子。
本发明具有以下有益效果:
本发明通过倾侧回路开环实验确定章动频率随转速的变化规律,在已知转速测量值存在噪声的情况下,改进了变增益章动阻尼器的实现策略,设计了中心频率在一定转速变化范围内稳定的变增益章动阻尼器,其对章动失稳的抑制效果如表1所示(表1数据均为变速倾侧动量轮在转速3200rpm下运行时得到)。由表1可知,当转速处于稳态且测量存在噪声时,章动阻尼器转速量化步长大于噪声幅值上限的陷波器效果最佳,优于转速量化步长小于噪声幅值上限的陷波器,而中心频率随转速实时变化的陷波器效果最差。另一方面,由图5可知,当转速动态变化时,使用转速量化步长越小的陷波器,其章动抑制效果越好。因此,在实际系统中需要综合考虑上述两方面的影响,设定转速量化步长,得到效果更优的变增益章动阻尼器。
附图说明
图1为本发明的流程图;
图2为变速倾侧动量轮从2700rpm至3600rpm稳速运行时的交轴开环频率特性,其中转速间隔为300rpm;
图3为变速倾侧动量轮转速-章动频率实际数据和拟合曲线图,其中实线为拟合后的章动-转速关系曲线,离散点为图2中章动频率-转速原始数据;
图4为转速为3200±3rpm情况下使用不同转速量化步长的闭环频谱及章动频率附近局部放大图;其中,采用陷波中心频率随测量转速实时变化的章动阻尼器的章动残余幅值最大,阻尼器效果最差;采用转速量化步长小于测量转速噪声幅值的章动阻尼器的章动残余幅值居中,阻尼器效果一般;采用转速量化步长大于转速测量噪声的章动阻尼器的章动残余幅值最小,阻尼效果最好;
图5为转速从3100rpm均匀升速到3300rpm过程中使用不同转速量化步长章动阻尼器的x轴倾侧角时域响应曲线;在转速变化情况下,选取转速量化步长越大,则相应的陷波效果越差;
图6为不使用章动阻尼器系统发散和使用章动阻尼器系统稳定的倾侧角响应曲线,其中,实线表示使用变增益章动阻尼器时倾侧角稳定在0°,虚线表示不使用变增益章动阻尼器时倾侧角发散。
具体实施方式
步骤一、通过倾侧回路开环实验,对得到的数据进行线性拟合,得到系统章动频率与动量轮转速的关系。
(a)令变速倾侧动量轮以转速nzi稳定运行,在倾侧控制回路开环的条件下,给倾侧轴x轴方向分别施加幅值为Ai,频率为fi的正弦电压,同时保持与其正交的y轴方向施加电压为0,测量正交于激励的y轴方向的开环稳态响应,记录在转速nzi下的幅值比ki,其中,i=1,2,…;
(b)得到转速nzi下的幅值特性的谐振峰对应的频率为fri,此频率即确定为nzi转速下的动量轮章动频率;
(c)改变设定的运行转速nz,重复上述步骤(a)、(b),得到尽可能多的章动频率-转速数据,对所得数据进行线性拟合,得到章动频率随转速变化的函数关系fr=f(nz)。
步骤二、通过转速闭环实验,记录测量转速的稳定情况,判断应设定的转速量化步长并得到章动阻尼器中心频率计算式。
(a)令变速倾侧动量轮在z轴(旋转轴)方向闭环运行,记录测得的飞轮转速信息,通过对数据的分析,得到噪声幅值上限Anoise;
(b)根据表1、图4得到的结论,为使章动阻尼器效果最优,应取转速量化步长nth≥Anoise;
(c)根据图5,应进行转速存在动态过程的实验,综合步骤(b)中的信息,得到阻尼器的转速量化步长nth;
步骤三、设计变增益章动阻尼器的参数,最终使用的章动阻尼器可由一个或多个串联的二阶陷波器构成,再对连续的章动阻尼器进行离散化处理。
(a)一个二阶陷波器的传递函数如式(1)所示
其中,s为拉普拉斯算子,a/b表示陷波器的陷波深度,b表示陷波器的陷波宽度。
(b)将设计的章动阻尼器进行离散化,常用的离散化方式为双线性变换法,可利用下式实现
对本发明的技术效果进行如下验证:
表1为转速量化步长在不同测量噪声条件下(转速分别取3200±3rmp和3200±6rmp)的章动抑制效果对比。表中数据为在不同转速噪声幅值下,使用不同转速量化步长的陷波器的剩余章动幅值,剩余章动幅值越小,效果越好。由表1可以看出,采用陷波中心频率随测量转速实时变化的章动阻尼器的章动残余幅值最大,阻尼器效果最差;采用转速量化步长小于测量转速噪声幅值的章动阻尼器的章动残余幅值居中,阻尼器效果一般;采用转速量化步长大于转速测量噪声的章动阻尼器的章动残余幅值最小,阻尼效果最好。
表1转速量化步长在不同测量噪声条件下的章动抑制效果对比
Claims (1)
1.一种变增益章动阻尼器的实现方法,其特征在于,所述方法具体为:
步骤一、通过动量轮倾侧回路开环扫频实验,确定变速倾侧动量轮在不同转速下的章动频率点,对转速-章动频率进行拟合,得到章动频率与动量轮转速的函数关系;
在步骤一中,通过倾侧回路开环实验,对得到的数据进行线性拟合,得到系统章动频率与动量轮转速的关系,具体为:
(a)令变速倾侧动量轮以转速nzi稳定运行,在倾侧控制回路开环的条件下,给倾侧轴x轴方向分别施加幅值为Ai,频率为fi的正弦电压,同时保持与x轴正交的y轴方向施加电压为0,测量正交于激励的y轴方向的开环稳态响应,记录在转速nzi下的幅值比ki,其中,i=1,2,...,用于区别不同转速;下角标z表示z轴;
(b)得到转速nzi下的幅值特性的谐振峰对应的频率为fri,此频率即确定为nzi转速下的动量轮章动频率;
(c)改变设定的运行转速nz,重复上述步骤(a)、(b),得到尽可能多的章动频率-转速数据,对所得数据进行线性拟合,得到章动频率随转速变化的函数关系fn=f(nz);
步骤二、根据变速倾侧动量轮运行时的实测转速,确定测速过程中引入的随机噪声变化范围,进而通过在转速时变和转速稳定条件下章动阻尼器效果的对比,确定章动阻尼器转速量化步长并得到章动阻尼器中心频率,章动阻尼器的转速在所述量化步长范围内不受噪声影响而改变;
在步骤二中,确定章动阻尼器转速量化步长并得到章动阻尼器中心频率的过程为:
(a)令变速倾侧动量轮在旋转轴z轴方向闭环运行,记录测得的变速倾侧动量轮转速,通过对转速数据的分析,得到噪声幅值上限Anoise;
(b)在变速倾侧动量轮转速稳定的情况下,为使章动阻尼器效果最优取转速量化步长nth1≥Anoise;
(c)在变速倾侧动量轮转速存在动态过程的情况下,应取章动阻尼器的转速量化步长nth2小于转速可控制步长;所述转速可控制步长是指变速倾侧动量轮能够实现的最小转速变化量;
(d)综合变速倾侧动量轮转速稳定和转速存在动态过程的情况,得到阻尼器的最终转速量化步长nth;nth的取值介于[nth2,nth1]之间;
步骤三、根据章动阻尼器转速量化步长和章动谐振峰幅值大小确定章动阻尼器的参数:串联陷波器的数量、中心陷波深度和陷波带宽,完成章动阻尼器的设计;将设计的章动阻尼器进行离散化处理以对变速倾侧动量轮进行章动抑制;
在步骤三中,设计变增益章动阻尼器的参数,章动阻尼器由一个或多个串联的二阶陷波器构成,再对连续的章动阻尼器进行离散化处理;
(a)一个二阶陷波器的传递函数如式(1)所示
其中,s为拉普拉斯算子,a/b表示陷波器的中心陷波深度,陷波深度应大于所述章动谐振峰幅值;b表示陷波器的陷波宽度,陷波宽度应大于章动谐振峰的频谱宽度;fr表示章动阻尼器中心频率;
(b)将设计的章动阻尼器进行离散化,常用的离散化方式为双线性变换法,可利用下式实现
式中,T表示离散控制系统的采样时间,z表示离散变换算子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810889031.5A CN108803646B (zh) | 2018-08-03 | 2018-08-03 | 一种变增益章动阻尼器的实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810889031.5A CN108803646B (zh) | 2018-08-03 | 2018-08-03 | 一种变增益章动阻尼器的实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108803646A CN108803646A (zh) | 2018-11-13 |
CN108803646B true CN108803646B (zh) | 2021-03-23 |
Family
ID=64079497
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810889031.5A Active CN108803646B (zh) | 2018-08-03 | 2018-08-03 | 一种变增益章动阻尼器的实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108803646B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6439509B1 (en) * | 1995-07-03 | 2002-08-27 | Space Systems/Loral, Inc. | Pre-bias scheme for roll momentum unloads on pitch momentum bias/yaw reaction wheel spacecraft |
CN101247097A (zh) * | 2008-03-03 | 2008-08-20 | 北京航空航天大学 | 一种设计磁悬浮扁平高速转子系统陷波器参数的方法 |
CN101554926A (zh) * | 2009-05-20 | 2009-10-14 | 上海微小卫星工程中心 | 航天器的姿态控制系统及方法 |
CN102425554A (zh) * | 2011-11-10 | 2012-04-25 | 北京中科科仪技术发展有限责任公司 | 一种磁悬浮分子泵增益调度控制方法 |
US8321076B2 (en) * | 2009-12-18 | 2012-11-27 | The Boeing Company | On-line inertia estimation for use in controlling an aerospace vehicle |
CN103072701A (zh) * | 2013-01-30 | 2013-05-01 | 北京控制工程研究所 | 一种欠驱动的卫星消旋控制方法 |
CN103072703A (zh) * | 2013-01-30 | 2013-05-01 | 北京控制工程研究所 | 一种欠驱动的卫星进动控制方法 |
US8620496B2 (en) * | 2008-07-23 | 2013-12-31 | The Boeing Company | Systems and method of controlling a spacecraft using attitude sensors |
CN106873611A (zh) * | 2017-01-22 | 2017-06-20 | 北京航空航天大学 | 一种多通道线性自抗扰控制器的设计方法 |
CN107134999A (zh) * | 2017-05-09 | 2017-09-05 | 徐州润泽电气有限公司 | 一种自适应变增益快速响应锁相环改进方法 |
CN107992063A (zh) * | 2017-12-29 | 2018-05-04 | 哈尔滨工业大学 | 基于变参数章动阻尼的变速倾侧动量轮进动控制方法 |
-
2018
- 2018-08-03 CN CN201810889031.5A patent/CN108803646B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6439509B1 (en) * | 1995-07-03 | 2002-08-27 | Space Systems/Loral, Inc. | Pre-bias scheme for roll momentum unloads on pitch momentum bias/yaw reaction wheel spacecraft |
CN101247097A (zh) * | 2008-03-03 | 2008-08-20 | 北京航空航天大学 | 一种设计磁悬浮扁平高速转子系统陷波器参数的方法 |
US8620496B2 (en) * | 2008-07-23 | 2013-12-31 | The Boeing Company | Systems and method of controlling a spacecraft using attitude sensors |
CN101554926A (zh) * | 2009-05-20 | 2009-10-14 | 上海微小卫星工程中心 | 航天器的姿态控制系统及方法 |
US8321076B2 (en) * | 2009-12-18 | 2012-11-27 | The Boeing Company | On-line inertia estimation for use in controlling an aerospace vehicle |
CN102425554A (zh) * | 2011-11-10 | 2012-04-25 | 北京中科科仪技术发展有限责任公司 | 一种磁悬浮分子泵增益调度控制方法 |
CN103072701A (zh) * | 2013-01-30 | 2013-05-01 | 北京控制工程研究所 | 一种欠驱动的卫星消旋控制方法 |
CN103072703A (zh) * | 2013-01-30 | 2013-05-01 | 北京控制工程研究所 | 一种欠驱动的卫星进动控制方法 |
CN106873611A (zh) * | 2017-01-22 | 2017-06-20 | 北京航空航天大学 | 一种多通道线性自抗扰控制器的设计方法 |
CN107134999A (zh) * | 2017-05-09 | 2017-09-05 | 徐州润泽电气有限公司 | 一种自适应变增益快速响应锁相环改进方法 |
CN107992063A (zh) * | 2017-12-29 | 2018-05-04 | 哈尔滨工业大学 | 基于变参数章动阻尼的变速倾侧动量轮进动控制方法 |
Non-Patent Citations (1)
Title |
---|
带章动阻尼器航天器的非线性动力学分析;孔令勇;《商丘师范学院学报》;20070630;第49-53页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108803646A (zh) | 2018-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2807332B1 (en) | Method for reducing drillstring oscillations | |
CN101262199B (zh) | 电动机控制装置以及电动机控制系统 | |
CN109995300B (zh) | 用于伺服系统谐振抑制、陷波参数优化方法、系统及介质 | |
CN107992063B (zh) | 基于变参数章动阻尼的变速倾侧动量轮进动控制方法 | |
CN107631738B (zh) | 一种陀螺飞轮偶不平衡辨识方法 | |
CN111251900B (zh) | 一种能降低扭矩损失的抑抖方法 | |
WO2007096993A1 (ja) | モータ制御装置 | |
CN109347391B (zh) | 一种考虑系统噪声的朗道自适应转动惯量辨识方法 | |
CN108351621B (zh) | 用于工业涡轮机的自适应多输入多输出pid控制系统 | |
EP3791061A1 (en) | Rotor control system for reducing structural vibrations based on m-blade transformation | |
CN111752153B (zh) | 一种基于1.5阶混合重复控制器的谐波电流抑制方法 | |
CN108803646B (zh) | 一种变增益章动阻尼器的实现方法 | |
CN108170027B (zh) | 云台控制器控制参数调整方法、装置、介质和云台控制器 | |
US20160243916A1 (en) | Damper control device | |
JP6481792B2 (ja) | 動力系の試験装置 | |
CN110058525B (zh) | 分数阶离心飞轮调速系统有限时间镇定控制方法 | |
Li et al. | Complex signal processing for Coriolis mass flow metering in two-phase flow | |
CN110535397A (zh) | 一种基于参数已知的电机控制方法及系统 | |
KR101251285B1 (ko) | 풍력 발전기의 토크 제어 방법 | |
CN113879566B (zh) | 一种控制力矩陀螺框架系统高频扰动的自适应补偿方法 | |
CN108809205A (zh) | 通过有源的标称转速调制来低噪声地操控电动机 | |
CN110299868B (zh) | 一种超声波电机转速的控制方法及装置 | |
Estupiñán et al. | Performance evaluation of inverted pendulum control: Linear and nonlinear techniques | |
KR100794893B1 (ko) | 모터 제어 장치 | |
JP2001005502A (ja) | 重み同定方法及びフィードバック制御方法 |
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 |