孟祥吉 李彬彬

（中國(guó)直升機(jī)設(shè)計(jì)研究所，江西景德鎮(zhèn) 333000）

0.引言

復(fù)合材料吸能結(jié)構(gòu)的吸能問題是一個(gè)多個(gè)參數(shù)變量共同影響決定的復(fù)雜問題，影響參數(shù)既包括材料的剪切強(qiáng)度、斷裂韌性等材料參數(shù)，還包括了結(jié)構(gòu)的厚度、鋪層方向等結(jié)構(gòu)參數(shù)。此外，還包括了材料剪切與拉伸之間的作用系數(shù)、最大破壞應(yīng)變等有關(guān)損傷演化規(guī)律的模型參數(shù)。對(duì)于試驗(yàn)中需要改變的輸入?yún)?shù)，我們通常稱之為因子，把因子的不同的數(shù)值定義為因子的水平。對(duì)于試驗(yàn)中所得到的結(jié)果，我們稱之為響應(yīng)。一個(gè)試驗(yàn)可以有多個(gè)因子，也可以有多個(gè)響應(yīng)。在復(fù)合材料吸能結(jié)構(gòu)吸能性能研究中，吸能結(jié)構(gòu)吸收的總能量與吸能結(jié)構(gòu)單位重量吸收的能量（SEA）作為試驗(yàn)的2個(gè)響應(yīng)。

通常參數(shù)試驗(yàn)方法只能研究單個(gè)因子對(duì)響應(yīng)的影響，不能研究多個(gè)因子共同作用對(duì)響應(yīng)的影響。拉丁超立方抽樣方法就是為了解決多因子共同對(duì)響應(yīng)的作用規(guī)律，同時(shí)又能較正交設(shè)計(jì)明顯的減少試驗(yàn)點(diǎn)的重疊，有效地降低了樣本數(shù)量，這對(duì)復(fù)合材料吸能結(jié)構(gòu)沖擊模擬分析是非常重要的，大大減少了計(jì)算時(shí)間。并且由于吸能性能響應(yīng)與復(fù)合材料參數(shù)因子之間具有較大非線性關(guān)系，因此采用拉丁超立方設(shè)計(jì)能更有效地對(duì)非線性關(guān)系進(jìn)行擬合分析。

1.復(fù)合材料模型

本文采用文獻(xiàn)[1]、[2]、[3]中提出的復(fù)合材料模型，包括剛度退化模型、改進(jìn)的Hashin損傷判定準(zhǔn)則和指數(shù)型損傷演化規(guī)律。

1.1 剛度退化模型

在模型中，定義有效應(yīng)力σ～與名義應(yīng)力σ有如下關(guān)系：

其中M是損傷算子，有如下的形式：

df、dm、ds分別為對(duì)應(yīng)纖維、基體和剪切失效的損傷變量。

1.2 改進(jìn)的Hashin判定準(zhǔn)則

判定材料損傷須定義損傷判定準(zhǔn)則。對(duì)Hashin判定準(zhǔn)則進(jìn)行了改進(jìn)，使得材料的損傷判定準(zhǔn)則在橫向和縱向上一致，并且考慮了面內(nèi)剪切損傷判定準(zhǔn)則。

1.3 指數(shù)型損傷演化規(guī)律

改進(jìn)的指數(shù)型損傷演化規(guī)律為：

2.采用拉丁超立方試驗(yàn)的Isight分析模型

在20世紀(jì)80年代，隨著對(duì)多學(xué)科優(yōu)化設(shè)計(jì)的需求，提出了“計(jì)算機(jī)輔助設(shè)計(jì)優(yōu)化技術(shù)”。Isight是計(jì)算機(jī)輔助優(yōu)化軟件的典型代表，由唐兆成在參與通用電氣公司對(duì)波音777發(fā)動(dòng)機(jī)進(jìn)行優(yōu)化設(shè)計(jì)時(shí)開發(fā)出來，并在之后不斷進(jìn)行完善，逐漸發(fā)展成了廣泛應(yīng)用于航空、航天以及電子等領(lǐng)域的多學(xué)科優(yōu)化設(shè)計(jì)軟件。Isight具有非常好的自動(dòng)化和可視化，能夠自動(dòng)執(zhí)行仿真流程并進(jìn)行比較，其過程可以通過軟件圖表顯示出來。此外還能對(duì)計(jì)算的數(shù)據(jù)進(jìn)行分析以及后處理，幫助用戶做出正確的方案決策。Isight繼承了各種試驗(yàn)設(shè)計(jì)方法，能夠方便調(diào)用，非常適合作為復(fù)合材料吸能結(jié)構(gòu)、吸能性能研究的手段。

