CN106441368A - 一种光纤陀螺光纤环温变特性测量方法及装置 - Google Patents

一种光纤陀螺光纤环温变特性测量方法及装置 Download PDF

Info

Publication number
CN106441368A
CN106441368A CN201610933044.9A CN201610933044A CN106441368A CN 106441368 A CN106441368 A CN 106441368A CN 201610933044 A CN201610933044 A CN 201610933044A CN 106441368 A CN106441368 A CN 106441368A
Authority
CN
China
Prior art keywords
fiber
temperature
fiber optic
optic loop
optical fibre
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
CN201610933044.9A
Other languages
English (en)
Other versions
CN106441368B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201610933044.9A priority Critical patent/CN106441368B/zh
Publication of CN106441368A publication Critical patent/CN106441368A/zh
Application granted granted Critical
Publication of CN106441368B publication Critical patent/CN106441368B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Gyroscopes (AREA)

Abstract

本发明公开了一种光纤陀螺光纤环温变特性测量方法及装置。将光纤陀螺的光纤环置于温控环境的转台上,将光纤陀螺非光纤环部分置于外部转台,在温控环境温度变化控制的同时控制转台进行特定角速度控制,采集光纤陀螺的耦合器输出的数据利用数据处理算法进行处理,获得光纤环温变特性;所述测量装置连接在对光纤陀螺的光学系统上,包括探测器、信号放大器、模数转换器、数字信号处理芯片和置于温控箱下的转台。本发明方法可对光纤陀螺光纤环结构温度性能进行测试评估,方法简单可靠,测量时间短,为光纤陀螺标度因数高温度稳定性的提高提供基础。

Description

