CN107515982B - 一种三维力学有限元模态分析中的接触分析方法 - Google Patents
一种三维力学有限元模态分析中的接触分析方法 Download PDFInfo
- Publication number
- CN107515982B CN107515982B CN201710724845.9A CN201710724845A CN107515982B CN 107515982 B CN107515982 B CN 107515982B CN 201710724845 A CN201710724845 A CN 201710724845A CN 107515982 B CN107515982 B CN 107515982B
- Authority
- CN
- China
- Prior art keywords
- contact
- finite element
- displacement
- grid
- equation
- 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
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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/06—Power analysis or power optimisation
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)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于三维结构力学分析数值求解技术领域,涉及一种三维力学有限元模态分析中的接触分析方法。本发明首先对目标电子器件结构进行建模,引入位移边界条件或者应力边界条件建立对应的几何结构模型,并根据各部件之间的相互关系,生成各部件之间的接触关系,然后对所建几何结构模型采用四面体网格进行剖分,并根据仿真区域接触面的面网格生成联合面网格,最后利用有限元法,建立考虑接触问题的目标电子器件的有限元广义本征方程,并求解该方程获得特征值和特征向量,进行后处理获得振动模态频率和振动振型,从而实现了模态分析中的接触分析,且获得高精度的数值计算结果。
Description
技术领域
本发明属于三维结构力学分析数值求解技术领域,涉及一种三维力学有限元模态分析中的接触分析方法。
背景技术
电子器件的使用环境往往十分恶劣,例如在崎岖上路上运输时的振动、飞机起飞、坦克高速行进、卫星和导弹上升阶段的重力加速等对电子器件的机械强度提出了十分严格的要求。机械性能却又是电子器件的可靠性和稳定性的重要组成部分,这直接影响到器件能否正常工作。因此对电子器件的机械性能进行优化设计是有必要的,而模态分析可以获得电子器件的振动特性,是其机械性能设计的重要环节,因此模态分析中高精度的获得器件的振动特性具有极其重要的意义。然而大部分的器件振动是一种多体运动,部件之间存在接触关系,适当的接触分析方法是能否高精度获得器件的振动特性的关键,因此研究模态分析中的接触分析方法具有重要意义。
目前,利用各种力学计算方法对电子器件结构模态进行仿真分析时,都是采用的有限元本征分析方法。有限元分析一般包括,单元划分、单元分析、系统综合、引入条件、求解方程组和后处理等几个步骤,这是一个线性分析过程。有关接触问题的分析,大量的文献资料上描述的都是一种非线性分析过程,因此对于模态分析中的接触问题分析仅有美国ANSYS公司开发的ANSYS软件等一些商业软件做到了这一点,但是其牵涉到的技术内容由于商业机密并没有对外公布。然而商业软件由于其广泛的适应性,其在算法的效率上并不具有优势,随着分析的结构复杂化,在一定程度上并不能满足专业用户的精度和效率需求。而且由于其高度的封装性,在用户二次开发时并不能透明化的使用接触分析这个功能,这也将影响专业用户的需求。因此我们需要一种模态分析中的接触分析方法,来高精度获得器件的振动特性,同时也方便用户二次开发。
发明内容
针对上述存在问题或不足,为解决构造有限元模态分析中的接触分析,从而获得高精度的数值模拟结果;本发明提供了一种三维力学有限元模态分析中的接触分析方法。
该三维力学有限元模态分析中的接触分析方法,包括以下步骤:
A.将目标电子器件结构进行建模,引入位移边界条件或者应力边界条件建立对应的几何结构模型;
B.根据电子器件各部件之间的相互关系,生成几何结构各部件之间的接触关系;
C.对所建几何结构模型采用四面体网格进行剖分,将连续的几何结构空间转化为离散空间;
D.根据仿真区域接触面的面网格生成联合面网格;
E.利用有限元法,将平衡微分方程、几何方程和物理方程等效的结构力学边值问题在步骤C建立的网格空间进行离散,结合步骤D建立的联合面网格,建立考虑接触问题的目标电子器件的有限元方程;
F.引入几何结构的惯性力,得到其自由振动有限元广义本征方程;
H.对G步骤获得的特征值和对应特征向量进行后处理获得振动模态频率和对应振动振型。
本发明通过在弹性体刚度矩阵中添加与接触约束相应的矩阵元素来实现接触分析,在求解本征方程时,矩阵方程求解只需要一次求解,而不需要非线性迭代多次求解矩阵方程,从而能够高效的获得高精度的数值计算结果。
附图说明
图1是本发明的流程图;
图2接触面联合面网格示意图。
具体实施方式
下面结合附图和实施例来详细说明本发明的技术方案。
参照附图1,一种三维力学有限元模态分析中的接触分析方法,包括以下步骤:
A.将目标电子器件结构进行建模,引入位移边界条件或者应力边界条件建立对应的几何结构模型。
建立目标电子器件的几何模型,根据电子器件的特性,引入位移边界条件来仿真整个结构的振动特性。具体的结构建模是结构力学数值计算中的一种公知过程,因此本步骤不再详细描述。
B.根据电子器件各部件之间的相互关系,生成几何结构各部件之间的接触关系。
根据电子器件各部件的间的连接关系,在几何结构中生成相应的接触对。一般的前处理软件都能够实现该功能,这是一种公知的过程,因此本步骤不再详细阐述。
C.对所建几何结构模型采用四面体网格进行剖分,将连续的几何结构空间转化为离散空间。
采用四面体网格剖分仿真区域,剖分后的仿真区域被人为分割为多个三维四面体网格,从而将连续的几何结构空间转化为离散的网格空间。此时,接触面由一系列离散的面网格组成。由于四面体网格剖分是有限元方法中的一种公知过程,因此本步骤不再详细描述。
D.根据仿真区域接触面的面网格生成联合面网格。
在上述步骤C中生成的四面体网格在接触面上会存在一个公共面,由于接触面两边属于不同的部件,各部件具有不同的属性,因此接触面在不同的部件网格中具有不同的属性,需要建立一套联合的接触面网格。下面给出联合面网格的生成方法与步骤。
1)确定接触面网格的位置和所属四面体单元
根据步骤C生成的四面体网格信息和步骤B生成的接触对关系,确定每一个接触对上的面网格以及其所属四面体单元。
2)复制接触单元面网格
对上述D步骤中的1)步骤所确定的接触对上的面网格进行共形复制,如附图2所示面网格A'B'C'是接触单元面网格ABC的一个复制,附图2中面网格A'B'C'和面网格ABC在空间几何位置上是重叠的,之间是没有缝隙的,图中只是为了示意方便,描述成那样。这样在接触面上就形成了一套共形网格,其中面网格ABC属于四面体ABCD,面网格A'B'C'属于四面体A'B'C'D',以分别携带不同的材料属性。
E.利用有限元法,将平衡微分方程、几何方程和物理方程等效的结构力学边值问题在步骤C建立的网格空间进行离散,结合步骤D建立的联合面网格,建立考虑接触问题的目标电子器件的有限元方程。
对于空间边值问题,在结构(弹性体)内部我们要考虑静力学、几何学、物理学三方面条件,分别建立三套方程;并在给定约束或面力的边界上,建立位移边界条件或应力边界条件。具体如下:
平衡微分方程
几何方程
物理方程
位移边界条件
应力边界条件
上述(1)(2)(3)(4)(5)式中,σx,σy,σz,τxy=τyx,τyz=τzy,τzx=τxz表示求解区域中6个应力分量,εx,εy,εz,γxy,γyz,γzx表示求解区域中6个形变应力分量,u,v,w表示求解区域中3个位移分量。E是求解区域中结构的杨氏弹性模量,μ是求解区域中结构的泊松比,Su表示位移边界面,Sσ表示应力边界面。a=cos(n′,x),b=cos(n′,y),c=cos(n′,z),表示应力边界面Sσ上的方向余弦,其中n′为应力边界面Sσ的外法线,x,y,z为应力边界面Sσ上三个方向的坐标值。为位移边界面Su上的位移值,fx,fy,fz为求解区域内结构受到的各个方向的体力,为应力边界面Sσ上受到的各个方向的面力,具体推导过程为一种公知过程,这里不再阐述。
经过有限元推导我们得到如下每一个四面体单元中的有限元方程,有限元法是一种众所周知的近似求解数理边值问题的数值技术,在结构力学中的应用已经很成熟,这里不再具体赘述。
Kmαm=Fm (6)
其中m为第几个单元,Km为第m个单元的刚度矩阵,αm为第m个单元的位移向量,Fm为第m个单元的外载荷量。具体表达式如下
Km=∫∫∫ΩNTATDANdV (7)
N为体插值基函数的矩阵形式
N=[N1 N2 … Ni … Nl] (12)
(13)式中Ni为体插值基函数,下标表示第几个插值点。对所有的四面体单元进行编号,同时对四面体单元内的插值点进行编号,最后去除重复的插值点,得到一组全局编号,该编号的个数即为整体系统的自由度,然后通过有限元系统装配以得到如下整体结构有限元方程,具体装配过程是一种公知的过程,这里不再描述。
Kα=F (14)
其中K为弹性体的刚度矩阵,α为结构位移向量,F为外载荷量。
其中nf为系统总自由度。
其中
式(18)、(19)、(20)、(21)是在总体坐标系中定义,为方便引入接触条件,需要将其转换到局部坐标系中,即
其中T是两种坐标系之间的转换矩阵,T是矩阵转置符号
在黏结接触状态下局部坐标中一个接触点对的接触力引起的等效节点力向量为
进一步以得到第h个联合面网格单元整体坐标系下接触力等效节点力向量
或者写为
对所有联合面网格单元的接触节点计算,并按照联合面网格所属体单元的编号,以及面插值点在体单元的位置组装到式(14),则可得到系统的运动方程,即
(K+Kc)α=F (27)
其中Kc为弹性体的接触刚度矩阵。
F.引入结构的惯性力,得到结构的自由振动有限元广义本征方程。
当研究结构振动问题时,上述E步骤的α位移向量为时间的函数,我们重新定义时间函数的结构位移向量
其中ui(t),vi(t),wi(t)表示是时间函数的三个位移分量。则根据E步骤讨论得到的有限元方程(27),引入物体的惯性力得到
其中
M=∫∫∫ΩρNTNdΩ (30)
当物体自由振动时,此时F(t)=0方程(29)退化为
其振动形式叫做自由振动,该方程有解的形式
消去ejωt后,有
该方程有非零解的条件是
|(K+Kc)-ω2M|=0 (35)
这就是接触问题模态分析广义本征方程。
H.对G步骤获得的特征值和对应特征向量进行后处理获得振动模态频率和对应振动振型。
对G步骤获得特征值λj进行处理,对应的振动模态频率为
Claims (3)
1.一种三维力学有限元模态分析中的接触分析方法,包括以下步骤:
A.将目标电子器件结构进行建模,引入位移边界条件或者应力边界条件建立对应的几何结构模型;
B.根据电子器件各部件之间的相互关系,生成几何结构各部件之间的接触对关系;
C.对所建几何结构模型采用四面体网格进行剖分,将连续的几何结构空间转化为离散空间;
采用四面体网格剖分仿真区域,剖分后的仿真区域被人为分割为多个三维四面体网格,从而将连续的几何结构空间转化为离散的网格空间;此时,接触面由一系列离散的面网格组成;
D.根据仿真区域接触面的面网格生成联合面网格;
E.利用有限元法,将平衡微分方程、几何方程和物理方程等效的结构力学边值问题在步骤C建立的网格空间进行离散,结合步骤D建立的联合面网格,建立考虑接触问题的目标电子器件的有限元方程;
F.引入几何结构的惯性力,得到其自由振动有限元广义本征方程;
H.对G步骤获得的特征值和对应特征向量进行后处理获得振动模态频率和对应振动振型。
2.如权利要求1所述三维力学有限元模态分析中的接触分析方法,其特征在于:
所述步骤D中联合面网格的步骤具体为:
1)确定接触面网格的位置和所属四面体单元;
根据步骤C生成的四面体网格信息和步骤B生成的接触对关系,确定每一个接触对上的面网格以及其所属四面体单元;
2)复制接触单元面网格;
对步骤1)中的步骤所确定的接触对上的面网格进行共形复制,在接触面上形成一套共形网格,以分别携带不同的材料属性。
3.如权利要求1所述三维力学有限元模态分析中的接触分析方法,其特征在于:
所述步骤E具体如下:
对于空间边值问题,在结构即弹性体的内部考虑静力学、几何学、物理学三方面条件,分别建立三套方程;并在给定约束或面力的边界上,建立位移边界条件或应力边界条件;
平衡微分方程:
几何方程:
物理方程:
位移边界条件:
应力边界条件:
上述(1)(2)(3)(4)(5)式中,σx,σy,σz,τxy=τyx,τyz=τzy,τzx=τxz表示求解区域中6个应力分量,εx,εy,εz,γxy,γyz,γzx表示求解区域中6个形变应力分量,u,v,w表示求解区域中3个位移分量;E是求解区域中结构的杨氏弹性模量,μ是求解区域中结构的泊松比,Su表示位移边界面,Sσ表示应力边界面;a=cos(n′,x),b=cos(n′,y),c=cos(n′,z),表示应力边界面Sσ上的方向余弦,其中n′为应力边界面Sσ的外法线,x,y,z为应力边界面Sσ上三个方向的坐标值;为位移边界面Su上的位移值,fx,fy,fz为求解区域内结构受到的各个方向的体力,为应力边界面Sσ上受到的各个方向的面力;
经过有限元推导得到如下每一个四面体单元中的有限元方程;
Kmαm=Fm (6)
其中m为第几个单元,Km为第m个单元的刚度矩阵,αm为第m个单元的位移向量,Fm为第m个单元的外载荷量;
Km=∫∫∫ΩNTATDANdV (7)
N为体插值基函数的矩阵形式
N=[N1 N2…Ni…Nl] (12)
(13)式中Ni为体插值基函数,下标表示第几个插值点,对所有的四面体单元进行编号,同时对四面体单元内的插值点进行编号,最后去除重复的插值点,得到一组全局编号,该编号的个数即为整体系统的自由度,然后通过有限元系统装配以得到如下整体结构有限元方程;
Kα=F (14)
其中K为弹性体的刚度矩阵,α为结构位移向量,F为外载荷量;
其中nf为系统总自由度;
其中
式(18)、(19)、(20)、(21)是在总体坐标系中定义,为方便引入接触条件,需要将其转换到局部坐标系中,即
其中T是两种坐标系之间的转换矩阵,T是矩阵转置符号
在黏结接触状态下局部坐标中一个接触点对的接触力引起的等效节点力向量为
进一步以得到第h个联合面网格单元整体坐标系下接触力等效节点力向量
或者写为
对所有联合面网格单元的接触节点计算,并按照联合面网格所属体单元的编号,以及面插值点在体单元的位置组装到式(14),则可得到系统的运动方程,即
(K+Kc)α=F (27)
其中Kc为弹性体的接触刚度矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710724845.9A CN107515982B (zh) | 2017-08-22 | 2017-08-22 | 一种三维力学有限元模态分析中的接触分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710724845.9A CN107515982B (zh) | 2017-08-22 | 2017-08-22 | 一种三维力学有限元模态分析中的接触分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107515982A CN107515982A (zh) | 2017-12-26 |
CN107515982B true CN107515982B (zh) | 2020-08-11 |
Family
ID=60723336
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710724845.9A Active CN107515982B (zh) | 2017-08-22 | 2017-08-22 | 一种三维力学有限元模态分析中的接触分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107515982B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108415267B (zh) * | 2018-03-02 | 2021-03-26 | 南京理工大学 | 基于apdl语言用于机床等效结合面快速定义的方法 |
CN110321571B (zh) * | 2018-03-29 | 2021-09-28 | 中国科学院沈阳自动化研究所 | 一种蜂窝板壳结构的力学参数数值提取方法 |
CN108694280B (zh) * | 2018-05-14 | 2021-02-26 | 电子科技大学 | 基于新型杂交应力四面体单元的冲击响应仿真模拟方法 |
CN108595893B (zh) * | 2018-05-16 | 2021-06-01 | 电子科技大学 | 一种基于三层预处理子的三维力学模态仿真模拟方法 |
CN108984829B (zh) * | 2018-06-06 | 2021-02-02 | 中国农业大学 | 堆石混凝土堆石体堆积过程的计算方法和系统 |
CN109033501A (zh) * | 2018-06-08 | 2018-12-18 | 昆明理工大学 | 一种包含刚体运动的求解域动态演化的几何模型建模方法 |
CN108875195B (zh) * | 2018-06-13 | 2022-03-15 | 电子科技大学 | 一种考虑接触的三维力学随机振动仿真模拟方法 |
CN109241596A (zh) * | 2018-08-24 | 2019-01-18 | 南京航空航天大学 | 一种三维层合结构动力学分析方法 |
CN109308382B (zh) * | 2018-08-31 | 2023-10-24 | 北京航空航天大学 | 一种针对奇异点附近接触应力的分析方法 |
CN109359360B (zh) * | 2018-09-30 | 2022-11-11 | 国家超级计算天津中心 | 一种基于局部特征的结构应力处理方法 |
CN109636912B (zh) * | 2018-11-27 | 2020-06-26 | 中国地质大学(武汉) | 应用于三维声呐图像重构的四面体剖分有限元插值方法 |
CN109726433B (zh) * | 2018-11-30 | 2022-10-14 | 电子科技大学 | 基于曲面边界条件的三维无粘低速绕流的数值模拟方法 |
CN109684723B (zh) * | 2018-12-24 | 2022-05-10 | 哈尔滨工程大学 | 一种二维结构内部声学性能分析方法 |
CN109783946A (zh) * | 2019-01-21 | 2019-05-21 | 河北工业大学 | 一种声子晶体带隙仿真的节点积分算法 |
CN111209675B (zh) * | 2020-01-10 | 2023-04-28 | 南方电网科学研究院有限责任公司 | 电力电子器件的仿真方法、装置、终端设备及存储介质 |
CN111832205B (zh) * | 2020-07-02 | 2022-12-09 | 西安交通大学 | 一种适应波群传递界面动态变化的能量有限元分析方法 |
CN112329303B (zh) * | 2020-11-03 | 2022-09-06 | 西安电子科技大学 | 基于有限元区域分解的阵列天线电磁特性求解方法 |
CN112665774A (zh) * | 2021-01-11 | 2021-04-16 | 哈尔滨工业大学 | 一种基于dic的液压管路接头侧向力测试方法 |
CN113032889B (zh) * | 2021-05-31 | 2022-03-01 | 北京盈建科软件股份有限公司 | 由基础结构和上部结构拼接成联合建筑模型的方法及其装置 |
CN114896710B (zh) * | 2022-07-15 | 2022-10-28 | 中国飞机强度研究所 | 用于航空结构冲击分析的接触碰撞分析方法 |
CN116306326B (zh) * | 2023-05-25 | 2023-09-15 | 南方医科大学珠江医院 | 一种关节接触力学仿真模型建立方法、装置及电子设备 |
CN118171547B (zh) * | 2024-05-16 | 2024-08-13 | 深圳十沣科技有限公司 | 自动化网格生成方法、装置、设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930071A (zh) * | 2012-08-29 | 2013-02-13 | 电子科技大学 | 基于非匹配网格的周期结构的三维电磁场仿真模拟方法 |
CN106354954A (zh) * | 2016-08-31 | 2017-01-25 | 电子科技大学 | 一种基于叠层基函数的三维力学模态仿真模拟方法 |
-
2017
- 2017-08-22 CN CN201710724845.9A patent/CN107515982B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102930071A (zh) * | 2012-08-29 | 2013-02-13 | 电子科技大学 | 基于非匹配网格的周期结构的三维电磁场仿真模拟方法 |
CN106354954A (zh) * | 2016-08-31 | 2017-01-25 | 电子科技大学 | 一种基于叠层基函数的三维力学模态仿真模拟方法 |
Non-Patent Citations (2)
Title |
---|
Structural damage detection based on static and modal analysis;P. Hajela 等;《AIAA JOURNAL》;20020517;第28卷(第6期);第1110-1115页 * |
基于参数化的齿轮传动接触有限元分析;刘晖;《中国优秀博硕士学位论文全文数据库—工程科技Ⅱ辑》;20060215(第2期);第12-40页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107515982A (zh) | 2017-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107515982B (zh) | 一种三维力学有限元模态分析中的接触分析方法 | |
Rao | The Finite Element Method in Engineering: Pergamon International Library of Science, Technology, Engineering and Social Studies | |
Fried | Numerical solution of differential equations | |
Park et al. | A variational principle for the formulation of partitioned structural systems | |
Atluri et al. | New concepts in meshless methods | |
Duarte et al. | Generalized finite element methods for three-dimensional structural mechanics problems | |
Nguyen-Thoi et al. | A cell-based smoothed three-node Mindlin plate element (CS-MIN3) for static and free vibration analyses of plates | |
CN108875195B (zh) | 一种考虑接触的三维力学随机振动仿真模拟方法 | |
Farhat et al. | Robust and provably second‐order explicit–explicit and implicit–explicit staggered time‐integrators for highly non‐linear compressible fluid–structure interaction problems | |
Johnson et al. | Parallel computation of incompressible flows with complex geometries | |
Lin | A mesh-based partition of unity method for discontinuity modeling | |
KR20130098306A (ko) | 계산용 데이터 생성 장치, 계산용 데이터 생성 방법 및 계산용 데이터 생성 프로그램 | |
CN106354954B (zh) | 一种基于叠层基函数的三维力学模态仿真模拟方法 | |
Bos et al. | Formula manipulation in the bond graph modelling and simulation of large mechanical systems | |
CN106295028A (zh) | 一种局部结构动力学建模方法及装置 | |
He et al. | Augmented numerical manifold method with implementation of flat-top partition of unity | |
Craig et al. | William R. Sherman | |
Sukumar et al. | Meshless methods and partition of unity finite elements | |
Marjanović et al. | Object-oriented framework for 3D bending and free vibration analysis of multilayer plates: Application to cross-laminated timber and soft-core sandwich panels | |
US7773814B2 (en) | Multi-scale analysis device | |
Montano et al. | Numerical simulation of tethered buoy dynamics using mixed finite elements | |
Frisch | A vector-dyadic development of the equations of motion for N-coupled flexible bodies and point masses | |
de la Fuente | An efficient procedure to obtain exact solutions in random vibration analysis of linear structures | |
Winkler et al. | A projection-based approach for the derivation of the floating frame of reference formulation for multibody systems | |
Han et al. | A Truly-Meshless Galerkin Method, through the MLPG" Mixed" Approach |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 |