CN113984276B - 检测怠速工况下发动机各气缸压力不均匀性的方法 - Google Patents

检测怠速工况下发动机各气缸压力不均匀性的方法 Download PDF

Info

Publication number
CN113984276B
CN113984276B CN202111232152.0A CN202111232152A CN113984276B CN 113984276 B CN113984276 B CN 113984276B CN 202111232152 A CN202111232152 A CN 202111232152A CN 113984276 B CN113984276 B CN 113984276B
Authority
CN
China
Prior art keywords
cylinder
engine
pressure
cycle
value
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
CN202111232152.0A
Other languages
English (en)
Other versions
CN113984276A (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.)
Chongqing Changan Automobile Co Ltd
Original Assignee
Chongqing Changan Automobile 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 Chongqing Changan Automobile Co Ltd filed Critical Chongqing Changan Automobile Co Ltd
Priority to CN202111232152.0A priority Critical patent/CN113984276B/zh
Publication of CN113984276A publication Critical patent/CN113984276A/zh
Application granted granted Critical
Publication of CN113984276B publication Critical patent/CN113984276B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L23/00Devices or apparatus for measuring or indicating or recording rapid changes, such as oscillations, in the pressure of steam, gas, or liquid; Indicators for determining work or energy of steam, internal-combustion, or other fluid-pressure engines from the condition of the working fluid
    • G01L23/08Devices or apparatus for measuring or indicating or recording rapid changes, such as oscillations, in the pressure of steam, gas, or liquid; Indicators for determining work or energy of steam, internal-combustion, or other fluid-pressure engines from the condition of the working fluid operated electrically
    • 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/04Testing internal-combustion engines
    • G01M15/08Testing internal-combustion engines by monitoring pressure in cylinders
    • 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
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)

Abstract

本发明公开的一种检测怠速工况下发动机各气缸压力不均匀性的方法,包括以下步骤:在待测发动机气缸内布置缸压传感器,并连接发动机曲轴转角传感器信号;采集发动机在稳定怠速工况下传感器数据信息;根据所测试燃烧压力与曲轴转角数据,计算出发动机各气缸每个工作循环下的平均指示有效压力;计算出每一循环下各气缸有效压力的最小值与平均值的比值,并将该值定义为抖振参数(LNV);根据计算得到的每一循环下抖振参数来绘制出该值随循环数的变化趋势图;根据Chebyshen不等式计算气缸压力在一定范围内波动的概率,并将该概率值作为抖振参数的检测标准;根据检测标准来客观衡量怠速发动机各气缸压力不均匀性。

Description

