CN107577843A - 碰撞波形与约束***特性耦合关系的评价方法 - Google Patents
碰撞波形与约束***特性耦合关系的评价方法 Download PDFInfo
- Publication number
- CN107577843A CN107577843A CN201710649471.9A CN201710649471A CN107577843A CN 107577843 A CN107577843 A CN 107577843A CN 201710649471 A CN201710649471 A CN 201710649471A CN 107577843 A CN107577843 A CN 107577843A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- occupant
- constrained system
- collision waveform
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 111
- 238000004458 analytical method Methods 0.000 claims abstract description 26
- 238000012935 Averaging Methods 0.000 claims abstract description 12
- 238000012417 linear regression Methods 0.000 claims abstract description 5
- 230000001133 acceleration Effects 0.000 claims description 63
- 238000000034 method Methods 0.000 claims description 32
- 238000013461 design Methods 0.000 claims description 26
- 238000002474 experimental method Methods 0.000 claims description 15
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 8
- 238000000611 regression analysis Methods 0.000 claims description 8
- 238000010521 absorption reaction Methods 0.000 claims description 7
- 230000003399 chemotactic effect Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 238000000205 computational method Methods 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 230000001419 dependent effect Effects 0.000 claims description 5
- 238000007689 inspection Methods 0.000 claims description 5
- 239000013589 supplement Substances 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000003786 synthesis reaction Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 8
- 238000010276 construction Methods 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 8
- 238000004088 simulation Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 5
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 208000027418 Wounds and injury Diseases 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 208000014674 injury Diseases 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000002633 protecting effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Landscapes
- Air Bags (AREA)
Abstract
本发明公开了碰撞波形与约束***特性耦合关系的评价方法，克服了CAE仿真技术在研究车体结构和乘员约束***特性耦合关系时计算量庞大且费时的问题，步骤：1.基于单自由度模型的乘员响应面的建立：1)定义双梯形波和三线性约束刚度曲线的形状特征参数:2)单自由度模型求解乘员响应面；3)对响应面XY向分别求平均；2.乘员响应与碰撞波形参数相关性分析：1)碰撞波形参数定义补充；2)从乘员响应面中提取出碰撞波形基本参数与Av对应关系，进行线性回归分析；3.碰撞波形与约束***特性耦合关系的评价：1)碰撞波形综合评价指标α的建立；2)约束***综合评价指标β的建立；3)碰撞波形与约束***特性综合评价指标的ao的建立。
Description
技术领域
本发明涉及一种汽车车身结构与约束***耦合作用效果的评价方法，更具体地说，本发明涉及一种碰撞波形与约束***特性耦合关系的评价方法。
背景技术
汽车安全性由车体结构的抗撞性和乘员约束***两方面共同决定。目前整车安全性开发流程主要是在整车结构的抗撞性设计后再进行约束***匹配。随着这两个领域研究的深入，企业陆续从各个方面进行抗撞性结构开发和约束***装置开发。车体抗撞性设计由于处于安全性开发流程的前期阶段，尚不能与乘员伤害建立直接关系，所以主要以较大吸能特性和良好的前端结构刚度等作为设计依据。约束***设计处于安全性开发的后期阶段，有时会出现无论怎样匹配约束***参数乘员伤害总是居高不下的现象。出现这种现象的主要原因是在概念设计阶段对车体结构和乘员约束***特性耦合关系考虑欠缺，致使车体抗撞性设计与约束***设计之间出现脱节。如果在安全性设计初期从耦合关系角度出发，提出车体结构的抗撞性和乘员约束***设计要求，就可以从宏观上对匹配效果进行整体把控。
由于车体结构和乘员约束***特性耦合关系影响参数众多，并且各个参数相互干扰，如果将所有参数水平放入仿真模型，计算量将非常庞大，所以不利于通过仿真方法进行耦合关系研究。国内外学者从碰撞力学理论层面展开了耦合关系研究，得到了一系列重要结论，但尚未有明确的评价指标。进一步加强耦合关系理论研究与工程实践的联系，对于指导车体结构抗撞性设计以及约束***匹配设计有重要意义。
发明内容
本发明所要解决的技术问题是克服了CAE仿真技术在研究车体结构和乘员约束***特性耦合关系时计算量庞大且耗费时间的问题，提供了一种碰撞波形与约束***特性耦合关系的评价方法。
为解决上述技术问题，本发明是采用如下技术方案实现的：所述的碰撞波形与约束***特性耦合关系的评价方法的步骤如下：
1)基于单自由度模型的乘员响应面的建立：
(1)定义双梯形波和三线性约束刚度曲线的形状特征参数:
(2)单自由度模型求解乘员响应面；
(3)对响应面XY向分别求平均：
对乘员响应面向X坐标方向求平均，得到每个碰撞波形对应的所有乘员响应的平均值，用符号Av表示；
对乘员响应面向Y坐标方向求平均，即可得到每个约束***刚度对应的所有乘员响应的平均值，用符号Ar表示；
2)乘员响应与碰撞波形参数相关性分析：
(1)碰撞波形参数的定义补充；
(2)从乘员响应面中提取出碰撞波形基本参数与Av对应关系，进行线性回归分析，并将相关性系数汇总；
3)碰撞波形与约束***特性耦合关系的评价：
(1)碰撞波形综合评价指标α的建立；
(2)约束***综合评价指标β的建立；
(3)碰撞波形与约束***特性综合评价指标的ao的建立。
技术方案中所述的定义双梯形波和三线性约束刚度曲线的形状特征参数是指:
G1为第一阶加速度值，即BC段所对应的加速度值；第一阶加速度的持续时间为tBC，单位为g；
G2为第二阶加速度值，即DE段对应的加速度值,第二阶加速度的持续时间为tDE，单位为g；
s1、s2、s3分别为直线AB、CD以及回弹段EF的斜率；
C1为发动机前端的压溃量，单位为m；
Cmax为车体前端结构的最大动态压溃量，单位为m；
k1为OP段斜率为安全带比刚度，OP段对应安全带线性伸长阶段；
GL为限力加速度；
k2为QR段斜率为安全气囊刚度；
D1为安全带作用力达到限力值时乘员相对于车体的运动位移；
D2为安全气囊开始作用时乘员相对于车体的运动位移。
技术方案中所述的碰撞波形参数的定义补充是指：
1)阶梯比i：定义等效双梯形波的两个阶高度G2与G1之比为阶梯比；
2)宽度比w：宽度比为发动机前端的压溃量C1与最大动态压溃量Cmax之比；
3)平均刚度KAE：平均刚度为等效双梯形波的第二阶高度G2与最大动态压溃量Cmax之比；
4)能量密度比a：能量密度比为发动机前端吸能量与车体前端结构总吸能量的比；
式中:v0为碰撞速度，单位m/s；G1为第一台阶的高度，单位为g；tC为发动机碰撞时刻，单位s；s1为直线AB段的斜率，单位g/s。
5)波形形心(to，Go)：波形形心是指加速度-时间曲线与时间轴围成的几何形状的形心，形心的横坐标称为形心时刻to，形心纵坐标称为形心加速度Go；
式中:t为时间，单位s；为车体加速度，单位g。
技术方案中所述的单自由度模型求解乘员响应面是指：
(1)求解该响应面的单自由度模型将车辆与乘员分别简化为集中质量块M和m,车体结构在碰撞中的压溃过程简化为弹簧刚度K的压缩过程，并近似为等效双梯形波；约束***在碰撞中的变形过程简化为弹簧刚度k的压缩过程，并近似为三线性约束刚度；乘员在两个弹簧振动***的共同振动作用下做向前的减速运动，因此乘员响应的大小由两个弹簧振动***的耦合关系决定；
(2)根据设定好的形状参数范围和步长，见表1，可得到3204个不同的双梯形波形与3136个三线性约束刚度；将获得的双梯形波形与约束***刚度进行一对一组合，利用已有的单自由度模型迭代算法快速对每个组合进行乘员响应求解，得到共计约1000万个不同双梯形波形与不同约束***刚度的匹配结果，如表2，表中a0为乘员加速度峰值；
表1形状参数的范围和步长
表2不同碰撞波形参数与约束***刚度参数的匹配结果
(3)以不同双梯形波形从1到3204编号作为X坐标，以不同约束***刚度从1到3136编号作为Y坐标，以乘员加速度峰值作为Z坐标，可绘制在不同碰撞波形和不同约束***刚度组合情况下的乘员加速度峰值响应面，简称为乘员响应面。
技术方案中所述的碰撞波形综合评价指标α的建立是指：
本步骤主要建立波形平均响应Av与碰撞波形参数之间的直接量化关系，确定函数P，实现碰撞波形初步的单独评价，本技术方案采用多元回归的方法近似构造Av与碰撞波形参数的函数关系P；
碰撞波形参数当中Go、tE、KAE、G2与乘员响应的相关性较高，但考虑到碰撞波形基本参数的重要性，所以本技术方案选取G1、G2、tC、tE、Go、KAE这六个碰撞波形参数进行多元回归分析；考虑到波形参数的复杂性以及相互之间的关系，本技术方案构建了多元二次回归模型，为方便描述将这六个参数分别用x1、x2……x6六个变量表示，对应因变量为Av，多元二次回归模型表达式如下：
其中a0为常数项，a1、a2……a6为一次项系数，a(1,1),a(1,2)……a(6,6)为二次项系数。通过回归分析剔除影响较小的项目，其中tC项被完全剔除，剩下对回归结果影响较大的7个项目，各个项目以及对应的系数如表4所示。
表4碰撞波形参数二次回归对应项目和系数
该多元二次回归模型的分析结果如表5所示，其中判定系数为接近1，F统计量较大，检验的P值为0，说明回归结果较好；
表5碰撞波形参数二次回归模型的分析结果
Av代表了碰撞波形的平均水平，通过对比不同碰撞波形的Av大小可以实现碰撞波形的优劣评价，为了使Av无量纲化并且与Ar之间同趋化，本技术方案采用min-max标准化方法，对Av作标准化处理，得到无量纲的碰撞波形评价指标，定义为“碰撞波形综合评价指标”，用符号α表示，碰撞波形综合评价指标α的计算方法为：
根据乘员响应面确定Av最大值为79.28g，最小值为36.70g，根据公式(7)和(8)，得到α计算公式为：
α分布在0到1之间，从乘员加速度角度考虑，该值越小对应的碰撞波形越好。然而当碰撞波形参数超出本技术方案确定的常用范围时会出现α大于1或者小于0的情况。
技术方案中所述的约束***综合评价指标β的建立是指：
本步骤主要建立Ar与约束***刚度参数之间的直接量化关系，确定函数Q，实现约束***特性初步的单独评价；本步骤同样采用多元回归的方法近似构造Ar与约束***刚度参数的函数关系Q；
由于约束***刚度参数较少，本技术方案选取所有约束***刚度参数k1、k2、GL、D2构建多元二次回归模型，为方便描述将这四个参数分别用y1、y2、y3、y4四个变量表示，对应因变量为Ar，多元二次回归模型表达式如下：
其中b0为常数项，b1、b2、b3、b4为一次项系数，b(1,1)、b(1,2)……b(4,4)为二次项系数；通过回归分析剔除影响较小的项目，其中参数D2被完全剔除，剩下对回归模型影响较大的7个项目，各个项目以及对应的系数如表6所示；
表6约束***刚度参数二次回归对应项目和系数
该多元二次回归模型的分析结果如表7所示，其中判定系数接近1，F统计量较大，检验的P值为0，说明回归结果较好；
表7约束***刚度参数二次回归模型的分析结果
Ar代表了约束***刚度的平均水平，通过对比不同约束***刚度的Ar大小可以实现约束***刚度的优劣评价，为了使Ar无量纲化并且与Av之间同趋化，本文采用min-max标准化方法，对Ar作标准化处理，得到无量纲约束***特性的评价指标，定义为“约束***综合评价指标”，用符号β表示，其计算方法为：
根据乘员响应面确定Ar最大值为69.12g，最小值为44.60g。根据公式(10)和(11)，得到β计算为：
β分布在0到1之间，从乘员加速度角度考虑，该值越小对应的约束***刚度越好；然而当约束***刚度参数超出本技术方案确定的常用范围时可能会出现β大于1或者小于0的情况。
技术方案中所述的碰撞波形与约束***特性综合评价指标的ao的建立是指：
本步骤主要建立乘员加速度响应ao和碰撞波形综合评价指标α与约束***综合评价指标β的直接量化关系，从而根据碰撞波形以及约束***特性初步的单独评价结果实现整体耦合关系评价，由于α和β分别由Av和Ar标准化处理得到，所以ao和Av与Ar之间的递增函数R可等价于ao和α与β之间的递增函数关系R；可将乘员响应面按照α与β的对应关系，绘制出α-β-ao响应面，为了便于观察，用网格代替响应面，整个曲面呈现出较好的递增趋势；
本技术方案采用曲面拟合的方法近似构造ao和α与β之间的递增函数关系R，以α与β为自变量，对α-β-ao响应面进行二次曲面拟合，曲面拟合如公式(13)所示，其中c0、c1、c2、c3、c4、c5为各项系数，通过最小二乘法得到各项系数值如表8所示；
ao＝c0+c1α+c2β+c3α2+c4αβ+c5β2 (13)
表8α-β-ao响应面二次曲面拟合结果
由于α与β的平方项系数接近0，所以曲面拟合公式可简化为：
ao＝280.5+335.1α+161.8β+168αβ (14)
该拟合结果的决定系数R2达到0.986，拟合精度较高。
与现有技术相比本发明的有益效果是：
1.本发明所述的一种碰撞波形与约束***特性耦合关系的评价方法可以通过碰撞波形综合评价指标α和约束***综合评价指标β可实现碰撞波形和约束***特性耦合关系评价，并直观地找到碰撞波形及约束***参数设计过程当中出现问题的环节。
2.本发明所述的一种碰撞波形与约束***特性耦合关系的评价方法建立α-β-ao响应面拟合曲面，可以通过碰撞波形综合评价指标α和约束***综合评价指标β估算乘员加速度峰值，简单、方便且与单自由度模型计算的误差控制在10％。
附图说明
下面结合附图对本发明作进一步的说明：
图1为本发明所述的碰撞波形与约束***特性耦合关系的评价方法的流程图；
图2-a为本发明所述的碰撞波形与约束***特性耦合关系的评价方法中车体-乘员单自由度求解模型示意图；
图2-b为本发明所述的碰撞波形与约束***特性耦合关系的评价方法中车体-乘员单自由度求解模型采用的双梯形波示意图；
图2-c为本发明所述的碰撞波形与约束***特性耦合关系的评价方法中车体-乘员单自由度求解模型采用的三线性刚度曲线示意图；
图3-a为本发明所述的碰撞波形与约束***特性耦合关系的评价方法中时间域等效双梯形波形参数示意图；
图3-b为本发明所述的碰撞波形与约束***特性耦合关系的评价方法中位移域等效双梯形波形参数示意图；
图4为本发明所述的碰撞波形与约束***特性耦合关系的评价方法中三线性约束刚度特征参数示意图；
图5为本发明所述的碰撞波形与约束***特性耦合关系的评价方法的原理图；
图6为本发明所述的碰撞波形与约束***特性耦合关系的评价方法的α-β-ao响应面图；
图7为本发明所述的碰撞波形与约束***特性耦合关系的评价方法中提取六个双梯形波形；
图8为本发明所述的碰撞波形与约束***特性耦合关系的评价方法中碰撞波形综合评价指标α与乘员胸部加速度峰值对应关系；
图9为本发明所述的碰撞波形与约束***特性耦合关系的评价方法中约束***综合评价指标β与乘员胸部加速度峰值对应关系；
图10为本发明所述的碰撞波形与约束***特性耦合关系的评价方法中乘员胸部加速度峰值对比关系图。
具体实施方式
下面结合附图对本发明作详细的描述：
参阅图1,本发明所述的碰撞波形与约束***特性耦合关系的评价方法的步骤如下:
1.基于单自由度模型的乘员响应面的建立
1)定义双梯形波和三线性约束刚度曲线的形状特征参数
参阅图3-a至图3-b，本发明所述的碰撞波形与约束***特性耦合关系的评价方法中使用的碰撞波形均为等效双梯形波,图3-a中A、B、C、D、E、F为等效双梯形波的特征点；其中：A是碰撞发生点，其对应的时刻为0；B为前纵梁碰到壁障的点；C为发动机与壁障接触点，其对应的时刻为tC；D为上边梁碰到壁障的点；E为车***移达到最大动态压溃量点，其对应的时刻为tE，tE与最大动态压溃时刻tmax相等；F为碰撞结束点；
G1为第一阶加速度值，即BC段所对应的加速度值；第一阶加速度的持续时间为tBC，单位为g；
G2为第二阶加速度值，即DE段对应的加速度值,第二阶加速度的持续时间为tDE，单位为g；
s1、s2、s3分别为直线AB、CD以及回弹段EF的斜率；
C1为发动机前端的压溃量，单位为m；
Cmax为车体前端结构的最大动态压溃量，单位为m。
参阅图4，本发明所述的碰撞波形与约束***特性耦合关系的评价方法中使用的约束***刚度曲线均为三线性约束刚度曲线。图4中的O、P、Q、R是“三线性约束刚度”的特征点，其中O为碰撞发生点；P点为安全带限力装置作用点；Q点为安全气囊作用开始点；R点为最大相对位移点。
k1为OP段斜率为安全带比刚度，OP段对应安全带线性伸长阶段；
GL为限力加速度；
k2为QR段斜率为安全气囊刚度；
D1为安全带作用力达到限力值时乘员相对于车体的运动位移；
D2为安全气囊开始作用时乘员相对于车体的运动位移。
2)单自由度模型求解乘员响应面
本发明所述碰撞波形与约束***特性耦合关系的评价方法以乘员响应面为数据基础。该响应面是利用单自由度模型迭代算法进行乘员加速度响应快速、准确的求解，得到的不同碰撞波形参数与约束***刚度参数的匹配结果。
(1)参阅图2-a至图2-c，求解该响应面的单自由度模型将车辆与乘员分别简化为集中质量块M和m,车体结构在碰撞中的压溃过程简化为弹簧刚度K的压缩过程，并近似为等效双梯形波；约束***在碰撞中的变形过程简化为弹簧刚度k的压缩过程，并近似为三线性约束刚度；乘员在两个弹簧振动***的共同振动作用下做向前的减速运动，因此乘员响应的大小由两个弹簧振动***的耦合关系决定。
(2)根据设定好的形状参数范围和步长，见表1，可得到3204个不同的双梯形波形与3136个三线性约束刚度。将获得的双梯形波形与约束***刚度进行一对一组合，利用已有的单自由度模型迭代算法快速对每个组合进行乘员响应求解，得到共计约1000万个不同双梯形波形与不同约束***刚度的匹配结果，如表2，表中a0为乘员加速度峰值。
表1形状参数的范围和步长
表2不同碰撞波形参数与约束***刚度参数的匹配结果
(3)以不同双梯形波形从1到3204编号作为X坐标，以不同约束***刚度从1到3136编号作为Y坐标，以乘员加速度峰值作为Z坐标，可绘制在不同碰撞波形和不同约束***刚度组合情况下的乘员加速度峰值响应面，简称为乘员响应面。
3)对响应面XY向分别求平均
对乘员响应面向X坐标方向求平均，得到每个碰撞波形对应的所有乘员响应的平均值，用符号Av表示；对乘员响应面向Y坐标方向求平均，即可得到每个约束***刚度对应的所有乘员响应的平均值，用符号Ar表示。
2.乘员响应与碰撞波形参数相关性分析
1)碰撞波形参数的定义补充如下：
(1)阶梯比i：定义等效双梯形波的两个阶高度G2与G1之比为阶梯比；
(2)宽度比w：宽度比为发动机前端的压溃量C1与最大动态压溃量Cmax之比；
(3)平均刚度KAE：平均刚度为等效双梯形波的第二阶高度G2与最大动态压溃量Cmax之比；
(4)能量密度比a：能量密度比为发动机前端吸能量与车体前端结构总吸能量的比；
式中:v0为碰撞速度，单位m/s；G1为第一台阶的高度，单位为g；tC为发动机碰撞时刻，单位s；s1为直线AB段的斜率，单位g/s。
(5)波形形心(to，Go)：波形形心是指加速度-时间曲线与时间轴围成的几何形状的形心，形心的横坐标称为形心时刻to，形心纵坐标称为形心加速度Go；
式中:t为时间，单位s；为车体加速度，单位g。
2)从乘员响应面中提取出碰撞波形基本参数与Av对应关系，利用Origin作图软件进行线性回归分析，并将相关性系数汇总，如表3所示；
表3碰撞波形参数与Av相关性
由表3可知，Go、tE、KAE、G2与乘员响应的相关性较高。
3.碰撞波形与约束***特性耦合关系的评价
在碰撞波形以及约束***刚度参数相关性分析的基础上，研究碰撞波形、约束***特性以及两者之间耦合关系的量化评价方法。Av虽然可以近似表征碰撞波形的好坏，但其计算过程仍需要约束***参数的介入。建立Av与碰撞波形参数之间的直接量化关系，如公式1函数P所示形式，可以脱离约束***参数仅仅根据碰撞波形参数实现碰撞波形初步的单独量化评价；建立Ar与约束***刚度参数之间的量化关系，如公式2函数Q所示形式，可以脱离碰撞波形仅根据约束***刚度参数实现约束***特性初步的单独量化评价；
Av＝P(G1,G2,tE,tC,...) (4)
Ar＝Q(k1,k2,GL,D1,...) (5)
根据乘员响应面，当Av和Ar值确定时，有唯一的乘员加速度响应ao相对应。由于排序后的乘员响应面整体呈现规律的递增，可认为乘员响应ao与Av和Ar之间存在一种递增函数关系，如公式3函数R所示形式，这样就可以建立起碰撞波形与约束***特性单独的量化评价与总体耦合之间的关系；
ao＝R(Av,Ar) (6)
参阅图5，确定以上函数关系便可建立整个耦合关系评价方法，该评价方法步骤主要包括：
(1)采用多元回归的方法确定函数P，从而根据碰撞波形参数实现碰撞波形初步的单独评价；
(2)采用多元回归的方法确定函数Q，从而根据约束***刚度参数实现约束***特性初步的单独评价；
(3)通过MATLAB软件的曲面拟合功能，确定函数R，从而根据碰撞波形以及约束***特性初步的单独评价结果实现整体耦合关系评价。
1)碰撞波形综合评价指标α的建立
本步骤主要建立波形平均响应Av与碰撞波形参数之间的直接量化关系，确定函数P，实现碰撞波形初步的单独评价。本技术方案采用多元回归的方法近似构造Av与碰撞波形参数的函数关系P。
碰撞波形参数当中Go、tE、KAE、G2与乘员响应的相关性较高，但考虑到碰撞波形基本参数的重要性，所以本技术方案选取G1、G2、tC、tE、Go、KAE这六个碰撞波形参数进行多元回归分析。考虑到波形参数的复杂性以及相互之间的关系，本技术方案构建了多元二次回归模型。为方便描述将这六个参数分别用x1、x2……x6六个变量表示，对应因变量为Av，多元二次回归模型表达式如下：
其中a0为常数项，a1、a2……a6为一次项系数，a(1,1),a(1,2)……a(6,6)为二次项系数。通过回归分析剔除影响较小的项目，其中tC项被完全剔除，剩下对回归结果影响较大的7个项目，各个项目以及对应的系数如表4所示。
表4碰撞波形参数二次回归对应项目和系数
该多元二次回归模型的分析结果如表5所示。其中判定系数为接近1，F统计量较大，检验的P值为0，说明回归结果较好。
表5碰撞波形参数二次回归模型的分析结果
Av代表了碰撞波形的平均水平，通过对比不同碰撞波形的Av大小可以实现碰撞波形的优劣评价。为了使Av无量纲化并且与Ar之间同趋化，本技术方案采用min-max标准化方法，对Av作标准化处理，得到无量纲的碰撞波形评价指标，定义为“碰撞波形综合评价指标”，用符号α表示。碰撞波形综合评价指标α的计算方法为：
根据乘员响应面确定Av最大值为79.28g，最小值为36.70g。根据公式(7)和(8)，得到α计算公式为：
α分布在0到1之间，仅从乘员加速度角度考虑，该值越小对应的碰撞波形越好。然而当碰撞波形参数超出本技术方案确定的常用范围时可能会出现α大于1或者小于0的情况。
2)约束***综合评价指标β的建立
本步骤主要建立Ar与约束***刚度参数之间的直接量化关系，确定函数Q，实现约束***特性初步的单独评价；本步骤同样采用多元回归的方法近似构造Ar与约束***刚度参数的函数关系Q。
由于约束***刚度参数较少，本技术方案选取所有约束***刚度参数k1、k2、GL、D2构建多元二次回归模型。为方便描述将这四个参数分别用y1、y2、y3、y4四个变量表示，对应因变量为Ar。多元二次回归模型表达式如下：
其中b0为常数项，b1、b2、b3、b4为一次项系数，b(1,1)、b(1,2)……b(4,4)为二次项系数。通过回归分析剔除影响较小的项目，其中参数D2被完全剔除，剩下对回归模型影响较大的7个项目，各个项目以及对应的系数如表6所示。
表6约束***刚度参数二次回归对应项目和系数
该多元二次回归模型的分析结果如表7所示。其中判定系数接近1，F统计量较大，检验的P值为0，说明回归结果较好。
表7约束***刚度参数二次回归模型的分析结果
Ar代表了约束***刚度的平均水平，通过对比不同约束***刚度的Ar大小可以实现约束***刚度的优劣评价。为了使Ar无量纲化并且与Av之间同趋化，本文采用min-max标准化方法，对Ar作标准化处理，得到无量纲约束***特性的评价指标，定义为“约束***综合评价指标”，用符号β表示，其计算方法为：
根据乘员响应面确定Ar最大值为69.12g，最小值为44.60g。根据公式(10)和(11)，得到β计算为：
β分布在0到1之间，仅从乘员加速度角度考虑，该值越小对应的约束***刚度越好。然而当约束***刚度参数超出本技术方案确定的常用范围时可能会出现β大于1或者小于0的情况。
3)碰撞波形与约束***特性综合评价指标的ao的建立
通过碰撞波形综合评价指标α以及约束***综合评价指标β可以分别对碰撞波形以及约束***特性作出初步的单独评价。然而乘员最终的保护效果是由碰撞波形以及约束***特性两者共同决定。即使碰撞波形设计很好，如果对应较差的约束***特性，很可能总体达不到较好的乘员保护效果，只能说良好的碰撞波形是保证约束***相容匹配的一个前提条件。同样良好的约束***特性也只是碰撞波形相容匹配的一个前提条件。
本步骤主要建立乘员加速度响应ao与碰撞波形综合评价指标α和约束***综合评价指标β的直接量化关系，从而根据碰撞波形以及约束***特性初步的单独评价结果实现整体耦合关系评价。由于α和β分别由Av和Ar标准化处理得到，所以ao与Av和Ar之间的递增函数R可等价于ao与α和β之间的递增函数关系R。可将乘员响应面按照α与β的对应关系，绘制出α-β-ao响应面，为了便于观察，用网格代替响应面，如图6所示，整个曲面呈现出较好的递增趋势。
本技术方案采用曲面拟合的方法近似构造ao与α和β之间的递增函数关系R。以α和β为自变量，对α-β-ao响应面进行二次曲面拟合，曲面拟合如公式(13)所示。其中c0、c1、c2、c3、c4、c5为各项系数。通过最小二乘法得到各项系数值如表8所示。
ao＝c0+c1α+c2β+c3α2+c4αβ+c5β2 (13)
表8α-β-ao响应面二次曲面拟合结果
由于α以及β的平方项系数接近0，所以曲面拟合公式可简化为：
ao＝280.5+335.1α+161.8β+168αβ (14)
该拟合结果的决定系数R2达到0.986，拟合精度较高。
实施例1
M6车型碰撞波形及约束***特性耦合关系评价
本实施例根据碰撞波形及约束***特性耦合关系评价方法对M6车型的碰撞波形、约束***特性以及两者耦合关系做出评价。M6车型的碰撞波形参数以及约束***刚度参数见表9。评价方法步骤如下：
表9M6车型碰撞波形及约束***刚度参数
1.建立基于单自由度模型的乘员响应面
根据表1中设定好的等效双梯形波和三线性刚度曲线的形状参数的取值范围和步长，可得到3204个不同的双梯形波形与3136个三线性约束刚度。将获得的双梯形波形与约束***刚度进行一对一组合，利用已有的单自由度模型迭代算法快速对每个组合进行乘员响应求解，得到共计约1000万个不同双梯形波形与不同约束***刚度的匹配结果，如表2，表中a0为乘员加速度峰值。
以不同双梯形波形从1到3204编号作为X坐标，以不同约束***刚度从1到3136编号作为Y坐标，以乘员加速度峰值作为Z坐标，可绘制在不同碰撞波形和不同约束***刚度组合情况下的乘员加速度峰值响应面，简称为乘员响应面
对乘员响应面向X坐标方向求平均，得到每个碰撞波形对应的所有乘员响应的平均值，用符号Av表示。对该响应面对Y坐标方向求平均，即可得到每个约束***刚度对应的所有响应的平均值，用符号Ar表示。
2.乘员响应与碰撞波形参数相关性分析
从乘员响应面中提取出碰撞波形基本参数(G1、G2、C1、Cmax、tC、tE、i、w、KAE、a、Go、to)与Av对应关系，进行线性回归分析，并将相关性系数汇总，如表3所示。由表3可知，Go、tE、KAE、G2与乘员响应的相关性较高。
3.碰撞波形与约束***特性耦合关系的评价
1)构建Av与碰撞波形参数(G1、G2、tE、GO、KAE)的函数关系P(多元二次回归模型)，分析结果如表5：
表5碰撞波形参数二次回归模型的分析结果
采用min-max标准化方法，利用公式(5)对Av作标准化处理，得到无量纲的碰撞波形综合评价指标α。根据乘员响应面确定Av最大值为79.28g，最小值为36.70g。根据公式(7)和(8)，得到α计算式为公式(9)；
根据表9中参数利用公式(9)可以计算得到M6车型的碰撞波形综合评价指标α：
2)同样构建Ar与约束***刚度参数(k1、k2、GL)的函数关系Q(多元二次回归模型)，分析结果如表6：
表6约束***刚度参数二次回归对应项目和系数
采用min-max标准化方法，根据公式(7)，对Ar作标准化处理，得到无量纲的约束***综合评价指标β。根据乘员响应面确定Ar最大值为69.12g，最小值为44.60g。根据公式(10)和(11)，得到β计算式为公式(12)：
同样根据表9中参数利用公式(12)可以计算得到M6车型的约束***综合评价指标β：
3)采用曲面拟合的方法近似构造乘员加速度响应ao与碰撞波形综合评价指标α和约束***综合评价指标β之间的递增函数关系R。以α和β为自变量，对α-β-ao响应面进行二次曲面拟合，曲面拟合公式如(13)所示，其中c0、c1、c2、c3、c4、c5为各项系数。通过最小二乘法得到各项系数值如表8所示；
表8α-β-ao响应面二次曲面拟合结果
舍掉α和β的平方项，曲面拟合公式为：
ao＝280.5+335.1α+161.8β+168αβ (14)
利用α-β-ao的曲面拟合公式(14)，可以计算得到M6车型的乘员加速度响应ao：
ao＝280.5+335.1α+161.8β+168αβ
＝280.5+335.1×0.56+161.8×0.44+168×0.56×0.44
＝580.74m/s2
已经求得的α和β的计算式，以及α-β-ao的曲面拟合公式可以直接应用于任何车型的评价；
由M6车型的碰撞波形和约束***特性的耦合关系评价可知，乘员加速度峰值为580.74m/s2超过50g，表明该车型的碰撞波形和约束***特性之间匹配的不好。β比α小，相对而言碰撞波形更差一些，后续可以优化碰撞波形来降低乘员加速度峰值。
实施例2
本实施例本文采用M6乘员约束***仿真模型，在Madymo软件中对本发明所述的三个评级指标做进一步验证。
从3204个碰撞波形中选取了α分别为0、0.2、0.4、0.6、0.8、1的六个双梯形波形，如图8所示。本文从3136个约束***刚度中选取了β分别为0、0.2、0.4、0.6、0.8、1时对应的六个约束***刚度，其约束***参数，如表9所示。
表9约束***刚度参数转化为约束***参数
将选好的六个波形依次作为M6乘员约束***仿真模型的加速度场，约束***不变，通过仿真计算得到对应的乘员胸部加速度响应曲线。根据α与乘员胸部加速度峰值关系，如图8所示，两者具有较高的一致性。
将M6乘员约束***仿真模型的约束***参数依次调整为表9中六组约束***参数值，碰撞波形不变，通过仿真计算得到对应的乘员胸部加速度响应曲线。根据β与乘员胸部加速度峰值关系，如图9所示，两者具有较高的一致性。
将6个碰撞波形和6组约束***参数交叉组合成36个耦合情况，依次带入到M6乘员约束***仿真模型，通过仿真计算得到36组乘员胸部加速度曲线。提取每个乘员胸部加速度曲线峰值，绘制成响应面，如图10所示，与拟合公式得到的响应面整体趋势基本一致，乘员胸部加速度峰值误差控制在10％。
Claims (7)
1.一种碰撞波形与约束***特性耦合关系的评价方法,其特征在于，所述的碰撞波形与约束***特性耦合关系的评价方法的步骤如下：
1)基于单自由度模型的乘员响应面的建立：
(1)定义双梯形波和三线性约束刚度曲线的形状特征参数:
(2)单自由度模型求解乘员响应面；
(3)对响应面XY向分别求平均：
对乘员响应面向X坐标方向求平均，得到每个碰撞波形对应的所有乘员响应的平均值，用符号Av表示；
对乘员响应面向Y坐标方向求平均，即可得到每个约束***刚度对应的所有乘员响应的平均值，用符号Ar表示；
2)乘员响应与碰撞波形参数相关性分析：
(1)碰撞波形参数的定义补充；
(2)从乘员响应面中提取出碰撞波形基本参数与Av对应关系，进行线性回归分析，并将相关性系数汇总；
3)碰撞波形与约束***特性耦合关系的评价：
(1)碰撞波形综合评价指标α的建立；
(2)约束***综合评价指标β的建立；
(3)碰撞波形与约束***特性综合评价指标的ao的建立。
2.按照权利要求1所述的碰撞波形与约束***特性耦合关系的评价方法,其特征在于，所述的定义双梯形波和三线性约束刚度曲线的形状特征参数是指:
G1为第一阶加速度值，即BC段所对应的加速度值；第一阶加速度的持续时间为tBC，单位为g；
G2为第二阶加速度值，即DE段对应的加速度值,第二阶加速度的持续时间为tDE，单位为g；
s1、s2、s3分别为直线AB、CD以及回弹段EF的斜率；
C1为发动机前端的压溃量，单位为m；
Cmax为车体前端结构的最大动态压溃量，单位为m；
k1为OP段斜率为安全带比刚度，OP段对应安全带线性伸长阶段；
GL为限力加速度；
k2为QR段斜率为安全气囊刚度；
D1为安全带作用力达到限力值时乘员相对于车体的运动位移；
D2为安全气囊开始作用时乘员相对于车体的运动位移。
3.按照权利要求1所述的碰撞波形与约束***特性耦合关系的评价方法,其特征在于，所述的碰撞波形参数的定义补充是指：
1)阶梯比i：定义等效双梯形波的两个阶高度G2与G1之比为阶梯比；
2)宽度比w：宽度比为发动机前端的压溃量C1与最大动态压溃量Cmax之比；
3)平均刚度KAE：平均刚度为等效双梯形波的第二阶高度G2与最大动态压溃量Cmax之比；
4)能量密度比a：能量密度比为发动机前端吸能量与车体前端结构总吸能量的比；
<mrow>
<mi>a</mi>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mi>C</mi>
</msub>
<msub>
<mi>G</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msubsup>
<mi>G</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>/</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
</mrow>
<msub>
<mi>v</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:v0为碰撞速度，单位m/s；G1为第一台阶的高度，单位为g；tC为发动机碰撞时刻，单位s；s1为直线AB段的斜率，单位g/s。
5)波形形心(to，Go)：波形形心是指加速度-时间曲线与时间轴围成的几何形状的形心，形心的横坐标称为形心时刻to，形心纵坐标称为形心加速度Go；
<mrow>
<msub>
<mi>t</mi>
<mi>o</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<mi>t</mi>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>v</mi>
</msub>
<mi>d</mi>
<mi>t</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>v</mi>
</msub>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>G</mi>
<mi>o</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<msup>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>v</mi>
</msub>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>t</mi>
</mrow>
<mrow>
<mn>2</mn>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>v</mi>
</msub>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:t为时间，单位s；为车体加速度，单位g。
4.按照权利要求1所述的碰撞波形与约束***特性耦合关系的评价方法,其特征在于，所述的单自由度模型求解乘员响应面是指：
(1)求解该响应面的单自由度模型将车辆与乘员分别简化为集中质量块M和m,车体结构在碰撞中的压溃过程简化为弹簧刚度K的压缩过程，并近似为等效双梯形波；约束***在碰撞中的变形过程简化为弹簧刚度k的压缩过程，并近似为三线性约束刚度；乘员在两个弹簧振动***的共同振动作用下做向前的减速运动，因此乘员响应的大小由两个弹簧振动***的耦合关系决定；
(2)根据设定好的形状参数范围和步长，见表1，可得到3204个不同的双梯形波形与3136个三线性约束刚度；将获得的双梯形波形与约束***刚度进行一对一组合，利用已有的单自由度模型迭代算法快速对每个组合进行乘员响应求解，得到共计约1000万个不同双梯形波形与不同约束***刚度的匹配结果，如表2，表中a0为乘员加速度峰值；
表1 形状参数的范围和步长
表2 不同碰撞波形参数与约束***刚度参数的匹配结果
(3)以不同双梯形波形从1到3204编号作为X坐标，以不同约束***刚度从1到3136编号作为Y坐标，以乘员加速度峰值作为Z坐标，可绘制在不同碰撞波形和不同约束***刚度组合情况下的乘员加速度峰值响应面，简称为乘员响应面。
5.按照权利要求1所述的碰撞波形与约束***特性耦合关系的评价方法,其特征在于，所述的碰撞波形综合评价指标α的建立是指：
本步骤主要建立波形平均响应Av与碰撞波形参数之间的直接量化关系，确定函数P，实现碰撞波形初步的单独评价，本技术方案采用多元回归的方法近似构造Av与碰撞波形参数的函数关系P；
碰撞波形参数当中Go、tE、KAE、G2与乘员响应的相关性较高，但考虑到碰撞波形基本参数的重要性，所以本技术方案选取G1、G2、tC、tE、Go、KAE这六个碰撞波形参数进行多元回归分析；考虑到波形参数的复杂性以及相互之间的关系，本技术方案构建了多元二次回归模型，为方便描述将这六个参数分别用x1、x2……x6六个变量表示，对应因变量为Av，多元二次回归模型表达式如下：
<mrow>
<msub>
<mi>A</mi>
<mi>v</mi>
</msub>
<mo>=</mo>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>6</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&ap;</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>6</mn>
</munderover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>6</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mi>m</mi>
</mrow>
<mn>6</mn>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mi>m</mi>
</msub>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中a0为常数项，a1、a2……a6为一次项系数，a(1,1),a(1,2)……a(6,6)为二次项系数。通过回归分析剔除影响较小的项目，其中tC项被完全剔除，剩下对回归结果影响较大的7个项目，各个项目以及对应的系数如表4所示。
表4 碰撞波形参数二次回归对应项目和系数
该多元二次回归模型的分析结果如表5所示，其中判定系数为接近1，F统计量较大，检验的P值为0，说明回归结果较好；
表5 碰撞波形参数二次回归模型的分析结果
Av代表了碰撞波形的平均水平，通过对比不同碰撞波形的Av大小可以实现碰撞波形的优劣评价，为了使Av无量纲化并且与Ar之间同趋化，本技术方案采用min-max标准化方法，对Av作标准化处理，得到无量纲的碰撞波形评价指标，定义为“碰撞波形综合评价指标”，用符号α表示，碰撞波形综合评价指标α的计算方法为：
<mrow>
<mi>&alpha;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>v</mi>
</msub>
<mo>-</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
根据乘员响应面确定Av最大值为79.28g，最小值为36.70g，根据公式(7)和(8)，得到α计算公式为：
<mrow>
<mi>&alpha;</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>42.58</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>76.2</mn>
<mo>-</mo>
<mn>0.7</mn>
<msub>
<mi>G</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mn>0.9</mn>
<msub>
<mi>G</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>0.01</mn>
<msub>
<mi>G</mi>
<mn>1</mn>
</msub>
<msub>
<mi>G</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>1875</mn>
<msub>
<mi>t</mi>
<mi>E</mi>
</msub>
<mo>-</mo>
<mn>14447</mn>
<msup>
<msub>
<mi>t</mi>
<mi>E</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>3.8</mn>
<msub>
<mi>G</mi>
<mi>o</mi>
</msub>
<mo>+</mo>
<mn>0.03</mn>
<msub>
<mi>K</mi>
<mrow>
<mi>A</mi>
<mi>E</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
α分布在0到1之间，从乘员加速度角度考虑，该值越小对应的碰撞波形越好。然而当碰撞波形参数超出本技术方案确定的常用范围时会出现α大于1或者小于0的情况。
6.按照权利要求1所述的碰撞波形与约束***特性耦合关系的评价方法,其特征在于，所述的约束***综合评价指标β的建立是指：
本步骤主要建立Ar与约束***刚度参数之间的直接量化关系，确定函数Q，实现约束***特性初步的单独评价；本步骤同样采用多元回归的方法近似构造Ar与约束***刚度参数的函数关系Q；
由于约束***刚度参数较少，本技术方案选取所有约束***刚度参数k1、k2、GL、D2构建多元二次回归模型，为方便描述将这四个参数分别用y1、y2、y3、y4四个变量表示，对应因变量为Ar，多元二次回归模型表达式如下：
<mrow>
<msub>
<mi>A</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>4</mn>
</munderover>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>4</mn>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mi>m</mi>
</mrow>
<mn>4</mn>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msub>
<msub>
<mi>y</mi>
<mi>m</mi>
</msub>
<msub>
<mi>y</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中b0为常数项，b1、b2、b3、b4为一次项系数，b(1,1)、b(1,2)……b(4,4)为二次项系数；通过回归分析剔除影响较小的项目，其中参数D2被完全剔除，剩下对回归模型影响较大的7个项目，各个项目以及对应的系数如表6所示；
表6 约束***刚度参数二次回归对应项目和系数
该多元二次回归模型的分析结果如表7所示，其中判定系数接近1，F统计量较大，检验的P值为0，说明回归结果较好；
表7 约束***刚度参数二次回归模型的分析结果
Ar代表了约束***刚度的平均水平，通过对比不同约束***刚度的Ar大小可以实现约束***刚度的优劣评价，为了使Ar无量纲化并且与Av之间同趋化，本文采用min-max标准化方法，对Ar作标准化处理，得到无量纲约束***特性的评价指标，定义为“约束***综合评价指标”，用符号β表示，其计算方法为：
<mrow>
<mi>&beta;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>r</mi>
</msub>
<mo>-</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
根据乘员响应面确定Ar最大值为69.12g，最小值为44.60g。根据公式(10)和(11)，得到β计算为：
<mrow>
<mi>&beta;</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2452</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>1295</mn>
<mo>+</mo>
<mn>21</mn>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>0.006</mn>
<msup>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>0.01</mn>
<msup>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</msup>
<mo>-</mo>
<mn>0.005</mn>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mn>0.3</mn>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>G</mi>
<mi>L</mi>
</msub>
<mo>-</mo>
<mn>0.2</mn>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<msub>
<mi>G</mi>
<mi>L</mi>
</msub>
<mo>+</mo>
<mn>2</mn>
<msup>
<msub>
<mi>G</mi>
<mi>L</mi>
</msub>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
β分布在0到1之间，从乘员加速度角度考虑，该值越小对应的约束***刚度越好；然而当约束***刚度参数超出本技术方案确定的常用范围时可能会出现β大于1或者小于0的情况。
7.按照权利要求1所述的碰撞波形与约束***特性耦合关系的评价方法,其特征在于，所述的碰撞波形与约束***特性综合评价指标的ao的建立是指：
本步骤主要建立乘员加速度响应ao和碰撞波形综合评价指标α与约束***综合评价指标β的直接量化关系，从而根据碰撞波形以及约束***特性初步的单独评价结果实现整体耦合关系评价，由于α和β分别由Av和Ar标准化处理得到，所以ao和Av与Ar之间的递增函数R可等价于ao和α与β之间的递增函数关系R；可将乘员响应面按照α与β的对应关系，绘制出α-β-ao响应面，为了便于观察，用网格代替响应面，整个曲面呈现出较好的递增趋势；
本技术方案采用曲面拟合的方法近似构造ao和α与β之间的递增函数关系R，以α与β为自变量，对α-β-ao响应面进行二次曲面拟合，曲面拟合如公式(13)所示，其中c0、c1、c2、c3、c4、c5为各项系数，通过最小二乘法得到各项系数值如表8所示；
ao＝c0+c1α+c2β+c3α2+c4αβ+c5β2 (13)
表8 α-β-ao响应面二次曲面拟合结果
由于α与β的平方项系数接近0，所以曲面拟合公式可简化为：
ao＝280.5+335.1α+161.8β+168αβ (14)
该拟合结果的决定系数R2达到0.986，拟合精度较高。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710649471.9A CN107577843B (zh) | 2017-08-02 | 2017-08-02 | 碰撞波形与约束***特性耦合关系的评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710649471.9A CN107577843B (zh) | 2017-08-02 | 2017-08-02 | 碰撞波形与约束***特性耦合关系的评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107577843A true CN107577843A (zh) | 2018-01-12 |
| CN107577843B CN107577843B (zh) | 2021-10-01 |
Family
ID=61034131
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710649471.9A Expired - Fee Related CN107577843B (zh) | 2017-08-02 | 2017-08-02 | 碰撞波形与约束***特性耦合关系的评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107577843B (zh) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108229079A (zh) * | 2018-03-20 | 2018-06-29 | 北京海纳川汽车部件股份有限公司 | 车辆的车体和约束***匹配的设计方法和车辆 |
| CN108647464A (zh) * | 2018-05-18 | 2018-10-12 | 吉林大学 | 概念设计阶段约束***的设计方法 |
| CN108932364A (zh) * | 2018-05-18 | 2018-12-04 | 吉林大学 | 参数化的约束***刚度设计方法 |
| CN109543259A (zh) * | 2018-11-09 | 2019-03-29 | 中国汽车技术研究中心有限公司 | 一种构建等效实车碰撞波形的方法 |
| CN109596365A (zh) * | 2018-12-27 | 2019-04-09 | 重庆长安汽车股份有限公司 | 一种模拟偏置碰撞的台车试验方法 |
| CN111580500A (zh) * | 2020-05-11 | 2020-08-25 | 吉林大学 | 一种针对自动驾驶汽车安全性的评价方法 |
| CN112347665A (zh) * | 2019-07-22 | 2021-02-09 | 广州汽车集团股份有限公司 | 构建车体结构与乘员损伤评价的关联模型的方法、装置和设备 |
| CN112382820A (zh) * | 2020-11-12 | 2021-02-19 | 上海理工大学 | 一种主动控制电池保护装置及其控制方法 |
| CN112818473A (zh) * | 2021-03-03 | 2021-05-18 | 吉林大学 | 一种用于求解汽车mpdb碰撞工况***动力学响应的解析方法 |
| CN113954991A (zh) * | 2021-12-06 | 2022-01-21 | 中国第一汽车股份有限公司 | 一种约束***的参数获取方法及装置 |
| CN117828310A (zh) * | 2024-03-04 | 2024-04-05 | 中汽研汽车检验中心（天津）有限公司 | 基于传递函数矩阵的胸部冲击损伤预测方法、设备及介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5345402A (en) * | 1992-02-25 | 1994-09-06 | Automotive Systems Laboratory, Inc. | Vehicle crash simulator system for testing crash sensors |
| CN102214256A (zh) * | 2011-05-20 | 2011-10-12 | 中国汽车技术研究中心 | 一种汽车碰撞波形特征参数提取与梯形波构建方法 |
| CN106441941A (zh) * | 2016-11-25 | 2017-02-22 | 北京汽车股份有限公司 | 汽车正面碰撞性能的评价方法和装置 |
| CN106599430A (zh) * | 2016-12-07 | 2017-04-26 | 江苏大学 | 基于能量分析的乘员约束***优化方法 |
-
2017
- 2017-08-02 CN CN201710649471.9A patent/CN107577843B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5345402A (en) * | 1992-02-25 | 1994-09-06 | Automotive Systems Laboratory, Inc. | Vehicle crash simulator system for testing crash sensors |
| CN102214256A (zh) * | 2011-05-20 | 2011-10-12 | 中国汽车技术研究中心 | 一种汽车碰撞波形特征参数提取与梯形波构建方法 |
| CN106441941A (zh) * | 2016-11-25 | 2017-02-22 | 北京汽车股份有限公司 | 汽车正面碰撞性能的评价方法和装置 |
| CN106599430A (zh) * | 2016-12-07 | 2017-04-26 | 江苏大学 | 基于能量分析的乘员约束***优化方法 |
Non-Patent Citations (1)
| Title |
|---|
| 张滕滕 等: "正面碰撞耐撞性的波形评价研究及相关性分析", 《上海汽车》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108229079A (zh) * | 2018-03-20 | 2018-06-29 | 北京海纳川汽车部件股份有限公司 | 车辆的车体和约束***匹配的设计方法和车辆 |
| CN108647464A (zh) * | 2018-05-18 | 2018-10-12 | 吉林大学 | 概念设计阶段约束***的设计方法 |
| CN108932364A (zh) * | 2018-05-18 | 2018-12-04 | 吉林大学 | 参数化的约束***刚度设计方法 |
| CN108647464B (zh) * | 2018-05-18 | 2019-08-09 | 吉林大学 | 概念设计阶段约束***的设计方法 |
| CN108932364B (zh) * | 2018-05-18 | 2019-09-03 | 吉林大学 | 参数化的约束***刚度设计方法 |
| CN109543259A (zh) * | 2018-11-09 | 2019-03-29 | 中国汽车技术研究中心有限公司 | 一种构建等效实车碰撞波形的方法 |
| CN109543259B (zh) * | 2018-11-09 | 2023-03-31 | 中国汽车技术研究中心有限公司 | 一种构建等效实车碰撞波形的方法 |
| CN109596365A (zh) * | 2018-12-27 | 2019-04-09 | 重庆长安汽车股份有限公司 | 一种模拟偏置碰撞的台车试验方法 |
| CN112347665A (zh) * | 2019-07-22 | 2021-02-09 | 广州汽车集团股份有限公司 | 构建车体结构与乘员损伤评价的关联模型的方法、装置和设备 |
| CN112347665B (zh) * | 2019-07-22 | 2024-04-23 | 广州汽车集团股份有限公司 | 构建车体结构与乘员损伤评价的关联模型的方法、装置和设备 |
| CN111580500B (zh) * | 2020-05-11 | 2022-04-12 | 吉林大学 | 一种针对自动驾驶汽车安全性的评价方法 |
| CN111580500A (zh) * | 2020-05-11 | 2020-08-25 | 吉林大学 | 一种针对自动驾驶汽车安全性的评价方法 |
| CN112382820A (zh) * | 2020-11-12 | 2021-02-19 | 上海理工大学 | 一种主动控制电池保护装置及其控制方法 |
| CN112818473A (zh) * | 2021-03-03 | 2021-05-18 | 吉林大学 | 一种用于求解汽车mpdb碰撞工况***动力学响应的解析方法 |
| CN113954991A (zh) * | 2021-12-06 | 2022-01-21 | 中国第一汽车股份有限公司 | 一种约束***的参数获取方法及装置 |
| CN117828310A (zh) * | 2024-03-04 | 2024-04-05 | 中汽研汽车检验中心（天津）有限公司 | 基于传递函数矩阵的胸部冲击损伤预测方法、设备及介质 |
| CN117828310B (zh) * | 2024-03-04 | 2024-05-31 | 中汽研汽车检验中心（天津）有限公司 | 基于传递函数矩阵的胸部冲击损伤预测方法、设备及介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107577843B (zh) | 2021-10-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107577843A (zh) | 碰撞波形与约束***特性耦合关系的评价方法 | |
| CN106599430B (zh) | 基于能量分析的乘员约束***优化方法 | |
| Filipi et al. | Engine-in-the-loop testing for evaluating hybrid propulsion concepts and transient emissions–HMMWV case study | |
| Oh | Evaluation of motor characteristics for hybrid electric vehicles using the hardware-in-the-loop concept | |
| CN107256289A (zh) | 汽车碰撞简化参数化有限元模型的建立方法 | |
| Fajri et al. | Emulating on-road operating conditions for electric-drive propulsion systems | |
| CN102339350B (zh) | 一种基于参数化设计的整车碰撞仿真方法 | |
| CN107063718A (zh) | 汽车正面碰撞波形参数化评价方法 | |
| CN112948983B (zh) | 一种协同正面碰撞多工况的汽车前端结构能量管理方法 | |
| Vairavel et al. | Modeling, validation and simulation of electric vehicles using MATLAB | |
| Zhang et al. | Identification of key design parameters of high-speed train for optimal design | |
| Hou et al. | Modeling and simulation of hybrid electric vehicles using HEVSIM and ADVISOR | |
| CN107256656A (zh) | 一种随动视景汽车驾驶模拟器滞后综合校正方法 | |
| Van Mierlo et al. | Vehicle simulation program: a tool to evaluate hybrid power management strategies based on an innovative iteration algorithm | |
| CN103279598B (zh) | 一种汽车车身多变量抗撞性优化设计的变量筛选方法 | |
| CN103398833B (zh) | 一种反映驾驶员肌肉动态特性的碰撞假人模型设计方法 | |
| Xie et al. | Multi-objective optimisation of a vehicle energy absorption structure based on surrogate model | |
| Zhang et al. | A two degrees of freedom model–based optimization method for occupant restraint systems in vehicle crash | |
| Harmati et al. | Energy distribution modeling of car body deformation using LPV representations and fuzzy reasoning | |
| Czapnik et al. | Conceptual design of battery electric vehicle powertrains | |
| Ahmad et al. | Applications of Hardware-in-the-Loop Simulation in Automotive Embedded Systems | |
| Jie et al. | Research on frontal collision reconstruction model based on coupling of PCCRASH and MADYMO | |
| CN109635484A (zh) | 基于多学科优化的混合动力车辆动力***优化测试方法 | |
| CN112818473B (zh) | 一种用于求解汽车mpdb碰撞工况***动力学响应的解析方法 | |
| Harmati et al. | Identification of energy distribution for crash deformational processes of road vehicles |
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 |
Granted publication date: 20211001 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |