李 盼，原梅妮，趙興成，李 瑤

(中北大學(xué) 機電工程學(xué)院， 太原 030051)

【航空和航海工程】

基于正交設(shè)計、組合權(quán)重-灰色關(guān)聯(lián)的Ω夾層蒙皮抗鳥撞優(yōu)化設(shè)計

李 盼，原梅妮，趙興成，李 瑤

(中北大學(xué) 機電工程學(xué)院， 太原 030051)

為提高Ω形夾層蒙皮的抗鳥撞特性，作者以Ω形夾層蒙皮的Ω半徑、外蒙皮厚度、Ω夾層厚度為影響因素，以蒙皮吸能特性、外蒙皮最大位移、外蒙皮最大應(yīng)力為評價指標(biāo)，設(shè)計了三因素四水平的正交試驗。采用組合賦權(quán)法和灰色關(guān)聯(lián)法處理有限元試驗的數(shù)據(jù)，獲得了優(yōu)化設(shè)計的最優(yōu)參數(shù)組合。Ω半徑10 mm、外蒙皮厚度2 mm、Ω夾層厚度3 mm、內(nèi)蒙皮厚度4.22 mm組合為最優(yōu)方案。將最優(yōu)組合與4 mm單蒙皮對比，其蒙皮吸能特性提高24.1%；外蒙皮最大位移減少37.9%；外蒙皮最大應(yīng)力減少1.08%。。

鳥撞，Ω形夾層，鳥撞，數(shù)值模擬，正交設(shè)計法

自1903年飛機進(jìn)入人類世界以來，飛機的鳥撞事故不斷發(fā)生。據(jù)統(tǒng)計，近年來全世界每年大約發(fā)生1萬次“鳥撞”事故[1-2]，這對飛機和機載人員存在著很大的危害。因此抗鳥撞研究對飛機設(shè)計至關(guān)重要。

對于機翼的鳥撞問題，已有學(xué)者做了一些研究[3-6]。石霄鵬等[3]在蜂窩夾芯結(jié)構(gòu)后增加了支撐，有效提高了該結(jié)構(gòu)的抗鳥撞性能。王文智等[4-5]針對機翼前緣提出了內(nèi)嵌斜支板的結(jié)構(gòu)，該結(jié)構(gòu)通過將鳥體動能分散，提高了機翼的抗鳥撞性能。本文采用正交試驗法對Ω夾層蒙皮的抗鳥撞性能進(jìn)行優(yōu)化。正交試驗法對于安排多因素試驗效率極高，已在很多領(lǐng)域得以應(yīng)用[6-15]。文廣等[6-13]學(xué)者通過正交試驗法，對各自研究的問題進(jìn)行優(yōu)化分析，均得到了最優(yōu)的結(jié)構(gòu)組合，經(jīng)驗證結(jié)果良好。劉思峰[14]，劉世豪[15]，孫曉東[16]等將灰色關(guān)聯(lián)法與正交試驗法組合應(yīng)用，彌補了正交試驗法只能針對一個指標(biāo)優(yōu)化的不足，有效解決了多指標(biāo)優(yōu)化問題。Pahange H[17]通過正交設(shè)計法+灰色關(guān)聯(lián)法，得出了蒙皮厚度、肋厚度、肋間距的最優(yōu)組合，有效提高了蒙皮抗鳥撞性能。

基于上文研究和Ω夾層蒙皮設(shè)計，本文通過正交優(yōu)化設(shè)計+組合賦權(quán)法和灰色關(guān)聯(lián)法對其三因素四水平三指標(biāo)進(jìn)行優(yōu)化分析，推斷出Ω形夾層蒙皮的最優(yōu)結(jié)構(gòu)組合。并通過Ls-dyna軟件對最優(yōu)組合進(jìn)行驗證分析。

1 實驗方法

1.1 正交試驗的原理及特點

正交試驗是采取部分試驗來代替全面試驗,挑選出有代表性的試驗點來進(jìn)行試驗,通過對有代表性的試驗的結(jié)果分析,了解全面試驗的情況。因此,挑選有代表性的試驗點是正交試驗的關(guān)鍵。自1951年日本統(tǒng)計學(xué)家田口玄一根據(jù)試驗的優(yōu)化規(guī)律提出了正交表[18]，正交表成為正交試驗設(shè)計的基本工具,使得正交試驗具備了分散性和整齊可比性。用正交表安排多因素設(shè)計的方法,稱為正交設(shè)計法。它具有均衡分散性與整齊可比性,即在數(shù)學(xué)上稱為正交性,是正交表安排設(shè)計試驗最重要的特點之一。正交設(shè)計法對于安排多因素試驗效率高,往往能用較少的設(shè)計次數(shù)獲得或推斷出最佳設(shè)計結(jié)果,這也是Ω夾層蒙皮鳥撞優(yōu)化所需要的效果。

1.2 灰色關(guān)聯(lián)法

灰色關(guān)聯(lián)度分析是灰色系統(tǒng)的基本內(nèi)容，其基本思想是: 從原始數(shù)據(jù)矩陣中找出最佳決策方案，并根據(jù)線性比例變換法對評判矩陣進(jìn)行標(biāo)準(zhǔn)化，之后根據(jù)權(quán)重系數(shù)計算判斷矩陣[14-16]，獲得各個方案的投影值，以此作為評價方案優(yōu)劣的依據(jù)。因此本文通過組合賦權(quán)法和灰色關(guān)聯(lián)法對Ω夾層蒙皮三因素四水平三指標(biāo)進(jìn)行優(yōu)化分析，推斷出Ω形夾層蒙皮的最優(yōu)結(jié)構(gòu)組合。

1) 按照灰色關(guān)聯(lián)分析法，假設(shè)評價系統(tǒng)中有m個評價指標(biāo)，n個試驗方案，則評價指標(biāo)矩陣A為:

(1)

2) 對指標(biāo)矩陣歸一化

在系統(tǒng)的各評價指標(biāo)中，有的是越大越優(yōu)的指標(biāo)，有的是越小越優(yōu)的指標(biāo)，針對這一現(xiàn)象，對相應(yīng)的數(shù)據(jù)作區(qū)分處理:

對越大越優(yōu)的指標(biāo):

(2)

對越小越優(yōu)的指標(biāo)：

(3)

其中:i=1,2,…,n;j=1,2,…,m。

將各指標(biāo)中的最大值作為參考方案，

參考方案矩陣為:K=k1,k2…,km，

其中kj=max(r1j,r2j,…,rnj),j=1,2,…,m。

3) 關(guān)聯(lián)系數(shù)矩陣計算

(4)

式中: 分辨系數(shù)ρ∈0,1，本文取ρ=0.5。

1.3 組合賦權(quán)法

權(quán)重賦值的方法一般有兩種: 主觀賦權(quán)法和客觀賦權(quán)法。主觀賦權(quán)法的權(quán)重由主觀判斷得到，決策結(jié)果具有一定的主觀隨意性?？陀^賦權(quán)法的權(quán)重是依據(jù)原始數(shù)據(jù)的關(guān)系確定的，主要方法有均方差法，熵值法、離差法等[19]。

為了使最終的決策結(jié)果更加真實、可靠，從而實現(xiàn)Ω夾層蒙皮的優(yōu)化設(shè)計，本文結(jié)合了主觀賦權(quán)法和熵值法確定評估指標(biāo)的組合權(quán)重。其中主觀賦權(quán)法：蒙皮吸能量χ1=0.6；外蒙皮最大位移χ2=0.25；外蒙皮最大應(yīng)力χ3=0.15。

1) 熵值法

對指標(biāo)矩陣A進(jìn)行歸一化處理:

(5)

各個指標(biāo)的信息熵:

(6)

熵值法權(quán)重:

(7)

2) 組合權(quán)重

主觀權(quán)重和客觀權(quán)重進(jìn)行組合，組合方法按照文獻(xiàn)[11]中組合權(quán)重計算方法進(jìn)行處理，具體公式為:

