CN103678775B - 一种输送管动强度分析方法 - Google Patents
一种输送管动强度分析方法 Download PDFInfo
- Publication number
- CN103678775B CN103678775B CN201310577106.3A CN201310577106A CN103678775B CN 103678775 B CN103678775 B CN 103678775B CN 201310577106 A CN201310577106 A CN 201310577106A CN 103678775 B CN103678775 B CN 103678775B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- conveying pipe
- mfrac
- msup
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000003068 static effect Effects 0.000 claims abstract description 34
- 230000006378 damage Effects 0.000 claims abstract description 28
- 230000003595 spectral effect Effects 0.000 claims abstract description 18
- 238000001228 spectrum Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 239000000463 material Substances 0.000 claims description 22
- 239000003380 propellant Substances 0.000 claims description 14
- 230000005284 excitation Effects 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 abstract description 5
- 208000027418 Wounds and injury Diseases 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 2
- 208000014674 injury Diseases 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 235000015842 Hesperis Nutrition 0.000 description 1
- 235000012633 Iberis amara Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Pipeline Systems (AREA)
Abstract
一种输送管动强度分析方法,本发明由静强度分析、模态分析、随机响应分析和疲劳损伤计算评估六部分组成,首先通过有限元软件Abaqus的static general模块、Frequency模块、Random reaponse模块依次实现输送管的静强度分析、模态分析和随机响应分析,然后采用Mises应力谱密度估算方法计算出Mises应力谱密度,基于随机振动的功率谱密度法,得到输送管的定量疲劳损伤后进行评估。本发明能够评价输送管动强度是否满足要求,在管路产品设计时提前揭示风险,提高了输送管动强度评价效率。
Description
技术领域
本发明涉及一种适用于输送管的动强度分析方法,属于运载火箭力学环境技术领域。
背景技术
疲劳失效是弹(箭)体结构件最常见的破坏模式之一,随着在研运载型号推力的提升,输送管特别是发动机舱内的输送管振动量级相应增加,载荷环境比较恶劣,动强度成为输送管是否满足使用要求的关键指标。
传统的疲劳寿命估算方法是通过试验获得危险点的应力—时间历程,采用不同的计数方法,将不规则应力时间历程转化为不同应力幅下的循环次数,再利用疲劳累积损伤理论估算疲劳损伤及寿命。但对于输送管,由于缺乏试验条件,无法获得危险点的应力时间历程,因此只能采用基于频域的统计方法来估算疲劳寿命。
目前基于频域的统计方法仅针对振动载荷下的动强度进行分析,未考虑输送管承受内压、温度及位移等静载荷对动强度的影响;仅针对实体结构的动强度进行分析,未考虑输送管内部携带推进剂重量的影响;采用应力分量功率谱密度作为参量进行疲劳寿命估算,则输送管处于多轴应力状态下时,忽略了其他应力分量的影响。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种输送管动强度分析方法,本方法通过采集频率域内的主要信息-功率谱密度,根据频率域内载荷过程的统计量分析完成输送管路的强度分析,本发明解决了管路疲劳寿命的评估标准问题,填补了管路设计领域在该专业方面的空白。
本发明的技术解决方案是:
一种输送管动强度分析方法包括步骤如下:
(1)首先将推进剂重量等效处理后,在有限元软件Abaqus中对输送管进行建模,所述的输送管道建模包括输送管的几何属性和材料属性;
(2)判断输送管模型几何和材料的正确性,若不正确返回步骤(1),若正确则对模型施加内压、位移和温度载荷,利用有限元软件Abaqus(采用static general模块)进行静强度分析;
(3)判断输送管静强度是否满足要求,若不满足要求,改进输送管道结构返回步骤(1),若满足要求,则利用有限元软件Abaqus(采用Frequency模块)进行模态分析;
(4)利用有限元软件Abaqus(采用Random response模块),结合步骤(3)中分析得到的模态以及对模型施加随机激励(例如:加速度功率谱密度)进行随机响应分析得到的应力分量均方根及功率谱密度;
(5)对步骤(4)中分析得到的应力分量均方根及功率谱密度,根据Mises应力功率谱密度估算方法计算Mises应力功率谱密度;
(6)根据步骤(5)中计算得到的Mises应力功率谱密度,由Dirlik公式计算出输送管的疲劳损伤;
(7)判断输送管疲劳损伤是否满足要求,若不满足要求,改进输送管道结构后返回步骤(1),若满足要求,则结束。
所述步骤(1)中推进剂重量等效处理包括如下步骤:
输送管等效密度ρ通过公式得到:
其中,D为输送管路内径,δ为输送管路壁厚,ρ1为输送管材密度,ρ2为推进剂密度。
所述步骤(3)中判断输送管静强度是否满足要求包括如下步骤:
(a)输送管静强度安全因子η通过如下公式得到:
其中,σ为静强度分析得到的输送管Mises应力,σs为管路材料在设计温度下的屈服强度;
(b)若η≥1,则输送管静强度满足要求,反之若η<1不满足要求。
所述步骤(5)中Mises应力功率谱密度估算方法如下:
输送管Mises应力谱密度G(σM)通过如下公式得到:
其中,G(σ11)、G(σ22)、G(σ12)为随机响应分析得到的输送管各应力分量的谱密度。
所述步骤(6)中疲劳损伤计算方法如下:
输送管疲劳损伤DL通过如下公式得到:
其中,C、m为材料疲劳曲线参数,σb为管路材料在设计温度下的强度极限,T为输送管振动时间,p(S)为Mises应力概率密度,E[p]为单位之间内的峰值数。
所述单位之间内的峰值数E[p]通过如下公式得到:
其中,f1为随机激励的频率下限,f2为随机激励的频率上限。
所述Mises应力概率密度p(S)通过公式:
其中,
所述步骤(7)中判断输送管疲劳损伤是否满足要求方法如下:
若输送管疲劳损伤DL≥0.1,则输送管疲劳损伤不满足要求,反之若DL<0.1满足要求。
本发明与现有技术相比的有益效果是:
(1)本发明给出了由静强度分析、模态分析、随机响应分析和疲劳损伤计算评估六部分依次组成的动强度分析流程,采用有限元软件Abaqus的static general模块、Frequency模块、Random reaponse模块实现输送管的静强度分析、模态分析和随机响应分析;通过该流程,内压、温度及位移等静载荷对动强度的影响将通过静强度分析向模态分析、随机响应分析和疲劳损伤计算评估传递,避免了单纯振动载荷下动强度分析的片面性;
(2)本发明对推进剂重量进行等效,可以在动强度分析中引入推进剂重量的影响。等效处理的好处是,既灵活的考虑了推进剂重量,又避免了需要对推进剂进行额外的建模、计算工作,使计算效率得到了保证。
(3)本发明给出了Mises应力谱密度与各应力分量应力谱密度的换算关系,采用Mises应力谱密度作为疲劳寿命估算参量,较之单一应力分量谱密度,充分考虑了输送管多轴应力状态的影响,可以更全面、更精确的评估输送管的疲劳寿命。
附图说明
图1为本发明方法流程图。
具体实施方式
下面结合附图对本发明的具体实施方式进行进一步的详细描述。
如图1所示,本发明一种输送管动强度分析方法,由静强度分析、模态分析、随机响应分析和疲劳损伤计算评估六部分组成,具体步骤如下:
(1)首先将推进剂重量等效处理后,在有限元软件Abaqus中对输送管进行建模,所述的输送管道建模包括输送管的几何属性和材料属性;
推进剂重量等效处理包括如下步骤:
输送管等效密度ρ通过公式得到:
其中,D为输送管路内径,δ为输送管路壁厚,ρ1为输送管材密度,ρ2为推进剂密度。
(2)判断输送管模型几何属性和材料属性的正确性(几何属性和材料属性的要求根据用户实际的需求相一致),若不正确返回步骤(1),若正确则对模型施加内压、位移和温度载荷,利用有限元软件Abaqus(采用有限元软件Abaqus中的static general模块)进行静强度分析;
(3)判断输送管静强度是否满足要求,若不满足要求,改进输送管道结构返回步骤(1),若满足要求,则利用有限元软件Abaqus(采用有限元软件Abaqus中的Frequency模块,截断频率阶数建议取值为20,小于等于20阶,采用subspace求解器;大于20阶,采用Lanczos求解器)进行模态分析;
判断输送管静强度是否满足要求包括如下步骤:
(a)输送管静强度安全因子η通过如下公式得到:
其中,σ为静强度分析得到的输送管Mises应力,σs为管路材料在设计温度下的屈服强度;
(b)若η≥1,则输送管静强度满足要求,反之若η<1不满足要求。
(4)利用有限元软件Abaqus(采用有限元软件Abaqus中的Random response模块,阻尼系数建议取值为0.03),结合步骤(3)中分析得到的模态以及对模型施加随机激励(例如:加速度功率谱密度)进行输送管道的随机响应分析得到应力分量均方根及功率谱密度;
(5)对步骤(4)中分析得到的应力分量均方根及功率谱密度,根据Mises应力功率谱密度估算方法计算Mises应力功率谱密度;
Mises应力功率谱密度估算方法如下:
输送管Mises应力谱密度G(σM)通过如下公式得到:
其中,G(σ11)、G(σ22)、G(σ12)为随机响应分析得到的输送管各应力分量的谱密度。
(6)根据步骤(5)中计算得到的Mises应力功率谱密度,由Dirlik公式计算出输送管的疲劳损伤;
疲劳损伤计算方法如下:
输送管疲劳损伤DL通过如下公式得到:
其中,C、m为材料疲劳曲线参数,σb为管路材料在设计温度下的强度极限,T为输送管振动时间,p(S)为Mises应力概率密度,E[p]为单位之间内的峰值数。
单位之间内的峰值数E[p]通过如下公式得到:
其中,f1为随机激励的频率下限,f2为随机激励的频率上限。
Mises应力概率密度p(S)通过公式:
其中,
(7)判断输送管疲劳损伤是否满足要求,若不满足要求,改进输送管道结构后返回步骤(1),若满足要求,则结束。
判断输送管疲劳损伤是否满足要求方法如下:
若输送管疲劳损伤DL≥0.1,则输送管疲劳损伤不满足要求,反之若DL<0.1满足要求。
本发明说明书中未作详细描述的内容属于本领域技术人员的公知技术。
Claims (1)
1.一种输送管动强度分析方法,其特征在于步骤如下:
(1)首先将推进剂重量等效处理后,在有限元软件Abaqus中对输送管进行建模,所述的输送管道建模包括输送管的几何属性和材料属性;
(2)判断输送管模型几何和材料的正确性,若不正确返回步骤(1),若正确则对模型施加内压、位移和温度载荷,利用有限元软件Abaqus进行静强度分析;
(3)判断输送管静强度是否满足要求,若不满足要求,改进输送管道结构返回步骤(1),若满足要求,则利用有限元软件Abaqus进行模态分析;
(4)利用有限元软件Abaqus,并结合步骤(3)中分析得到的模态以及对模型施加随机激励进行随机响应分析得到的应力分量均方根及功率谱密度;
(5)对步骤(4)中分析得到的应力分量均方根及功率谱密度,根据Mises应力功率谱密度估算方法计算Mises应力功率谱密度;
(6)根据步骤(5)中计算得到的Mises应力功率谱密度,由Dirlik公式计算出输送管的疲劳损伤;
(7)判断输送管疲劳损伤是否满足要求,若不满足要求,改进输送管道结构后返回步骤(1),若满足要求,则结束;
所述步骤(1)中推进剂重量等效处理包括如下步骤:
输送管等效密度ρ通过公式得到:
<mrow>
<mi>&rho;</mi>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mfrac>
<mi>D</mi>
<mrow>
<mn>4</mn>
<mi>&delta;</mi>
</mrow>
</mfrac>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
</mrow>
其中,D为输送管路内径,δ为输送管路壁厚,ρ1为输送管材密度,ρ2为推进剂密度;
所述步骤(3)中判断输送管静强度是否满足要求包括如下步骤:
(a)输送管静强度安全因子η通过如下公式得到:
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>1.67</mn>
<msub>
<mi>&sigma;</mi>
<mi>s</mi>
</msub>
</mrow>
<mi>&sigma;</mi>
</mfrac>
</mrow>
其中,σ为静强度分析得到的输送管Mises应力,σs为管路材料在设计温度下的屈服强度;
(b)若η≥1,则输送管静强度满足要求,反之若η<1不满足要求;
所述步骤(5)中Mises应力功率谱密度估算方法如下:
输送管Mises应力谱密度G(σM)通过如下公式得到:
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mi>M</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>11</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>22</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msqrt>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>11</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>22</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>+</mo>
<mn>3</mn>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>12</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中,G(σ11)、G(σ22)、G(σ12)为随机响应分析得到的输送管各应力分量的谱密度;
所述步骤(6)中疲劳损伤计算方法如下:
输送管疲劳损伤DL通过如下公式得到:
<mrow>
<msub>
<mi>D</mi>
<mi>L</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>p</mi>
<mo>&rsqb;</mo>
<mi>T</mi>
</mrow>
<mi>C</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&sigma;</mi>
<mi>b</mi>
</msub>
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>b</mi>
</msub>
<mo>-</mo>
<mi>&sigma;</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>&infin;</mi>
</msubsup>
<msup>
<mi>S</mi>
<mi>m</mi>
</msup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>S</mi>
</mrow>
其中,C、m为材料疲劳曲线参数,σb为管路材料在设计温度下的强度极限,T为输送管振动时间,p(S)为Mises应力概率密度,E[p]为单位之间内的峰值数;
所述单位之间内的峰值数E[p]通过如下公式得到:
<mrow>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>p</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<msub>
<mi>M</mi>
<mn>4</mn>
</msub>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
</mrow>
其中,n=2,4,f1为随机激励的频率下限,f2为随机激励的频率上限;
所述Mises应力概率密度p(S)通过公式:
<mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
<mi>q</mi>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mi>z</mi>
<mi>q</mi>
</mfrac>
</mrow>
</msup>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>D</mi>
<mn>2</mn>
</msub>
<mi>z</mi>
</mrow>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
<msup>
<mi>ze</mi>
<mrow>
<mo>-</mo>
<mfrac>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
</mfrac>
</mrow>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msqrt>
<msub>
<mi>M</mi>
<mn>0</mn>
</msub>
</msqrt>
</mrow>
</mfrac>
</mrow>
其中, D3=1-D1-D2,
所述步骤(7)中判断输送管疲劳损伤是否满足要求方法如下:
若输送管疲劳损伤DL≥0.1,则输送管疲劳损伤不满足要求,反之若DL<0.1满足要求;
所述的σ为静强度分析得到的输送管Mises应力;
所述的σ11、σ22、σ12为随机响应分析得到的输送管各应力分量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310577106.3A CN103678775B (zh) | 2013-11-18 | 2013-11-18 | 一种输送管动强度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310577106.3A CN103678775B (zh) | 2013-11-18 | 2013-11-18 | 一种输送管动强度分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103678775A CN103678775A (zh) | 2014-03-26 |
CN103678775B true CN103678775B (zh) | 2017-08-29 |
Family
ID=50316311
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310577106.3A Active CN103678775B (zh) | 2013-11-18 | 2013-11-18 | 一种输送管动强度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103678775B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105651478A (zh) * | 2015-12-15 | 2016-06-08 | 西安交通大学青岛研究院 | 一种基于振动信号测试零部件疲劳寿命的分析方法 |
WO2018116787A1 (ja) * | 2016-12-20 | 2018-06-28 | Jfeスチール株式会社 | 鋼管の変形性能評価方法、鋼管の製造方法 |
CN116432335B (zh) * | 2023-03-03 | 2024-06-25 | 江阴市南方管件制造有限公司 | 一种基于有限元分析的化工用非标管件的制备工艺 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101509855B (zh) * | 2009-03-17 | 2012-01-25 | 华东理工大学 | 一种高温管道损伤及寿命的在线预测方法 |
CN102567567B (zh) * | 2011-11-15 | 2014-05-28 | 北京宇航系统工程研究所 | 一种基于有限元分析的管路随机振动疲劳寿命分析方法 |
-
2013
- 2013-11-18 CN CN201310577106.3A patent/CN103678775B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103678775A (zh) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104699976B (zh) | 一种包含平均应力影响的金属材料多轴高周疲劳失效预测方法 | |
CN103926152A (zh) | 一种高温多轴谱载下低周蠕变-疲劳寿命评估方法 | |
CN103942418B (zh) | 一种多轴载荷条件下细节疲劳强度额定值的确定方法 | |
CN102567567B (zh) | 一种基于有限元分析的管路随机振动疲劳寿命分析方法 | |
CN107145641B (zh) | 叶片振动疲劳概率寿命预测方法 | |
CN104850713B (zh) | 机械结构随机振动动态应力高精度计算方法 | |
CN107644149B (zh) | 考虑载荷分布及损伤一致性的多通道均衡关联方法 | |
CN104316388B (zh) | 一种对各向异性材料结构件进行疲劳寿命测定的方法 | |
CN103678775B (zh) | 一种输送管动强度分析方法 | |
CN106777457B (zh) | 一种固体发动机药柱结构可靠性评估软件系统 | |
CN102562239A (zh) | 一种航空发动机排气温度的监测方法 | |
CN109726411B (zh) | 风力机机舱结构疲劳强度的计算方法 | |
CN102539315B (zh) | 一种快速准确确定金属塑性成形过程摩擦系数的方法 | |
CN105260574A (zh) | 一种基于临界面法疲劳破坏准则的高周多轴疲劳寿命预测方法 | |
CN109726414A (zh) | 一种风力发电机轮毂结构疲劳计算方法及疲劳分析装置 | |
CN111090957B (zh) | 一种高温结构危险点应力-应变计算方法 | |
CN106934144B (zh) | 汽车悬置支架设计方法 | |
CN113919078A (zh) | 一种尺寸效应下耦合应力梯度的缺口结构疲劳分析方法 | |
CN106446503B (zh) | 遗忘自协方差矩阵递推主元的时变工作模态识别方法 | |
CN106407620B (zh) | 一种基于abaqus的工程结构响应面随机有限元分析处理方法 | |
Gan et al. | Fatigue reliability analysis of turbine disk alloy using saddlepoint approximation | |
CN107944116B (zh) | 一种针对时变结构的瞬态能量响应高效预示方法 | |
CN115221607B (zh) | 疲劳分析方法、终端及存储介质 | |
CN114936494A (zh) | 数据物理融合驱动的高温部件可靠性评定方法和系统 | |
CN109883709B (zh) | 一种基于相对等效应变的随机多轴热机计数方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |