馬智博， 李海杰， 殷建偉， 黃文斌

（1.北京應(yīng)用物理與計算數(shù)學(xué)研究所，北京 100088；2.長沙機電產(chǎn)品研究開發(fā)中心，湖南長沙 410100；3.流體物理研究所，四川綿陽 621900）

可靠性數(shù)值模擬的不確定度量化

馬智博1， 李海杰2， 殷建偉1， 黃文斌3

（1.北京應(yīng)用物理與計算數(shù)學(xué)研究所，北京 100088；2.長沙機電產(chǎn)品研究開發(fā)中心，湖南長沙 410100；3.流體物理研究所，四川綿陽 621900）

根據(jù)可靠性認(rèn)證對不確定度量化的技術(shù)要求，結(jié)合復(fù)雜工程系統(tǒng)的層級結(jié)構(gòu)和數(shù)值模擬從校準(zhǔn)、驗證與確認(rèn)到最終具有預(yù)測能力的歷程，對可靠性數(shù)值模擬不確定度量化的要求和方法進(jìn)行探索，并結(jié)合爆轟算例對這些方法進(jìn)行演示和驗證.

數(shù)值模擬；不確定度量化；驗證與確認(rèn)；可靠性認(rèn)證

0 引言

可靠性問題源自產(chǎn)品本身和任務(wù)環(huán)境中A、B兩類不確定度的存在，其中前者對應(yīng)于客觀隨機性引起的偶然不確定度（aleatory uncertainty），后者對應(yīng)于主觀認(rèn)知缺陷帶來的認(rèn)知不確定度（epistemic uncertainty）.根據(jù)對不確定問題的存在是否預(yù)知，認(rèn)知不確定度又分為已知的未知（known unknowns）和未知的未知（unknown unknowns）［1］.為了保證產(chǎn)品工作可靠，必須在綜合考慮兩類不確定度的前提下，為定型的批量產(chǎn)品留下足夠的設(shè)計裕度.

實物試驗和統(tǒng)計分析可用于量化這兩類不確定度，也可用于消減認(rèn)知不確定度并發(fā)現(xiàn)其中未知的未知，但偶然不確定度只能通過改進(jìn)設(shè)計方案和制造工藝來實現(xiàn)其消減.

數(shù)值模擬包含建模和模擬（modeling and simulation，M＆S）兩個技術(shù)環(huán)節(jié)，建模和模擬均受人類主觀意志的主導(dǎo)，只能在人的認(rèn)識基礎(chǔ)上借助計算機完成對客觀現(xiàn)象的分析，所以數(shù)值模擬主要用于研究偶然不確定度以及認(rèn)知不確定度中已知的未知，一般不直接用于研究認(rèn)知不確定度中未知的未知.

當(dāng)受到成本因素或其它條件制約時，很多大型復(fù)雜產(chǎn)品往往難以實施系統(tǒng)級可靠性試驗，此時，需要用M＆S技術(shù)對可靠性特征量進(jìn)行工程預(yù)測，而基于裕度和不確定度量化的QMU認(rèn)證方法（quantification of margins and uncertainties，QMU）為M＆S在可靠性認(rèn)證中發(fā)揮巨大作用提供了科學(xué)的技術(shù)平臺，M＆S預(yù)測結(jié)果以及預(yù)測結(jié)果的不確定度將是QMU認(rèn)證的重要信息來源［2－4］.

可靠性數(shù)值模擬的任務(wù)主要表現(xiàn)在兩個方面：其一是研究產(chǎn)品狀態(tài)和任務(wù)環(huán)境與產(chǎn)品性能之間復(fù)雜而又確定的物理關(guān)系（對隨機過程而言，統(tǒng)計意義下的物理關(guān)系是確定的），發(fā)現(xiàn)失效模式并量化可靠性特征量，其二是研究產(chǎn)品狀態(tài)和任務(wù)環(huán)境的不確定度所引起的產(chǎn)品性能的不確定度，促成不確定度的傳遞分析.因為后者的不確定度起源于工程因素而非數(shù)值模擬本身，所以數(shù)值模擬的不確定度主要體現(xiàn)在反映前者對應(yīng)的確定性物理關(guān)系時所帶來的模擬結(jié)果的不確定度.

