風洞MDOE的形式實驗設(shè)計方法研究

張 江，秦永明，馬漢東

（中國航天空氣動力技術(shù)研究院，北京 100074）

風洞MDOE的形式實驗設(shè)計方法研究

張 江＊，秦永明，馬漢東

（中國航天空氣動力技術(shù)研究院，北京 100074）

MDOE風洞實驗方法能夠用相對于傳統(tǒng)實驗方法更少的吹風次數(shù)，獲得更高精準度的數(shù)據(jù)。為了解決現(xiàn)有基于參數(shù)模型的MDOE方法獲取較強非線性氣動規(guī)律能力的不足，需要發(fā)展基于非參數(shù)模型的MDOE方法。本文對基于非參數(shù)模型的MDOE的形式實驗設(shè)計方法進行研究。通過“虛擬”風洞實驗的方法，對兩種常用的“空間填充設(shè)計”——拉丁超立方設(shè)計和均勻設(shè)計應(yīng)用于風洞實驗的適用性進行對比，并在此基礎(chǔ)上發(fā)展了幾種對均勻設(shè)計的優(yōu)化改進方法，可以進一步提高樣本點設(shè)計質(zhì)量，使其滿足風洞實驗的要求。研究表明：均勻設(shè)計較拉丁超立方設(shè)計更為穩(wěn)健、均勻，更適合基于非參數(shù)模型的風洞MDOE方法；在均勻設(shè)計方法基礎(chǔ)上，根據(jù)風洞實驗的特點發(fā)展了優(yōu)化方法，包括邊界點補充、樣本點密度調(diào)整和重復(fù)點設(shè)計，能夠?qū)⒁延械摹跋闰炐畔ⅰ睉?yīng)用于實驗設(shè)計中；所發(fā)展的形式實驗設(shè)計方法所需的測量點要少于OFAT方法的測量點（如示例中所用的測量點數(shù)僅為OFAT方法的66.7%），且能夠充分和準確地對較劇烈的非線性變化規(guī)律進行采樣。本文對風洞MDOE的形式實驗設(shè)計方法的研究結(jié)果，為后續(xù)發(fā)展基于非參數(shù)模型的風洞MDOE方法奠定了基礎(chǔ)。

MDOE；非參數(shù)模型；均勻設(shè)計；拉丁超立方；風洞實驗

0 引 言

現(xiàn)代風洞實驗測試系統(tǒng)的構(gòu)成越來越復(fù)雜，增加了實驗中各環(huán)節(jié)的微弱變化所導致的系統(tǒng)誤差來源，使得風洞實驗誤差的控制難度加大，同飛行器設(shè)計不斷提高的數(shù)據(jù)質(zhì)量要求之間的矛盾越來越突出。傳統(tǒng)的單因子實驗方法（OFAT方法）對實驗環(huán)境變化導致的誤差缺少控制和消除能力［1］，且生產(chǎn)效率過低［2］。現(xiàn)代實驗設(shè)計（MDOE）領(lǐng)域的研究表明通過改進實驗設(shè)計，能夠有效地減小甚至消除實驗環(huán)境中不可控因素對實驗結(jié)果的影響。風洞MDOE實驗方法就是將現(xiàn)代實驗設(shè)計理論應(yīng)用于風洞實驗產(chǎn)生的新技術(shù)，它能夠用相對少量的實驗次數(shù)，獲得更高質(zhì)量的數(shù)據(jù)。

MDOE方法通過建立統(tǒng)計模型來獲取氣動規(guī)律。統(tǒng)計模型有方差分析模型、參數(shù)回歸模型、非參數(shù)回歸模型和穩(wěn)健回歸模型等。形式實驗設(shè)計方法是MDOE的重要組成，是對測量樣本點數(shù)量和取值進行設(shè)計的方法。不同的統(tǒng)計模型需要不同的形式實驗設(shè)計方法。方差模型適用于因素線性效應(yīng)的研究，采用具有合適分辨力的部分因子設(shè)計［3］，可以通過較全因子設(shè)計更少的實驗次數(shù)，計算出主因素的效應(yīng)和交互效應(yīng)，并且有效避免主因素與交互作用及交互作用間出現(xiàn)混淆（Confounded）現(xiàn)象。參數(shù)回歸模型能夠給出氣動力系數(shù)隨姿態(tài)角、馬赫數(shù)、控制面偏轉(zhuǎn)角等自變量變化的非線性規(guī)律，在風洞MDOE中被廣泛應(yīng)用［4］?；趨?shù)回歸模型的MDOE的形式實驗設(shè)計方法有中心復(fù)合設(shè)計（Center Composite Design，CCD）、面心復(fù)合設(shè)計（Face－centered Central Composite Design，F(xiàn)CD）和最優(yōu)回歸設(shè)計（Optimum Design）等［5－6］。這些設(shè)計方法都屬于回歸設(shè)計，在進行設(shè)計前要確定氣動模型的結(jié)構(gòu)形式。

隨著現(xiàn)代戰(zhàn)斗機和導彈對機動性要求的不斷提高，飛行器氣動力規(guī)律通常是非線性和多變量耦合的，基于參數(shù)回歸模型的MDOE方法已經(jīng)不能滿足要求。當氣動規(guī)律的非線性較強時，基于多項式的響應(yīng)面模型難以很好地對全局氣動規(guī)律進行擬合［7－8］。由于非參數(shù)模型較多項式模型有更強的非線性建模能力，所以基于非參數(shù)模型的MDOE方法逐漸被人們所關(guān)注［9］。非參數(shù)模型的形式實驗設(shè)計一般采用拉丁超立方設(shè)計（Latin Hypercube Sampling，LHS）和均勻設(shè)計（Uniform Design），在風洞MDOE中還沒有相關(guān)的研究。本文針對風洞非參數(shù)模型MDOE的形式實驗設(shè)計方法開展研究，采用Kriging方法進行建模，對拉丁超立方設(shè)計和均勻設(shè)計應(yīng)用于風洞實驗的適用性進行比較，并根據(jù)風洞實驗的特點發(fā)展相應(yīng)的優(yōu)化方法。

1 研究方法

為了深入分析研究不同的實驗設(shè)計方法的優(yōu)缺點，本文通過建立標模實驗氣動模型，進行“虛擬”風洞實驗?！疤摂M”實驗是實驗設(shè)計方法研究中常用的方法，實驗設(shè)計方法研究中需要頻繁地對各種設(shè)計進行比較和優(yōu)化，真實的風洞實驗周期長、成本高，而虛擬風洞實驗可以及時地獲得結(jié)果，采用這種方法能夠高效、低成本地進行實驗設(shè)計方法的研究［10］。本文所采用的實驗氣動模型是利用在中國航天十一院FD－06風洞已有的大量標模實驗數(shù)據(jù)構(gòu)建的。輸入一組新的設(shè)計樣本點，就可以獲得該樣本對應(yīng)的六分量氣動數(shù)據(jù)和不確定度，完成一次“虛擬”的風洞實驗過程。然后根據(jù)虛擬風洞實驗結(jié)果，進行建模和數(shù)據(jù)分析，通過和實驗氣動模型（即標準響應(yīng)面）相比較，考核其還原氣動模型的能力。

自變量的設(shè)計空間如下：

迎角：0°～40°；

滾轉(zhuǎn)角：0°～90°。

響應(yīng)變量：法向力系數(shù)CN，前體軸向力系數(shù)CAF，側(cè)向力系數(shù)CZ，俯仰力矩系數(shù)mzg，偏航力矩系數(shù)myg和滾轉(zhuǎn)力矩系數(shù)mx。

檢驗設(shè)計優(yōu)劣除了通過數(shù)據(jù)和繪圖直接對比外，還利用統(tǒng)計方法對各種設(shè)計方法建立的響應(yīng)面模型的近似精度進行比較。選擇一定量均布于設(shè)計空間的檢驗點，分別代入標準響應(yīng)面和樣本點建立的響應(yīng)面模型得到標準值yi、預(yù)測值y＾i和預(yù)測值的不確定度Ui，采用整體均方根誤差RMSE和判定系數(shù)R2對新響應(yīng)面模型的擬合程度進行評估，采用整體平均不確定度M＿U 對預(yù)測誤差進行評估，公式分別如下：

