CN114065572A - 一种厚度偏差与叶片厚度分布相关的叶片模型构建方法 - Google Patents
一种厚度偏差与叶片厚度分布相关的叶片模型构建方法 Download PDFInfo
- Publication number
- CN114065572A CN114065572A CN202111253781.1A CN202111253781A CN114065572A CN 114065572 A CN114065572 A CN 114065572A CN 202111253781 A CN202111253781 A CN 202111253781A CN 114065572 A CN114065572 A CN 114065572A
- Authority
- CN
- China
- Prior art keywords
- blade
- thickness
- model
- thickness deviation
- distribution
- 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
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Abstract
本发明公开了一种厚度偏差与叶片厚度分布相关的叶片模型构建方法,通过叶片不同位置控制点处叶片厚度分布确定叶片厚度偏差在叶片不同位置的分布情况,以实现考虑叶片厚度偏差几何模型的重新构建。本发明采用参数化方法重建叶片模型并将代理模型方法用于不确定性分析,既能保证在叶片前、尾缘半径、叶片弦长等参数不变的情况下改变叶片厚度分布,同时保证了叶片型线的连续性。
Description
技术领域
本发明属于航空技术领域,具体涉及一种叶片模型构建方法及一种构造叶片厚度偏差影响叶片气动性能的代理模型方法。
背景技术
压气机叶片是航空发动机中最重要的部件之一,其性能的优劣直接影响着发动机的经济性、气动性能和可靠性。随着航空发动机对压气机提出的性能要求越来越高,压气机叶片的气动热力学设计越来越多的突破了传统技术的约束,采取了越来越多的如弯、掠叶片的新技术。而压气机叶片载荷高、扭曲度大、数量多,其生产量占到了整台发动机的30%-40%,这对加工工艺和日常维护提出了很高的要求。
在叶片的制造加工过程中,因加工振动、工件内部应变等因素很容易造成叶片变形,实际加工得到的叶片与设计叶片厚度分布之间不可避免的会出现一些偏差。虽然叶片各位置处厚度偏差均较小,但其对叶片气动性能的综合影响是不容忽视的。文献(Wu CY.Arbitrary Surface Flank Milling and Flank SAM in the Design andManufacturing of Jet Engine Fan and Compressor Airfoils,ASME Turbo Expo 2012,2012)发现某大型风扇叶片中某些区域0.0254mm-0.0762mm的加工误差足以改变叶片的服役性能和寿命。文献(高丽敏,蔡宇桐,曾瑞慧,等.叶片加工误差度压气机叶栅气动性能的影响,推进技术,2017)表明0.1mm的叶片轮廓度偏差最高可导致20.4%的叶片性能恶化。因此,对叶片厚度偏差影响叶片性能的不确定性研究具有重要的实际意义。
以往对于叶片厚度偏差的研究往往依托于叶片整体轮廓度偏差对叶片气动性能带来的影响,而叶片轮廓度偏差则包含了叶片厚度,弦长和前、尾缘半径偏差,无法得知叶片厚度偏差对叶片性能影响程度的大小。同时,在以往大部分针对叶片厚度偏差研究的叶片模型重建时,往往将叶片吸力面和压力面整体平移,即认为叶片各区域厚度偏差值相同,这对数值研究起到了一定的简化效果,但无法保证叶片前、尾缘附近叶片型线的连续性。
发明内容
为了克服现有技术的不足,本发明提供了一种厚度偏差与叶片厚度分布相关的叶片模型构建方法,通过叶片不同位置控制点处叶片厚度分布确定叶片厚度偏差在叶片不同位置的分布情况,以实现考虑叶片厚度偏差几何模型的重新构建。本发明采用参数化方法重建叶片模型并将代理模型方法用于不确定性分析,既能保证在叶片前、尾缘半径、叶片弦长等参数不变的情况下改变叶片厚度分布,同时保证了叶片型线的连续性。
本发明解决其技术问题所采用的技术方案包括如下步骤:
步骤1:叶片参数化;
根据叶片的几何型线构造参数化控制点,获得不同叶高截面、不同弦长位置处叶片厚度分布;
步骤2:获得厚度偏差分布情况;
根据给定叶片厚度偏差值,以及各控制点处叶片厚度与叶片最大厚度比,计算得到各控制点处厚度偏差值;
步骤3:修改控制点并生成叶片几何模型;
通过厚度叠加法修改各控制点处的厚度分布,获得一组新的叶片型面控制点组合;
步骤4:代理模型初始化;
采用4阶非嵌入式混沌多项式方法,通过叶片厚度偏差概率分布形式以及叶片厚度偏差所满足的公差带生成代理模型;求得相应采样点并通过步骤3中方法生成叶片几何模型;
采用Numeca进行叶片网格划分及性能计算得到叶片性能参数;
通过控制点叶片厚度偏差值与叶片性能参数初次训练代理模型;
步骤5:代理模型验证;
在叶片厚度公差带内随机抽取若干偏差值,并进行求解,将求解结果与代理模型计算的性能进行对比,若精度满足要求则认为代理模型收敛;若精度不够,则提升非嵌入式混沌多项式方法阶数以训练更高精度代理模型,直到满足精度要求;
步骤6:按照NIPC方法中对均值和标准差的计算方法,获得输入概率密度分布时叶片性能的统计参数,并获得对应的最终代理模型。
本发明的有益效果如下:
本发明可以快速的对带有厚度偏差的压气机叶片进行模型重建,所得到的叶片模型既能保证叶片型线连续,同时保证叶片其他几何参数不受影响,从而完成对叶片厚度偏差影响叶片性能的不确定性研究。
附图说明
图1为本发明方法实施流程图。
图2为本发明实施例拟合截面位置示意图。
图3为本发明实施例拟合优化过程。
图4为本发明实施例厚度改变前后沿径向各截面最大厚度分布对比。
图5为本发明实施例厚度改变前后S6截面吸、压力面厚度分布对比。
图6为本发明实施例代理模型的搭建过程。
图7为本发明实施例代理模型精度测试数值模拟结果;其中图(a)是代理模型与数值模拟对叶片峰值效率工况质量流量计算结果对比图;图(b)是代理模型与数值模拟对叶片峰值效率工况绝热效率计算结果对比图;图(c)是代理模型与数值模拟对叶片峰值效率工况总压比计算结果对比图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
一种厚度偏差与叶片厚度分布相关的叶片模型构建方法,包括如下步骤:
步骤1:叶片参数化;
根据叶片的几何型线构造参数化控制点,获得不同叶高截面、不同弦长位置处叶片厚度分布;
步骤2:获得厚度偏差分布情况;
根据给定叶片厚度偏差值,以及各控制点处叶片厚度与叶片最大厚度比,计算得到各控制点处厚度偏差值;
步骤3:修改控制点并生成叶片几何模型;
通过厚度叠加法修改各控制点处的厚度分布,获得一组新的叶片型面控制点组合;
步骤4:代理模型初始化;
采用4阶非嵌入式混沌多项式方法,通过叶片厚度偏差概率分布形式以及叶片厚度偏差所满足的公差带生成代理模型;求得相应采样点并通过步骤3中方法生成叶片几何模型;
采用Numeca进行叶片网格划分及性能计算得到叶片性能参数;
通过控制点叶片厚度偏差值与叶片性能参数初次训练代理模型;
步骤5:代理模型验证;
在叶片厚度公差带内随机抽取若干偏差值,并进行求解,将求解结果与代理模型计算的性能进行对比,若精度满足要求则认为代理模型收敛;若精度不够,则提升非嵌入式混沌多项式方法阶数以训练更高精度代理模型,直到满足精度要求;
步骤6:按照NIPC方法中对均值和标准差的计算方法,获得输入概率密度分布时叶片性能的统计参数,并获得对应的最终代理模型。
具体实施例:
下面结合NASA Rotor 37叶片几何模型重构实例以及CFD计算结果,详细说明本发明的具体实施过程。
1:如图2所示,沿叶片径向均匀取11个截面,从叶根到叶顶分别命名为S1~S11截面。叶型积叠方式采用中弧线积叠,将中弧线40%弦长位置视为积叠点,该积叠方式在改变叶型厚度分布时,截面不会产生周向和轴向偏移。在中弧线的基础上采用厚度叠加法获取截面吸、压力面型线。采用B样条曲线(B-Spline Curve)对叶型各截面进行参数化拟合,选择模拟退火法作为拟合优化算法,拟合优化过程如图3所示。
每个截面的吸、压力面各设置30个控制点。其中,1~26号控制点均匀分布于叶型表面,用于获得截面厚度分布;27和28号控制点用于获得截面弦长参数;29和30号控制点用于获取叶片前、尾缘半径大小。
将Si截面上用于获取吸力面(SS)和压力面(PS)厚度分布的控制点沿前缘到尾缘的方向分别命名为Di,s,1…Di,s,26和Di,p,1…Di,p,26,将各控制点处截面厚度分别记为di,j,k,其中,i=1…11,表示S1~S11截面;j=s,p表示截面吸力面和压力面;k=1…26,表示1~26号控制点。
2:如图4和图5所示所示,通过拟合得到的叶型厚度分布得到叶片最大厚度dmax,并分别找出各截面吸、压力面最大厚度,记为Mdi,j,i和j的含义与上述相同。
对于给定的厚度偏差值E,通过Mdi,j获得各截面吸、压力面最大厚度偏差值Ei,j:
在计算沿各截面吸、压力面型线控制点厚度偏差值时,考虑到尽量保持叶型的连续性,则将Di,j,1和Di,j,26视为吸、压力面与前、尾缘的交点,厚度偏差为0。由于型线上各控制点为均匀分布,则可通过此两点获得型线上各控制点处名义厚度ndi,j,k:
式中,k=1…26,表示1~26号控制点。通过比较型线上26个控制点处名义厚度大小,得到该型线最大名义厚度,记为Mndi,j。则各控制点处厚度偏差量ei,j,k可通过最大厚度偏差值以及名义厚度比计算获得:
修改控制点位置处叶片厚度分布,并结合原始叶片前、尾缘及弦长等参数,生成新的叶片几何模型。
4:通过叶片厚度偏差概率分布形式以及所满足的公差带,确定代理模型需要用到的正交多项式组,首先通过4阶NIPC方法采样点生成方法确定输入样本。
通过第2步和第3步进行样本叶片几何模型重建,调用CFD软件Numeca进行网格划分和流场计算,获得所有样本的性能参数。将样本用于初次训练代理模型,代理模型搭建过程见图6。代理模型初始化完成后,输入厚度偏差值即可输出对应的性能参数。
5:在公差带内随机抽取若干偏差值,并调用CFD程序进行求解,将CFD求解结果与代理模型输出的性能参数值进行对比,若有效精确数字位数足够,则认为代理模型精度足够;若有效精确数字位数不足,则提升NIPC方法阶数,重复步骤4以训练更高精度代理模型。
本发明以叶片厚度偏差满足公差带为[-0.16mm,0.16mm]且概率分布为正态分布为例构建了代理模型,输入量为叶片厚度偏差,输出量为峰值效率工况时的质量流量、绝热效率和总压比。图7将CFD对随机抽样样本的计算结果与代理模型输出结果进行了对比,其中纵坐标为代理模型的输出结果,横坐标为通过CFD数值模拟得到的结果,三张图分别通过压气机的质量流量、绝热效率和增压比对代理模型进行了精度测试,图(a)是代理模型与数值模拟对叶片峰值效率工况质量流量计算结果对比图,图(b)是代理模型与数值模拟对叶片峰值效率工况绝热效率计算结果对比图,图(c)是代理模型与数值模拟对叶片峰值效率工况总压比计算结果对比图。三张图中性能分布对比均近似服从y=x,发现4阶NIPC代理模型得到的质量流量结果已具有至少6位有效数字精度,绝热效率和总压比具有至少5位有效数字精度,认为该模型精度足够。
6:按照NIPC方法中对均值和标准差的计算方法,获得输入概率密度分布时叶片性能的统计参数,并获得对应的代理模型。
Claims (1)
1.一种厚度偏差与叶片厚度分布相关的叶片模型构建方法,其特征在于,包括如下步骤:
步骤1:叶片参数化;
根据叶片的几何型线构造参数化控制点,获得不同叶高截面、不同弦长位置处叶片厚度分布;
步骤2:获得厚度偏差分布情况;
根据给定叶片厚度偏差值,以及各控制点处叶片厚度与叶片最大厚度比,计算得到各控制点处厚度偏差值;
步骤3:修改控制点并生成叶片几何模型;
通过厚度叠加法修改各控制点处的厚度分布,获得一组新的叶片型面控制点组合;
步骤4:代理模型初始化;
采用4阶非嵌入式混沌多项式方法,通过叶片厚度偏差概率分布形式以及叶片厚度偏差所满足的公差带生成代理模型;求得相应采样点并通过步骤3中方法生成叶片几何模型;
采用Numeca进行叶片网格划分及性能计算得到叶片性能参数;
通过控制点叶片厚度偏差值与叶片性能参数初次训练代理模型;
步骤5:代理模型验证;
在叶片厚度公差带内随机抽取若干偏差值,并进行求解,将求解结果与代理模型计算的性能进行对比,若精度满足要求则认为代理模型收敛;若精度不够,则提升非嵌入式混沌多项式方法阶数以训练更高精度代理模型,直到满足精度要求;
步骤6:按照NIPC方法中对均值和标准差的计算方法,获得输入概率密度分布时叶片性能的统计参数,并获得对应的最终代理模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111253781.1A CN114065572B (zh) | 2021-10-27 | 2021-10-27 | 一种厚度偏差与叶片厚度分布相关的叶片模型构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111253781.1A CN114065572B (zh) | 2021-10-27 | 2021-10-27 | 一种厚度偏差与叶片厚度分布相关的叶片模型构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114065572A true CN114065572A (zh) | 2022-02-18 |
CN114065572B CN114065572B (zh) | 2023-04-07 |
Family
ID=80235904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111253781.1A Active CN114065572B (zh) | 2021-10-27 | 2021-10-27 | 一种厚度偏差与叶片厚度分布相关的叶片模型构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114065572B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116160291A (zh) * | 2023-04-24 | 2023-05-26 | 沈阳强航时代精密科技有限公司 | 一种基于在机测量的叶片型面加工刀长补偿方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3009104A1 (fr) * | 2013-07-29 | 2015-01-30 | Snecma | Procede de modelisation d'une pale d'une helice non-carenee |
US20150241211A1 (en) * | 2014-02-26 | 2015-08-27 | Hans Joachim Pehle | Method and program for analyzing pipe wall-thickness variations |
WO2016101182A1 (zh) * | 2014-12-23 | 2016-06-30 | 清华大学 | 基于贝叶斯网络和极限学习机的区间型指标预报方法 |
CN108090275A (zh) * | 2017-12-14 | 2018-05-29 | 中国航发沈阳发动机研究所 | 一种面向参数化模型的涡轮气冷叶片构型方法 |
CN109409019A (zh) * | 2018-12-15 | 2019-03-01 | 江西理工大学南昌校区 | 一种风力机叶片气动外形的预弯优化方法 |
CN110555214A (zh) * | 2018-05-30 | 2019-12-10 | 中国航发商用航空发动机有限责任公司 | 压气机叶片叶型构建方法及压气机叶片 |
CN110727995A (zh) * | 2018-06-29 | 2020-01-24 | 中国航发商用航空发动机有限责任公司 | 叶片形状的构造方法、叶片以及计算机设备 |
CN112347579A (zh) * | 2020-11-04 | 2021-02-09 | 中国科学院工程热物理研究所 | 一种压气机叶型设计方法及压气机叶型 |
-
2021
- 2021-10-27 CN CN202111253781.1A patent/CN114065572B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3009104A1 (fr) * | 2013-07-29 | 2015-01-30 | Snecma | Procede de modelisation d'une pale d'une helice non-carenee |
US20150241211A1 (en) * | 2014-02-26 | 2015-08-27 | Hans Joachim Pehle | Method and program for analyzing pipe wall-thickness variations |
WO2016101182A1 (zh) * | 2014-12-23 | 2016-06-30 | 清华大学 | 基于贝叶斯网络和极限学习机的区间型指标预报方法 |
CN108090275A (zh) * | 2017-12-14 | 2018-05-29 | 中国航发沈阳发动机研究所 | 一种面向参数化模型的涡轮气冷叶片构型方法 |
CN110555214A (zh) * | 2018-05-30 | 2019-12-10 | 中国航发商用航空发动机有限责任公司 | 压气机叶片叶型构建方法及压气机叶片 |
CN110727995A (zh) * | 2018-06-29 | 2020-01-24 | 中国航发商用航空发动机有限责任公司 | 叶片形状的构造方法、叶片以及计算机设备 |
CN109409019A (zh) * | 2018-12-15 | 2019-03-01 | 江西理工大学南昌校区 | 一种风力机叶片气动外形的预弯优化方法 |
CN112347579A (zh) * | 2020-11-04 | 2021-02-09 | 中国科学院工程热物理研究所 | 一种压气机叶型设计方法及压气机叶型 |
Non-Patent Citations (5)
Title |
---|
TIANYUAN JI等: "Mechanism Analysis of the Influence of Blade Thickness Deviation on The Performance of Axial Flow Compressor", 《TURBO EXPO: POWER FOR LAND, SEA, AND AIR》 * |
余锐等: "基于自适应Kriging代理模型的叶型气动优化设计", 《动力工程学报》 * |
刘润泽等: "非均匀有理B样条曲线及节点插入算法在透平叶片优化设计中的应用", 《航空动力学报》 * |
李丹等: "GFRP风机叶片结构设计的二级优化方法", 《南京航空航天大学学报》 * |
李超东等: "加工误差对2维叶型损失系数影响的数值分析", 《航空发动机》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116160291A (zh) * | 2023-04-24 | 2023-05-26 | 沈阳强航时代精密科技有限公司 | 一种基于在机测量的叶片型面加工刀长补偿方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114065572B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lange et al. | Introduction of a parameter based compressor blade model for considering measured geometry uncertainties in numerical simulation | |
CN110929357A (zh) | 一种高性能舰船燃机压气机气动设计方法 | |
CN112685968B (zh) | 基于空间负荷定制思想的轴流压气机气动设计方法 | |
CN112765732B (zh) | 一种基于选区激光熔化工艺的航空叶片拓扑优化设计方法 | |
CN114065572B (zh) | 一种厚度偏差与叶片厚度分布相关的叶片模型构建方法 | |
CN112287580A (zh) | 一种基于全三维数值模拟的轴流压气机喘振边界计算方法 | |
CN115081130B (zh) | 基于动态支持向量回归的叶片与端壁联合气动优化方法 | |
CN112594064A (zh) | 一种基于轴流压气机级间测量参数的s2流场诊断方法 | |
CN113283024A (zh) | 一种基于有限元建模的拓扑优化和参数优化的风机轮毂结构优化方法 | |
CN111027148B (zh) | 损失落后角模型自动标定及工业轴流压缩机性能计算方法 | |
CN112270139A (zh) | 基于母型库的燃料电池离心压缩机气动优化设计方法 | |
CN114429090A (zh) | 一种基于数据驱动的压气机叶片鲁棒性设计方法 | |
CN109595040A (zh) | 一种燃气轮机叶片扭曲规律设计方法 | |
CN112487570A (zh) | 基于自由变形技术的离心压缩机分流叶片形状优化方法 | |
EP4258074A1 (en) | Systems and methods for aerodynamic analysis for inspected bladed rotors | |
CN115859536B (zh) | 压气机转子叶片非同步振动锁频数值模拟方法 | |
CN115994991A (zh) | 一种涡轮叶片模型冷热态快速转换方法 | |
CN114861315A (zh) | 一种基于机器学习的叶轮机二维叶型优化方法 | |
CN115680784B (zh) | 层合结构树脂基复合材料叶片建模及计算方法 | |
CN113959691B (zh) | 一种中介机匣气动性能试验件的设计方法 | |
CN113742861B (zh) | 一种适用于风力发电机风洞试验的叶片模型优化设计方法 | |
CN114510798B (zh) | 一种航空发动机气动轴向力防错分析方法 | |
CN115270343A (zh) | 轴流涡轮一维反问题高可靠度气动设计方法 | |
CN114547786A (zh) | 一种考虑加工变形的弯扭薄壁叶片测点规划方法 | |
Li et al. | The optimization of a centrifugal impeller based on a new multi-objective evolutionary strategy |
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 |