對M＆S實施不確定度量化（uncertainty quantification，UQ），本身也是一項復(fù)雜的系統(tǒng)工程，相關(guān)技術(shù)正處在發(fā)展中［5－7］.本文根據(jù)可靠性認(rèn)證的要求，對M＆S不確定度的基本概念、組成結(jié)構(gòu)、量化原則、量化方法和方法檢驗進(jìn)行探索，并以爆轟計算為例展示這些方法的思想基礎(chǔ)和應(yīng)用效果.

1 數(shù)值模擬不確定度量化的基本原則

不確定度量化的源頭信息往往是誤差數(shù)據(jù)，綜合不同行業(yè)對誤差和不確定度的定義［8－13］，可以從廣義的角度上總結(jié)出兩者的普適性內(nèi)涵.本文認(rèn)為，誤差是行為結(jié)果與期望目標(biāo)之間可以辨識的差別，它來自有期望目標(biāo)的行為，沒有期望目標(biāo)的行為不產(chǎn)生誤差，不能辨識的差別形不成誤差，而不確定度是行為結(jié)果的分散性或背離期望目標(biāo)的程度，對無期望目標(biāo)的行為，不確定度主要表現(xiàn)為行為結(jié)果的分散性，對有期望目標(biāo)的行為，不確定度主要表現(xiàn)為行為結(jié)果背離期望目標(biāo)的程度.

可靠性數(shù)值模擬的期望目標(biāo)是準(zhǔn)確預(yù)測客觀物理過程中的確定性規(guī)律，所以既產(chǎn)生誤差，也產(chǎn)生不確定度，其中誤差為模擬結(jié)果減去參照解所得的差值，基于可靠性工程的需要，不確定度應(yīng)被理解為模擬結(jié)果背離真值的程度.

無論可靠性認(rèn)證（針對抽象的型號）還是可靠性評估（針對具體的產(chǎn)品），都需要對可靠性特征量的真值作出估計，根據(jù)QMU方法對應(yīng)的數(shù)據(jù)格式，該估計結(jié)果包含兩個部分，其一是最佳估計結(jié)果，其二是針對該最佳估計的不確定度估計結(jié)果，兩者共同形成關(guān)于可靠性特征量估計的不確定度區(qū)間.

再考慮到其它方面的因素，數(shù)值模擬的不確定度量化應(yīng)該遵循以下原則：

1）UQ的前提是M＆S的技術(shù)狀態(tài)得到相對固化且計算過程不依賴人為干預(yù)；

2）UQ的目標(biāo)是量化對應(yīng)于工程預(yù)測且無參照解時數(shù)值模擬的不確定度；

3）可靠性數(shù)值模擬的不確定度應(yīng)該理解為模擬結(jié)果背離真值的程度；

4）模擬結(jié)果和不確定度估計值所形成的不確定區(qū)間要立足于覆蓋真值；

5）基于給定輸入信息，方法上要盡量使不確定度的量化結(jié)果達(dá)到最小；

6）量化方法要有利于診斷和消除M＆S存在的問題以促進(jìn)不確定度消減；

7）量化方法要有利于保證量化信息的充分性和量化結(jié)果的可信性.

2 數(shù)值模擬不確定度的組成結(jié)構(gòu)

對大部分工程問題，M＆S的本質(zhì)是求解物理系統(tǒng)初始狀態(tài)和邊界條件與其最終性能的映射關(guān)系，該映射關(guān)系可表示為一個泛函［14］

其中x為物質(zhì)點的空間坐標(biāo)；t為時間變量，下標(biāo)“0”和“e”分別表示初始和最終時刻；ξ為物質(zhì)點的物理狀態(tài)參數(shù)組，反映材料的物理狀態(tài)；ω為隨空間和時間而變化的邊界條件函數(shù)組，反映產(chǎn)品的工作環(huán)境；η為由物質(zhì)行為屬性決定的描述物質(zhì)行為規(guī)律的函數(shù)組，是物理狀態(tài)ξ的函數(shù)；ζ為隨時間變化的系統(tǒng)性能參數(shù)組；P為由普遍物理規(guī)律決定的聯(lián)結(jié)系統(tǒng)性能和眾多參數(shù)或函數(shù)的泛函關(guān)系式.

建模與模擬的不確定度，指的是用模擬計算結(jié)果來預(yù)測對應(yīng)實體模型的真實行為所帶來的不確定度，根據(jù)式（1），該不確定度可分解為三個組成部分：U1為實體建模不確定度U01所引起的模擬結(jié)果不確定度，實體建模即構(gòu)造ξ（x，t0）和ω（x，t）；U2為物理建模不確定度U02所引起的模擬結(jié)果不確定度，物理建模即構(gòu)造η（ξ）和P［·］；U3為數(shù)值計算所引起的模擬結(jié)果不確定度，數(shù)值計算即用數(shù)值方法求解泛函P［·］.

因此，M＆S的總不確定度可表示為

其中符號“⊕”表示總不確定度由三個部分組成，在某些情況下，總不確定度可近似表示為右端三項簡單相加所得的總和.

對建模引起的不確定度U1和U2，需要首先量化建模的不確定度U01和U02，然后通過傳遞分析得到U1和U2，U01和U02均由兩部分構(gòu)成，其一是模型形式的不確定度Uform，其二是模型輸入?yún)?shù)的不確定度Upara．?dāng)?shù)值計算引起的不確定度U3，主要源自數(shù)值離散、迭代計算和四舍五入等.

3 不確定度量化的技術(shù)環(huán)節(jié)

數(shù)值模擬不確定度量化所涉及到的技術(shù)環(huán)節(jié)及其對應(yīng)方法主要有以下幾個方面.

1）物理系統(tǒng)的層級劃分

可根據(jù)物理過程從簡單到復(fù)雜的邏輯關(guān)系或從開始到結(jié)束的先后順序，將物理系統(tǒng)劃分為不同的層級，層級劃分有助于隔離和診斷M＆S存在的問題，也有利于獲取更多的不確定度信息.

2）建模的不確定度量化

主要目的是量化U01和U02．對不同特性的模型，可基于相應(yīng)的測量技術(shù)、科學(xué)理論和實驗統(tǒng)計方法來獲取U01和U02的量化信息.

3）不確定度的傳遞分析

常用于對建模的不確定度U01和U02進(jìn)行傳遞以得到U1和U2，該過程可借助數(shù)值模擬程序用變動分析的方法來完成，兩者分別用變分簡記為

4）模擬結(jié)果與參照解的對比分析

該方法主要基于數(shù)值模擬的驗證與確認(rèn)（verification＆validation，V＆V）技術(shù)，通過與高可信度參照解的對比來量化數(shù)值模擬的不確定度［1］，基準(zhǔn)模型的設(shè)計及其參照解的構(gòu)造是V＆V的重要技術(shù)支撐.

5）不確定度的預(yù)測推斷

數(shù)值模擬的實際價值體現(xiàn)在對新模型的工程預(yù)測，但此時的不確定度不能通過對比分析來獲取，需要基于已有的對比信息以及對不確定度隨模型參數(shù)變化規(guī)律的認(rèn)識來綜合推斷預(yù)測結(jié)果的不確定度.

6）不確定度信息的整合

對一個新系統(tǒng)進(jìn)行數(shù)值模擬，其不確定度信息至少有兩個來源，其一是建模不確定度傳遞到系統(tǒng)級模擬結(jié)果上的不確定度信息，其二是基于各層級對比分析和預(yù)測推斷得到的系統(tǒng)級模擬結(jié)果不確定度信息，另外還可能有不同程序之間計算結(jié)果的對比信息，這些信息可通過廣義信息理論進(jìn)行整合.

4 應(yīng)用域內(nèi)不確定度的預(yù)測推斷

根據(jù)文獻(xiàn)［1］，可將實體模型設(shè)計參數(shù)（包括環(huán)境參數(shù)）所對應(yīng)的多維參數(shù)空間劃分為確認(rèn)域和應(yīng)用域，確認(rèn)域即已經(jīng)實施的物理試驗所對應(yīng)的參數(shù)空間，應(yīng)用域即數(shù)值模擬預(yù)期應(yīng)用所對應(yīng)的參數(shù)空間，由于應(yīng)用域內(nèi)一般沒有試驗數(shù)據(jù)用于對比分析，所以其不確定度信息需要根據(jù)確認(rèn)域內(nèi)的不確定度以及兩個區(qū)域之間不確定度的某種聯(lián)系來進(jìn)行推斷.

為了揭示確認(rèn)域和應(yīng)用域之間不確定度的聯(lián)系方式，文獻(xiàn)［14］提出了不確定度與實體模型的設(shè)計參數(shù)存在函數(shù)關(guān)系的假設(shè)，并結(jié)合有理論解的黎曼問題對基于該假設(shè)的應(yīng)用域不確定度預(yù)測量化方法進(jìn)行了驗證分析，本文擬對該方法的思想基礎(chǔ)做進(jìn)一步完善，并結(jié)合試驗?zāi)Ｐ蛯υ摲椒ǖ挠行赃M(jìn)行檢驗.

本文就應(yīng)用域內(nèi)不確定度預(yù)測推斷的科學(xué)依據(jù)提出如下觀點：

1）用規(guī)定的工具和方法對物理量（如尺寸，質(zhì)量，密度，溫度，速度等）進(jìn)行測量，測量結(jié)果的不確定度是被測量的函數(shù)；

2）用規(guī)定的運算法則由觀測量推算其它導(dǎo)出量，導(dǎo)出量的不確定度是觀測量的函數(shù)；

3）計算程序固化后，計算環(huán)節(jié)引起的不確定度是實體模型設(shè)計參數(shù)的函數(shù)；

4）基于以上觀點，數(shù)值模擬不確定度可表示為實體模型設(shè)計參數(shù)的函數(shù)：

5）該函數(shù)既適用于確認(rèn)域也適用于應(yīng)用域.

根據(jù)以上五點，如果在確認(rèn)域獲取了不確定度隨實體模型設(shè)計參數(shù)的變化規(guī)律，就可將該函數(shù)內(nèi)推或外推以獲得應(yīng)用域內(nèi)M＆S的不確定度.

5 數(shù)值模擬不確定度量化的步驟

根據(jù)數(shù)值模擬不確定度量化的原則和V＆V的技術(shù)思想，建議按如下步驟實施數(shù)值模擬的不確定度量化：

1）根據(jù)可靠性認(rèn)證的工程需求，確認(rèn)數(shù)值模擬的不確定度要求；

2）根據(jù)物理過程從簡單到復(fù)雜的邏輯關(guān)系或從開始到結(jié)束的先后順序，將物理系統(tǒng)劃分為不同的層級，低層級的模型有利于診斷和解決M＆S中的問題以及不確定度信息的獲??；

3）根據(jù)不同層級中物理現(xiàn)象的特點和耦合程度，分析對應(yīng)的算法中需要重點考核的技術(shù)要素，基于這些考核要素為代碼驗證逐層建立基準(zhǔn)模型并獲取基準(zhǔn)模型的高可信度參照解（比如理論解、人為解或由其它程序得到的被公認(rèn)為比較可信的數(shù)值解）；

4）結(jié)合基準(zhǔn)模型開展代碼驗證（code verification），發(fā)現(xiàn)并隔離算法設(shè)計和代碼編制中存在的錯誤，驗證離散方程的相容性以及數(shù)值求解的收斂性和穩(wěn)定性，保證計算程序能夠忠實地表現(xiàn)數(shù)學(xué)模型.代碼驗證是解算驗證的前提與基礎(chǔ)，其中計算誤差來自基準(zhǔn)模型數(shù)值解與參照解的對比，可通過對該誤差進(jìn)行統(tǒng)計分析和區(qū)間估計以實現(xiàn)對計算不確定度的評估；

5）根據(jù)M＆S的應(yīng)用需求在不同層級上開展解算驗證（solution verification），以驗證數(shù)值計算是否滿足給定的準(zhǔn)確度要求.解算驗證也要估計輸入數(shù)據(jù)和后處理結(jié)果的不確定度，它是確認(rèn)活動的前提和技術(shù)基礎(chǔ)，所對應(yīng)的模型較復(fù)雜，一般沒有用于對比分析的參照解，可通過數(shù)值解隨離散尺度的收斂趨勢來估計其收斂解——如Richardson外推方法等，并借助帶安全系數(shù)的網(wǎng)格收斂指標(biāo)（grid convergence index，GCI）來估計數(shù)值計算的不確定度，并使該結(jié)果滿足不確定度量化的真值覆蓋原則；

6）針對所需的實體建模和物理建模，設(shè)計并實施各個層級的校準(zhǔn)試驗，對建模和模擬的模型參數(shù)、計算參數(shù)和可調(diào)參數(shù)進(jìn)行校準(zhǔn).可調(diào)參數(shù)指的是無物理意義或雖然有物理意義但存在較大認(rèn)知缺陷的對M＆S結(jié)果有明顯影響的建?；蛴嬎銋?shù)；

7）開展各個層級的確認(rèn)試驗，用校準(zhǔn)后的模型參數(shù)進(jìn)行數(shù)值模擬，通過數(shù)值模擬結(jié)果與試驗結(jié)果的對比信息，量化各層級數(shù)值模擬的不確定度，并根據(jù)確認(rèn)過程發(fā)現(xiàn)的問題，完成對建模和模擬的技術(shù)改進(jìn)，以縮減數(shù)值模擬不確定度；

8）通過預(yù)測推斷由系統(tǒng)級確認(rèn)域內(nèi)的不確定度獲得應(yīng)用域內(nèi)的不確定度；

9）通過傳遞分析獲得由建模不確定度所引起的系統(tǒng)級數(shù)值模擬不確定度；

10）將預(yù)測推斷和傳遞分析得到的不確定度信息進(jìn)行整合，實現(xiàn)應(yīng)用域內(nèi)數(shù)值模擬總體不確定度的量化.

6 實例演示

6.1 實體模型

該實體模型主要包含一個起爆器和一塊炸藥，在軸線上炸藥的厚度為δ，該模型及其試驗裝置主要研究炸藥起爆初期的爆轟波發(fā)展過程，用以建立該階段的起爆和傳爆物理模型.測量的參數(shù)為爆轟波傳到炸藥右端面時在軸線上貼近LiF窗口的炸藥峰值速度，該速度的相對測量誤差為±3.0%，如圖1所示.

圖1 實體模型及試驗設(shè)備Fig.1 The entitymodel and test facility

6.2 層級劃分

將該實體模型視為一個系統(tǒng)，并劃分為兩個層級：

1）起爆器層級，僅有起爆器，沒有炸藥件；

2）全系統(tǒng)層級，包含起爆器和炸藥件.

起爆器層級要先于系統(tǒng)層級開展校準(zhǔn)和確認(rèn)活動，用于建立和評價起爆器動作過程的數(shù)值模擬可信度，在此基礎(chǔ)上，再開展全系統(tǒng)層級的校準(zhǔn)和確認(rèn)活動.

6.3 模型校準(zhǔn)