2 MDOE風洞實驗設(shè)計方法的選擇

拉丁超立方設(shè)計和均勻設(shè)計，二者也被稱為是“空間填充設(shè)計”［11］或者“空間抽樣設(shè)計”［12］。“沒有最好的設(shè)計，只有最適合的設(shè)計”［13］，本文根據(jù)風洞實驗的特點，對這兩種方法的適用性進行研究。

2.1 拉丁超立方設(shè)計

拉丁超立方法是應(yīng)用廣泛的現(xiàn)代實驗設(shè)計技術(shù)之一，取樣覆蓋均勻，適用于影響因素較多的情況，可顯著減少實驗規(guī)模［14］，實驗次數(shù)等于水平數(shù)。拉丁超立方抽樣是按如下方法選取設(shè)計點的：

設(shè)有n個設(shè)計變量，每個變量p個水平，拉丁超立方抽樣將每個設(shè)計變量分為p個區(qū)間，因均勻分布，p個區(qū)間等間隔。這樣整個變量空間分成pn個了區(qū)域。遵循下列兩個原則取實驗點：

1）樣本點在每個了區(qū)域隨機選??；

2）在任一維，即任一變量上的投影有p個區(qū)間，每個區(qū)間有且僅有一個樣本點。

由于拉丁超立方設(shè)計是一種隨機抽樣設(shè)計方法，容易產(chǎn)生分布不均勻的較差設(shè)計，所以通常加以一定的準則以保證其均勻性，常用的準則有最大最小距離（Maximize minimum distance）、最小最大距離（Minimize maximum distance）、最小差異準則（Minimize discrepancy）、相對理想累積分布最小均方根差異準則（Minimize RMSvariation from CDF）和相對理想累積分布最小最大差異準則（Minimize maximum variation from CDF）。本文對以上幾種準則條件下的LHS設(shè)計進行對比。設(shè)計樣本點為40個，通過“虛擬”風洞實驗獲得樣本數(shù)據(jù)，然后用樣本數(shù)據(jù)建立樣本點響應(yīng)面和標準響應(yīng)面進行比較。

對于模型的6個氣動力和力矩分量，需要建立6個響應(yīng)面。通?？v向氣動系數(shù)的響應(yīng)面模型非線性程度要比橫向的弱，如果樣本點能夠滿足橫向氣動系數(shù)響應(yīng)面模型的建立，縱向也能夠滿足。圖1～圖5給出了基于上述5種不同準則下的LHS樣本設(shè)計及其建立的滾轉(zhuǎn)力矩系數(shù)mx響應(yīng)面。表1給出了檢驗點的統(tǒng)計分析結(jié)果。

圖1 LHS樣本點及響應(yīng)面（最大最小距離準則）Fig.1 Sample points of LHS design and the response surface（Maximize minimum distance）

圖2 LHS樣本點及響應(yīng)面（最小最大距離準則）Fig.2 Sample points of LHS design and the response surface（Minimize maximum distance）

圖3 LHS樣本點及響應(yīng)面（最小差異準則）Fig.3 Sample points of LHS design and the response surface（Minimize discrepancy）

圖4 LHS樣本點及響應(yīng)面（相對理想累積分布最小均方根差異準則）Fig.4 Sample points of LHS design and the response surface（Minimize RMS variation from CDF）

圖5 LHS樣本點及響應(yīng)面（相對理想累積分布最小最大差異準則）Fig.5 Sample points of LHS design and the response surface（Minimize maximum variation from CDF）

表1 不同LHS設(shè)計方法的響應(yīng)面檢驗結(jié)果Table 1 Test results of the response surfaces of different LHS designs

從檢驗參數(shù)看，不同準則下的LHS設(shè)計的整體模型精度均較高，差異不明顯。從樣本點響應(yīng)面與標準響應(yīng)面對比圖可以看出，如果在某局部區(qū)域缺少樣本點，就容易出現(xiàn)局部預(yù)測誤差較大的問題，如滾轉(zhuǎn)角0°、迎角接近40°附近，預(yù)測響應(yīng)面的mx普遍偏大，而標準響應(yīng)面的mx在此處為小量。只有采用了相對理想累積分布最小最大差異準則的LHS樣本點的響應(yīng)面（圖5），其原因是該設(shè)計在此位置有一個樣本點。

