CN113688508A - 一种用于多相区板带轧制的变形抗力预测方法 - Google Patents

一种用于多相区板带轧制的变形抗力预测方法 Download PDF

Info

Publication number
CN113688508A
CN113688508A CN202110874036.2A CN202110874036A CN113688508A CN 113688508 A CN113688508 A CN 113688508A CN 202110874036 A CN202110874036 A CN 202110874036A CN 113688508 A CN113688508 A CN 113688508A
Authority
CN
China
Prior art keywords
deformation resistance
phase
temperature
dislocation density
deformation
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.)
Withdrawn
Application number
CN202110874036.2A
Other languages
English (en)
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.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
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 University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN202110874036.2A priority Critical patent/CN113688508A/zh
Publication of CN113688508A publication Critical patent/CN113688508A/zh
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00—Investigating or analyzing materials by the use of thermal means
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Immunology (AREA)
  • Control Of Metal Rolling (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Probability & Statistics with Applications (AREA)
  • Computer Hardware Design (AREA)
  • Pathology (AREA)
  • Geometry (AREA)
  • Analytical Chemistry (AREA)

Abstract

本发明实施例公开一种用于多相区板带轧制的变形抗力预测方法，属于冶金机械、自动化及轧制的技术领域。包括制备试样，利用热膨胀实验获得冷却相变的开始温度和结束温度，通过插值方法确定两相区给定温度下的各相体积分数；然后基于位错密度和流变应力得到流变应力与变形条件之间的关系表达式，并通过高温压缩实验数据回归得到奥氏体区和铁素体区不同温度下的材料常数、硬化系数和软化系数，两相区的变形抗力则基于混合法则计算得出。本发明考虑了不同相区的体积分数、位错密度和流变应力等对热轧带钢变形抗力的影响，从而确定了两相区的变形抗力与前述因素之间的数值关系，实现了精准预测，确保了热精轧过程中无取向钢的热轧稳定性。

Description

一种用于多相区板带轧制的变形抗力预测方法
技术领域
本发明属于冶金机械、自动化及轧制技术的技术领域，涉及一种用于多相区板带轧制的变形抗力预测方法。
背景技术
在钢铁工业中，制定合理的轧制规程是获得高质量产品的关键。如何准确确定不同变形条件下的金属变形抗力，则是制定合理轧制规程的必要条件。
尤其是对于生产工艺窗口窄、相变温度较高的钢种，如电工钢等，在热精轧过程中不可避免地存在相变，往往上游机架位于奥氏体区，中间机架落入两相区，下游机架则变成铁素体区。由于各相区组织结构不同，材料变形抗力随温度的变化特点也因而存在差异，并不能一概而论。
然而，现有的热轧模型主要针对奥氏体单相区轧制，鲜有考虑相变对热轧板形的影响，故而在实际生产中经常出现质量问题，严重影响热精轧过程中无取向钢的热轧稳定性。所以需要针对多相区轧制建立专门的变形抗力预测方法，以保证轧制力的预测精度，提高厚度控制精度和轧制稳定性。
申请公布号CN105290119A的发明专利申请公开了基于相变的热轧带钢变形抗力预测方法，是直接针对钢种对实绩数据进行回归分析，是仅考虑某温度段与变形抗力的影响关系，并未考虑相变即不同相区(包括奥氏体区、铁素体区)对变形抗力的影响；故而多相区板带轧制的变形抗力预测并不精准和具体。
综上，现有技术并不存在准确确定不同相区(包括奥氏体区、铁素体区)对变形抗力影响的方法，至于具体不同相区的流变应力与变形条件关系的表达式和两相区的变形抗力的精准预测公式的获得以及二者之间关系的调整也不存在。
发明内容
本发明解决的技术问题是现有的基于相变的热轧带钢变形抗力预测方法仅考虑某温度段与变形抗力的影响关系，并未考虑奥氏体区、铁素体区的不同体积分数、位错密度和流变应力等对热轧带钢变形抗力的影响，预测方法并不精确。
为解决上述技术问题，本发明提供如下技术方案：
一种用于多相区板带轧制的变形抗力预测方法，所述变形抗力预测方法具体包括如下步骤：
步骤一：通过热膨胀实验获得膨胀量-温度曲线，在膨胀量-温度曲线选择冷却相变的开始温度Fs和结束温度Ff；
步骤二：通过插值方法获得两相区给定温度下的奥氏体体积分数1-X和铁素体体积分数X；
步骤三：采用位错密度ρ表征高温下奥氏体区和铁素体区的微观组织演变，流变应力σ表征宏观力学变形；
步骤四：获得流变应力σ与变形条件之间的关系；
步骤五：两相区的变形抗力σA+F根据奥氏体区和铁素体区的流变应力的混合法则计算得到，从而实现多相区板带轧制的变形抗力预测。
优选地，所述步骤一中的通过热膨胀实验获得膨胀量-温度曲线为制备试样、热处理、热膨胀仪自动记录得到膨胀量-温度曲线。
优选地，所述步骤一中的制备试样为制备Φ4mm×10mm的圆柱试样。
优选地，所述步骤一中的热处理为以5-10℃/s的升温速率加热至奥氏体化温度1120-1200℃，保温10min，然以10℃/s的冷却速率冷却至室温。
优选地，所述步骤一中的冷却相变的开始温度和结束温度是通过选择热膨胀仪自动记录的膨胀量-温度曲线冷却段直线部分的延长线与曲线部分的分离点对应的温度作获得的。
优选地，所述步骤三中的位错密度包括加工硬化引起位错增殖的加工硬化位错密度和动态软化引起位错湮灭的动态软化位错密度。两种位错密度综合作用，共同决定热变形过程中流变应力的大小及变化特点。
优选地，所述步骤三中的位错密度ρ采用以下微分方程来表示，
Figure BDA0003189694110000031
代表加工硬化效果，
Figure BDA0003189694110000032
代表动态软化效果；
Figure BDA0003189694110000033
其中：
加工硬化过程中，位错密度随变形量的增加而线性增加，表示为：
Figure BDA0003189694110000034
动态软化过程中，位错密度随时间的增加而减少，表示为：
Figure BDA0003189694110000035
式中，h为硬化系数，表示位错密度的增殖速率；s为软化系数，表示位错密度的湮灭速率。
优选地，所述步骤三中的位错密度的表达式：
Figure BDA0003189694110000036
式中，
Figure BDA0003189694110000037
为应变速率；ε为应变量；ρ0为变形开始前的位错密度值，是为了满足t＝0时刻ρ＝ρ0的初值条件。
优选地，所述步骤三中的流变应力的表达式：
σ＝G·ρ0.5，
Figure BDA0003189694110000039
G＝αμb，
式中，σ0为屈服强度，MPa；α为材料常数；μ为剪切模量，取8×104MPa；b为伯格矢量，取2.6×10-10m；
所述步骤四中的流变应力σ与变形条件之间的关系如下：
Figure BDA0003189694110000038
优选地，基于高温压缩实验测得不同温度下的应力-应变数据，采用多元回归分析，确定奥氏体区A和铁素体区F不同温度下各自对应的材料常数α、硬化系数h和软化系数s值。
优选地，所述步骤五中的两相区的变形抗力的表达式：
σA+F＝σA(1-X)+σFX，
式中，X为铁素体的体积分数。
本发明实施例提供的上述技术方案，至少具有如下有益效果：
上述方案中，用于多相区板带轧制的变形抗力预测方法考虑了奥氏体区、铁素体区的不同体积分数、位错密度和流变应力等对热轧带钢变形抗力的影响，从而确定了两相区的变形抗力与前述因素之间的数值关系，实现了精准预测，确保了热精轧过程中无取向钢的热轧稳定性，利于工业多相区板带轧制，成品率高。
综上，用于多相区板带轧制的变形抗力预测方法能够预测存在多相区轧制的特殊钢种或薄规格板带的变形抗力，提高质量控制精度和轧制稳定性。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的多相区板带轧制的变形抗力预测方法在冷却速率10℃/s的膨胀量-温度曲线；
图2为本发明的多相区板带轧制的变形抗力预测方法在不同温度下的真应力-真应变曲线；
图3为本发明的多相区板带轧制的变形抗力预测方法的两相区高温变形的实验应力应变数据与模型预测值的真应力-真应变曲线。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。
如图1-3所示，本发明实施例提供了一种用于多相区板带轧制的变形抗力预测方法，包括如下步骤：
步骤一：通过热膨胀实验获得膨胀量-温度曲线，在膨胀量-温度曲线选择冷却相变的开始温度Fs和结束温度Ff；
所述膨胀量-温度曲线是由制备试样、热处理、热膨胀仪自动记录而得到。
具体步骤如下：
制备试样：试样选择板带轧制的无取向电工钢，制备Φ4mm×10mm的圆柱试样。
热处理：以10℃/s的升温速率加热至奥氏体化温度1120℃，保温10min，然以10℃/s的冷却速率冷却至室温。
如图1所示，冷却结束后，取热膨胀仪自动记录的膨胀量-温度曲线冷却段直线部分的延长线与曲线部分的分离点对应的温度作为冷却相变的开始温度Fs＝965℃和结束温度Ff＝940℃。
步骤二：通过插值方法获得两相区给定温度下的奥氏体体积分数1-X和铁素体体积分数X；
利用膨胀量-温度曲线插值的方法，可知960℃时奥氏体、铁素体的体积分数分别为80％、20％；故而X＝0.2。
步骤三：采用位错密度ρ表征高温下奥氏体区和铁素体区的微观组织演变，流变应力σ表征宏观力学变形。
所述位错密度包括加工硬化引起位错增殖的加工硬化位错密度和动态软化引起位错湮灭的动态软化位错密度。两种位错密度综合作用，共同决定热变形过程中流变应力的大小及变化特点。
具体地：基于位错密度理论建立统一变形抗力模型。
位错密度ρ可采用如下微分方程来表示，以分别体现两种微观演变行为的作用效果。其中，
Figure BDA0003189694110000051
代表加工硬化效果，
Figure BDA0003189694110000052
代表动态软化效果。
Figure BDA0003189694110000061
其中：
加工硬化过程中，位错密度随变形量的增加而线性增加，表示为：
Figure BDA0003189694110000062
动态软化过程中，位错密度随时间的增加而减少，表示为：
Figure BDA0003189694110000063
其中，h为硬化系数，表示位错密度的增殖速率；s为软化系数，表示位错密度的湮灭速率。
将式(2)和(3)代入式(1)中，积分后得到位错密度的表达式：
Figure BDA0003189694110000064
式中，ρ0为变形开始前的位错密度值，是为了满足t＝0时刻ρ＝ρ0的初值条件；ε为应变量；
Figure BDA0003189694110000065
为应变速率。
另外，流变应力σ可表示为：
σ＝G·ρ0.5 (5)
Figure BDA0003189694110000066
G＝αμb (7)
式中σ0为屈服强度，MPa；α为材料常数；μ为剪切模量，取8×104MPa；b为伯格矢量，取2.6×10-10m。
步骤四：获得流变应力σ与变形条件之间的关系；
将式(5)、式(6)和式(7)代入式(4)中，得到表征流变应力与变形条件关系的表达式：
Figure BDA0003189694110000067
如图2所示，基于高温压缩实验测得不同温度下的应力-应变数据，实验采用试样为Φ10mm×15mm的小圆柱，以10℃/s的速率加热到1200℃，保温5min以实现完全奥氏体化；然后以5℃/s的冷却速率冷却至变形温度(800℃、850℃、900℃、960℃、1000℃、1050℃、1100℃)，再以10s-1的变形速率进行压缩变形(为了使实验效果更接近实际轧制工况，应变速率
Figure BDA0003189694110000073
定为10s-1)，应变量ε＝60％，变形完成后水冷至室温。
采用多元回归分析，确定不同相区、不同温度下的材料常数α、硬化系数h和软化系数s值，具体如下表1所述。
表1不同相区、不同温度下的特征参数
Figure BDA0003189694110000071
结合上述公式，采用线性回归方法得到各特征参数与温度的函数关系。
对于奥氏体区：
Figure BDA0003189694110000072
对于铁素体区：
Figure BDA0003189694110000081
将T＝960℃带入以上(9)和(10)，得：
σA＝208.79，σF＝57.58。
步骤五：两相区的变形抗力σA+F根据奥氏体区和铁素体区的流变应力的混合法则计算得到，从而实现多相区板带轧制的变形抗力预测。
根据混合法则，两相区的变形抗力可表示为：
σA+F＝σA(1-X)+σFX (11)
式中，X为铁素体的体积分数，根据上述计算，X＝0.2。
将以上计算的σA，σF，X带入上式，得到960℃的带钢变形抗力：σA+F＝178.54。
如图3所示，将两相区高温变形的实验应力应变数据与模型预测值作对比，结果发现计算值与实验值吻合较好，相关系数平方值R2达到0.95。
上述方案中，用于多相区板带轧制的变形抗力预测方法考虑了奥氏体区、铁素体区的不同体积分数、位错密度和流变应力等对热轧带钢变形抗力的影响，从而确定了两相区的变形抗力与前述因素之间的数值关系，实现了精准预测，确保了热精轧过程中无取向钢的热轧稳定性，利于工业多相区板带轧制，成品率高。
综上，用于多相区板带轧制的变形抗力预测方法能够预测存在多相区轧制的特殊钢种或薄规格板带的变形抗力，提高质量控制精度和轧制稳定性。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种用于多相区板带轧制的变形抗力预测方法，其特征在于，所述变形抗力预测方法具体包括如下步骤：
步骤一：通过热膨胀实验获得膨胀量-温度曲线，在膨胀量-温度曲线选择冷却相变的开始温度Fs和结束温度Ff；
步骤二：通过插值方法获得两相区给定温度下的奥氏体体积分数1-X和铁素体体积分数X；
步骤三：采用位错密度ρ表征高温下奥氏体区和铁素体区的微观组织演变，流变应力σ表征宏观力学变形；
步骤四：获得流变应力σ与变形条件之间的关系；
步骤五：两相区的变形抗力σA+F根据奥氏体区和铁素体区的流变应力的混合法则计算得到，从而实现多相区板带轧制的变形抗力预测。
2.根据权利要求1所述的用于多相区板带轧制的变形抗力预测方法，其特征在于，所述步骤一中的热膨胀实验为制备试样、热处理、热膨胀仪自动记录得到膨胀量-温度曲线。
3.根据权利要求2所述的用于多相区板带轧制的变形抗力预测方法，其特征在于，所述步骤一中的制备试样为制备Φ4mm×10mm的圆柱试样。
4.根据权利要求2所述的用于多相区板带轧制的变形抗力预测方法，其特征在于，所述步骤一中的热处理为以5-10℃/s的升温速率加热至奥氏体化温度1120-1200℃，保温10min，然以10℃/s的冷却速率冷却至室温。
5.根据权利要求1所述的用于多相区板带轧制的变形抗力预测方法，其特征在于，所述步骤一中的冷却相变的开始温度和结束温度是通过选择热膨胀仪自动记录的膨胀量-温度曲线冷却段直线部分的延长线与曲线部分的分离点对应的温度作获得的。
6.根据权利要求1所述的用于多相区板带轧制的变形抗力预测方法，其特征在于，所述步骤三中的位错密度包括加工硬化引起位错增殖的加工硬化位错密度和动态软化引起位错湮灭的动态软化位错密度。
7.根据权利要求6所述的用于多相区板带轧制的变形抗力预测方法，其特征在于，所述步骤三中的位错密度ρ采用以下微分方程来表示，
Figure FDA0003189694100000021
代表加工硬化效果，
Figure FDA0003189694100000022
代表动态软化效果；
Figure FDA0003189694100000023
其中：
加工硬化过程中，位错密度随变形量的增加而线性增加，表示为：
Figure FDA0003189694100000024
动态软化过程中，位错密度随时间的增加而减少，表示为：
Figure FDA0003189694100000025
式中，h为硬化系数，表示位错密度的增殖速率；s为软化系数，表示位错密度的湮灭速率。
8.根据权利要求7所述的用于多相区板带轧制的变形抗力预测方法，其特征在于，所述步骤三中的位错密度的表达式：
Figure FDA0003189694100000026
式中，
Figure FDA0003189694100000027
为应变速率；ε为应变量；ρ0为变形开始前的位错密度值，是为了满足t＝0时刻ρ＝ρ0的初值条件。
9.根据权利要求8所述的用于多相区板带轧制的变形抗力预测方法，其特征在于，所述步骤三中的流变应力σ的表达式：
σ＝G·ρ0.5，
Figure FDA0003189694100000028
G＝αμb，
式中，σ0为屈服强度，MPa；α为材料常数；μ为剪切模量，取8×104MPa；b为伯格矢量，取2.6×10-10m；
所述步骤四中的流变应力σ与变形条件之间的关系如下：
Figure FDA0003189694100000031
10.根据权利要求9所述的用于多相区板带轧制的变形抗力预测方法，其特征在于，所述步骤五中的两相区的变形抗力的表达式：
σA+F＝σA(1-X)+σFX，
式中，X为铁素体的体积分数。
CN202110874036.2A 2021-07-30 2021-07-30 一种用于多相区板带轧制的变形抗力预测方法 Withdrawn CN113688508A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110874036.2A CN113688508A (zh) 2021-07-30 2021-07-30 一种用于多相区板带轧制的变形抗力预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110874036.2A CN113688508A (zh) 2021-07-30 2021-07-30 一种用于多相区板带轧制的变形抗力预测方法

Publications (1)

Publication Number Publication Date
CN113688508A true CN113688508A (zh) 2021-11-23

Family

ID=78578446

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110874036.2A Withdrawn CN113688508A (zh) 2021-07-30 2021-07-30 一种用于多相区板带轧制的变形抗力预测方法

Country Status (1)

Country Link
CN (1) CN113688508A (zh)

Similar Documents

Publication Publication Date Title
CN103143563B (zh) 一种利用相变热膨胀防止65Mn热轧板卷塌卷的方法
Li et al. Experimental study and numerical simulation of dynamic recrystallization behavior of a micro-alloyed plastic mold steel
CN109977442B (zh) 一种超厚板坯多道次轧制工艺的数值模拟方法
CN1330930C (zh) 轧制过程钢板内部组织晶粒尺寸的软测量方法
JPWO2008012881A1 (ja) 圧延ラインの材質予測および材質制御装置
Montheillet et al. A critical assessment of three usual equations for strain hardening and dynamic recovery
CN104238498A (zh) 热连轧管线钢组织及力学性能预测的方法
CN115815345A (zh) 预测全流程热轧带钢力学性能的机理协同预报方法及系统
Sun et al. Thermo-mechanical coupled analysis of hot ring rolling process
Ganguly et al. Modeling and simulation of mass flow during hot rolling low carbon steel I-beam
Qiao et al. Static recrystallization behavior of SA508Gr. 4N reactor pressure vessel steel during hot compressive deformation
Li et al. Study on hot deformation behavior and bending forging process of 7075 aluminum alloy
CN109856180A (zh) 一种含硼微合金钢热成形工艺参数的确定方法
CN107818184B (zh) 一种构建材料变形抗力模型的方法
CN113688508A (zh) 一种用于多相区板带轧制的变形抗力预测方法
CN107423469B (zh) 一种06Cr19Ni9NbN钢锻透的判定方法
CN102750453B (zh) 一种用于轧机设计的材料变形抗力统一模型的构造方法
CN115502233A (zh) 板坯加热参数的获取方法、板坯加热方法及钢板板坯
Ashtiani et al. Flow Behavior and Metallurgical Phenomena of Micro-alloy Steel Under Elevated Temperature Conditions
Li et al. Constitutive modeling of flow behavior and processing maps of a low-carbon steel
Kruglova et al. Improvement of chemical composition and production regimes for manufacture of K65–K70 (X80–X90) strip based on simulation
KR100340544B1 (ko) 페라이트계스테인레스강의열간변형저항예측방법
Wang et al. Characteristics and prevention of thermomechanical controlled process plate deflection resulting from uneven cooling
Serajzadeh et al. Effects of deformation parameters on the final microstructure and mechanical properties in warm rolling of a low-carbon steel
CN110119542B (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
WW01 Invention patent application withdrawn after publication

Application publication date: 20211123

WW01 Invention patent application withdrawn after publication