進(jìn)行模型確認(rèn)之前，要先進(jìn)行模型校準(zhǔn)，即通過試驗對物理模型的參數(shù)進(jìn)行校準(zhǔn)，在校準(zhǔn)之前，要首先選定物理模型的形式，模型形式不在此細(xì)述.該模型的校準(zhǔn)活動分兩步：

1）起爆器層級校準(zhǔn)

通過起爆器外半球形金屬薄膜中心位置的峰值速度測量值來校準(zhǔn)起爆器內(nèi)數(shù)值模擬的物理模型和計算參數(shù)；

2）全系統(tǒng)層級校準(zhǔn)

選擇不同的炸藥厚度，通過炸藥右端面峰值速度測量值來校準(zhǔn)炸藥對應(yīng)起爆和傳爆過程的物理模型和計算參數(shù).

6.4 模型確認(rèn)

該實體模型只有一個設(shè)計變量，即炸藥的厚度δ，炸藥密度不作為設(shè)計變量.

該例主要關(guān)心炸藥峰值速度V的數(shù)值模擬不確定度隨炸藥厚度的變化規(guī)律，炸藥厚度δ和炸藥密度均具有加工隨機性，數(shù)值模擬時要輸入它們的實測值.

確認(rèn)試驗分4組，δ的設(shè)計值分別為5 mm、10 mm、15 mm和20 mm，每組共有2次重復(fù)性試驗，由于加工隨機性，即使重復(fù)性試驗，炸藥樣品的厚度和密度也不會嚴(yán)格等于其設(shè)計值.

在確認(rèn)環(huán)節(jié)，固定校準(zhǔn)后的模型參數(shù)（比如JWL狀態(tài)方程中的5個模型參數(shù)）不變，分別獲得8個實體模型的峰值速度模擬結(jié)果，如表1所示，對應(yīng)的試驗結(jié)果如表2所示.

表1 確認(rèn)域內(nèi)峰值速度V的模擬結(jié)果（δ／mm，V／（m·s－1））Table 1 Simulation results of V in the validated domain（δ／mm，V／（m·s－1））

表2 確認(rèn)域內(nèi)峰值速度V的試驗結(jié)果（δ／mm，V／（m·s－1））Table 2 Test results of V in validated domain（δ／mm，V／（m·s－1））

本文將模擬誤差定義為模擬結(jié)果減去試驗結(jié)果：

由于每組的重復(fù)試驗只有兩次，我們根據(jù)這兩次試驗中模擬結(jié)果對試驗結(jié)果的最大背離來近似提取確認(rèn)域內(nèi)建模與模擬的不確定度信息，如表3所示.

6.5 應(yīng)用域內(nèi)的不確定度量化

確認(rèn)域內(nèi)炸藥厚度δ的設(shè)計范圍是5 mm～20 mm，本文將炸藥厚度δ分別為2 mm、3 mm、4 mm和25 mm的實體模型劃入應(yīng)用域內(nèi)，根據(jù)δ的取值范圍，該例的確認(rèn)域和應(yīng)用域不重合.

在現(xiàn)實情況下，應(yīng)用域內(nèi)沒有試驗值可以參照，不可能通過與試驗結(jié)果的對比來量化M＆S不確定度，必須建立不確定度的預(yù)測模型.

基于式（5）所表達(dá)的基本思想，擬用二次多項式構(gòu)造不確定度預(yù)測模型，獲取不確定度隨炸藥厚度δ的變化規(guī)律

表3 確認(rèn)域內(nèi)峰值速度V的模擬誤差和不確定度（δ／mm，EV／（m·s－1），UV／（m·s－1））Table 3 Simulation errors＆uncertainties of V in validated domain（δ／mm，EV／（m·s－1），UV／（m·s－1））

方程（8）一般為超定方程，根據(jù)最小二范數(shù)原則進(jìn)行求解，得到

至此，式（7）所對應(yīng)的函數(shù)形式被完全確定下來，并由此得到應(yīng)用域內(nèi)的不確定度預(yù)測結(jié)果，如表 4所示.

