CN106649954B - 一种基于扩展pbee2理论框架的地震易损性分析方法 - Google Patents

一种基于扩展pbee2理论框架的地震易损性分析方法 Download PDF

Info

Publication number
CN106649954B
CN106649954B CN201610877011.7A CN201610877011A CN106649954B CN 106649954 B CN106649954 B CN 106649954B CN 201610877011 A CN201610877011 A CN 201610877011A CN 106649954 B CN106649954 B CN 106649954B
Authority
CN
China
Prior art keywords
earthquake
vulnerability
seismic
model
establishing
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
CN201610877011.7A
Other languages
English (en)
Other versions
CN106649954A (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.)
Huatian Engineering and Technology Corp MCC
Original Assignee
Huatian Engineering and Technology Corp MCC
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 Huatian Engineering and Technology Corp MCC filed Critical Huatian Engineering and Technology Corp MCC
Priority to CN201610877011.7A priority Critical patent/CN106649954B/zh
Publication of CN106649954A publication Critical patent/CN106649954A/zh
Application granted granted Critical
Publication of CN106649954B publication Critical patent/CN106649954B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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]

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)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于扩展PBEE2理论框架的地震易损性分析方法，属于地震风险及震害预测分析技术领域。本发明的步骤为：(1)建立典型结构空间：(1.1)选取抗震设防烈度、高度和宽度比作为影响结构抗震性能的主导参数；(1.2)基于逻辑树原理，建立典型结构空间；(1.3)采用PKPM设计各典型结构；(2)建立主导参数与其概率地震需求模型参数间的相关性矩阵；(3)基于对数正态分布的一类结构的概率地震易损性函数。该方法实现了一类结构解析地震易损性函数的建立，操作方便可靠，为城市建筑群的区域震害快速预测提供了理论基础，更加适用于城市建筑群的区域震害预测。

Description

