CN106599430B - 基于能量分析的乘员约束系统优化方法 - Google Patents
基于能量分析的乘员约束系统优化方法 Download PDFInfo
- Publication number
- CN106599430B CN106599430B CN201611113305.9A CN201611113305A CN106599430B CN 106599430 B CN106599430 B CN 106599430B CN 201611113305 A CN201611113305 A CN 201611113305A CN 106599430 B CN106599430 B CN 106599430B
- Authority
- CN
- China
- Prior art keywords
- value
- energy
- occupant
- energy absorption
- relationship
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000005457 optimization Methods 0.000 title claims abstract description 26
- 238000004458 analytical method Methods 0.000 title claims abstract description 21
- 238000010521 absorption reaction Methods 0.000 claims abstract description 87
- 238000012360 testing method Methods 0.000 claims abstract description 36
- 230000006378 damage Effects 0.000 claims abstract description 33
- 208000027418 Wounds and injury Diseases 0.000 claims abstract description 8
- 208000014674 injury Diseases 0.000 claims abstract description 8
- 230000001133 acceleration Effects 0.000 claims description 24
- 238000004088 simulation Methods 0.000 claims description 18
- 210000000689 upper leg Anatomy 0.000 claims description 13
- 210000003127 knee Anatomy 0.000 claims description 10
- 210000002414 leg Anatomy 0.000 claims description 10
- 206010019196 Head injury Diseases 0.000 claims description 8
- 238000000418 atomic force spectrum Methods 0.000 claims description 7
- 239000004744 fabric Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 238000011156 evaluation Methods 0.000 claims description 5
- 244000309466 calf Species 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 4
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 4
- 230000004888 barrier function Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000010206 sensitivity analysis Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 210000001364 upper extremity Anatomy 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 2
- 230000000875 corresponding effect Effects 0.000 claims 2
- 239000002131 composite material Substances 0.000 claims 1
- 230000002596 correlated effect Effects 0.000 claims 1
- 238000005516 engineering process Methods 0.000 claims 1
- 239000000284 extract Substances 0.000 claims 1
- 230000008676 import Effects 0.000 claims 1
- 230000005540 biological transmission Effects 0.000 abstract 1
- 238000004904 shortening Methods 0.000 abstract 1
- 210000000038 chest Anatomy 0.000 description 14
- 230000000694 effects Effects 0.000 description 4
- 210000001699 lower leg Anatomy 0.000 description 2
- 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 description 1
- 206010039203 Road traffic accident Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 238000009863 impact test Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 210000000115 thoracic cavity Anatomy 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
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)
- Air Bags (AREA)
Abstract
本发明公开了基于能量分析的乘员约束系统优化方法，该方法首先建立某A级车的正面碰撞模型，通过仿真和试验的方法验证其有效性。然后结合能量理论，从乘员的能量传递方式出发，分析各个约束子系统的吸能情况，获得约束系统能量分配曲线。结合正交试验的多组数据，比较优化前后约束系统能量的分配差异，寻找到约束系统能量分配与乘员综合损伤值(WIC)的关系。再分析各个约束子系统的参数，构建参数与约束系统能量分配的关系，从而可调节相关参数，实现约束系统对乘员保护的性能优化。本发明对于提高乘员约束系统保护性能，缩短优化周期具有重要意义。
Description
技术领域
本发明属于汽车安全领域，尤其涉及一种基于能量分析的乘员约束系统优化方法。
背景技术
车辆正面碰撞是交通事故中最常见的碰撞形式。在正碰过程中，乘员约束系统对乘员的保护作用十分显著，对其性能的研究非常有必要。目前，主要以乘员损伤值作为评价依据，对约束系统进行优化。然而，乘员自身存在很多不确定因素，以此来间接地评价约束系统，步骤繁琐而且并不精确。因此可以从能量的角度，分析各个约束子系统的吸能情况，寻找出约束系统能量分配与其对乘员保护之间的关系。根据约束系统参数与其吸能的关系，优化乘员约束系统的性能，使其达到更好的保护效果。
发明内容
针对乘员约束系统常规的优化方法中以乘员损伤值作为评价依据，步骤繁琐并存在不确定性因素，本发明提出一种基于能量分析的乘员约束系统优化方法，旨在更加便捷与精准地优化乘员约束系统，提高其保护性能。
根据本发明的一个方面，提出一种基于能量分析的乘员约束系统优化方法，所述方法包括以下步骤：
步骤1：根据中国汽车技术研究中心2015年版的C-NCAP管理规则，进行车辆正面100％重叠刚性壁障碰撞试验。并在试验过程中测量车辆、约束系统以及乘员的相关曲线。
步骤2：根据车辆正面碰撞试验数据，建立正碰仿真模型。
步骤3：车辆正面碰撞仿真模型的验证。判断与实车碰撞试验是否相符。如果相符则转入步骤4；如果不相符，则转入步骤2。直至符合试验情形，即模型有效。
步骤4：采用常规方式优化车辆正碰模型，在优化过程中获得乘员损伤值，作为第一组数据。
步骤5：按照C-NCAP法规对乘员损伤值的评价以及乘员各部位间的连接关系，将乘员划分为头部、胸部、髋部和腿部四个部分。
步骤6：将汽车正碰仿真模型划分为5个子模型：Hybrid_III_50th假人多刚体模型，汽车车身多刚体模型，安全气囊有限元模型，安全肩带和安全腰带混合模型，汽车座椅有限元模型。根据能量分析理论，在正碰过程中，子模型之间主要产生四种形式的力，即惯性力、阻尼和摩擦力、弹性变形力、外界载荷作用力。由运动学方程的力学平衡可知能量也平衡。因此对应四种形式的能量，即动能、耗散能、内能、外能，依据如下所示的能量公式，计算并输出乘员各部位的能量曲线。
式中，WTotal为总能量；v为碰撞速度，单位为m/s；M为模型质量,单位为kg；a为碰撞加速度，单位为m/s2；C为阻尼系数；K为弹性模量；l为变形量，单位为m；Fext为外界载荷作用力，单位为N。
步骤7：去除乘员各部位间内部的能量传递，获得各约束子系统的吸能曲线，并得到各约束子系统的吸能峰值，作为第二组数据。
步骤8：分析第一组数据和第二组数据，利用回归模型建立各个约束子系统吸能峰值与乘员损伤值之间的关系，得到关系A。
步骤9：选取与乘员损伤值存在较强相关关系(在相关系数R2大于0.6时，两者间存在较强的相关关系，见图6)的约束子系统，根据实际情况(实际情况是现实情形下每个参数都有一定的取值范围，比如安全带预紧器点火时间为8ms到28ms)，选择系统的主要参数及其取值范围，作为第三组数据。
步骤10：分别计算各个约束子系统在不同参数下的吸能峰值，作为第四组数据。
步骤11：分析第三组数据与第四组数据，建立约束系统参数与其吸能的关系，得出关系B。
步骤12：将关系A与关系B建立联系，利用此方法即可对乘员约束系统进行优化。具体地：从关系A中可知，安全肩带、安全腰带和汽车座椅吸能越多，防火墙及地板吸能越少时，乘员损伤值WIC越小，即约束系统对乘员的保护作用越好。所以在关系B得到的公式中，分别在参数范围内取对应的最大最小值，此时的参数组合，即为约束系统性能最优的组合。
优选地，在所述步骤1中，测量相关曲线包括车辆的碰撞减速度曲线，乘员头部加速度曲线、胸部加速度曲线、髋部加速度曲线、腿部径向力曲线，以及安全肩带力与安全腰带力对应曲线。
优选地，在所述步骤2中，建立车辆正碰仿真模型。包括构建车身模型，假人的导入与定位，建立乘员约束系统模型，定义各部件的载荷特性与接触，以及载入碰撞波形。
优选地，在所述步骤2中，构建车身模型，包括挡风玻璃、仪表板、转向系统、汽车座椅、油门踏板的建模。
优选地，在所述步骤2中，选用HybridШ第50百分位男性假人，采用预模拟和试验数据相结合的方式进行假人定位。
优选地，在所述步骤2中，建立乘员约束系统模型中安全肩带与腰带采用有限元和多刚体结合的方式建模。安全气囊首先三维建模，进行网格划分后，得到有限元模型。
优选地，在所述步骤2中，定义与实车试验数据一致的各部件载荷和接触特性。
优选地，在所述步骤2中，将实车试验中，车辆B柱下方测得的加速度曲线及重力场，赋予仿真模型。
优选地，在所述步骤3中，正面碰撞仿真模型的验证遵循从下至上原则，获取软件输出的大腿力曲线，髋部、胸部和头部的加速度曲线，将其与试验采集曲线进行对比。主要验证曲线的起始时刻，峰值大小，峰值时刻，波形形状是否基本相符。
优选地，在所述步骤4中，通过灵敏度分析的方式筛选出主要约束系统参数，作为设计变量。用头部损伤值HIC36ms、胸部3毫秒加速度C3ms、胸部压缩量D、大腿的径向力FF以及小腿的径向力Fc作为约束变量，以乘员综合损伤指标WIC值作为优化目标。通过正交试验的方法获得各组试验中的乘员损伤值。其中，头部损伤值HIC36ms的计算公式如下所示：
式中a(t)是在T0≤t≤TE时间段内，头部质心处合成加速度；T0是碰撞开始时刻；TE是碰撞结束时刻；法规中规定将计算前排乘员HIC的时间间隔(t2-t1)定为36ms。
乘员综合损伤指标WIC值的计算公式如下所示：
式中，0.6、0.35和0.05分别代表损伤指标加权系数；HIC36ms为头部伤害评价指标；C3ms为胸部3MS准则的值，单位为g；D为胸部压缩量，单位为m；Fl为左大腿骨轴向力，单位为kN；Fr为右大腿骨轴向力，单位为kN。
优选地，在所述步骤5中，乘员部位划分时，头部包括头部和颈部，胸部包括胸部和上肢，腿部包括大腿、小腿和脚。
优选地，所述步骤6的实现包括：用命令语句在MADYMO中定义多刚体、有限元及相互间的接触，并定义节点力和能量输出。经过MADYMO运算，利用HyperGraph，选取有限元能量和多刚体能量，计算和输出相应曲线。
优选地，所述步骤7的实现：根据乘员各部位的连接关系以及试验录像和仿真动画，分析乘员各部位的能量耗散方式。提取出各个约束子系统吸收乘员的能量曲线。
优选地，在所述步骤8中，采用直线拟合与二次多项式拟合回归模型构建约束系统吸能与乘员损伤值的关系。
优选地，在所述步骤8中，所得具体关系式如下所示：
安全肩带吸能峰值与WIC值关系：Y1＝-1322X+2516。
式中，X为WIC值；Y1为安全肩带吸能峰值，单位为J。
安全腰带吸能峰值与WIC值关系：Y2＝-705.5X+2256。
式中，X为WIC值；Y2为安全腰带吸能峰值，单位为J。
汽车座椅吸能峰值与WIC值关系：Y3＝1721X2-1697X+766.5。
式中，X为WIC值；Y3为汽车座椅吸能峰值，单位为J。
防火墙及地板吸能峰值与WIC值关系：Y4＝388.5X+516。
式中，X为WIC值；Y4为防火墙及地板吸能峰值，单位为J。
优选地，在所述步骤8中，所得关系A即为WIC值与约束系统能量峰值相关性较强、与安全肩带、安全腰带以及汽车座椅吸能峰值成负相关、与防火墙及地板的吸能峰值成正相关。当约束系统能量分配给安全肩带、安全腰带和汽车座椅的能量越多，分配给防火墙的能量越少时，WIC值越低，约束系统对乘员的保护性能越好。
优选地，在所述步骤11中，所得具体关系式如下所示：
安全肩带参数与其吸能峰值的关系式：
Y1＝1967.9-0.0078X1+0.1365X2-466.9286X3-0.0045X2 2+1293.9X3 2
式中，Y1为安全肩带吸能峰值；X1为限力器限力值，单位为N；X2为预紧器点火时间，单位为ms；X3为织物延伸率。
安全腰带参数与其吸能峰值的关系式：
Y2＝1881.1-0.0032X1+0.2371X2+1011.5X3-0.0067X2 2-2060.7X3 2
式中，Y2为安全腰带吸能峰值；X1为限力器限力值，单位为N；X2为预紧器点火时间，单位为ms；X3为织物延伸率。
汽车座椅参数与其吸能峰值的关系式：
Y3＝342.879+5.2778X1+3.8992X2-0.0976X1 2-0.3029X1X2-0.5556X2 2
式中，Y3为汽车座椅吸能峰值；X1为座椅倾角，单位为°；X2为座椅刚度，单位为N/m。
汽车防火墙及地板参数与其吸能峰值的关系式：
Y4＝683.3349-0.4925X1+50.0686X2-34.2943X3+0.0143X1 2+7.7143X2 2+18.3571X3 2
式中，Y4为汽车防火墙及地板吸能峰值；X1为膝垫倾角，单位为°；X2为地板摩擦系数；X3为膝垫刚度，单位为N/m。
优选地，在所述步骤11中，所得关系B即为根据上述四个关系式，调节各个约束子系统参数，实现安全肩带、安全腰带以及汽车座椅的吸能峰值达最大值，地板与膝垫的吸能峰值达最小值，此时约束系统吸能达最优分布。
本发明创新性地提出了一种根据正面碰撞中约束系统的吸能情况来实现对其性能的优化。本发明的有益效果：
发明所提出的基于能量分析的乘员约束系统优化方法，从约束系统自身出发，分析其在正面碰撞过程中对乘员能量的吸收，调节约束系统参数，达到最优的保护性能。该方法可更加便捷地实现对乘员约束系统的优化，进而提高其性能。
附图说明
图1是本发明所述方法的流程图。
图2是本发明中车辆正面碰撞的仿真模型图。
图3是本发明中乘员部位划分示意图。
图4是本发明中乘员各部位的能量耗散图。
其中：(a)表示头部能量分布曲线；(b)表示胸部能量分布曲线；(c)表示髋部能量分布曲线；(d)表示腿部能量分布曲线。
图5是本发明中三种典型状态下各个约束子系统的吸能曲线分析图。
图6是本发明中乘员损伤值与各个约束子系统吸能峰值的关系图。
其中：(a)表示WIC值与安全肩带吸能峰值关系；(b)表示WIC值与安全腰带吸能峰值关系；(c)表示WIC值与汽车座椅吸能峰值关系；(d)表示WIC值与防火墙及地板吸能峰值关系。
具体实施方式
下面结合附图，对本发明的具体实施方式进行详细说明。
本发明针对现有的乘员约束系统优化方法，创新性地提出从能量分析的角度，选择合适的约束系统参数对其优化，提高其保护性能。
图1表示出了本发明具体实施方式中基于能量分析的乘员约束系统优化方法的流程步骤，具体实施步骤如下：
步骤1：根据中国汽车技术研究中心2015年版的C-NCAP管理规则，进行车辆正面100％重叠刚性壁障碰撞试验。并在试验过程中测量车辆、约束系统以及乘员的相关曲线。
在本发明具体实施例中，上述相关曲线包括车辆的碰撞减速度曲线，乘员头部加速度曲线、胸部加速度曲线、髋部加速度曲线、腿部径向力曲线，以及安全肩带力与安全腰带力对应曲线。
步骤2：根据车辆正面碰撞试验的相关数据，利用MADYMO软件建立正碰仿真模型，如图2所示。
在本发明具体实施例中，建立车辆正面碰撞模型包括构建车身模型，假人的导入与定位，建立乘员约束系统模型，定义各部件的载荷特性与接触，以及载入碰撞波形。
其中构建车身模型，包括挡风玻璃、仪表板、转向系统、汽车座椅、油门踏板的建模。
上述假人的导入与定位，即调用MADYMO模型库中的HybridШ第50百分位男性假人，采用预模拟和试验数据相结合的方式进行假人定位。
建立乘员约束系统模型时，安全肩带与腰带采用有限元和多刚体结合的方式建模。安全气囊建立CAD模型，利用HYPERMESH进行网格划分，最后在MADYMO的FOLDER模块下自动折叠。
其中，定义各部件的载荷特性与接触时，定义与实车试验数据一致的各部件载荷和接触特性。
载入碰撞波形时，将实车试验中，将车辆B柱下方测得的加速度曲线及重力场赋予仿真模型。
步骤3：车辆正面碰撞仿真模型的验证。判断与实车试验是否相符。如果相符则转入步骤4；如果不相符，则转入步骤2。直至符合试验情形，即模型有效。
上述正面碰撞仿真模型的验证遵循从下至上原则，获取软件输出的大腿力曲线，髋部、胸部和头部的加速度曲线，将其与试验采集曲线进行对比。主要验证曲线的起始时刻，峰值大小，峰值时刻，波形形状是否基本相符，两者间误差应控制在5％以内。
步骤4：采用常规方式优化车辆正碰模型，在优化过程中获得乘员损伤值，作为第一组数据。
在本发明具体实施例中，用常规方式优化车辆正碰模型即通过灵敏度分析的方式筛选出主要约束系统参数，作为设计变量。用头部损伤值HIC36ms、胸部3毫秒加速度C3ms、胸部压缩量D、大腿的径向力FF以及小腿的径向力Fc作为约束变量，以乘员综合损伤指标WIC值作为优化目标。通过正交试验的方法获得各组试验中的乘员损伤值。其中，头部损伤值HIC36ms的计算公式如下所示：
式(1)中a(t)是在T0≤t≤TE时间段内，头部质心处合成加速度；T0是碰撞开始时刻；TE是碰撞结束时刻；法规中规定将计算前排乘员HIC的时间间隔(t2-t1)定为36ms。
乘员综合损伤指标WIC值的计算公式如下所示：
式(2)中，0.6、0.35和0.05分别代表损伤指标加权系数；HIC36ms为头部伤害评价指标；C3ms为胸部3MS准则的值，单位为g；D为胸部压缩量，单位为m；Fl为左大腿骨轴向力，单位为kN；Fr为右大腿骨轴向力，单位为kN。
步骤5：按照C-NCAP法规对乘员损伤值的评价以及乘员各部位间的连接关系，将乘员划分为头部、胸部、髋部和腿部四个部分，如图3所示。
上述乘员部位划分时，是进行广义的划分，即头部包括头部和颈部，胸部包括胸部和上肢，腿部包括大腿、小腿和脚。
步骤6：根据能量分析理论，以MADYMO和HyperGraph软件辅助，计算并输出乘员各部位的能量曲线，如图4所示。
在本发明具体实施例中，利用能量分析公式进行能量计算，公式如下所示：
式(3)中，WTotal为总能量；v为碰撞速度，单位为m/s；M为模型质量,单位为kg；a为碰撞加速度，单位为m/s2；C为阻尼系数；K为弹性模量；l为变形量，单位为m；Fext为外界载荷作用力，单位为N。
在本发明具体实施例中，用命令语句在MADYMO中定义多刚体、有限元及相互间的接触，并定义节点力和能量输出。经过MADYMO运算，利用HyperGraph，选取有限元能量和多刚体能量，计算和输出相应曲线。
步骤7：去除乘员各部位间内部的能量传递，获得各约束子系统的吸能曲线，并得到各约束子系统的吸能峰值，作为第二组数据。
根据乘员各部位的连接关系以及试验录像和仿真动画，分析乘员各部位的能量耗散方式。提取出各个约束子系统吸收乘员的能量曲线，如图5所示。
步骤8：分析第一组数据和第二组数据，利用MATLAB拟合回归模型建立各个约束子系统吸能峰值与乘员损伤值之间的关系，如图6所示。
在本发明具体实施例中，主要采用直线拟合与二次多项式拟合回归模型构建约束系统吸能与乘员损伤值的关系，具体关系式如下所示：
安全肩带吸能峰值与WIC值关系：
Y1＝-1322X+2516 (4)
式(4)中，X为WIC值；Y1为安全肩带吸能峰值，单位为J。
安全腰带吸能峰值与WIC值关系：
Y2＝-705.5X+2256 (5)
式(5)中，X为WIC值；Y2为安全腰带吸能峰值，单位为J。
汽车座椅吸能峰值与WIC值关系：
Y3＝1721X2-1697X+766.5 (6)
式(6)中，X为WIC值；Y3为汽车座椅吸能峰值，单位为J。
防火墙及地板吸能峰值与WIC值关系：
Y4＝388.5X+516 (7)
式(7)中，X为WIC值；Y4为防火墙及地板吸能峰值，单位为J。
在本发明具体实施例中，所得关系A即为WIC值与约束系统能量峰值相关性较强，与安全肩带、安全腰带以及汽车座椅吸能峰值成负相关，与防火墙及地板的吸能峰值成正相关。当约束系统能量分配给安全肩带、安全腰带和汽车座椅的能量越多，分配给防火墙的能量越少时，WIC值越低，约束系统对乘员的保护性能越好。
步骤9：选取与乘员损伤值存在较强相关关系的约束子系统，根据实际情况，选择系统的主要参数及其取值范围，作为第三组数据。
在本发明具体实施例中，与乘员损伤值存在较强相关关系的约束子系统为安全肩带、安全腰带、座椅、地板与膝垫。
其中，影响安全肩带与安全腰带性能的主要参数是安全带限力器限力值、安全带预紧器点火时间和安全带织物刚度；影响座椅性能的主要参数是座椅倾角和座椅刚度；影响地板的主要参数是摩擦系数；影响膝垫性能的主要参数是膝垫刚度和安装角度。
步骤10：分别计算各个子约束系统在不同参数下的吸能峰值，作为第四组数据。
步骤11：分析第三组数据与第四组数据，建立约束系统参数与其吸能的关系，得出关系B。
在本发明具体实施例中，所得具体关系式如下所示：
安全肩带参数与其吸能峰值的关系式：
Y1＝1967.9-0.0078X1+0.1365X2-466.9286X3-0.0045X2 2+1293.9X3 2 (8)
式(8)中，Y1为安全肩带吸能峰值；X1为限力器限力值，单位为N；X2为预紧器点火时间，单位为ms；X3为织物延伸率。
安全腰带参数与其吸能峰值的关系式：
Y2＝1881.1-0.0032X1+0.2371X2+1011.5X3-0.0067X2 2-2060.7X3 2 (9)
式(9)中，Y2为安全腰带吸能峰值；X1为限力器限力值，单位为N；X2为预紧器点火时间，单位为ms；X3为织物延伸率。
汽车座椅参数与其吸能峰值的关系式：
Y3＝342.879+5.2778X1+3.8992X2-0.0976X1 2-0.3029X1X2-0.5556X2 2 (10)
式(10)中，Y3为汽车座椅吸能峰值；X1为座椅倾角，单位为°；X2为座椅刚度，单位为N/m。
汽车防火墙及地板参数与其吸能峰值的关系式：
Y4＝683.3349-0.4925X1+50.0686X2-34.2943X3+0.0143X1 2+7.7143X2 2+18.3571X3 2 (11)
式(11)中，Y4为汽车防火墙及地板吸能峰值；X1为膝垫倾角，单位为°；X2为地板摩擦系数；X3为膝垫刚度，单位为N/m。
在本发明具体实施例中，所得关系B即为根据上述四个关系式(8)～(11)，调节各个约束子系统参数，实现安全肩带、安全腰带以及汽车座椅的吸能峰值达最大值，地板与膝垫的吸能峰值达最小值，此时约束系统吸能达最优分布。
步骤12：将关系A与关系B建立联系，利用此方法即可对乘员约束系统进行优化。具体地：从关系A中可知，安全肩带、安全腰带和汽车座椅吸能越多，防火墙及地板吸能越少时，乘员损伤值WIC越小，即约束系统对乘员的保护作用越好。所以在关系B得到的公式中，分别在参数范围内取对应的最大最小值，此时的参数组合，即为约束系统性能最优的组合。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明，它们并非用以限制本发明的保护范围，凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
Claims (6)
1.基于能量分析的乘员约束系统优化方法，其特征在于，包括以下步骤：
步骤1：根据中国汽车技术研究中心2015年版的C-NCAP管理规则，进行车辆正面100％重叠刚性壁障碰撞试验，并在试验过程中测量车辆、约束系统以及乘员的相关曲线；
步骤2：根据车辆正面碰撞试验数据，建立车辆正面碰撞模型；
步骤3：验证车辆正面碰撞模型，判断与实车碰撞试验是否相符，如果相符则转入步骤4；如果不相符，则转入步骤2，直至符合试验情形，即模型有效；
步骤4：采用常规方式优化车辆正面碰撞模型，在优化过程中获得乘员损伤值，作为第一组数据；
步骤5：按照C-NCAP法规对乘员损伤值的评价以及乘员各部位间的连接关系，将乘员划分为头部、胸部、髋部和腿部四个部分；
步骤6：根据能量分析理论，计算并输出乘员各部位的能量曲线；
步骤7：去除乘员各部位间内部的能量传递，获得各约束子系统的吸能曲线，并得到各约束子系统的吸能峰值，作为第二组数据；
步骤8：分析第一组数据和第二组数据，利用拟合回归模型建立各个约束子系统吸能峰值与乘员损伤值之间的关系，得出关系A；
步骤9：选取与乘员损伤值存在较强相关关系的约束子系统，根据实际情况，选择各个约束子系统的参数及其取值范围，作为第三组数据；
步骤10：分别计算各个子约束系统在不同参数下的吸能峰值，作为第四组数据；
步骤11：分析第三组数据与第四组数据，建立约束系统参数与其吸能的关系，得出关系B；
步骤12：将关系A与关系B建立关系，利用此关系即可对乘员约束系统进行优化；
所述步骤4中的具体实现：通过灵敏度分析的方式筛选出约束系统参数，作为设计变量；用头部损伤值HIC36ms、胸部3毫秒加速度C3ms、胸部压缩量D、大腿的径向力FF以及小腿的径向力Fc作为约束变量，以乘员综合损伤指标WIC值作为优化目标；
通过正交试验的方法获得各组试验中的乘员损伤值；其中，头部损伤值HIC36ms的计算公式如下所示：
式中：a(t)是在T0≤t≤TE时间段内，头部质心处合成加速度；T0是碰撞开始时刻；TE是碰撞结束时刻；规定将计算前排乘员HIC的时间间隔设为36ms；
乘员综合损伤指标WIC值的计算公式如下所示：
式中，0.6、0.35和0.05分别代表损伤指标加权系数；HIC36ms为头部损伤值；C3ms为胸部3MS准则的值，单位为g；D为胸部压缩量，单位为m；Fl为左大腿骨轴向力，单位为kN；Fr为右大腿骨轴向力，单位为kN；
所述步骤6中利用能量分析的计算公式如下所示：
式中WTotal为总能量；v为碰撞速度，单位为m/s；M为模型质量，单位为kg；a为碰撞加速度，单位为m/s2；C为阻尼系数；K为弹性模量；l为变形量，单位为m；Fext为外界载荷作用力，单位为N；
用命令语句在MADYMO中定义多刚体、有限元及相互间的接触，并定义节点力和能量输出；经过MADYMO运算，利用HyperGraph，选取有限元能量和多刚体能量，计算和输出相应曲线；
在所述步骤8中，采用直线拟合与二次多项式拟合回归模型构建约束系统吸能与乘员损伤值的关系；所得关系具体公式如下：
WIC值与安全肩带吸能峰值关系：Y1＝-1322X+2516；
WIC值与安全腰带吸能峰值关系：Y2＝-705.5X+2256；
WIC值与汽车座椅吸能峰值关系：Y3＝1721X2-1697X+766.5；
WIC值与汽车防火墙及地板吸能峰值关系：Y4＝388.5X+516；
所述关系A即为WIC值与约束系统能量峰值相关性较强，与安全肩带、安全腰带以及汽车座椅吸能峰值成负相关，与防火墙及地板的吸能峰值成正相关；
所述步骤11中建立约束系统参数与其吸能的关系具体如下：
安全肩带参数与其吸能峰值的关系式为：
Y1＝1967.9-0.0078X1+0.1365X2-466.9286X3-0.0045X2 2+1293.9X3 2
式中，Y1为安全肩带吸能峰值；X1为限力器限力值，单位为N；X2为预紧器点火时间，单位为ms；X3为织物延伸率；
安全腰带参数与其吸能峰值的关系式为：
Y2＝1881.1-0.0032X1+0.2371X2+1011.5X3-0.0067X2 2-2060.7X3 2
式中，Y2为安全腰带吸能峰值；X1为限力器限力值，单位为N；X2为预紧器点火时间，单位为ms；X3为织物延伸率；
汽车座椅参数与其吸能峰值的关系式为：
Y3＝342.879+5.2778X1+3.8992X2-0.0976X1 2-0.3029X1X2-0.5556X2 2
式中，Y3为汽车座椅吸能峰值；X1为座椅倾角，单位为°；X2为座椅刚度，单位为N/m；
汽车防火墙及地板参数与其吸能峰值的关系式如下所示：
Y4＝683.3349-0.4925X1+50.0686X2-34.2943X3+0.0143X1 2+7.7143X2 2+18.3571X3 2
式中，Y4为汽车防火墙及地板吸能峰值；X1为膝垫倾角，单位为°；X2为地板摩擦系数；X3为膝垫刚度，单位为N/m；
所述关系B即为根据上述四个关系式，调节各个约束子系统参数，实现安全肩带、安全腰带以及汽车座椅的吸能峰值达最大值，地板与膝垫的吸能峰值达最小值，此时约束系统吸能达最优分布。
2.根据权利要求1所述的基于能量分析的乘员约束系统优化方法，其特征在于，步骤1中所述相关曲线包括：测量车辆的碰撞减速度曲线、乘员头部加速度曲线、胸部加速度曲线、髋部加速度曲线、腿部径向力曲线、以及安全肩带力与安全腰带力对应曲线。
3.根据权利要求1所述的基于能量分析的乘员约束系统优化方法，其特征在于，所述步骤2中建立车辆正面碰撞模型包括：构建车身模型，假人的导入与定位，建立乘员约束系统模型，定义各部件的载荷特性与接触，以及载入碰撞波形；
所述构建车身模型包括挡风玻璃、仪表板、转向系统、汽车座椅、油门踏板的建模；
所述假人的导入与定位：调用MADYMO模型库中的Hybrid Ш第50百分位男性假人，采用预模拟和试验数据相结合的方式进行假人定位；
所述建立乘员约束系统模型时，安全肩带与腰带采用有限元和多刚体结合的方式建模；安全气囊建立CAD模型，利用HYPERMESH进行网格划分，最后在MADYMO的FOLDER模块下自动折叠；
所述定义各部件的载荷特性与接触是定义与实车试验数据一致的各部件载荷和接触特性；
所述载入碰撞波形是将实车碰撞试验中，车辆B柱下方测得的加速度曲线及重力场，赋予仿真模型。
4.根据权利要求1所述的基于能量分析的乘员约束系统优化方法，其特征在于，所述步骤3中正面碰撞模型的验证遵循从下至上原则，获取软件输出的大腿力曲线，髋部、胸部和头部的加速度曲线，将其与试验采集曲线进行对比，验证曲线的起始时刻，峰值大小，峰值时刻，波形形状是否相符。
5.根据权利要求1所述的基于能量分析的乘员约束系统优化方法，其特征在于，在所述步骤5中乘员部位划分时，其中头部包括头部和颈部，胸部包括胸部和上肢，腿部包括大腿、小腿和脚。
6.根据权利要求1所述的基于能量分析的乘员约束系统优化方法，其特征在于，所述步骤7的具体实现：根据乘员各部位的连接关系以及试验录像和仿真动画，分析乘员各部位的能量耗散方式，提取出各个约束子系统吸收乘员的能量曲线。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611113305.9A CN106599430B (zh) | 2016-12-07 | 2016-12-07 | 基于能量分析的乘员约束系统优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611113305.9A CN106599430B (zh) | 2016-12-07 | 2016-12-07 | 基于能量分析的乘员约束系统优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106599430A CN106599430A (zh) | 2017-04-26 |
| CN106599430B true CN106599430B (zh) | 2020-03-31 |
Family
ID=58597213
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611113305.9A Expired - Fee Related CN106599430B (zh) | 2016-12-07 | 2016-12-07 | 基于能量分析的乘员约束系统优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106599430B (zh) |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107063718A (zh) * | 2017-06-16 | 2017-08-18 | 吉林大学 | 汽车正面碰撞波形参数化评价方法 |
| CN107577843B (zh) * | 2017-08-02 | 2021-10-01 | 吉林大学 | 碰撞波形与约束系统特性耦合关系的评价方法 |
| CN108229079A (zh) * | 2018-03-20 | 2018-06-29 | 北京海纳川汽车部件股份有限公司 | 车辆的车体和约束系统匹配的设计方法和车辆 |
| CN108625636A (zh) * | 2018-04-30 | 2018-10-09 | 罗杰敏 | 一种智慧城市安全交通系统 |
| CN108664910A (zh) * | 2018-04-30 | 2018-10-16 | 许晓山 | 校园学生接送管理系统 |
| CN108608982B (zh) * | 2018-04-30 | 2020-12-18 | 苏州天浩汽车科技股份有限公司 | 车辆安全系统 |
| CN109506944B (zh) * | 2018-10-16 | 2020-05-19 | 中国汽车技术研究中心有限公司 | 一种用于测试车辆正面对碰中攻击性的方法 |
| CN109543259B (zh) * | 2018-11-09 | 2023-03-31 | 中国汽车技术研究中心有限公司 | 一种构建等效实车碰撞波形的方法 |
| CN109522647A (zh) * | 2018-11-15 | 2019-03-26 | 中国直升机设计研究所 | 一种直升机乘员座椅约束系统符合性验证方法 |
| CN109766579A (zh) * | 2018-12-10 | 2019-05-17 | 山东国金汽车制造有限公司 | 一种基于cae技术的汽车乘员约束系统优化方法 |
| CN111695189A (zh) * | 2019-03-13 | 2020-09-22 | 东汉新能源汽车技术有限公司 | 一种乘员约束系统的碰撞优化设计方法和装置 |
| CN110309556B (zh) * | 2019-06-13 | 2023-07-28 | 上汽大众汽车有限公司 | 汽车碰撞使用的头部碰撞模型的有限元建模方法及其模型 |
| CN110377988B (zh) * | 2019-07-08 | 2020-05-22 | 清华大学 | 乘员道路交通事故损伤预测方法 |
| CN110766213A (zh) * | 2019-10-15 | 2020-02-07 | 清华大学 | 汽车乘员碰撞损伤预测方法及计算机设备 |
| CN110826275B (zh) * | 2019-10-31 | 2020-07-14 | 清华大学 | 车辆前舱罩盖设计参数的优化方法 |
| CN110955992A (zh) * | 2019-11-21 | 2020-04-03 | 天津科技大学 | 六岁儿童乘员的有限元模型的构建方法及系统、儿童约束系统评价方法及系统 |
| CN111353220A (zh) * | 2020-02-24 | 2020-06-30 | 西南交通大学 | 一种列车碰撞能量分配快速分析方法 |
| CN117574539B (zh) * | 2023-11-22 | 2024-07-26 | 苏州郅荣软件有限公司 | 一种车辆约束系统优化参数验证系统 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1843804A (zh) * | 2005-12-07 | 2006-10-11 | 江苏大学 | 机动车碰撞安全座椅 |
| CN202703345U (zh) * | 2012-03-13 | 2013-01-30 | 江苏大学 | 弹簧驱动式主动响应头枕 |
| CN103310119A (zh) * | 2013-07-05 | 2013-09-18 | 湖南大学 | 基于损伤评价的车辆碰撞乘员约束系统摩擦特性的反分析方法 |
-
2016
- 2016-12-07 CN CN201611113305.9A patent/CN106599430B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN106599430A (zh) | 2017-04-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106599430B (zh) | 基于能量分析的乘员约束系统优化方法 | |
| CN103310119A (zh) | 基于损伤评价的车辆碰撞乘员约束系统摩擦特性的反分析方法 | |
| CN110851958A (zh) | 碰撞严重性的预测方法 | |
| CN107577843B (zh) | 碰撞波形与约束系统特性耦合关系的评价方法 | |
| Yang et al. | Design and injury analysis of the seated occupant protection posture in train collision | |
| Huang et al. | Optimisation study of occupant restraint system concerning variations in occupant size and crash severity in frontal collisions | |
| Hongling et al. | Combined with finite element analysis of car seat safety performance improvement | |
| CN117990389A (zh) | 基于波形等效的台车侧面碰撞方法 | |
| Kuznetcov et al. | Influence of non-standard seating postures on safety of motorcoach passengers in frontal collisions: a numerical study | |
| Kuznetcov | Effect of acceleration time history in motorcoach frontal collision on passenger safety under uncertainty of seating posture | |
| Ji et al. | Influence of school bus seat spacing on crash protection of six–year olds | |
| Soliman | Finite element modelling of vehicles and human body models for passive safety | |
| Ye et al. | Automated analysis of driver response in a finite element crash test reconstruction | |
| Lin et al. | Usage of LSTC_NCAC hybrid III 50th dummy in frontal occupant simulation | |
| Martínez et al. | Analysis of coaches rows seats distance influence on the passengers comfort and safety | |
| Kuwahara et al. | Research of Occupant kinematics and Injury values of Hybrid III, THOR, and human FE model in Oblique Frontal Impact | |
| He et al. | A study on rear seat occupant injuries in Side impact | |
| Fung et al. | Main effects of driving posture changes in frontal crashes due to ageing | |
| Sankarasubramanian et al. | Optimization of vehicle front for safety of pedestrians | |
| Zhang et al. | Development of a biofidelic rollover dummy-part II: validation of the kinematic response of THOR multi-body and finite element models relative to response of the physical THOR dummy under laboratory rollover conditions | |
| Yu | Finite element analysis of passenger multiple belt restraint configurations | |
| Ge et al. | Research of front seat design parameters on rear seat passenger safety | |
| CN113954991A (zh) | 一种约束系统的参数获取方法及装置 | |
| CN118010378A (zh) | 一种用于膝部保护的测试方法及测试系统 | |
| Naughton et al. | An approach to front-end system design for pedestrian safety |
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: 20200331 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |