CN112948983A - 一种协同正面碰撞多工况的汽车前端结构能量管理方法 - Google Patents
一种协同正面碰撞多工况的汽车前端结构能量管理方法 Download PDFInfo
- Publication number
- CN112948983A CN112948983A CN202110422464.1A CN202110422464A CN112948983A CN 112948983 A CN112948983 A CN 112948983A CN 202110422464 A CN202110422464 A CN 202110422464A CN 112948983 A CN112948983 A CN 112948983A
- Authority
- CN
- China
- Prior art keywords
- vehicle
- barrier
- collision
- energy
- end structure
- 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
- 238000007726 management method Methods 0.000 title claims abstract description 12
- 230000004888 barrier function Effects 0.000 claims abstract description 133
- 238000010521 absorption reaction Methods 0.000 claims abstract description 66
- 238000000034 method Methods 0.000 claims abstract description 49
- 230000033001 locomotion Effects 0.000 claims abstract description 37
- 230000004044 response Effects 0.000 claims abstract description 29
- 230000008569 process Effects 0.000 claims abstract description 24
- 230000001133 acceleration Effects 0.000 claims abstract description 16
- 230000015556 catabolic process Effects 0.000 claims abstract description 5
- 238000006073 displacement reaction Methods 0.000 claims description 34
- 230000000452 restraining effect Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 238000000354 decomposition reaction Methods 0.000 description 38
- 238000011156 evaluation Methods 0.000 description 25
- 238000004364 calculation method Methods 0.000 description 20
- 238000013461 design Methods 0.000 description 13
- 238000011160 research Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 238000012360 testing method Methods 0.000 description 6
- 230000006378 damage Effects 0.000 description 5
- 230000009545 invasion Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 238000012795 verification Methods 0.000 description 5
- 208000027418 Wounds and injury Diseases 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 208000014674 injury Diseases 0.000 description 4
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 3
- 229910052782 aluminium Inorganic materials 0.000 description 3
- 230000001413 cellular effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 238000012356 Product development Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 241001391944 Commicarpus scandens Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004141 dimensional analysis Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000009467 reduction Effects 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/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Body Structure For Vehicles (AREA)
Abstract
本发明公开了一种协同正面碰撞多工况的汽车前端结构能量管理方法,包括:步骤一、将车辆前端吸能总空间划分为多个吸能子空间,并将FRB工况碰撞原波形分解为与吸能子空间一一对应的子波形,建立车辆前端结构三维解析模型;每个子波形的幅值为FRB工况碰撞原波形的幅值与其对应的吸能子空间的吸收能量占比的乘积;步骤二、基于MPDB工况,对前端结构三维解析模型进行求解,得到车辆和壁障的运动响应,确定碰撞过程中壁障和车辆的变形、壁障和车辆的加速度及壁障和车辆的速度;并根据运动响应得到乘员负载指数、壁障表面击穿情况、壁障的均匀性指标和兼容性罚分;以及基于SOB工况,对前端结构三维解析模型进行求解,得到碰撞结束时车辆乘员舱的形变量。
Description
技术领域
本发明属于汽车碰撞解析模型的建模及求解技术领域,特别涉及一种协同正面碰撞多工况的汽车前端结构能量管理方法。
背景技术
汽车正面碰撞安全性指汽车在发生正面碰撞事故时,车体结构按照既定的方式变形,并与乘员约束系统(包括安全带、安全气囊、吸能式转向柱、座椅和内饰件等)协同作用,有效地吸收碰撞能量,最大限度地减轻车内、外人员碰撞损伤的能力。汽车正面碰撞安全性根据保护对象的不同可以进一步划分为,以保护自身车内乘员安全为目标的汽车耐撞性设计,和以保护对方车内乘员安全为目标的汽车兼容性设计。目前,关于车辆耐撞性的研究日趋成熟,利用汽车正面碰撞概念设计方法已经可以实现车辆结构和约束系统的部分正向设计,法律规范也较为成熟,但主要适用于单工况设计,且倾向于对自车乘员的保护,忽略了汽车碰撞兼容性问题。而关于车辆兼容性的研究,直到2020年E-NCAP以及2021年C-NCAP启用MPDB碰撞工况之前,仍未被法规明确提出并且引起广泛重视,在此前关于车辆兼容性的研究仅涉及到个别自主法规。新型壁障对车碰撞试验,从原来的车与固定壁障的碰撞升级为车与带壁障的移动小车之间的相互碰撞,更加真实地模拟车辆前部碰撞安全事故,即车对车碰撞,而不是车与固定物碰撞,车辆与壁障碰撞速度分别为50km/h,重叠量为50%,台车质量为1400kg采用新型MPDB蜂窝铝,与此前的ODB工况采用的蜂窝铝相比不易被击穿,更接近于真实工况,离地高度150mm,在主驾驶位置放置THRO 50th男性假人,在副驾驶放置Hybrid III50th男性假人,在后排放置Q10、Q6儿童假人考察车辆的乘员保护性能。与现有评价指标相对比,除了对乘员进行伤害值评价,还增加了兼容性评价,与伤害值评价不同的是,兼容性评价为罚分制,包含了三个评价参数:壁障的OLC值、壁障穿透和壁障均匀性系数h。新法规规定车辆需同时满足耐撞性和兼容性,然而如今现有研究中,关于车辆兼容性的研究主要是利用仿真技术,缺少相关的理论研究以及数据作指导,通常需要在结构设计完毕之后通过试验或CAE(Computer Aided Engineering)方法获得车辆与壁障前端变形数据以及吸收能量,再根据MPDB评价指标进行评价,如果评价结果达不到法规要求,则需要对车辆前端结构进行反复修改直至满足评价标准。
此外,在所有正面碰撞死亡事故中,小偏置正面碰撞约占正面碰撞总量的四分之一。且根据已有研究表明,在正面全宽、正面40%偏置以及正面25%重叠率碰撞工况中车辆在IIHS正面25%重叠率碰撞工况的满分率是最低的。并且该工况下乘员舱侵入量较大、乘员伤害较严重。但是直到2012年,各国安全性法规及评价规程中,仅IIHS在2012年对小偏置这种正面碰撞类型做出了评价,随着国内消费者对汽车安全性重视程度和国内主机厂被动安全设计能力的提升,2017年中国保险汽车安全指数发布,引入25%偏置碰撞测试,受到媒体和社会的广泛关注,并逐渐受到各主机厂的重视,成为业内研究的新热点。但是,目前各企业及研究机构对该工况的研究主要集中在车辆结构修改上,通常是根据某款车在该工况试验中出现的吸能不足、乘员舱侵入的情况,对结构进行修改,再通过试验或CAE方法对改进后的车辆结构进行耐撞性的验证。这种结构修改方法,很难把握乘员舱侵入量减小的程度,只能通过乘员伤害预测反复修改结构来达到乘员保护的目的,增加了设计的难度和项目开发的风险。
发明内容
本发明的目的是针对MPDB及SOB偏置碰撞工况的动力学响应,主要是利用仿真技术获得的研究现状,提供了一种协同正面碰撞多工况的汽车前端结构能量管理方法,实现不同工况条件下车辆响应的快速求解,避免了大量重复的试算过程,缩短产品开发周期。
本发明提供的技术方案为:
一种协同正面碰撞多工况的汽车前端结构能量管理方法,包括如下步骤:
步骤一、根据总布置阶段的车辆硬点尺寸确定车辆前端吸能总空间;将所述车辆前端吸能总空间划分为多个吸能子空间,并将FRB工况碰撞原波形分解为与所述吸能子空间一一对应的子波形,建立车辆前端结构三维解析模型;
其中,每个所述子波形的幅值为FRB工况碰撞原波形的幅值与其对应的吸能子空间的吸收能量占比的乘积;
步骤二、基于MPDB工况,对所述前端结构三维解析模型进行求解,得到车辆和壁障的运动响应,确定碰撞过程中壁障和车辆的变形、壁障和车辆的加速度及壁障和车辆的速度;并根据所述运动响应得到乘员负载指数、壁障表面击穿情况、壁障的均匀性指标和兼容性罚分;以及
基于SOB工况,对所述前端结构三维解析模型进行求解,得到碰撞结束时车辆乘员舱的形变量。
优选的是,在所述步骤一中,所述车辆前端解析模型的建立方法为:
将车辆前端每个吸能子空间都简化为刚度为kij的弹簧,将车辆简化为质量块;
其中,
式中,qij为吸能子空间的吸收能量占比;av为碰撞波形,dv为车辆位移,M为整车质量,i为车辆前端吸能总空间沿车辆的横向分解的个数,j为车辆前端吸能总空间沿车辆的垂向分解的个数。
优选的是,在所述步骤二中,对所述前端结构三维解析模型进行求解之前,还包括:
基于MPDB工况,将壁障在横向和垂向划分为与车辆前端的吸能子空间相对应的多个区域,将所述多个区域分别简化为弹簧;并根据机械原理将壁障和车辆的碰撞模型简化为等效单自由度模型。
优选的是,在所述步骤二中,乘员负载指数通过如下公式进行计算:
式中,OLC表示乘员负载指数,t1为乘员运动阶段向前移动65mm的时间,t2为乘员约束阶段向前移动235mm时的时间,vB(t)为壁障的速度,v0为壁障初速度。
优选的是,在所述步骤二中,兼容性罚分的计算公式为:
当OLC<25g时,Mcompat=-2·h-MBO;
当25g≤OLC≤40g时,
Mcompat=-2·OLC/15+10/3-h·((4·OLC/10-8)-(2·OLC/15-10/3))-MBO;
当OLC>40g时,Mcompat=-2-6·h-MBO;
式中,h为均匀系数,MBO为触底罚分。
优选的是,所述均匀系数,根据壁障的均匀性指标确定;其中:
当SD<50时,h=0;
当50≤SD≤150时,h=(SD-50)/100;
当SD≥150时,h=100%;
式中,SD为壁障的均匀性指标。
优选的是,所述壁障均匀性指标的计算公式为:
优选的是,在所述步骤二中,车辆乘员舱的形变量通过如下公式进行计算:
D25%max=max(xv)-L1;
式中,D25%max为车辆乘员舱的形变量,max(xv)为碰撞过程结束时车辆的纵向最大形变量,L1为防撞梁前端到防火墙的距离。
优选的是,将通过所述三维解析模型计算得到的车辆的最大纵向位移作为车辆的纵向最大形变量。
优选的是,在所述步骤二中,基于SOB工况判断碰撞过程是否结束的方法为:
当vv≤0或者y≥25%B时,判定碰撞过程结束;
其中,vv为车辆速度,y为车辆的横向位移,B为车辆宽度;
y=rsinθ-25%×B;
式中,θ为车辆的旋转角度,r为车辆的旋转半径。
本发明的有益效果是:
(1)本发明基于三维分解方法划分吸能空间和FRB工况碰撞波形,在概念设计阶段,缺少三维数模的情况下,能够根据不同工况条件改变模型参数,适用于多工况条件下车辆响应的快速求解。
(2)本发明提供的三维解析模型,揭示了MPDB碰撞工况下壁障、车辆的动力学关系,还有二者在碰撞过程的相互作用机理;以及SOB碰撞工况的车辆响应和碰撞结束条件;
(3)本发明建立了车辆兼容性和乘员舱安全性的评价指标,以车辆基本信息、碰撞工况和碰撞波形的分解方案为输入,利用振动方程实现不同工况条件下车辆响应的求解,从而评价车辆是否满足耐撞性,兼容性和乘员舱的安全性要求。
(4)本发明提供的三维解析模型适用于汽车碰撞安全性的概念设计和碰撞能量管理阶段,可以通过求解车辆在MPDB以及SOB碰撞工况下的评价指标对车辆前端结构的能量管理进行优化,使车辆结构的设计更具有针对性,避免了大量重复的有限元试算过程,节约人力、财力,缩短产品开发周期。
附图说明
图1为本发明所述的协同正面碰撞多工况的汽车前端结构能量管理方法的原理图。
图2a-2b为本发明所述的正面碰撞能量的初步分解方案。
图3为本发明所述的车辆前端结构三维解析模型图。
图4为本发明所述的MPDB碰撞工况的三维模型(模型在车身坐标系XZ面投影)示意图。
图5为本发明所述的根据机械原理简化的壁障与车辆的等效单自由度模型。
图6为本发明所述的OLC的计算图像。
图7为本发明所述的SOB碰撞工况的三维模型(模型在车身坐标系XY面投影)示意图。
图8为本发明所述的车辆在SOB碰撞工况旋转过程中的受力分析图。
图9为本发明所述的在MPDB工况模型验证环节所选的三个车型在FRB工况下的碰撞波形图。
图10为本发明所述的在SOB工况模型验证环节所选的某车型在FRB工况下的碰撞波形图。
具体实施方式
下面结合附图对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
如图1所示,本发明提供了一种协同正面碰撞多工况的汽车前端结构能量管理方法,基于三维分解方法划分吸能空间和碰撞波形,以吸能空间和碰撞波形的分解方案为输入,利用振动方程实现不同工况条件下车辆响应的求解;在概念设计阶段,没有三维数模的情况下,根据总布置参数确定汽车正面碰撞的能量管理方案,使车辆在满足耐撞性的同时,还能兼顾车辆的兼容性和乘员舱的安全性要求。本发明是基于自然基金项目(项目批准号:51775228,名称:汽车车身多材料复杂断面薄壁梁结构抗撞性理论模型研究)实施的。
具体求解过程如下:
1、输入参数
1)根据经验以及车身前端结构特点对FRB工况碰撞波形进行分解
本次所选取的车辆前端结构的三维分解为横向50%+25%区(即在在横向划分四个区,每个区的宽度为25%车宽),垂向三层,纵向二段的分解方案,其中,纵向第一段即为车辆前端吸能总空间,纵向第二段布置发动机等刚性部件,不参与变形吸能。
由于横向两个50%区域的吸能结构完全相同,吸能量各占比均为50%。考虑SOB工况的吸能效果,25%区的吸能量最好不要少于总能量的10%,同时考虑车轮的占用空间和轴距的设置要求,该区能用于布置吸能结构的空间有限,吸能量很难达到总能量的30%。垂向第二层作为正面碰撞的主要吸能区,其能量占比不少于40%。同时考虑MPDB工况中前舱吸能结构布置的均匀性要求,第一层和第三层的吸能量占比分别大于10%。即25%区吸能量占比为[10%,30%],第一层和第三层的吸能量占比分别为[10%,30%],第二层的吸能量占比为[40%,80%]。假设初步设定正面碰撞能量的分解方案如图2a所示,记为Q(见公式(1))。则每个吸能子空间的吸能量占比qij需要满足公式(2)的条件,其中,i为前舱吸能总空间横向分解的个数,j为吸能总空间垂向分解的个数。其中,QY2为25%区占比,QZ1为第一层占比,QZ2为第二层占比,QZ3为第三层占比。
碰撞波形的纵向分解初步是根据发动机前后的吸能空间划分为两段,在确定波形的横向和垂向分解比例后,需要根据吸能构件的设置和变形模式进一步做详细分解。
2)吸能空间、工况偏置条件、碰撞速度
吸能空间需要根据所研究的车辆的结构特点进行带入研究计算;工况偏置条件以及碰撞速度与选取的碰撞工况实验条件有关,如MPDB工况为50%重叠率偏置,车与带壁障的移动小车之间以64km/h的速度相互碰撞。
2、车辆前端结构三维解析模型的建立
车辆正面碰撞吸能总空间经过三维分解之后,每一个吸能子空间都对应一段分解后的碰撞波形,为了便于区分,记为“子波形”。在形状上与原波形相同,但是波形幅值不同,为原波形幅值与各子空间吸能量占比的乘积。位移域的碰撞波形可以看作是车辆前端结构的等效比刚度,同理,分解到各吸能子空间的子波形则可以看作是该空间的等效比刚度。
在此基础上搭建车辆前端结构三维解析模型,如图3所示,以公式(1)中给出的分解方案为例,将车辆前端每一个吸能子空间都简化为刚度为kij的弹簧,分别与图2a中的子空间吸能量占比qij相对应,计算公式为(3),将车辆简化为质量块。
公式(3)中,av为碰撞波形,dv为车辆位移,M为整车质量,i为前舱吸能总空间横向分解的个数,j为吸能总空间垂向分解的个数。
相比于将车身前部结构刚度直接简化为一根弹簧的简化方式,本发明所提出的三维解析模型更加接近真实情况。并且,本发明所提出的三维解析模型除了可以根据振动方程求解出MPDB工况的OLC和MD两个评价指标外,还可以求解壁障均匀性指标SD以及SOB工况车辆响应。
3、基于MPDB工况的模型求解
MPDB碰撞工况中的壁障在纵向上呈现线性刚度特性,而在横向和垂向上壁障的材料是均匀分布的。因此壁障的吸能空间划分可以在横向和垂向与车辆相对应,纵向不再分段且每个空间的刚度相同。搭建MPDB碰撞工况的三维模型如图4所示(计算时忽略车辆与壁障在碰撞过程中的旋转作用)。
图4中,将车辆和壁障简化为两个质量块并做相对运动,M和MB分别为车辆和壁障的质量,车辆前端的吸能子空间分别简化为刚度为kij的数根弹簧,以同样的划分形式对壁障前端结构进行划分,简化为同等数量(与车辆前端的吸能子空间的数量相同)刚度为kBij的弹簧,dBV为二者的相对位移,dij为车辆局部变形,dBij为壁障局部变形。其中,kBij的值等于壁障的总刚度除以弹簧(吸能子空间)的数量。图4根据机械原理可以简化为图5所示的等效单自由度模型,其中,veq为等效速度,Meq为等效质量,Keq为等效刚度。在碰撞过程中,车辆和壁障的局部刚度kij和kBij先串联获得每个吸能子空间的等效刚度keqij,之后每个吸能子空间的等效刚度并联,获得车辆和壁障整体的等效刚度Keq,具体计算方式如公式(4)所示。
模型的等效质量为Meq,等效速度为veq,如公式(5)所示。
等效单自由度模型的振动方程为:
式中,t为时间,f(t)为等效质量Meq的振动响应,ωeq为系统固有频率。
系统的固有频率ωeq为:
求解振动方程(7)的振动响应如下所示:
壁障和车辆在碰撞过程中做协同振动,振动频率与公式(10)相同,但是振动幅值不同。壁障和车辆的振动响应如下:
式中,dB和dv分别为壁障和车辆碰撞过程中的变形,aB和av分别为壁障和车辆的加速度,vB和vv分别为壁障和车辆的速度。
根据车辆和壁障的运动响应,进一步推导MPDB碰撞工况的三个兼容性评价指标:
壁障的OLC值计算如图6所示。图中实线为壁障的速度曲线,是已知的曲线;虚线为等效的乘员速度曲线,是待求解的曲线。OLC定义为从乘员相对于车辆移动0.065m的时间到乘员相对于车辆移动总计0.3m的时间为止的乘员速度的斜率。
在图6中,虚线与实线之间的面积即为乘员与壁障的相对位移,碰撞开始时,乘员不受约束以匀速运动到t1时刻,乘员与壁障的相对位移S1达到0.065m。此后乘员受到约束作用,开始以恒定的加速度做减速运动,直至t2时刻乘员与壁障速度相同,二者以相同速度继续运动,此时乘员与壁障的相对位移为S1+S2等于0.3m,即S2=0.235m。则乘员减速过程的斜率为OLC值,相当于乘员的加速度。其中,乘员载荷标准OLC的高性能指标为25g,低性能指标为40g,g为一个重力加速度。
OLC通过以下方程计算求得:
v0-OLC×(t2-t1)=vB(t2) (17)
式中:v0为壁障的初速度,t1为乘员运动阶段向前移动65mm的时间,t2为乘员约束阶段向前移动235mm时的时间(在运动阶段之后),即300mm为乘员的运动距离,v0为壁障初速度,也是乘员的初速度;vB(t)为壁障的速度,它可以通过壁障加速度aB(t)积分获得(这里描述的乘员运动是假设壁障驾驶位有乘员存在,便于计算乘员负载指数,用以评价兼容性)。
计算OLC、t1、t2等数值:
(1)依据公式(15)求出时刻t1
公式(15)的含义就是图6中阴影部分的面积为0.065时,求这时的t1值。
(2)求解时刻t2
将公式(15)+(16)结果如下:
公式(18)可变为
根据公式(19)及t1的值求出t2。
车辆和壁障的相对运动位移dBv为:
dBv=dB+dv (20)
壁障的均匀性指标SD的计算公式为:
式中m为前舱吸能子空间的个数。
壁障的最大变形MD为:
dBmax=max(dBij) (23)
获得以上相关评估数值后,相容性修正分数Mcompat计算式如下:
Mcompat=-2·h-MBO(OLC<25g) (24)
Mcompat=-2·OLC/15+10/3-h·((4·OLC/10-8)-(2·OLC/15-10/3))-MBO(25g≤OLC≤40g,Mcompat的分数极限为-8分) (25)
Mcompat=-2-6·h-MBO(当OLC>40g时,Mcompat的分数极限为-8分) (26)
其中,h为均匀系数,当SD<50时,h=0;当50≤SD≤150时,h=(SD-50)/100;当SD≥150时,h=100%。MBO为触底罚分,最终计算得到的Mcompat直接从MPDB正面碰撞工况的总分中扣除。
4、基于SOB工况的模型求解
在SOB工况中,车辆与壁障的解析模型如图7所示。车辆在与刚性壁障接触后沿纵向做减速运动,同时绕其与壁障的接触点做旋转运动。在研究中定义发生以下两种情况之一,则碰撞过程结束:1)车辆的速度减为0的时刻,车辆不再发生结构变形;2)车辆在Y向位移≥25%车宽,车辆将擦过壁障,同样不发生结构变形。因此车辆在SOB工况碰撞过程中的运动响应可以从两方面考虑,即车辆的减速运动和旋转运动。
在图7中,O为车辆与壁障的接触点,也是车辆的旋转中心,v0为碰撞速度,M为整车质量,kij为车辆刚度。由于该工况车辆与壁障是25%重叠,所以主要是25%区的结构发生变形在碰撞过程中吸能,该区域等效刚度为k25%,计算方法如公式(27)所示。
k25%=∑ki1=∑ki4 (27)
车辆在X轴方向做单自由度自由振动,振动方程如下所示:
将公式(29)带入公式(28)中,有:
SOB工况的初始条件为:
t=0,x=0
解得车辆的运动响应为:
车辆在SOB碰撞工况中的受力分析如图8所示。车辆会在碰撞反力F的作用下绕O点做旋转运动,车辆质心到旋转中心的距离为旋转半径r。在车辆做减速运动时,可以求得其纵向位移为xv,则质心到旋转中心的纵向向距离为1/2L-xv,车辆在做旋转运动时,增大质心的横向位移y,则质心到旋转中心的横向距离为25%B+y。由此车辆的旋转半径r的计算公式为(33)。
车辆的旋转角度θ计算公式为(34)。
碰撞反力F垂直于旋转半径的分力为F',车辆受到的转矩为Me,计算公式如(35)所示。
Me=F'·r=F·sin(θ)·r=F·(y+25%B) (35)
已知车辆的转动惯量为J,进一步求得车辆旋转的角加速度β为:
Me=F'·r=J·β (36)
根据车辆的旋转角加速度积分可以获得旋转角度θ,进一步可以求得车辆的横向位移y为:
y=rsinθ-25%×B (37)
到此,车辆在SOB工况中的纵向、横向位移和车辆的加速度、速度响应都已获得。根据碰撞过程结束的判别条件,如果vv≤0或者y≥25%B,则碰撞过程结束。该过程中车辆纵向位移的最大值即为车辆在碰撞过程中的结构变形量,已知防撞梁前端到防火墙的距离为L1,则乘员舱的侵入量D25%max为:
D25%max=max(xv)-L1 (38)
至此,所要求得的车辆MPDB工况下的兼容性评价指标:乘员载荷标准OLC、壁障表面击穿情况、壁障的均匀性指标和兼容性罚分;SOB工况下的乘员舱侵入量全部求解完成。若MPDB工况兼容性指标与SOB工况乘员舱侵入量指标均满足目标需求,则输出对应的能量分解方案作为最终分解方案,若不满足或不完全满足,则需要重新进行能量分解并计算,直至所有指标均满足要求,输出最终分解方案。
实施例
本实施例利用上述方法对所选取的车型进行计算,步骤如下:
1、输入参数
1)根据经验以及车身前端结构特点对FRB工况碰撞波形进行分解
本实施例所选取的车辆前端结构的三维分解方案如图2a-2b所示,为横向50%+25%区,垂向三层,纵向二段。
初步设定正面碰撞能量的分解方案如图2a所示,记为Q(见公式(1))。则每个吸能子空间的吸能量占比qij需要满足公式(2)的条件,其中i为前舱吸能总空间横向分解的个数,j为吸能总空间垂向分解的个数。其中,QY2为25%区占比,QZ1为第一层占比,QZ2为第二层占比,QZ3为第三层占比。
碰撞波形的纵向分解初步是根据发动机前后的吸能空间划分为两段(如图2b所示),在确定波形的横向和垂向分解比例后,需要根据吸能构件的设置和变形模式进一步做详细分解。
MPDB工况:本次选取的车型为中型车,质量M为1400kg;车辆前舱变形空间为0.644m;车辆前端结构的等效刚度K是根据FRB工况的碰撞波形计算的,波形图如图9所示;该工况下的碰撞速度v0为50km/h其前端结构的刚度分解方案为横向50%+25%区,垂向三层,记为“4×3”分解,如图2a所示。该车型的分解方案Q如公式(39)所示。
SOB:本次选取的某车型的质量为1365kg,前舱吸能空间为0.68m。车辆在FRB工况的碰撞波形如图10所示。车辆前端结构的刚度分解方案如公式(40)所示。
2)吸能空间、工况偏置条件、碰撞速度
MPDB工况为车辆与壁障以50km/h的碰撞速度相撞,壁障与车的重叠量为50%,台车质量为1400kg。
SOB工况为车辆以64km/h的速度撞击刚性壁障,壁障与车的重叠率为25%,壁障的圆角半径R=150mm。
2、车辆前端结构三维解析模型的建立
车辆正面碰撞吸能总空间经过三维分解之后,每一个吸能子空间都对应一段分解后的碰撞波形,为了便于区分,记为“子波形”。在形状上与原波形相同,但是波形幅值不同,为原波形与各子空间吸能量占比的乘积。位移域的碰撞波形可以看作是车辆前端结构的等效比刚度,同理,分解到各吸能子空间的子波形则可以看作是该空间的等效比刚度。
在此基础上搭建车辆前端结构三维解析模型。模型中各吸能子空间的等效比刚度用弹簧刚度表示,车辆用质量块表示。以公式(1)中给出的分解方案为例,搭建的车辆前端结构三维解析模型如图3所示。
图3中,弹簧刚度记为kij,与图2a中的子空间吸能量占比qij相对应,计算公式为(3)。
相比于将车身前部结构刚度直接简化为一根弹簧的简化方式,本发明所提出的三维解析模型更加接近真实情况。并且,本发明所提出的三维解析模型除了可以根据振动方程求解出MPDB工况的OLC和MD两个评价指标外,还可以求解壁障均匀性指标SD以及SOB工况车辆响应。
3、基于MPDB工况的模型求解
MPDB碰撞工况中的壁障在纵向上呈现线性刚度特性,而在横向和垂向上壁障的材料是均匀分布的。因此壁障的吸能空间划分可以在横向和垂向与车辆相对应,纵向不再分段且每个空间的刚度相同。搭建MPDB碰撞工况的三维模型如图4所示(计算时忽略车辆与壁障在碰撞过程中的旋转作用)。
图4中,车辆和壁障做相对运动,M和MB分别为车辆和壁障的质量,dBV为二者的相对位移,dij为车辆局部变形,dBij为壁障局部变形,kij为车辆的局部刚度,kBij为壁障局部刚度。图4根据机械原理可以简化为图5所示的等效单自由度模型。其中,veq为等效速度,Meq为等效质量,Keq为等效刚度。在碰撞过程中,车辆和壁障的局部刚度kij和kBij先串联获得每个吸能子空间的等效刚度keqij,之后每个吸能子空间的等效刚度并联,获得车辆和壁障整体的等效刚度Keq,具体计算方式如公式(4)所示。
模型的等效质量为Meq,等效速度为veq,由公式(5)计算得。
将Keq和Meq带入公式(8)得到系统的固有频率ωeq为:
由公式(7)求得简化模型的振动响应:
由公式(11)、(12)、(13)和(14)得到壁障和车辆的振幅系数u1、u2,壁障和车辆的位移dB和dv,壁障和车辆的加速度aB和av,壁障和车辆的速度vB和vv:
式中,dB和dv分别为壁障和车辆碰撞过程中的变形,aB和av分别为壁障和车辆的加速度,vB和vv分别为壁障和车辆的速度.
根据车辆和壁障的运动响应,进一步推导MPDB碰撞工况的三个兼容性评价指标:
壁障的OLC值计算如图6所示。图中实线为壁障的速度曲线,是已知的曲线;虚线为等效的乘员速度曲线,是待求解的曲线。OLC定义为从乘员相对于车辆移动0.065m的时间到乘员相对于车辆移动总计0.3m的时间为止的乘员速度的斜率。
在图6中,虚线与实线之间的面积即为乘员与壁障的相对位移,碰撞开始时,乘员不受约束以匀速运动到t1时刻,乘员与壁障的相对位移S1达到0.065m。此后乘员受到约束作用,开始以恒定的加速度做减速运动,直至t2时刻乘员与壁障速度相同,二者以相同速度继续运动,此时乘员与壁障的相对位移为S1+S2等于0.3m,即S2=0.235m。则乘员减速过程的斜率为OLC值,相当于乘员的加速度。其中,乘员载荷标准OLC的高性能指标为25g,低性能指标为40g,g为一个重力加速度。
OLC通过以下方程计算求得:
v0-OLC×(t2-t1)=vB(t2)
式中,v0为壁障的初速度,t1为乘员运动阶段向前移动65mm的时间,t2为乘员约束阶段向前移动235mm时的时间(在运动阶段之后),即300mm为乘员的运动距离,v0为壁障初速度,也是乘员的初速度;vB(t)为壁障的速度,它可以通过壁障加速度aB(t)积分获得(这里描述的乘员运动是假设壁障驾驶位有乘员存在,便于计算乘员负载指数,用以评价兼容性)。
根据求得的壁障响应以及公式(15)~(19)联立求解可得到OLC的值,本例中OLC的值为27.33g。
由公式(20)、(21)、(23)求解得到壁障的最大变形MD为:
dBmax=max(dBij)=0.62m,小于车辆前舱变形空间为0.644m;所以,壁障未触底。
由公式(22)求解得到壁障的均匀性指标SD为:
获得以上相关评估数值后,再计算相容性修正分数Mcompat,因为SD=98.3mm,则均匀系数h取0.483,由公式(25)计算可得:
Mcompat=-2·OLC/15+10/3-h·((4·OLC/10-8)-(2·OLC/15-10/3))-MBO
=-2·27.33/15+10/3-0.483·((4·27.33/10-8)-(2·27.33/15-10/3))
=-1.577
所以该车的计算罚分为1.577分。
到此,所要求得的MPDB工况兼容性评价指标:试验后壁障变形情况、乘员载荷标准OLC和壁障表面击穿情况以及罚分情况全部求解完成。
由于整车质量对碰撞系统的动力学响应影响较大,所以本次模型验证选取了不同质量的三种车型作为基础数据,三种车型的质量分别为大型车质量1800kg、中型车质量1400kg、小型车质量1100kg,加速度曲线分别如图9中曲线所示,分别计算三种车型的均匀性指标、壁障的OLC值和最大变形量MD以及罚分情况,以验证本发明提供的汽车碰撞解析方法的准确性。本实施例已公开了质量M=1400kg的中型车的计算过程,其他两个车型的计算同上,分解方案如公式(41)所示
并将计算结果与本次验证选取的三个不同质量的车型的仿真实验数据(实际数据)进行对比,其结果与实际数据对比如表1所示:
表1 MPDB工况三个评价指标的验证
将公式(41)中的分解方案输入到解析模型中计算MPDB工况的兼容性指标并计算兼容性罚分。如表1所示,模型计算结果与原数据的误差均在10%左右。
4、基于SOB工况的模型求解
在SOB工况中,车辆与壁障的解析模型如图7所示。车辆在与刚性壁障接触后沿纵向做减速运动,同时绕其与壁障的接触点做旋转运动。在研究中定义发生以下两种情况之一,则碰撞过程结束:1)车辆的速度减为0的时刻,车辆不再发生结构变形;2)车辆在Y向位移≥25%车宽,车辆将擦过壁障,同样不发生结构变形。因此车辆在SOB工况碰撞过程中的运动响应可以从两方面考虑,即车辆的减速运动和旋转运动。
本实施例以质量为1365kg的某车型在SOB工况的有限元计算结果为基础数据,验证该解析模型计算25%碰撞工况中乘员舱上部最大侵入量的准确性,对其进行FRB与SOB碰撞工况的仿真,并取乘员舱上部4个侵入量测量点的最大侵入量值为34.2cm。
在图7中,O为车辆与壁障的接触点,也是车辆的旋转中心,v0为碰撞速度,M为整车质量,kij为车辆刚度。由于该工况车辆与壁障是25%重叠,所以主要是25%区的结构发生变形在碰撞过程中吸能,该区域等效刚度为k25%。
系统的固有频率ω为:
由公式(32)可求得简化模型的振动响应:
将图10中的碰撞波形转化为位移域的10台阶波,并与公式(40)中的分解方案一起带入解析模型中计算,获得车辆纵向的最大位移为1.0019m。
根据求得的简化模型的振动响应,利用公式(32)~(37)判断碰撞结束时间,并计算碰撞过程中的乘员舱的侵入量,由公式(38)计算可得:
D25%max=max(xv)-L1=1.0019-0.68=0.3219m
即解析模型计算得到的乘员舱最大侵入量值为32.19cm,与该车型的乘员舱上部侵入量仿真结果的最大值对比的误差在5.8%左右。
由此可见,该模型计算MPDB工况的兼容性指标以及SOB工况的乘员舱侵入量指标准确性均较好。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域,对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的图例。
Claims (10)
1.一种协同正面碰撞多工况的汽车前端结构能量管理方法,其特征在于,包括如下步骤:
步骤一、根据总布置阶段的车辆硬点尺寸确定车辆前端吸能总空间;将所述车辆前端吸能总空间划分为多个吸能子空间,并将FRB工况碰撞原波形分解为与所述吸能子空间一一对应的子波形,建立车辆前端结构三维解析模型;
其中,每个所述子波形的幅值为FRB工况碰撞原波形的幅值与其对应的吸能子空间的吸收能量占比的乘积;
步骤二、基于MPDB工况,对所述前端结构三维解析模型进行求解,得到车辆和壁障的运动响应,确定碰撞过程中壁障和车辆的变形、壁障和车辆的加速度及壁障和车辆的速度;并根据所述运动响应得到乘员负载指数、壁障表面击穿情况、壁障的均匀性指标和兼容性罚分;以及
基于SOB工况,对所述前端结构三维解析模型进行求解,得到碰撞结束时车辆乘员舱的形变量。
3.根据权利要求2所述的协同正面碰撞多工况的汽车前端结构能量管理方法,其特征在于,在所述步骤二中,对所述前端结构三维解析模型进行求解之前,还包括:
基于MPDB工况,将壁障在横向和垂向划分为与车辆前端的吸能子空间相对应的多个区域,将所述多个区域分别简化为弹簧;并根据机械原理将壁障和车辆的碰撞模型简化为等效单自由度模型。
5.根据权利要求4所述的协同正面碰撞多工况的汽车前端结构能量管理方法,其特征在于,在所述步骤二中,兼容性罚分的计算公式为:
当OLC<25g时,Mcompat=-2·h-MBO;
当25g≤OLC≤40g时,
Mcompat=-2·OLC/15+10/3-h·((4·OLC/10-8)-(2·OLC/15-10/3))-MBO;
当OLC>40g时,Mcompat=-2-6·h-MBO;
式中,h为均匀系数,MBO为触底罚分。
6.根据权利要求5所述的协同正面碰撞多工况的汽车前端结构能量管理方法,其特征在于,所述均匀系数,根据壁障的均匀性指标确定;其中:
当SD<50时,h=0;
当50≤SD≤150时,h=(SD-50)/100;
当SD≥150时,h=100%;
式中,SD为壁障的均匀性指标。
8.根据权利要求6或7所述的协同正面碰撞多工况的汽车前端结构能量管理方法,其特征在于,在所述步骤二中,车辆乘员舱的形变量通过如下公式进行计算:
D25%max=max(xv)-L1;
式中,D25%max为车辆乘员舱的形变量,max(xv)为碰撞过程结束时车辆的纵向最大形变量,L1为防撞梁前端到防火墙的距离。
9.根据权利要求8所述的协同正面碰撞多工况的汽车前端结构能量管理方法,其特征在于,将通过所述三维解析模型计算得到的车辆的最大纵向位移作为车辆的纵向最大形变量。
10.根据权利要求9所述的协同正面碰撞多工况的汽车前端结构能量管理方法,其特征在于,在所述步骤二中,基于SOB工况判断碰撞过程是否结束的方法为:
当vv≤0或者y≥25%B时,判定碰撞过程结束;
其中,vv为车辆速度,y为车辆的横向位移,B为车辆宽度;
y=rsinθ-25%×B;
式中,θ为车辆的旋转角度,r为车辆的旋转半径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110422464.1A CN112948983B (zh) | 2021-04-20 | 2021-04-20 | 一种协同正面碰撞多工况的汽车前端结构能量管理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110422464.1A CN112948983B (zh) | 2021-04-20 | 2021-04-20 | 一种协同正面碰撞多工况的汽车前端结构能量管理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112948983A true CN112948983A (zh) | 2021-06-11 |
CN112948983B CN112948983B (zh) | 2022-04-05 |
Family
ID=76233015
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110422464.1A Expired - Fee Related CN112948983B (zh) | 2021-04-20 | 2021-04-20 | 一种协同正面碰撞多工况的汽车前端结构能量管理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112948983B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113361019A (zh) * | 2021-07-19 | 2021-09-07 | 中国第一汽车股份有限公司 | 一种基于正面耐撞的前端结构多工况设计方法 |
CN114001974A (zh) * | 2021-09-23 | 2022-02-01 | 中汽研汽车检验中心(天津)有限公司 | 一种评价车辆前部构件对乘员舱碰撞响应贡献的方法 |
CN114047000A (zh) * | 2021-10-28 | 2022-02-15 | 中汽研汽车检验中心(天津)有限公司 | 侧面碰撞蜂窝铝壁障前端刚度的确定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106709151A (zh) * | 2016-11-28 | 2017-05-24 | 吉林大学 | 汽车正面25%重叠率碰撞的车身前端结构设计方法 |
CN107256289A (zh) * | 2017-05-18 | 2017-10-17 | 吉林大学 | 汽车碰撞简化参数化有限元模型的建立方法 |
CN107992668A (zh) * | 2017-11-28 | 2018-05-04 | 吉林大学 | 一种基于双梯形波的碰撞波形概念设计方法 |
JP2019191041A (ja) * | 2018-04-26 | 2019-10-31 | 三菱重工機械システム株式会社 | 自動車衝突模擬試験の試験条件設定方法及び自動車衝突模擬試験の試験条件設定システム |
CN111855234A (zh) * | 2020-09-01 | 2020-10-30 | 东风汽车集团有限公司 | 一种车辆前端结构子系统评价试验系统及方法 |
CN112519710A (zh) * | 2020-12-23 | 2021-03-19 | 凌云工业股份有限公司 | 提高车辆正面碰撞兼容性和小偏置碰撞安全性防撞梁结构 |
-
2021
- 2021-04-20 CN CN202110422464.1A patent/CN112948983B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106709151A (zh) * | 2016-11-28 | 2017-05-24 | 吉林大学 | 汽车正面25%重叠率碰撞的车身前端结构设计方法 |
CN107256289A (zh) * | 2017-05-18 | 2017-10-17 | 吉林大学 | 汽车碰撞简化参数化有限元模型的建立方法 |
CN107992668A (zh) * | 2017-11-28 | 2018-05-04 | 吉林大学 | 一种基于双梯形波的碰撞波形概念设计方法 |
JP2019191041A (ja) * | 2018-04-26 | 2019-10-31 | 三菱重工機械システム株式会社 | 自動車衝突模擬試験の試験条件設定方法及び自動車衝突模擬試験の試験条件設定システム |
CN111855234A (zh) * | 2020-09-01 | 2020-10-30 | 东风汽车集团有限公司 | 一种车辆前端结构子系统评价试验系统及方法 |
CN112519710A (zh) * | 2020-12-23 | 2021-03-19 | 凌云工业股份有限公司 | 提高车辆正面碰撞兼容性和小偏置碰撞安全性防撞梁结构 |
Non-Patent Citations (2)
Title |
---|
JUNYUAN ZHANG 等: "The vehicle frontal impact safety evaluation and design method based on the coupling effect between crash pulse and restraint system characteristic", 《INTERNATIONAL JOURNAL OF VEHICLE DESIGN (IJVD)》 * |
李林峰 等: "基于25%小偏置正面碰撞的某乘用车前端结构改进设计", 《中国机械工程》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113361019A (zh) * | 2021-07-19 | 2021-09-07 | 中国第一汽车股份有限公司 | 一种基于正面耐撞的前端结构多工况设计方法 |
CN114001974A (zh) * | 2021-09-23 | 2022-02-01 | 中汽研汽车检验中心(天津)有限公司 | 一种评价车辆前部构件对乘员舱碰撞响应贡献的方法 |
CN114001974B (zh) * | 2021-09-23 | 2023-09-22 | 中汽研汽车检验中心(天津)有限公司 | 一种评价车辆前部构件对乘员舱碰撞响应贡献的方法 |
CN114047000A (zh) * | 2021-10-28 | 2022-02-15 | 中汽研汽车检验中心(天津)有限公司 | 侧面碰撞蜂窝铝壁障前端刚度的确定方法 |
CN114047000B (zh) * | 2021-10-28 | 2024-01-02 | 中汽研汽车检验中心(天津)有限公司 | 侧面碰撞蜂窝铝壁障前端刚度的确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112948983B (zh) | 2022-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112948983B (zh) | 一种协同正面碰撞多工况的汽车前端结构能量管理方法 | |
CN106709151B (zh) | 汽车正面25%重叠率碰撞的车身前端结构设计方法 | |
Christensen et al. | Nonlinear optimization of vehicle safety structures: Modeling of structures subjected to large deformations | |
Reichert et al. | Validation of a toyota camry finite element model for multiple impact configurations | |
Yadav et al. | Investigations into dynamic response of automobile components during crash simulation | |
Alekseev et al. | Validation of euroncap frontal impact of frame off-road vehicle–road traffic accident simulation | |
Yıldızhan et al. | Design of improved energy absorbing pads to reduce occupant injuries in vehicle side impact | |
Geigl et al. | Reconstruction of occupant kinematics and kinetics for real world accidents | |
Aşkar et al. | Crash analysis and size optimization of a vehicle’s front bumper system | |
Tahan et al. | The Effect of Reduced Mass on Frontal Crashworthiness | |
Yu et al. | Optimum design of vehicle frontal structure and occupant restraint system for crashworthiness (a multilevel approach using SDSS) | |
ZHANG et al. | Lightweight design of automotive front end material-structure based on frontal collision | |
Jenefeldt et al. | A methodology to assess frontal stiffness to improve crash compatibility | |
CN112818473B (zh) | 一种用于求解汽车mpdb碰撞工况系统动力学响应的解析方法 | |
Choi et al. | Design methodology to reduce the chest deflection in US NCAP and EURO NCAP tests | |
Woitsch et al. | Influence of pitching and yawing during frontal passenger vehicle crash tests on driver occupant's kinematics and injury | |
Trella et al. | Application of derived characteristics from dynamic test data for simulation of car-to-car side impacts using a lumped mass approach | |
Brewer | Effects of angles and offsets in crash simulations of automobiles with light trucks | |
Noorsumar et al. | Modeling of Modified Vehicle Crashworthiness using a Double Compound Pendulum | |
Wang et al. | Crash Energy Management of Vehicle Front-end Structures Considering Multiple Conditions: Modelling and Solution Method | |
Li et al. | The effects of curtain airbag on occupant kinematics and injury index in rollover crash | |
Barbat et al. | A comparative analysis of vehicle-to-vehicle and vehicle-to-rigid fixed barrier frontal impacts | |
Jeyakumar et al. | Improvement of the frontal structure of a bus for crash accidents | |
Vysochinskiy et al. | Development and Extended Validation of a Lumped Parameter Prediction Model for Analysing Injury Parameters in a Vehicle Crash | |
Reichert et al. | Development of a 2015 Mid-Size Sedan Vehicle Model |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220405 |