CN103486996B - 未知cad模型的航空发动机叶片型面测量方法 - Google Patents
未知cad模型的航空发动机叶片型面测量方法 Download PDFInfo
- Publication number
- CN103486996B CN103486996B CN201310354360.7A CN201310354360A CN103486996B CN 103486996 B CN103486996 B CN 103486996B CN 201310354360 A CN201310354360 A CN 201310354360A CN 103486996 B CN103486996 B CN 103486996B
- Authority
- CN
- China
- Prior art keywords
- point
- blade
- curvature
- curve
- measurement
- 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.)
- Expired - Fee Related
Links
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种未知CAD模型的航空发动机叶片型面测量方法,用于解决现有航空发动机叶片型面三坐标测量方法精度低的技术问题。技术方案是基于对若干条叶片型面上反应主轴方向曲率分布状态的曲线测量与分析,基于曲率及容差控制求出曲线的控制点,进而确定叶片测量截面;基于曲率连续预测法,规划每条测量截面的测量点集,实现对CAD模型未知的航空发动机叶片型面的三坐标测量。由于该方法根据叶片型面曲率分布特点,选取叶片测量截面;并对测量截面基于连续曲率预测法实现测点规划,连续曲率预测法具有比较好的鲁棒性,既减少了型线冗余点的提取,又保证了提取点较好的贴合测量截面型线曲率变化情况。既提高了测量效率又可保证了测量精度。
Description
技术领域
本发明涉及一种航空发动机叶片型面三坐标测量方法,特别涉及一种未知CAD模型的航空发动机叶片型面测量方法。
背景技术
在航空发动机领域,叶片是透平机械的最核心部件,叶片型面的气动外形将直接影响航空发动机的性能。为了提高发动机推重比,叶片的型面通常设计成扭曲的变截面曲线,而精确的叶片型面数据,是叶片实现CAD/CAM、逆向制造数字化、产品检测及质量控制等的关键和前提。为了得到叶片型面数据,需要对叶片型面进行测量,叶片型面的测量总体可规划为两种模式,第一种模式叫检测误差模式,在该模式中被测曲面有设计数据或CAD模型,测量该曲面的目的是为了得到被测曲面与理论数据或数字模型的误差,判断该曲面是否合格。第二种模式叫重构模式,在该模式中,被测曲面既没有设计数据也没有CAD模型,对之进行测量的目的是为了重构该曲面,进而实现对该曲面气动性能的分析与评价,最后实现对该曲面的逆向制造。
目前叶片类零件常用测量方法主要有光学测量、CT断层扫描、三坐标测量,其中三坐标测量是接触式测量仪器的典型代表,具有精度高、通用性强、效率高、测量稳定的优点而被广泛使用;叶片型面属于自由曲面,测量时一般按照“曲面-曲线-点集-测点集”的方式;在叶片测量的重构模式上,现有方法一般根据叶片主轴尺寸,基于等距或分段等距方法规划出叶片若干测量截面,并对每个测量截面基于等步长或分段等步长方法规划出测量点,由于叶片型面曲率分布不均,故现有方法存在明显不足。
一般而言,测量截面数越多、每个测量截面的测量点越密,测量精度就越高;然而过多的测量截面和测量点将严重影响测量效率;不但加重测量工作量,而且难以快速高效的得到叶片型面的测量数据。因此,叶片类自由曲面零件型面的测量应该根据叶片曲率分布特点,在保证整体测量结果精度的基础上,尽可能减少测量截面数量和每条测量截面的测点数。
发明内容
为了克服现有航空发动机叶片型面三坐标测量方法精度低的不足,本发明提供一种未知CAD模型的航空发动机叶片型面测量方法。该方法基于对若干条叶片型面上反应主轴方向曲率分布状态的曲线测量与分析,基于曲率及容差控制求出曲线的控制点,进而确定叶片测量截面;并基于曲率连续预测法,规划每条测量截面的测量点集,实现对CAD模型未知的航空发动机叶片型面精确、高效的三坐标测量。由于该方法根据叶片型面曲率分布特点,选取叶片测量截面;并对测量截面基于连续曲率预测法实现测点规划,连续曲率预测法具有比较好的鲁棒性,既可减少型线冗余点的提取,又能保证提取点较好的贴合测量截面型线曲率变化情况。该方法可避免叶片型面传统测量中根据经验规划测量截面及截面测量点的缺陷,既可提高测量效率又可保证测量精度。
本发明解决其技术问题所采用的技术方案:一种未知CAD模型的航空发动机叶片型面测量方法,其特点是包括以下步骤:
步骤一、沿叶片主轴方向在叶盆、叶背型面划出若干条曲线,曲线在叶盆、叶背曲面呈等间距分布,用三坐标测量机采集各规划曲线上若干点,得到测点集信息。
步骤二、将步骤一得到的测点集信息导入UG,并拟合曲线,将叶盆、叶背曲面上对应的两条曲线拟合曲面,并提取该曲面中分线。
步骤三、基于曲率分布及容差控制的中分线最小控制点集提取。
第一步,任取一条中分线,以曲线的任一端点为起点,向前搜索一微小弧长,计算搜索始末点连线的中点到该段搜索曲线的距离d,设h为曲线弦公差,若|d-h|<ε,h=nc,n为大于0的自然数,c为曲线弦公差单位变化值,考虑计算精度和效率,c设为0.005;ε为设定阈值,则该点为控制点;若h-d>ε,则继续向前搜索。搜索完毕,将搜索到的点拟合曲线,记原始中分线容差值为S;将拟合曲线进行一定数量的点离散,并计算出离散点到原始曲线的最大距离dmax,若|S-dmax|>ε,则调整n,重复上述操作,直到|S-dmax|<ε。
第二步,取其中一条中分线为参考线,顺次过其上控制点作叶片主轴垂面与其余等参数曲线交于一点,计算其余曲线上离交点最近的控制点,如果该控制点与交点Z向坐标差值<nr,n为自然数,r为测头半径,ΔZmin为第一小步所求每条曲线各控制点Z向坐标差值的最小值,则将该控制点以交点代替;如果两点Z向坐标差值>nr,则将交点加入该曲线的控制点集;第一条中分线上控制点遍历结束,再以第二条中分线为参考线重复上述操作,直到所有中分线遍历结束,初步实现全体中分线的控制点集统一。
第三步,以得到的控制点集拟合B样条曲线,设定原始中分线容差值为S;将拟合曲线进行点离散,并计算出离散点到原始中分线的最大距离dmax,若|S-dmax|>ε,则调整第二小步中n,直到|S-dmax|<ε。
步骤四、基于连续曲率预测法提取测量截面测点集。
根据自由曲面多阶连续的特性,利用曲率连续预测法实现测点的自适应规划,建立拟合多项式Z=AX5+BX4+CX3+DX2+EX+F。
利用三坐标测量机在求出的测量截面上初测6个点,并指定搜索方向。
利用前6点的坐标信息拟合5次多项式。
求解第6点一阶、二阶导数和曲率半径;
根据第6点切线方向和曲率半径,预测第7点坐标值。
判断预测点坐标是否回到起始点,未到达则根据第7点预测值引导三坐标测量机测量并获取第7点真实坐标,若到达或超出起始点则搜索结束。
步骤五、将叶片型面测点集数据导入UG,根据各测量截面测点信息拟合出测量截面曲线,并由各测量截面曲线拟合出叶身型面。
步骤六、在重构模型上提取曲率代表线,并与第一步中测量点拟合的曲率代表线进行比较,将重构模型曲率代表线以等弧长方式离散为若干点,并计算各离散点到对应测量点拟合曲率代表线的距离,若大于给定容差值,则增加步骤一中规划的曲率代表线条数,重复上述操作,直到重构模型曲率代表线各离散点到原始曲率代表线的距离小于给定容差值。
本发明的有益效果是:该方法基于对若干条叶片型面上反应主轴方向曲率分布状态的曲线测量与分析,基于曲率及容差控制求出曲线的控制点,进而确定叶片测量截面;并基于曲率连续预测法,规划每条测量截面的测量点集,实现对CAD模型未知的航空发动机叶片型面精确、高效的三坐标测量。由于该方法根据叶片型面曲率分布特点,选取叶片测量截面;并对测量截面基于连续曲率预测法实现测点规划,连续曲率预测法具有比较好的鲁棒性,既减少了型线冗余点的提取,又保证了提取点较好的贴合测量截面型线曲率变化情况。该方法避免了叶片型面传统测量中根据经验规划测量截面及截面测量点的缺陷,既提高了测量效率又可保证了测量精度。
下面结合附图和实施例对本发明作详细说明。
附图说明
图1是本发明未知CAD模型的航空发动机叶片型面测量方法的流程图。
图2是曲率代表线测点集。
图3是曲率代表线拟合曲面及中分线求取。
图4是中分线最小控制点集求取。
图5是Z=72.92mm处叶片截面测量点集。
图6是叶片型面测量点集。
具体实施方式
参照图1-6。以某型航空发动机压气机叶片型面为例详细说明本发明。
步骤1:分析叶片结构特点,沿叶片主轴方向在叶盆、叶背曲面各规划五条曲线,曲线在叶盆叶背曲面呈等间距分布。
步骤2:对规划的曲线进行数据点采集,根据主轴方向叶片尺寸,采用手动打点的方式,在叶盆、叶背型面上各段曲线均采集30个点,结果如图2所示。
步骤3:将上述曲线测点集数据导入UG软件中并拟合曲线,将叶盆、叶背曲面上对应曲线进行曲面拟合,并提取曲面中分线,结果如图3所示。
步骤4:对五条中分线进行基于曲率分布及容差控制的最小控制点集提取,提取的步骤如下;
[1]任取一条中分线,以曲线的任一端点为起点,向前搜索一微小弧长,计算搜索始末点连线的中点到该段搜索曲线的距离d,设h为曲线弦公差,若|d-h|<ε,(h=nc,n为大于0的自然数,c为曲线弦公差单位变化值,考虑计算精度和效率,c设为0.005;ε为设定阈值,工程上一般设置为10-3),则该点为控制点;若h-d>ε,则继续向前搜索。搜索完毕,将搜索到的点拟合曲线,记原始中分线容差值为S;将拟合曲线进行一定数量的点离散,并计算出离散点到原始曲线的最大距离dmax,若|S-dmax|>ε,则调整n,重复上述操作,直到|S-dmax|<ε。
[2]取其中一条中分线为参考线,顺次过其上控制点作叶片主轴垂面与其余等参数曲线交于一点,计算其余曲线上离交点最近的控制点,如果该控制点与交点Z向坐标差值<nr(n为自然数,r为测头半径,ΔZmin为第一小步所求每条曲线各控制点Z向坐标差值的最小值),则将该控制点以交点代替;如果两点Z向坐标差值>nr,则将交点加入该曲线的控制点集;如此,第一条中分线上控制点遍历结束,再以第二条中分线为参考线重复上述操作,直到所有中分线遍历结束,即可初步实现全体中分线的控制点集统一。
[3]以得到的控制点集拟合B样条曲线,设定原始中分线容差值为S;将拟合曲线进行一定数量的点离散,并计算出离散点到原始中分线的最大距离dmax,若|S-dmax|>ε(ε为设计阈值,此处取10-3),则调整第二小步中n,直到|S-dmax|<ε。
中分线弦公差取值依次为(图5中从左到右),第一条为0.065,第二条为0.08,第三条为0.09,第四条为0.07,第五条为0.06,所求控制点个数为19;控制点提取结果如图4所示。
步骤5:由上一步骤可确定叶片测量截面。以Z=72.92mm处叶型曲线为例,其控制点的提取步骤如下。
[1]利用三坐标测量机在测量截面上初测6个点,并指定搜索方向。
[2]利用前6点的坐标信息拟合5次多项式。
[3]求解第6点一阶、二阶导数和曲率半径;
[4]根据第6点切线方向和曲率半径,预测第7点坐标值。
[5]判断预测点坐标是否回到起始点,未到达则根据第7点预测值引导三坐标测量机测量并获取第7点真实坐标,若到达或超出起始点则搜索结束;测量结果如图5所示。
步骤5:对所有测量截面实施步骤4操作,结果如图6所示;将测点集数据导入UG,根据各测量截面测点信息拟合出测量截面曲线,并由各测量截面曲线拟合出叶身型面。
步骤6:将原始规划的叶身曲线与重构模型上相应曲线进行比较,将重构模型曲线以等弧长方式离散为1000个点,离散点到原始规划曲线最大距离均小于0.005mm,符合工程精度要求。
Claims (1)
1.一种未知CAD模型的航空发动机叶片型面测量方法,其特征在于包括以下步骤:
步骤一、沿叶片主轴方向在叶盆、叶背型面划出若干条曲线,曲线在叶盆、叶背曲面呈等间距分布,用三坐标测量机采集各规划曲线上若干点,得到测点集信息;
步骤二、将步骤一得到的测点集信息导入UG,并拟合曲线,将叶盆、叶背曲面上对应的两条曲线拟合曲面,并提取该曲面中分线;
步骤三、基于曲率分布及容差控制的中分线最小控制点集提取;
第一步,任取一条中分线,以曲线的任一端点为起点,向前搜索一微小弧长,计算搜索始末点连线的中点到该段搜索曲线的距离d,设h为曲线弦公差,若|d-h|<ε,h=nc,n为大于0的自然数,c为曲线弦公差单位变化值,考虑计算精度和效率,c设为0.005;ε为设定阈值,则该点为控制点;若h-d>ε,则继续向前搜索;搜索完毕,将搜索到的点拟合曲线,记原始中分线容差值为S;将拟合曲线进行一定数量的点离散,并计算出离散点到原始曲线的最大距离dmax,若|S-dmax|>ε,则调整n,重复第一步操作,直到|S-dmax|<ε;
第二步,取其中一条中分线为参考线,顺次过其上控制点作叶片主轴垂面与其余参数曲线交于一点,计算其余曲线上离交点最近的控制点,如果该控制点与交点Z向坐标差值<nr,n为自然数,r为测头半径,ΔZmin为第一步所求每条曲线各控制点Z向坐标差值的最小值,则将该控制点以交点代替;如果两点Z向坐标差值>nr,则将交点加入该曲线的控制点集;第一条中分线上控制点遍历结束,再以第二条中分线为参考线重复第二步操作,直到所有中分线遍历结束,初步实现全体中分线的控制点集统一;
第三步,以得到的控制点集拟合B样条曲线,设定原始中分线容差值为S;将拟合曲线进行点离散,并计算出离散点到原始中分线的最大距离dmax,若|S-dmax|>ε,则调整第二步中n,直到|S-dmax|<ε;
步骤四、基于连续曲率预测法提取测量截面测点集;
根据自由曲面多阶连续的特性,利用曲率连续预测法实现测点的自适应规划,建立拟合多项式Z=AX5+BX4+CX3+DX2+EX+F;
利用三坐标测量机在求出的测量截面上初测6个点,并指定搜索方向;
利用前6点的坐标信息拟合5次多项式;
求解第6点一阶、二阶导数和曲率半径;
根据第6点切线方向和曲率半径,预测第7点坐标值;
判断预测点坐标是否回到起始点,未到达则根据第7点预测值引导三坐标测量机测量并获取第7点真实坐标,若到达或超出起始点则搜索结束;
步骤五、将叶片型面测点集数据导入UG,根据各测量截面测点信息拟合出测量截面曲线,并由各测量截面曲线拟合出叶身型面;
步骤六、在重构模型上提取曲率代表线,并与第一步中测量点拟合的曲率代表线进行比较,将重构模型曲率代表线以等弧长方式离散为若干点,并计算各离散点到对应测量点拟合曲率代表线的距离,若大于给定容差值,则增加步骤一中规划的曲率代表线条数,重复上述操作,直到重构模型曲率代表线各离散点到原始曲率代表线的距离小于给定容差值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310354360.7A CN103486996B (zh) | 2013-08-14 | 2013-08-14 | 未知cad模型的航空发动机叶片型面测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310354360.7A CN103486996B (zh) | 2013-08-14 | 2013-08-14 | 未知cad模型的航空发动机叶片型面测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103486996A CN103486996A (zh) | 2014-01-01 |
CN103486996B true CN103486996B (zh) | 2016-01-20 |
Family
ID=49827412
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310354360.7A Expired - Fee Related CN103486996B (zh) | 2013-08-14 | 2013-08-14 | 未知cad模型的航空发动机叶片型面测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103486996B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110928326A (zh) * | 2019-11-26 | 2020-03-27 | 南京航空航天大学 | 一种面向飞行器外形的测量点差异性规划方法 |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104180789B (zh) * | 2014-09-12 | 2017-02-22 | 北京航空航天大学 | 基于图形匹配算法的叶片检测方法 |
CN104462807B (zh) * | 2014-12-05 | 2017-10-27 | 武汉船用机械有限责任公司 | 一种精密铸造式喷推叶轮的叶片型值检验方法 |
CN104864829B (zh) * | 2015-06-14 | 2017-06-27 | 吉林大学 | 一种叶片曲面的快速测量方法 |
CN105069312B (zh) * | 2015-08-24 | 2018-02-09 | 北京星航机电装备有限公司 | 自由曲面测量点分布的规划方法 |
CN107806836A (zh) * | 2017-09-28 | 2018-03-16 | 中国航发动力股份有限公司 | 一种基于五轴联动点光源的精锻叶片型面检测方法 |
CN110738728B (zh) * | 2019-10-09 | 2022-05-03 | 西北工业大学 | 基于线性组合过渡算法的叶片修复模型重建方法 |
CN110703686B (zh) * | 2019-11-11 | 2022-12-27 | 苏州千机智能技术有限公司 | 一种整体叶盘叶片截面在线测量路径规划方法 |
CN110969697B (zh) * | 2019-12-20 | 2023-09-19 | 杭州汽轮控股有限公司 | 一种燃气轮机涡轮动叶片关键结构检测装置和处理方法 |
CN110966969B (zh) * | 2019-12-20 | 2021-05-14 | 杭州汽轮动力集团有限公司 | 一种燃气轮机涡轮静叶片关键结构检测装置和处理方法 |
CN111400667B (zh) * | 2020-03-31 | 2021-11-02 | 华中科技大学 | 一种基于变公差带约束的航空叶片型面检测方法和*** |
CN111907727B (zh) * | 2020-06-17 | 2022-04-08 | 成都飞机工业(集团)有限责任公司 | 一种飞机机翼骨架外形测量点选取方法 |
CN111678477A (zh) * | 2020-06-20 | 2020-09-18 | 贵阳航发精密铸造有限公司 | 一种涡轮工作叶片最终检验自动检测测量方法 |
CN112414354A (zh) * | 2020-12-15 | 2021-02-26 | 西安飞机工业(集团)有限责任公司 | 一种对曲面零件的外形测量点规划方法 |
CN112729129A (zh) * | 2021-02-03 | 2021-04-30 | 西安工业大学 | 一种四轴叶片激光测量平台的点激光测量装置和方法 |
CN113867258B (zh) * | 2021-09-18 | 2023-09-01 | 华中科技大学 | 一种基于在机测量的航空叶片加工定位方法 |
CN113983991B (zh) * | 2021-10-11 | 2024-07-19 | 中国航发南方工业有限公司 | 一种径向扩压器的叶尖圆弧测量方法 |
CN113970311A (zh) * | 2021-10-12 | 2022-01-25 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种航空发动机叶片矢量逼近迭代测量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7379193B2 (en) * | 2002-02-09 | 2008-05-27 | Lang Liu | Sensing device for measuring the three-dimension shape and its measuring method |
CN102155928A (zh) * | 2011-03-03 | 2011-08-17 | 西北工业大学 | 用于叶片三坐标测量机的夹具及建立测量坐标系的方法 |
CN102323962A (zh) * | 2011-05-20 | 2012-01-18 | 西北工业大学 | 航空发动机叶片型面流道线测量路径规划方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4491687B2 (ja) * | 2005-04-21 | 2010-06-30 | パルステック工業株式会社 | 座標変換関数の補正方法 |
-
2013
- 2013-08-14 CN CN201310354360.7A patent/CN103486996B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7379193B2 (en) * | 2002-02-09 | 2008-05-27 | Lang Liu | Sensing device for measuring the three-dimension shape and its measuring method |
CN102155928A (zh) * | 2011-03-03 | 2011-08-17 | 西北工业大学 | 用于叶片三坐标测量机的夹具及建立测量坐标系的方法 |
CN102323962A (zh) * | 2011-05-20 | 2012-01-18 | 西北工业大学 | 航空发动机叶片型面流道线测量路径规划方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110928326A (zh) * | 2019-11-26 | 2020-03-27 | 南京航空航天大学 | 一种面向飞行器外形的测量点差异性规划方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103486996A (zh) | 2014-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103486996B (zh) | 未知cad模型的航空发动机叶片型面测量方法 | |
CN103411574B (zh) | 航空发动机叶片型面三坐标测量方法 | |
CN104316016B (zh) | 一种整体叶盘复杂曲面叶片的纵向测量方法 | |
CN101672637B (zh) | 一种复杂曲面的数字化检测方法 | |
CN104864829B (zh) | 一种叶片曲面的快速测量方法 | |
CN102735204B (zh) | 一种基于弦线的航空薄壁叶片加工扭曲度误差测量方法 | |
CN108682043A (zh) | 一种基于参数映射的复杂曲面测量规划方法 | |
CN102607484B (zh) | 一种变间距自适应搜寻测量点的飞机壁板法矢确定方法 | |
CN102622479B (zh) | 一种基于三维草图的逆向工程cad建模方法 | |
CN102049731B (zh) | 对盘形凸轮进行精确轮廓测量及精确加工的方法 | |
CN109683552B (zh) | 一种基面曲线导向的复杂点云模型上的数控加工路径生成方法 | |
CN107526875A (zh) | 一种航空叶片型面加工误差可视化方法 | |
CN106446343B (zh) | 一种径流叶轮直纹叶片参数化型线的自动提取方法 | |
CN103942837B (zh) | 基于序列线性规划的叶片点云模型截面曲线直接构造方法 | |
CN112033338B (zh) | 一种叶片类曲面接触式扫描测量测头半径面补偿方法 | |
CN101281403A (zh) | 异型玻璃加工图形定位方法及*** | |
CN105446264A (zh) | 基于特征的机床精度优化设计方法 | |
CN103900520A (zh) | 一种整体壁板板坯模型几何尺寸检测方法 | |
CN105045973A (zh) | 一种变弧长自适应采样方法 | |
CN104280000A (zh) | 凸轮型线的在线检测方法 | |
CN114912159A (zh) | 一种轨道交通线路平面几何线形的拟合方法 | |
CN110928326A (zh) | 一种面向飞行器外形的测量点差异性规划方法 | |
CN107103136A (zh) | 一种基于距离场的轨道车辆轮轨接触点计算方法 | |
WO2021196407A1 (zh) | 基于割线旋转迭代的航空叶片叶型弦长检测方法和*** | |
CN111806720B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 Termination date: 20160814 |
|
CF01 | Termination of patent right due to non-payment of annual fee |