表4 應(yīng)用域內(nèi)峰值速度V的M＆S不確定度預(yù)測值（δ／mm，UV／（m·s－1））Table 4 Predicted M＆S uncertainties of V in the applied domain（δ／mm，UV／（m·s－1））

7 對不確定度量化方法的檢驗

本文僅對基于預(yù)測推斷的應(yīng)用域內(nèi)不確定度量化方法進(jìn)行檢驗.

基于上述算例模型，本文安排以下步驟進(jìn)行檢驗：

1）根據(jù)預(yù)測模型，推斷應(yīng)用域內(nèi)對應(yīng)不同炸藥厚度δ的M＆S不確定度，如表4所示；

2）對應(yīng)用域內(nèi)每個δ設(shè)計值，加工2個試驗?zāi)Ｐ?，?個試驗?zāi)Ｐ?，制定試驗計劃?/p>

3）對8個試驗?zāi)Ｐ瓦M(jìn)行初態(tài)參數(shù)的測量，在試驗前完成數(shù)值模擬并凍結(jié)模擬結(jié)果；

4）完成8個試驗?zāi)Ｐ偷脑囼?，記錄試驗的測試結(jié)果；

5）將模擬結(jié)果減去測試結(jié)果，得到8次試驗的模擬誤差；

6）對比模擬誤差和推斷結(jié)果以檢驗應(yīng)用域內(nèi)不確定度量化方法的準(zhǔn)確性.

對比結(jié)果列于表5.將表5中每個δ設(shè)計值下兩次重復(fù)試驗的最大誤差絕對值與表4中由預(yù)測模型得到的不確定度進(jìn)行對比，發(fā)現(xiàn)預(yù)測模型能基本上反映出不確定度隨δ的變化規(guī)律，具有一定的預(yù)測推斷效果.由于試驗數(shù)太少，統(tǒng)計結(jié)果的隨機性較大，應(yīng)用域內(nèi)的量化結(jié)果還不能完全達(dá)到真值覆蓋的要求，并且隨著設(shè)計點遠(yuǎn)離確認(rèn)域，預(yù)測的準(zhǔn)確性也會逐漸降低.

表5 應(yīng)用域內(nèi)峰值速度V的模擬結(jié)果和模擬誤差（δ／mm，V／（m·s－1））Table 5 Simulation and test results and errors of V in applied domain（δ／mm，V／（m·s－1））

8 結(jié)論

在缺乏試驗的前提下，全系統(tǒng)級的可靠性認(rèn)證主要依靠數(shù)值模擬結(jié)果以及該結(jié)果對應(yīng)的不確定度，對不確定度過大或過小的估計將顯著增加產(chǎn)品認(rèn)證決策中“棄真”和“納偽”的錯誤，為供需雙方帶來損失.本文根據(jù)復(fù)雜工程系統(tǒng)的可靠性認(rèn)證將越來越多地依賴數(shù)值模擬這一發(fā)展趨勢，對可靠性數(shù)值模擬不確定度量化的概念、原則、思想方法與實施步驟進(jìn)行了分析和歸納，重點對數(shù)值模擬用于預(yù)測任務(wù)時應(yīng)用域內(nèi)的不確定度量化方法進(jìn)行了探索，并結(jié)合算例演示了不確定度量化所經(jīng)歷的層級劃分、模型校準(zhǔn)、驗證與確認(rèn)、應(yīng)用域內(nèi)不確定度量化及其方法檢驗等過程.檢驗結(jié)果表明，以實體模型參數(shù)為自變量而建立的不確定度函數(shù)，能夠基本上反映數(shù)值模擬不確定度在該參數(shù)空間中的變化規(guī)律，并能夠用于應(yīng)用域內(nèi)數(shù)值模擬不確定度的推斷，但隨著實體模型遠(yuǎn)離確認(rèn)域，推斷的準(zhǔn)確性會逐漸降低，另外，由于該方法的思想基礎(chǔ)源自對科學(xué)規(guī)律的主觀認(rèn)識，所以還需要有更多的實踐檢驗.

