CN111241653A - 基于蒙特卡洛法的多维度叶片误差性能影响仿真实现方法 - Google Patents
基于蒙特卡洛法的多维度叶片误差性能影响仿真实现方法 Download PDFInfo
- Publication number
- CN111241653A CN111241653A CN201910318468.8A CN201910318468A CN111241653A CN 111241653 A CN111241653 A CN 111241653A CN 201910318468 A CN201910318468 A CN 201910318468A CN 111241653 A CN111241653 A CN 111241653A
- Authority
- CN
- China
- Prior art keywords
- geometric
- error
- blade
- grid
- module
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000004088 simulation Methods 0.000 title claims abstract description 31
- 238000000342 Monte Carlo simulation Methods 0.000 title claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 20
- 238000005070 sampling Methods 0.000 claims abstract description 20
- 238000009826 distribution Methods 0.000 claims abstract description 14
- 238000013515 script Methods 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 12
- 239000012530 fluid Substances 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000005315 distribution function Methods 0.000 claims description 2
- 230000003068 static effect Effects 0.000 claims description 2
- 238000010923 batch production Methods 0.000 claims 2
- 238000005094 computer simulation Methods 0.000 abstract description 2
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Classifications
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
一种基于蒙特卡洛法的多维度叶片几何误差性能影响仿真实现方法,针对几何误差类型设置高斯分布概率密度函数,然后在选定的叶片控制截面上抽样加入几何误差以批量生成新叶片几何样本,用于自动生成网格后通过仿真计算获得压气机性能参数及其概率估计量。本发明通过批量生成几何样本、网格自动生成、计算仿真及计算结果自动获取实现了多维度叶片几何误差性能影响的自动化仿真流程,以减少仿真所耗费的人力及时间成本。
Description
技术领域
本发明涉及的是一种航空发动机制造领域的技术,具体是一种基于蒙特卡洛法的多维度叶片几何误差性能影响仿真实现方法,可应用于航空发动机等各种叶轮机械领域。
背景技术
航空发动机中的轴流压气机叶片作为高负荷做功部件,其截面位置度、扭转度及轮廓度等参数对于压气机性能有不可忽视的影响。由于实际加工中误差种类及其分布位置具有高度随机性,为得到叶片几何误差与压气机性能间的关系,所需要的叶片几何样本数量巨大,实验成本很高。即使应用计算流体力学对每个样本进行逐一几何生成和计算仿真,也需要大量的人力成本与工作时间。
发明内容
本发明针对现有叶片几何误差的性能影响研究多数仍采用单维度的误差叶片仿真方法无法适用于多维度叶片几何误差的性能影响分析的不足,提出一种基于蒙特卡洛法的多维度叶片几何误差性能影响仿真实现方法,通过批量生成几何样本、网格自动生成、计算仿真及计算结果自动获取实现了多维度叶片几何误差性能影响的自动化仿真流程,以减少仿真所耗费的人力及时间成本。
本发明是通过以下技术方案实现的:
本发明针对几何误差类型设置高斯分布概率密度函数,然后在选定的叶片控制截面上抽样加入几何误差以批量生成新叶片几何样本,用于自动生成网格后通过仿真计算获得压气机性能参数及其概率估计量。
所述的几何误差类型包括:叶片控制截面叶型的轴向位置度、周向位置度、扭转度及轮廓度。
所述的高斯分布概率密度函数为:其中:μ为几何误差的平均值,取值为0,σ为标准差,取值由误差类型对应的公差范围而定;使几何偏差的分布满足:偏离原型越远,其概率密度越低,且到达公差边界时概率密度接近于0的状态。
所述的叶片控制截面是指:垂直于径向的叶片横截面,其数量与位置可依需要确定。
所述的抽样加入是指:对单类几何误差的样本值采用但不限于拉丁超立方抽样方法经随机组合后得到带多维度误差的叶片几何样本。
所述的网格,采用但不限于商业计算流体力学软件(Numeca)自动生成,其对应的网格绘制模块为Autogrid,对应的计算模块为Fine Turbo。
所述的压气机性能参数包括但不限于一级压气机的折合流量、总压比、等熵效率、轴向推力及扭矩;所述的概率估计量包括但不限于期望值、均方差及拟合概率密度分布函数。
本发明涉及一种实现上述方法的系统,包括:样本抽样模块、几何生成模块、网格生成模块、性能仿真模块及结果读取分析模块,其中:样本抽样模块与几何生成模块相连并传输几何误差样本值用于生成误差模型,几何生成模块与网格生成模块相连并传输误差模型用于生成相应网格,网格生成模块与性能仿真模块相连并传输网格数据用以计算压气机级性能参数,性能仿真模块与结果读取分析模块相连并传输性能参数计算结果,最后得到结果报告及性能参数概率估计量。
技术效果
本发明实现了叶片多维度几何误差性能影响仿真,与现有的单维度叶片误差构造相比,此方法构造的叶片几何更贴近实际生产情况,仿真结果更具工程应用价值。同时本发明通过自动化流程设计,减少了仿真消耗的人力成本及时间成本,提高了仿真效率,使得更大样本的批量计算成为可能。
附图说明
图1为实施例的压气机转、静子示意图;
图2为实施例的压气机转子叶片轴向位置度产生偏差时的几何示意图;
图中:ΔZ表示轴向误差,数字为轴向误差大小示例;
图3为实施例的压气机转子叶片周向位置度产生偏差时的几何示意图;
图中:ΔY表示周向误差,数字为周向误差大小示例;
图4为实施例的压气机转子叶片扭转度定义及其产生偏差时的几何示意图;
图中:γ表示扭转度,C表示弦长,Δγ表示扭转度误差,数字为周向误差大小示例;
图5为实施例的压气机转子叶片轮廓度产生偏差时的几何示意图;
图中:ΔP表示周向误差,数字为周向误差大小示例;
图6为实施例计算生成的概率密度分布拟合图;
图中:横轴为新叶片性能参数与原型参数的相对误差,纵轴为概率密度,直方图为样本概率密度分布情况,其面积等于样本出现频率,曲线为概率密度分布正态拟合结果。
具体实施方式
如图1所示,本实施例为某轴流压气机的转子叶片几何误差对压气机级性能的影响分析,具体包括以下步骤:
步骤1、叶片几何误差的概率密度函数设置及多截面组合抽样:
1.1如图2至图5所示,所述的叶片几何误差为:轴向位置度误差ΔZ、周向位置度误差ΔY、扭转度误差Δγ及轮廓度ΔP。
1.2所述的概率密度函数为:对于本压气机叶片,确定轴、周向位置度、扭转度及轮廓度的高斯分布概率密度函数标准差σ分别为0.5mm、0.5°及0.05mm。
1.3所述的多截面组合抽样为:将动叶沿径向15%、50%、85%叶高处横截面作为3个控制截面,分别进行上述4类几何误差抽样,使得叶片共具有12个维度的几何误差。
1.4所述的抽样所用的方法具体为:确定叶片样本数量,此例中样本数量为200片,随后进行拉丁超立方抽样,将每个维度的误差概率密度函数分为等概率的200个误差区间,从每个区间中随机抽取一个误差样本值,再进行多维度误差样本值的随机组合,最终得到200个叶片的几何误差样本,记录为表格形式误差文件。
步骤2、基于原型几何模型数据及几何误差抽样数据,在径向控制截面内加入几何误差,批量生成新叶片的几何模型数据:
2.1所述的几何模型数据为叶片三个控制截面的轮廓控制点坐标集合,本例中每个截面各有190个轮廓控制点。
2.2所述的叶片径向截面内的坐标系是以叶轮机械转轴为中心轴的圆柱坐标系。
2.3所述的批量生成的方法具体为:通过python脚本读取误差文件内的叶片几何误差值,在原型叶片控制截面轮廓点的基础上进行平移、旋转及缩放坐标变换,生成具有位置度、扭转度、轮廓度误差的截面轮廓点坐标,再将三个截面组合形成新叶片的几何模型。
此例中的200个算例使用英特尔酷睿i7-5500U CPU处理器生成,耗时7秒。
步骤3、执行脚本,将新叶片几何模型导入商业计算流体力学软件Numeca中的网格模块Autogrid,使用与原型网格相同的参数自动批量生成网格并保存。
3.1所述的脚本文件为Autogrid性能影响仿真模块支持的Python语言脚本,用以自动设置并生成网格数据文件。
3.2所述的原型网格的参数具体为:转子流道径向网格点数为105,静子流道径向网格点数为89,转子和静子叶片径向网格点数均为33。
3.3所述的批量生成的方法为:将多个新叶片几何对应的网格生成指令写入一个代码文件中并运行,即可批量获得所有网格文件。
此例中网格生成使用英特尔酷睿i7-5500U CPU处理器生成,每个样本耗时约2.5分钟,200个样本网格生成共耗时约8.5小时。
步骤4、执行脚本,将网格文件导入Numeca中的计算模块Fine Turbo中,使用与原型计算文件相同的设置和参数进行批量计算。
4.1所述的脚本文件为Fine Turbo模块支持的Python语言脚本,用以自动执行计算步骤。
4.2所述的原型计算文件的设置和参数具体为:流体为空气,流动模型设为定常,湍流模型使用Spalart-Allmaras;给定转子旋转速度-1.6×104rpm、进口给定总温600K、总压9.5×105Pa及流动与轴向夹角arctg(Vt/Vz)=-0.1rad,arctg(Vr/Vz)=-0.1rad;出口给定静压1.08×103kPa,满足径向平衡方程。
4.3所述的批量生成的方法为:将多个文件的设置计算指令写入一个代码文件中运行,即可批量进行运算。
此例中的每个样本在英特尔CPU处理器上进行本地四核并行计算,每个样本平均耗时约1小时,且该CPU可支持四个样本并行计算,本例中200个样本共耗时约50小时。
步骤5、如图6所示,自动获取计算结果并得到概率估计量。
5.1所述的自动获取计算结果的方法为:运行Python语言脚本批量读取计算结果数据,本例中读取压气机级进出口质量流量、总压比、等熵效率、扭矩及轴向推力数据,记录为结果文件。
5.2所述的得到概率估计量的方法为:运行Python语言脚本从结果文件中读取结果数值大小,并导入MATLAB运行概率统计分析脚本计算得到概率估计量的值,本例中对计算所得样本无量纲化后的总压比参数进行概率估计量计算及正态分布拟合,得到参数数学期望值为1.0000,标准差为1.081×10-3,拟合得到正态分布数学期望值为1.0000,标准差为1.091×10-3。
本方法的几何生成模块对系统最终技术效果给出了最大贡献,显著提高了叶片几何误差维度并使模型生成耗时显著降低。在本实施例中,以叶片几何误差维度和仿真中各步骤耗时为性能指标进行本方法的评估,对比现有的单维度叶片几何误差性能仿真方法与本方法的性能指标,如表1所示。由表1可见,相比于现有方法,叶片几何误差参数由1个提升至12个,各仿真步骤耗时均有降低,本实施例中整个系统可于3日内自动完成200个样本的生成与计算,有效降低了时间成本及人力成本。
表1各性能指标对比
本方法实现了叶片多维度几何误差性能影响仿真,与现有的单维度叶片误差构造相比,此方法构造的叶片几何更贴近实际生产情况,仿真结果更具工程应用价值。同时本方法通过自动化流程设计,减少了仿真消耗的人力成本及时间成本,提高了仿真效率,使得更大样本的批量计算成为可能。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。
Claims (10)
1.一种基于蒙特卡洛法的多维度叶片几何误差性能影响仿真实现方法,其特征在于,针对几何误差类型设置高斯分布概率密度函数,然后在选定的叶片控制截面上抽样加入几何误差以批量生成新叶片几何样本,用于自动生成网格后通过仿真计算获得压气机性能参数及其概率估计量;
所述的几何误差类型包括:叶片控制截面叶型的轴向位置度、周向位置度、扭转度及轮廓度;
所述的抽样加入是指:对单类几何误差的样本值采用但不限于拉丁超立方抽样方法经随机组合后得到带多维度误差的叶片几何样本。
3.根据权利要求1所述的方法,其特征是,所述的叶片控制截面是指:垂直于径向的叶片横截面。
4.根据权利要求1所述的方法,其特征是,所述的网格,采用Numeca自动生成,其对应的网格绘制模块为Autogrid,对应的计算模块为Fine Turbo。
5.根据权利要求1所述的方法,其特征是,所述的压气机性能参数包括:一级压气机的折合流量、总压比、等熵效率、轴向推力及扭矩;所述的概率估计量包括:期望值、均方差及拟合概率密度分布函数。
6.根据权利要求1所述的方法,其特征是,具体包括:
步骤1、叶片几何误差的概率密度函数设置及多截面组合抽样;
步骤2、基于原型几何模型数据及几何误差抽样数据,在径向控制截面内加入几何误差,批量生成新叶片的几何模型数据;
步骤3、执行脚本,将新叶片几何模型导入商业计算流体力学软件Numeca中的网格模块Autogrid,使用与原型网格相同的参数自动批量生成网格并保存;
步骤4、执行脚本,将网格文件导入Numeca中的计算模块Fine Turbo中,使用与原型计算文件相同的设置和参数进行批量计算;
步骤5、自动获取计算结果并得到概率估计量。
7.根据权利要求6所述的方法,其特征是,所述的步骤1包括:
1.1设置叶片几何误差为:轴向位置度误差ΔZ、周向位置度误差ΔY、扭转度误差Δγ及轮廓度ΔP;
1.2设置概率密度函数为:对于本压气机叶片,确定轴、周向位置度、扭转度及轮廓度的高斯分布概率密度函数标准差σ分别为0.5mm、0.5°及0.05mm;
1.3设置多截面组合抽样为:将动叶沿径向15%、50%、85%叶高处横截面作为3个控制截面,分别进行上述4类几何误差抽样,使得叶片共具有12个维度的几何误差;
1.4抽样:确定叶片样本数量,此例中样本数量为200片,随后进行拉丁超立方抽样,将每个维度的误差概率密度函数分为等概率的200个误差区间,从每个区间中随机抽取一个误差样本值,再进行多维度误差样本值的随机组合,最终得到200个叶片的几何误差样本,记录为表格形式的误差文件。
8.根据权利要求6所述的方法,其特征是,所述的步骤3包括:
3.1设置脚本文件为Autogrid性能影响仿真模块支持的Python语言脚本,用以自动设置并生成网格数据文件;
3.2设置原型网格的参数具体为:转子流道径向网格点数为105,静子流道径向网格点数为89,转子和静子叶片径向网格点数均为33;
3.3批量生成:将多个新叶片几何对应的网格生成指令写入一个代码文件中并运行,即可批量获得所有网格文件。
9.根据权利要求6所述的方法,其特征是,所述的步骤4包括:
4.1设置脚本文件为Fine Turbo模块支持的Python语言脚本,用以自动执行计算步骤;
4.2设置原型计算文件的设置和参数具体为:流体为空气,流动模型设为定常,湍流模型使用Spalart-Allmaras;转子旋转速度-1.6×104rpm、进口给定总温600K、总压9.5×105Pa及流动与轴向夹角arctg(Vt/Vz)=-0.1rad,arctg(Vr/Vz)=-0.1rad;出口给定静压1.08×103kPa,满足径向平衡方程;
4.3批量生成:将多个文件的设置计算指令写入一个代码文件中运行,即可批量进行运算。
10.一种实现上述任一权利要求所述方法的系统,其特征在于,包括:样本抽样模块、几何生成模块、网格生成模块、性能仿真模块及结果读取分析模块,其中:样本抽样模块与几何生成模块相连并传输几何误差样本值用于生成误差模型,几何生成模块与网格生成模块相连并传输误差模型用于生成相应网格,网格生成模块与性能仿真模块相连并传输网格数据用以计算压气机级性能参数,性能仿真模块与结果读取分析模块相连并传输性能参数计算结果,最后得到结果报告及性能参数概率估计量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910318468.8A CN111241653B (zh) | 2019-04-19 | 2019-04-19 | 基于蒙特卡洛法的多维度叶片误差性能影响仿真实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910318468.8A CN111241653B (zh) | 2019-04-19 | 2019-04-19 | 基于蒙特卡洛法的多维度叶片误差性能影响仿真实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111241653A true CN111241653A (zh) | 2020-06-05 |
CN111241653B CN111241653B (zh) | 2024-03-19 |
Family
ID=70872616
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910318468.8A Active CN111241653B (zh) | 2019-04-19 | 2019-04-19 | 基于蒙特卡洛法的多维度叶片误差性能影响仿真实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111241653B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114091278A (zh) * | 2021-11-26 | 2022-02-25 | 清华大学无锡应用技术研究院 | 一种评估压气机叶片局部制造误差影响的仿真实现方法 |
CN114152552A (zh) * | 2021-11-26 | 2022-03-08 | 国家能源投资集团有限责任公司 | 煤矿地下水库储水系数测算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5462105A (en) * | 1992-08-07 | 1995-10-31 | Supernak; Janusz | Adjustments for window shades |
CN102519407A (zh) * | 2011-12-05 | 2012-06-27 | 西北工业大学 | 一种建立叶片三维公差模型的方法 |
WO2016026355A1 (zh) * | 2014-08-18 | 2016-02-25 | 国家电网公司 | 一种有源配电网电压跌落仿真与评估方法 |
CN108388747A (zh) * | 2018-03-12 | 2018-08-10 | 上海交通大学 | 安装角误差叶片的多通道周向类正弦分布样本实现方法 |
-
2019
- 2019-04-19 CN CN201910318468.8A patent/CN111241653B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5462105A (en) * | 1992-08-07 | 1995-10-31 | Supernak; Janusz | Adjustments for window shades |
CN102519407A (zh) * | 2011-12-05 | 2012-06-27 | 西北工业大学 | 一种建立叶片三维公差模型的方法 |
WO2016026355A1 (zh) * | 2014-08-18 | 2016-02-25 | 国家电网公司 | 一种有源配电网电压跌落仿真与评估方法 |
CN108388747A (zh) * | 2018-03-12 | 2018-08-10 | 上海交通大学 | 安装角误差叶片的多通道周向类正弦分布样本实现方法 |
Non-Patent Citations (2)
Title |
---|
蔡凯武;: "飞机叶片双重螺旋法加工的调整参数优化", 机械与电子, no. 08, 24 August 2018 (2018-08-24) * |
郑似玉;滕金芳;羌晓青;: "位置度超差对轴流压气机流场性能影响的数值研究", 流体机械, no. 11, 30 November 2016 (2016-11-30) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114091278A (zh) * | 2021-11-26 | 2022-02-25 | 清华大学无锡应用技术研究院 | 一种评估压气机叶片局部制造误差影响的仿真实现方法 |
CN114152552A (zh) * | 2021-11-26 | 2022-03-08 | 国家能源投资集团有限责任公司 | 煤矿地下水库储水系数测算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111241653B (zh) | 2024-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Righi et al. | Three-dimensional through-flow modelling of axial flow compressor rotating stall and surge | |
US20210209264A1 (en) | Modeling and calculation aerodynamic performances of multi-stage transonic axial compressors | |
CN109598081B (zh) | 基于数据降维及多二维流面的径流式透平气动优化方法 | |
CN109815624B (zh) | 一种考虑进气总压畸变影响的压气机稳定边界判断方法 | |
CN111241653B (zh) | 基于蒙特卡洛法的多维度叶片误差性能影响仿真实现方法 | |
US10168306B2 (en) | Time reference derivation from time of arrival measurements | |
Popov et al. | Multicriteria optimization of axial low pressure compressor of gas turbine power plant | |
Winhart et al. | Application of the time transformation method for a detailed analysis of multistage blade row interactions in a shrouded turbine | |
CN112268011B (zh) | 多级轴流压气机气动性能预估方法和装置 | |
CN112182802B (zh) | 计入随机失谐的叶轮机械气动弹性优化设计实现方法 | |
Guo et al. | Statistical evaluation of stability margin of a multi-stage compressor with geometric variability using adaptive polynomial chaos-Kriging model | |
CN116663202A (zh) | 一种多级轴流压气机性能仿真模型的校核方法 | |
CN116595874A (zh) | 叶轮机械性能预测模型参数优化方法及装置、存储介质 | |
Fei et al. | Application of new empirical models based on mathematical statistics in the through-flow analysis | |
Simpson et al. | Numerical and experimental study of the performance effects of varying vaneless space and vane solidity in radial inflow turbine stators | |
Kuzmenko et al. | Optimization of the gas turbine engine parts using methods of numerical simulation | |
CN114662218A (zh) | 基于体积力模型的叶轮机模拟、失稳预测方法和装置 | |
Lecheler et al. | Experimental and numerical investigation of the flow in a 5-stage transonic compressor rig | |
CN114091278A (zh) | 一种评估压气机叶片局部制造误差影响的仿真实现方法 | |
Kunte et al. | Experimental investigation of a truncated pipe diffuser with a tandem deswirler in a centrifugal compressor stage | |
CN113673060B (zh) | 多级压气机模化方法及系统 | |
Schemmann et al. | A Multi-Fidelity Sampling Method for Efficient Design and Optimization of Centrifugal Compressor Impellers | |
Righi et al. | Body-force and mean-line models for the generation of axial compressor sub-idle characteristics | |
CN114576177B (zh) | 一种径流式叶轮轮背间隙摩擦损失高精度预测方法 | |
Ziegler et al. | Development of a Novel Axial Compressor Generation for Industrial Applications: Part 1—Compressor Design and Performance |
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 |