所以，LHS設(shè)計盡管已經(jīng)采用了一些促進其“均勻化”的準則，但仍然容易出現(xiàn)某個局部沒有樣本點的情況，這會導致較大的預(yù)測偏差和不確定度，是風洞實驗難以接受的。其次，LHS的樣本點生成有隨機性，帶來兩個問題：一是每次設(shè)計的質(zhì)量會有差異，穩(wěn)定性不高；二是每次設(shè)計中樣本點分布不佳的區(qū)域和程度不同，不便于進行優(yōu)化。

2.2 均勻設(shè)計

均勻設(shè)計［15］由我國方開泰和王元院士提出，是一種空間填充（Space Filling）設(shè)計。均勻設(shè)計要求實驗點均勻的分布在實驗區(qū)域內(nèi)，均勻設(shè)計基于偏差來衡量樣本點在設(shè)計空間分布的均勻性，使得設(shè)計空間內(nèi)各區(qū)域能夠均等的“獲得”樣本點信息，每個樣本點所承載的信息也能夠在建模中被充分利用，這一優(yōu)點非常適合MDOE風洞實驗設(shè)計。

均勻性準則要求設(shè)計具有最好的均勻性，對于一個設(shè)計，若沒有其他設(shè)計比它有更小的偏差，則稱該設(shè)計為均勻設(shè)計［16］。為了衡量設(shè)計的均勻性，人們提出了偏差這個概念來衡量實驗點集在實驗區(qū)域內(nèi)散布均勻程度。不失一般性，假設(shè)實驗域為Cs＝［0，1］s。記P＝｛x1，…，xn｝＝（xij）n×s為Cs上的一個有n個實驗點的設(shè)計，它可表為n×s的矩陣，其中n表示實驗點數(shù)，s表示因子的個數(shù)。Hickernell［17］用泛函分析中再生核Hilbert空間的方法將古典的諸偏差表成統(tǒng)一形式，并且定義了許多新偏差，其中應(yīng)用得最廣泛的是如下的中心化L2－偏差CD（P）：

和可卷型L2－偏差WD（P）：

進行樣本點數(shù)量為40個的均勻設(shè)計，圖6（a）給出了一個設(shè)計結(jié)果，表2給出了樣本點的偏差計算和響應(yīng)面檢驗結(jié)果。均勻設(shè)計的樣本分布整體上非常均勻，沒有較明顯的“團聚”和“缺失”區(qū)域，所建立的響應(yīng)面檢驗結(jié)果也優(yōu)于LHS設(shè)計結(jié)果。從圖6（c）和6（d）可以看出，除了設(shè)計空間的邊緣位置外，樣本點響應(yīng)面和標準響應(yīng)面吻合良好，沒有差異或不確定度特別大的局部區(qū)域。

所以，對于風洞實驗均勻設(shè)計較LHS更為穩(wěn)健、均勻，更適合在MDOE風洞實驗中應(yīng)用。

表2 均勻設(shè)計樣本點的偏差計算和響應(yīng)面檢驗結(jié)果Table 2 Discrepancies of Uniform Design and test results of the response surface

圖6 均勻設(shè)計樣本點及響應(yīng)面Fig.6 Sample points of Uniform Design and the response surface

3 均勻設(shè)計應(yīng)用于MDOE風洞實驗的優(yōu)化與改進

均勻設(shè)計得到初始實驗設(shè)計結(jié)果還無法滿足風洞實驗的要求，還要根據(jù)風洞實驗的具體特點加以優(yōu)化。本文發(fā)展了幾種對均勻設(shè)計的優(yōu)化改進方法，可以進一步提高樣本點設(shè)計質(zhì)量，使其滿足風洞實驗的要求。

3.1 邊界點補充

對于飛行器氣動模型，邊界位置的數(shù)據(jù)質(zhì)量往往是非常重要的，而且當飛行器處于最大姿態(tài)角位置時，氣動規(guī)律出現(xiàn)非線性的可能性增大，要求數(shù)據(jù)盡可能內(nèi)插，以防止外插帶來的風險。邊界點的增加包括角點和邊線點，補充邊界點后的設(shè)計如圖7所示。補充邊界點時先將設(shè)計空間進行了標準化，即將各自變量范圍轉(zhuǎn)換到［0，1］區(qū)域內(nèi)，補充方法如下：

1）角點直接在設(shè)計空間角點增加，即［0，0］，［0，1］，［1，0］，［1，1］；

2）計算原始設(shè)計樣本間的最大最小距離Dmaxmin，即每個樣本點和其他點的最小距離中的最大值，本例中Dmaxmin為0.190 39，出現(xiàn)Dmaxmin的兩個點在圖7中用“☆”標記；

3）確定邊界影響區(qū)域范圍，這是給邊界提供建模信息的區(qū)域，通過判斷該區(qū)域點的分布來決定邊線點增加的位置。如本例取0.1作為劃分界限，距離四個邊界為0.1的點視為邊界影響區(qū)域樣本點，在圖7中用“＋”標出。

4）對邊界影響區(qū)域樣本點在其邊界方向的距離進行判斷，如果距離大于Dmaxmin，則在邊界以等分方式增加樣本點，直至邊界影響區(qū)域樣本點在其邊界方向的距離都小于Dmaxmin，結(jié)果如圖7所示。

圖7 均勻設(shè)計邊界點補充示意圖Fig.7 Boundary points filling of Uniform Design

3.2 樣本點密度調(diào)整

為了提高實驗數(shù)據(jù)的精確度，風洞實驗中常常在會發(fā)生氣動規(guī)律劇烈變化的自變量區(qū)域，或者研究者特別關(guān)注的實驗狀態(tài)附近，減小實驗點間隔，增加樣本點。均勻設(shè)計的樣本點在設(shè)計空間內(nèi)“均勻散布”，適合于研究中各區(qū)域的信息密度需求相同的情況。當有先驗信息時，在其基礎(chǔ)上改變樣本點分布密度，更符合風洞實驗的實際。圖8給出了迎角方向改變密度分布前后的對比，對于迎角0°～12°范圍的樣本點，迎角方向間隔增加到1.2倍，多出的點移入12°～40°的區(qū)間；迎角12°～40°范圍的樣本點，迎角方向間隔減小為原來的0.914 3倍。由于迎角大于12°后氣動規(guī)律隨迎角變化更為強烈，改變密度分布后迎角12°到40°的區(qū)間樣本點增加了2個，提高了大迎角區(qū)域的樣本點密度，更為合理。

圖8 樣本點密度調(diào)整前后對比Fig.8 Comparison the distributions of pre－and－after adjusting of the density of sample points

3.3 重復(fù)點設(shè)計

重復(fù)點數(shù)量由實驗的不確定度目標和實驗環(huán)境所決定，實驗環(huán)境誤差大，要獲得高質(zhì)量水平的數(shù)據(jù)，重復(fù)的次數(shù)需要多些。重復(fù)點的設(shè)計可以根據(jù)先驗來決定，如在規(guī)律變化劇烈或評價檢測風洞系統(tǒng)穩(wěn)定性的某些關(guān)鍵位置。除此之外，角點和邊線點位置也是設(shè)置重復(fù)點的重點區(qū)域。對設(shè)計空間的內(nèi)部點，重復(fù)點設(shè)置可以隨機選擇，也可以按照均勻設(shè)計方法進行選擇，圖9給出內(nèi)點重復(fù)設(shè)計的結(jié)果。

圖9 設(shè)計空間內(nèi)重復(fù)點的設(shè)計Fig.9 Design of the repetitive points in the design space

4 應(yīng)用示例

本節(jié)將均勻設(shè)計及其優(yōu)化方法應(yīng)用于基于Kriging模型MDOE實驗，設(shè)計變量為迎角和滾轉(zhuǎn)角。

4.1 實驗條件

風洞模型為導彈類標模，外形為舵身組合體，四片尾舵在0°滾轉(zhuǎn)角時呈“X”字布局。實驗馬赫數(shù)為1.79，自變量的設(shè)計空間為：迎角0°～40°；滾轉(zhuǎn)角0°～90°。

OFAT實驗的實驗測量點如圖10所示，共有108個測量點。表3給出了對六分量測力數(shù)據(jù)的實驗數(shù)據(jù)不確定度。