一种光纤陀螺光纤环温变特性测量方法及装置
技术领域
本发明涉及光纤陀螺光纤环温变特性测量,主要用于光纤环温度性能的提升方案的筛选,与光纤陀螺光学标度因数相关的温度效应测试,以及光纤陀螺标度因数温度稳定性提升,方法简单可靠,测量时间短。
背景技术
光纤陀螺是一种敏感角速率的光纤传感器,具有无运动部件、工艺简单、精度覆盖面广、动态范围大、启动快、寿命长、抗冲击、耐过载等优点。光纤陀螺在航空、航天、航海和兵器等军用领域具有广阔的发展前景,引起世界各国的关注。进入21世纪后,光纤陀螺已成为惯性测量和制导领域的主流仪表之一。
光纤环是光纤陀螺中用于敏感输入角速度引起的两相向传播光波之间Sagnac相位差的部件,因而是影响光纤陀螺性能的主要因素。光纤环可以分为脱骨架的光纤环和有骨架的光纤环。
由Sagnac效应可知,光纤环的光纤长度越长、平均直径越大,则标度因数越大,光纤陀螺灵敏度越高,动态范围越小。故通过改变光纤环的结构尺寸可以改变光纤陀螺的灵敏度和动态范围。
光纤环的光纤长度L、平均直径D影响光纤陀螺的测量范围、灵敏度和陀螺的外形尺寸。L和D的乘积越大,光纤陀螺越敏感。
光纤环本身的热胀冷缩和骨架的热膨胀会引起L和D的乘积的变化,直接影响光纤陀螺的标度因数,因此温度的变化将引起光纤陀螺标度因数的变化。
光纤环的直径和光纤的长度对光纤陀螺的性能影响在Sagnac效应表达式中有明确的体现。因而光纤陀螺在温变环境下标度因数会发生变化。国外可以将该变化控制在1ppm/℃,相比之下,国内在这方面仍存在差距。为改善标度因数温度效应,其中最重要方案之一的是改善光纤环温变特性,而因此本发明提出了一种对光纤陀螺光纤环温变特性进行测试的专用方法和装置。
发明内容
针对光纤陀螺光纤环温变特性,本发明的目的在于提供了一种光纤陀螺光纤环温变特性测量方法及装置。
本发明的技术方案是:
一、一种光纤陀螺光纤环温变特性测量方法:
将光纤陀螺的光纤环置于温控环境的转台上,将光纤陀螺非光纤环部分置于外部,在温控环境温度变化控制的同时控制转台进行特定角速度控制,采集光纤陀螺的耦合器输出的数据利用数据处理算法进行处理,获得光纤环温变特性。
所述采集光纤陀螺的耦合器输出的数据的时刻是温控环境内温度稳定的时刻。温控环境的温度变化一般可以是间隔20摄氏度依次递增变化,再递减,以此循环。
所述的特定角速度控制是将转台从静止(角速度为0)开始,以恒定加速度a加速直至角速度为ω0,随后再以恒定加速度a减速直到角速度为-ω0,接着以加速度a加速直到角速度为0,以此重复循环50次以上;
上述角速度ω0满足:ω0>λc/LD,式中λ为光源平均波长,波动忽略不计,c为真空中光速,D为光纤环的线圈平均直径,L为光纤环的光纤长度。
所述的数据处理算法具体是:针对采集获得的耦合器输出数据,通过时域平均获得角速度从-ω0到ω0的时域平均数据,用时域平均数据进行拟合获得数据时域周期Tt,再采用以下公式计算获得光纤长度L和平均直径D乘积随温度变化的光纤环特征参数:
式中,a为转台角加速度,Tt为探测光强的时域周期,λ为光源平均波长,波动忽略不计,c为真空中光速,D为光纤环的线圈平均直径,L为光纤环的光纤长度,T为测试温度,[LD]T为测试温度T下L与D的乘积,即光纤环特征参数。
二、一种光纤陀螺光纤环温变特性测量装置:
所述测量装置连接在对光纤陀螺的光学系统(包括光源、控温电路、耦合器、Y波导相位调制器、光纤环)上,对光学系统进行测试,包括探测器、信号放大器、模数转换器、数字信号处理芯片和置于温控箱下的转台,光纤陀螺光学系统中的耦合器的输出端依次经探测器、信号放大器、模数转换器后与数字信号处理芯片连接,将光纤陀螺光学系统的光纤环置于温控环境的转台上,转台连接电机,由电机控制转台以特定角速度旋转控制,将光纤陀螺光学系统的非光纤环部分置于温控环境转台的外部。
所述的探测器探测光纤陀螺光学系统中耦合器的输出端光强,探测器信号经信号放大器后输入到模数转换器转换,数字信号处理芯片采集模数转换器信号,并连接传输到上位机。
本发明原理为,根据Sagnac干涉仪的原理可知,由角速度引起的两束光波相位差为:
式中,λ为光源的平均波长,c为真空中光速,D为线圈平均直径,L为光纤长度,△φR为由Ω转速产生的相位差。
干涉式光纤陀螺的干涉光强表示为:
I=I0(1+cos△φR)
若需产生2π相位差,则输入转速信号为:
由上可知在某一温度下,只需获得产生2π相位差的角速度输入即可获得该温度下光纤环特性,由此获得光纤环温变特性。
本发明的有益效果是:
本发明能测试获得光纤长度L、平均直径D乘积随温度变化的特性,可用于光纤陀螺光纤环温度性能评价及筛选,方法简单可靠,避免了非光纤环部分的干扰,为提高光纤陀螺标度因数温度稳定性提供基础。
本发明方法可对光纤陀螺光纤环结构温度性能进行测试评估,方法简单可靠,测量时间短,为光纤陀螺标度因数高温度稳定性的提高提供基础。
附图说明
图1是本发明测试装置的转台结构图。
图2是本发明测试装置的电气结构示意图。
图3是本发明测试方法的温度控制图。
图4是本发明测试方法的时序控制图。
图5是本发明测试方法测量两光纤环结构获得其光纤环参数随温度变化曲线。
图1中:1、温控箱,2、温箱内转台,3、光纤环,4、温控转台的电机,5、光纤陀螺光学系统的非光纤环部分,6、温箱外转台。
具体实施方式
下面结合附图和实例对本发明作进一步说明。
本发明的实施例如下:
实施的测量装置放置于温控转台中,如图1所示,数据采集及供电通过温控转台内部滑环进行连接。将光纤陀螺光学系统的光纤环3置于温控箱1中温控环境的转台2上,转台2连接电机4,由电机4控制转台2以特定角速度旋转控制,将光纤陀螺光学系统的非光纤环部分5置于温控环境外部的转台6,与转台2相同的转动情况,数据采集及供电通过温控转台内部滑环进行连接。
测量装置包括探测器、信号放大器、模数转换器、数字信号处理芯片和置于温控箱1下的转台2,如图2所示,光纤陀螺光学系统中的耦合器的输出端依次经探测器、信号放大器、模数转换器后与数字信号处理芯片连接,如图2所示。
如图2所示,探测器探测光纤陀螺光学系统中耦合器的输出端光强,探测器信号经信号放大器后输入到模数转换器转换,数字信号处理芯片采集模数转换器信号,并连接传输到上位机。
因需测试不同温度下光纤环特性,故采用如图3所示的温度控制,从室温降为至-40℃后开始第一次测试,随后分别测试-20℃、0℃、20℃、40℃、60℃条件下光纤环特性,如需获得更详细的光纤环温度特性,可将温度点设置的更为密集,测试时间也相应的更长。
在某一温度下进行测试时,转台由静止(角速度为0)开始以恒定加速度a加速,直至角速度为ω00>λc/LD,式中λ为光源平均波长,波动忽略不计,c为真空中光速,D为光纤环的线圈平均直径,L为光纤环的光纤长度。),随后恒定加速度a减速,直到角速度为-ω0;再以加速度a加速,直到角速度为0,并重复以上循环50次以上。
例如,针对光源平均波长为1310nm,光纤环长1km,平均直径0.1m的待测设备,设置以下参数:
ω0=300°/s(>3.93rad/s)
a=10°/s2
式中,ω0为最大角速度,a为设定的角加速度,T为执行一个循环所用时间。
根据转台精度可设置更小的角加速度,但会稍微增加测试时间。根据该测试流程,获得角加速度及角速度变化趋势如图4所示,用公式可以表示如下:
在某一固定温度下,可认为光纤环性能不变,故其产生的Sagnac相位差为
式中λ为光源平均波长,在采用相应温控电路的情况下,可认为波动较小,忽略不计,c为真空中的光速,L、D分别为光纤环长度和光纤环平均直径,△φR(t)为在设定测试流程下,Sagnac相位差变化。
如图4所示,△φR(t)峰峰值大于4π,以获得完整的周期信号。
根据干涉式光纤陀螺干涉光强公式可得,探测器信号呈现如图4所示。
设探测器信号中一个余弦相位的周期为Tt,则
ωT=aTt
式中T为测试环境温度,ωT即为该测试环境温度T下光强随角速度输入变化的周期,也即Sagnac相位差为2π对应的角速度输入。也即:
式中T表示不同测试环境温度,[LD]T为该测试环境温度T下光纤陀螺光纤环参数,λ为光源平均波长,在采用相应温控电路的情况下,可认为波动较小,忽略不计,c为真空中的光速。
对不同温度点进行测试即可获得光纤环温度特性。以下是对两种不同结构进行测试仿真计算流程及数据。
测试过程中,以加速度a=10°/s2,最大角速度为ω0=300°/s(>3.93rad/s)进行测试,分别测试-40℃、-20℃、0℃、20℃、40℃、60℃六个温度点。得到以下表格,(因使用的采样频率为2000Hz,故时间精度可以到万分之一左右)
表1
根据上述表格数据及前文所述公式进行计算得,两种不同结构光纤环参数如下表。
表2
将光纤陀螺光纤环参数[LD]T与温度的函数关系绘制成图,如图5所示。由图5及上表可知,光纤环结构2光纤环参数的温度温度性明显优于光纤环结构1
通过对不同结构或材料的光纤陀螺光纤环进行测试,可从中验证并获取最佳结构,为更高精度光纤陀螺打下基础。

