廖卓　趙橄培　羅龍健

摘 要：仿真模型是人們開展各類科學(xué)研究的重要工具，其自身的可信度將直接影響仿真運(yùn)行結(jié)果的可靠性，而模型驗(yàn)證則是評(píng)價(jià)仿真模型可信度的主要方法。基于此，本文綜述了汽車動(dòng)態(tài)系統(tǒng)仿真模型相關(guān)概念及模型驗(yàn)證方法發(fā)展歷程，根據(jù)汽車安全領(lǐng)域仿真模型的特點(diǎn)，重點(diǎn)對(duì)時(shí)域分析方向已有的模型驗(yàn)證方法進(jìn)行了分類，并對(duì)相應(yīng)的模型驗(yàn)證輔助工具存在的問題進(jìn)行了梳理，提出了未來的重點(diǎn)研究方向。

關(guān)鍵詞：汽車安全 仿真技術(shù) 動(dòng)態(tài)系統(tǒng) 模型驗(yàn)證

1 前言

數(shù)字化建模與仿真技術(shù)已完成為人們開展各類研究的重要手段，尤其在汽車被動(dòng)安全研究領(lǐng)域的地位更為突出，正因?yàn)榉抡婕夹g(shù)的重要作用，仿真模型的可信度必然備受用戶關(guān)注，模型的可信度將直接影響仿真結(jié)果的準(zhǔn)確性和實(shí)用價(jià)值，因此，針對(duì)仿真模型，尤其是復(fù)雜動(dòng)態(tài)系統(tǒng)仿真模型的驗(yàn)證方法及相關(guān)輔助工具的研究逐漸成為國內(nèi)外研究熱點(diǎn)。本文旨在綜述汽車動(dòng)態(tài)系統(tǒng)仿真模型驗(yàn)證方法的研究進(jìn)展，梳理現(xiàn)階段模型驗(yàn)證工具存在的問題，指出其未來發(fā)展方向。

2 汽車動(dòng)態(tài)系統(tǒng)仿真模型及模型驗(yàn)證

2.1 汽車動(dòng)態(tài)系統(tǒng)仿真模型

數(shù)字化仿真已經(jīng)成為新車開發(fā)過程中必不可少的一部分，更是汽車被動(dòng)安全研究領(lǐng)域的主要方式之一。汽車被動(dòng)安全領(lǐng)域的汽車碰撞過程（如圖1所示），是一個(gè)具有高度非線性和復(fù)雜性的動(dòng)態(tài)過程，該過程中包含的車、人和物（道路安全設(shè)施）等相關(guān)環(huán)境共同組成汽車動(dòng)態(tài)系統(tǒng)。在針對(duì)此類汽車動(dòng)態(tài)系統(tǒng)進(jìn)行碰撞研究分析時(shí)，一方面需要關(guān)注動(dòng)態(tài)系統(tǒng)輸出響應(yīng)數(shù)據(jù)之間高度的時(shí)間相關(guān)性和時(shí)變性，同時(shí)需要將其相位、幅度及頻率等曲線特征差異也納入監(jiān)控范圍，另一方面動(dòng)態(tài)系統(tǒng)通常包含多個(gè)輸出響應(yīng)，且響應(yīng)間存在著很強(qiáng)的相關(guān)性，因而需要建立可信的汽車動(dòng)態(tài)系統(tǒng)仿真模型進(jìn)行研究分析。

2.2 仿真模型驗(yàn)證方法發(fā)展歷程

仿真模型本身的可信度會(huì)直接影響應(yīng)用該模型研究所得結(jié)果的準(zhǔn)確性，尤其是汽車動(dòng)態(tài)系統(tǒng)等具有復(fù)雜運(yùn)行機(jī)理和交互關(guān)系的仿真模型，其內(nèi)部缺陷和錯(cuò)誤不易被發(fā)現(xiàn)，若不對(duì)其進(jìn)行驗(yàn)證，將會(huì)給模型應(yīng)用帶來巨大的風(fēng)險(xiǎn)。模型驗(yàn)證技術(shù)和仿真技術(shù)的發(fā)展密不可分，總的來看，仿真技術(shù)經(jīng)歷了探索、發(fā)展、成熟和高級(jí)等四個(gè)階段。在仿真技術(shù)應(yīng)用需求的驅(qū)動(dòng)下，模型驗(yàn)證技術(shù)也從最初的概念提出發(fā)展到現(xiàn)代的復(fù)雜仿真模型驗(yàn)證階段，取得了高速的發(fā)展。

在上世紀(jì)50年代末期，國外學(xué)者開始對(duì)模型驗(yàn)證相關(guān)理論和方法進(jìn)行研究并取得大量研究成果，其具體發(fā)展歷程如圖2所示。仿真模型的可信度問題由Conway等研究人員于1959年首次提出[1]；Fishman等則在1967年闡述了模型驗(yàn)證的相關(guān)概念[2]；進(jìn)入20世紀(jì)70年代，發(fā)展技術(shù)逐漸成熟，相應(yīng)的模型驗(yàn)證方法也步入初步發(fā)展階段，在這一時(shí)期，相關(guān)研究人員提出了一些定性和定量模型驗(yàn)證分析方法，例如圖示比較法、Turing測試法、頻譜分析法、假設(shè)檢驗(yàn)法、Theil不等式系數(shù)法等，但這些方法主要應(yīng)用于簡單仿真模型的驗(yàn)證分析，其研究重點(diǎn)主要集中在模型驗(yàn)證的基本理論和分析方法。從1990年開始，仿真技術(shù)進(jìn)入高速發(fā)展階段，仿真模型的復(fù)雜程度大幅提升，為滿足復(fù)雜仿真模型的驗(yàn)證需求，相應(yīng)的仿真模型驗(yàn)證方法和相關(guān)技術(shù)也得到進(jìn)一步的深入研究。Balci等根據(jù)仿真模型的特點(diǎn)，將仿真模型驗(yàn)證分為動(dòng)態(tài)模型驗(yàn)證和靜態(tài)模型驗(yàn)證兩類[3]。針對(duì)復(fù)雜仿真模型驗(yàn)證過程中面臨的影響因素不確定性、交互關(guān)系復(fù)雜性、輸出變量多元性等問題，相關(guān)學(xué)者分別提出了貝葉斯網(wǎng)絡(luò)、多元概率積分變換、概率包絡(luò)等驗(yàn)證方法。此外，Wu等針對(duì)復(fù)雜仿真模型的多元輸出等問題，提出了基于多元輸出面積差法的模型驗(yàn)證方法；Claudia等則提出了一種可重用仿真模型的語義驗(yàn)證方法；Lizcano等則提出通過Web終端軟件來進(jìn)行仿真模型的自動(dòng)校核與驗(yàn)證等。

國內(nèi)學(xué)者對(duì)仿真模型驗(yàn)證的相關(guān)研究起步相對(duì)較晚，部分學(xué)者于上世紀(jì)80年代末，開始對(duì)仿真模型驗(yàn)證相關(guān)理論及方法開展研究，并對(duì)國外仿真模型的可信度評(píng)估及VV&A（Verification， Validation and Accreditation）的相關(guān)理論和方法開展了綜述分析。另有部分學(xué)者則對(duì)應(yīng)用于仿真實(shí)驗(yàn)結(jié)果驗(yàn)證的頻譜分析法進(jìn)行了進(jìn)一步研究。在國外已有相關(guān)研究的基礎(chǔ)上，有學(xué)者在針對(duì)導(dǎo)彈系統(tǒng)仿真模型驗(yàn)證的相關(guān)研究中，提出了灰色關(guān)聯(lián)分析法。進(jìn)入21世紀(jì)后，則重點(diǎn)針對(duì)復(fù)雜仿真模型的驗(yàn)證方法開展了相關(guān)研究工作，研究了基于動(dòng)態(tài)系統(tǒng)仿真模型的多元響應(yīng)分析方法MRA（Multivariate Responses Analysis），基于數(shù)據(jù)特征的評(píng)估方法、考慮不確定性的多元?jiǎng)討B(tài)輸出結(jié)果驗(yàn)證方法、時(shí)域分析和頻域分析相結(jié)合的可信度評(píng)估方法等。

3 動(dòng)態(tài)系統(tǒng)仿真模型驗(yàn)證方法及工具

3.1 動(dòng)態(tài)系統(tǒng)仿真模型驗(yàn)證方法