圖10 OFAT實驗的實驗測量點分布Fig.10 Distribution of the sample points in OFAT test

表3 實驗數(shù)據(jù)不確定度要求Table 3 Requirement of the uncertainty of the test data

4.2 形式試驗設(shè)計結(jié)果

初始樣本點為35點的均勻設(shè)計，增補了邊界點和重復(fù)點后最終設(shè)計為72個。樣本點組成為：初始均勻設(shè)計35個；角點4個；邊線點9個；重復(fù)點24個。

MDOE方法測量點要比OFAT方法的少36個點，測量點為OFAT方法的66.7%。為了增加大迎角下的測試點密度，迎角0°～12°范圍的樣本點在迎角方向間隔增加到1.2倍，最終設(shè)計如圖11所示。

圖11 MDOE方法最終樣本點設(shè)計Fig.11 Final design of the sample points in MDOE test

4.2.1 實驗結(jié)果

通過“虛擬”風洞實驗獲得72個樣本點的六分量氣動力和力矩系數(shù)后對其進行建模，獲得圖12所示的響應(yīng)面模型和不確定度評估?？梢钥闯?，利用均勻設(shè)計為基礎(chǔ)進行優(yōu)化所得到的樣本點設(shè)計，擬合出的氣動響應(yīng)面規(guī)律合理。盡管該模型的橫向氣動力和力矩系數(shù)表現(xiàn)出隨自變量強烈的非線性變化規(guī)律，樣本點對其非線性特征的采樣非常充分和準確。通過非參數(shù)建模的方式，能夠以非常直觀的形式表示出氣動特性隨迎角和滾轉(zhuǎn)角變化的規(guī)律，有利于分析和研究；而且以響應(yīng)面模型給出風洞實驗結(jié)果，任意自變量組合的值都可以直接獲得，具有信息量大的優(yōu)點。圖12右側(cè)一列給出了MDOE結(jié)果的不確定度，滿足實驗要求。

圖12 MDOE方法的六分量氣動力和力矩系數(shù)的響應(yīng)面及不確定度Fig.12 Response surfaces and uncertainty of aerodynamic parameters in MDOE test

5 結(jié) 論

本文對風洞非參數(shù)模型MDOE的形式實驗設(shè)計方法進行了研究，對拉丁超立方設(shè)計和均勻設(shè)計應(yīng)用于風洞實驗的適用性進行比較，并根據(jù)風洞實驗的特點發(fā)展幾種優(yōu)化方法，主要結(jié)論有：

（1）對于風洞實驗均勻設(shè)計較拉丁超立方設(shè)計更為穩(wěn)健、均勻，更適合基于非參數(shù)模型的風洞MDOE方法。

（2）在均勻設(shè)計方法基礎(chǔ)上，根據(jù)風洞實驗的特點發(fā)展了優(yōu)化的方法，包括邊界點補充、樣本點密度調(diào)整和重復(fù)點設(shè)計，能夠?qū)⒁延械摹跋闰炐畔ⅰ睉?yīng)用于實驗設(shè)計中。

（3）利用均勻設(shè)計為基礎(chǔ)進行優(yōu)化所得到的樣本點設(shè)計，所需的測量點要少于OFAT方法的（如示例中所用的測量點數(shù)僅為OFAT方法的66.7%），且能夠充分和準確地對較劇烈的非線性變化規(guī)律進行采樣。

［1］ DeLoach R.Applications of modern experiment design to wind tunnel testing at NASA Langley research center［R］.AIAA 98－0713，36th AIAA Aerospace Sciences Meeting and Exhibit，Reno，NV，Jan.1998.

［2］ DeLoach R，Micol J R.Comparison of resource requirements for a wind tunnel test designed with conventional vs.modern design of experiments methods［R］.AIAA 2011－1260，49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition，Orlando，F(xiàn)lorida，January 4－7，2011.

［3］ Box G E，Hunter W G，Hunter J S.Statistics for experimenters：an introduction to design，data analysis and model building［M］.1978.

［4］Drew Landman，Jim Simpson，Dan Vicroy，et al.Efficient methods for complex aircraft configuration aerodynamic characterization using response surface methodologies［J］.Journal of Aircraft，2007，44（4）：1189－1195.