Claims (6)

1.一种光纤陀螺光纤环温变特性测量方法,其特征在于该方法如下:将光纤陀螺的光纤环置于温控环境的转台上,将光纤陀螺非光纤环部分置于外部,在温控环境温度变化控制的同时控制转台进行特定角速度控制,采集光纤陀螺的耦合器输出的数据利用数据处理算法进行处理,获得光纤环温变特性。
2.根据权利要求1所述的一种光纤陀螺光纤环温变特性测量方法,其特征在于:所述采集光纤陀螺的耦合器输出的数据的时刻是温控环境内温度稳定的时刻。
3.根据权利要求1所述的一种光纤陀螺光纤环温变特性测量方法,其特征在于:所述的特定角速度控制是将转台从静止(角速度为0)开始,以恒定加速度a加速直至角速度为ω0,随后再以恒定加速度a减速直到角速度为-ω0,接着以加速度a加速直到角速度为0,以此重复循环50次以上;
上述角速度ω0满足:ω0>λc/LD,式中λ为光源平均波长,c为真空中光速,D为光纤环的线圈平均直径,L为光纤环的光纤长度。
4.根据权利要求1所述的一种光纤陀螺光纤环温变特性测量方法,其特征在于:所述的数据处理算法具体是:针对采集获得的耦合器输出数据,通过时域平均获得角速度从-ω0到ω0的时域平均数据,用时域平均数据进行拟合获得数据时域周期Tt,再采用以下公式计算获得光纤长度L和平均直径D乘积随温度变化的光纤环特征参数:
[ L D ] T = λ c aT t
式中,a为转台角加速度,Tt为探测光强的时域周期,λ为光源平均波长,波动忽略不计,c为真空中光速,D为光纤环的线圈平均直径,L为光纤环的光纤长度,T为测试温度,[LD]T为测试温度T下L与D的乘积,即光纤环特征参数。
5.用于实施权利要求1-4任一所述的一种光纤陀螺光纤环温变特性测量装置,其特征在于:所述测量装置连接在对光纤陀螺的光学系统上,包括探测器、信号放大器、模数转换器、数字信号处理芯片和置于温控箱(1)下的转台(2),光纤陀螺光学系统中的耦合器的输出端依次经探测器、信号放大器、模数转换器后与数字信号处理芯片连接,将光纤陀螺光学系统的光纤环(3)置于温控环境的转台(2)上,转台(2)连接电机(4),将光纤陀螺光学系统的非光纤环部分(5)置于温控环境外部的转台(6)。
6.根据权利要求5所述的一种光纤陀螺光纤环温变特性测量装置,其特征在于:所述的探测器探测光纤陀螺光学系统中耦合器的输出端光强,探测器信号经信号放大器后输入到模数转换器转换,数字信号处理芯片采集模数转换器信号,并连接传输到上位机。
CN201610933044.9A 2016-10-25 2016-10-25 一种光纤陀螺光纤环温变特性测量方法及装置 Active CN106441368B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610933044.9A CN106441368B (zh) 2016-10-25 2016-10-25 一种光纤陀螺光纤环温变特性测量方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610933044.9A CN106441368B (zh) 2016-10-25 2016-10-25 一种光纤陀螺光纤环温变特性测量方法及装置

