CN103196666A - 一种基于等强度原则的水平轴风力机叶片缩比模型静力测试方法 - Google Patents

一种基于等强度原则的水平轴风力机叶片缩比模型静力测试方法 Download PDF

Info

Publication number
CN103196666A
CN103196666A CN2013101004255A CN201310100425A CN103196666A CN 103196666 A CN103196666 A CN 103196666A CN 2013101004255 A CN2013101004255 A CN 2013101004255A CN 201310100425 A CN201310100425 A CN 201310100425A CN 103196666 A CN103196666 A CN 103196666A
Authority
CN
China
Prior art keywords
blade
scale
load
model
full
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
Application number
CN2013101004255A
Other languages
English (en)
Other versions
CN103196666B (zh
Inventor
石可重
徐建中
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Engineering Thermophysics of CAS
Original Assignee
Institute of Engineering Thermophysics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Engineering Thermophysics of CAS filed Critical Institute of Engineering Thermophysics of CAS
Priority to CN201310100425.5A priority Critical patent/CN103196666B/zh
Publication of CN103196666A publication Critical patent/CN103196666A/zh
Application granted granted Critical
Publication of CN103196666B publication Critical patent/CN103196666B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种水平轴风力机叶片缩比模型的静力测试方法,其目的在于通过采用低成本、高效率的缩比模型试验检测与计算分析手段,替代全尺寸叶片静力测试试验。该方法特点在于,采用尺寸等比例缩小的叶片模型静力试验替代传统的全尺寸叶片静力试验,并依据缩比模型与原型叶片满足结构等强度这一要求,用以确定缩比模型的试验载荷。该测试方法既可采用行业现有使用的叶片静力测试设备,也可采用低成本的小型试验装置和设备,进行静力测试。通过准确的设计与数值计算,可实现缩比模型与全尺寸叶片有着比较准确的结构尺寸对应关系与载荷对应关系,以及测试结果的准确性。

Description

一种基于等强度原则的水平轴风力机叶片缩比模型静力测试方法
技术领域
本发明涉及机械行业风力发电技术领域,尤其涉及一种水平轴风力机叶片在静态载荷作用下结构性能的测试评估方法。
背景技术
风电叶片在极限载荷作用下发生损伤破坏问题是风电叶片设计、制造、运行中所关注的一项重要内容。在风电叶片设计中,一般都要求风电叶片能满足20年的使用寿命。在如此长的服役期间和恶劣的工作环境下,风电机组经常要在多种极限载荷作用下运行,严重影响了风电机组安全运行的可靠性和使用寿命。而恶劣的工作环境,特殊的材料性质,以及结构与工艺所带来的种种问题,给叶片极限载荷下性能评估带来相当大的难度。
由于风电叶片结构及运行情况的特殊性,其它行业已建立的工程实践经验和知识不能有效的应用到风电机组叶片的性能评估中来。进行叶片在极限载荷作用下的静力测试可以提供对设计的可靠确认。对于新设计、新工艺、新材料的风电叶片,静力测试是保证叶片质量关键的一环。通过测试,可以将测量数据与设计的数据进行比较,以掌握叶片的结构性能。并且可以利用测试的结果改进与优化结构设计。
然而,目前所开展的叶片静力测试方法大都为全尺寸叶片测试,试验叶片的制作与试验周期很长,而且大都采用大型试验设备进行,成本很高。降低叶片测试成本,缩短测试周期,对于叶片行业具有重要的现实意义。
发明内容
针对现有叶片静力测试的以上问题,本专利提出一种利用叶片缩比模型静力测试试验,代替全尺寸叶片静力测试的叶片静力性能测试评估方法。该方法具有周期短、成本低等优点。此外,本方法中缩比模型试验所需的设备与实验步骤,与现有全尺寸叶片静力测试方法相似。该方法即可利用现有国内外已有的叶片静态检测设施开展测试,也可采用规模形式较小的试验设备开展,因此在设备投资上也是比较经济的。
(一)要解决的技术问题
针对目前风力机叶片静力测试测试周期长、成本高等缺点,本发明提出了一种缩比模型静力测试方法,通过对缩比模型结构与试验载荷准确的设计与计算,可保证缩比模型测试具有比较准确的测试结果。采用该方法进行静力测试,可实现周期短、成本低、能够获得比较准确的叶片静态加载情况下结构性能数据等优点。
(二)技术方案
根据本发明的一个方面,提供了一种利用缩比模型试验替代全尺寸风力机叶片静力测试的方法,其特征在于,该方法包括以下步骤:
1)依据所要测试评估的原型全尺寸叶片结构特征,设计并制作一个结构尺寸等比例缩小的缩尺寸模型叶片,并以该叶片作为试验叶片;
2)依据等强度原则,(即缩尺寸模型叶片与原型全尺寸叶片在各自的试验载荷作用下,应保持相同的结构应力分布),将原型全尺寸叶片的试验载荷转化为缩尺寸模型叶片的试验载荷。其机理与方法如下:
为了实现用缩比模型试验替代全尺寸叶片静力试验,首先要确定两个试验的相似准则。由于静力试验的核心目标是研究叶片在外部极限载荷作用下,叶片是否发生损坏,以及一旦发生损坏,损坏将以何种模式存在。因此,静力试验的要求是要能够实现对外部极限载荷作用下,叶片结构应力分布情况的准确描述。根据这一要求,在本发明中,缩比模型试验载荷确定的原则为,缩比模型与原型叶片在各自的试验载荷作用下,应保持相同的结构应力分布。
然而,要实现这一目标,首要解决的是外部静力载荷与叶片应力分布关系的建立问题。由于风电叶片是由外部壳体与内部大梁组成的复合材料结构体,沿叶片展向,外部壳体的形状、尺寸与铺层结构都在发生变化,而内部大梁的结构形状、尺寸与铺层结构也都随之发生变化。这种复杂的结构形式,难以采用解析方法建立载荷与叶片应力分布的数学关系。这也成为实现缩比模型试验替代全尺寸叶片试验的难点。
本发明则提出利用数值分析的方法,确定缩比模型试验的载荷。其方法如下:
(a)采用有限元方法,分别建立全尺寸叶片与缩比模型的有限元模型;
(b)对两个模型分别定义一致的边界条件;例如,由于开展台架疲劳试验时,叶片的结构边界条件只包括叶根处的约束,因此有限元模型也只需定义叶根约束即可。
(c)根据设计部门所提供的全尺寸叶片试验载荷,施加到全尺寸叶片有限元模型上,一般为沿叶片展向布置的几处集中载荷;并据此计算在该载荷作用下叶片结构的应力分布情况。
(d)将(c)中全尺寸叶片各加载点处的载荷,按照一个相同的比例缩小,(例如按照缩比模型与原型叶片的尺寸缩小比例值),并将缩小后的载荷施加到缩比模型的有限元模型上,采取与(c)中一致的计算方法,计算在该载荷作用下叶片结构应力分布情况。
(e)分别在全尺寸叶片和缩比模型的对应位置选取若干位置点,作为应力对比点,选取的原则为,利用这些点处的强度计算结果,可以较好反映或识别叶片的损伤情况。
(f)比较全尺寸叶片和缩比模型对比点处的应力情况,如应力不一致,可调整缩比模型的载荷,(调整时须保证缩比模型上各点处的载荷按照相同比例,同步地增减)。
(g)经过若干次调整计算,可以找到实现对比点处应力一致(或近似)的缩比模型施加载荷,以该载荷作为缩比模型静力试验的试验载荷。
3)完成以上步骤后,即可进行缩比模型的静力试验。包括将叶片模型固定在试验台架上,进行加载,测量等步骤。所采用的加载装置与试验过程可与现有通用的叶片静力试验方法相同或类似,这里不再表述;
4)根据缩比模型静力试验获得的应力测量结果,对原型叶片的结构性能进行评估。
优选的,采用该方法进行叶片静力测试,其前提在于依据原型叶片设计并制作等比例缩小的试验叶片。
优选的,采用该方法进行叶片静力性能评估,应依据等强度原则,采用有限元数值模拟方法,对原型叶片与缩比模型进行准确的强度计算,用以确定缩比模型的试验载荷。
优选的,采用该方法进行叶片静力性能评估,试验步骤及试验装置可利用现有通用的叶片静态测试步骤。
(三)有益效果
本发明与现有技术相比,具有如下明显的实质特点和显著优点:
1)与现有采用全尺寸叶片进行静力载荷作用下结构性能测试评估相比,该方法采用缩比模型进行静力载荷作用下性能测试评估,可实现周期短、成本低、能够比较准确获得叶片静态结构性能等优点。
2)该方法可利用现有国内外叶片静态检测设施,或更小型的试验设备,因此在设备投资上也是比较经济的。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,对本发明进一步详细说明。
本实施例中,以当前主流的MW级风力机叶片为例,采用以下步骤进行检测:
1)依据所要评估的原型叶片结构特征,设计并制作一个结构尺寸等比例缩小的叶片,并以该叶片作为试验叶片,本例中缩小比例取为10;
2)依据等强度原则,(即缩尺寸模型叶片与原型叶片在各自的试验载荷作用下,应保持相同的结构应力分布),将原型叶片的试验载荷转化为缩尺寸模型叶片的试验载荷。其步骤如下:
(a)采用有限元方法,分别建立全尺寸叶片与缩比模型的有限元模型;
(b)对两个模型分别定义一致的边界条件;
(c)根据设计部门所提供的全尺寸叶片试验载荷,施加到全尺寸叶片有限元模型上,一般为沿叶片展向布置的几处集中载荷;并据此计算在该载荷作用下叶片结构的应力分布情况。
(d)将(c)中全尺寸叶片各加载点处的载荷,按照一个相同的比例缩小,本例中缩小比例取为10,并将缩小后的载荷施加到缩比模型的有限元模型上,采取与(c)中一致的计算方法,计算在该载荷作用下叶片结构应力分布情况。
(e)分别在全尺寸叶片和缩比模型的对应位置选取若干位置点,作为应力对比点,选取的原则为,利用这些点处的强度计算结果,可以较好反映或识别叶片的损伤情况。
(f)比较全尺寸叶片和缩比模型对比点处的应力情况,如应力不一致,可调整缩比模型的载荷,(调整时须保证缩比模型上各点处的载荷按照相同比例,同步地增减)。
(g)经过若干次调整计算,可以找到实现对比点处应力一致(或近似)的缩比模型施加载荷,以该载荷作为缩比模型静力试验的试验载荷。
3)将该缩比模型叶片固定在叶片实验台架上,在叶片表面及内部设置应变片,本实例中传感器的数量为30片,布置位置尽可能地覆盖叶片在极限载荷下的危险发生位置(危险发生位置根据叶片设计得到);
4)依据2)中所确定的试验载荷,对叶片进行静态加载实验,实验过程与现有通用的全尺寸静力试验过程相似,这里不介绍。
5)对实验所采集的应变等结果进行分析,可获得原型叶片的相关静态结构力学性能数据。
通过以上步骤的实施,完成了对该叶片静态加载性能的测试与评估。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明。所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种水平轴风力机叶片的静力测试方法,其特征在于,所述方法包括以下步骤:
1)依据所要测试评估的原型全尺寸叶片结构特征,设计并制作一个材料相同而结构尺寸等比例缩小的缩尺寸模型叶片,并以该叶片作为试验叶片;
2)依据等强度原则,(即缩尺寸模型叶片与原型全尺寸叶片在各自的试验载荷作用下,应保持相同的结构应力分布),将原型全尺寸叶片的试验载荷转化为缩尺寸模型叶片的试验载荷,具体步骤如下:
(a)采用有限元数值模拟方法,分别建立原型全尺寸叶片与缩尺寸模型叶片的有限元模型;
(b)对上述两个有限元模型分别定义一致的边界条件;
(c)根据设计部门所提供的原型全尺寸叶片试验载荷,施加到原型全尺寸叶片有限元模型上,并据此计算在该载荷作用下原型全尺寸叶片的结构应力分布情况;
(d)将步骤(c)中原型全尺寸叶片各加载点处的载荷,按照一个相同的比例缩小,并将缩小后的载荷施加到缩尺寸模型叶片的有限元模型的对应各加载点处,采取与步骤(c)中一致的计算方法,计算在该载荷作用下缩尺寸模型叶片的结构应力分布情况;
(e)分别在原型全尺寸叶片和缩尺寸模型叶片的对应位置选取若干位置点,作为应力对比点;
(f)比较原型全尺寸叶片和缩尺寸模型叶片在上述应力对比点处的应力分布情况,如应力分布不一致,则调整施加到缩尺寸模型叶片有限元模型上的载荷;
(g)经过若干次调整计算,找到实现原型全尺寸叶片和缩尺寸模型叶片在应力对比点处应力分布一致或近似一致的缩尺寸模型叶片的施加载荷,以该载荷作为缩尺寸模型叶片静力试验的试验载荷;
3)按照步骤2)所确定的施加到缩尺寸模型叶片的试验载荷,对缩尺寸模型叶片进行静力试验;
4)根据步骤3)中缩尺寸模型叶片静力试验获得的应力测量结果,对原型全尺寸叶片的结构性能进行评估。
2.根据权利要求1所述的水平轴风力机叶片缩比模型静力测试方法,其特征在于,所述步骤2)的子步骤(c)中,施加到原型全尺寸叶片有限元模型上的试验载荷,为沿叶片展向布置的几处集中载荷。
3.根据权利要求1或2所述的水平轴风力机叶片的缩比模型静力测试方法,其特征在于,所述步骤2)的步骤(d)中,按照原型全尺寸叶片与缩尺寸模型叶片的尺寸缩小比例值,缩小步骤(c)中原型全尺寸叶片各加载点处的载荷,得到施加到缩尺寸模型叶片有限元模型上的载荷。
4.根据上述任一项权利要求所述的水平轴风力机叶片的缩比模型静力测试方法,其特征在于,所述步骤2)的步骤(e)中,依据可较好反映或识别叶片损伤情况的强度计算结果,分别在原型全尺寸叶片和缩尺寸模型叶片的对应位置选取若干位置点,作为应力对比点。
5.根据上述任一项权利要求所述的水平轴风力机叶片的缩比模型静力测试方法,其特征在于,所述步骤2)的步骤(f)中,调整施加到缩尺寸模型叶片上的载荷时,须保证施加到各位置点处的载荷按照相同比例同步增减。
6.根据上述任一项权利要求所述的水平轴风力机叶片的缩比模型静力测试方法,其特征在于,所述步骤2)中,对原型全尺寸叶片与缩尺寸模型叶片的有限元模型采用一致或相似的计算网格以及相同的边界条件进行数值计算,以使二者具备一致或相似的计算误差。
7.根据上述任一项权利要求所述的水平轴风力机叶片的缩比模型静力测试方法,其特征在于,采用该方法进行叶片静力性能评估,试验步骤及试验装置可利用现有通用的叶片静态测试步骤。
CN201310100425.5A 2013-03-26 2013-03-26 一种基于等强度原则的水平轴风力机叶片缩比模型静力测试方法 Active CN103196666B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310100425.5A CN103196666B (zh) 2013-03-26 2013-03-26 一种基于等强度原则的水平轴风力机叶片缩比模型静力测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310100425.5A CN103196666B (zh) 2013-03-26 2013-03-26 一种基于等强度原则的水平轴风力机叶片缩比模型静力测试方法