通過綜述仿真模型驗(yàn)證方法發(fā)現(xiàn)，總體可將其分為定性分析、定量分析和綜合分析三類，定性分析方法主要包括：Turing測試、圖示比較、區(qū)間估值等。根據(jù)具體分析對(duì)象特點(diǎn)，定量分析方法又可分為靜態(tài)系統(tǒng)分析方法和動(dòng)態(tài)系統(tǒng)分析方法兩種類型，而動(dòng)態(tài)系統(tǒng)的輸出數(shù)據(jù)種包含了時(shí)域數(shù)據(jù)和頻域數(shù)據(jù)，因此，可將其分析方法分為時(shí)域分析法和頻域分析法，針對(duì)頻域輸出數(shù)據(jù)的分析，主要是通過利用試驗(yàn)結(jié)果和仿真模型模擬數(shù)據(jù)的相互關(guān)系數(shù)來衡量兩個(gè)響應(yīng)之間的近似度。現(xiàn)階段使用的頻域分析方法主要有：基于相關(guān)系數(shù)的頻率響應(yīng)函數(shù)（FRF）方法、傅里葉變換估算功率譜密度等，但此類頻域分析方法的工作原理，是通過從響應(yīng)數(shù)據(jù)中提取能夠表征動(dòng)態(tài)仿真系統(tǒng)特征的變量，并不能給出實(shí)際試驗(yàn)數(shù)據(jù)和仿真數(shù)據(jù)之間的標(biāo)準(zhǔn)化差值，因此，還需有其他驗(yàn)證方法來進(jìn)行性定量分析比較。

汽車碰撞領(lǐng)域的動(dòng)態(tài)系統(tǒng)輸出數(shù)據(jù)具有高度的時(shí)間相關(guān)性及時(shí)變性，主要表現(xiàn)為時(shí)域輸出響應(yīng)特性，因此，本節(jié)將重點(diǎn)從時(shí)域分析的角度進(jìn)行模型驗(yàn)證方法的相關(guān)研究。

Sebdur等提出的AVASIM模型驗(yàn)證方法，已經(jīng)在汽車領(lǐng)域得到應(yīng)用，該方法是基于動(dòng)態(tài)系統(tǒng)響應(yīng)的統(tǒng)計(jì)特性提出，以范數(shù)為理論計(jì)算基礎(chǔ)，此類方法雖然在計(jì)算結(jié)果上存在較大的主觀性，其應(yīng)用也具有一定的局限，但任可為其他模型驗(yàn)證方法提供理論基礎(chǔ)[4]。Liu等則提出通過互相關(guān)系數(shù)對(duì)輸出響應(yīng)進(jìn)行相似度比較，相較于其他方法，互相關(guān)系數(shù)可進(jìn)一步減小相位誤差，但其對(duì)響應(yīng)數(shù)據(jù)的幅度和頻率等特征差異的敏感度較低，其應(yīng)用范圍也因此受限。Geers等提出了時(shí)間序列誤差度量方法，而后Sprague等又對(duì)其進(jìn)行了優(yōu)化，該方法能夠在對(duì)響應(yīng)數(shù)據(jù)進(jìn)行綜合度量時(shí)，有效的將幅度誤差和相位誤差進(jìn)行區(qū)分，但由于該方法不具有對(duì)稱性，從而導(dǎo)致計(jì)算結(jié)果對(duì)比較基準(zhǔn)的依賴性較高，同時(shí)也忽略了形狀誤差對(duì)驗(yàn)證結(jié)果的影響。Russe等通過對(duì)幅度誤差因子進(jìn)行修正，并和相位誤差組合成類似于S&G指標(biāo)的綜合誤差，在一定程度上解決了上述的不對(duì)稱問題，但形狀誤差的影響還是未能解決。Gehre等利用走廊評(píng)分法和互相關(guān)評(píng)分法的結(jié)合，提出了一種響應(yīng)曲線評(píng)分法，可以有效彌補(bǔ)上述方法中忽略了形狀誤差的不足，但其中需要人為設(shè)置的參數(shù)過多，試驗(yàn)人員主觀因素對(duì)實(shí)驗(yàn)結(jié)果影響較大，因此其可信度還有待進(jìn)一步提高。