Publications (2)

Publication Number Publication Date
CN106441368A true CN106441368A (zh) 2017-02-22
CN106441368B CN106441368B (zh) 2019-02-22

Family

ID=58177493

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610933044.9A Active CN106441368B (zh) 2016-10-25 2016-10-25 一种光纤陀螺光纤环温变特性测量方法及装置

Country Status (1)

Country Link
CN (1) CN106441368B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108151761A (zh) * 2017-12-05 2018-06-12 浙江大学 一种光纤环有效面积的实时闭环测量装置和方法
CN108630330A (zh) * 2018-05-29 2018-10-09 岭东核电有限公司 压水堆核电站仪表系统探测器试验处理方法、装置及系统
CN108692741A (zh) * 2018-04-04 2018-10-23 中航捷锐(北京)光电技术有限公司 一种闭环光纤陀螺标度因数快速稳定方法及装置
CN110595505A (zh) * 2019-09-18 2019-12-20 哈尔滨工程大学 一种测量光纤环温度特性的多参数仿真平台
CN111964659A (zh) * 2020-06-28 2020-11-20 北京航天时代光电科技有限公司 一种光纤陀螺光纤环温度测试与评价系统
CN113124899A (zh) * 2021-03-23 2021-07-16 西安航天精密机电研究所 一种基于仿真技术获取光纤陀螺变温标度因数的方法
CN113865574A (zh) * 2021-09-28 2021-12-31 中国船舶重工集团公司第七0七研究所 一种改善光纤陀螺标度因数线性度的系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101008569A (zh) * 2006-12-31 2007-08-01 北京航天控制仪器研究所 采用低偏和保偏混合光路的光纤陀螺
CN101387524A (zh) * 2008-10-09 2009-03-18 北京航空航天大学 一种适用于光纤陀螺的偏置温度误差测试与补偿系统
CN101285690B (zh) * 2008-05-26 2010-06-02 浙江大学 光纤陀螺磁场-温度灵敏度测试方法
CN101840529A (zh) * 2010-03-26 2010-09-22 东南大学 基于局部可变集成神经网络的光纤陀螺随机漂移建模方法
JP2012185046A (ja) * 2011-03-07 2012-09-27 Japan Aviation Electronics Industry Ltd 光ファイバジャイロ
CN104180798A (zh) * 2014-09-16 2014-12-03 中国科学院光电技术研究所 一种多光纤环串联的单轴光纤陀螺仪及多光纤环串联方法
CN105136165A (zh) * 2015-08-13 2015-12-09 北京航空航天大学 一种适用于光纤陀螺的定向梯度温度场响应特性测试装置
CN106017511A (zh) * 2016-07-28 2016-10-12 中国船舶重工集团公司第七0七研究所 一种光纤陀螺温度系数测试和补偿方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101008569A (zh) * 2006-12-31 2007-08-01 北京航天控制仪器研究所 采用低偏和保偏混合光路的光纤陀螺
CN101285690B (zh) * 2008-05-26 2010-06-02 浙江大学 光纤陀螺磁场-温度灵敏度测试方法
CN101387524A (zh) * 2008-10-09 2009-03-18 北京航空航天大学 一种适用于光纤陀螺的偏置温度误差测试与补偿系统
CN101840529A (zh) * 2010-03-26 2010-09-22 东南大学 基于局部可变集成神经网络的光纤陀螺随机漂移建模方法
JP2012185046A (ja) * 2011-03-07 2012-09-27 Japan Aviation Electronics Industry Ltd 光ファイバジャイロ
CN104180798A (zh) * 2014-09-16 2014-12-03 中国科学院光电技术研究所 一种多光纤环串联的单轴光纤陀螺仪及多光纤环串联方法
CN105136165A (zh) * 2015-08-13 2015-12-09 北京航空航天大学 一种适用于光纤陀螺的定向梯度温度场响应特性测试装置
CN106017511A (zh) * 2016-07-28 2016-10-12 中国船舶重工集团公司第七0七研究所 一种光纤陀螺温度系数测试和补偿方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108151761A (zh) * 2017-12-05 2018-06-12 浙江大学 一种光纤环有效面积的实时闭环测量装置和方法
CN108151761B (zh) * 2017-12-05 2020-05-19 浙江大学 一种光纤环有效面积的实时闭环测量装置及方法
CN108692741A (zh) * 2018-04-04 2018-10-23 中航捷锐(北京)光电技术有限公司 一种闭环光纤陀螺标度因数快速稳定方法及装置
CN108692741B (zh) * 2018-04-04 2020-09-25 中航捷锐(北京)光电技术有限公司 一种闭环光纤陀螺标度因数快速稳定方法及装置
CN108630330A (zh) * 2018-05-29 2018-10-09 岭东核电有限公司 压水堆核电站仪表系统探测器试验处理方法、装置及系统
CN110595505A (zh) * 2019-09-18 2019-12-20 哈尔滨工程大学 一种测量光纤环温度特性的多参数仿真平台
CN111964659A (zh) * 2020-06-28 2020-11-20 北京航天时代光电科技有限公司 一种光纤陀螺光纤环温度测试与评价系统
CN113124899A (zh) * 2021-03-23 2021-07-16 西安航天精密机电研究所 一种基于仿真技术获取光纤陀螺变温标度因数的方法
CN113124899B (zh) * 2021-03-23 2022-09-16 西安航天精密机电研究所 一种基于仿真技术获取光纤陀螺变温标度因数的方法
CN113865574A (zh) * 2021-09-28 2021-12-31 中国船舶重工集团公司第七0七研究所 一种改善光纤陀螺标度因数线性度的系统