圖1為本文研究復(fù)合材料吸能結(jié)構(gòu)性能參數(shù)的Isight分析模型。模型為一個(gè)循環(huán)的流程，它主要由4個(gè)模塊構(gòu)成。第一部分為實(shí)驗(yàn)設(shè)計(jì)方法模塊，它決定著整個(gè)分析模型所采用的方法。本文采用的是拉丁超立方方法，需要在此模塊中定義需要的樣本數(shù)量，選取20個(gè)樣本數(shù)量。此外還需要定義所需研究的參數(shù)以及對(duì)應(yīng)該參數(shù)的水平。

圖1 拉丁超立方Isight分析模型

通過設(shè)置研究參數(shù)與樣本數(shù)量，可以隨機(jī)生成拉丁超立方式樣設(shè)計(jì)方案的設(shè)計(jì)矩陣。由于每次生成的設(shè)計(jì)矩陣都不一樣，可采用保存設(shè)計(jì)矩陣的方法，下次直接導(dǎo)入設(shè)計(jì)矩陣即可，設(shè)計(jì)矩陣如圖2所示。

圖2 Isight參數(shù)設(shè)計(jì)矩陣

第二部分為文件的讀寫組件，在該流程中主要負(fù)責(zé)將試驗(yàn)設(shè)計(jì)組件中定義的參數(shù)寫入模板文件中，生成后續(xù)程序的計(jì)算輸入文件。其原理過程如圖3所示，寫入?yún)?shù)的方式可以分為通過查找特征字符與指定行號(hào)來進(jìn)行，并且一個(gè)文件讀寫模塊可以讀寫多個(gè)文件，反之可以從文件中讀取參數(shù)數(shù)值。本文通過該組件將參數(shù)寫入了MSC.Dytran的計(jì)算文件與子程序文件中。

圖3 文件讀寫組件工作原理

第三部分為命令行組件，命令行組件主要用來啟動(dòng)外部的程序，執(zhí)行各種操作系統(tǒng)的shell腳本，本文這一部分主要包括3個(gè)OS Command組件。第一個(gè)OS Command組件用來執(zhí)行子程序的編譯，通過把編譯命令寫成compile.bat文件，再通過OS Command組件來執(zhí)行compile.bat文件。第二個(gè)OS Command組件的目的是用來執(zhí)行編譯好的程序來進(jìn)行有限元計(jì)算，在此組件中需要的文件為上一步中編譯好的程序與讀寫組件中生成的MSC.Dytran計(jì)算文件。第三個(gè)OS Command組件的目的是將MSC.Dytran的計(jì)算結(jié)果THS文件轉(zhuǎn)換為TXT文件格式，使得下一部分的Execl組件能夠讀取TXT文件中的計(jì)算結(jié)果。

第四部分為Excel組件部分，其主要作用為讀寫Excel中的單元值，并且執(zhí)行Excel中定義的宏。預(yù)先在模板Excel文件中輸入宏命令，該宏命令能將上一步生成的計(jì)算結(jié)果的txt文件導(dǎo)入Excel，并計(jì)算出平均載荷與最大載荷。Excel組件首先載入模板Excel，之后調(diào)用宏命令計(jì)算出結(jié)果，最后將結(jié)果讀取出來存儲(chǔ)到結(jié)果參數(shù)中。

復(fù)合材料性能參數(shù)與吸能性能評(píng)估參數(shù)：

評(píng)價(jià)復(fù)合材料結(jié)構(gòu)吸能性能的一個(gè)重要指標(biāo)是比吸能SEA[4-5]，其定義為：

即結(jié)構(gòu)有效破壞長(zhǎng)度內(nèi)單位質(zhì)量吸收的能量，其中：M為結(jié)構(gòu)壓潰破壞部分的質(zhì)量，F(xiàn)為接觸力，ρ為結(jié)構(gòu)平均密度，A為橫截面積，h為結(jié)構(gòu)有效壓潰位移。在穩(wěn)定壓潰模式下，平均壓潰載荷做功等于沖擊能量，由此得到平均壓潰載荷和比吸能有如下關(guān)系：

本節(jié)選取G827/5224復(fù)合材料圓管作為基礎(chǔ)材料，利用拉丁超立方的實(shí)驗(yàn)設(shè)計(jì)方法，通過改變材料的屬性得到復(fù)合材料圓管的吸能性能與材料屬性參數(shù)之間的關(guān)系，并研究材料參數(shù)之間的相互作用。在復(fù)合材料吸能結(jié)構(gòu)的軸向沖擊過程中，復(fù)合材料的剪切性能對(duì)其吸能性能具有較大的影響，因此將復(fù)合材料的剪切性能參數(shù)作為拉丁超立方試驗(yàn)的變量，以平均壓潰載荷作為目標(biāo)變量來衡量復(fù)合材料吸能結(jié)構(gòu)的吸能性能。

