CN113076670A - 一种多物理场耦合的相控阵制导微系统协同优化方法 - Google Patents

一种多物理场耦合的相控阵制导微系统协同优化方法 Download PDF

Info

Publication number
CN113076670A
CN113076670A CN202110325777.5A CN202110325777A CN113076670A CN 113076670 A CN113076670 A CN 113076670A CN 202110325777 A CN202110325777 A CN 202110325777A CN 113076670 A CN113076670 A CN 113076670A
Authority
CN
China
Prior art keywords
phased array
guidance
physical
micro
analysis model
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
CN202110325777.5A
Other languages
English (en)
Other versions
CN113076670B (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.)
Northwestern Polytechnical University
Guizhou Aerospace Electronic Technology Co Ltd
Original Assignee
Northwestern Polytechnical University
Guizhou Aerospace Electronic 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 Northwestern Polytechnical University, Guizhou Aerospace Electronic Technology Co Ltd filed Critical Northwestern Polytechnical University
Priority to CN202110325777.5A priority Critical patent/CN113076670B/zh
Publication of CN113076670A publication Critical patent/CN113076670A/zh
Application granted granted Critical
Publication of CN113076670B publication Critical patent/CN113076670B/zh
Active 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
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/06Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]

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)
  • Radar Systems Or Details Thereof (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明提供了一种多物理场耦合的相控阵制导微系统协同优化方法,分别建立相控阵制导微系统的结构有限元模型、热分析模型和力分析模型;建立相应的结构温差载荷;得到结构位移场;确定辐射单元的空间相位误差;获得相控阵制导微系统天线单元的和差方向图;计算相控阵天线的电性能参数;得到相控阵制导系统的指向精度;判断指向精度,如果满足要求则系统协同优化完成,否则,修改结构有限元模型、热分析模型和力分析模型的设计参数,直到满足要求。本发明提供了一种多物理场耦合的相控阵制导微系统协同优化方法,通过对热场、力场、电磁场和电路结构四类物理量耦合,优化系统指向精度,来保证制导系统复杂应用环境中的探测性能和可靠性。

Description

一种多物理场耦合的相控阵制导微系统协同优化方法
技术领域
本发明属于天线技术领域,具体涉及一种多物理场耦合的相控阵制导微系统协同优化方法。
背景技术
相控阵制导微系统是利用先进相控阵技术,通过波控计算机控制电移相器的相移量实现雷达波束的扫描的一种新型雷达制导系统。在整个相控阵制导微系统工作过程中,由于其自身结构特点和运行环境,热场、力场、电磁场同时存在,且各场并不是独立对“电路”产生影响的,它是一个叠加的效果,多物理场协同作用使得波束指向误差加剧,这也将严重影响制导系统的性能。
传统的相控阵制导微系统由于存在功能结构复杂、技术性能先进、信息密集、高度集成、协同复杂等诸多特征,在优化分析中一般考虑单一的热仿真或力学仿真对电磁场的影响,利用电磁场对天线方向图的影响,对系统热学或力学参数进行优化。如公开号为CN108920831A的中国专利,公开了高速飞行器罩体高温烧蚀对天线电性能影响的快速计算方法,仅仅是对热场的影响进行了分析,该方法未考虑多载荷应力耦合的工程应用环境,也不能考虑热力之间的相互影响,导致分析结果不准确。此外,仅考虑天线方向图对系统的影响,未考虑捷联去耦技术,也会导致相控阵制导系统的指向误差模型不准确。通过上述分析可知,传统的多物理场耦合协同优化方法难以在满足相控阵制导微系统可靠性设计要求。
发明内容
为解决上述技术问题,本发明提供了一种多物理场耦合的相控阵制导微系统协同优化方法,通过对热场、力场、电磁场和电路结构四类物理量耦合,优化系统指向精度,来保证制导系统复杂应用环境中的探测性能和可靠性。
本发明通过以下技术方案得以实现:
一种多物理场耦合的相控阵制导微系统协同优化方法,包括以下步骤:
步骤一、分别建立相控阵制导微系统的结构有限元模型、热分析模型和力分析模型。
步骤二、根据热分析模型建立相应的结构温差载荷。
步骤三、利用步骤一中得到的三种模型,共同建立相控阵制导微系统的约束和边界条件,并根据应力条件和步骤二中的结构温差载荷,得到包括阵面辐射单元的位置偏移量在内的结构位移场。
步骤四、根据阵面辐射单元的位置偏移,确定辐射单元的位置变化在远区目标处导致的空间相位误差。
步骤五、根据步骤四中的空间相位误差,结合结构温差载荷获得TR芯片激励电流的幅度变化量和相位变化量,进而获得相控阵制导微系统天线单元的和差方向图。
步骤六、利用步骤五中和差方向图的变化,结合阵面辐射函数以及幅相分布、天线阵面辐射单元的排列形式,计算相控阵天线的电性能参数。
步骤七、根据步骤六中的电性能参数的变化,得到相控阵制导系统的指向精度。
步骤八、判断相控阵制导系统的指向精度是否满足要求,如果满足要求则系统协同优化完成,否则,修改结构有限元模型、热分析模型和力分析模型的设计参数,并重复步骤一至八,直到满足要求。
所述结构有限元模型的设计参数包括元器件类型及其位置、尺寸、管脚和功耗,电路板的层数、厚度、镀孔信息。
所述力分析模型的设计参数至少包括机械振动。
所述热分析模型的设计参数包括环境温度。
所述步骤一中,根据相控阵制导微系统的结构参数,确定相控阵制导微系统的有限元模型,得到阵面每一个辐射单元的理论坐标。
所述步骤二中,根据热分析模型建立相应的结构温差载荷包括以下步骤:
步骤二一、建立相控阵制导微系统的热分析模型,确定散热设计参数,并对该模型施加系统工作要求参数,得到相控阵制导微系统、相控阵天线和相控阵芯片的温度分布和温升曲线。
步骤二二、根据相控阵收发芯片的温度分布,以及实测得到的相控阵收发芯片的温度电流曲线图,得到相控阵收发芯片激励电流的幅度变化量和相位变化量。
步骤二三、设定相控阵制导微系统的参考温度,并将相控阵制导微系统、相控阵天线辐射单元和相控阵芯片的温度分布与该参考温度比较,确定相应的结构温差载荷。
所述步骤六中的电性能参数包括副瓣电平、波束指向和增益。
所述步骤七中和差方向图变化通过系统的和差测角方法和捷联去耦技术计算得到相控阵制导系统的指向精度。
所述相控阵为有源相控阵。
本发明的有益效果在于:
与现有技术相比,通过对热场、力场、电磁场和电路结构四类物理量耦合,利用系统工作过程中的温升曲线、温度梯度和振动谱线,深入研究热场、力场、电磁场耦合关系,结合测角原理和捷联去耦方法建立热应力影响下的系统指向误差模型,并通过热场、力场、电路结构参数优化,优化系统指向精度,来保证相控阵末制导微系统在复杂应用环境中的探测性能和可靠性,使其角度测量精度满足系统性能要求,并提供了一种多物理场耦合条件下相控阵制导微系统协同优化方法,该方法相比现有技术影响因素考虑更加全面,与实际测量误差小,精度更高。
附图说明
图1是本发明中多物理场耦合的相控阵制导微系统协同优化方法的路线图。
具体实施方式
下面进一步描述本发明的技术方案,但要求保护的范围并不局限于所述。
如图1所示,一种多物理场耦合的相控阵制导微系统协同优化方法,包括以下步骤:
步骤一、分别建立有源相控阵制导微系统的结构有限元模型、热分析模型和力分析模型,并根据相控阵制导微系统的结构参数,确定相控阵制导微系统的有限元模型,得到阵面每一个辐射单元的理论坐标。
所述结构有限元模型的设计参数包括元器件类型及其位置、尺寸、管脚和功耗,电路板的层数、厚度、镀孔信息;所述力分析模型的设计参数至少包括机械振动;所述热分析模型的设计参数包括环境温度。
步骤二、根据热分析模型建立相应的结构温差载荷,具体包括以下步骤:
步骤二一、建立相控阵制导微系统的热分析模型,确定散热设计参数,并对该模型施加系统工作要求参数,得到相控阵制导微系统、相控阵天线和相控阵芯片的温度分布和温升曲线。
步骤二二、根据相控阵收发芯片的温度分布,以及实测得到的相控阵收发芯片的温度电流曲线图,得到相控阵收发芯片激励电流的幅度变化量和相位变化量。
步骤二三、设定相控阵制导微系统的参考温度,并将相控阵制导微系统、相控阵天线辐射单元和相控阵芯片的温度分布与该参考温度比较,确定相应的结构温差载荷。
步骤三、利用步骤一中得到的三种模型,共同建立相控阵制导微系统的约束和边界条件,并根据应力条件和步骤二中的结构温差载荷,得到包括阵面辐射单元的位置偏移量在内的结构位移场。
步骤四、根据阵面辐射单元的位置偏移,确定辐射单元的位置变化在远区目标处导致的空间相位误差。
步骤五、根据步骤四中的空间相位误差,结合结构温差载荷获得TR芯片激励电流的幅度变化量和相位变化量,进而获得相控阵制导微系统天线单元的和差方向图。
步骤六、利用步骤五中和差方向图的变化,结合阵面辐射函数以及幅相分布、天线阵面辐射单元的排列形式,计算相控阵天线的副瓣电平、波束指向和增益等电性能参数。
步骤七、根据步骤六中的电性能参数的变化,通过系统的和差测角方法和捷联去耦技术计算得到相控阵制导系统的指向精度。
步骤八、判断相控阵制导系统的指向精度是否满足要求,如果满足要求则系统协同优化完成,否则,修改结构有限元模型、热分析模型和力分析模型的设计参数,并重复步骤一至八,直到满足要求。
本发明提供的一种多物理场耦合的相控阵制导微系统协同优化方法,通过对热场、力场、电磁场和电路结构四类物理量耦合,利用系统工作过程中的温升曲线、温度梯度和振动谱线,深入研究热场、力场、电磁场耦合关系,结合测角原理和捷联去耦方法建立热应力影响下的系统指向误差模型,并通过热场、力场、电路结构参数优化,优化系统指向精度,来保证相控阵末制导微系统在复杂应用环境中的探测性能和可靠性,使其角度测量精度满足系统性能要求,并提供了一种多物理场耦合条件下相控阵制导微系统协同优化方法,该方法相比现有技术影响因素考虑更加全面,与实际测量误差小,精度更高。

Claims (9)

1.一种多物理场耦合的相控阵制导微系统协同优化方法,其特征在于,包括以下步骤:
步骤一、分别建立相控阵制导微系统的结构有限元模型、热分析模型和力分析模型;
步骤二、根据热分析模型建立相应的结构温差载荷;
步骤三、利用步骤一中得到的三种模型,共同建立相控阵制导微系统的约束和边界条件,并根据应力条件和步骤二中的结构温差载荷,得到包括阵面辐射单元的位置偏移量在内的结构位移场;
步骤四、根据阵面辐射单元的位置偏移,确定辐射单元的位置变化在远区目标处导致的空间相位误差;
步骤五、根据步骤四中的空间相位误差,结合结构温差载荷获得TR芯片激励电流的幅度变化量和相位变化量,进而获得相控阵制导微系统天线单元的和差方向图;
步骤六、利用步骤五中和差方向图的变化,结合阵面辐射函数以及幅相分布、天线阵面辐射单元的排列形式,计算相控阵天线的电性能参数;
步骤七、根据步骤六中的电性能参数的变化,得到相控阵制导系统的指向精度;
步骤八、判断相控阵制导系统的指向精度是否满足要求,如果满足要求则系统协同优化完成,否则,修改结构有限元模型、热分析模型和力分析模型的设计参数,并重复步骤一至八,直到满足要求。
2.如权利要求1所述的一种多物理场耦合的相控阵制导微系统协同优化方法,其特征在于:所述结构有限元模型的设计参数包括元器件类型及其位置、尺寸、管脚和功耗,电路板的层数、厚度、镀孔信息。
3.如权利要求1所述的一种多物理场耦合的相控阵制导微系统协同优化方法,其特征在于:所述力分析模型的设计参数至少包括机械振动。
4.如权利要求1所述的一种多物理场耦合的相控阵制导微系统协同优化方法,其特征在于:所述热分析模型的设计参数包括环境温度。
5.如权利要求1所述的一种多物理场耦合的相控阵制导微系统协同优化方法,其特征在于:所述步骤一中,根据相控阵制导微系统的结构参数,确定相控阵制导微系统的有限元模型,得到阵面每一个辐射单元的理论坐标。
6.如权利要求1所述的一种多物理场耦合的相控阵制导微系统协同优化方法,其特征在于:所述步骤二中,根据热分析模型建立相应的结构温差载荷包括以下步骤:
步骤二一、建立相控阵制导微系统的热分析模型,确定散热设计参数,并对该模型施加系统工作要求参数,得到相控阵制导微系统、相控阵天线和相控阵芯片的温度分布和温升曲线;
步骤二二、根据相控阵收发芯片的温度分布,以及实测得到的相控阵收发芯片的温度电流曲线图,得到相控阵收发芯片激励电流的幅度变化量和相位变化量;
步骤二三、设定相控阵制导微系统的参考温度,并将相控阵制导微系统、相控阵天线辐射单元和相控阵芯片的温度分布与该参考温度比较,确定相应的结构温差载荷。
7.如权利要求1所述的一种多物理场耦合的相控阵制导微系统协同优化方法,其特征在于:所述步骤六中的电性能参数包括副瓣电平、波束指向和增益。
8.如权利要求1所述的一种多物理场耦合的相控阵制导微系统协同优化方法,其特征在于:所述步骤七中和差方向图变化通过系统的和差测角方法和捷联去耦技术计算得到相控阵制导系统的指向精度。
9.如权利要求1所述的一种多物理场耦合的相控阵制导微系统协同优化方法,其特征在于:所述相控阵为有源相控阵。
CN202110325777.5A 2021-03-26 2021-03-26 一种多物理场耦合的相控阵制导微系统协同优化方法 Active CN113076670B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110325777.5A CN113076670B (zh) 2021-03-26 2021-03-26 一种多物理场耦合的相控阵制导微系统协同优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110325777.5A CN113076670B (zh) 2021-03-26 2021-03-26 一种多物理场耦合的相控阵制导微系统协同优化方法

Publications (2)

Publication Number Publication Date
CN113076670A true CN113076670A (zh) 2021-07-06
CN113076670B CN113076670B (zh) 2022-10-18

Family

ID=76610481

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110325777.5A Active CN113076670B (zh) 2021-03-26 2021-03-26 一种多物理场耦合的相控阵制导微系统协同优化方法

Country Status (1)

Country Link
CN (1) CN113076670B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113901663A (zh) * 2021-10-14 2022-01-07 中国电子科技集团公司第三十八研究所 一种瓦片式相控阵天线机电热耦合设计方法

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101344564A (zh) * 2008-08-14 2009-01-14 西安电子科技大学 基于机电热三场耦合的有源相控阵天线电性能预测方法
CN102998540A (zh) * 2012-10-22 2013-03-27 西安电子科技大学 共形承载微带天线阵面形貌对电性能影响的预测方法
CN104101786A (zh) * 2014-06-24 2014-10-15 中国电子科技集团公司第十研究所 全空域有源多波束球面相控阵天线方向图测量系统
CN105095603A (zh) * 2015-09-09 2015-11-25 哈尔滨工业大学 一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法
CN105757888A (zh) * 2016-03-09 2016-07-13 华为技术有限公司 一种精密空调内外机联动的控制方法及装置
CN105787160A (zh) * 2016-02-17 2016-07-20 西安电子科技大学 基于机电耦合的星载有源相控阵天线结构轻量化设计方法
CN106484984A (zh) * 2016-09-29 2017-03-08 北京大学 星载柔性附件热致微振动响应仿真分析平台
CN106549235A (zh) * 2016-10-21 2017-03-29 四川九洲电器集团有限责任公司 一种相控阵列天线方向图性能参数优化方法及电子设备
CN106898395A (zh) * 2017-03-14 2017-06-27 北京科瑞华安科技有限公司 一种乏核燃料格架地震安全分析方法
CN107608942A (zh) * 2017-09-08 2018-01-19 北京理工大学 温度和弯矩载荷下聚合物基微带天线变形量的计算方法
CN108381742A (zh) * 2018-01-16 2018-08-10 宁波工程学院 一种长混凝土预制构件横向裂缝的控制方法
JP2019123004A (ja) * 2018-01-18 2019-07-25 Jfeスチール株式会社 熱間圧延の粗圧延方法、熱間圧延の粗圧延装置、熱延鋼板の製造方法及び熱延鋼板の製造装置
CN110414164A (zh) * 2019-08-01 2019-11-05 西安电子科技大学 一种基于机电耦合的6g通信天线信道容量预测方法
CN110781630A (zh) * 2019-11-28 2020-02-11 青海电研科技有限责任公司 一种gis设备管线结构温差应力变形在线监测方法
CN111693979A (zh) * 2020-05-25 2020-09-22 北京理工大学 一种基于泰勒展开的数字阵列单脉冲测角方法

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101344564A (zh) * 2008-08-14 2009-01-14 西安电子科技大学 基于机电热三场耦合的有源相控阵天线电性能预测方法
CN102998540A (zh) * 2012-10-22 2013-03-27 西安电子科技大学 共形承载微带天线阵面形貌对电性能影响的预测方法
CN104101786A (zh) * 2014-06-24 2014-10-15 中国电子科技集团公司第十研究所 全空域有源多波束球面相控阵天线方向图测量系统
CN105095603A (zh) * 2015-09-09 2015-11-25 哈尔滨工业大学 一种高超声速流动-传热与结构响应的多场耦合瞬态数值的方法
CN105787160A (zh) * 2016-02-17 2016-07-20 西安电子科技大学 基于机电耦合的星载有源相控阵天线结构轻量化设计方法
CN105757888A (zh) * 2016-03-09 2016-07-13 华为技术有限公司 一种精密空调内外机联动的控制方法及装置
CN106484984A (zh) * 2016-09-29 2017-03-08 北京大学 星载柔性附件热致微振动响应仿真分析平台
CN106549235A (zh) * 2016-10-21 2017-03-29 四川九洲电器集团有限责任公司 一种相控阵列天线方向图性能参数优化方法及电子设备
CN106898395A (zh) * 2017-03-14 2017-06-27 北京科瑞华安科技有限公司 一种乏核燃料格架地震安全分析方法
CN107608942A (zh) * 2017-09-08 2018-01-19 北京理工大学 温度和弯矩载荷下聚合物基微带天线变形量的计算方法
CN108381742A (zh) * 2018-01-16 2018-08-10 宁波工程学院 一种长混凝土预制构件横向裂缝的控制方法
JP2019123004A (ja) * 2018-01-18 2019-07-25 Jfeスチール株式会社 熱間圧延の粗圧延方法、熱間圧延の粗圧延装置、熱延鋼板の製造方法及び熱延鋼板の製造装置
CN110414164A (zh) * 2019-08-01 2019-11-05 西安电子科技大学 一种基于机电耦合的6g通信天线信道容量预测方法
CN110781630A (zh) * 2019-11-28 2020-02-11 青海电研科技有限责任公司 一种gis设备管线结构温差应力变形在线监测方法
CN111693979A (zh) * 2020-05-25 2020-09-22 北京理工大学 一种基于泰勒展开的数字阵列单脉冲测角方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王从思 等: "星载智能天线结构的机电热耦合优化分析", 《电波科学学报》 *
肖瑞 等: "相控阵雷达有源子阵结构设计研究", 《电子机械工程》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113901663A (zh) * 2021-10-14 2022-01-07 中国电子科技集团公司第三十八研究所 一种瓦片式相控阵天线机电热耦合设计方法

Also Published As

Publication number Publication date
CN113076670B (zh) 2022-10-18

Similar Documents

Publication Publication Date Title
CN101344564B (zh) 基于机电热三场耦合的有源相控阵天线电性能预测方法
Wang et al. Coupled structural–electromagnetic–thermal modelling and analysis of active phased array antennas
CN107038299B (zh) 一种考虑互耦效应的变形阵列天线远场方向图补偿方法
CN106021764B (zh) 面向机电耦合的有源相控阵天线性能仿真置信度的计算方法
CN104036114B (zh) 一种基于机电耦合的六边形有源相控阵天线结构公差的快速确定方法
CN109031226B (zh) 基于机电耦合的变形有源相控阵雷达探测性能快速评估方法
CN104615836A (zh) 星载微带阵列天线热变形对电性能影响的快速预测方法
CN107103124B (zh) 基于机电耦合理论的变形阵列天线远场方向图分析方法
CN102708257B (zh) 一种有源相控阵天线结构公差的快速确定方法
CN102253290B (zh) 基于机电耦合模型的变形对数周期天线电性能预测方法
CN104701637B (zh) 基于机电耦合与最小二乘法的变形阵列天线电性能补偿方法
CN113076670B (zh) 一种多物理场耦合的相控阵制导微系统协同优化方法
CN104036093A (zh) 一种基于阵元互耦的大型变形阵列天线副瓣性能预测方法
CN106991210B (zh) 一种基于机电耦合模型的赋形反射面天线电性能预测方法
CN111539143B (zh) 一种基于应变电磁耦合的有源相控阵天线电性能补偿方法
CN106886619A (zh) 相控阵天线阵面综合实现方法
Luther et al. A Low-Cost 2$\times $2 Planar Array of Three-Element Microstrip Electrically Steerable Parasitic Array Radiator (ESPAR) Subcells
CN105718662B (zh) 基于机电耦合的圆柱共形阵列天线结构公差快速确定方法
Wang et al. Electromechanical coupling based performance evaluation of distorted phased array antennas with random position errors
Rengarajan et al. Design, analysis, and development of a large Ka-band slot array for digital beam-forming application
CN115021833A (zh) 相控阵天线阵元通道一致性多模并行处理标校方法
CN110729821B (zh) 用于多目标无线能量传输的准无衍射波束形成方法
Fan et al. Optimal synthesis of sum and difference beam patterns with a common weight vector for symmetric and asymmetric antenna arrays
US8134494B1 (en) Simulating the mutual performance of an antenna array coupled to an electrical drive circuit
CN111539141B (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