一种基于扩展PBEE2理论框架的地震易损性分析方法
技术领域
本发明涉及地震风险及震害预测分析技术领域，更具体地说，涉及一种基于扩展PBEE2理论框架的地震易损性分析方法。
背景技术
我国是一个地震多发的国家，地震给国家及人民带来了不可挽回的破坏与损失，如2008年的汶川地震造成了巨大人员伤亡与经济损失，由于准确地震预测的不可能性，需要对地震灾害进行风险分析，特别是对土木工程结构进行地震易损性分析就成为当前工程中主要的防灾减灾措施。我国对于地震易损性的研究多数停留于经验地震易损性，而由于震害资料的匮乏，仅有少数结构类型的建筑物，比如砌体结构，而RC结构的震害资料则尤为稀少，随着我国城镇化的发展，使得RC结构得以广泛应用。针对目前震害预测中缺乏RC框架、RC框架-剪力墙和RC剪力墙等诸多结构形式的经验及解析地震易损性函数，且考虑各单体建筑物间的变异性和城市建筑群中建筑物数量的海量性，提出一种适用于一类结构的解析地震易损性分析方法，成为城市建筑群大尺度震害预测的面临巨大挑战与亟待解决问题。
根据PBEE理论框架，对于单体结构的地震风险评估，其准确性取决于地震危险性模型、地震易损性模型(包括结构分析和破坏分析)和损失模型。因经验易损性曲线具有很大的地域局限性、数据有限且相关理论研究进展缓慢，因此对于解析易损性的研究成为提供单体风险评估精度的有效途径，比如采用更为精细化的结构数值模型和更为先进的非线性时程分析方法，均可一定程度提高预测结果的精度，降低预测结果的不确定性。
根据扩展的PBEE-2理论框架，相比于单体建筑物，影响一类建筑物地震风险评估精度的因素更为复杂。理论上可对评估区域中每个单体结构均进行的风险评估并求和，所得评估结果才为精确解，但由于城市建筑群中建筑物数量的海量性，因而基于所有单体逐个分析的研究思路并不可取。因此，区域性地震风险评估中，如何将无限多个建筑物简化为少量代表性的建筑物，而不失群体建筑物特性的方法，则成为了控制评估结果准确度的关键问题。
基于矩匹配(MM)抽样获取典型建筑物样本属性的方法，虽在一定程度上反映了样本的统计特性，但目前尚未给出样本模型的建立方法。显然，若要建立建筑物的数值模型，而仅已确定三个建筑物属性，比如层数、竖向不规则度和平面不规则度，模型的建立必然会存在很大的不确定，导致评估精度降低。
发明内容
1、要解决的问题
本发明的目的在于解决现有技术存在的不足和缺陷等问题，提供了一种基于扩展PBEE2理论框架的地震易损性分析方法。该方法给出了基于有限单体结构(典型结构)的一类结构解析地震易损性函数的建立方法，通过建立主导参数与地震易损性模型参数间的相关性矩阵，分析各结构设计参数对于地震易损性模型的影响规律，进而实现了一类结构解析地震易损性函数的建立，操作方便可靠，为城市建筑群的区域震害快速预测提供了理论基础，更加适用于城市建筑群的区域震害预测。
2、技术方案
为达到上述目的，本发明提供的技术方案为：
本发明的一种基于扩展PBEE2理论框架的地震易损性分析方法，其步骤为：
(1)建立典型结构空间：
(1.1)选取抗震设防烈度、高度和宽度比作为影响结构抗震性能的主导参数，其中层数、层高及跨度的选取，依照该类结构的大量数据统计分析得到；
(1.2)基于逻辑树原理，建立RC框架结构的典型结构空间；
(1.3)采用PKPM设计各典型结构，其恒载、活载取值依照该类结构的大量数据统计分析得到；
(2)建立主导参数与其概率地震需求模型参数间的相关性矩阵：
(2.1)基于OpenSEEs有限元软件平台,采用纤维模型和塑性铰模型，分别建立结构的非倒塌状态及倒塌状态数值模型；
(2.2)选取在未来地震下可能遭受的具有各种地震动特性的地震动记录，作为该类结构地震需求模型建立的输入地震动记录集，以考虑地震动的不确定性；
(2.3)采用基于位移和基于加速度的地震易损性模型，建立主导参数与其地震易损性模型参数间的相关性矩阵；
(3)基于对数正态分布的一类结构的概率地震易损性函数：
(3.1)根据典型结构设计参数与地震易损性函数参数间的数学模型关系，采用线性插值方法；计算其他跨度、层数非典型单体结构的地震易损性函数参数，从而建立一类结构中每个层数、跨度的结构地震易损性模型；
(3.2)以一类结构中不同典型结构的地震响应为总样本，按照单体结构的地震易损性建立方法，直接建立此类结构的地震易损性模型；
对于(3.2)中总不确定，按照下式考虑：
Figure GDA0002216427260000021
式中，βD为考虑地震动变异性、跨度变化的对数标准差，βM为建模不确定性。
作为本发明的进一步改进，步骤(2.3)中基于典型结构建立主导参数与其地震易损性模型参数间的相关性矩阵。
作为本发明的进一步改进，步骤(3)中给出了两种方法实现一类结构地震易损性模型的建立。
3、有益效果
相比于现有技术，本发明的有益效果为：
本发明的一种基于扩展PBEE2理论框架的地震易损性分析方法，体现了建筑物设计参数变化对于其地震易损性模型参数的影响规律，能在充分考虑我国各类建筑物诸主要设计参数变异性并兼顾有限数量典型结构计算量的同时，提出了基于有限单体(典型结构)到一类结构地震易损性模型的建立方法，搭建起了“单体”到“群体”建筑物(从特例到一般)地震易损性分析间的桥梁，将单体结构的精细化建模分析方法，应用于一类结构中，对于传统的基于经验的地震易损矩阵及基于能力谱法的简化方法，进行了完善与发展，为城市建筑群的区域震害快速预测提供了理论基础，更加适用于城市建筑群的区域震害预测。
附图说明
图1为本发明的框架流程图；
图2为典型RC框架结构逐层细化与典型结构空间划分示意图；
图3为典型RC框架结构层数、层高及跨度示意图；
图4为典型结构地震易损性评估流程图；
图5为F-5-Ⅶ-6m概率需求模型、IDA曲线及倒塌地震易损性曲线图，(a)为F-5-Ⅶ-6m概率需求模型图，(b)为IDA曲线图，(c)为倒塌地震易损性曲线图；
图6为对数标准差βD|IM与层数、跨度、设计基本加速度值和基本周期间的相关性图，(a)为对数标准差βD|IM与层数相关性图，(b)为对数标准差βD|IM与跨度相关性图，(c)为对数标准差βD|IM与设计基本加速度值相关性图，(d)为对数标准差βD|IM与基本周期间的相关性图；
图7为对数线性回归一次项系数β1与层数、跨度、设计基本加速度值和基本周期间的相关性图，(a)为对数线性回归一次项系数β1与层数相关性图，(b)为对数线性回归一次项系数β1与跨度相关性图，(c)为对数线性回归一次项系数β1与设计基本加速度值相关性图，(d)为对数线性回归一次项系数β1与基本周期间的相关性图；
图8为对数线性回归常数项β0与层数、跨度、设计基本加速度值和基本周期间的相关性图，(a)为对数线性回归常数项β0与层数相关性图，(b)为对数线性回归常数项β0与跨度相关性图，(c)为对数线性回归常数项β0与设计基本加速度值相关性图，(d)为对数线性回归常数项β0与基本周期间的相关性图；
图9为各烈度RC框架结构地震易损性曲线图，(a)为6度设防地震易损性曲线图，(b)为7度设防地震易损性曲线图，(c)为7.5度设防地震易损性曲线图，(d)为8度设防地震易损性曲线图，(e)为8.5度设防地震易损性曲线图，(f)为9度设防地震易损性曲线图。
具体实施方式
下面结合具体实施例和附图对本发明进一步进行描述：
实施例1
如图1所示，本实施例的一种基于扩展PBEE2理论框架的地震易损性分析方法，考虑了一类结构的主要设计参数并基于逻辑树原理，建立了典型结构空间；基于OpenSEEs有限元软件平台和第二代基于性能的理论框架，建立了一类结构的主要设计参数与其概率地震需求模型参数间的相关性矩阵；以一类结构概率地震需求模型参数为统计分析对象，建立了适用于该类结构的概率地震需求参数及地震易损性模型，其具体步骤为：
1、RC框架结构的典型结构空间的建立：
(1.1)选取设防烈度、层数、跨度为主导参数，基于逻辑树原理建立RC框架结构的典型结构空间，如图2所示；
(1.2)确定典型结构层数、层高及跨度的取值，如图3所示；
(1.3)根据(1)和(2)，采用PKPM设计各典型结构；
2、主导参数与其概率地震需求模型参数间的相关性分析：
(2.1)基于OpenSEEs有限元软件平台，采用纤维模型和塑性铰模型，分别建立结构的非倒塌状态及倒塌状态数值模型；
(2.2)因地震动记录的选取原则与结构和场地特性无关，选取ATC 63推荐的未经调幅的22条远场地震动记录作为输入地震动记录集；
(2.3)采用基于位移和基于加速度的地震易损性模型，如图4所示，以8度设防、6m跨度、5层RC框架结构(F-5-Ⅶ-6m)为例，所得概率地震需求模型与IDA曲线，如图5所示；
(2.4)对于非倒塌状态，假定结构响应与地震动强度指标之间为对数线性关系，建立主导参数与其地震易损性模型参数间的相关性矩阵，如图6～8所示；
3、基于对数正态分布的一类结构的概率地震易损性曲线的建立分析：
对于基于位移的地震易损性模型，概率地震能力分析中分别采用0.002、0.005、0.015、0.04作为轻微、中等、严重和倒塌极限状态的层间位移角限值；
基于对数正态分布的地震易损性模型：
Figure GDA0002216427260000051
式中，基位移的地震易损性模型时
Figure GDA0002216427260000052
基于加速度的地震易损性模型时，mR为IDA分析所得倒塌点的中位值；
Figure GDA0002216427260000053
βR为易损性函数的对数标准差，βRR为偶然的不确定性，βRU为知识的不确定性或βM建模的不确定性，βD|IM和βC分别为地震动的不确定性和结构能力的不确定性。所建立各烈度、各层数RC框架结构地震易损性曲线，如图9所示。
本实施例的一种基于扩展PBEE2理论框架的地震易损性分析方法，第一部分中典型结构空间的建立，是根据我国建筑物特性统计资料及设计规范，基于逻辑树原理的建筑物分类标准并采用PKPM设计工具，设计出的用于数值分析的典型结构集。第二部分中在进行概率地震需求分析和地震易损性分析中，考虑了结构的非倒塌数值模型和倒塌数值模型，分别采用了基于位移和基于加速度的地震易损性模型，实现了各破坏极限状态下地震易损性函数的准确建立。第三部分中在形成一类结构的概率地震易损性函数时，根据第二部分的相关性分析结果，考虑了建筑物与建筑物间变异性所造成的不确定性(即层数、跨度等参数的变异性导致的结构抗震性能变异性)，使用了对数正态分布对易损性曲线进行拟合。
本实施例的一种基于扩展PBEE2理论框架的地震易损性分析方法，该方法给出了基于有限单体结构(典型结构)的一类结构解析地震易损性函数的建立方法，通过建立主导参数与地震易损性模型参数间的相关性矩阵，分析各结构设计参数对于地震易损性模型的影响规律，进而实现了一类结构解析地震易损性函数的建立，操作方便可靠，为城市建筑群的区域震害快速预测提供了理论基础，更加适用于城市建筑群的区域震害预测。
以上示意性的对本发明及其实施方式进行了描述，该描述没有限制性，附图中所示的也只是本发明的实施方式之一，实际的结构并不局限于此。所以，如果本领域的普通技术人员受其启示，在不脱离本发明创造宗旨的情况下，不经创造性的设计出与该技术方案相似的结构方式及实施例，均应属于本发明的保护范围。