(8)

3) 理想方案的關(guān)聯(lián)度計算

駱駝城井灌區(qū)位于高臺縣城西南25 km處，東與南華鎮(zhèn)接壤，西至黑泉鄉(xiāng)，南依戈壁荒漠，北至巷道鎮(zhèn)；轄1個鎮(zhèn)、13個行政村、113個農(nóng)場總?cè)丝谶_(dá)1.48萬人；總耕地面積0.793萬hm2，其中純機井灌溉面積0.667萬hm2，配套機電井728眼。灌區(qū)以提取地下水灌溉為主，引黑河水補充灌溉為輔,是一個水資源嚴(yán)重短缺的灌區(qū)。多年來由于地下水過度開采，致使地下水位逐年下降（年均下降0.67～1.19 m），灌區(qū)內(nèi)地下水形成了降落漏斗，雖然西干渠從2000年開始輸水，但對地下水位影響甚微，被確定為水資源管理紅區(qū)。

r=ξw

(9)

其中w為wj(j=(1,2,…m))組成的組合權(quán)重系數(shù)矩陣，r為目標(biāo)函數(shù)灰色關(guān)聯(lián)度矩陣。

1.4 鳥撞機翼模型

為優(yōu)化Ω形夾層蒙皮結(jié)構(gòu)以提高其抗鳥撞能力，本文采用LS-DYNA對Ω夾層蒙皮的鳥撞過程進(jìn)行模擬。鳥體采用SOLID164單元，其材料采用各向同性彈塑性材料模型。機翼采用SHELL163單元，其材料(2024-Al)采用JC材料模型，機翼四邊采用固接。材料參數(shù)見表1。Ω夾層蒙皮見圖1。

表1 材料參數(shù)

圖1 Ω夾層蒙皮簡圖

由圖1可知，Ω夾層蒙皮由外蒙皮，Ω夾層，內(nèi)蒙皮三部分組成，易知影響Ω夾層蒙皮抗鳥撞性能的因素有Ω半徑r，外蒙皮厚度h1，Ω厚度h2，內(nèi)蒙皮厚度h3。以蒙皮厚度為4 mm的單層蒙皮為基準(zhǔn)，鋁合金單層蒙皮結(jié)構(gòu)重量為13.2 kg。為保證Ω夾層蒙皮的結(jié)構(gòu)重量不超過單蒙皮。根據(jù)網(wǎng)格劃分單元數(shù)可得，在5 mm、10 mm、15 mm、20 mm四種Ω半徑下：外蒙皮面積最大為Ω形夾層的3.5倍，最大為內(nèi)蒙皮面積的3.7倍。由式(10)可知，當(dāng)外蒙皮厚度與Ω形夾層厚度確定時，內(nèi)蒙皮厚度也隨之確定。

由上所述，可得結(jié)構(gòu)重量等式：

2 780×4×蒙皮面積=2 780×h1×蒙皮面積+

(10)

化簡可得：

式中：h1為外蒙皮厚度；h2為Ω夾層蒙皮厚度；h3為內(nèi)蒙皮厚度。

2 Ω夾層蒙皮的優(yōu)化設(shè)計。

2.1 創(chuàng)建正交表

本文是在鳥速為370 km/h(飛機最大起飛速度)撞擊Ω形夾層蒙皮的情況上進(jìn)行優(yōu)化的。由于在鳥撞過程中由于在鳥撞過程中,蒙皮對于鳥體動能的吸收情況是影響蒙皮結(jié)構(gòu)的抗鳥撞性能的主要因素之一，因此把蒙皮的吸能特性最大作為主要優(yōu)化目標(biāo)，外蒙皮最大位移，外蒙皮最大應(yīng)力作為次要優(yōu)化目標(biāo)。在滿足此目標(biāo)的前提下, Ω形夾層蒙皮的結(jié)構(gòu)重量不能大于蒙皮厚度為4mm時的單蒙皮的結(jié)構(gòu)重量。