检测怠速工况下发动机各气缸压力不均匀性的方法
技术领域
本发明实施例涉及发动机技术领域,尤其涉及一种检测怠速工况下发动机各气缸压力不均匀性的方法。
背景技术
怠速工况作为汽车发动机最典型也是最重要的运行工况之一,要消耗约25%的运行时间和30%的燃料,在该工况下发动机气缸内混合气体燃烧产生稳定压力对汽车动力性、经济性以及排放要求都十分重要;同时,气缸燃烧过程中气缸压不均匀性,也使得振动及噪声增大,直接影响对汽车NVH性能的评价。因此,客观检测怠速发动机气缸压力不均匀性对优化怠速控制策略,改善发动机怠速稳定性能,实现发动机节能减排具有重要意义。发动机气缸压不均匀性的产生与气缸内混合气体不规则燃烧过程密切联系,而气缸不规则燃烧主要包括各循环之间的燃烧变动和各气缸之间的燃烧差异,前者称为循环波动,后者称为各气缸工作不均匀。对于多气缸汽油发动机,除了单气缸燃烧压力循环波动的影响外,各气缸工作时燃烧压力的不均匀也是影响发动机性能的主要原因。
目前普遍适用的衡量发动机气缸循环波动的参数是循环变动系数(Cov),该方法是通过采集发动机工作时的各循环最高燃烧压力pmax,并计算最高燃烧压力pmax的循环波动率δp来评价循环波动程度,其定义如下:
式中,σp为pmax的标准方差,为pmax的平均值。但是这种方法只能通过单一数值检测单个气缸各循环间的燃烧压力波动情况,而不能检测各个气缸之间每个循环下的气缸压不均匀性,这样无法有效衡量各气缸之间的工作不均匀,使得动力性、经济性、排放性以及噪声与振动等整机指标难以优化。
因此,有必要开发一种检测怠速工况下发动机各气缸压力不均匀性的方法。
发明内容
本发明提供一种检测怠速工况下发动机各气缸压力不均匀性的方法,能客观衡量发动机在怠速工作过程中各气缸工作不均匀性。
本发明所述的一种检测怠速工况下发动机各气缸压力不均匀性的方法,包括以下步骤:
(1)布置测试系统:在待测发动机的每个气缸内布置缸压传感器,各缸压传感器经电荷放大器后连接至数据采集前端模块;同时,将曲轴转角传感器连接至数据采集前端模块;
(2)怠速工况数据采集:将发动机点火,并充分热机至怠速转速趋于稳定,通过数据采集前端采集N个工作循环数下的缸压传感器测试数据;同时收集曲轴转角传感器的测试数据;
(3)计算单气缸一个工作循环下的平均有效指示压力IMEP;
(4)计算抖振参数LNV:计算每个循环下各气缸有效压力的最小值与平均值的比值,即LNVj值;
(5)抖振参数LNV的检测标准:计算每个气缸在N个循环内的方差,计算各气缸压力控制在ε范围内波动的概率P,将计算出的各气缸压力在一定范围ε内波动的概率值作为各循环下抖振参数的检测标准;
(6)发动机各气缸压力不均匀性的检测:根据计算得到的各循环下LNVj值,绘制出该LNVj值随循环数的变化趋势图,并标注检测标准;若检测出各循环LNVj值均未超出检测标准,则表明发动机在工作过程中各气缸压力波动均控制在一定范围内,各气缸工作不均匀性较小;若检测出某些循环的LNVj值超出检测标准,则表明发动机在这些循环下各气缸压力波动较大,各气缸工作不均匀性较大。
可选地,计算单气缸一个工作循环下的平均有效指示压力,具体为:
其中,Pj (k)表示第k个气缸在第j个循环的IMEP值,k=1,2,…,m,m为发动机的气缸总数,j=1,2,…,N,N为循环总次数;D为发动机活塞有效直径;Vs为单个气缸的排量;n表示单个循环的总采样点数;pi表示在第j个循环下第i个采样点所采集的缸压传感器测试数据;r为发动机曲轴半径;θi为第i个采样点时的发动机曲轴转角;θi-1是第i-1个采样点时的发动机曲轴转角;λ为发动机曲轴半径与活塞连杆长度之比。
可选地,计算LNVj的公式如下:
其中,Min(P′j)为第j个循环下m个气缸内IMEP的最小值;Mean(P′j)为第j个循环下m个气缸内IMEP的平均值,即
可选地,计算每个气缸在N个循环内的方差,具体为:
其中:D(P)k表示第k个气缸在N个循环内的方差,为第k个气缸在N个循环内IMEP的平均值,即/>
可选地,根据Chebyshev不等式,计算各气缸压力控制在ε范围内波动的概率P:
其中,为各气缸计算得到的IMEP方差平均值,即/>
本发明具有以下优点:通过检测各个气缸每一循环下气缸压力波动的大小,来客观衡量发动机在怠速工作过程中各气缸工作不均匀性,该方法为优化发动机怠速控制策略提供依据。
附图说明
图1为发动机气缸压力及曲轴转角测试系统示意图;
图2为抖振参数LNV随循环数的变化趋势图之一;
图3为抖振参数LNV随循环数的变化趋势图之二;
图4为本实施例的流程图;
图5为本实施例中怠速工况气缸压不均匀性检测示意图;
图中:1、缸压传感器,2、电荷放大器,3、数据采集前端模块,4、曲轴转角传感器。
具体实施方式
下面结合附图对本发明作进一步说明。
如图4所示,本实施例中,一种检测怠速工况下发动机各气缸压力不均匀性的方法,在待测发动机气缸内布置缸压传感器1,并连接发动机曲轴转角传感器4信号;采集发动机在稳定怠速工况下传感器数据信息;根据所测试燃烧压力与曲轴转角数据,计算出发动机各气缸每个工作循环下的平均指示有效压力;计算出每一循环下各气缸有效压力的最小值与平均值的比值,并将该值定义为抖振参数(LNV);根据计算得到的每一循环下抖振参数来绘制出该值随循环数的变化趋势图;根据Chebyshen不等式计算气缸压力在一定范围内波动的概率,并将该概率值作为抖振参数的检测标准;根据检测标准来客观衡量怠速发动机各气缸压力不均匀性。
如图4和图5所示,本实施例中,一种检测怠速工况下发动机各气缸压力不均匀性的方法,包括以下步骤:
(1)布置测试系统:在待测发动机的每个气缸内布置缸压传感器1,各缸压传感器1经电荷放大器2后连接至数据采集前端模块3;同时,将曲轴转角传感器4连接至数据采集前端模块3,参见图1。
(2)怠速工况数据采集:将发动机点火,并充分热机至怠速转速趋于稳定,通过数据采集前端采集N个工作循环数下的缸压传感器1测试数据;同时收集曲轴转角传感器4的测试数据。本实施例中,N一般取值为50~200次,若气缸压波动较大时应进一步加大采样循环数;采样频率一般为曲轴转角每转过1°采集一次燃烧压力。
(3)计算平均指示有效压力IMEP:平均指示有效压力Pimep是指单位气缸工作容积所作的循环指示有效功,也是一个作用于活塞上的假想平均压力,当其作用于活塞上使其运动一个冲程所作之功,正好为循环指示功W,因此:
Pimep=W/Vs (1);
其中:W为单个气缸循环有效功,单位为KJ;Vs为单气缸排量(同一个发动机的各气缸的排量是一样的),单位为L。当缸内系统经历微小的准静态过程使活塞移动一微小的距离Δs时,气体所做的功为dW=FΔs=PiAds=PidV,其中A为活塞等效面积(即D为活塞的有效直径);如果缸内的体积经过一个准静态过程由V1变化为V2,则在该过程中,气体所做的功为/>
根据定义所得的平均指示有效压力的表达式为:
根据活塞-曲轴-连杆结构运动学,活塞位移Δs的近似公式为: 当活塞移动微小的距离Δs,则平均有效指示压力表示为:
其中,Pimep是指单位气缸工作容积所作的循环指示有效功,D为发动机活塞有效直径,r为发动机曲轴半径,l为发动机活塞连杆长度,λ为发动机曲轴半径与活塞连杆长度之比,θi为第i个采样点时的发动机曲轴转角,pi为该转角下所测得的燃烧压力。
计算单气缸一个工作循环下的平均有效指示压力为:
其中,Pj (k)表示第k个气缸在第j个循环的IMEP值,k=1,2,…,m,m为发动机的气缸总数,j=1,2,…,N,N为循环总次数;D为发动机活塞有效直径;Vs为单个气缸的排量;n表示单个循环的总采样点数;pi表示在第j个循环下第i个采样点所采集的缸压传感器测试数据;r为发动机曲轴半径;θi为第i个采样点时的发动机曲轴转角;θi-1是第i-1个采样点时的发动机曲轴转角;λ为发动机曲轴半径与活塞连杆长度之比。
本实施例中,某一发动机各气缸每个循环下的平均有效指示压力表:
(4)计算抖振参数LNV:计算每个循环下各气缸有效压力的最小值与平均值的比值,即LNVj值;
其中,Min(P′j)为第j个循环下m个气缸内IMEP的最小值,比如:当m为4,j为1,Min(P′1)即为和/>中的最小值;Mean(P′j)为第j个循环下m个气缸内IMEP的平均值,即/>
(5)抖振参数LNV的检测标准:计算每个气缸在N个循环内的方差:
其中:D(P)k表示第k个气缸在N个循环内的方差,为第k个气缸在N个循环内IMEP的平均值,即/>
根据Chebyshev不等式,计算各气缸压力控制在一定范围ε内波动的概率P:
其中,表示计算得出的有效压力Pj (k)偏离平均值/>的程度小于ε的概率。/>为各气缸计算得到的IMEP方差平均值,即/>一般地将ε控制在/>以内时,/>
将计算出的各气缸压力在一定范围ε内波动的概率值作为各循环下抖振参数的检测标准。
(6)发动机各气缸压力不均匀性的检测:根据计算得到的各循环下LNVj值,绘制出该LNVj值随循环数的变化趋势图,并标注检测标准;若检测出各循环LNVj值均未超出检测标准(即),则表明发动机在工作过程中各气缸压力波动均控制在一定范围内,各气缸工作不均匀性较小,参见图2。若检测出某些循环的LNVj值超出检测标准,则表明发动机在这些循环下各气缸压力波动较大,各气缸工作不均匀性较大,参见图3。通过本方法就能够客观检测出各气缸在每个循环下的气缸压不均匀情况,而当检测出发动机各气缸工作压力不均匀性较大时,就需要对发动机怠速控制策略进行优化,以提升发动机工作性能。
本方法对多气缸发动机各气缸工作不均匀性均可检测,同时针对不同发动机的性能要求,可通过调整ε值来修正检测标准以满足不同性能需求。

