高開(kāi)展,羅 巧,張志飛,徐中明
(重慶大學(xué)機(jī)械與運(yùn)載工程學(xué)院,重慶 400030)
汽車座椅作為與駕乘人員接觸的主要人機(jī)交互界面,是影響汽車乘坐舒適性的重要因素[1]。對(duì)于座椅舒適性的評(píng)價(jià)主要分為兩個(gè)方面:①主觀評(píng)價(jià)方法,如語(yǔ)義差分法[2]、舒適性量表法[3]、舒適性問(wèn)卷法[4]等;②客觀評(píng)價(jià)方法,如體壓分布測(cè)量[5]、加速度測(cè)量[6]、生物肌電力測(cè)量[7]以及姿勢(shì)分析[8]等。其中,加權(quán)加速度均方根值作為量化全身振動(dòng)對(duì)人類健康和舒適性影響的指標(biāo)廣泛應(yīng)用于坐姿人體振動(dòng)評(píng)價(jià)[9]和汽車座椅舒適性評(píng)價(jià)[10]。
近年來(lái),隨著測(cè)試技術(shù)的發(fā)展,人體與座椅界面間的壓力測(cè)量作為量化人-椅界面間相互作用的方法,在汽車座椅的舒適性評(píng)價(jià)和設(shè)計(jì)過(guò)程中逐漸得到應(yīng)用[11]。通過(guò)對(duì)靜態(tài)工況下人-椅接觸面的壓力分布進(jìn)行研究,發(fā)現(xiàn)坐姿人體的不舒適度與平均壓力[12-14]、最大壓力[15]、壓力梯度[16]或坐骨結(jié)節(jié)附近的最大壓力有關(guān)[17-18]。Uenishi等[19]通過(guò)駕駛狀態(tài)下的體壓分布測(cè)量對(duì)比3種不同的泡沫座椅,發(fā)現(xiàn)平均壓力變化率均方根是評(píng)價(jià)振動(dòng)不舒適較為合適的指標(biāo)。Ahmadian等[20]根據(jù)人-椅界面間壓力大小對(duì)人體細(xì)胞的生理影響不同,對(duì)不同范圍內(nèi)的壓力進(jìn)行加權(quán),提出區(qū)域壓力均方根的概念,以評(píng)估長(zhǎng)時(shí)間乘坐對(duì)人體健康造成的影響。Na等[21]發(fā)現(xiàn)駕駛工況下的體壓分布能夠反映駕駛員姿勢(shì)變化,并依此來(lái)預(yù)測(cè)駕駛員的移動(dòng)頻率。Tahir等[22]通過(guò)對(duì)駕駛工況下不同身高和體質(zhì)量的駕駛員壓力分布進(jìn)行定量分析,用于指導(dǎo)符合人機(jī)工程學(xué)的座椅舒適性設(shè)計(jì)。Lantoine等[5]基于對(duì)駕駛工況的模擬,探究了不同座椅軟硬程度和道路類型對(duì)壓力分布和不適感的影響,以幫助改進(jìn)座椅靠背和坐墊的舒適性設(shè)計(jì)。
盡管體壓分布能更全面地描述人-椅界面的相互作用,但其應(yīng)用主要集中在靜態(tài)工況,對(duì)于振動(dòng)激勵(lì)這一動(dòng)態(tài)工況下的研究進(jìn)展緩慢。Wu等[23]通過(guò)對(duì)1-10 Hz正弦振動(dòng)下人體與剛性座椅接觸面的壓力分布進(jìn)行研究,發(fā)現(xiàn)在整個(gè)頻率范圍內(nèi)壓力幅值隨著振動(dòng)激勵(lì)強(qiáng)度的增加而增加。此外,Wu等[17]還發(fā)現(xiàn)正弦振動(dòng)下人體與泡沫座椅接觸面間坐骨結(jié)節(jié)處的最大壓力同樣隨著振動(dòng)激勵(lì)強(qiáng)度的增加而增加。Inagaki等[24]將壓力隨時(shí)間的變化率納入1 Hz及以下的座椅振動(dòng)舒適的評(píng)價(jià)中,但其應(yīng)用較有限。Hinz等[25]對(duì)3種振動(dòng)激勵(lì)強(qiáng)度下的壓力分布進(jìn)行測(cè)量,發(fā)現(xiàn)無(wú)論是平均壓力還是最大壓力都隨激勵(lì)強(qiáng)度的增大而增大。在振動(dòng)舒適性評(píng)價(jià)中,振動(dòng)激勵(lì)強(qiáng)度是影響人體振動(dòng)響應(yīng)特性、舒適性感受的重要因素[26-27],同樣也會(huì)顯著改變?nèi)?椅接觸面的體壓分布[23]。盡管振動(dòng)激勵(lì)下的壓力分布會(huì)對(duì)乘坐舒適性產(chǎn)生影響[28],但目前對(duì)于振動(dòng)舒適性的評(píng)價(jià)主要依賴于加速度測(cè)試[6,29-30],體壓分布并未充分發(fā)揮其作用,值得進(jìn)一步探究。
為探究如何借助體壓分布測(cè)量來(lái)表征汽車座椅的振動(dòng)舒適性,設(shè)計(jì)汽車座椅振動(dòng)舒適性主客觀評(píng)價(jià)試驗(yàn)。以垂向振動(dòng)為激勵(lì),在不同振動(dòng)激勵(lì)強(qiáng)度下開(kāi)展振動(dòng)加速度測(cè)試、體壓分布測(cè)試和主觀評(píng)價(jià)。利用非參數(shù)統(tǒng)計(jì)分析方法和史蒂文斯冪定律分析振動(dòng)加速度和體壓分布指標(biāo)與主觀評(píng)分之間的相關(guān)程度和關(guān)聯(lián)性大小,為汽車座椅振動(dòng)舒適性評(píng)價(jià)提供參考。
選用某款汽車座椅固定在振動(dòng)臺(tái)臺(tái)面,如圖1所示。試驗(yàn)采用的振動(dòng)臺(tái)為北京思齊致新科技公司生產(chǎn)的電磁激勵(lì)單軸振動(dòng)臺(tái),臺(tái)面最大負(fù)載3 kN,最大行程±150 mm,工作頻率范圍為0.1-50 Hz。利用SIControl分析系統(tǒng)對(duì)試驗(yàn)所用的6種垂向激勵(lì)信號(hào)進(jìn)行迭代重現(xiàn),頻率帶寬0.5-20 Hz,加速度均方根值為0.2、0.4、0.6、0.8、1.0、1.2 m/s2,保證激勵(lì)信號(hào)的重現(xiàn)誤差小于1%。
圖1 振動(dòng)舒適性試驗(yàn)
利用IMC數(shù)據(jù)采集器、B&K 4515-B三向加速度傳感器分別采集座椅椅面和靠背的加速度,采樣頻率設(shè)置為200 Hz,采樣時(shí)間為60 s。采用XSENSOR體壓測(cè)量系統(tǒng)對(duì)人-椅接觸面的壓力進(jìn)行采集,如圖2所示。其中靠背用體壓墊型號(hào)為PX100:40.64,有效測(cè)量面積為50.8×81.2 cm2,量程1.4-27 kPa;坐墊用體壓墊型號(hào)為L(zhǎng)X210:40.40,有效測(cè)量面積為50.8×50.8 cm2,量程0.69-103.43 kPa;體壓分布測(cè)試的采樣頻率為32 Hz。
圖2 體壓分布測(cè)量試驗(yàn)
本次試驗(yàn)招募了12名受試者,受試者特征如表1所示。要求全體受試者近期無(wú)身體不適,且近3年內(nèi)無(wú)腰背痛或肌肉骨骼疾病,所有受試者均簽署了知情同意書(shū)。試驗(yàn)開(kāi)始前對(duì)每位受試者進(jìn)行簡(jiǎn)單指導(dǎo),并對(duì)不同振動(dòng)激勵(lì)逐個(gè)體驗(yàn),以幫助受試者熟悉試驗(yàn)內(nèi)容和流程。此外,試驗(yàn)開(kāi)始前每位受試者須在汽車座椅上靜坐3 min,以確保泡沫坐墊的松弛特性達(dá)到相對(duì)穩(wěn)定狀態(tài)。
表1 受試者特征
開(kāi)展振動(dòng)加速度測(cè)試時(shí),調(diào)整并固定座椅靠背與坐墊夾角為110°,然后正確放置、固定坐墊加速度傳感器。受試者以放松姿勢(shì)坐在座椅上,背靠椅背,手自然放置于大腿上,將腳放在振動(dòng)臺(tái)的腳踏板上并保持大腿水平,如圖1所示。振動(dòng)臺(tái)隨機(jī)播放(隨機(jī)序號(hào)在Matlab中生成)6種不同加速度均方根值的垂向激勵(lì)信號(hào)。為保證測(cè)量數(shù)據(jù)的一致性,每次測(cè)量重復(fù)進(jìn)行3次,且保證3次測(cè)量誤差不超過(guò)5%。此外,要求受試者全程佩戴安全帶,如有不適可按下急停按鈕以終止試驗(yàn)。
體壓分布測(cè)試與振動(dòng)加速度測(cè)試分開(kāi)進(jìn)行,要求正確放置、固定體壓墊,保證表面平整無(wú)褶皺。受試者身穿較柔軟的衣服,采用與振動(dòng)加速度測(cè)試時(shí)相同的坐姿,如圖2所示。體壓數(shù)據(jù)采集前先根據(jù)體壓圖像進(jìn)行調(diào)試,待體壓數(shù)據(jù)穩(wěn)定后對(duì)靠背和坐墊壓力數(shù)據(jù)同時(shí)進(jìn)行采集,每次測(cè)量重復(fù)進(jìn)行3次。
主觀評(píng)價(jià)與體壓測(cè)試同時(shí)進(jìn)行,采用5分制不舒適性評(píng)價(jià)量表(見(jiàn)表2),對(duì)振動(dòng)引起的不適進(jìn)行打分。測(cè)試時(shí),隨機(jī)播放6種不同強(qiáng)度的垂向激勵(lì)信號(hào),受試者根據(jù)體驗(yàn)對(duì)指定振動(dòng)信號(hào)給出對(duì)應(yīng)不舒適度分值,并做記錄。為避免視覺(jué)和噪聲對(duì)受試者主觀感受造成影響,要求受試者在試驗(yàn)過(guò)程中帶上耳機(jī)(播放白噪聲),同時(shí)閉上雙眼。若受試者對(duì)自己的判斷不確定,可以要求重新進(jìn)行測(cè)試。
表2 不舒適性評(píng)價(jià)量表
1.3.1 振動(dòng)加速度
ISO 2631-1:1997將加權(quán)加速度均方根作為基本的評(píng)價(jià)方法,用來(lái)評(píng)價(jià)振動(dòng)對(duì)人體舒適性的影響[29]。該指標(biāo)綜合考慮了振動(dòng)強(qiáng)度、頻率和振動(dòng)方向?qū)θ梭w舒適性的影響。加權(quán)加速度均方根的計(jì)算方法如下。
首先將采集到的加速度時(shí)間歷程x(n)進(jìn)行頻譜分析得到功率譜密度函數(shù)G(af),其頻譜X(k)計(jì)算的離散表達(dá)式為
式中:N為傅里葉變換的點(diǎn)數(shù);k代表離散的頻率點(diǎn)。
功率譜密度函數(shù)G(af)為
式中Fs為采樣頻率。
則各軸的加權(quán)加速度均方根為
式中w(f)為頻率加權(quán)函數(shù),其中xs、ys、yb、zb方向的加權(quán)函數(shù)為w(df),zs的加權(quán)函數(shù)為w(kf),xb的加權(quán)函數(shù)為w(cf)。
加權(quán)加速度均方根值的計(jì)算公式為
式中i=1,2。當(dāng)i=1時(shí),av1代表座椅椅面處的加權(quán)加速度均方根值,此時(shí)軸加權(quán)系數(shù)kx、ky、kz均為1.0;當(dāng)i=2時(shí),av2代表座椅靠背處的加權(quán)加速度均方根值,此時(shí)軸加權(quán)系數(shù)kx、ky、kz分別為0.8、0.5、0.4[29]。
則同時(shí)考慮椅面和靠背處的加權(quán)加速度均方根值的計(jì)算公式為
1.3.2 體壓分布
平均壓力均值指人-椅接觸面間不同時(shí)刻全部接觸點(diǎn)壓力在時(shí)間和空間上的算術(shù)平均值,即
最大壓力均值指人-椅接觸面間不同時(shí)刻所有接觸點(diǎn)中最大壓力在時(shí)間上的算術(shù)平均值,即
式中pip,j=max(pi,j),pip,j代表測(cè)量周期內(nèi)第j幀的最大壓力。
壓力變化率均方根值是指振動(dòng)條件下壓力在指定時(shí)間歷程內(nèi)變化率的均方根[19],即
式中:T代表總的時(shí)間周期;p(t)代表平均壓力時(shí)間歷程。
文獻(xiàn)[31]表明,振動(dòng)下坐姿人體的不舒適感與人-椅界面間的力有關(guān)。由于XSENSOR體壓測(cè)量系統(tǒng)所測(cè)得的壓力為法向壓力,本文中結(jié)合人-椅接觸面的接觸面積,提出法向力變化率均方根指標(biāo),即法向壓力和接觸面積的乘積在指定時(shí)間歷程內(nèi)變化率的均方根為
式中F(t)=p(t)×A(t),其中A(t)為接觸面積時(shí)間歷程。
根據(jù)式(6)~式(9)分別在靠背面、坐墊面計(jì)算平均壓力均值、最大壓力均值、壓力變化率均方根和法向力變化率均方根,再計(jì)算人-椅接觸面的總體指標(biāo),即
式中:SPcrms代表坐墊處的平均壓力變化率均方根值;BPcrms代表靠背處的平均壓力變化率均方根值;代表人-椅接觸面的平均壓力變化率均方根值。
式中:SFcrms代表坐墊處的法向力變化率均方根值;BFcrms代表靠背處的法向力變化率均方根值;代表人-椅接觸面的法向力變化率均方根值。
使用非參數(shù)檢驗(yàn)對(duì)多個(gè)相關(guān)樣本進(jìn)行差異顯著性檢驗(yàn)。Friedman檢驗(yàn)對(duì)分析樣本無(wú)正態(tài)分布的要求,故采用Friedman雙向方差分析對(duì)12組樣本進(jìn)行非參數(shù)統(tǒng)計(jì)學(xué)分析,根據(jù)顯著性因子p是否小于0.05來(lái)判斷測(cè)試樣本是否存在顯著差異。使用Spearman秩相關(guān)分析來(lái)確定客觀指標(biāo)和不舒適評(píng)分之間的相關(guān)程度,將顯著性因子p與0.05比較,判斷測(cè)試變量是否顯著相關(guān),根據(jù)相關(guān)系數(shù)β的大小判斷相關(guān)性的強(qiáng)弱。
12名受試者在不同振動(dòng)激勵(lì)強(qiáng)度下的主觀評(píng)分結(jié)果如圖3所示。從圖中可以看出,隨著振動(dòng)激勵(lì)強(qiáng)度的逐漸增大,受試者的主觀不舒適性評(píng)分也隨之增加。對(duì)于同一振動(dòng)加速度均方根值而言,受試者之間的主觀不舒適性評(píng)分有所不同,說(shuō)明個(gè)體間的差異可能導(dǎo)致不同受試者對(duì)同一振動(dòng)激勵(lì)的主觀感受不同。差異顯著性水平統(tǒng)計(jì)學(xué)分析結(jié)果表明,振動(dòng)激勵(lì)強(qiáng)度的變化對(duì)受試者主觀不舒適度有顯著影響(Friedman,p<0.01)。因此,隨著振動(dòng)激勵(lì)強(qiáng)度的增大,經(jīng)由汽車座椅傳遞給人體的振動(dòng)隨之增大,對(duì)于同一個(gè)座椅和固定坐姿而言,振動(dòng)輸入激勵(lì)強(qiáng)度的改變是導(dǎo)致舒適性變化的根本原因。
圖3 不同振動(dòng)激勵(lì)強(qiáng)度下的主觀不舒適性評(píng)分
12名受試者在不同振動(dòng)激勵(lì)強(qiáng)度下加權(quán)加速度均方根值分析結(jié)果如圖4所示。同樣地,隨著振動(dòng)激勵(lì)強(qiáng)度的逐漸增大,人-椅接觸面的加權(quán)加速度均方根值隨之增加。差異顯著性水平統(tǒng)計(jì)學(xué)分析結(jié)果表明,輸入激勵(lì)強(qiáng)度的變化對(duì)人-椅接觸面的加權(quán)加速度均方根值有顯著影響(Friedman,p<0.01)。
圖4 不同振動(dòng)激勵(lì)強(qiáng)度下加權(quán)加速度均方根值
圖5展示了某一受試者在垂向激勵(lì)工況下10、20、30 s時(shí)刻(從左至右)體壓分布的測(cè)試結(jié)果。從圖中可以看出,靠背最大壓力均出現(xiàn)在肩胛骨附近,坐墊最大壓力均出現(xiàn)在人體坐骨結(jié)節(jié)附近,以坐骨結(jié)節(jié)為中心向外逐漸減小,從坐骨結(jié)節(jié)至大腿前沿,壓力逐漸減小。對(duì)不同時(shí)刻肩胛骨和坐骨結(jié)節(jié)區(qū)域的壓力分布進(jìn)行局部放大,從中可以看出靠背壓力和坐墊處壓力在不同時(shí)刻有明顯變化,其中坐墊處的壓力分布形式與文獻(xiàn)[17]報(bào)道的結(jié)果一致。
圖5 某一受試者不同時(shí)刻的體壓分布云圖
對(duì)12名受試者體壓分布時(shí)間歷程中的50 s數(shù)據(jù)進(jìn)行處理,分別得到不同振動(dòng)激勵(lì)強(qiáng)度下各個(gè)受試者與座椅接觸面的平均壓力均值、最大壓力均值、平均壓力變化率均方根值和法向力變化率均方根值,如圖6~圖9所示。從圖中可以看出,隨著振動(dòng)激勵(lì)強(qiáng)度的增大,人-椅接觸面的平均壓力均值和最大壓力均值都有所變化,但規(guī)律不明顯,而平均壓力變化率均方根值和法向力變化率均方根值則逐漸增大。差異顯著性水平統(tǒng)計(jì)學(xué)分析結(jié)果表明,輸入激勵(lì)強(qiáng)度的變化對(duì)人-椅接觸面的平均壓力均值和最大壓力均值的變化無(wú)顯著影響(Friedman,p>0.05),而對(duì)平均壓力變化率均方根值和法向力變化率均方根值均有顯著影響(Friedman,p<0.01)。
圖6 不同振動(dòng)激勵(lì)強(qiáng)度下的平均壓力均值
圖9 不同振動(dòng)激勵(lì)強(qiáng)度下的法向力變化率均方根值
從以上分析結(jié)果可以看出,不同的振動(dòng)激勵(lì)強(qiáng)度對(duì)人-椅接觸面的平均壓力變化率均方根值和法向力變化率均方根值的影響與對(duì)主觀不舒適性評(píng)分的影響趨勢(shì)一致,而平均壓力均值、最大壓力均值既無(wú)法反映振動(dòng)激勵(lì)強(qiáng)度的變化,也沒(méi)有呈現(xiàn)與主觀不舒適性評(píng)分類似的變化趨勢(shì)。
圖7 不同振動(dòng)激勵(lì)強(qiáng)度下的最大壓力均值
圖8 不同振動(dòng)激勵(lì)強(qiáng)度下的平均壓力變化率均方根值
為進(jìn)一步探究動(dòng)態(tài)工況下12名受試者的主觀不舒適性評(píng)分與客觀測(cè)試量的相關(guān)關(guān)系,對(duì)主觀評(píng)分、加速度、體壓分布等主客觀指標(biāo)進(jìn)行相關(guān)性分析,結(jié)果如表3所示。從表中可以看出,主觀不舒適性評(píng)分與人-椅接觸面的加權(quán)加速度均方根值、平均壓力變化率均方根值和法向力變化率均方根值均有顯著相關(guān)關(guān)系(Spearman,p<0.01),且相關(guān)系數(shù)均大于0.8。雖然不舒適評(píng)分與平均壓力均值在0.05級(jí)別(雙尾)存在相關(guān)關(guān)系,但相關(guān)系數(shù)β僅為0.26,說(shuō)明相關(guān)性并不強(qiáng)。而不舒適性評(píng)分與最大壓力之間不存在顯著相關(guān)性(Spearman,p>0.05)。對(duì)比表3中的分析結(jié)果可以發(fā)現(xiàn),受試者主觀不舒適性評(píng)分與人-椅界面間的平均壓力變化率均方根值以及法向力變化率均方根值的聯(lián)系更為緊密。
表3 主客觀指標(biāo)相關(guān)性分析
對(duì)于振動(dòng)加速度、體壓分布等客觀指標(biāo)和主觀不舒適性評(píng)分,采用史蒂文斯冪定律確定其關(guān)聯(lián)性。史蒂文斯冪定律描述了刺激的客觀物理量與受試者產(chǎn)生的感覺(jué)之間的關(guān)系[32],即
式中:φ代表振動(dòng)激勵(lì)的客觀量;ψ代表對(duì)應(yīng)激勵(lì)下的心理主觀物理量;k為一個(gè)常數(shù);n反映了人體不適隨客觀量的增長(zhǎng)率。
利用對(duì)數(shù)坐標(biāo)系將其轉(zhuǎn)化為線性關(guān)系:
用12名受試者的客觀測(cè)試量中位數(shù)代替客觀量φ,主觀不舒適性評(píng)分的中位數(shù)代替心理主觀物理量ψ,利用式(14)對(duì)主客觀物理量進(jìn)行線性擬合。擬合的效果通過(guò)擬合優(yōu)度進(jìn)行判斷,其表達(dá)式為
式中:yi表示第i個(gè)試驗(yàn)樣本;yˉ表示試驗(yàn)樣本的平均值;表示通過(guò)線性回歸得到的預(yù)測(cè)變量,一般情況下,R2取值在[0,1]之間,該值越大則表示模型擬合精度越好。
采用史蒂文斯冪定律對(duì)主客觀測(cè)試量進(jìn)行線性回歸擬合,相應(yīng)的擬合參數(shù)如表4所示。受試者不舒適評(píng)分與人-椅接觸面的加權(quán)加速度均方根值、平均壓力變化率均方根值和法向力變化率均方根值的關(guān)聯(lián)性分析結(jié)果如圖10~圖12所示??梢钥闯?,主觀不舒適性評(píng)分與加權(quán)加速度和體壓分布指標(biāo)均呈正相關(guān)關(guān)系,擬合優(yōu)度R2均大于98.0%,說(shuō)明人-椅接觸面的加權(quán)加速度均方根值、平均壓力變化率均方根值和法向力變化率均方根值與主觀不舒適性評(píng)分均有較高的關(guān)聯(lián)性。
表4 主客觀關(guān)聯(lián)性擬合參數(shù)
圖10 加權(quán)加速度均方根值與不舒適評(píng)分關(guān)聯(lián)性曲線
圖11 平均壓力變化率均方根與不舒適評(píng)分關(guān)聯(lián)性曲線
圖12 法向力變化率均方根與不舒適評(píng)分關(guān)聯(lián)性曲線
綜上所述,類似于利用人-椅界面間的速度對(duì)時(shí)間的變化率(即加速度)作為指標(biāo)描述人體在振動(dòng)下的不舒適感[29],由振動(dòng)導(dǎo)致的人-椅界面間的壓力或法向力對(duì)時(shí)間的變化率也可以作為汽車座椅振動(dòng)不舒適性評(píng)價(jià)的客觀度量。
以12名受試者為對(duì)象,設(shè)計(jì)6種不同的垂向激勵(lì)信號(hào),在振動(dòng)臺(tái)上對(duì)某款汽車座椅進(jìn)行了振動(dòng)舒適性主客觀評(píng)價(jià)試驗(yàn),并對(duì)加速度、體壓分布和主觀不舒適性評(píng)分等主客觀指標(biāo)進(jìn)行了相關(guān)性和關(guān)聯(lián)性分析,得出以下結(jié)論。
(1)靜態(tài)工況下常用的平均壓力和最大壓力等指標(biāo)與主觀不舒適性評(píng)分間的相關(guān)性較弱(β均小于0.3),因而不能用來(lái)反映振動(dòng)舒適性。
(2)人-椅接觸面的平均壓力變化率均方根值、法向力變化率均方根值與主觀不舒適性評(píng)分間的相關(guān)系數(shù)β均大于0.8、關(guān)聯(lián)性分析結(jié)果中擬合優(yōu)度R2均大于99.0%,與人-椅接觸面的加權(quán)加速度均方根值指標(biāo)相近,可以作為表征汽車座椅振動(dòng)舒適性的客觀評(píng)價(jià)指標(biāo)。