CN109256995A - 感应电机定子电阻启动辨识方法 - Google Patents
感应电机定子电阻启动辨识方法 Download PDFInfo
- Publication number
- CN109256995A CN109256995A CN201811221258.9A CN201811221258A CN109256995A CN 109256995 A CN109256995 A CN 109256995A CN 201811221258 A CN201811221258 A CN 201811221258A CN 109256995 A CN109256995 A CN 109256995A
- Authority
- CN
- China
- Prior art keywords
- current
- stator resistance
- voltage
- follows
- resistance
- 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
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
- Motor And Converter Starters (AREA)
Abstract
本发明公开了感应电机定子电阻启动辨识方法,分别通过系统建模、PI参数整定、计算电流相应和电压响应计算出每次启动时定子电阻相对于离线辨识定子电阻变化的增量ΔRs再通过最终公式Rs1=Rs+ΔRs,得出每次启动时的定子电阻Rs1,本方法提供的感应电机定子电阻启动辨识方法简单快速,且在最大程度上消除由于温度变化带来的低频出力不足和稳定性问题。
Description
技术领域
本发明涉及电机控制技术领域,特别涉及感应电机定子电阻启动辨识方法。
背景技术
异步电机无速度传感器控制系统需要精确的磁链信息,电压型磁链观测器是实现电机定转子磁链估算最简单的方法,其公式为:ψs=∫(us-Ris)dt。
定子电阻是该模型中唯一涉及的电机参数,直接影响到磁场定向的准确性,关系到低频时电机带载能力和矢量控制系统的稳定性。然而,在电机运转中,由于电机绕组发热以及环境温度的变化,定子电阻将发生变化,从而影响控制性能。
目前工业场合变频器驱动异步电机常用的定子电阻辨识方法有两种,离线辨识和在线辨识。离线辨识一般包含在参数自学习功能中,该功能会对所有电机参数进行自学习。其中学习定子电阻时采用电流闭环阶跃响应测试,在电压达到稳态时用伏安法计算定子电阻,电压稳定时间和电机功率相关。一般功率越大,所需时间越长。对于大功率机器可能需要数秒的时间。而在工业应用中,由于参数自学习功能耗时较长,因此不会在每次运行前都进行参数自学习。因此,由于环境温度的变化导致的定子电阻变化的问题依旧没有得到解决。如果只将定子电阻学习这一部分照搬到电机启动前进行,由于需要数秒的时间,往往难以忍受。并且对于一些特殊的场合,比如升降机设备,这个过程可能导致设备的降落。
在线辨识一般是基于电机模型的方法,需进行观测器设计,算法复杂,计算量大,工程实现较困难。
发明内容
本发明解决的技术问题是提供一种能最大程度消除由于温度变化带来的低频出力不足和稳定性问题的感应电机定子电阻启动辨识方法。
本发明解决其技术问题所采用的技术方案是:感应电机定子电阻启动辨识方法,其步骤为:
A、系统建模:在静止坐标系下,建立电流闭环控制定子电阻辨识,其中包含PI环节,逆变器环节,电机环节以及电流采样环节,各环节模型如下所示:
(1)建立感应电机静止状态下的模型,包括依次连接的定子电阻Rs,电机定子漏感Lls,电机转子漏感Llr和电阻Rr,还包括连接在电感Llr和电阻Rr之间的电机互感Lm,还包括施加的电阻Rs输入端和电阻Rr输出端之间的电压u,其输入电流为i;
其中可近似认为近似认为Lls=Llr,则定转子全电感均为L=Lm+Lls,定子时间常数转子时间常数漏感系数
定子电流和电压之间的传递函数为:
所述s代表微分算子;
(2)考虑逆变器延迟和死区的SVPWM逆变器传递函数为:
其中T∑i为跟载波频率相关的小时间常数,一般小于1ms;
(3)电流采样采用硬件有源低通滤波器,滤波时间远小于T∑i,因此环节等效为1;
B、PI参数整定:
从到isα的开环传函为:
所述为给定α轴参考电流,isα为反馈α轴电流;
该系统为I型系统,因此阶跃响应无稳态误差,
电机模型的极点:
考虑到(Ts+Tr)2>>4σTrTs,因此很显然,s1>>s2。
电机模型为带零点的二阶系统,三个转折频率分别为:
显然ω1<<ω3,ω2<<ω3,ω1和ω2构成偶极子且不十分靠近坐标极点,因此对系统的瞬态响应很小,忽略后,从到isα的开环传函近似为:
由于用PI零点对消掉大惯性环节极点,按照二阶“最优”系统整定,得到
C、电流响应和电压响应:
按照上述方法整定出PI参数后,下面来推导在电流阶跃指令给定下,实际电流和电压表达式从而获得定子电阻启动辨识方法:
整定PI参数时,需要考虑高频段,忽略了偶极子下面推导低频段,也即稳态时电压响应表达式,此时可忽略掉项,系统开环传函从(2)变为:
从而系统闭环传递函数为:
式(8)特征方程明显有两个实根。考虑到TΣi<<Ts,TΣi<<Tr,则:
其中,表示是和同阶的数,一般小于1ms2了;
从而,在电流单位阶跃指令给定下,isα(t)时域表达式为:
从(9)可以看出,第一项i1(t)=1为稳态分量,第二项的时间常数近似为Tr,第三项的时间常数近似为2TΣi;
电流阶跃响应时,暂态分量前几个PWM周期由i3(t)主导,幅值大,但持续时间短,经过10TΣi此分量就衰减到0.7%,后面主要由i2(t)主导,幅值小但持续时间长,此分量幅值最大为因此可以认为i2(t)主导暂态响应时电流已经稳定,也即电流阶跃响应在10TΣi达到稳定,
接着可以得到电压在电流单位阶跃响应下时域表达式为:
从(10)可以看出,第一项u1(t)=Rs为稳态分量,第二项的时间常数近似为Tr,第三项的时间常数近似为2TΣi;
与电流阶跃响应表达式对比可知,电压暂态响应与电流具有相同的时间常数,电压暂态分量前几个PWM周期由u3(t)主导,幅值大,但持续时间短,经过10TΣi此分量衰减到0.7%Kp,经过20TΣi此分量衰减到与u2(t)相比已经可忽略,后面主要由u2(t)主导,幅值最大为Rr,以转子时间常数Tr指数衰减,相比稳态分量u1(t)不可忽略;
对电机离线辨识结束后以及每次启动前两个不同时刻进行电流阶跃控制,20TΣi后的采样电压响应曲线如图7,其表达式为:
则由式(12)与式(11)相减得:
Δu=Rs1-Rs=ΔRs
ΔRs即为每次启动时定子电阻相对于离线辨识定子电阻变化的增量,只要知道离线辨识时的定子电阻值Rs,每次启动时的定子电阻Rs1可通过下述公式计算得到:
Rs1=Rs+ΔRs(13)。
本发明的有益效果是:本发明提供的感应电机定子电阻启动辨识方法简单快速,且在最大程度上消除由于温度变化带来的低频出力不足和稳定性问题。
附图说明
图1为基于静止坐标系下的电流闭环控制定子电阻辨识控制框图;
图2感应电机电机静止状态下的等效电路;
图3为图1模型的传递函数框图;
图4电压响应曲线;
图5电压响应曲线中u3(t)占主导部分;
图6电压响应曲线中u2(t)占主导部分;
图7不同时刻电流阶跃时电压响应曲线;
具体实施方式
下面结合附图和具体实施方式对本发明进一步说明。
感应电机定子电阻启动辨识方法,其步骤为:
A、系统建模:在静止坐标系下,建立电流闭环控制定子电阻辨识框图如图1所示,其中包含PI环节,逆变器环节,电机环节以及电流采样环节,各环节模型如下所示:
(1)建立感应电机静止状态下的模型,包括依次连接的定子电阻Rs,电机定子漏感Lls,电机转子漏感Llr和电阻Rr,还包括连接在电感Llr和电阻Rr之间的电机互感Lm,还包括施加的电阻Rs输入端和电阻Rr输出端之间的电压u,其输入电流为i;
其中可近似认为近似认为Lls=Llr,则定转子全电感均为L=Lm+Lls,定子时间常数转子时间常数漏感系数
定子电流和电压之间的传递函数为:
(2)考虑逆变器延迟和死区的SVPWM逆变器传递函数为:
其中T∑i为跟载波频率相关的小时间常数,一般小于1ms;
(3)电流采样采用硬件有源低通滤波器,滤波时间远小于T∑i,因此环节等效为1;
综合以上环节,可得框图1模型的传递函数框图如图3
B、PI参数整定:
从到isα的开环传函为:
该系统为I型系统,因此阶跃响应无稳态误差,
电机模型的极点:
考虑到(Ts+Tr)2>>4σTrTs,因此很显然,s1>>s2。
电机模型为带零点的二阶系统,三个转折频率分别为:
显然ω1<<ω3,ω2<<ω3,ω1和ω2构成偶极子且不十分靠近坐标极点,因此对系统的瞬态响应很小,忽略后,从到isα的开环传函近似为:
由于用PI零点对消掉大惯性环节极点,按照二阶“最优”系统整定,得到
C、电流响应和电压响应:
按照上述方法整定出PI参数后,下面来推导在电流阶跃指令给定下,实际电流和电压表达式从而获得定子电阻启动辨识方法:
整定PI参数时,需要考虑高频段,忽略了偶极子下面推导低频段,也即稳态时电压响应表达式,此时可忽略掉项,系统开环传函从(2)变为:
从而系统闭环传递函数为:
式(8)特征方程明显有两个实根。考虑到TΣi<<Ts,TΣi<<Tr,则:
其中,表示是和同阶的数,一般小于1ms2了;
从而,在电流单位阶跃指令给定下,isα(t)时域表达式为:
从(9)可以看出,第一项i1(t)=1为稳态分量,第二项的时间常数近似为Tr,第三项的时间常数近似为2TΣi;
电流阶跃响应时,暂态分量前几个PWM周期由i3(t)主导,幅值大,但持续时间短,经过10TΣi此分量就衰减到0.7%,后面主要由i2(t)主导,幅值小但持续时间长,此分量幅值最大为因此可以认为i2(t)主导暂态响应时电流已经稳定,也即电流阶跃响应在10TΣi达到稳定,
接着可以得到电压在电流单位阶跃响应下时域表达式为:
从(10)可以看出,第一项u1(t)=Rs为稳态分量,第二项的时间常数近似为Tr,如图6所示,第三项的时间常数近似为2TΣi,如图5所示;
与电流阶跃响应表达式对比可知,电压暂态响应与电流具有相同的时间常数,电压暂态分量前几个PWM周期由u3(t)主导,幅值大,但持续时间短,经过10TΣi此分量衰减到0.7%Kp,经过20TΣi此分量衰减到与u2(t)相比已经可忽略,后面主要由u2(t)主导,幅值最大为Rr,以转子时间常数Tr指数衰减,相比稳态分量u1(t)不可忽略;
对电机离线辨识结束后以及每次启动前两个不同时刻进行电流阶跃控制,20TΣi后的采样电压响应曲线如图7所示,其表达式为:
则由式(12)与式(11)相减得:
Δu=Rs1-Rs=ΔRs
ΔRs即为每次启动时定子电阻相对于离线辨识定子电阻变化的增量,只要知道离线辨识时的定子电阻值Rs,每次启动时的定子电阻Rs1可通过下述公式计算得到:
Rs1=Rs+ΔRs (13)。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.感应电机定子电阻启动辨识方法,其步骤为:
A、系统建模:在静止坐标系下,建立电流闭环控制定子电阻辨识,其中包含PI环节,逆变器环节,电机环节以及电流采样环节,各环节模型如下所示:
(1)建立感应电机静止状态下的模型,包括依次连接的定子电阻Rs,电机定子漏感Lls,电机转子漏感Llr和电阻Rr,还包括连接在电感Llr和电阻Rr之间的电机互感Lm,还包括施加的电阻Rs输入端和电阻Rr输出端之间的电压u,其输入电流为i;
其中可近似认为近似认为Lls=Llr,则定转子全电感均为L=Lm+Lls,定子时间常数转子时间常数漏感系数
定子电流和电压之间的传递函数为:
所述s代表微分算子;
(2)考虑逆变器延迟和死区的SVPWM逆变器传递函数为:
其中T∑i为跟载波频率相关的小时间常数,一般小于1ms;
(3)电流采样采用硬件有源低通滤波器,滤波时间远小于T∑i,因此环节等效为1;B、PI参数整定:
从到isα的开环传函为:
所述为给定α轴参考电流,isα为反馈α轴电流;
该系统为I型系统,因此阶跃响应无稳态误差,
电机模型的极点:
考虑到(Ts+Tr)2>>4σTrTs,因此很显然,s1>>s2。
电机模型为带零点的二阶系统,三个转折频率分别为:
显然ω1<<ω3,ω2<<ω3,ω1和ω2构成偶极子且不十分靠近坐标极点,因此对系统的瞬态响应很小,忽略后,从到isα的开环传函近似为:
由于用PI零点对消掉大惯性环节极点,按照二阶“最优”系统整定,得到
C、电流响应和电压响应:
按照上述方法整定出PI参数后,下面来推导在电流阶跃指令给定下,实际电流和电压表达式从而获得定子电阻启动辨识方法:
整定PI参数时,需要考虑高频段,忽略了偶极子下面推导低频段,也即稳态时电压响应表达式,此时可忽略掉项,系统开环传函从(2)变为:
从而系统闭环传递函数为:
式(8)特征方程明显有两个实根。考虑到TΣi<<Ts,TΣi<<Tr,则:
其中,表示是和同阶的数,一般小于1ms2了;
从而,在电流单位阶跃指令给定下,isα(t)时域表达式为:
从(9)可以看出,第一项i1(t)=1为稳态分量,第二项的时间常数近似为Tr,第三项的时间常数近似为2TΣi;
电流阶跃响应时,暂态分量前几个PWM周期由i3(t)主导,幅值大,但持续时间短,经过10TΣi此分量就衰减到0.7%,后面主要由i2(t)主导,幅值小但持续时间长,此分量幅值最大为因此可以认为i2(t)主导暂态响应时电流已经稳定,也即电流阶跃响应在10TΣi达到稳定,
接着可以得到电压在电流单位阶跃响应下时域表达式为:
从(10)可以看出,第一项u1(t)=Rs为稳态分量,第二项的时间常数近似为Tr,第三项的时间常数近似为2TΣi;
与电流阶跃响应表达式对比可知,电压暂态响应与电流具有相同的时间常数,电压暂态分量前几个PWM周期由u3(t)主导,幅值大,但持续时间短,经过10TΣi此分量衰减到0.7%Kp,经过20TΣi此分量衰减到0.45%00Kp,与u2(t)相比已经可忽略,后面主要由u2(t)主导,幅值最大为Rr,以转子时间常数Tr指数衰减,相比稳态分量u1(t)不可忽略;
对电机离线辨识结束后以及每次启动前两个不同时刻进行电流阶跃控制,20TΣi后的采样电压响应曲线表达式为:
则由式(12)与式(11)相减得:
Δu=Rs1-Rs=ΔRs
ΔRs即为每次启动时定子电阻相对于离线辨识定子电阻变化的增量,只要知道离线辨识时的定子电阻值Rs,每次启动时的定子电阻Rs1可通过下述公式计算得到:
Rs1=Rs+ΔRs (13)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811221258.9A CN109256995B (zh) | 2018-10-19 | 2018-10-19 | 感应电机定子电阻启动辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811221258.9A CN109256995B (zh) | 2018-10-19 | 2018-10-19 | 感应电机定子电阻启动辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109256995A true CN109256995A (zh) | 2019-01-22 |
CN109256995B CN109256995B (zh) | 2021-09-24 |
Family
ID=65046050
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811221258.9A Active CN109256995B (zh) | 2018-10-19 | 2018-10-19 | 感应电机定子电阻启动辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109256995B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110635738A (zh) * | 2019-10-12 | 2019-12-31 | 东南大学 | 一种永磁同步电机定子电阻及电机温度的实时辨识方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393507A (zh) * | 2011-09-01 | 2012-03-28 | 北京配天大富精密机械有限公司 | 一种电机参数检测方法及电机参数检测装置 |
CN103825524A (zh) * | 2014-03-14 | 2014-05-28 | 中冶南方(武汉)自动化有限公司 | 一种永磁同步电机基本电器参数的离线辨识方法 |
CN103869172A (zh) * | 2012-12-10 | 2014-06-18 | 台达电子企业管理(上海)有限公司 | 感应电机的转子电阻测量方法及测量装置 |
US20140241032A1 (en) * | 2013-02-26 | 2014-08-28 | Seagate Technology Llc | Methods and apparatuses using a transfer function to predict resistance shifts and/or noise of resistance-based memory |
CN104037795A (zh) * | 2014-06-27 | 2014-09-10 | 武汉大学 | 一种vsc-hvdc控制系统的pi参数整定方法 |
CN106773710A (zh) * | 2017-01-13 | 2017-05-31 | 东南大学 | 一种直流电机电流控制器参数自整定系统及方法 |
CN107196569A (zh) * | 2017-02-28 | 2017-09-22 | 常州联力自动化科技有限公司 | 一种基于dsp的转速估计环节pi参数定量整定方法 |
-
2018
- 2018-10-19 CN CN201811221258.9A patent/CN109256995B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393507A (zh) * | 2011-09-01 | 2012-03-28 | 北京配天大富精密机械有限公司 | 一种电机参数检测方法及电机参数检测装置 |
CN103869172A (zh) * | 2012-12-10 | 2014-06-18 | 台达电子企业管理(上海)有限公司 | 感应电机的转子电阻测量方法及测量装置 |
US20140241032A1 (en) * | 2013-02-26 | 2014-08-28 | Seagate Technology Llc | Methods and apparatuses using a transfer function to predict resistance shifts and/or noise of resistance-based memory |
CN103825524A (zh) * | 2014-03-14 | 2014-05-28 | 中冶南方(武汉)自动化有限公司 | 一种永磁同步电机基本电器参数的离线辨识方法 |
CN104037795A (zh) * | 2014-06-27 | 2014-09-10 | 武汉大学 | 一种vsc-hvdc控制系统的pi参数整定方法 |
CN106773710A (zh) * | 2017-01-13 | 2017-05-31 | 东南大学 | 一种直流电机电流控制器参数自整定系统及方法 |
CN107196569A (zh) * | 2017-02-28 | 2017-09-22 | 常州联力自动化科技有限公司 | 一种基于dsp的转速估计环节pi参数定量整定方法 |
Non-Patent Citations (2)
Title |
---|
R.N. CONTE ETAL.: "Parameters identification of induction machines based on frequency response and optimization techniques", 《IECON"03. 29TH ANNUAL CONFERENCE OF THE IEEE INDUSTRIAL ELECTRONICS SOCIETY (IEEE CAT. NO.03CH37468)》 * |
杨宗军 等: "表贴式永磁同步电机的多参数在线辨识", 《电工技术学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110635738A (zh) * | 2019-10-12 | 2019-12-31 | 东南大学 | 一种永磁同步电机定子电阻及电机温度的实时辨识方法 |
CN110635738B (zh) * | 2019-10-12 | 2020-12-18 | 东南大学溧阳研究院 | 一种永磁同步电机定子电阻及电机温度的实时辨识方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109256995B (zh) | 2021-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108183648B (zh) | 一种基于逆变器非线性补偿的永磁同步电机参数辨识方法 | |
Soto et al. | Reduced-order observers for rotor flux, rotor resistance and speed estimation for vector controlled induction motor drives using the extended Kalman filter technique | |
Solvar et al. | Sensorless second order sliding mode observer for induction motor | |
CN107482982A (zh) | 一种基于铁损模型的异步电机矢量控制方法 | |
CN110609230A (zh) | 一种永磁同步电机的功率级模拟驱动系统及方法 | |
CN103501150A (zh) | 一种内嵌式永磁同步电机参数辨识装置及方法 | |
CN111585488B (zh) | 一种永磁电机无速度传感器控制方法及系统 | |
Jannati et al. | Speed sensorless fault-tolerant drive system of 3-phase induction motor using switching extended kalman filter | |
CN113422550B (zh) | 基于复矢量解耦及延时补偿的高速电机低载波比控制方法 | |
CN109256995A (zh) | 感应电机定子电阻启动辨识方法 | |
CN106817056A (zh) | 基于谐波观测器的电机驱动系统 | |
Antonello et al. | Theory and implementation of a MTPA tracking controller for anisotropic PM motor drives | |
Xu et al. | Online estimation of the stator resistances of a PMSM with consideration of magnetic saturation | |
Wang et al. | Simplified rotor and stator resistance estimation method based on direct rotor flux identification | |
Nikmaram et al. | Comparison study of active flux based sliding-mode observer and pll based sliding-mode observer sensorless control of synrm | |
CN108683376A (zh) | 一种基于滤波原理的交流电机电流测量偏移误差补偿方法 | |
Sun et al. | Hybrid excitation synchronous motor feedback linearization decoupling sliding mode control | |
Hu et al. | On-line Parameter Identification of Permanent Magnet Synchronous Motor based on Extended Kalman Filter | |
CN109742992B (zh) | 一种无位置传感器的新型混合励磁同步电机控制方法 | |
Renukrishna et al. | Sensorless vector control of induction motor drives using rotor flux observer | |
CN109842337B (zh) | 一种磁链差异量处理方法及电机控制装置 | |
Devadason et al. | Analysis of Free Acceleration and Steady State Characteristics of Line-Started Induction Motors with Shunt Capacitor Compensation | |
CN117856691B (zh) | 一种永磁同步电机鲁棒模型预测速度控制方法 | |
Fang et al. | Online identification of flux linkage of PMSM with unbalanced stator resistance | |
CN115333424A (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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 215000 No. 1000 Songhua Road, Guo Xiang street, Wuzhong Economic and Technological Development Zone, Suzhou, Jiangsu Applicant after: Suzhou Weichuang Electrical Technology Co., Ltd. Address before: 215000 No. 1000 Songhua Road, Guo Xiang street, Wuzhong Economic and Technological Development Zone, Suzhou, Jiangsu Applicant before: SUZHOU WEICHUANG ELECTRICAL EQUIPMENT TECHNOLOGY CO., LTD. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |