張紹偉，朱大煒，翟光照

（1.法國(guó)ESI集團(tuán)上海分公司，上海 200000；2.上汽大眾汽車有限公司產(chǎn)品研發(fā)車身研發(fā)部，上海 200000）

前言

在汽車追尾碰撞事故中，78.2%的損傷發(fā)生在頸部［1］，稱之為揮鞭傷。揮鞭傷易造成頸部的骨骼或軟組織損傷，這種損傷并發(fā)癥多、潛伏期長(zhǎng)，有些甚至是不可治愈的永久傷害。正確合理的座椅設(shè)計(jì)直接影響到追尾碰撞中乘員的頸部保護(hù)效果［2-3］。目前C-NCAP、IIHS、C-IASI、Euro-NCAP等均采用BioRIDII假人來(lái)考察乘員頸部傷害。但BioRIDII假人的物理結(jié)構(gòu)導(dǎo)致鞭打?qū)嶒?yàn)重復(fù)性差。目前，現(xiàn)有的仿真方法多用于分析座椅各參數(shù)對(duì)座椅抗揮鞭傷性能（下簡(jiǎn)稱揮鞭性能或鞭打性能）的影響趨勢(shì)分析［4］，并不能在座椅設(shè)計(jì)初期對(duì)揮鞭性能有較好的預(yù)測(cè)。所以，設(shè)計(jì)一種穩(wěn)健、良好的驗(yàn)證方法來(lái)保證仿真模型的精確性就顯得尤為重要。

人工智能（artificial intelligence，AI）是計(jì)算機(jī)科學(xué)領(lǐng)域、涉及自動(dòng)化智能行為和機(jī)器學(xué)習(xí)的方法。有限元分析（FEA）是使用數(shù)值方法對(duì)給定的物理現(xiàn)象進(jìn)行仿真，以獲得由偏微分方程控制的一組問(wèn)題的近似解決方案。在機(jī)器學(xué)習(xí)領(lǐng)域，人工神經(jīng)網(wǎng)絡(luò)用于映射某些輸入和輸出變量之間的相互關(guān)系和依存關(guān)系。本文中利用CAE仿真結(jié)合深度學(xué)習(xí)算法，建立了一種長(zhǎng)短期記憶（long short-term memory，LSTM）網(wǎng)絡(luò)模型來(lái)預(yù)測(cè)假人響應(yīng)曲線，通過(guò)該方法可以快速獲得改變變量后假人的響應(yīng)曲線，為快速查看假人傷害響應(yīng)提供了一種有效的方法。

1 C-NCAP鞭打?qū)嶒?yàn)評(píng)分規(guī)則

C-NCAP2018規(guī)則對(duì)于鞭打試驗(yàn)的考核只針對(duì)前排駕駛員座椅，總分5分［5］；而C-NCAP2021規(guī)則增加了后排座椅的鞭打性能考核［6］，其評(píng)價(jià)指標(biāo)與駕駛員座椅鞭打試驗(yàn)相同，將結(jié)果乘以0.4，作為第二排座椅鞭打試驗(yàn)最終得分。新規(guī)則中，鞭打試驗(yàn)總成績(jī)?yōu)?分（駕駛員座椅5分，第二排座椅2分），在總分5分制的規(guī)則下，假人各部分的傷害值評(píng)價(jià)細(xì)則如表1所示［6］。

表1 鞭打?qū)嶒?yàn)總體評(píng)分準(zhǔn)則

其中，NIC（neck injury criteria）是評(píng)價(jià)揮鞭傷的重要指標(biāo)。NIC的計(jì)算方法綜合考慮了假人頭部加速度和胸部T1加速度，NIC的計(jì)算公式［5-6］為

其中

NIC在一定程度上反映了座椅頭枕與座椅靠背向后移動(dòng)的同步運(yùn)動(dòng)程度，座椅頭枕與座椅靠背的同步程度越高，對(duì)假人的頸部拉伸作用就越小，假人上頸部與下頸部的受力也就越小，反之亦然。

假人頭部X向加速度與胸部T1加速度的時(shí)間歷程如圖1所示，提高NIC可以通過(guò)降低圖1中所示的灰色區(qū)域的面積來(lái)實(shí)現(xiàn)，主要的實(shí)現(xiàn)方法有兩種。

（1）將假人與頭枕的接觸時(shí)間提前，主要依靠?jī)煞N方式來(lái)實(shí)現(xiàn)：

①在設(shè)計(jì)上，盡可能減小假人頭部與座椅頭枕間隙（backset）；

②采用主動(dòng)式頭枕，在發(fā)生碰撞初始，通過(guò)頭枕主動(dòng)式前移，減少假人頭部與座椅頭枕間隙。

（2）將假人胸部T1椎骨與座椅靠背的接觸時(shí)間后移，主要依靠?jī)煞N方式來(lái)實(shí)現(xiàn)：

①大幅減少靠背剛度，使假人能較快陷入座椅之中；

②采用可活動(dòng)式靠背［7-8］，在假人向后運(yùn)動(dòng)的同時(shí)，座椅靠背在保持頭枕位置相對(duì)不變的情況下，下沉或后陷。

一般情況下，提高鞭打成績(jī)的優(yōu)化方案與座椅舒適性、頭枕靜態(tài)打擊等相關(guān)要求相矛盾，為能夠在座椅設(shè)計(jì)初期階段就考慮揮鞭傷性能，須對(duì)可能影響揮鞭傷性能的所有因素進(jìn)行考慮，但多因子分析須進(jìn)行多輪實(shí)驗(yàn)，整個(gè)過(guò)程周期長(zhǎng)、費(fèi)用高。因此，如何在設(shè)計(jì)初期，尚無(wú)物理實(shí)驗(yàn)條件的情況下快速而精準(zhǔn)地預(yù)測(cè)鞭打性能，顯得尤為重要。

2 座椅零部件物理實(shí)驗(yàn)與仿真對(duì)標(biāo)

為能夠在設(shè)計(jì)初期就了解座椅的鞭打性能，傳統(tǒng)方法是對(duì)鞭打現(xiàn)象進(jìn)行有限元仿真，指導(dǎo)設(shè)計(jì)開(kāi)發(fā)。后期通過(guò)物理實(shí)驗(yàn)對(duì)基于仿真的設(shè)計(jì)進(jìn)行驗(yàn)證。在此過(guò)程中，核心工作就是保證仿真響應(yīng)與物理實(shí)驗(yàn)響應(yīng)盡量一致。

2.1 座椅零部件實(shí)驗(yàn)

座椅鞭打是一個(gè)具有高度敏感性的非線性系統(tǒng)，座椅骨架、發(fā)泡、面套、頭枕?xiàng)U、頭枕等多個(gè)設(shè)計(jì)因素對(duì)鞭打結(jié)果有直接影響，因此須對(duì)整個(gè)座椅模型進(jìn)行精細(xì)化建模，并通過(guò)一系列的材料級(jí)、零部件級(jí)、分總成級(jí)對(duì)座椅仿真模型進(jìn)行標(biāo)定。在確保仿真模型精度的情況下，采用標(biāo)定后的仿真模型對(duì)影響鞭打性能的因素進(jìn)行分析，找出這些因素對(duì)汽車座椅揮鞭傷的影響程度，以便在座椅開(kāi)發(fā)過(guò)程中對(duì)這些因素進(jìn)行優(yōu)化［3］。

本文中以上汽大眾某平臺(tái)座椅作為研究對(duì)象，對(duì)人與座椅相互作用的重點(diǎn)接觸區(qū)域，例如頭枕與靠背，取樣同一批產(chǎn)線樣件，分別進(jìn)行材料級(jí)、零部件級(jí)、分總成級(jí)和整椅級(jí)的一系列實(shí)驗(yàn)，部分實(shí)驗(yàn)設(shè)置如圖2所示。

2.2 曲線相似度評(píng)價(jià)方法

采用Cora等級(jí)方法評(píng)價(jià)兩條曲線之間的區(qū)別。Cora等級(jí)方法是一個(gè)計(jì)算兩個(gè)非模糊信號(hào)（如時(shí)間歷程曲線）之間相關(guān)性水平的軟件［9］。利用在目標(biāo)曲線以外生成兩個(gè)通道并形成3個(gè)區(qū)域的方法來(lái)計(jì)算偏差，并將評(píng)價(jià)結(jié)果歸一化至"0"和"1"之間，具體取決于用戶自定義的匹配質(zhì)量。其中，"1"代表評(píng)價(jià)曲線完全位于內(nèi)通道內(nèi)，即在定義的容差范圍內(nèi)100%匹配，"0"代表曲線位于外通道以外，即在定義的容差范圍內(nèi)0%匹配，中間部分線性插值［9］，如圖3所示。

2.3 座椅零部件仿真對(duì)標(biāo)

為保證實(shí)驗(yàn)的一致性，并驗(yàn)證實(shí)驗(yàn)的重復(fù)性，材料級(jí)實(shí)驗(yàn)分別重復(fù)5次，分總成級(jí)至少重復(fù)3次。同時(shí)，分別建立了相同實(shí)驗(yàn)設(shè)置的仿真模型，結(jié)合實(shí)驗(yàn)結(jié)果，對(duì)仿真模型進(jìn)行標(biāo)定，材料與零部件實(shí)驗(yàn)的矩陣及其仿真與實(shí)驗(yàn)的Cora等級(jí)對(duì)比如表2所示。

表2 材料、零部件實(shí)驗(yàn)矩陣及其對(duì)標(biāo)情況

其中分總成級(jí)實(shí)驗(yàn)中的座椅整體骨架剛度、腰托剛度、頭枕動(dòng)態(tài)打擊的仿真與實(shí)驗(yàn)結(jié)果對(duì)比如圖4所示。整體Cora等級(jí)大于98%，由此可以看出，仿真模型能夠較好地表征物理樣件的真實(shí)性能。

3 整椅鞭打?qū)嶒?yàn)與仿真標(biāo)定

鞭打系統(tǒng)仿真預(yù)測(cè)工作涉及到有限元假人的精度、座椅模型的精度和假人與座椅的交互響應(yīng)，傳統(tǒng)的實(shí)驗(yàn)與仿真對(duì)標(biāo)方法是將仿真假人與座椅的位置調(diào)整到與實(shí)驗(yàn)設(shè)置一致的情況，并對(duì)座椅的相關(guān)材料屬性等進(jìn)行微調(diào)來(lái)達(dá)到實(shí)驗(yàn)與仿真假人的動(dòng)態(tài)響應(yīng)趨勢(shì)吻合［10］，這種對(duì)標(biāo)方式是基于物理假人與仿真假人100%一致的情況下進(jìn)行的，忽略了物理假人與仿真假人本身的不同。本文中基于第2節(jié)對(duì)標(biāo)良好的整椅模型與上汽大眾鞭打仿真“最佳實(shí)踐”方法，建立了如圖5所示的鞭打仿真模型進(jìn)行整椅級(jí)別的仿真對(duì)標(biāo)。

對(duì)標(biāo)完成后的實(shí)驗(yàn)與仿真結(jié)果一致性非常高，假人傷害相關(guān)通道的實(shí)驗(yàn)與仿真對(duì)比Cora等級(jí)評(píng)價(jià)如圖6所示。

假人頭部加速度曲線、胸部T1平均加速度曲線和NIC曲線如圖7所示。由圖可見(jiàn)，仿真預(yù)測(cè)結(jié)果與實(shí)驗(yàn)結(jié)果吻合良好，其中頭部加速度曲線的Cora等級(jí)達(dá)到了0.974。

至此，已經(jīng)得到一個(gè)仿真精度良好的整椅計(jì)算模型，將作為下一步深度學(xué)習(xí)生成數(shù)據(jù)的基礎(chǔ)。根據(jù)深度學(xué)習(xí)方法的特點(diǎn)，必須選取一些對(duì)于揮鞭傷結(jié)果有重要影響的因素，同時(shí)為使數(shù)據(jù)集生成過(guò)程更高效，也有必要對(duì)此步驟進(jìn)行流程化開(kāi)發(fā)。