表2 Ω形夾層蒙皮正交優(yōu)化設(shè)計的影響因素與因素的水平 mm

將外蒙皮厚度h1、Ω夾層厚度h2代入式(10)得出內(nèi)蒙皮厚度h3，其中實驗4、8、12(見表3)中因外蒙皮厚度h1、Ω夾層厚度h2過大，超出蒙皮厚度為4 mm的單蒙皮結(jié)構(gòu)重量，為保證實驗完整性，在這三個實驗中取內(nèi)蒙皮厚度0.1 mm完成計算。計算結(jié)果見表3。試驗結(jié)果有鳥體回彈、鳥體穿透兩種情況。鳥體穿透時：蒙皮吸能為(鳥體初始動能-鳥體末時動能)；鳥體反彈時：蒙皮吸能為(鳥體初始動能+鳥體末時動能)。

Ω形夾層蒙皮正交優(yōu)化時的影響因素與影響水平數(shù)如表2所示.容易看出這是一個三因素四水平的問題,因此選取L16(43)有空列的正交表合適,此時n=16,即一共要做16次計算,r=4,m=3,d=16。

表3 Ω形夾層蒙皮正交優(yōu)化設(shè)計的計算結(jié)果

2.2 灰色關(guān)聯(lián)法+組合賦權(quán)法分析

將表3中Ω夾層蒙皮的3個評價指標(biāo)各仿真數(shù)據(jù)化為矩陣A。

根據(jù)式(2)、式(3)將矩陣A處理可得矩陣R；

根據(jù)式(4)可得矩陣ξ：

由式(5)、式(6)和式(7)，可得客觀權(quán)重為：

由上文可得主觀權(quán)重為：

由式(8)可得組合權(quán)重為：

wj=0.245 7, 0.151 2, 0.603 1

由式(9)計算的目標(biāo)函數(shù)的灰色關(guān)聯(lián)度，如表4所示。

表4 目標(biāo)函數(shù)關(guān)聯(lián)度

表5 各水平平均關(guān)聯(lián)度

其中各水平的平均關(guān)聯(lián)度，如表5所示。取各因素關(guān)聯(lián)度最大的，得最佳組合方案是：Ω半徑為10 mm；外蒙皮厚度2 mm；Ω夾層半徑3 mm。

3 Ω夾層蒙皮優(yōu)化設(shè)計驗證分析

由表6可知Ω夾層蒙皮的最優(yōu)組合為：Ω半徑為10 mm、外蒙皮厚度為2 mm、內(nèi)蒙皮厚度為3 mm。為了確定該組合就是最佳參數(shù)組合，對該方案進(jìn)行仿真試驗分析，與單蒙皮結(jié)果對比，如表6、圖2、圖3和圖4所示。

表6 三種指標(biāo)下兩種結(jié)構(gòu)對比

圖2 鳥體動能-時間曲線

圖3 外蒙皮最大位移-時間曲線

圖4 外蒙皮最大應(yīng)力-時間曲線

由圖2可知在370 km/h鳥速撞擊下，鋁合金單蒙皮的鳥體動能先下降然后穩(wěn)定在933 J(因鳥體穿透未出現(xiàn)回彈現(xiàn)象)。而Ω形夾層蒙皮鳥體動能下降至0 J，然后動能上升且穩(wěn)定在455 J，鳥體出現(xiàn)回彈現(xiàn)象。單蒙皮吸收了5 781 J的鳥體動能，Ω形夾層蒙皮吸收(鳥體初始動能+鳥體末時動能)7 175 J的鳥體動能。Ω夾層蒙皮可有效的防止鳥體穿透和吸收更多的鳥體動能，保護(hù)機翼內(nèi)部結(jié)構(gòu)。蒙皮吸能率相比于單蒙皮優(yōu)化率為24.1%。

由圖3可知在370 km/h鳥速撞擊下，鋁合金單蒙皮外蒙皮最大位移為104 mm，外蒙皮形變量較大。而Ω形夾層蒙皮外蒙皮最大位移為64.6 mm，外蒙皮形變量得到有效改善。相比于單蒙皮優(yōu)化率為37.9%。