Claims (4)

1.一种检测怠速工况下发动机各气缸压力不均匀性的方法,其特征在于,包括以下步骤:
(1)布置测试系统:在待测发动机的每个气缸内布置缸压传感器,各缸压传感器经电荷放大器后连接至数据采集前端模块;同时,将曲轴转角传感器连接至数据采集前端模块;
(2)怠速工况数据采集:将发动机点火,并充分热机至怠速转速趋于稳定,通过数据采集前端采集N个工作循环数下的缸压传感器测试数据;同时收集曲轴转角传感器的测试数据;
(3)计算单气缸一个工作循环下的平均有效指示压力IMEP,具体为:
其中,Pj (k)表示第k个气缸在第j个循环的IMEP值,k=1,2,…,m,m为发动机的气缸总数,j=1,2,…,N,N为循环总次数;D为发动机活塞有效直径;Vs为单个气缸的排量;n表示单个循环的总采样点数;pi表示在第j个循环下第i个采样点所采集的缸压传感器测试数据;r为发动机曲轴半径;θi为第i个采样点时的发动机曲轴转角;θi-1是第i-1个采样点时的发动机曲轴转角;λ为发动机曲轴半径与活塞连杆长度之比;
(4)计算抖振参数LNV:计算每个循环下各气缸有效压力的最小值与平均值的比值,即LNVj值:
(5)抖振参数LNV的检测标准:计算每个气缸在N个循环内的方差,计算各气缸压力控制在ε范围内波动的概率P,将计算出的各气缸压力在一定范围ε内波动的概率值作为各循环下抖振参数的检测标准;
(6)发动机各气缸压力不均匀性的检测:根据计算得到的各循环下LNVj值,绘制出该LNVj值随循环数的变化趋势图,并标注检测标准;若检测出各循环LNVj值均未超出检测标准,则表明发动机在工作过程中各气缸压力波动均控制在一定范围内,各气缸工作不均匀性较小;若检测出某些循环的LNVj值超出检测标准,则表明发动机在这些循环下各气缸压力波动较大,各气缸工作不均匀性较大。
2.根据权利要求1所述的检测怠速工况下发动机各气缸压力不均匀性的方法,其特征在于:计算LNVj的公式如下:
其中,Min(P′j)为第j个循环下m个气缸内IMEP的最小值;Mean(P′j)为第j个循环下m个气缸内IMEP的平均值,即
3.根据权利要求2所述的检测怠速工况下发动机各气缸压力不均匀性的方法,其特征在于:计算每个气缸在N个循环内的方差,具体为:
其中:D(P)k表示第k个气缸在N个循环内的方差,为第k个气缸在N个循环内IMEP的平均值,即/>
4.根据权利要求3所述的检测怠速工况下发动机各气缸压力不均匀性的方法,其特征在于:
根据Chebyshev不等式,计算各气缸压力控制在ε范围内波动的概率P:
其中,为各气缸计算得到的IMEP方差平均值,即/>
CN202111232152.0A 2021-10-22 2021-10-22 检测怠速工况下发动机各气缸压力不均匀性的方法 Active CN113984276B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111232152.0A CN113984276B (zh) 2021-10-22 2021-10-22 检测怠速工况下发动机各气缸压力不均匀性的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111232152.0A CN113984276B (zh) 2021-10-22 2021-10-22 检测怠速工况下发动机各气缸压力不均匀性的方法