［1］OberkampfWL，Roy C J.Verification and validation in science computing［M］.Cambridge University Press，2010：1－767.

［2］Ma Zhibo，Ying Yangjun，Zhu Jianshi.QMU certifyingmethod and its implementation［J］.Chinese Journal of Nuclear Science and Engineering，2009，29（1）：1－9.

［3］Pilch M，Trucano T G，Helton J C.Ideas underlying the quantification of margins and uncertainties［J］.Reliability Engineering and System Safety，2011，96：965－975.

［4］Helton JC.Quantification ofmargins and uncertainties：Conceptual and computational basis［J］.Reliability Engineering and System Safety，2011，96：976－1013.

［5］Doostan A，Iaccarino G.A least-squares approximation of partial differential equations with high-dimensional random inputs［J］.Journal of Computational Physics，2009，228：4332－4342.

［6］Helton JC，Johnson JD，Sallaberry C J，Storlie C B.Survey of sampling-based method for uncertainty and sensitivity analysis［J］.Reliability Engineering and System Safety，2006，91（10－11）：1175－1209.

［7］Roy C J，OberkampfWL.A comprehensive framework for verification，validation，and uncertainty quantification in scientific computing［J］.Computer Methods in Applied Mechanics and Engineering，2011，200：2131－2144.

［8］Reston，VA，Guide for verification and validation of computational fluid dynamics［S］.American Institute of Aeronautics and Astronautics，AIAA G－077－1998.

［9］Ye Depei，Song Zhenguo，Wang Xianzhi，Expression and evaluation of uncertainty inmeasurement［S］.Military Standard of the People's Republic of China，1999，GJB3756－99.

［10］The American Society ofMechanical Engineers.Testuncertainty［S］.American National Standard，2005，ASME PTC 19.1－2005.

［11］The American Society of Mechanical Engineers.Guide for verification and validation in computational solid mechanics［S］. American National Standard，2006，ASME V＆V 10－2006.

［12］The American Society of Mechanical Engineers.Standard for verification and validation in computational fluid dynamics and heat transfer［S］.American National Standard，2009，ASME V＆V 20－2009.

［13］The American Society of Mechanical Engineers.An illustration of the concepts of verification and validation in computational solid mechanics［S］.American National Standard，2012，ASME V＆V 10.1－2012.

［14］Ma Zhibo，Zheng Miao，Yin Jianwei，et al.Quantification of uncertainties in detonation simulations［J］.Chinese Journal of Computational Physics，2011，28（1）：66－74.

Uncertainty Quantification of Numerical Simulation for Reliability Analysis

MA Zhibo1，LIHaijie2，YIN Jinwei1，HUANGWenbin3
（1.Institute of Applied Physics and Computational Mathematics，Beijing 100088，China；2.Electromechanical Product Research and Development Center，Changsha 410100，China；3.Institute of Fluid Physics，Mianyang 621900，China）

According to demands of uncertainty quantification for reliability certification，in hierarchy of complex engineering system and development from calibration，verification＆validation to prediction ability ofmodeling＆simulation（M＆S），the principle and method of uncertainty quantification of M＆S are investigated.An example on detonation process is shown in which the ideas are demonstrated.

numerical simulation；uncertainty quantification；verification＆validation；reliability certification

date：2013－07－12；Revised date：2013－12－07

0242.1

A

1001-246X（2014）04-0424-07

2013－07－12；

2013－12－07

國家自然科學(xué)基金（11371066）、中國工程物理研究院科學(xué)技術(shù)發(fā)展基金（2013A0101004）及裝備預(yù)先研究（51319010207）資助項目

馬智博（1963－），男，河南，博士，研究員，主要從事計算力學(xué)和可靠性工程研究，E-mail：mazhibo＠iapcm.ac.cn