Sarin等為能夠全面考慮仿真模型輸出響應(yīng)的曲線特征之間的差異，提出了一種基于時(shí)間響應(yīng)的誤差評(píng)估方法（CEARTH），該方法對(duì)相位、幅度、頻率三個(gè)獨(dú)立的誤差進(jìn)行測量，并利用互相關(guān)系數(shù)等方法，對(duì)三類誤差進(jìn)行分離和獨(dú)立測量，從而進(jìn)一步提高了汽車安全領(lǐng)域的模型驗(yàn)證的科學(xué)性和準(zhǔn)確性，但由于在誤差分離時(shí)還存在一定缺陷，導(dǎo)致最終評(píng)分結(jié)果還是具有一定誤差。為進(jìn)一步提升該方法度量結(jié)果的直觀性，詹振飛等在此方法的基礎(chǔ)上進(jìn)行了改進(jìn)，提出了一種增強(qiáng)的時(shí)間響應(yīng)誤差評(píng)估（EEARTH）方法，但評(píng)分結(jié)果的誤差問題還是未能解決[5]。雷正保等提出的多元響應(yīng)分析（MRA）方法，則是專門針對(duì)模型驗(yàn)證中輸出響應(yīng)間強(qiáng)相關(guān)性問題的，同時(shí)開發(fā)了配套的驗(yàn)證工具，該方法有效提升了汽車安全領(lǐng)域仿真模型的驗(yàn)證效率和魯棒性[6]。

近些年來，各類仿真模型的復(fù)雜程度大幅提升，有學(xué)者開始探索將人工智能、機(jī)器學(xué)習(xí)、數(shù)據(jù)挖掘、神經(jīng)網(wǎng)絡(luò)等新方法引入到模型驗(yàn)證中來，希望能夠給復(fù)雜動(dòng)態(tài)系統(tǒng)的仿真模型驗(yàn)證問題帶來新的突破。Birta等借助信息化手段，對(duì)專家知識(shí)或觀測數(shù)據(jù)進(jìn)行形象化描述，從而構(gòu)建模型驗(yàn)證知識(shí)庫；Min等同樣是基于知識(shí)開展模型驗(yàn)證，但對(duì)行為劃分更加合理，評(píng)分指標(biāo)更加詳細(xì)；Mamen等嘗試將神經(jīng)網(wǎng)絡(luò)算法引入模型驗(yàn)證當(dāng)中；DONAHUE等則在將數(shù)據(jù)挖掘、集成學(xué)習(xí)等方法應(yīng)用到模型驗(yàn)證中做了有效探索，上述新方法是國內(nèi)外學(xué)者將現(xiàn)代信息技術(shù)充分融入模型驗(yàn)證方法做出的探索實(shí)踐，其在在驗(yàn)提高驗(yàn)證效率等方面具有一定優(yōu)勢，同時(shí)在驗(yàn)證結(jié)果的可信度，驗(yàn)證工具的智能化程度等方面也還有較大提升空間。

3.2 仿真模型驗(yàn)證工具

仿真模型、尤其是復(fù)雜的動(dòng)態(tài)系統(tǒng)仿真模型的驗(yàn)證工作涉及到大量的數(shù)據(jù)分析與計(jì)算，工作量巨大且繁瑣，僅依靠人力很難完成，因此學(xué)者們?cè)谘芯糠抡婺Ｐ万?yàn)證方法的同時(shí)，也會(huì)針對(duì)相應(yīng)的模型驗(yàn)證工具開展研究。

國外已有模型驗(yàn)證工具主要包括：Holms等針對(duì)導(dǎo)彈系統(tǒng)開發(fā)的模型驗(yàn)證工具（Random Tool Box， RTB ），該工具箱可提供多種模型驗(yàn)證方法給用戶進(jìn)行選擇；Tsai等針對(duì)大型分布式仿真系統(tǒng)的開發(fā)和評(píng)估，設(shè)計(jì)了一個(gè)面向服務(wù)的建模與仿真框架；以及Ferson等開發(fā)的不確定性量化與分析工具，RAMAS Risk等開發(fā)的基于區(qū)間的概率邊界自動(dòng)化計(jì)算工具、美國海軍開發(fā)的VDT工具等。