Publications (2)

Publication Number Publication Date
CN113984276A CN113984276A (zh) 2022-01-28
CN113984276B true CN113984276B (zh) 2024-01-12

Family

ID=79740290

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111232152.0A Active CN113984276B (zh) 2021-10-22 2021-10-22 检测怠速工况下发动机各气缸压力不均匀性的方法

Country Status (1)

Country Link
CN (1) CN113984276B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4525781A (en) * 1980-10-01 1985-06-25 Toyota Jidosha Kogyo Kabushiki Kaisha Method of determining fluctuations in indicated mean effective pressure of engine and apparatus therefor
JPH07180598A (ja) * 1993-12-21 1995-07-18 Mitsubishi Motors Corp 内燃機関の燃焼状態判定方法および燃焼状態制御装置
JPH11183330A (ja) * 1997-12-19 1999-07-09 Honda Motor Co Ltd 内燃エンジンの失火判別方法
FR2862711A1 (fr) * 2003-11-24 2005-05-27 Peugeot Citroen Automobiles Sa Systeme de calibrage d'une chaine d'acquisition de la pression dans un cylindre de moteur diesel de vehicule automobile
JP2007211722A (ja) * 2006-02-10 2007-08-23 Toyota Motor Corp 内燃機関の制御装置
WO2010050544A1 (ja) * 2008-10-30 2010-05-06 日立オートモティブシステムズ株式会社 筒内圧測定装置
CN102808694A (zh) * 2011-05-31 2012-12-05 通用汽车环球科技运作有限责任公司 用于估计发动机气缸的平均指示有效压力的系统和方法
CN105571871A (zh) * 2016-01-04 2016-05-11 哈尔滨工程大学 一种在线诊断柴油机工作不均匀性的方法
CN106150719A (zh) * 2015-04-28 2016-11-23 长城汽车股份有限公司 发动机的控制方法、系统及车辆
CN106687676A (zh) * 2014-09-24 2017-05-17 瓦锡兰芬兰有限公司 用于起动双燃料发动机的方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4525781A (en) * 1980-10-01 1985-06-25 Toyota Jidosha Kogyo Kabushiki Kaisha Method of determining fluctuations in indicated mean effective pressure of engine and apparatus therefor
JPH07180598A (ja) * 1993-12-21 1995-07-18 Mitsubishi Motors Corp 内燃機関の燃焼状態判定方法および燃焼状態制御装置
JPH11183330A (ja) * 1997-12-19 1999-07-09 Honda Motor Co Ltd 内燃エンジンの失火判別方法
FR2862711A1 (fr) * 2003-11-24 2005-05-27 Peugeot Citroen Automobiles Sa Systeme de calibrage d'une chaine d'acquisition de la pression dans un cylindre de moteur diesel de vehicule automobile
JP2007211722A (ja) * 2006-02-10 2007-08-23 Toyota Motor Corp 内燃機関の制御装置
WO2010050544A1 (ja) * 2008-10-30 2010-05-06 日立オートモティブシステムズ株式会社 筒内圧測定装置
CN102808694A (zh) * 2011-05-31 2012-12-05 通用汽车环球科技运作有限责任公司 用于估计发动机气缸的平均指示有效压力的系统和方法
CN106687676A (zh) * 2014-09-24 2017-05-17 瓦锡兰芬兰有限公司 用于起动双燃料发动机的方法
CN106150719A (zh) * 2015-04-28 2016-11-23 长城汽车股份有限公司 发动机的控制方法、系统及车辆
CN105571871A (zh) * 2016-01-04 2016-05-11 哈尔滨工程大学 一种在线诊断柴油机工作不均匀性的方法