4 獲取深度學(xué)習(xí)所需的關(guān)鍵數(shù)據(jù)

4.1 影響鞭打的關(guān)鍵因素

如前文所述，影響座椅鞭打成績(jī)的因素有很多，陳長(zhǎng)亮等［3］對(duì)靠背剛度、頭枕剛度等進(jìn)行了研究。任海波等［11］研究了座椅面套與假人的摩擦因數(shù)對(duì)座椅防揮鞭性能的影響。李永攀等［12］對(duì)面套鋼絲、防異響無(wú)紡布、靠背骨架彈簧等進(jìn)行研究。周澤斌等［13］通過(guò)改進(jìn)頭枕的結(jié)構(gòu)并優(yōu)化頭枕強(qiáng)度和發(fā)泡厚度等措施，來(lái)提高座椅的防揮鞭傷性能。王晉等［14］對(duì)假人的坐姿、座椅的試驗(yàn)狀態(tài)（坐墊高度調(diào)節(jié)裝置）和座椅發(fā)泡蒙皮的光滑程度等進(jìn)行研究。寧成業(yè)等［15］通過(guò)調(diào)整座椅頭枕頭后間隙和座椅靠背骨架焊接的方式提升座椅整體剛度，改進(jìn)了鞭打性能。

雖然國(guó)內(nèi)外學(xué)者對(duì)座椅鞭打做了大量研究，但都基于某個(gè)或某幾個(gè)關(guān)鍵因素進(jìn)行，整椅的鞭打性能與各因素相關(guān)性很大，所以有必要嘗試同時(shí)對(duì)所有的關(guān)鍵因素進(jìn)行研究。本文中基于上汽大眾某平臺(tái)座椅項(xiàng)目經(jīng)驗(yàn)與各學(xué)者的研究，盡可能考慮所有影響因素，探尋各影響因素之間的相關(guān)性以及與假人傷害值之間的關(guān)系，并將所得到的對(duì)應(yīng)關(guān)系作為下一步深度學(xué)習(xí)的基礎(chǔ)數(shù)據(jù)，具體的影響因素如表3所示。

表3 影響鞭打的關(guān)鍵因素

4.2 分析流程自動(dòng)化

考慮到同時(shí)計(jì)算所有影響因素，數(shù)據(jù)樣本非常龐大，數(shù)據(jù)準(zhǔn)備過(guò)程與處理生成過(guò)程冗繁，因此，須對(duì)整個(gè)數(shù)據(jù)準(zhǔn)備過(guò)程自動(dòng)化，把標(biāo)定良好的鞭打模型作為基礎(chǔ)模型，并將所有的變量通過(guò)PAMCRASH的PYVAR關(guān)鍵字設(shè)置為變量，賦予每個(gè)變量一定的變動(dòng)范圍，利用Python建立自動(dòng)化腳本實(shí)現(xiàn)自動(dòng)化因子更改、自動(dòng)化和HPC（高性能計(jì)算機(jī)）交互。另外有：

（1）假人位置的改變會(huì)導(dǎo)致假人與座椅泡沫之間產(chǎn)生干涉，同時(shí)引起假人與安全帶不匹配，因此須對(duì)座椅泡沫重新壓縮，安全帶重新進(jìn)行纏繞。這整個(gè)過(guò)程也需要自動(dòng)化，因此開(kāi)發(fā)了針對(duì)假人定位、座椅調(diào)節(jié)、座椅泡沫自動(dòng)壓縮、替換、安全帶自動(dòng)纏繞的程序，并應(yīng)用于本項(xiàng)研究，程序界面如圖8所示。

（2）傳統(tǒng)的商用后處理軟件一般僅能同時(shí)承載5～10個(gè)大型計(jì)算結(jié)果文件，因此整個(gè)計(jì)算結(jié)果的后處理應(yīng)不依賴于商用后處理軟件。為此，利用Python重新編程，對(duì)PAMCRASH計(jì)算后的二進(jìn)制結(jié)果文件.erfh5進(jìn)行快速讀取、分類、統(tǒng)計(jì)和呈現(xiàn)。最終，實(shí)現(xiàn)了大量計(jì)算結(jié)果的快速自動(dòng)化處理，極大提高了效率。整個(gè)自動(dòng)化流程如圖9所示。

5 深度學(xué)習(xí)工作流程

5.1 訓(xùn)練數(shù)據(jù)準(zhǔn)備

將表3中影響揮鞭傷的關(guān)鍵因素A～M按照其允許的最大波動(dòng)范圍建立全因子輸入?yún)?shù)集，并利用碰撞分析有限元軟件Pamcrash進(jìn)行仿真，共完成15 000次仿真，獲得了輸入數(shù)據(jù)集。

按照機(jī)器學(xué)習(xí)的術(shù)語(yǔ)，將影響揮鞭傷的關(guān)鍵因素A～M稱之為特征值，由于各特征值之間存在量綱差異，為消除這種單位以及尺度差異等的影響，需要對(duì)每個(gè)特征值進(jìn)行歸一化［16-17］。

本文中采用標(biāo)準(zhǔn)差方法作為歸一化方法：

式中xˉ和σ分別為對(duì)應(yīng)傷害值的均值和標(biāo)準(zhǔn)差。歸一化后的特征將分布在［-1，1］區(qū)間且服從正態(tài)分布。

雖然每個(gè)變量都提供了一定信息，但變量太多增加了計(jì)算的復(fù)雜性，也給合理性分析和解釋帶來(lái)困難，每個(gè)變量所提供的信息在一定程度上有所重疊，因此須設(shè)法將原來(lái)變量重新組合成一組相互無(wú)關(guān)的新綜合變量，進(jìn)行維度降低，同時(shí)根據(jù)實(shí)際需要，從中選出盡量少的變量盡可能多地表達(dá)原來(lái)變量的信息，本文中利用主成分分析法進(jìn)行數(shù)據(jù)降維，將降維后的數(shù)據(jù)作為全新的數(shù)據(jù)樣本，用于深度學(xué)習(xí)的模型訓(xùn)練，整個(gè)方法流程如圖10所示。

5.2 模型訓(xùn)練

BioRIDII假人傷害響應(yīng)是假人向后運(yùn)動(dòng)過(guò)程中與座椅接觸，假人頭部與座椅頭枕碰撞以及座椅變形過(guò)程中（0.1 s內(nèi)）采集到的假人體內(nèi)若干部位的傳感器信號(hào)，其響應(yīng)曲線是時(shí)間序列的增量。較為傳統(tǒng)的方式是可以使用前饋神經(jīng)網(wǎng)絡(luò)（feedforward neural network，F(xiàn)FNN）來(lái)表示整個(gè)過(guò)程假人的響應(yīng)，但這個(gè)過(guò)程需要大量的FFNN以及一些不確定變量，另外，傳統(tǒng)的FFNN方法無(wú)法表征響應(yīng)對(duì)于時(shí)間的依賴性［17-18］。

循環(huán)神經(jīng)網(wǎng)絡(luò)（recurrent neural network，RNN）是一種可以預(yù)測(cè)數(shù)據(jù)序列（如時(shí)間序列）的體系結(jié)構(gòu)。與前饋神經(jīng)網(wǎng)絡(luò)不同，RNN使用序列中先前評(píng)估的信息預(yù)測(cè)下一個(gè)值。本文采用一種特殊的RNN結(jié)構(gòu)，即LSTM網(wǎng)絡(luò)。圖11示出該神經(jīng)網(wǎng)絡(luò)的主要架構(gòu)。通過(guò)LSTM與FFNN的結(jié)合來(lái)預(yù)測(cè)假人響應(yīng)，同時(shí)解決整個(gè)響應(yīng)過(guò)程中對(duì)時(shí)間的依賴問(wèn)題［19-21］。模型由開(kāi)源的Python深度學(xué)習(xí)框架PyTorch進(jìn)行訓(xùn)練［22］。

輸入數(shù)據(jù)集被隨機(jī)分為兩部分，75%作為訓(xùn)練數(shù)據(jù)集，25%作為測(cè)試數(shù)據(jù)集，用于評(píng)估“訓(xùn)練后的”人工神經(jīng)網(wǎng)絡(luò)的泛化能力。

模型采用Adam優(yōu)化器［22］，初始學(xué)習(xí)率為0.001，激活函數(shù)采用tanh函數(shù)，具體的模型架構(gòu)如表4所示。

表4 LSTM-FFNN模型架構(gòu)

同時(shí)，使用均方誤差（MSE）函數(shù)度量模型性能：

當(dāng)ypredict=ytest時(shí)，誤差降為0，同時(shí)計(jì)算每條曲線的均方誤差（MSE），訓(xùn)練集與測(cè)試集的均方誤差變化過(guò)程如圖12所示。

5.3 模型預(yù)測(cè)

為進(jìn)一步驗(yàn)證訓(xùn)練后的神經(jīng)網(wǎng)絡(luò)的泛化能力，重新生成了30個(gè)全新數(shù)據(jù)集，新數(shù)據(jù)集與原訓(xùn)練數(shù)據(jù)集及測(cè)試數(shù)據(jù)集無(wú)交叉。為此自由組合不同的變種，同時(shí)利用PAMCRASH計(jì)算出每個(gè)變種的假人傷害值曲線，并利用訓(xùn)練后的神經(jīng)網(wǎng)絡(luò)模型對(duì)假人的傷害曲線進(jìn)行預(yù)測(cè)，再對(duì)經(jīng)過(guò)計(jì)算的假人傷害值曲線與預(yù)測(cè)的假人傷害值曲線進(jìn)行對(duì)比，其中一組特征值的對(duì)比結(jié)果如圖13所示。預(yù)測(cè)曲線在很大程度上表征了仿真曲線的峰值與脈寬，能夠較好地反映假人響應(yīng)。

6 總結(jié)與展望

汽車座椅是一個(gè)涉及金屬和多種軟質(zhì)材料的復(fù)雜系統(tǒng)，尤其在揮鞭傷工況中，各種影響因素交織在一起形成了最終的假人傷害響應(yīng)?；趥鹘y(tǒng)的有限元法，如未分層次對(duì)其中的各因素進(jìn)行對(duì)標(biāo)，仿真預(yù)測(cè)效果會(huì)非常一般。本文提出了一系列材料級(jí)、零部件級(jí)、子系統(tǒng)級(jí)和整椅級(jí)的實(shí)驗(yàn)方法，通過(guò)這些實(shí)驗(yàn)與仿真的對(duì)標(biāo)，獲得仿真精度較高的計(jì)算模型。然后，從眾多影響因素中選取若干關(guān)鍵因素，通過(guò)自動(dòng)化的CAE分析流程獲得大量數(shù)據(jù)集，借助于深度學(xué)習(xí)技術(shù)長(zhǎng)短期記憶網(wǎng)絡(luò)構(gòu)建了一種可以快速預(yù)測(cè)響應(yīng)值，即假人傷害曲線的神經(jīng)網(wǎng)絡(luò)架構(gòu)，可以快速地預(yù)測(cè)假人在不同設(shè)計(jì)變量下的響應(yīng)，使工程師在更改設(shè)計(jì)變量的同時(shí)，即刻預(yù)測(cè)假人響應(yīng)。由于不必再進(jìn)行時(shí)間冗長(zhǎng)的CAE計(jì)算，所以這種快速獲得響應(yīng)的優(yōu)勢(shì)可以幫助工程技術(shù)人員在較短的時(shí)間內(nèi)獲得最優(yōu)解，有效地提高了開(kāi)發(fā)效率。

目前該方法還在進(jìn)一步拓展之中，在預(yù)測(cè)假人時(shí)間歷程響應(yīng)的同時(shí)，同步獲得座椅與假人隨時(shí)間變化的變形情況，將是后續(xù)研究的重點(diǎn)之一。