Publications (2)

Publication Number Publication Date
CN103196666A true CN103196666A (zh) 2013-07-10
CN103196666B CN103196666B (zh) 2015-04-22

Family

ID=48719360

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310100425.5A Active CN103196666B (zh) 2013-03-26 2013-03-26 一种基于等强度原则的水平轴风力机叶片缩比模型静力测试方法

Country Status (1)

Country Link
CN (1) CN103196666B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104134013A (zh) * 2014-08-16 2014-11-05 中国科学院工程热物理研究所 一种风力机叶片模态分析方法
CN104155092A (zh) * 2014-08-16 2014-11-19 中国科学院工程热物理研究所 一种风力机叶片静态分析方法
CN105203310A (zh) * 2015-09-17 2015-12-30 中国飞机强度研究所 一种航空发动机的轴类部件载荷标定方法及系统
CN105628356A (zh) * 2015-12-22 2016-06-01 中国航空工业集团公司济南特种结构研究所 大曲率雷达罩静力试验加载帆布带精确定位方法
CN106840721A (zh) * 2016-12-15 2017-06-13 中国航空工业集团公司西安飞机设计研究所 一种飞行器结构传载特性模型试验设计方法
CN107084844A (zh) * 2017-04-19 2017-08-22 西北工业大学 一种叶片结构模拟件的试验方法
CN112434456A (zh) * 2020-11-26 2021-03-02 国网江西省电力有限公司建设分公司 一种输电线路塔型力学薄弱环节计算及加固方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101634604A (zh) * 2009-09-07 2010-01-27 中国人民解放军国防科学技术大学 风电用大型叶片的静强度测试方法及其测试系统
CN201945447U (zh) * 2010-12-27 2011-08-24 中航惠腾风电设备股份有限公司 一种大型风轮叶片静力试验的多点加载侧拉试验装置
CN102175511A (zh) * 2010-12-24 2011-09-07 深圳市大族激光科技股份有限公司 材料性能评估方法和系统
CN102455264A (zh) * 2010-12-27 2012-05-16 中航惠腾风电设备股份有限公司 一种大型风轮叶片静力试验的多点加载侧拉试验装置
CN202274981U (zh) * 2011-10-12 2012-06-13 华锐风电科技(集团)股份有限公司 大功率叶片试验台静力加载装置
CN102607831A (zh) * 2012-02-25 2012-07-25 中国科学院工程热物理研究所 一种水平轴风力机叶片疲劳损伤与寿命评估方法
EP2522975A2 (en) * 2011-03-07 2012-11-14 Vestas Wind Systems A/S A wind turbine blade tester

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101634604A (zh) * 2009-09-07 2010-01-27 中国人民解放军国防科学技术大学 风电用大型叶片的静强度测试方法及其测试系统
CN102175511A (zh) * 2010-12-24 2011-09-07 深圳市大族激光科技股份有限公司 材料性能评估方法和系统
CN201945447U (zh) * 2010-12-27 2011-08-24 中航惠腾风电设备股份有限公司 一种大型风轮叶片静力试验的多点加载侧拉试验装置
CN102455264A (zh) * 2010-12-27 2012-05-16 中航惠腾风电设备股份有限公司 一种大型风轮叶片静力试验的多点加载侧拉试验装置
EP2522975A2 (en) * 2011-03-07 2012-11-14 Vestas Wind Systems A/S A wind turbine blade tester
CN202274981U (zh) * 2011-10-12 2012-06-13 华锐风电科技(集团)股份有限公司 大功率叶片试验台静力加载装置
CN102607831A (zh) * 2012-02-25 2012-07-25 中国科学院工程热物理研究所 一种水平轴风力机叶片疲劳损伤与寿命评估方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
毛火军等: "基于CFD和BEM方法的风电叶片强度分析比较", 《太阳能学报》 *
毛火军等: "大型风电叶片的模态测试与数值模拟", 《工程热物理学报》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104134013A (zh) * 2014-08-16 2014-11-05 中国科学院工程热物理研究所 一种风力机叶片模态分析方法
CN104155092A (zh) * 2014-08-16 2014-11-19 中国科学院工程热物理研究所 一种风力机叶片静态分析方法
CN104134013B (zh) * 2014-08-16 2017-02-08 中国科学院工程热物理研究所 一种风力机叶片模态分析方法
CN104155092B (zh) * 2014-08-16 2017-05-24 中国科学院工程热物理研究所 一种风力机叶片静态分析方法
CN105203310A (zh) * 2015-09-17 2015-12-30 中国飞机强度研究所 一种航空发动机的轴类部件载荷标定方法及系统
CN105203310B (zh) * 2015-09-17 2018-04-13 中国飞机强度研究所 一种航空发动机的轴类部件载荷标定方法及系统
CN105628356A (zh) * 2015-12-22 2016-06-01 中国航空工业集团公司济南特种结构研究所 大曲率雷达罩静力试验加载帆布带精确定位方法
CN105628356B (zh) * 2015-12-22 2018-08-31 中国航空工业集团公司济南特种结构研究所 大曲率雷达罩静力试验加载帆布带精确定位方法
CN106840721A (zh) * 2016-12-15 2017-06-13 中国航空工业集团公司西安飞机设计研究所 一种飞行器结构传载特性模型试验设计方法
CN106840721B (zh) * 2016-12-15 2019-06-11 中国航空工业集团公司西安飞机设计研究所 一种飞行器结构传载特性模型试验设计方法
CN107084844A (zh) * 2017-04-19 2017-08-22 西北工业大学 一种叶片结构模拟件的试验方法
CN112434456A (zh) * 2020-11-26 2021-03-02 国网江西省电力有限公司建设分公司 一种输电线路塔型力学薄弱环节计算及加固方法