由圖4可知在370 km/h鳥速撞擊下，鋁合金單蒙皮外蒙皮最大應(yīng)力為748 MPa。而Ω形夾層蒙皮外蒙皮最大位移為740 MPa。兩者差別不大，優(yōu)化效果不明顯。相比于單蒙皮優(yōu)化率為1.08%。但在6 ms后，由于Ω形夾層良好的抗沖擊性能，Ω形夾層蒙皮的外蒙皮最大應(yīng)力急速下降。

4 結(jié)論

1) 基于正交優(yōu)化設(shè)計方法，通過灰色關(guān)聯(lián)法與組合賦權(quán)法得出114結(jié)構(gòu)Ω形夾層蒙皮抗鳥撞能力最優(yōu)，其Ω形半徑為10 mm，外蒙皮厚度為2 mm，Ω形夾層厚度為3 mm，內(nèi)蒙皮厚度為4.22 mm。

2) 相對于4 mm單蒙皮，Ω形夾層蒙皮的蒙皮吸能提高24.1%；外蒙皮最大位移提高37.9%；外蒙皮最大應(yīng)力提高1.08%。Ω夾層蒙皮在外蒙皮最大位移的優(yōu)化較為明顯，蒙皮吸能特性次之，外蒙皮最大應(yīng)力優(yōu)化結(jié)果欠佳。

[1] 龔倫.飛機結(jié)構(gòu)的抗鳥撞設(shè)計與分析[D].西安：西北工業(yè)大學(xué), 2007.

[2] 李玉龍, 石霄鵬.民用飛機鳥撞研究現(xiàn)狀[J].航空學(xué)報, 2012, 33(2): 189-198.

[3] 石霄鵬, 李玉龍, 劉軍,等.某夾芯結(jié)構(gòu)抗鳥撞分析與設(shè)計[J].航空學(xué)報, 2012, 33(1):68-76.

[4] 王文智, 萬小朋, 郭葳.一種機翼前緣抗鳥撞結(jié)構(gòu)設(shè)計與分析[J].機械設(shè)計與制造, 2009(12):33-35.

[5] 楊立文, 韓全民, 張彥軍.隔板支持剛度對機翼前緣鳥撞性能的影響分析[J].力學(xué)與實踐, 2016, 38(1):33.

[6] 文廣, 周桂源, 何成剛,等.基于正交設(shè)計試驗的試驗機轉(zhuǎn)軸靜態(tài)靈敏度分析[J].組合機床與自動化加工技術(shù), 2015(4):52-54.

[7] 胡盛德, 郭志旺, 肖頤,等.利用正交試驗設(shè)計優(yōu)化冷彎成形工藝參數(shù)[J].軋鋼, 2015, 32(3):44-47.

[8] 叢小青, 周日, 韓泳濤,等.基于CFX正交試驗的深井離心泵導(dǎo)葉的優(yōu)化設(shè)計[J].流體機械, 2015(9):22-25.

[9] 王欣, 楊洪海, 張總輝,等.小型吸收式制冷熱虹吸泵的正交試驗設(shè)計[J].流體機械, 2015(7):63-66.

[10] 張旭, 吉衛(wèi)喜, 朱明清,等.基于虛擬正交試驗的TA2彎頭冷推制工藝參數(shù)優(yōu)化設(shè)計[J].鍛壓技術(shù), 2015, 40(2):46-51.

[11] SHUANG G, QIU B L, ZHU C Q, et al.Effects of comprehensive function of factors on retention behavior of microparticles in gravitational field-flow fractionation[J].Journal of Chromatography B, 2016, 1031:1-7.

[12] 王曉, 劉星榮, 葛如海.波紋管在汽車碰撞吸能中的正交優(yōu)化設(shè)計[J].江蘇大學(xué)學(xué)報(自然科學(xué)版), 2001, 22(3):29-32.