復(fù)合材料的剪切性能參數(shù)包括復(fù)合材料的剪切強(qiáng)度S（簡(jiǎn)稱強(qiáng)度S）、復(fù)合材料的剪切對(duì)拉伸的作用系數(shù)α（簡(jiǎn)稱剪切系數(shù)α）以及復(fù)合材料剪切損傷演化系數(shù)AM（簡(jiǎn)稱演化系數(shù)AM）。材料的剪切強(qiáng)度可以通過dat文件的MAT8A卡片進(jìn)行設(shè)置。根據(jù)一般復(fù)合材料的剪切性能假設(shè)其強(qiáng)度變化區(qū)間為60Mpa～150Mpa，而剪切對(duì)拉伸的作用系數(shù)α的取值范圍在0～1，復(fù)合材料剪切損傷演化系數(shù)AM的取值范圍則可根據(jù)材料性能與網(wǎng)格尺度進(jìn)行估算，一般材料取值區(qū)間為[0.02，0.2]。根據(jù)以上變量的取值范圍，通過拉丁超立方試驗(yàn)設(shè)計(jì)方法，選取10個(gè)試驗(yàn)樣本，樣本參數(shù)如表1所示。

表1 拉丁超立方實(shí)驗(yàn)設(shè)計(jì)樣本點(diǎn)參數(shù)

3.材料參數(shù)對(duì)復(fù)合材料吸能結(jié)構(gòu)軸向沖擊平均載荷的影響

復(fù)合材料吸能結(jié)構(gòu)的沖擊平均載荷一般決定著復(fù)合材料吸能結(jié)構(gòu)吸能能力的大小。通過對(duì)計(jì)算計(jì)算結(jié)果進(jìn)行分析可以得到復(fù)合材料的剪切強(qiáng)度S、復(fù)合材料的剪切對(duì)拉伸的作用系數(shù)α以及復(fù)合材料剪切損傷演化系數(shù)AM關(guān)于平均載荷的Pareto圖，如圖4所示，藍(lán)色表示平均載荷隨著變量的增加而增加，紅色表示平均載荷隨變量增加而減小。結(jié)果表明，剪切強(qiáng)度S對(duì)平均載荷的影響最大，其次為復(fù)合材料剪切損傷演化系數(shù)AM，剪切對(duì)拉伸的作用系數(shù)α對(duì)平均載荷的影響較小。

圖4 軸向沖擊平均載荷的Pareto圖

復(fù)合材料參數(shù)之間還存在相互作用，圖5材料剪切損傷演化系數(shù)AM與剪切對(duì)拉伸的作用系數(shù)α的交互作用圖，圖6材料剪切損傷演化系數(shù)AM與剪切強(qiáng)度S的交互作用圖。從圖中可以看出，2幅圖中的2條曲線基本保持趨勢(shì)一致，沒有發(fā)生交叉作用，表明2個(gè)參數(shù)間的交互作用較弱。

圖5 演化系數(shù)AM與剪切系數(shù)α交互作用圖

圖6 演化系數(shù)AM與強(qiáng)度S交互作用圖

圖7為剪切損傷剪切系數(shù)α與剪切強(qiáng)度S的交互作用圖。與前面2幅圖不同，圖中2條曲線趨勢(shì)相差較為明顯。在較高強(qiáng)度水平時(shí)，平均載荷隨著α的增加而逐漸減小，載荷變化區(qū)間較大。而在較低強(qiáng)度水平時(shí)，平均載荷隨著α的增加先逐漸減小再逐漸增大，載荷變化區(qū)間較小。結(jié)果表明剪切系數(shù)α與剪切強(qiáng)度S間存在較強(qiáng)的耦合效應(yīng)。進(jìn)行復(fù)合材料設(shè)計(jì)選材時(shí)，應(yīng)盡量選取強(qiáng)度較高、剪切損傷演化系數(shù)較小的材料。

圖7 剪切系數(shù)α與強(qiáng)度S交互作用圖

4.結(jié)論

（1）剪切強(qiáng)度S對(duì)平均載荷的影響最大，其次為復(fù)合材料剪切損傷演化系數(shù)AM，剪切對(duì)拉伸的作用系數(shù)α對(duì)平均載荷的影響較小。

（2）復(fù)合材料參數(shù)之間還存在相互作用，材料剪切損傷演化系數(shù)AM與剪切對(duì)拉伸的作用系數(shù)α的交互作用較弱，剪切系數(shù)α與剪切強(qiáng)度S之間存在較強(qiáng)的耦合效應(yīng)。

（3）進(jìn)行復(fù)合材料吸能結(jié)構(gòu)設(shè)計(jì)選材時(shí)應(yīng)盡量選取強(qiáng)度較高、剪切損傷演化系數(shù)較小的材料。