CN113704901B - 发动机频率响应载荷谱自动处理方法 - Google Patents

发动机频率响应载荷谱自动处理方法 Download PDF

Info

Publication number
CN113704901B
CN113704901B CN202110842717.0A CN202110842717A CN113704901B CN 113704901 B CN113704901 B CN 113704901B CN 202110842717 A CN202110842717 A CN 202110842717A CN 113704901 B CN113704901 B CN 113704901B
Authority
CN
China
Prior art keywords
acceleration
curve
engine
load spectrum
frequency
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.)
Active
Application number
CN202110842717.0A
Other languages
English (en)
Other versions
CN113704901A (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.)
Jiangling Motors Corp Ltd
Original Assignee
Jiangling Motors Corp 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 Jiangling Motors Corp Ltd filed Critical Jiangling Motors Corp Ltd
Priority to CN202110842717.0A priority Critical patent/CN113704901B/zh
Publication of CN113704901A publication Critical patent/CN113704901A/zh
Application granted granted Critical
Publication of CN113704901B publication Critical patent/CN113704901B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G01H1/12Measuring characteristics of vibrations in solids by using direct conduction to the detector of longitudinal or not specified vibrations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G01H1/12Measuring characteristics of vibrations in solids by using direct conduction to the detector of longitudinal or not specified vibrations
    • G01H1/14Frequency
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M15/00Testing of engines
    • G01M15/02Details or accessories of testing apparatus
    • 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
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/18Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being spectral information of each sub-band

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Optimization (AREA)
  • Health & Medical Sciences (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Computational Linguistics (AREA)
  • Signal Processing (AREA)
  • Pure & Applied Mathematics (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种发动机频率响应载荷谱自动处理方法,自动处理发动机加速度载荷谱,获得全转速范围内各个频率的最大加速度值;获得全转速下各个频率最大加速度峰值的包络线作为振动响应和声学响应分析载荷谱;使用一个加速度谱在发动机全转速范围内考核零部件振动响应和声学响应;运用该自动化处理载荷谱技术能够自动处理全转速下发动机加速度载荷谱,保证载荷谱处理数据完整性,降低了工程师的劳动强度,提升工作效力;能够获得零部件在发动机全转速范围内的振动响应和声学响应结果,节省分析资源,为零部件NVH优化提供更可靠的保障。

Description

发动机频率响应载荷谱自动处理方法
技术领域
本发明涉及汽车发动机噪音控制方法,具体的说,是涉及一种发动机频率响应载荷谱自动处理方法。
背景技术
发动机振动噪声问题对整车NVH性能影响较大,有效控制发动机振动噪声成为了汽车NVH开发过程中不可缺少的一部分。发动机主要结构件及薄壁件对发动机NVH性能起决定性作用。现有发动机结构件NVH以控制模态为主,同时分析零部件振动响应,响应载荷谱通常有两种,一是额定功率点的测试加速度载荷,二是发动机最高转速的测试加速度载荷。这两种方式都只能体现发动机在特定转速下的振动激励,不能代表发动机全转速范围内的振动激励,以这两种振动载荷作为零件相应分析的输入无法反映零部件在全转速范围发动机振动激励下的NVH特性。
发明内容
针对上述现有技术中的不足,本发明提供一种能够综合考虑发动机全转速范围振动激励对零部件的影响,为发动机零部件NVH设计及优化提供可靠保障的发动机频率响应载荷谱自动处理方法。
本发明所采取的技术方案是:
一种发动机频率响应载荷谱自动处理方法,包括如下步骤:
步骤S101,NVH测试设备测得发动机零部件激励点的加速度载荷,获得发动机不同转速下激励点的时域加速度值;
步骤S102,通过傅里叶变换将时域加速度值转换成频域加速度数据;
步骤S103,将频域加速度数据导出为txt格式的文件;使用编写的程序读取频域加速度数据;
将频域加速度数据导出为txt格式的文件;使用编写的MATLAB程序读取频域加速度数据,每一个转速分别读入频率及对应的加速度值,如读取1000转的加速度数据,软件自动找到txt数据中1000转的数据区域,从0开始间隔1Hz读取频率值直到4000Hz为止并保存,后读取每个频率对应的加速度值并保存为1000转的频谱曲线,按此顺序可读取不同转速的频谱曲线。
步骤S104,以怠速工况加速度频谱曲线作为初始曲线,记录不同频率下的加速度值,该曲线定义为L1;
如步骤S103所述,起始转速为怠速转速,保存下来的曲线L1是过程数据,不会输出的。
步骤S105,用L1曲线的加速度值与其他转速的加速度值在频域范围内进行两两对比,加速度取大值保存下来形成频域加速度曲线,该曲线定义为L2;
每个转速的频率范围都是1-4000Hz,每个频率都会有一个加速度值,通过对比取大值软件后台会自动保存,继续与下一个转速相同频率下的加速度值进行对比。
步骤S106,判断发动机所有转速加速度频域数据参与对比是否完成,完成跳转步骤S107;未完成跳转步骤S105;
步骤S107,经过n轮加速度对比,直至发动机所有转速加速度频域数据参与对比计算,形成加速度频域曲线Ln;
步骤S108,对形成的加速度载荷曲线Ln峰值连线形成载荷谱的包络线,该包络线能够完全覆盖发动机全转速下不同频率的加速度峰值,该加速度载荷包络线命名为Lmax;加速度载荷谱曲线Lmax为目标加速度载荷谱曲线;
Ln曲线频率1Hz的加速度值为A1,依次规律频率为nHz的加速度为An,4000Hz对应的加速度为A4000,频率按从低到高的顺序相邻频率的加速度值进行两两对比,若An>An-1且An>An+1,则该点的加速度及频率值输出作为Lmax曲线上的一个点记录为P1,按此程序第二个输出点记录为P2,依次可获得n个点输出,这些数据作为Lmax曲线上的点形成Lmax曲线,及目标加速度载荷曲线。
步骤S109,导出加速度载荷谱曲线Lmax,加速度载荷谱曲线Lmax作为发动机零部件振声响应的加速度载荷谱。
步骤S110,加速度载荷谱曲线Lmax进行零部件振动声学响应分析。
优选的,所述加速度频域曲线Ln为发动机全转速范围内各频率加速度最大值曲线。
优选的,所述加速度载荷谱曲线Lmax进行零部件振动声学响应分析的步骤如下:
步骤S301,建立需要分析零件的有限元网格模型;
步骤S302,按实际情况建立零部件之间的约束关系;
步骤S303,不同零件按实际情况给定不同的属性参数;
步骤S304,零部件安装孔位置进行全约束;
步骤S305,约束点加载Lmax载荷谱进行振动声学响应分析;
步骤S306,输出零部件的振动声学响应结果,对结果进行分析,与零部件振动响应分析数据库进行对比,确实是否满足NVH要求,不满足则需要根据分析结果对零部件进行结构优化重新计算。
本发明相对现有技术的有益效果:
本发明发动机频率响应载荷谱自动处理方法,使用一个载荷谱综合考虑发动机零部件在全转速振动激励下的振动响应,为发动机零部件NVH优化提供准确载荷。
本发明发动机频率响应载荷谱自动处理方法,克服了现有技术存在发动机零部件响应分析中载荷考虑不完善,分析效率低等问题,提供一种全新的载荷谱自动处理技术,运用该自动化处理载荷谱技术能够自动处理全转速下发动机加速度载荷谱,保证载荷谱处理数据完整性,降低了工程师的劳动强度,提升工作效力并完善分析结果。
本发明发动机频率响应载荷谱自动处理方法,自动处理发动机加速度载荷谱,获得全转速范围内各个频率的最大加速度值;获得全转速下各个频率最大加速度峰值的包络线作为振动响应和声学响应分析载荷谱;使用一个加速度谱在发动机全转速范围内考核零部件振动响应和声学响应;获得零部件在发动机全转速范围内的振动响应和声学响应结果,节省分析资源,为零部件NVH优化提供更可靠的保障。
附图说明
图1是发动机频率响应载荷谱自动处理方法的流程图;
图2是发动机频率响应载荷谱自动处理方法的处理后的加速度载荷谱对比图(一);
图3是发动机频率响应载荷谱自动处理方法的处理后的加速度载荷谱对比图(二);
图4是发动机频率响应载荷谱自动处理方法的薄壁件振动响应结果对比图;
图5是发动机频率响应载荷谱自动处理方法的薄壁件声学响应结果对比图。
具体实施方式
以下参照附图及实施例对本发明进行详细的说明:
图2是频率范围0-1000赫兹的加速度载荷谱图,图3是频率范围0-2000赫兹的加速度载荷谱图。
附图1-5可知,一种发动机频率响应载荷谱自动处理方法,包括如下步骤:
步骤S101,NVH测试设备测得发动机零部件激励点的加速度载荷,获得发动机不同转速下激励点的时域加速度值;
步骤S102,通过傅里叶变换将时域加速度值转换成频域加速度数据;
步骤S103,将频域加速度数据导出为txt格式的文件;使用编写的程序读取频域加速度数据;
将频域加速度数据导出为txt格式的文件;使用编写的MATLAB程序读取频域加速度数据,每一个转速分别读入频率及对应的加速度值,如读取1000转的加速度数据,软件自动找到txt数据中1000转的数据区域,从0开始间隔1Hz读取频率值直到4000Hz为止并保存,后读取每个频率对应的加速度值并保存为1000转的频谱曲线,按此顺序可读取不同转速的频谱曲线。
步骤S104,以怠速工况加速度频谱曲线作为初始曲线,记录不同频率下的加速度值,该曲线定义为L1;
如步骤S103所述,起始转速为怠速转速,保存下来的曲线L1是过程数据,不会输出的。
步骤S105,用L1曲线的加速度值与其他转速的加速度值在频域范围内进行两两对比,加速度取大值保存下来形成频域加速度曲线,该曲线定义为L2;
每个转速的频率范围都是1-4000Hz,每个频率都会有一个加速度值,通过对比取大值软件后台会自动保存,继续与下一个转速相同频率下的加速度值进行对比。
步骤S106,判断发动机所有转速加速度频域数据参与对比是否完成,完成跳转步骤S107;未完成跳转步骤S105;
步骤S107,经过n轮加速度对比,直至发动机所有转速加速度频域数据参与对比计算,形成加速度频域曲线Ln;
步骤S108,对形成的加速度载荷曲线Ln峰值连线形成载荷谱的包络线,该包络线能够完全覆盖发动机全转速下不同频率的加速度峰值,该加速度载荷包络线命名为Lmax;加速度载荷谱曲线Lmax为目标加速度载荷谱曲线;
Ln曲线频率1Hz的加速度值为A1,依次规律频率为nHz的加速度为An,4000Hz对应的加速度为A4000,频率按从低到高的顺序相邻频率的加速度值进行两两对比,若An>An-1且An>An+1,则该点的加速度及频率值输出作为Lmax曲线上的一个点记录为P1,按此程序第二个输出点记录为P2,依次可获得n个点输出,这些数据作为Lmax曲线上的点形成Lmax曲线,及目标加速度载荷曲线。
步骤S109,导出加速度载荷谱曲线Lmax,加速度载荷谱曲线Lmax作为发动机零部件振声响应的加速度载荷谱。
步骤S110,加速度载荷谱曲线Lmax进行零部件振动声学响应分析。
所述加速度频域曲线Ln为发动机全转速范围内各频率加速度最大值曲线。
所述加速度载荷谱曲线Lmax进行零部件振动声学响应分析的步骤如下:
步骤S301,建立需要分析零件的有限元网格模型;
步骤S302,按实际情况建立零部件之间的约束关系;
步骤S303,不同零件按实际情况给定不同的属性参数;
步骤S304,零部件安装孔位置进行全约束;
步骤S305,约束点加载Lmax载荷谱进行振动声学响应分析;
步骤S306,输出零部件的振动声学响应结果,对结果进行分析,与零部件振动响应分析数据库进行对比,确实是否满足NVH要求,不满足则需要根据分析结果对零部件进行结构优化重新计算。
本发明发动机频率响应载荷谱自动处理方法,使用一个载荷谱综合考虑发动机零部件在全转速振动激励下的振动响应,为发动机零部件NVH优化提供准确载荷。
本发明发动机频率响应载荷谱自动处理方法,克服了现有技术存在发动机零部件响应分析中载荷考虑不完善,分析效率低等问题,提供一种全新的载荷谱自动处理技术,运用该自动化处理载荷谱技术能够自动处理全转速下发动机加速度载荷谱,保证载荷谱处理数据完整性,降低了工程师的劳动强度,提升工作效力并完善分析结果。
现有技术进行发动机零部件振动声学响应分析,大多以发动机功率最大点的转速加速度频谱曲线或者是最大转速的加速度频谱曲线进行响应分析,但单一转速的载荷谱不能反应发动机全部工作转速的振动声学响应,若需要计算发动机所有工作转速下的响应结果则需要以各个转速的加速度载荷谱进行多轮的响应分析,这个劳动强度是较大的,除了计算建模需要花很多时间外,后续的结果处理也需要较多的工作时间,所以通常情况下只用额定转速或者最大转速下的加速度载荷谱进行计算分析。本发明通过处理发动机不同工作转速下加速度载荷谱,获得全转速范围内各频率加速度的最大值,形成Lmax载荷谱,可以减少分析轮次,降低分析及结果后处理时间。也方便形成数据库,对后续分析结果评判提供有效数据。如图4是额定功率点转速下加速度载荷谱及用Lmax载荷谱计算得到的零部件振动响应分析结果,用Lmax载荷谱计算的响应结果比额定转速载荷谱计算的结果大,能够反应零部件在发动机工作转速下的振动水平,而不只是额定转速的响应情况,为零部件NVH优化提供准确的方向。
本发明发动机频率响应载荷谱自动处理方法,自动处理发动机加速度载荷谱,获得全转速范围内各个频率的最大加速度值;获得全转速下各个频率最大加速度峰值的包络线作为振动响应和声学响应分析载荷谱;使用一个加速度谱在发动机全转速范围内考核零部件振动响应和声学响应;获得零部件在发动机全转速范围内的振动响应和声学响应结果,节省分析资源,为零部件NVH优化提供更可靠的保障。
以上所述,仅是本发明的较佳实施例而已,并非对本发明的结构作任何形式上的限制。凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均属于本发明的技术方案范围内。

Claims (2)

1.一种发动机频率响应载荷谱自动处理方法,其特征在于,
步骤S101,NVH测试设备测得发动机零部件激励点的加速度载荷,获得发动机不同转速下激励点的时域加速度值;
步骤S102,通过傅里叶变换将时域加速度值转换成频域加速度数据;
步骤S103,将频域加速度数据导出为txt格式的文件;读取频域加速度数据;
步骤S104,以怠速工况加速度频谱曲线作为初始曲线,记录不同频率下的加速度值,该曲线定义为L1;
步骤S105,用L1曲线的加速度值与其他转速的加速度值在频域范围内进行两两对比,加速度取大值保存下来形成频域加速度曲线,该曲线定义为L2;
步骤S106,判断发动机所有转速加速度频域数据参与对比是否完成,完成跳转步骤S107;未完成跳转步骤S105;
步骤S107,经过n轮加速度对比,直至发动机所有转速加速度频域数据参与对比计算,形成加速度频域曲线Ln;
步骤S108,对形成的加速度载荷曲线Ln峰值连线形成载荷谱的包络线,该包络线能够完全覆盖发动机全转速下不同频率的加速度峰值,该加速度载荷包络线命名为Lmax;加速度载荷谱曲线Lmax为目标加速度载荷谱曲线;
步骤S109,导出加速度载荷谱曲线Lmax,加速度载荷谱曲线Lmax作为发动机零部件振声响应的加速度载荷谱;
步骤S110,加速度载荷谱曲线Lmax进行零部件振动声学响应分析;
所述加速度载荷谱曲线Lmax进行零部件振动声学响应分析的步骤如下:
步骤S301,建立需要分析零件的有限元网格模型;
步骤S302,按实际情况建立零部件之间的约束关系;
步骤S303,不同零件按实际情况给定不同的属性参数;
步骤S304,零部件安装孔位置进行全约束;
步骤S305,约束点加载Lmax载荷谱进行振动声学响应分析;
步骤S306,输出零部件的振动声学响应结果,对结果进行分析,与零部件振动响应分析数据库进行对比,确实是否满足NVH要求,不满足则需要根据分析结果对零部件进行结构优化重新计算。
2.根据权利要求1所述发动机频率响应载荷谱自动处理方法,其特征在于:所述加速度频域曲线Ln为发动机全转速范围内各频率加速度最大值曲线。
CN202110842717.0A 2021-07-26 2021-07-26 发动机频率响应载荷谱自动处理方法 Active CN113704901B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110842717.0A CN113704901B (zh) 2021-07-26 2021-07-26 发动机频率响应载荷谱自动处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110842717.0A CN113704901B (zh) 2021-07-26 2021-07-26 发动机频率响应载荷谱自动处理方法

Publications (2)

Publication Number Publication Date
CN113704901A CN113704901A (zh) 2021-11-26
CN113704901B true CN113704901B (zh) 2023-10-20

Family

ID=78650462

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110842717.0A Active CN113704901B (zh) 2021-07-26 2021-07-26 发动机频率响应载荷谱自动处理方法

Country Status (1)

Country Link
CN (1) CN113704901B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6389887B1 (en) * 1999-12-23 2002-05-21 Snecma Moteurs Process for the detection of damage to components of an engine
CN102023074A (zh) * 2010-11-10 2011-04-20 中国第一汽车集团公司 汽车发动机机油盘振动噪声性能分析方法
CN104834763A (zh) * 2014-10-31 2015-08-12 北汽福田汽车股份有限公司 获取发动机辐射噪声的方法及发动机设计优化方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6389887B1 (en) * 1999-12-23 2002-05-21 Snecma Moteurs Process for the detection of damage to components of an engine
CN102023074A (zh) * 2010-11-10 2011-04-20 中国第一汽车集团公司 汽车发动机机油盘振动噪声性能分析方法
CN104834763A (zh) * 2014-10-31 2015-08-12 北汽福田汽车股份有限公司 获取发动机辐射噪声的方法及发动机设计优化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
曾利奎 ; 杨平 ; .电控自动变速器载荷谱的采集及分析.现代制造工程.2017,(第04期),第61-65页. *

Also Published As

Publication number Publication date
CN113704901A (zh) 2021-11-26

Similar Documents

Publication Publication Date Title
US7505949B2 (en) Process model error correction method and system
CN111766080B (zh) 汽车工况筛选方法、发动机标定方法及系统和可读存储介质
JP2010530179A (ja) 仮想センサ・システムおよび方法
CN101532915A (zh) 汽车尾气排放及车载自诊断系统的标定方法
Jiang et al. Implementation of model-based calibration for a gasoline engine
CN113704901B (zh) 发动机频率响应载荷谱自动处理方法
CN116341312A (zh) 一种发动机分体式动力学振动响应结果评价方法
CN116220870A (zh) 车辆尾气排放的控制方法、装置、电子设备和存储介质
CN113094913A (zh) 排放预估方法和装置
CN114625107B (zh) 一种基于stars测控系统的ecu自动标定方法及系统
CN116577104A (zh) 一种发动机故障检测方法及设备
CN113378304B (zh) 一种车辆性能目标的确定方法、装置、存储介质及设备
CN113408560B (zh) 发动机测试数据分类方法、电子设备和可读存储介质
CN112084656A (zh) 一种叶片振动疲劳概率寿命预测系统及预测方法
CN110943485A (zh) 双馈风电场等值模型仿真可信度的指标评价方法
CN115638900B (zh) 一种排气管温度的确定方法、系统、存储介质及电子设备
CN114757117B (zh) 柴油机进气道性能评估方法、系统、终端设备及存储介质
CN114091540B (zh) 柴油机冷试智能检测模型构建方法、检测方法及系统
CN112729849B (zh) 用于测试至少一个驱动系构件的方法和系统
CN114818259A (zh) 一种验证阶段rde水平预测模型建立方法、存储介质、电子设备以及预测方法
US20240176849A1 (en) Method for Characterizing a Test Bench and Associated Method for Testing and Producing a Component
CN117688294A (zh) 车辆起动异常诊断方法、装置、设备及存储介质
CN115219220A (zh) 一种车辆nvh性能衰减分析方法、装置及存储介质
CN117516949A (zh) 发动机爆震诊断标定方法、装置、电子设备和存储介质
CN116337453A (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