CN114741938B - 一种斜拉桥调索非弹性收缩量智能监测方法 - Google Patents
一种斜拉桥调索非弹性收缩量智能监测方法 Download PDFInfo
- Publication number
- CN114741938B CN114741938B CN202210658962.0A CN202210658962A CN114741938B CN 114741938 B CN114741938 B CN 114741938B CN 202210658962 A CN202210658962 A CN 202210658962A CN 114741938 B CN114741938 B CN 114741938B
- Authority
- CN
- China
- Prior art keywords
- cable
- stay
- adjusting
- stayed bridge
- shrinkage
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000012544 monitoring process Methods 0.000 title claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 41
- 238000006073 displacement reaction Methods 0.000 claims abstract description 20
- 239000013598 vector Substances 0.000 claims description 64
- 230000001105 regulatory effect Effects 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 235000017166 Bambusa arundinacea Nutrition 0.000 description 1
- 235000017491 Bambusa tulda Nutrition 0.000 description 1
- 241001330002 Bambuseae Species 0.000 description 1
- 235000015334 Phyllostachys viridis Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000011425 bamboo Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
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
- 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
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- Bridges Or Land Bridges (AREA)
Abstract
本发明公开了一种斜拉桥调索非弹性收缩量智能监测方法,包括步骤一、根据斜拉桥的结构信息建立有限元模型,并根据有限元模型求解得到各斜拉索单元的平衡方程;步骤二、将非弹性约束变量施加于有限元模型中得到斜拉索状态模型;步骤三、根据斜拉索状态模型计算斜拉索张拉力对斜拉桥各控制节点位移的影响矩阵;步骤四、根据影响矩阵和所述斜拉索状态模型建立矩阵方程,实时获取斜拉桥各控制节点的位移量,并根据所述矩阵方程,依次求解得到斜拉索张拉力;步骤五、根据斜拉索张拉力和目标索力,基于所述影响矩阵求解得到调索收缩量,实现了对斜拉桥调索非弹性收缩量智能监测,对斜拉桥的安全稳定性监测以及进一步制定索力调整方案具有现实意义。
Description
技术领域
本发明涉及斜拉桥安全监测技术领域,特别涉及一种斜拉桥调索非弹性收缩量智能监测方法。
背景技术
随着社会经济建设的迅猛发展,以及土木工程理论与技术的不断成熟,交通基础设施的建设日新月异,桥梁作为跨越障碍的主要结构,各类桥型如雨后春笋般层出不穷,斜拉桥作为其中一类因其跨越能力大而备受关注与重视,并得到广泛的应用,随着新材料(如碳纤维增强聚合物索) 和新施工工艺的发展,斜拉桥的跨越能力不断得到提高,从而在桥梁建设中更具有竞争力。但是,由于斜拉桥跨度大,桥面梁和斜拉索相对于桥塔为柔性结构,在交通载荷、强风和地震等复杂环境载荷作用下,其动力学行为相当复杂,在斜拉桥投入使用及运营中,需要实时监测斜拉桥的安全性能,并及时进行维修或调整。
斜拉桥一般用斜拉索作为主梁的支撑,主梁的受力与多跨连续梁的受力状态相似,斜拉索可以理解为对斜拉桥施加了一个体外预应力。斜拉桥在力学上属于为高次超静定结构,影响其成桥内力的主要因素有恒载分布、索力和边界支承反力等。斜拉桥经营过程中可以通过索力调整来改变成桥内力,从而使斜拉桥结构的线形、索力、结构受力等达到设计成桥理想状态。调整斜拉桥索力对斜拉桥的内力状态影响较大,故实时检测斜拉桥的调索收缩量以判断斜拉桥的安全稳定性,及其后续调索方案的制定和实施具有重要的现实意义。
发明内容
本发明提供了一种斜拉桥调索非弹性收缩量智能监测方法,通过建立斜拉桥的有限元模型,采用影响矩阵法,基于实时获取的斜拉桥各控制节点的位移量,求解得到斜拉索张拉力,进而得到斜拉索索力调整值及调索收缩量。
本发明的技术方案为:
一种斜拉桥调索非弹性收缩量智能监测方法,包括:
步骤一、根据斜拉桥的结构信息建立有限元模型,并根据有限元模型求解得到各斜拉索单元的平衡方程;
步骤二、将非弹性约束变量施加于有限元模型中得到斜拉索状态模型;
步骤三、根据斜拉索状态模型计算斜拉索张拉力对斜拉桥各控制节点位移的影响矩阵;
步骤四、根据影响矩阵和斜拉索状态模型建立矩阵方程,实时获取斜拉桥各控制节点的向变量,并根据所述矩阵方程,依次求解得到斜拉索张拉力;
步骤五、根据斜拉索张拉力和目标索力,基于影响矩阵求解得到调索收缩量。
优选的是,结构信息包括材料、弹性模量和泊松比。
优选的是,斜拉索平衡方程为:
其中,表示斜拉索单元的弹性模量,表示斜拉索横截面面积,表示建模长度,表示斜拉索角度余弦,表示斜拉索角度正弦,分别表示一组斜拉索的轴向位移,分别表示一组斜拉索的纵向位移,表示斜拉索非弹性收缩量,分别表示一组斜拉索的轴向荷载向量,分别表示一组斜拉索的轴向荷载向量,分别表示一组斜拉索的纵向荷载向量,表示斜拉索拉力。
优选的是,非弹性约束变量包括斜拉索位移和垂度、桥面线形、主缆线形、桥塔偏位和索鞍位置。
优选的是,步骤三包括:
自定义受调向量和施调向量;
将施调向量施加于斜拉索状态模型,得到受调向量的应变量,以生成施调向量对所述受调向量的影响矩阵。
优选的是,受调向量的受调元为各关键控制节点的向变量,向变量包括位移、截面应力、内力为受调向量中的几种;施调向量的施调元为各斜拉索张拉力。
优选的是,影响矩阵的生成过程为:
计算调值阶段整体刚度矩阵;
令施调向量中的任意施调元为单位调值量,得到相应的结构荷载阵列;
计算荷载阵列对应的任意施调元与各受调元的影响向量;
遍历施调元,得到影响矩阵。
优选的是,矩阵方程为:
其中,表示第个施调元对第个受调元的竖向影响向量,表示第个施调元对第个受调元的横向影响向量,表示第个施调元对第个受调元的旋转影响向量,表示斜拉索张拉力,表示恒载对第个受调元的竖向影响向量,表示恒载对第个受调元的横向影响向量,表示恒载对第个受调元的竖向影响向量,表示第个受调元的竖向控制值,表示第个受调元的横向控制值,表示第个受调元的旋转控制值。
优选的是,目标索力为斜拉桥各斜拉索的成桥索力。
优选的是,步骤五包括:
将所述斜拉索张拉力和目标索力进行对比,计算各斜拉索索力调整值;
将各斜拉索的所述索力调整值代入所述影响矩阵,计算得到各斜拉索的调索收缩量。
本发明的有益效果是:
本发明通过建立斜拉桥的有限元模型,采用影响矩阵法,基于实时获取的斜拉桥各控制节点的位移量,求解得到斜拉索张拉力,进而得到斜拉索索力调整值及调索收缩量,通过对调索收缩量的实时监测有利于及时对斜拉桥的安全稳定性做出预判、制定调索方案。
附图说明
图1为本发明提供的实施例中一种斜拉桥调索非弹性收缩量智能监测方法流程图。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式,熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明的是,在本发明的描述中,术语“中”、“上”、“下”、“横”、“内”等指示的方向或位置关系的术语是基于附图所示的方向或位置关系,这仅仅是为了便于描述,而不是指示或暗示装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性。
此外,还需要说明的是,在本发明的描述中,除非另有明确的规定和限定,术语“设置”、“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域技术人员而言,可根据具体情况理解上述术语在本发明中的具体含义。
如图1所示,一种斜拉桥调索非弹性收缩量智能监测方法,包括:
S110、根据斜拉桥的结构信息建立有限元模型,并根据有限元模型求解得到各斜拉索单元的平衡方程。
采用通用有限元软件ANSYS进行建模,通过输入点,定义插值曲线,定义函数等多种方式建立有限元单元,通过输入结构信息,定义材料等方式进行模拟。
具体的,结构信息包括材料、弹性模量和泊松比。
斜拉索平衡方程为:
其中,表示斜拉索单元的弹性模量,表示斜拉索横截面面积,表示建模长度,表示斜拉索角度余弦,表示斜拉索角度正弦,分别表示一组斜拉索的轴向位移,分别表示一组斜拉索的纵向位移,表示斜拉索非弹性收缩量,分别表示一组斜拉索的轴向荷载向量,分别表示一组斜拉索的轴向荷载向量,分别表示一组斜拉索的纵向荷载向量,表示斜拉索拉力。
S120、将非弹性约束变量施加于有限元模型中得到斜拉索状态模型。
斜拉桥在力学上属于为高次超静定结构,约束边界和荷载是斜拉桥的斜拉索状态的主要影响因素,因此引入非弹性约束变量为约束边界和/或荷载,并加载到有限元模型中,可以对斜拉桥进行仿真模拟,得到斜拉桥拉索状态模型。
其中,非弹性约束变量包括斜拉索位移和垂度、桥面线形、主缆线形、桥塔偏位和索鞍位置。
S130、根据斜拉索状态模型计算斜拉索张拉力对斜拉桥各控制节点位移的影响矩阵。
自定义受调向量和施调向量;
一个优选的实施例是,受调向量的受调元为各关键控制节点的向变量,向变量包括位移、截面应力、内力为受调向量中的几种。
一个优选的实施例是,施调向量的施调元为各斜拉索张拉力。
将施调向量施加于斜拉索状态模型,得到受调向量的应变量,以生成施调向量对所述受调向量的影响矩阵。
具体的,影响矩阵的生成过程为:
计算调值阶段整体刚度矩阵;
令施调向量中的任意施调元为单位调值量,得到相应的结构荷载阵列;
计算荷载阵列对应的任意施调元与各受调元的影响向量;
遍历施调元,得到影响矩阵。
S140、根据影响矩阵和斜拉索状态模型建立矩阵方程,实时获取斜拉桥各控制节点的向变量,并根据所述矩阵方程,依次求解得到斜拉索张拉力。
具体的,矩阵方程为:
其中,表示第个施调元对第个受调元的竖向影响向量,表示第个施调元对第个受调元的横向影响向量,表示第个施调元对第个受调元的旋转影响向量,表示斜拉索张拉力,表示恒载对第个受调元的竖向影响向量,表示恒载对第个受调元的横向影响向量,表示恒载对第个受调元的竖向影响向量,表示第个受调元的竖向控制值,表示第个受调元的横向控制值,表示第个受调元的旋转控制值。
具体的,斜拉桥各控制节点的向变量通过在各关键控制节点设置应变仪、锚索计、索力传感器或弦振式索力仪等实时采集。
S150、根据斜拉索张拉力和目标索力,基于影响矩阵求解得到调索收缩量。
具体的,目标索力为斜拉桥各斜拉索的成桥索力。
首先,将斜拉索张拉力和目标索力进行对比,计算各斜拉索索力调整值;
将各斜拉索的索力调整值代入所述影响矩阵,计算得到各斜拉索的调索收缩量。
本发明通过建立斜拉桥的有限元模型,采用影响矩阵法,基于实时获取的斜拉桥各控制节点的位移量,求解得到斜拉索张拉力,进而得到斜拉索索力调整值及调索收缩量,通过对调索收缩量的实时监测有利于及时对斜拉桥的安全稳定性做出预判、制定调索方案。
以上内容仅是本发明的实施例,方案中公知的具体结构及特性等常识在此未作过多描述,本领域技术人员容易理解的是,本发明的保护范围显然不局限于这些具体实施方式。在不脱离本发明的前提下,还可以作出若干变形和改进,这些也应该视为本发明的保护范围,这些都不会影响本发明实施的效果和专利的实用性。
Claims (10)
1.一种斜拉桥调索非弹性收缩量智能监测方法,其特征在于,包括:
步骤一、根据斜拉桥的结构信息建立有限元模型,并根据所述有限元模型求解得到各斜拉索单元的平衡方程;
步骤二、将非弹性约束变量施加于所述有限元模型中得到斜拉索状态模型;
步骤三、根据所述斜拉索状态模型计算斜拉索张拉力对斜拉桥各控制节点位移的影响矩阵;
步骤四、根据所述影响矩阵和所述斜拉索状态模型建立矩阵方程,实时获取斜拉桥各控制节点的向变量,并根据所述矩阵方程,依次求解得到斜拉索张拉力;
步骤五、根据所述斜拉索张拉力和目标索力,基于所述影响矩阵求解得到调索收缩量。
2.如权利要求1所述的斜拉桥调索非弹性收缩量智能监测方法,其特征在于,所述结构信息包括材料、弹性模量和泊松比。
4.如权利要求3所述的斜拉桥调索非弹性收缩量智能监测方法,其特征在于,所述非弹性约束变量包括斜拉索位移和垂度、桥面线形、主缆线形、桥塔偏位和索鞍位置。
5.如权利要求4所述的斜拉桥调索非弹性收缩量智能监测方法,其特征在于,所述步骤三包括:
自定义受调向量和施调向量;
将所述施调向量施加于所述斜拉索状态模型,得到受调向量的应变量,以生成所述施调向量对所述受调向量的影响矩阵。
6.如权利要求5所述的斜拉桥调索非弹性收缩量智能监测方法,其特征在于,所述受调向量的受调元为各关键控制节点的向变量,所述向变量包括位移、截面应力、内力为受调向量中的几种;所述施调向量的施调元为各斜拉索张拉力。
7.如权利要求6所述的斜拉桥调索非弹性收缩量智能监测方法,其特征在于,所述影响矩阵的生成过程为:
计算调值阶段整体刚度矩阵;
令施调向量中的任意施调元为单位调值量,得到相应的结构荷载阵列;
计算所述荷载阵列对应的所述任意施调元与各受调元的影响向量;
遍历所述施调元,得到影响矩阵。
9.如权利要求8所述的斜拉桥调索非弹性收缩量智能监测方法,其特征在于,所述目标索力为斜拉桥各斜拉索的成桥索力。
10.如权利要求9所述的斜拉桥调索非弹性收缩量智能监测方法,其特征在于,所述步骤五包括:
将所述斜拉索张拉力和目标索力进行对比,计算各斜拉索索力调整值;
将各斜拉索的所述索力调整值代入所述影响矩阵,计算得到各斜拉索的调索收缩量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210658962.0A CN114741938B (zh) | 2022-06-13 | 2022-06-13 | 一种斜拉桥调索非弹性收缩量智能监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210658962.0A CN114741938B (zh) | 2022-06-13 | 2022-06-13 | 一种斜拉桥调索非弹性收缩量智能监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114741938A CN114741938A (zh) | 2022-07-12 |
CN114741938B true CN114741938B (zh) | 2022-08-19 |
Family
ID=82287476
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210658962.0A Active CN114741938B (zh) | 2022-06-13 | 2022-06-13 | 一种斜拉桥调索非弹性收缩量智能监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114741938B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115434247B (zh) * | 2022-08-30 | 2024-09-27 | 中铁大桥勘测设计院集团有限公司 | 一种斜拉索张拉两阶段控制方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08326013A (ja) * | 1995-06-02 | 1996-12-10 | Kajima Corp | 斜張橋斜材ケーブルの制振方法およびその装置 |
CN103870702A (zh) * | 2014-03-26 | 2014-06-18 | 重庆交通大学 | 一种含非弹性收缩量的预应力筋单元 |
CN104713673A (zh) * | 2015-03-11 | 2015-06-17 | 哈尔滨工业大学 | 基于扩展卡尔曼滤波器的拉索时变索力历程识别方法 |
WO2018001147A1 (zh) * | 2016-06-29 | 2018-01-04 | 深圳市智能机器人研究院 | 一种基于优化张紧弦模型的桥索监测方法及系统 |
CN111475883A (zh) * | 2020-04-16 | 2020-07-31 | 重庆交通大学 | 一种大跨度钢管混凝土拱桥拱肋线形控制方法 |
CN111611644A (zh) * | 2020-05-28 | 2020-09-01 | 广西新发展交通集团有限公司 | 施工过程索力最优的斜拉桥斜拉索张拉方法 |
-
2022
- 2022-06-13 CN CN202210658962.0A patent/CN114741938B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08326013A (ja) * | 1995-06-02 | 1996-12-10 | Kajima Corp | 斜張橋斜材ケーブルの制振方法およびその装置 |
CN103870702A (zh) * | 2014-03-26 | 2014-06-18 | 重庆交通大学 | 一种含非弹性收缩量的预应力筋单元 |
CN104713673A (zh) * | 2015-03-11 | 2015-06-17 | 哈尔滨工业大学 | 基于扩展卡尔曼滤波器的拉索时变索力历程识别方法 |
WO2018001147A1 (zh) * | 2016-06-29 | 2018-01-04 | 深圳市智能机器人研究院 | 一种基于优化张紧弦模型的桥索监测方法及系统 |
CN111475883A (zh) * | 2020-04-16 | 2020-07-31 | 重庆交通大学 | 一种大跨度钢管混凝土拱桥拱肋线形控制方法 |
CN111611644A (zh) * | 2020-05-28 | 2020-09-01 | 广西新发展交通集团有限公司 | 施工过程索力最优的斜拉桥斜拉索张拉方法 |
Non-Patent Citations (1)
Title |
---|
一种基于指定应力的斜拉桥成桥索力调整方法;王家林等;《公路交通科技》;20200715(第07期);75-81 * |
Also Published As
Publication number | Publication date |
---|---|
CN114741938A (zh) | 2022-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fleming et al. | Dynamic behaviour of a cable‐stayed bridge | |
CN111695188B (zh) | 一种地锚式悬索桥动力特性的快速精细分析方法 | |
CN114741938B (zh) | 一种斜拉桥调索非弹性收缩量智能监测方法 | |
Wahba et al. | Evaluation of non-linear analysis of guyed antenna towers | |
CN111783201A (zh) | 一种三跨自锚式悬索桥动力特性的快速分析方法 | |
CN103696356A (zh) | 一种设有双排支座体系的多塔斜拉桥 | |
CN109492293B (zh) | 一种倾斜悬索的静、动力作用刚度模型的构建方法 | |
CN105803941A (zh) | 一种斜拉桥双索同时张拉无应力状态施工控制方法 | |
Fei et al. | Vertical vibrations of suspension bridges: a review and a new method | |
CN110020498A (zh) | 曲线钢-混凝土组合箱型梁一维有限元模型的构建方法 | |
CN103913338A (zh) | 混合监测问题索载荷广义位移识别方法 | |
CN117271949B (zh) | 考虑弹性边界和大垂度影响的悬索自振分析方法及系统 | |
CN113468632B (zh) | 一种偏心活载作用下的悬索桥全桥响应的确定方法 | |
CN103852289A (zh) | 空间坐标监测问题索载荷广义位移识别方法 | |
CN105631169B (zh) | 一种索杆结构初始应变确定方法 | |
CN112211114A (zh) | 一种塔梁固结体系斜拉桥的悬臂拼装施工方法 | |
Ko et al. | Modal analysis of suspension bridge deck units in erection stage | |
CN111666626A (zh) | 索结构形态分析中的预应力水平调节方法 | |
CN103852329A (zh) | 混合监测受损索载荷广义位移识别方法 | |
CN115906333A (zh) | 一种桁架结构的几何非线性等效板动力学建模与响应分析方法 | |
CN109977565A (zh) | 一种基于fea的斜拉桥动力特性分析方法 | |
CN103852308A (zh) | 混合监测问题索载荷线位移识别方法 | |
CN103868733A (zh) | 线位移角度监测问题索载荷识别方法 | |
Raman et al. | Large displacement analysis of guyed towers | |
Xiaowei et al. | Numerical simulation on wind-induced responses of steel guyed transmission tower-line coupling systems |
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 |