Also Published As

Publication number Publication date
CN103196666B (zh) 2015-04-22

Similar Documents

Publication Publication Date Title
CN103196666B (zh) 一种基于等强度原则的水平轴风力机叶片缩比模型静力测试方法
CN103217282B (zh) 一种基于等寿命原则的水平轴风力机叶片缩比模型疲劳测试方法
CN102607831B (zh) 一种水平轴风力机叶片疲劳损伤与寿命测试方法
CN102410928B (zh) 一种水平轴风力机叶片的疲劳检测方法
CN107223206B (zh) 疲劳测试
CN107346357B (zh) 一种基于整体耦合模型的海上风机疲劳分析系统
Rašuo et al. Harmonization of new wind turbine rotor blades development process: A review
CN105205264B (zh) 预测抽水蓄能发电电动机转子构件使用寿命的试验方法
CN105701337A (zh) 风电机组的疲劳寿命预测方法和装置
CN110390146B (zh) 基于扇区载荷的风电机组塔筒焊缝疲劳损伤计算方法及产品
CN105893757A (zh) 一种用于电压暂降严重程度的综合评估方法
CN103196689B (zh) 一种水平轴风力机叶片静力测试结果修正方法
CN102564750A (zh) 一种水平轴风力机叶片模态测试结果修正方法
CN104155092B (zh) 一种风力机叶片静态分析方法
CN104134013B (zh) 一种风力机叶片模态分析方法
CN108333446A (zh) 基于风电变桨控制系统的超级电容检测方法
CN108869174A (zh) 一种非线性建模的风力发电机叶片固有频率工况补偿方法
CN205138791U (zh) 一种燃气涡轮发动机轴向力测试传感器装配状态检测装置
CN109798226B (zh) 一种风电机组塔筒载荷预测方法和系统
CN108009311A (zh) 一种针对蠕变试验的蠕变本构模型参数识别方法
CN105606391B (zh) 大型风机塔架强度试验装置及试验方法
WO2014207528A1 (en) Static testing of wind turbine blades
CN104573197B (zh) 一种水轮机模型参数的测试方法及系统
CN103514329A (zh) 变电站避雷针塔安全性评定方法
CN105067364A (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