Also Published As

Publication number Publication date
CN113984276A (zh) 2022-01-28

Similar Documents

Publication Publication Date Title
Naber et al. Analysis of combustion knock metrics in spark-ignition engines
US8078389B2 (en) Method and apparatus for determining a normal combustion characteristic for an internal combustion engine from an accelerometer signal
US7113861B2 (en) System and method for diagnosing and calibrating internal combustion engines
Brown et al. Determination of engine cylinder pressures from crankshaft speed fluctuations
EP2375038A1 (en) Diagnosis device and method using an in-cylinder pressure sensor in an internal combustion engine
BG63832B1 (bg) Метод за детектиране на неправилно запалване в двигател с вътрешно горене и система за осъществяване на метода
US8170777B2 (en) Indicating system and method for determining an engine parameter
CN1975363B (zh) 运行内燃机的方法和装置
US5020360A (en) Process for determining and evaluating the combustion pressure of an internal combustion engine
CN109973280A (zh) 一种发动机爆震自动标定方法
US9046048B2 (en) Method of fuel quality determination
CA2525020C (en) A combustion state detecting apparatus for an engine
CN109211572A (zh) 一种整车状态下发动机实时扭矩的在线检测装置及方法
CN113984276B (zh) 检测怠速工况下发动机各气缸压力不均匀性的方法
EP2681682B1 (en) Using ion current signal for soot and in-cylinder variable measuring techniques in internal combustion engines and methods for doing the same
US5758307A (en) Normalized misfire detection method
CN208595993U (zh) 一种整车状态下发动机实时扭矩的在线检测装置
CN103790720A (zh) 用于识别在内燃发动机中所实际使用的燃料类型的方法
US6595043B2 (en) Pressure measurement system
CN112267998B (zh) 一种多缸柴油机缸内最高温度及各缸不均匀性测试方法
CN114689324A (zh) 标定点火提前角的方法及台架试验设备和发动机控制系统
Roth et al. In-cylinder pressure measurements with optical fiber and piezoelectric pressure transducers
CN114542281B (zh) 一种基于多源数据融合的柴油机失火故障识别方法
Zhang et al. Diagnosis of the working unevenness of each cylinder by the transient crankshaft speed
KR20230163837A (ko) 불꽃 점화 엔진의 토크 모델 보정 장치 및 방법

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