［5］ 何為，薛衛(wèi)東，唐斌.優(yōu)化實驗設(shè)計方法及數(shù)據(jù)分析［M］.北京：化學工業(yè)出版社，2012.

［6］ Box G E，Hunter W G，Hunter J S.Statistics for experimenters：an introduction to design，data analysis and model building［M］.1978.

［7］ Richard DeLoach，Erickson Gary.Low－order response surface modeling of wind tunnel data over truncated inference subspaces［R］.AIAA 2003－0456，2003.

［8］ Erickson G E.Estimation of supersonic stage separation aerodynamics of winged－body launch vehicles using response surface methods［R］.NASA／TM－2010－216196，2010.

［9］ Ching F Lo，Zhao J L，DeLoach R.Application of neural networks to wind tunnel data response surface methods［R］.AIAA 2000－2639，2000.

［10］DeLoach R，Micol J R.Comparison of resource requirements for a wind tunnel test designed with conventional vs.modern design of experiments methods［R］.AIAA 2011－1260，49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition，Orlando，F(xiàn)lorida，January 4－7，2011.

［11］Kai－Tai Fang，Runze Li，Agus Sudjianto.Design and modeling for computer experiments［M］.London：Taylor ＆Francis Group，LLC.

［12］王勁峰，姜成晟，李連發(fā).空間抽樣與統(tǒng)計判斷［M］.北京：科學出版社，2009.

［13］方開泰.實驗設(shè)計與建模［M］.北京：高等教育出版社，2001.

［14］游海龍，賈新章，張小波，等.Kriging插值與拉丁超立方實驗相結(jié)合構(gòu)造電路元模型［J］.系統(tǒng)仿真學報，2005，17（11）：2752－2755.

［15］方開泰.均勻設(shè)計與均勻設(shè)計表［M］.北京：科學出版社，1994.

［16］Fang K T，Wang Y.Number－theoretic methods in statistics［M］.Chaprrzanand Hall，London，1994.

［17］Hickernell F J.A generalized discrepancy and quadrature error bound［J］.Mathematics of Computation，1998，67：299－332.

Formal design of experiments method of MDOE in wind tunnel tests

Zhang Jiang＊，Qin Yongming，Ma Handong

（China Academy of Aerospace Aerodynamics，Beijing 100074，China）

The existing MDOE methods based on the parametric model cannot meet the tunnel test requirements of advanced vehicle because of strong nonlinear aerodynamic behaviors.It is necessary to develop the MDOE method based on the non－parametric model instead.The formal design of experiments method of MDOE for wind tunnel tests is researched in this paper.The applicability of Latin Hypercube Sampling and Uniform Design in wind tunnel tests is compared through virtual experiments.And several optimizing methods based on Uniform Design are devoleped，which can improve the quality of design results to meet the requirement of wind tunnel tests.Studies show that，Uniform Design is more robust and uniform than Latin Hypercube Sampling for the MDOE method based on the non－parametric model，the optimization including boundary complement，density adjustment of sampling points and repetition points design can utilize the priori information to the design，the formal design of experiments method developed in the paper requires less sampling points than the OFAT，which are only 66.7%in the demonstration，and samples the intensive non－linear changing law adequately and exactly.The studies on the formal design of experiments method of MDOE in wind tunnel tests laid the foundation for further study of the MDOE method based on the non－parametric model.

MDOE；non－parametric model；Uniform Design；Latin Hypercube Sampling；wind tunnel test

V211.74

Adoi：10.7638／kqdlxxb－2015.0116

0258－1825（2016）01－0059－11

2015－07－21；

2015－10－25

張江＊（1978－），男，陜西綏德，高工，博士，研究方向：實驗流體力學.E－mail：13611319903＠163.com

張江，秦永明，馬漢東.風洞MDOE的形式實驗設(shè)計方法研究［J］.空氣動力學學報，2016，34（1）：59－69.

10.7638／kqdlxxb－2015.0116 Zhang J，Qin Y M，Ma H D.Formal design of experiments method of MDOE in wind tunnel tests［J］.Acta Aerodynamica Sinica，2016，34（1）：59－69.