國內(nèi)的相關(guān)研究主要集中在高校，例如，哈爾濱工業(yè)大學(xué)開發(fā)的基于層次評(píng)估開發(fā)了仿真系統(tǒng)可信度評(píng)估輔助工具，上海交通大學(xué)開發(fā)了模型驗(yàn)證系統(tǒng)，南京理工大學(xué)開發(fā)了W&A軟件包，北京理工大學(xué)開發(fā)了運(yùn)動(dòng)體系統(tǒng)仿真模型驗(yàn)證軟件，北京航空航天大學(xué)開發(fā)了仿真系統(tǒng)校核與驗(yàn)證的自測軟件，國防科學(xué)技術(shù)大學(xué)開發(fā)了制導(dǎo)仿真系統(tǒng)W&A軟件工具，長沙理工大學(xué)開發(fā)了動(dòng)態(tài)系統(tǒng)仿真模型輔助驗(yàn)證工具[7]，空軍工程大學(xué)開發(fā)了M&S VV&A工具等[8]，樹優(yōu)公司開發(fā)了應(yīng)用于汽車碰撞實(shí)驗(yàn)?zāi)Ｐ偷尿?yàn)證工具OntoTest等。

現(xiàn)有的模型驗(yàn)證工具還存在功能較為單一、多目標(biāo)協(xié)同評(píng)估性能較差、便利性和可拓展性不足等問題，隨著仿真技術(shù)的發(fā)展，各類復(fù)雜動(dòng)態(tài)系統(tǒng)仿真模型對(duì)支撐工具平臺(tái)的可拓展性、操作便利性、功能完備性、性能穩(wěn)定性等要求越來越高，如何提供一個(gè)更加開放、規(guī)范的在線協(xié)同環(huán)境，成為支撐復(fù)雜仿真模型驗(yàn)證工作的迫切需求。

4 結(jié)論

動(dòng)態(tài)系統(tǒng)仿真模型是汽車被動(dòng)安全領(lǐng)域研究非常重要的數(shù)字化工具，同時(shí)也是驅(qū)動(dòng)復(fù)雜仿真系統(tǒng)運(yùn)行的核心，仿真模型驗(yàn)證則是最終仿真結(jié)果是否可信的重要依據(jù)，因此，仿真模型的可信度會(huì)直接影響仿真結(jié)果的可信度。本文在分析總結(jié)汽車動(dòng)態(tài)系統(tǒng)仿真模型及其驗(yàn)證方法發(fā)展歷程的基礎(chǔ)上，重點(diǎn)針對(duì)從時(shí)域分析的角度對(duì)動(dòng)態(tài)系統(tǒng)仿真模型驗(yàn)證方法進(jìn)行了歸納分析和分類，同時(shí)，對(duì)國內(nèi)外主要模型驗(yàn)證工具存在的問題進(jìn)行了梳理，并指出了復(fù)雜動(dòng)態(tài)系統(tǒng)仿真模型驗(yàn)證方法及工具存在的若干問題和未來發(fā)展方向。

基金項(xiàng)目：湖南省教育廳科學(xué)研究項(xiàng)目（20C1004）：汽車復(fù)雜動(dòng)態(tài)系統(tǒng)模型驗(yàn)證的多元響應(yīng)分析方法研究。

參考文獻(xiàn)：

[1]Conway R W，Johnson B M，Maxwell W L. Some problems of digital systems simulation[J]. Management Science（S0025-1909），1959，6（1）：92-110.

[2]Fishman G S，Kiviat P J. The analysis of simulation generated time series[J]. Management Science.

（S0025-1909），1967，3（7）： 525-557.

[3]Balci O，Sargent R G. Some examples of simulation model validation using hypothesis testing[C]. Proceedings of the 1982 Winter Simulation Conference.

[4]Sendur P，Stein J L，Peng H，et al. A model accuracy and validation algorithm[J].Dynamic Systems and Control，2002，71：573-583.

[5]詹振飛，胡潔，符雁，等. 虛擬樣機(jī)環(huán)境下不確定性多元?jiǎng)討B(tài)系統(tǒng)模型驗(yàn)證[J]. 機(jī)械工程學(xué)報(bào)，2012，48（5）： 138-146.

[6]雷正保，邢歡，陳小勇.汽車多元?jiǎng)討B(tài)系統(tǒng)仿真模型響應(yīng)分析方法[J].汽車技術(shù)，2020（08）：41-45.

[7]陳小勇. 動(dòng)態(tài)系統(tǒng)仿真模型驗(yàn)證方法及工具研究[D].長沙理工大學(xué)，2019.

[8]馬震，吳曉燕，卜祥偉等.仿真模型VV&A工具研究[J].現(xiàn)代防御技術(shù)，2016，44（1）：225-229.