Also Published As

Publication number Publication date
CN106441368B (zh) 2019-02-22

Similar Documents

Publication Publication Date Title
CN106441368B (zh) 一种光纤陀螺光纤环温变特性测量方法及装置
CN101339093B (zh) 光纤陀螺用光纤环质量的测量方法及其装置
CN102538822B (zh) 光纤陀螺仪的动态特性快速测试标定方法
CN104990547B (zh) 一种稳定光纤陀螺保持标度因数的方法和装置
CN101387524B (zh) 一种适用于光纤陀螺的偏置温度误差测试与补偿系统
CN107270884A (zh) 一种光纤陀螺及抑制输出零偏随温度变化漂移的方法
CN103115628B (zh) 一种谐振式光学陀螺标度因数测试方法
CN106643836B (zh) 一种同时测量轴向加速度与水平旋转角速度的光纤传感装置
CN106017511A (zh) 一种光纤陀螺温度系数测试和补偿方法
CN107525946A (zh) 一种基于光学波导中原子干涉的加速度测量方法及装置
CN104713575A (zh) 一种闭环光纤陀螺频率特性的测试方法
CN103791917B (zh) 光纤陀螺振动-温度交联耦合度的测试方法
CN103439530A (zh) 光学加速度计
Zhang et al. Analysis for the thermal performance of a modified quadrupolar fiber coil
Celikel et al. Establishment of all digital closed-loop interferometric fiber-optic gyroscope and scale factor comparison for open-loop and all digital closed-loop configurations
CN209623725U (zh) 一种光纤线圈动态仿真测试装置
CN104457792B (zh) 一种在无机械转动条件下测量光纤陀螺标度因数的方法
CN106643791A (zh) 一种光纤陀螺反馈回路性能测试方法及装置
CN109946042A (zh) 一种y型波导耦合精度的测量方法及装置
CN105444777A (zh) 一种摇摆条件下光纤陀螺仪误差测试方法
CN104457791B (zh) 一种静态条件下测量光纤陀螺带宽的方法
CN114216450B (zh) 一种光纤陀螺的标度因数温度补偿方法及系统
CN115452196A (zh) 一种光纤敏感环高精度温度灵敏度系数测试的装置及方法
CN110186500A (zh) 一种采用绝对法的非平衡光纤干涉仪臂长差测量装置及测量方法
CN109489686A (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