[13] XU X M, JIA H, ZHU Y P, et al.Optimization Design of Muffler Based on the Theory of Orthogonal Test[J].Applied Mechanics & Materials, 2014, 556-562:1191-1195.

[14] 劉思峰，蔡華，楊英杰，等.灰色關(guān)聯(lián)分析模型研究進(jìn)展[J].系統(tǒng)工程理論與實踐，2013，33(8): 2041 - 2044.

[15] 劉世豪，葉文華，陳蔚芳，等.基于正交試驗和灰色關(guān)聯(lián)的機床主軸箱優(yōu)化設(shè)計[J].振動與沖擊，2011，30(7):127-132.

[16] 孫曉東, 焦玥, 胡勁松.基于組合權(quán)重的灰色關(guān)聯(lián)理想解法及其應(yīng)用[J].工業(yè)工程與管理, 2006, 11(1):62-66.

[17] PAHANGE H, ABOLBASHARI M H.Mass and performance optimization of an airplane wing leading edge structure against bird strike using Taguchi-based grey relational analysis[J].中國航空學(xué)報(英文版), 2016, 29(4):934-944.

[18] 劉瑞江, 張業(yè)旺, 聞崇煒,等.正交試驗設(shè)計和分析方法研究[J].實驗技術(shù)與管理, 2010, 27(9):52-55.

[19] ZHENG C, KIZIL M S, CHEN Z, et al.Effects of coal properties on ventilation air leakage into methane gas drainage boreholes: Application of the orthogonal design[J].Journal of Natural Gas Science & Engineering, 2017, 45:88-95.

OptimizationDesignForbird-StrikeResistantofΩ-ShapedSandwichSkinBasedonOrthogonalExperimentalMethodandGreyRelationalAnalysis

LI Pan, YUAN Meini, ZHAO Xingcheng, LI Yao

(College of Mechanical and Electrical Engineering, North University of China, Taiyuan 030051, China)

In order to improve the bird-strike resistant of Ω-shaped sandwich skin, a series numerical experiments were performed using four-factor and three-level orthogonal method in this paper. The effect of radius of Ω-shaped skin, the thickness of outer skin and Ω-shaped middle skin on the energy absorption ability, the maximum displacement and stress of the outer skin were detailed investigated. The optimum combination of parameters were obtained from the simulations results which were analyzed with combination weighting and gray correlation method. The results showed that the sandwich structures with 2 mm-thickness of outer skin, 4.22 mm-thickness of inner skin, 3 mm thickness and 10 mm-radius Ω-shaped middle skin performs the best. Compared with the typical single layer skin with a thickness of 4 mm, the sandwich structures with the optimum combination of parameters improve the energy absorption ability by 24.1%, reduce the maximum displacement and stress of the outer skin by 37.9% and 1.08%, respectively.

bird-strike; Ω-shaped sandwich skin;numerical simulation；orthogonal experimental method

2017-08-20；

2017-09-12

國家自然科學(xué)基金(51201155)；山西省自然科學(xué)基金(2012011019-1)；教育部新教師基金(20101420120006)

李盼(1994—)，男，碩士研究生，主要從事飛行器沖擊動力學(xué)。

原梅妮(1974—)，女，博士，教授，主要從事復(fù)合材料力學(xué)性能研究，E-mail：mnyuan@126.com。

10.11809/scbgxb2017.12.052

本文引用格式：李盼,原梅妮，趙興成，等.基于正交設(shè)計、組合權(quán)重-灰色關(guān)聯(lián)的Ω夾層蒙皮抗鳥撞優(yōu)化設(shè)計[J].兵器裝備工程學(xué)報，2017(12):236-241.

formatLI Pan, YUAN Meini, ZHAO Xingcheng, et al.Optimization Design Forbird-Strike Resistant of Ω-Shaped Sandwich Skin Based on Orthogonal Experimental Method and Grey Relational Analysis[J].Journal of Ordnance Equipment Engineering,2017(12):236-241.

V215.2;TJ85

A

2096-2304(2017)12-0236-06

(責(zé)任編輯楊繼森)