CN109492275B - 一种风机主轴制动器强度校核方法及系统 - Google Patents

一种风机主轴制动器强度校核方法及系统 Download PDF

Info

Publication number
CN109492275B
CN109492275B CN201811237654.0A CN201811237654A CN109492275B CN 109492275 B CN109492275 B CN 109492275B CN 201811237654 A CN201811237654 A CN 201811237654A CN 109492275 B CN109492275 B CN 109492275B
Authority
CN
China
Prior art keywords
main shaft
dimensional geometric
bolt
shaft brake
finite element
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
Application number
CN201811237654.0A
Other languages
English (en)
Other versions
CN109492275A (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.)
Xuji Group Co Ltd
Xuchang Xuji Wind Power Technology Co Ltd
Original Assignee
Xuji Group Co Ltd
Xuchang Xuji Wind Power Technology Co Ltd
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 Xuji Group Co Ltd, Xuchang Xuji Wind Power Technology Co Ltd filed Critical Xuji Group Co Ltd
Priority to CN201811237654.0A priority Critical patent/CN109492275B/zh
Publication of CN109492275A publication Critical patent/CN109492275A/zh
Application granted granted Critical
Publication of CN109492275B publication Critical patent/CN109492275B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02BCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
    • Y02B10/00Integration of renewable energy sources in buildings
    • Y02B10/30Wind power
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Braking Arrangements (AREA)

Abstract

本发明提供了一种风机主轴制动器强度校核方法及系统,首先建立风机主轴制动器各部件的三维几何模型,三维几何模型包括活动钳体、固定钳体、上挡板、下挡板及螺栓,不包括制动器的安装底座;然后将三维几何模型转换为有限元模型,进行有限元分析;再将螺栓预紧力施加到螺栓圆柱面上,并将液压力施加到活动钳体的活塞孔上,计算风机主轴制动器的应力变形情况。本发明的风机主轴制动器各部件的三维几何模型只考虑了关键因素,而未考虑非关键因素,降低了模型复杂性,提高了风机主轴制动器强度校核的效率与准确性,更加符合风电机组主轴制动器的实际运行情况。

Description

一种风机主轴制动器强度校核方法及系统
技术领域
本发明属于风力发电机组零件强度校核技术领域,特别涉及一种风机主轴制动器强度校核方法及系统。
背景技术
随着社会的发展,资源的不断枯竭,新能源行业越来越受到人民的青睐,这几年风力发电得到了快速发展。制动系统是风力发电机组重要的组成部分,在风力发电机组需要制动时起着极为关键的作用。风机主轴制动器是制动系统的一个重要部件,目前对于风机主轴制动器强度的校核方法主要采用有限元分析法,但是在采用有限元分析法对风机主轴制动器强度进行校核时,需要将主轴制动器相关部件全部建模分析,考虑的非关键因素过多,如公开号为“CN101893047A”,名称为“碟簧加载式液压盘式制动器”的中国发明专利申请,该专利在采用有限元结构对风机主轴制动器进行分析时考虑了很多因素,这将导致建立的有限元模型比较复杂,一些非关键因素被考虑进来,非关键因素导致的计算误差被进一步扩大,不能有效反映风机主轴制动器的情况。
发明内容
本发明的目的在于提供一种风机主轴制动器强度校核方法,用于解决现有技术中风机主轴制动器强度校核方法不能准确反映风机主轴制动器情况的问题;同时还提供了一种风机主轴制动器强度校核系统。
为实现上述目的,本发明提供了一种风机主轴制动器强度校核方法,包括如下步骤:
1)建立风机主轴制动器各部件的三维几何模型,所述三维几何模型包括活动钳体、固定钳体、上挡板、下挡板及螺栓,所述三维几何模型不包括制动器的安装底座;
2)将所述三维几何模型进行转换得到有限元模型,并对所述有限元模型进行有限元分析;
3)在所述有限元模型中对螺栓圆柱面施加螺栓预紧力,并在有限元模型中对活动钳体的活塞孔施加液压力,根据所述预紧力、液压力及由液压力产生的作用在上挡板、下挡板内侧面上的摩擦力,计算风机主轴制动器的等效应力。
本发明首先建立风机主轴制动器各部件的三维几何模型,三维几何模型包括活动钳体、固定钳体、上挡板、下挡板及螺栓,不包括制动器的安装底座;然后将三维几何模型转换为有限元模型,进行有限元分析;再将螺栓预紧力施加到螺栓圆柱面上,并将液压力施加到活动钳体的活塞孔上,计算风机主轴制动器的应力变形情况。本发明的风机主轴制动器各部件的三维几何模型只考虑了关键因素,而未考虑非关键因素,降低了模型复杂性,提高了风机主轴制动器强度校核的效率与准确性,更加符合风电机组主轴制动器的实际运行情况。
为了得到风机主轴制动器的安全系数及危险位置,根据风机主轴制动器的等效应力进一步计算风机主轴制动器的安全系数及危险位置,所述安全系数的计算公式为:
Figure BDA0001838575110000021
其中,rm为材料安全系数,取1.1;σeqv为等效应力值,单位为MPa;Re为材料屈服强度,单位为MPa。
为了使风机主轴制动器强度校核计算结果更加准确,在将所述三维几何模型转换为有限元模型时,需将所述三维几何模型中的圆角半径、螺纹孔去除,得到去除圆角半径、螺纹孔的三维几何模型。
为了使风机主轴制动器强度校核计算结果更加准确,在活动钳体、固定钳体、上挡板、下挡板及螺栓之间建立Bonded或Rough接触。
为了使风机主轴制动器强度校核计算结果更加准确,所述活动钳体、固定钳体之间通过螺栓连接,在所述活动钳体顶部端面施加Fixed Support约束,对所述固定钳体两侧圆柱孔施加圆柱约束。
本发明还提供了一种风机主轴制动器强度校核系统,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
1)建立风机主轴制动器各部件的三维几何模型,所述三维几何模型包括活动钳体、固定钳体、上挡板、下挡板及螺栓,所述三维几何模型不包括制动器的安装底座;
2)将所述三维几何模型进行转换得到有限元模型,并对所述有限元模型进行有限元分析;
3)在所述有限元模型中对螺栓圆柱面施加螺栓预紧力,并在有限元模型中对活动钳体的活塞孔施加液压力,根据所述预紧力、液压力及由液压力产生的作用在上挡板、下挡板内侧面上的摩擦力,计算风机主轴制动器的等效应力。
本发明首先建立风机主轴制动器各部件的三维几何模型,三维几何模型包括活动钳体、固定钳体、上挡板、下挡板及螺栓,不包括制动器的安装底座;然后将三维几何模型转换为有限元模型,进行有限元分析;再将螺栓预紧力施加到螺栓圆柱面上,并将液压力施加到活动钳体上,计算风机主轴制动器的应力变形情况。本发明的风机主轴制动器各部件的三维几何模型只考虑了关键因素,而未考虑非关键因素,降低了模型复杂性,提高了风机主轴制动器强度校核的效率与准确性,更加符合风电机组主轴制动器的实际运行情况。
为了得到风机主轴制动器的安全系数及危险位置,根据风机主轴制动器的等效应力进一步计算风机主轴制动器的安全系数及危险位置,所述安全系数的计算公式为:
Figure BDA0001838575110000031
其中,rm为材料安全系数,取1.1;σeqv为等效应力值,单位为MPa;Re为材料屈服强度,单位为MPa。
为了使风机主轴制动器强度校核计算结果更加准确,在将所述三维几何模型转换为有限元模型时,需将所述三维几何模型中的圆角半径、螺纹孔去除,得到去除圆角半径、螺纹孔的三维几何模型。
为了使风机主轴制动器强度校核计算结果更加准确,在活动钳体、固定钳体、上挡板、下挡板及螺栓之间建立Bonded或Rough接触。
为了使风机主轴制动器强度校核计算结果更加准确,所述活动钳体、固定钳体之间通过螺栓连接,在所述活动钳体顶部端面施加Fixed Support约束,对所述固定钳体两侧圆柱孔施加圆柱约束。
附图说明
图1是本发明的风机主轴制动器三维几何模型的示意图;
图2是本发明的零件间接触关系示意图;
图3是本发明的固定约束示意图;
图4是本发明的螺栓预紧示意图;
图5是本发明的液压力与水平作用力示意图;
图6是本发明的风机主轴制动器的三维视图;
图7为本发明的风机主轴制动器的活动钳体与固定钳体连接关系示意图;
图中:1-活动钳体,2-固定钳体,3-上挡板,4-下挡板,5-螺栓。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的说明:
本发明建立了风机主轴制动器的三维几何模型并进行装配,模型不包括摩擦片、活塞、密封圈等。将该模型导入有限元分析软件建立风机主轴制动器有限元模型。设置模型各部件的边界条件和接触,对模型施加工作载荷,求解得到风机主轴制动器的应力变形情况,进而得出风机主轴制动器的安全系数与危险位置。
具体的,本发明的风机主轴制动器强度校核方法,包括如下步骤:
1、如图1所示,建立风机主轴制动器各部件三维几何模型,该三维几何模型包括活动钳体1、固定钳体2、上挡板3、下挡板4、螺栓5,三维几何模型不包括固定制动器的安装底座,其约束通过边界条件施加。风机主轴制动器的各部件装配成风机主轴制动器,风机主轴制动器的各部件的真实连接关系如图6及图7所示。
2、将上述三维几何模型转换为有限元模型,并导入有限元分析软件进行有限元分析。为了简化网格、消除奇异网格以提高计算准确性,主要去除了三维几何模型中的个别细节,如圆角半径、螺纹孔。本实施例采用的有限元软件为ANSYS Workbench15.0,对风机主轴制动器的三维几何模型进行网络划分时,为了更好模拟螺栓5连接部位的受力传递情况,螺栓5采用实体建模,螺栓5及螺孔区域适当加密网格以提高网格质量。
3、在活动钳体1、固定钳体2、上挡板3、下挡板4、螺栓5之间建立Bonded或Rough接触,如图2所示。活动钳体1、固定钳体2之间通过螺栓5连接,由于真实情况下安装底座的固定作用和活动钳体1、固定钳体2之间螺栓5的连接作用,施加边界条件时只需对活动钳体1顶部端面施加Fixed Support约束Ux=Uy=Uz=0,对固定钳体2两侧圆柱孔施加圆柱约束Ux=Uz=0,如图3所示。
4、将螺栓预紧力通过预紧力选项Pretension直接施加到螺杆圆柱面上,如图4所示,并将此设为载荷步(1)。其中,螺栓预紧力的施加步骤为:在每个螺栓上都施加预紧力,螺栓预紧力通过Pretension载荷命令施加到螺杆圆柱面上,这一命令可以将螺杆在预紧力截面处断开,在这两部分一一对应的节点之间,创建一个Pre179单元,通过Pre179单元的被拉伸的情况来控制预紧力从而产生预紧力效果,根据作用力与反作用力原理,将压力通过螺栓帽与活动钳体1之间的接触传递到活动钳体1上。比如螺栓5受拉力,被螺栓5紧固的物体受到相同大小的压力,单个螺栓的预紧力为:
Figure BDA0001838575110000051
其中,d为螺纹部分危险剖面的计算直径;
σ0=(0.5~0.7)σs
其中,σs为螺栓材料的屈服极限。
5、在活动钳体1活塞孔的底部端面上施加液压力,并将此设为载荷步(2),具体的为:使用Force选项施加,其计算方法为液压力=工作压力×活塞孔底部作用面积。由于液压力的施加,液压力使得活动钳体1的上挡板3、下挡板4相向移动挤压制动盘产生摩擦力,并分别作用在左侧上、下挡板4的内侧面上,摩擦力=液压力×摩擦系数,因作用力与反作用力的关系,作用在固定钳体2上的等效集中力与活动钳体1上的等效应力相等,如图5所示。
6、分别求解载荷步(1)、载荷步(2),载荷步(1)计算的为预紧力作用下结构的应力与变形并将计算结果锁定之后参与载荷步(2)的计算,载荷步(2)计算的是液压力、摩擦力、预紧力共同作用下风机主轴制动器的应力与变形情况,计算等效应力,该等效应力指的是Von Mises等效应力,其是基于剪切应变能的一种等效应力,其表达式为:
Figure BDA0001838575110000052
其中,σx、σy、σz分别指正应力,τxy、τyz、τzx为对应的切应力。采用有限元计算方法,计算时在液压力、摩擦力、预紧力的共同作用下,制动器产生应力变形,经计算,有限元软件提取各个节点或单元上的各项正应力与切应力值,根据等效应力计算公式得到各个节点或单元上的等效应力值。根据最终得出的风机主轴制动器的等效应力值计算其安全系数,风机主轴制动器上每个位置上的等效应力都可以根据下述的公式计算出对应的安全系数值,风机主轴制动器上安全系数最小的位置即为等效应力最大的位置,也即为风机主轴制动器的危险位置,此安全系数最小的值即为风机主轴制动器的安全系数,风机主轴制动器的安全系数的表达式为:
Figure BDA0001838575110000061
其中,rm为材料安全系数,取1.1;σeqv为等效应力值,单位为MPa;Re为材料屈服强度,单位为MPa。
本发明还提供了一种风机主轴制动器强度校核系统,该校核系统是与上述校核方法相对应的进程和程序,由于校核方法的具体实施方式已经在上述实施例中进行了详细的说明,因此,在这里不再赘述。
以上给出了具体的实施方式,但本发明不局限于以上所描述的实施方式。本发明的基本思路在于上述基本方案,对本领域普通技术人员而言,根据本发明的教导,设计出各种变形的模型、公式、参数并不需要花费创造性劳动。在不脱离本发明的原理和精神的情况下对实施方式进行的变化、修改、替换和变形仍落入本发明的保护范围内。

Claims (8)

1.一种风机主轴制动器强度校核方法,其特征在于,包括如下步骤:
1)建立风机主轴制动器各部件的三维几何模型,所述三维几何模型包括活动钳体、固定钳体、上挡板、下挡板及螺栓,所述三维几何模型不包括制动器的安装底座;
2)将所述三维几何模型进行转换得到有限元模型,并对所述有限元模型进行有限元分析;
3)在所述有限元模型中对螺栓圆柱面施加螺栓预紧力,并在有限元模型中对活动钳体的活塞孔施加液压力,根据所述预紧力、液压力及由液压力产生的作用在上挡板、下挡板内侧面上的摩擦力,计算风机主轴制动器的等效应力;根据风机主轴制动器的等效应力进一步计算风机主轴制动器的安全系数及危险位置,所述安全系数的计算公式为:
Figure FDA0004060016290000011
其中,rm为材料安全系数,取1.1;σeqv为等效应力值,单位为MPa;Re为材料屈服强度,单位为MPa。
2.根据权利要求1所述的风机主轴制动器强度校核方法,其特征在于,在将所述三维几何模型转换为有限元模型时,需将所述三维几何模型中的圆角半径、螺纹孔去除,得到去除圆角半径、螺纹孔的三维几何模型。
3.根据权利要求2所述的风机主轴制动器强度校核方法,其特征在于,在活动钳体、固定钳体、上挡板、下挡板及螺栓之间建立Bonded或Rough接触。
4.根据权利要求3所述的风机主轴制动器强度校核方法,其特征在于,所述活动钳体、固定钳体之间通过螺栓连接,在所述活动钳体顶部端面施加Fixed Support约束,对所述固定钳体两侧圆柱孔施加圆柱约束。
5.一种风机主轴制动器强度校核系统,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现以下步骤:
1)建立风机主轴制动器各部件的三维几何模型,所述三维几何模型包括活动钳体、固定钳体、上挡板、下挡板及螺栓,所述三维几何模型不包括制动器的安装底座;
2)将所述三维几何模型进行转换得到有限元模型,并对所述有限元模型进行有限元分析;
3)在所述有限元模型中对螺栓圆柱面施加螺栓预紧力,并在有限元模型中对活动钳体的活塞孔施加液压力,根据所述预紧力、液压力及由液压力产生的作用在上挡板、下挡板内侧面上的摩擦力,计算风机主轴制动器的等效应力;根据主轴制动器的等效应力进一步计算主轴制动器的安全系数及危险位置,所述安全系数的计算公式为:
Figure FDA0004060016290000021
其中,rm为材料安全系数,取1.1;σeqv为等效应力值,单位为MPa;Re为材料屈服强度,单位为MPa。
6.根据权利要求5所述的风机主轴制动器强度校核系统,其特征在于,在将所述三维几何模型转换为有限元模型时,需将所述三维几何模型中的圆角半径、螺纹孔去除,得到去除圆角半径、螺纹孔的三维几何模型。
7.根据权利要求6所述的风机主轴制动器强度校核系统,其特征在于,在活动钳体、固定钳体、上挡板、下挡板及螺栓之间建立Bonded或Rough接触。
8.根据权利要求7所述的风机主轴制动器强度校核系统,其特征在于,所述活动钳体、固定钳体之间通过螺栓连接,在所述活动钳体顶部端面施加Fixed Support约束,对所述固定钳体两侧圆柱孔施加圆柱约束。
CN201811237654.0A 2018-10-23 2018-10-23 一种风机主轴制动器强度校核方法及系统 Expired - Fee Related CN109492275B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811237654.0A CN109492275B (zh) 2018-10-23 2018-10-23 一种风机主轴制动器强度校核方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811237654.0A CN109492275B (zh) 2018-10-23 2018-10-23 一种风机主轴制动器强度校核方法及系统

Publications (2)

Publication Number Publication Date
CN109492275A CN109492275A (zh) 2019-03-19
CN109492275B true CN109492275B (zh) 2023-04-18

Family

ID=65692525

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811237654.0A Expired - Fee Related CN109492275B (zh) 2018-10-23 2018-10-23 一种风机主轴制动器强度校核方法及系统

Country Status (1)

Country Link
CN (1) CN109492275B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111680363B (zh) * 2020-06-09 2023-04-18 南方天合底盘系统有限公司 一种用于浮动式鼓式制动器行车制动效能因数的计算方法
CN113627046B (zh) * 2021-07-12 2023-08-08 许昌许继风电科技有限公司 一种风机偏航制动盘建模与极限强度校核方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005085233A (ja) * 2003-09-11 2005-03-31 Sharp Corp 3次元cad解析システム、解析方法、解析プログラム及び解析プログラムを記録したコンピュータ読取可能な記録媒体
CN105022876B (zh) * 2015-07-15 2018-12-07 许继集团有限公司 一种风机风轮锁强度校核的计算方法

Also Published As

Publication number Publication date
CN109492275A (zh) 2019-03-19

Similar Documents

Publication Publication Date Title
CN109492275B (zh) 一种风机主轴制动器强度校核方法及系统
CN109299532B (zh) 一种风机主机架与后机架连接螺栓强度校核方法及系统
CN103528779B (zh) 一种电机整体机座振动特性的试验系统及试验方法
Babu et al. Stress analysis on steering knuckle of the automobile steering system
CN104573172A (zh) 风力发电机组中的结构件的疲劳分析方法和疲劳分析装置
CN105740504B (zh) 一种乘用车气瓶安装强度的建模方法、计算方法及装置
CN104077428A (zh) 一种服务产业联盟的远程有限元分析方法
CN103020377A (zh) 一种直驱风力发电机的螺栓强度分析方法
CN106339541A (zh) 大型风力发电机组中塔筒门框的焊缝疲劳强度分析方法
CN110472324A (zh) 一种螺栓选型方法及使用该方法获得的螺栓
CN109726435B (zh) 一种接头钉群载荷计算的方法
CN117390914A (zh) 一种考虑螺纹蠕变的连接系统松动有限元仿真方法
CN113127998A (zh) 考虑实际装配过程的螺栓应力场变化的有限元建模方法
CN103366094A (zh) 一种基于力载荷和力边界条件的确定复合材料结构强度的方法
CN103678753B (zh) 一种分离界面连接强度的有限元计算方法
CN110232218B (zh) 一种风力发电机组叶片锁定装置强度的计算方法及装置
CN114996877A (zh) 一种风力发电机组主轴承滚动体受力的计算方法及装置
CN102542118A (zh) 一种几何模型力学性能分析方法及装置
Piscan et al. Finite element analysis of bolted joint
Anand et al. Process automation of simulation using Toolkit/Tool Command Language (TK/TCL) scripting
CN110955930A (zh) 一种矿用工程车辆轻量化模型获取方法及装置
CN112084683A (zh) 一种输电线路用铝合金材料金具的塑性变形分析方法
CN112528393A (zh) 一种船舶轴系连接法兰动力学分析建模方法
CN110954330A (zh) 一种航空发动机转子紧固界面切向接触阻尼预测方法
CN113486562A (zh) 风力发电机组弹性支撑连接螺栓强度评估方法及装置

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20230418