CN111695268A - 智慧城市震灾推演快速动力时程分析方法 - Google Patents
智慧城市震灾推演快速动力时程分析方法 Download PDFInfo
- Publication number
- CN111695268A CN111695268A CN202010554048.2A CN202010554048A CN111695268A CN 111695268 A CN111695268 A CN 111695268A CN 202010554048 A CN202010554048 A CN 202010554048A CN 111695268 A CN111695268 A CN 111695268A
- Authority
- CN
- China
- Prior art keywords
- mode
- displacement
- vibration
- iteration
- order
- 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
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 15
- 238000006073 displacement reaction Methods 0.000 claims abstract description 79
- 238000010008 shearing Methods 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 23
- 239000000758 substrate Substances 0.000 claims abstract description 14
- 230000003595 spectral effect Effects 0.000 claims description 15
- 230000001133 acceleration Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000013016 damping Methods 0.000 claims description 6
- 230000002194 synthesizing effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 11
- 239000010410 layer Substances 0.000 description 7
- 238000005452 bending Methods 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000004088 simulation Methods 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
Abstract
本发明涉及一种地震分析方法,本发明的智慧城市震灾推演快速动力时程分析方法,包括如下步骤:建立智慧城市需要监控的建筑信息模型,将数据导入平台;计算各阶振动周期;获得不同振型下基底剪力与结构顶点位移的关系;选取振型数,获得不同模态下顶点位移与基底剪力曲线;计算第一振型的折减系数;计算单自由度动力体系弹塑性微分方程;计算多振型下结构的顶点位移;合成的多模态位移响应,确定不同楼层多模态位移响应比例,按照多模态位移响应比例,分配各楼层的位移响应;将各楼层的位移响应,实时动态加载建立的平台,实现智慧城市震灾的快速推演。本发明提高了计算精度和速度,满足智慧城市地震后震害快速推演的需求。
Description
技术领域
本发明涉及一种地震分析方法,尤其涉及一种地震动推演分析方法。
背景技术
地震是危及人类生命财产安全的突发式自然灾害,是对人类威胁最严重的自然灾害。全球每年会发生500多万次地震,地震几乎时时刻刻都有发生,不可避免。近一个世纪以来,世界范围内有20多座城市遭到地震灾害的毁灭性破坏,著名的地震包括日本关东(1923)、智利(1960)、阿拉斯加(1957、1964、1965)、唐山(1976)、墨西哥(1985)、中国台湾集集(1999)等大地震。近年来,世界各地地壳运动活跃,地震发生频繁,全球地震活动以2004年12月26日苏门答腊9.0级、2008年5月12日汶川9.0级、2011年3月11日日本宫城县9.0级、2018年8月19日斐济群岛等系列地震为标志,已经进入了8级地震活动的高发期,到目前为止共发生8级以上地震15次,9级以上地震2次,呈现了类似20 世纪前半叶全球大地震多发的状态,因此,研发防灾减灾关键技术对智慧城市的建设具有重大意义。
智慧城市在地震后能快速进行震害的推演,做出震灾研判和救灾部署,是智慧城市防灾减灾的重要内容,是智慧城市御震韧性的重要体现。基于时程动力分析的震害推演方法,主要有东京大学的区域建筑震害模拟方法Integrated Earthquake Simulation(2008)、清华大学的Lu方法(2014)和Xiong方法(2016),前二种方法基于多自由度剪切层模型,进行动力弹塑性时程分析,第三种方法考虑弯曲-剪切耦合变形,基于剪切/弯曲模型进行动力弹塑性时程分析。
以上三种方法存在以下问题:(1)结构自振周期是按经验公式预估出来,对计算结果影响大,计算准确性有待完善;(2)剪切或剪切/弯曲模型的动力弹塑性时程分析,计算时间长、计算效率差,难以满足地震后快速分析震害的需求;(3)剪切或剪切/弯曲模型难以充分考虑高振型的影响,计算精度差。
发明内容
为了解决以上问题,本发明提出了一种智慧城市震灾推演快速动力时程分析方法,该方法满足智慧城市地震后震害快速推演的需求。
为实现上述目的,本发明采用如下技术方案:包括如下步骤:
步骤1,建立智慧城市需要监控的建筑信息模型BIM,将数据导入智慧城市 GIS平台;
步骤2,计算步骤1中需要监控的建筑各阶振动周期;
步骤3,采用PUSHOVER推覆方法,获得不同振型下基底剪力与结构顶点位移的关系,指定不同动力状态,计算不同动力状态下基底剪力;
步骤4,选取振型数,获得不同模态下顶点位移与基底剪力曲线;
步骤5,根据步骤3确定的不同动力状态的基底剪力,乘以不同模态剪力调整系数ζj=Ηj/Η1,其中Η1为第一阶振型参与系数,Ηj为第一阶至最大阶振型数的振型参与系数,j=1,2,3,……,r,r为最大振型数;
步骤6,在二阶以上最大振型数以下模态PUSHOVER推覆曲线中,根据步骤 5,计算动力状态的基底剪力乘以不同模态剪力调整系数ζj后的基底剪力值i为第一次至n次的迭代数;Vi是第i次迭代时,根据上次迭代在第一阶振型下,顶点位移对应的基底剪力;Vip是第i次迭代时,p 阶振型下动力状态的基底剪力,基于二阶以上最大振型数以下模态PUSHOVER推覆曲线,根据Vip确定不同模态下的顶点位移uip,uip是第i次迭代时,p阶振型下动力状态的位移;
步骤8,第一次迭代时,根据步骤3计算的第一振型顶点位移u11,计算第一次迭代位移折减系数θ1=u11/U1,第一阶振型第二次迭代的位移值u21=θ1*u11;同理,计算第一次迭代后,第i次迭代位移折减系数θi=ui1/Ui,ui+11=θi*ui1,其中i=1,2,3,……,n,n为最后一次迭代次数;
步骤9,根据计算出的u21替代步骤8的u11,重复步骤5-7,依次类推,直到(Un-Un-1)/Un<δ,迭代终止,其中δ=1%;
步骤10,当迭代终止时,计算第一振型的折减系数,θi=ui1/Ui,计算第一振型位移占结构总动力状态位移的百分比;
步骤11,将步骤3转化为第一模态下的单自由度体系的方程:
D1(t)第一模态的动力响应,κ1第一模态阻尼比,ω1为第一模态频率,f1(t) 为第一模态下的谱加速度值和谱位移关系函数;
步骤14,按以下公式计算第一模态下结构顶点位移
步骤15,根据步骤10,计算多振型下结构的顶点位移
Utop(t)=u1(t)/θn;
步骤16,基于步骤3和4多模态PUSHOVER推覆结果,按SRSS规则,合成的多模态位移响应,确定不同楼层多模态位移响应比例,基于Utop(t)按照多模态位移响应比例,分配各楼层的位移响应;
步骤17,将各楼层的位移响应,实时动态加载基于步骤1BIM的建立的GIS 平台,实现智慧城市震灾的快速推演。
优选地,步骤1中,采用IFC数据导入智慧城市GIS平台。
优选地,步骤3中,第一振型下最大层间位移角与结构顶点位移的关系获得通过采用第一阶模态PUSHOVER推覆结构。
优选地,步骤3中,不同动力状态下的基底剪力,通过以下方法获得:
在地震力作用下多自由度体系结构的运动微分方程:
按照多模态PUSHOVER推覆结构,n阶模态推覆力fn为:
优选地,步骤4中,选取振型数按照振型数应使各振型参与质量之和不小于结构总质量的90%的原则;获得不同模态下顶点位移与基底剪力曲线采用二阶以上最大振型数以下的模态PUSHOVER推覆结构。
优选地,步骤11中,第一模态下的谱加速度值和谱位移关系函数f1(t),按下式计算:
优选地,步骤13中,计算单自由度动力体系弹塑性微分方程采用Wilson- θ法。
本发明智慧城市震灾推演快速动力时程分析方法,一是提高了计算精度,二是提高了计算速度,满足智慧城市地震后震害快速推演的需求。
附图说明
图1为本发明第一振型PUSHOVER推覆曲线;
图2为本发明第二振型PUSHOVER推覆曲线;
图3为本发明第三振型PUSHOVER推覆曲线。
具体实施方式
本发明的智慧城市震灾推演快速动力时程分析方法,包括如下步骤:
1、基于BIM技术,建立智慧城市需要监控的建筑信息模型,采用IFC数据导入智慧城市GIS平台,例如超图SuperMap GIS平台。
2、计算建筑各阶振动周期。
3、采用PUSHOVER推覆方法,获得不同振型下基底剪力与结构顶点位移的关系,指定不同动力状态,计算不同动力状态下基底剪力。并获得不同动力状态下的基底剪力。
在地震力作用下多自由度体系结构的运动微分方程:
假定按照多模态PUSHOVER推覆结构,n阶模态推覆力fn
4、按照计算振型数应使各振型参与质量之和不小于结构总质量的90%的原则,选取振型数r。采用二阶以上最大振型数以下的模态PUSHOVER推覆结构,获得不同模态下顶点位移与基底剪力曲线。
假定应考虑振型为3个,第二、三PUSHOVER推覆曲线见图2、图3。
5、根据步骤3确定的不同动力状态的基底剪力,乘以不同模态剪力调整系数ζj=Ηj/Η1,其中Η1为第一阶振型参与系数,Ηj为第一阶至最大阶振型数的振型参与系数,j=1,2,3,……,r,r为最大振型数。
6、在二阶以上最大振型数以下模态PUSHOVER推覆曲线中,根据步骤5,计算动力状态的基底剪力乘以不同模态剪力调整系数ζj后的基底剪力值i为第一次至n次的迭代数;Vi是第i次迭代时,根据上次迭代在第一阶振型下,顶点位移对应的基底剪力;Vip是第i次迭代时,p 阶振型下动力状态的基底剪力,基于二阶以上最大振型数以下模态PUSHOVER推覆曲线,根据Vip确定不同模态下的顶点位移uip,uip是第i次迭代时,p阶振型下动力状态的位移。
8、第一次迭代时,根据步骤3计算的第一振型顶点位移u11,第一阶振型第二次迭代的位移值u21=θ1*u11;同理,计算第一次迭代后,第i次迭代位移折减系数μi=ui1/Ui,ui+11=θi*ui1,其中i=1,2,3,……,n,n为最后一次迭代次数。
9、根据新计算的u21替代步骤8的u11,重复步骤5-7,依次类推,直到 (Un-Un-1)/Un<δ,迭代终止,其中δ=1%。
10、当迭代终止时,计算第一振型的折减系数,θi=ui1/Ui,计算第一振型位移占结构总动力状态位移的百分比。
11、将步骤3转化为第一模态下的单自由度体系的方程:
D1(t)第一模态的动力响应,κ1第一模态阻尼比,ω1为第一模态频率,f1(t)为第一模态下的谱加速度值和谱位移关系函数,第一模态下的谱加速度值和谱位移关系函数f1(t),按下式计算:
13、将代入方程(3),采用Wilson-θ法,快速计算单自由度动力体系弹塑性微分方程。该方法避免了已有技术需要计算多自由度动力体系弹塑性微分方程,计算时间数千倍节省,采用了准确的结构周期,单自由度动力体系弹塑性微分方程的计算精度得到明显提高。
14、按以下公式计算第一模态下结构顶点位移
15、根据步骤10,计算考虑多振型下结构的顶点位移
Utop(t)=u1(t)/θn (10)
这一步考虑了多阶振型的因素,因此,计算精度得到提高。
16、基于步骤3和4多模态PUSHOVER推覆结果,按SRSS规则,合成的多模态位移响应,确定不同楼层多模态位移响应比例,基于Utop(t)按照多模态位移响应比例,分配各楼层的位移响应。
17、将各楼层的位移响应,实时动态加载步骤1建立的BIM+GIS平台,实现智慧城市震灾的快速推演。
Claims (6)
1.一种智慧城市震灾推演快速动力时程分析方法,其特征在于,包括如下步骤:
步骤1,建立智慧城市需要监控的建筑信息模型BIM,将数据导入智慧城市GIS平台;
步骤2,计算步骤1中需要监控的建筑各阶振动周期;
步骤3,采用PUSHOVER推覆方法,获得不同振型下基底剪力与结构顶点位移的关系,指定不同动力状态,计算不同动力状态下基底剪力;
步骤4,选取振型数,获得不同模态下顶点位移与基底剪力曲线;
步骤5,根据步骤3确定的不同动力状态的基底剪力,乘以不同模态剪力调整系数βj=Γj/Γ1,其中Γ1为第一阶振型参与系数,Γj为第一阶至最大阶振型数的振型参与系数,j=1,2,3,……,r,r为最大振型数;
步骤6,在二阶以上最大振型数以下模态PUSHOVER推覆曲线中,根据步骤5,计算动力状态的基底剪力乘以不同模态剪力调整系数βj后的基底剪力值i为第一次至n次的迭代数;Vi是第i次迭代时,根据上次迭代在第一阶振型下,顶点位移对应的基底剪力;Vip是第i次迭代时,p阶振型下动力状态的基底剪力,基于二阶以上最大振型数以下模态PUSHOVER推覆曲线,根据Vip确定不同模态下的顶点位移uip,uip是第i次迭代时,p阶振型下动力状态的位移;
步骤8,第一次迭代时,根据步骤3计算的第一振型顶点位移u11,计算第一次迭代位移折减系数θ1=u11/U1,第一阶振型第二次迭代的位移值u21=θ1*u11;同理,计算第一次迭代后,第i次迭代位移折减系数θi=ui1/Ui,ui+1 1=θi*ui 1,其中i=1,2,3,……,n,n为最后一次迭代次数;
步骤9,根据计算的u21替代步骤8的u11,重复步骤5-7,依次类推,直到(Un-Un-1)/Un<δ,迭代终止,其中δ=1%;
步骤10,当迭代终止时,计算第一振型的折减系数,θi=ui1/Ui,计算第一振型位移占结构总动力状态位移的百分比;
步骤11,将步骤3转化为第一模态下的单自由度体系的方程:
D1(t)为第一模态的动力响应,ζ1为第一模态阻尼比,ω1为第一模态频率,f1(t)为第一模态下的谱加速度值和谱位移关系函数;
步骤14,按以下公式计算第一模态下结构顶点位移
步骤15,根据步骤10,计算多振型下结构的顶点位移
Utop(t)=u1(t)/θn;
步骤16,基于步骤3和4多模态PUSHOVER推覆结果,按SRSS规则,合成的多模态位移响应,确定不同楼层多模态位移响应比例,基于Utop(t)按照多模态位移响应比例,分配各楼层的位移响应;
步骤17,将各楼层的位移响应,实时动态加载基于步骤1的BIM建立的GIS平台,实现智慧城市震灾的快速推演。
2.根据权利要求1所述的智慧城市震灾推演快速动力时程分析方法,其特征在于,步骤1中,采用IFC数据导入智慧城市GIS平台。
4.根据权利要求1所述的智慧城市震灾推演快速动力时程分析方法,其特征在于,步骤4中,选取振型数按照振型数应使各振型参与质量之和不小于结构总质量的90%的原则;获得不同模态下顶点位移与基底剪力曲线采用二阶以上最大振型数以下的模态PUSHOVER推覆结构。
6.根据权利要求1所述的智慧城市震灾推演快速动力时程分析方法,其特征在于,步骤13中,计算单自由度动力体系弹塑性微分方程采用Wilson-θ法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010554048.2A CN111695268B (zh) | 2020-06-17 | 2020-06-17 | 智慧城市震灾推演快速动力时程分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010554048.2A CN111695268B (zh) | 2020-06-17 | 2020-06-17 | 智慧城市震灾推演快速动力时程分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111695268A true CN111695268A (zh) | 2020-09-22 |
CN111695268B CN111695268B (zh) | 2022-05-31 |
Family
ID=72481666
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010554048.2A Expired - Fee Related CN111695268B (zh) | 2020-06-17 | 2020-06-17 | 智慧城市震灾推演快速动力时程分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111695268B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113484911A (zh) * | 2021-04-30 | 2021-10-08 | 石家庄学院 | 一种基于地震动演化功率谱的地震位移反应计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102787677A (zh) * | 2011-05-20 | 2012-11-21 | 青岛理工大学 | 工程结构抗震烈度预测方法 |
US20120310607A1 (en) * | 2011-05-20 | 2012-12-06 | Qingdao Technological University | Engineering structure multi-objective performance-based seismic design |
CN106049951A (zh) * | 2016-05-25 | 2016-10-26 | 青岛理工大学 | 多级地震作用下工程结构抗震性能设计评估方法 |
CN108460200A (zh) * | 2018-02-12 | 2018-08-28 | 青岛理工大学 | 基于性态双谱的多模态性能化抗震设计方法 |
CN108647366A (zh) * | 2018-02-05 | 2018-10-12 | 清华大学 | 城市建筑群地震反应非线性历程分析方法及装置 |
-
2020
- 2020-06-17 CN CN202010554048.2A patent/CN111695268B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102787677A (zh) * | 2011-05-20 | 2012-11-21 | 青岛理工大学 | 工程结构抗震烈度预测方法 |
US20120310607A1 (en) * | 2011-05-20 | 2012-12-06 | Qingdao Technological University | Engineering structure multi-objective performance-based seismic design |
CN106049951A (zh) * | 2016-05-25 | 2016-10-26 | 青岛理工大学 | 多级地震作用下工程结构抗震性能设计评估方法 |
CN108647366A (zh) * | 2018-02-05 | 2018-10-12 | 清华大学 | 城市建筑群地震反应非线性历程分析方法及装置 |
CN108460200A (zh) * | 2018-02-12 | 2018-08-28 | 青岛理工大学 | 基于性态双谱的多模态性能化抗震设计方法 |
Non-Patent Citations (3)
Title |
---|
孙江伟: "多级地震作用下剪力墙结构的多模态Pushover分析", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
张曦: "高层建筑结构高宽比与剪重比问题分析", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
朱春明: "中、美规范在建筑结构抗震设计中的对比", 《钢结构》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113484911A (zh) * | 2021-04-30 | 2021-10-08 | 石家庄学院 | 一种基于地震动演化功率谱的地震位移反应计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111695268B (zh) | 2022-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shibata et al. | Substitute-structure method for seismic design in R/C | |
Chopra et al. | A modal pushover analysis procedure to estimate seismic demands for unsymmetric‐plan buildings | |
Kalkan et al. | Assessment of current nonlinear static procedures for seismic evaluation of buildings | |
CN105868477A (zh) | 基于层间位移角约束的结构抗震敏感性优化方法 | |
CN103774772B (zh) | 一种控制基础隔震结构扭转的方法 | |
CN111695268B (zh) | 智慧城市震灾推演快速动力时程分析方法 | |
Tarbali et al. | Story shear and torsional moment-based pushover procedure for asymmetric-plan buildings using an adaptive capacity spectrum method | |
CN108460200B (zh) | 基于性态双谱的多模态性能化抗震设计方法 | |
Kolay et al. | Multi-hazard real-time hybrid simulation of a tall building with damped outriggers | |
CN106372324B (zh) | 基于地震剪力系数约束的结构抗震敏感性优化方法 | |
CN111767596A (zh) | 基于性能水准的地震动评估方法 | |
Prasanth et al. | Estimation of hysteretic energy demand using concepts of modal pushover analysis | |
Nooranad | Analytical investigation on the performance of tube-in-tube structures subjected to lateral loads | |
Xiong et al. | The Multangular-Pyramid Concave Friction System (MPCFS) for seismic isolation: A preliminary numerical study | |
CN108427833B (zh) | 基于三维弹塑性地震反应谱的性能化抗震设计方法 | |
Manoukas et al. | Evaluation of a multimode pushover procedure for asymmetric in plan buildings under biaxial seismic excitation | |
Khy et al. | Application of conditional mean spectrum in nonlinear response history analysis of tall buildings on soft soil | |
Lepage et al. | Best-fit models for nonlinear seismic response of reinforced concrete frames | |
Pian et al. | A hybrid force/displacement seismic design method for reinforced concrete infilled-MRFs and wall-frame dual systems | |
Karimi et al. | A three-dimensional drift pushover method for unsymmetrical plan buildings | |
Kim et al. | An efficient structural analysis of super tall mega frame buildings using a multi-level condensation method | |
Manoukas | Evaluation of a multimode pushover procedure for asymmetric in plan and non-regular in elevation R/C buildings | |
CN110727026B (zh) | 一种基于地震-海啸连续作用的非弹性位移比谱模型 | |
Qu et al. | An improved multidimensional modal pushover analysis procedure for seismic evaluation of latticed arch‐type structures under lateral and vertical earthquakes | |
Numada et al. | Material point method for run-out analysis of earthquake-induced long-traveling soil flows |
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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Liu Wenfeng Inventor after: Han Shuangyang Inventor after: Zhou Zheng Inventor before: Liu Wenfeng Inventor before: Chen Guanjun Inventor before: Zhou Zheng |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220531 |
|
CF01 | Termination of patent right due to non-payment of annual fee |