Claims (3)

1.一种基于扩展PBEE2理论框架的地震易损性分析方法，其特征在于：其步骤为：
(1)建立典型结构空间：
(1.1)选取抗震设防烈度、高度和宽度比作为影响结构抗震性能的主导参数，其中层数、层高及跨度的选取，依照该类结构的大量数据统计分析得到；
(1.2)基于逻辑树原理，建立RC框架结构的典型结构空间；
(1.3)采用PKPM设计各典型结构，其恒载、活载取值依照该类结构的大量数据统计分析得到；
(2)建立主导参数与其概率地震需求模型参数间的相关性矩阵：
(2.1)基于OpenSEEs有限元软件平台,采用纤维模型和塑性铰模型，分别建立结构的非倒塌状态及倒塌状态数值模型；
(2.2)选取在未来地震下可能遭受的具有各种地震动特性的地震动记录，作为该类结构地震需求模型建立的输入地震动记录集，以考虑地震动的不确定性；
(2.3)采用基于位移和基于加速度的地震易损性模型，建立主导参数与其地震易损性模型参数间的相关性矩阵；
(3)基于对数正态分布的一类结构的概率地震易损性函数：
(3.1)根据典型结构设计参数与地震易损性函数参数间的数学模型关系，采用线性插值方法；计算其他跨度、层数非典型单体结构的地震易损性函数参数，从而建立一类结构中每个层数、跨度的结构地震易损性模型；
(3.2)以一类结构中不同典型结构的地震响应为总样本，按照单体结构的地震易损性建立方法，直接建立此类结构的地震易损性模型；
对于(3.2)中总不确定，按照下式考虑：
Figure FDA0002216427250000011
式中，βD为考虑地震动变异性、跨度变化的对数标准差，βM为建模不确定性。
2.根据权利要求1所述的一种基于扩展PBEE2理论框架的地震易损性分析方法，其特征在于：步骤(2.3)中基于典型结构建立主导参数与其地震易损性模型参数间的相关性矩阵。
3.根据权利要求2所述的一种基于扩展PBEE2理论框架的地震易损性分析方法，其特征在于：步骤(3)中给出了两种方法实现一类结构地震易损性模型的建立。
CN201610877011.7A 2016-10-08 2016-10-08 一种基于扩展pbee2理论框架的地震易损性分析方法 Active CN106649954B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610877011.7A CN106649954B (zh) 2016-10-08 2016-10-08 一种基于扩展pbee2理论框架的地震易损性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610877011.7A CN106649954B (zh) 2016-10-08 2016-10-08 一种基于扩展pbee2理论框架的地震易损性分析方法

Publications (2)

Publication Number Publication Date
CN106649954A CN106649954A (zh) 2017-05-10
CN106649954B true CN106649954B (zh) 2020-01-31

Family

ID=58853491

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610877011.7A Active CN106649954B (zh) 2016-10-08 2016-10-08 一种基于扩展pbee2理论框架的地震易损性分析方法

Country Status (1)

Country Link
CN (1) CN106649954B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107292545B (zh) * 2017-08-23 2020-05-26 中铁二院贵阳勘察设计研究院有限责任公司 一种基于地震风险评估的桥梁抗震分析方法
CN108256141B (zh) * 2017-12-11 2021-04-06 哈尔滨工业大学 一种基于Copula理论的主余震联合易损性分析方法
CN108595845B (zh) * 2018-04-26 2022-05-24 中冶华天工程技术有限公司 一种基于性能设计的桥梁地震风险概率分析方法
CN109241673A (zh) * 2018-10-11 2019-01-18 中国人民解放军63921部队 钢筋混凝土框架结构建筑物抗毁伤能力确定方法
CN109408965B (zh) * 2018-10-26 2023-04-28 中国地震局地球物理研究所 基于地震动参数的房屋震害矩阵曲线化分析系统及方法
CN111143931B (zh) * 2019-12-25 2022-11-22 上海市建筑科学研究院有限公司 一种基于增量静力方法的非迭代砌体结构易损性分析方法
CN111401747B (zh) * 2020-03-17 2023-05-19 上海市建筑科学研究院有限公司 序列型地震作用下砌体结构的易损性分析方法
CN112270125A (zh) * 2020-10-23 2021-01-26 同济大学 一种基于有向图逻辑模型和蒙特卡洛模拟的系统地震易损性分析方法
CN112329377A (zh) * 2020-11-02 2021-02-05 同济大学 用于量化多输入多输出变电站系统地震易损性的评估方法
CN112949149B (zh) * 2021-05-17 2021-07-23 西南石油大学 配气站设备工艺区管道地震易损性分析方法及分析系统
CN114117617B (zh) * 2021-12-08 2022-05-31 哈尔滨工业大学 一种三跨内廊式rc框架建筑地震响应的快速预测方法
CN114417465B (zh) * 2021-12-31 2022-10-11 清华大学 基于深度学习的隔震结构生成方法及装置
CN115270635B (zh) * 2022-08-04 2023-05-12 中国矿业大学 贝叶斯-神经网络高层建筑地震需求与易损性预测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102855219A (zh) * 2012-09-13 2013-01-02 同济大学 基于ann-mc-ud方法的公路桥梁地震易损性分析方法
CN105654414A (zh) * 2015-12-25 2016-06-08 浙江大学城市学院 基于开源体系架构和建筑物空间数据库的城市多灾害风险损失评测系统及方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102855219A (zh) * 2012-09-13 2013-01-02 同济大学 基于ann-mc-ud方法的公路桥梁地震易损性分析方法
CN105654414A (zh) * 2015-12-25 2016-06-08 浙江大学城市学院 基于开源体系架构和建筑物空间数据库的城市多灾害风险损失评测系统及方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Seismic fragility assessment of concrete bridge pier reinforced with Shape Memory Alloy considering residual displacement;Muntasir Billah 等;《Proceedings of SPIE - The International Society for Optical Engineering》;20120331;第83411F-1—83411F-13页 *
基于IDA方法的RC框架结构地震易损性分析研究;李文博;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20130215(第2期);第C038-463页 *
基于多元增量动力分析( MIDA) 方法的RC核心筒结构地震易损性分析;郑山锁 等;《振动与冲击》;20150115;第34卷(第1期);第117-123,189页 *
节点变形对在役混凝土框架地震易损性影响分析;郑山锁 等;《工业建筑》;20141220;第44卷(第12期);第63-67页 *
钢筋混凝土结构易损性分析与地震风险研究;朱健;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20101015(第10期);第C038-12页 *

Also Published As

Publication number Publication date
CN106649954A (zh) 2017-05-10

Similar Documents

Publication Publication Date Title
CN106649954B (zh) 一种基于扩展pbee2理论框架的地震易损性分析方法
CN111458748B (zh) 基于三层数据集神经网络的性能地震动危险性分析方法
CN105606063B (zh) 一种基于正交应变比的土层边坡稳定性测定方法
CN108256141B (zh) 一种基于Copula理论的主余震联合易损性分析方法
Ortega et al. A vulnerability index formulation for the seismic vulnerability assessment of vernacular architecture
WO2019148950A1 (zh) 城市建筑群地震反应非线性历程分析方法及装置
Sun et al. Study on vulnerability matrices of masonry buildings of mainland China
Bhatia et al. Analyzing architectural space: identifying salient regions by computing 3D isovists
CN110503267B (zh) 基于时空尺度自适应模型的城市侵财案件预测系统及预测方法
CN103134433B (zh) 一种利用位移监测鉴别边坡失稳致滑因子的方法
CN114792020B (zh) 一种基于机器学习的建筑抗震韧性快速评估方法及系统
Zhou et al. An area‐based intensity measure for incremental dynamic analysis under two‐dimensional ground motion input
Nale et al. Fragility functions for local failure mechanisms in unreinforced masonry buildings: a typological study in Ferrara, Italy
Olivito et al. Inventory and monitoring of historical cultural heritage buildings on a territorial scale: a preliminary study of structural health monitoring based on the CARTIS approach
Tsalouchidis et al. Amplitude scaling of ground motions as a potential source of bias: large‐scale investigations on structural drifts
CN116738553B (zh) 抗震结构的不确定参数敏感性分析方法
CN113657610A (zh) 一种基于随机森林的冰雹气候特征预测方法
Olivito et al. SHM systems applied to the built heritage inventory at the territorial scale. A preliminary study based on CARTIS approach
CN111582634A (zh) 一种地下大空间施工多因素安全分级方法及系统
CN107220705A (zh) 常减压装置常压塔顶干点预测方法
CN112883339B (zh) 一种地震有感范围确定方法及系统
CN113298426B (zh) 知识图谱驱动的大坝安全评判权重动态拟定方法及系统
CN104462739B (zh) 一种适用于大范围多尺度卫星遥感数据反演的生态环境参数地面采样方法
CN114021644A (zh) 一种基于K-means和去丛聚法的区域代表性地下水位计算方法
CN109241673A (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