王廣彬,任金東,劉洪浩
(吉林大學(xué)汽車工程學(xué)院,長春 130022)
座椅舒適性是消費者考慮的一個重要因素。人椅界面間的體壓分布是度量座椅靜態(tài)舒適性的重要指標(biāo)[1-7]。研究表明,體壓分布與座椅舒適性明顯相關(guān),可綜合人椅間平均壓強、峰值壓強和壓強梯度等來評價座椅舒適性[8-11]。體壓分布通常通過實驗獲取,但存在精度低、可重復(fù)性差、影響因素多且復(fù)雜、很難得到人體組織應(yīng)力、人椅界面間剪力和關(guān)節(jié)肌肉負荷等信息,以及容易受主觀因素影響等缺點。相比之下,仿真方法不僅能克服實驗的不足,還具有成本低、周期短、不受主觀因素影響和能在設(shè)計早期應(yīng)用等優(yōu)點。
關(guān)于體壓分布仿真,國內(nèi)外學(xué)者建立了人體和座椅有限元模型來進行研究[12-16];除了研究舒適性評價和座椅設(shè)計因素之外[17-18],還進行人體生物力學(xué)建模和體壓分布的準(zhǔn)確模擬[12]。但絕大多數(shù)研究還僅限于仿真人椅之間的相互作用;所建立的人體模型主要只針對一種身材(如50百分位身高),其原因是建模和仿真工作量非常大且復(fù)雜。對于確定的座椅和駕駛室布置參數(shù),不同身材駕駛員由于坐姿和人體尺度的差異,使體壓分布存在差別。因此,如何綜合不同身材駕駛員體壓分布對座椅舒適性進行評價和優(yōu)化很值得研究。同時,真實的姿勢是準(zhǔn)確模擬體壓分布的關(guān)鍵,也是目前許多研究欠缺之處。在輔助設(shè)計方面,少數(shù)學(xué)者通過實驗研究不同形狀、不同座椅參數(shù)(靠背和坐墊傾角、腰部支撐量等)和不同發(fā)泡材料對于體壓分布的影響,從而指導(dǎo)設(shè)計[12,18-19];乃至研究面向良好的體壓分布進行座椅設(shè)計的方法[20]。但是針對如何改進設(shè)計和輔助設(shè)計還須進行更深入的研究。
本文中建立了 H點裝置(H point machine,HPM)的有限元模型,用于通過仿真準(zhǔn)確測定座椅的布置參數(shù);建立了人體和參數(shù)化座椅有限元模型,對人體入座過程進行仿真,并分別對座椅設(shè)計時影響體壓分布的因素進行了仿真研究,得出了可用于指導(dǎo)座椅設(shè)計的規(guī)律。
1.1.1 幾何建模
通過測繪HPM實物建立幾何模型;HPM各組成部分幾何模型分別建立并組裝到一起,通過參數(shù)控制 HPM的姿勢,參見圖 1。運用 DMU(digital mockup)模塊驗證機構(gòu)的運動,尤其是腰部支撐量的變化能否體現(xiàn)在背板機構(gòu)上,尺寸的變化是否正確等。幾何模型的主要尺寸驗證如圖2所示。
軀干部機構(gòu)主要驗證腰部支撐量(L81),參見圖3。
1.1.2 有限元建模
基于幾何模型并考慮必要的幾何清理來建立有限元模型。HPM有限元模型如圖4所示。
圖1 HPM幾何模型
圖2 主要尺寸的驗證
圖3 軀干部機構(gòu)和腰部支撐量的驗證
圖4 HPM有限元模型
表1為校準(zhǔn)后的HPM模型各部分質(zhì)量,與SAE標(biāo)準(zhǔn)對比可知,HPM仿真模型的各部分質(zhì)量均在允許誤差范圍之內(nèi)。
表1 HPM仿真模型的各部分測量值
本文中利用HPM仿真方法測量駕駛室和座椅布置參數(shù)。為驗證其準(zhǔn)確性,制定了驗證方案,見表2,共25組實驗,對比測量和仿真結(jié)果。座椅設(shè)計坐墊角和靠背角分別為14°和22°,仿真結(jié)果分別為13.8°和21.8°。圖5為坐墊角為14°時改變靠背角和靠背角19°時改變坐墊角的結(jié)果??梢姡抡娼Y(jié)果和實測結(jié)果非常接近。
表2 HPM參數(shù)驗證方案
以美國ANSUR(anthropometry survey of US army)人體數(shù)據(jù)庫為數(shù)據(jù)源,利用人體切片[21]建立幾何模型,并通過縮放滿足人體尺寸要求。再根據(jù)駕駛室布置參數(shù)(表3),使用CPM(cascade prediction model)模型[22]計算和設(shè)置人體姿勢,如圖6所示。
圖5 仿真結(jié)果與實際測量結(jié)果的比較
表3 駕駛室布置參數(shù)
圖6 人體幾何模型
網(wǎng)格尺寸對仿真精度和收斂速度有重要影響。使用10mm大小的四面體單元對軟組織和骨骼建模;使用10mm大小的殼單元對皮膚建模,厚度為1mm。單元信息見表4,人體有限元模型如圖 7所示。
表4 網(wǎng)格信息
圖7 人體有限元模型
人體骨骼、皮膚采用各向同性的線彈性材料描述,利用彈性模量和泊松比控制其力學(xué)特性[23-24],具體參數(shù)如表5所示。
表5 皮膚和骨骼材料屬性
人體軟組織使用Mooney-Rivlin形式的各向同性超彈性材料進行描述,其應(yīng)變能函數(shù)為
式中:U為應(yīng)變能;C10和C01為剪切模量;D1為體積模量;I1和I2為Cauchy-Green等容形式張量的第一和第二不變量;Jel為變形前后的體積比。參數(shù)C10,C01和 D1取值為 1.65kPa,3.35kPa 和 3.653。
建立參數(shù)化的座椅幾何模型,包括坐墊、靠背、頭枕和腰托4部分,通過坐墊角、靠背角、腰托支撐量和支撐位置4個參數(shù)控制,見圖8(a)。在此基礎(chǔ)上建立有限元模型,用10mm大小的四面體單元劃分網(wǎng)格,如圖8(b)所示。
座椅發(fā)泡材料通常為聚氨酯泡沫,采用超彈性泡沫材料描述其力學(xué)特性[25],其材料特性可通過單軸壓縮實驗和剪切實驗數(shù)據(jù)來擬合[26],實驗數(shù)據(jù)如圖9所示。設(shè)置泊松比為0。
圖8 參數(shù)化座椅模型
圖9 座椅發(fā)泡材料實驗數(shù)據(jù)曲線
人體和座椅有限元模型組裝在一起,二者接觸面采用“罰函數(shù)”進行控制,以防止相互穿透。摩擦因數(shù)設(shè)置為0.3,采用“surface-to-surface”的形式定義接觸。軟組織和骨骼間的接觸面通過Tie連接定義。約束座椅發(fā)泡背面單元和手部與腳部單元各6個自由度,并施加重力場。
為方便分析,在座椅上定義了局部坐標(biāo)系和關(guān)鍵截面,見圖10和表6,坐標(biāo)原點與H點重合。
圖10 局部坐標(biāo)系和關(guān)鍵截面
建立了5,50和95 3個百分位男性駕駛員人體模型,并進行體壓分布仿真。圖11為斷面Ⅱ內(nèi)的人體軟組織應(yīng)力分布。圖12為座椅上的壓強分布云圖。可見,軟組織中最大應(yīng)力都出現(xiàn)在坐骨結(jié)節(jié)下方,說明坐骨結(jié)節(jié)處承受的壓強較大。坐墊上壓強以坐骨結(jié)節(jié)下方為中心向四周逐漸減小。從坐骨結(jié)節(jié)下方至大腿前下方,壓強也逐漸減小??勘成献畲髩簭姵霈F(xiàn)的位置因人而異,這是由于不同人體身材差異造成的。厚度建議為40mm左右,靠背其余部分厚度建議為30mm左右。
表6 4個關(guān)鍵截面的定義
圖11 人體軟組織應(yīng)力
圖12 座椅表面的壓強分布
圖13 不同海綿厚度座椅的壓強分布
圖14 座椅接觸面積和最大壓強隨厚度的變化趨勢
坐墊和靠背海綿發(fā)泡厚度從20變到70mm,步長10mm。限于篇幅,僅列出50百分位模型落座時的壓強分布(其他人的規(guī)律相同),見圖13??梢?,海綿發(fā)泡厚度的增加會增大人體與座椅的接觸面積,有利于減小最大壓強。圖14為各厚度時的接觸面積和最大壓強,同樣反映了上述規(guī)律。從圖15可見,軟組織最大應(yīng)力隨發(fā)泡厚度的增加明顯減小,大腿與坐墊前緣接觸區(qū)域的壓強和人體坐骨結(jié)節(jié)下方的軟組織應(yīng)力也明顯減小??紤]到成本等原因,坐墊臀部下方的發(fā)泡厚度建議為60mm左右,大腿前下方發(fā)泡厚度建議為30mm左右;靠背腰托處發(fā)泡
圖15 軟組織最大應(yīng)力隨厚度的變化趨勢
以座椅初始姿態(tài)下的腰部支撐位置為基礎(chǔ)進行腰部支撐量仿真,其中座椅初始姿態(tài)下的腰部支撐量和支撐位置都為0。圖16為腰托突出量的示意圖。以此分析了腰托凸出量0~40mm(步長10mm)時,腰托位置相對于初始位置從-20變化到20mm(步長10mm,由初始位置沿靠背線向下為負)時靠背的體壓分布。隨著腰部支撐量的增加靠背的接觸面積明顯減小,見圖17(b);而最大壓強則增大,見圖17(d);而坐墊接觸面積和最大壓強與腰部支撐量關(guān)系不大。肩胛骨處的最大壓強減小,而腰部最大壓強則增大。經(jīng)驗表明,適當(dāng)提高腰部的最大壓強有利于提高舒適性和保持舒適的坐姿。
圖16 腰托突出量示意圖
圖17 接觸面積和最大壓強隨腰部支撐量的變化趨勢
由于身材尺寸的差異,不同駕駛員對腰部支撐量的需求不同。從仿真結(jié)果看,腰部支撐量為20mm左右時3位駕駛員腰部能得到較為合理的支撐。圖18為腰部支撐量為20mm時人體與座椅的接觸面積和峰值壓強隨腰部支撐位置變化的趨勢。可見,接觸面積基本不受腰部支撐位置變化的影響,而靠背峰值壓強對于不同駕駛員呈現(xiàn)不規(guī)律的變化。
圖18 接觸面積和最大壓強隨支撐位置的變化趨勢
圖19 為單獨改變靠背角(坐墊角為15°,靠背角分別取20°,23°和 26°)時體壓分布的變化,每次坐墊角或靠背角的仿真調(diào)整都是通過HPM測量實現(xiàn)的。圖20給出了斷面Ⅲ和斷面Ⅰ處體壓分布隨靠背角增加的變化趨勢。由圖可見,隨著靠背角的增加,靠背腰部附近壓強增加,且越靠近下方增加幅度越大;上背部的最大壓強峰值降低。對于坐墊,雖然壓強水平在各處沒有明顯變化,但峰值位置略向后方移動。
圖19 不同靠背角時的壓強分布
圖20 斷面Ⅲ和斷面Ⅰ處的壓強分布
圖21 為單獨改變坐墊角(靠背角為23°,坐墊角分別取12°,15°和18°)時體壓分布的變化。圖22給出了斷面Ⅲ和斷面Ⅰ處體壓分布隨坐墊角增加的變化趨勢。由圖可見,坐墊角對靠背體壓分布幾乎沒影響;而坐墊角增大時,坐墊上大腿前下方的壓強明顯增加。
圖21 不同坐墊角時的壓強分布
圖22 斷面Ⅲ和斷面Ⅰ處的壓強分布
結(jié)果表明,靠背角過小不能有效地支撐人體腰部;靠背角過大時雖然能充分支撐人體,但易使人體過度后傾,容易造成頸部緊張。坐墊角過大會增加大腿下方與坐墊前緣接觸處的壓強,易造成局部血液不暢;坐墊角過小則易導(dǎo)致臀部下方壓強過大,產(chǎn)生不舒適的感覺。不同身材駕駛員對座椅靠背和坐墊角的需求不同,二者必須相互協(xié)調(diào)才能提高座椅舒適性。
(1)建立了H點裝置有限元模型,實現(xiàn)對座椅坐墊角、靠背角、H點位置和腰部支撐量的仿真測量方法。建立了人體和座椅有限元模型,對人椅間體壓分布進行了仿真研究。
(2)研究了座椅海綿發(fā)泡厚度、腰部支撐量與支撐位置、坐墊角和靠背角對體壓分布的影響。結(jié)果表明,適當(dāng)增加發(fā)泡厚度有利于降低人椅接觸面最大壓強,但要綜合考慮成本等因素確定。適當(dāng)增加腰部支撐能改善舒適性,支撐量為20mm時3位駕駛員舒適性最好。支撐位置應(yīng)兼顧不同身材駕駛員。坐墊角和靠背角的最佳組合,取決于不同身材駕駛員對這兩個角度的要求。
(3)本文中的仿真方法已經(jīng)經(jīng)過實驗驗證[12],非常適用于在座椅設(shè)計早期對駕駛員體壓分布進行評估,和面向理想的體壓分布來確定座椅的設(shè)計參數(shù)。但由于目前的仿真還存在一定的簡化,今后將考慮更多的實際因素(座椅面料、發(fā)泡支撐彈簧和著裝因素等),使仿真更接近實際情況,確保得出的規(guī)律具有更實用的價值。