1. 引言
座椅是非常重要的驾乘界面。人体–座椅之间的体压分布常用于检验座椅静态舒适性 [1] [2] [3]。一般认为,良好的体压分布特性是保证座椅为人体提供舒适的乘坐界面的必要条件。因此,座椅开发后期通常制作座椅样件,并测量样件的体压分布,据此评判座椅静态舒适性的优劣。然而,这种传统的座椅开发模式只能在设计后期才能根据检测到的体压分布来评价静态舒适性,一旦设计方案不理想,则由于优化设计方案造成的设计周期大大加长,势必增加了开发成本。
体压分布仿真技术是利用少量的座椅几何轮廓、海绵发泡轮廓和大致的骨架结构信息即能够对体压分布进行快速预估的仿真方法。采用这种方法能够在设计的早期就根据有限的设计资料进行仿真,从而能够很早就预估体压分布的情况,并评判静态舒适性的好坏;此外,还能方便地评估不同的海绵和蒙皮材料、海绵形状、骨架构造形式等因素可能对体压分布的影响,从而在设计早期就综合权衡多方面因素后选定最佳设计方案。
关于体压分布仿真,目前行业中普遍采用有限元方法来建立人体有限元模型,其中一般包括骨骼、体表、软组织三层来模拟人体 [4] [5] [6] [7]。人体骨骼和体表通常采用线弹性材料进行模拟,软组织则通常采用超弹性Mooney-Rivlin材料模拟。目前的研究还存在姿势不准确、人体材料不够准确、采用单一尺寸的人体进行仿真分析等问题。很明显,人体姿势是影响仿真准确性的重要因素,文献 [6] 考虑了姿势因素来模拟体压分布,但模拟休息状态驾驶员或者其他乘员还需要进一步建立相应的姿势预测模型来进行模拟。人体软组织材料不是均匀的。事实上,人体肌肉组织和一般的皮下软组织的力学特性差异很大,因此,在模拟的时候就不能够泛泛地将软组织统一赋予一种材料属性。虽然实际应用中需要测量座椅样件的体压分布,但通常也都是需要测量不同尺寸人体所产生的体压分布,并综合考虑来对座椅静态舒适性进行评价的,但目前的体压分布仿真往往还都采用单一尺寸(例如:95百分位身高男子的人体尺寸)人体模型进行仿真,其结果对于体压分布的评价说服力还相对欠缺 [4] [5] [6] [7]。在座椅建模方面,目前主要考虑海绵发泡的材料特性进行模拟,而蒙皮材料则考虑较少。蒙皮材料在座椅设计中经常更改材料类别,例如从织物材料改为PVC涂层面料,以获得更为理想的外观效果等。不同类别的蒙皮材料其力学特性有较大的差异,而且在座椅装配中为保证美观还施加了张紧力,这些都会显著影响体压分布。因此,在体压分布仿真中不能忽略蒙皮材料。
本文研究了考虑蒙皮材料、考虑不同尺寸人体的座椅体压分布仿真和评价,通过与实验数据的对比,验证了本文所建立的仿真分析方法的合理性、准确性、实用性。
2. 座椅蒙皮材料建模
座椅采用的蒙皮材料为PVC涂层面料,其力学特性主要取决于基布材料。通常基布材料经纬向力学特性存在差异,考虑到在正常使用范围内单一方向的力学特性近似可以视为线弹性,于是可采用正交各向异性线弹性(orthotropic linear elasticity)材料来进行模拟,其应力应变关系为 [8] [9]:
(1)
其中,
,为应力;
,为应变;而Q为刚度矩阵,且
(2)
其中,Ex、Ey分别为经纬向的杨氏模量;vxy、vyx为泊松比;
;Gxy为切变模量。经纬向杨氏模量、泊松比、切变模量可根据经纬向拉伸试验(图1)和相框剪切试验(图2)确定。图3为通过仿真得到的蒙皮材料拉伸曲线,他们都经过了试验验证,仿真的蒙皮材料特性是准确的。
(a)
(b)
(c)
Figure 1. Tension test of cover materials. (a) Test apparatus; (b) 0% elongation; (c) 40% elongation
图1. 蒙皮材料拉伸试验。(a) 试验设备;(b) 拉伸0%;(c) 拉伸40%
(a)
(b)
Figure 2. Trellising shearing test of cover materials. (a) Trellising shearing frame; (b) Test process
图2. 蒙皮材料相框剪切试验。(a) 相框剪切试验夹具;(b) 试验过程
(a)
(b)
Figure 3. Unidirectional tension curves of cover material. (a) Force-translation curve; (b) Stress-strain curve
图3. 蒙皮单向拉伸曲线。(a) 力–位移曲线;(b) 应力–应变曲线
3. 人体生物力学建模
3.1. 人体软组织材料
对于体压分布数值模拟,影响较大的因素之一是软组织材料的力学特性。人体软组织一般用非线性、超弹性、不可压缩的Mooney-Rivlin材料模拟 [4] [5] [6] [7],其应变能函数为:
(3)
其中,U为应变能;
和
分别为Cauchy-Green等容形式张量的第一和第二不变量;J为材料变形前后的体积比;C10、C01、D1为材料参数,应用时需要针对不同部位的软组织赋予不同的数值。参考前人的研究成果,对于一般的皮下软组织,三者取值分别为3.65 kPa、3.35 kPa和5.71,对于肌肉组织可以适当加大刚度来模拟。
3.2. 人体有限元建模
为了能反映目标人群中不同尺寸的人的体压分布特性,原则上应该建立不同身材的人体模型来模拟;但由于几何和有限元建模工作量非常巨大,并考虑到行业里常常用95、50、5百分位身高人体的体压分布来评价座椅静态舒适性,于是建立了P95男子、P50男子和女子、P5女子人体模型,见图4。人体数据来自最新中国标准化研究院2016年的人体尺寸普查数据。
根据身高、坐高和体重,利用人体尺寸预测模型(见文献 [10] )计算人体尺寸,参照购买的人体体表点云和人体切片数据,利用几何放缩的方法建立人体几何模型,并利用姿势模型预测和设定人体模型的姿势 [11]。基于几何模型建立有限元模型(图5),包括皮肤、骨骼和软组织三部分,在骨骼和皮肤之间的部分都简化为软组织。皮肤和骨骼材料的定义参见表1 [4] [5] [6] [7]。
(a)
(b)
(c)
Figure 5. Geometric and FE models of seated human body. (a) Geometric model of bone; (b) Geometric model of soft-tissue and skin; (c) FE model
图5. 坐姿人体几何和有限元模型。(a) 骨骼几何模型;(b) 软组织和体表几何模型;(c) 有限元模型
![](Images/Table_Tmp.jpg)
Table 1. Material properties of skin and bone
表1. 皮肤和骨骼材料属性
人体骨架的各部分肢体之间需要定义关节。此外,骨骼-软组织之间、软组织-体表之间采用绑定的方式装配。
4. 体压分布仿真
4.1. 座椅有限元建模
根据座椅几何模型来建立有限元模型,包括蒙皮、软海绵、主海绵、弹簧、骨架。蒙皮采用壳单元模拟,其他均采用体单元模拟。由于很多的细节部分包含很多的不规则形状,因此,均自动划分网格。上述这些部分的材料属性的设置参见表2。
![](Images/Table_Tmp.jpg)
Table 2. Settings of seat material parameters
表2. 座椅材料参数的设置
4.2. 边界条件
人体模型建立好之后,需要和座椅装配在一起。考虑到直接装配到落座后的位置会产生人体和座椅相互穿透的现象而导致计算出错,因此,也需要将人体模型事先相对于座椅产生一个位移(分解为前后向位移X_offset和垂向位移Z_offset),参见下图6所示。在仿真计算的时候,同样先施加前后向的强制位移,然后再施加垂向的重力场,使得人体正常下落,完成落座过程。人体和座椅接触面之间的切向属性采用“Penalty”(罚函数)进行设置,摩擦系数设置为0.3 (具体视蒙皮材料而有所变化);接触压力与穿透之间的关系通过法向属性定义,一般选用“Penalty”(罚函数)约束增强方法。为了尽可能避免穿透,在接触表面采用主从面权重算法,根据主-从面设置原则确定座椅表面和人体皮肤的主从面。为保证在模拟过程中座椅整体不产生运动,需要约束骨架的六个方向自由度;如果省略了骨架,则需要约束座垫、靠背和头枕背面六个方向的自由度。重力加载能够准确地模拟人体落座过程,但是对有限元模型的精确性要求较高,需要事先准确测量人体各部分的体积,并设置合理的密度,以保证人体模型以实际的重力进行落座。
![](//html.hanspub.org/file/4-2340967x27_hanspub.png?20220209164530754)
Figure 6. Assembly of human body and seat
图6. 人体与座椅的装配
在蒙皮缝线处施加朝向海绵方向的拉力,目的是使得蒙皮绷紧。
考虑到能够分析不同座椅(不同人体坐姿)的体压分布,将脚踝相对于胯点的位置参数化,并设置计算步(Step)专门用于调节脚踝相对于胯点的位置,从而使下肢实现不同的姿势。
为了能够分析座椅上背部区域和头枕设计是否适合头颈部舒适性要求,增加了头部靠在头枕上的计算步,通过施加扭矩在颈部使得头向后靠,该扭矩通过骨肌生物力学计算(图7)获得 [12] [13]。建立了头颈部姿势舒适性评价模型来评价头部靠在头枕上之后的姿势舒适性。
![](//html.hanspub.org/file/4-2340967x28_hanspub.png?20220209164530754)
Figure 7. Calculation of head-neck joint moments
图7. 头颈部关节扭矩计算
4.3. 仿真结果及验证
按照前述方法分别计算P95、P50男、P50女、P5女的体压分布。限于篇幅,仅给出P95男子的仿真和实测结果。仿真计算结果参见图8,峰值压强1.25 kPa左右,出现在坐垫上坐骨下面的区域;由于模拟中考虑的人体背部姿态是自然休息的姿态,腰背部都有压强分布。图9为对应的体压分布实测结果,表3为主要分区的压力比率和峰值压强(定义见下节)的仿真和试验得到的置信区间的对比结果,可见,仿真和实测结果很接近。
![](//html.hanspub.org/file/4-2340967x29_hanspub.png?20220209164530754)
Figure 8. Simulation result of the interface pressure distribution of man with the percentile-95 stature
图8. 身高95百分位男子的体压分布仿真结果
![](//html.hanspub.org/file/4-2340967x30_hanspub.png?20220209164530754)
Figure 9. Test result of interface pressure distribution
图9. 体压分布实测结果
![](Images/Table_Tmp.jpg)
Table 3. Comparison of the simulation and test results for load ratio and peak pressure
表3. 压力比率和峰值压强仿真和试验对比
4.4. 结果的评价
4.4.1. 体压分布的评价
目前主流的体压分布评价都是基于座椅表面分区的 [14]。典型的分区方法参见图10。
(a)
(b)
Figure 10. Regions of seat surface. (a) Seat cushion; (b) Seat back
图10. 座椅表面分区。(a) 座垫;(b) 靠背
根据座椅分区,统计得到座椅上确定部位的体压分布统计特征值,包括峰值压强、平均压强、压力比率。若Pkmax为区域k的峰值压强;Pij为第i行第j列压强值;MAX{…}为最大值函数;rk为区域k的压力比率;Lk为区域k的总压力;L为区域所处的座垫或靠背的总压力,则
(4)
(5)
(6)
(7)
对典型的身高百分位的一定数量的被试人员舒适驾驶姿势的体压分布进行测量,并计算座椅各分区的峰值压强、平均压强、压力比率;再统计上述参量的90%的置信区间,作为舒适座椅的评价标准。如果新设计方案仿真得到的上述参量落入该90%置信区间,就认为是舒适的压力分布。
4.4.2. 颈部姿势的评价
首先需要评价标准。选取一定数量被试人员,通过测量统计(图11)得到他们的舒适颈部姿势角度(图12)的90%置信区间,作为舒适的颈部姿势角度范围。
体压分布仿真得到人体落座稳定后的头颈和颈胸关节、以及胯点坐标,利用这三点坐标计算颈部姿势角度,并对比舒适的颈部姿势角度范围来作出评价。
![](//html.hanspub.org/file/4-2340967x37_hanspub.png?20220209164530754)
Figure 11. Test for comfort neck posture angles
图11. 颈部姿势角度实验
体压分布仿真得到人体落座稳定后的头颈和颈胸关节、以及胯点坐标,利用这三点坐标计算颈部姿势角度,并对比舒适的颈部姿势角度范围来作出评价。
5. 结论
1) 在体压分布仿真中考虑了座椅蒙皮的因素,并研究了蒙皮材料的有限建模方法,对于分析蒙皮材料的变化对体压分布的影响有帮助;
2) 根据最新的中国人体尺寸建立了典型尺寸的人体生物力学模型,可用于面向中国人群进行体压分布仿真和座椅静态舒适性分析;
3) 研究了适用于不同车型的体压分布仿真、分析和评价方法,其中评价包括体压分布本身的评价和人体颈部姿势的评价;
4) 本文所建立的体压分布仿真、评价方法自成体系,具有较高的应用价值。
NOTES
*通讯作者。