郭 明， 柴 山, 劉金釗

(山東理工大學(xué) 交通與車輛工程學(xué)院， 山東 淄博 255091)

1 預(yù)備知識(shí)

結(jié)構(gòu)優(yōu)化設(shè)計(jì)問(wèn)題的數(shù)學(xué)模型一般形式為

findX=[x1,x2,…，xn]T∈Rn

minf(X)=f(x1,x2,…,xn)

s.t.σil≤[σ]i

(i=1,2,…,n；l=1,2,…,N)

(1)

?

其中X為設(shè)計(jì)變量；f(X)為目標(biāo)函數(shù)；σil及δlk分別為應(yīng)力約束條件和位移約束條件.

在一個(gè)結(jié)構(gòu)設(shè)計(jì)方案中,全部變量可分為三種類型,即設(shè)計(jì)變量、性能變量和中間變量.

1)設(shè)計(jì)變量是結(jié)構(gòu)優(yōu)化設(shè)計(jì)中的自變量，通常由設(shè)計(jì)者主動(dòng)選擇.

2)性能變量是結(jié)構(gòu)的各種性態(tài)變量,例如應(yīng)力、位移、自振頻率等等,是設(shè)計(jì)變量的因變量,設(shè)計(jì)者不能直接選出所需要的性態(tài)變量,而只能靠結(jié)構(gòu)分析來(lái)描述.

3)中間變量是由設(shè)計(jì)變量求性態(tài)變量運(yùn)算中的一些量.例如單元的應(yīng)力是一個(gè)性態(tài)變量,求應(yīng)力時(shí)所需的內(nèi)力就是一個(gè)中間變量.

一般說(shuō)來(lái)，在結(jié)構(gòu)優(yōu)化設(shè)計(jì)數(shù)學(xué)模型中，目標(biāo)函數(shù)是設(shè)計(jì)者關(guān)心的目標(biāo)與設(shè)計(jì)變量相關(guān)的函數(shù)，而約束條件是對(duì)性能變量的某些限制，是性能變量的函數(shù)[1].在有限元方程中，剛度矩陣、質(zhì)量矩陣、阻尼矩陣是設(shè)計(jì)變量的函數(shù)，而應(yīng)力、位移這些性能變量是要通過(guò)求解有限元方程求得的，因而是設(shè)計(jì)變量的隱函數(shù).一般的數(shù)學(xué)規(guī)劃問(wèn)題的目標(biāo)函數(shù)和約束條件都是設(shè)計(jì)變量的顯函數(shù)，結(jié)構(gòu)優(yōu)化設(shè)計(jì)數(shù)學(xué)模型中約束條件是設(shè)計(jì)變量的隱函數(shù)，這也是結(jié)構(gòu)優(yōu)化設(shè)計(jì)問(wèn)題區(qū)別于一般的數(shù)學(xué)規(guī)劃問(wèn)題的一個(gè)重要特點(diǎn)，由于這一特點(diǎn)，大大增加了結(jié)構(gòu)優(yōu)化設(shè)計(jì)問(wèn)題的復(fù)雜性，增加了求解難度.

要解決結(jié)構(gòu)優(yōu)化設(shè)計(jì)數(shù)學(xué)模型中約束條件是設(shè)計(jì)變量的隱函數(shù)的問(wèn)題，就必須進(jìn)行約束條件的顯式化，響應(yīng)面方法是最常用的顯式化方法.響應(yīng)面方法最初由Box和Wilson于1951年提出[2]，其基本思想是通過(guò)一系列確定性實(shí)驗(yàn)，用多項(xiàng)式函數(shù)來(lái)近似隱式極限狀態(tài)函數(shù).通過(guò)合理地選取試驗(yàn)點(diǎn)和迭代策略，來(lái)保證多項(xiàng)式函數(shù)能夠逼近于真實(shí)的隱式極限狀態(tài)函數(shù).當(dāng)系統(tǒng)變量和系統(tǒng)輸出響應(yīng)以某種隱含的方式存在時(shí)，RSM無(wú)疑提供了一種近似表達(dá)這種隱含關(guān)系的合適手段.RSM方法在統(tǒng)計(jì)數(shù)據(jù)處理、工業(yè)過(guò)程控制和可靠度分析等領(lǐng)域已為人們所熟知，近幾年在結(jié)構(gòu)優(yōu)化中也逐步得到了廣泛的應(yīng)用.1995年，Myers和Montgomery對(duì)響應(yīng)面的及其應(yīng)用進(jìn)行了全面的闡述，并把響應(yīng)面方法定義為“一種用于開(kāi)發(fā)、改進(jìn)、優(yōu)化的統(tǒng)計(jì)和數(shù)學(xué)方法”[3].如今，響應(yīng)面方法廣泛用于優(yōu)化設(shè)計(jì)中，即通過(guò)合理的試驗(yàn)設(shè)計(jì)方法構(gòu)建目標(biāo)函數(shù)、約束函數(shù)和設(shè)計(jì)變量之間的近似函數(shù)[4].

由于傳統(tǒng)的多項(xiàng)式基響應(yīng)面是應(yīng)用最小二乘法建立的逼近函數(shù)，在某些問(wèn)題中其擬合精度偏低.本文在結(jié)合桁架結(jié)構(gòu)優(yōu)化設(shè)計(jì)問(wèn)題，分析傳統(tǒng)多項(xiàng)式響應(yīng)面法及其擬合精度的基礎(chǔ)上，提出了基于泰勒展開(kāi)的響應(yīng)面法，可以有效地提高在展開(kāi)點(diǎn)附近響應(yīng)面的擬合精度，這對(duì)于結(jié)構(gòu)優(yōu)化設(shè)計(jì)非常有利.

2 傳統(tǒng)多項(xiàng)式基響應(yīng)面方法及其擬合精度分析

2.1 傳統(tǒng)多項(xiàng)式基響應(yīng)面算法簡(jiǎn)介

常用的多項(xiàng)式基響應(yīng)面法的形式有下面幾種:

線性型

(2)

二次型

(3)

完整二次型(含交叉項(xiàng))

(4)

(5)

每次試驗(yàn)的表達(dá)式可統(tǒng)一寫成如下矩陣形式

Y=Xβ+ε

(6)

其中

系數(shù)向量的無(wú)偏估計(jì)β可由最小二乘法獲得，即令每次試驗(yàn)的誤差平方和δ為最小

δ=εTε=(Y-Xβ)T(Y-Xβ)→min

(7)

β=(XTX)-1XTY

(8)

2.2 多項(xiàng)式基響應(yīng)面算法的精度分析

利用多項(xiàng)式基響應(yīng)面方法進(jìn)行顯式化時(shí)，響應(yīng)面的擬合的精度是要考慮的，文中用幾個(gè)簡(jiǎn)單的算例對(duì)多項(xiàng)式基響應(yīng)面方法的擬合精度進(jìn)行探討.本文主要討論的是桁架結(jié)構(gòu)，選取桁架結(jié)構(gòu)的截面積為設(shè)計(jì)變量，截面積變化范圍為0.01m2～0.05m2，且截面均為圓截面，彈性模量E=2.1×1011Pa，泊松比為0.3，以節(jié)點(diǎn)的位移為結(jié)構(gòu)的響應(yīng).

算例1 以最簡(jiǎn)單的二桿桁架為例,二桿桁架結(jié)構(gòu)如圖1所示.

圖1 二桿桁架結(jié)構(gòu)

通過(guò)結(jié)構(gòu)力學(xué)求得的二桿桁架的節(jié)點(diǎn)A沿力F方向的位移解析解如下：

利用多項(xiàng)式基響應(yīng)面求解的A節(jié)點(diǎn)的近似響應(yīng)值與真實(shí)響應(yīng)值之間的相對(duì)誤差和以及響應(yīng)面形式的選擇情況見(jiàn)表1.

表1 二桿桁架的擬合精度及響應(yīng)面形式的選擇 %

算例2 分析一個(gè)超靜定的六桿桁架的例子，六桿桁架的結(jié)構(gòu)如圖2所示[6].

圖2 六桿桁架結(jié)構(gòu)

六桿桁架節(jié)點(diǎn)B沿力F方向的位移的解析解為

利用多項(xiàng)式基響應(yīng)面求解B節(jié)點(diǎn)的近似響應(yīng)值與真實(shí)響應(yīng)值之間的相對(duì)誤差和以及響應(yīng)面形式的選擇情況見(jiàn)表2.

表2 六桿桁架的擬合精度及響應(yīng)面形式的選擇 %

由表1和表2可以看出，利用多項(xiàng)式基響應(yīng)面方法對(duì)桁架結(jié)構(gòu)的位移進(jìn)行擬合時(shí)，線性形式的響應(yīng)面在檢驗(yàn)點(diǎn)處相對(duì)誤都較大，不能夠滿足足夠的逼近；二次形式的響應(yīng)面在檢驗(yàn)點(diǎn)處計(jì)算精度雖然有所提高，但是誤差也很大；靜定結(jié)構(gòu)增加交叉項(xiàng)對(duì)擬合的結(jié)果影響很小，超靜定結(jié)構(gòu)增加交叉項(xiàng)后精度有所提高，但是也不能滿足實(shí)際需求.

3 基于泰勒展開(kāi)的響應(yīng)面方法

(9)

若線性形式的基于泰勒展開(kāi)的響應(yīng)面方程不能滿足精度的要求，則可增加高階形式的差分來(lái)進(jìn)行補(bǔ)充構(gòu)建響應(yīng)面函數(shù)，根據(jù)工程經(jīng)驗(yàn)和大量的資料分析可知，多元的二階泰勒展開(kāi)形式就能很好的滿足實(shí)際需要，基于多元泰勒展開(kāi)形式的二次響應(yīng)面公式.

(10)

其中：若式(10)中不含混合的二階偏導(dǎo)數(shù)，則稱該展開(kāi)為二次形式展開(kāi)，若式(10)中含有混合的二階偏導(dǎo)數(shù)，則稱為完整二次形式展開(kāi).

寫成矩陣形式為

S=S0+[S0]T(x-x0)+

其中H為海森矩陣

在差分網(wǎng)格圖3中，Si,j，Si+1,j，…,Si-1,j-1,Si+1,j-1，分別表示節(jié)點(diǎn)(i,j)，(i+1,j)，…，(i-1,j-1)，(i+1,j-1)處性能函數(shù)S的響應(yīng)值，若以節(jié)點(diǎn)(i,j)為展開(kāi)點(diǎn)，則上述各式中偏導(dǎo)數(shù)的求解公式如下：

圖3 差分的網(wǎng)格圖

一階偏導(dǎo)數(shù)：

(11)

二階偏導(dǎo)數(shù)：

(12)

混合的二階偏導(dǎo)數(shù)：

(Si+1,j+1+Si-1,j-1)]

(13)

對(duì)超靜定結(jié)構(gòu)利用基于泰勒展開(kāi)的響應(yīng)面進(jìn)行分析時(shí)，超靜定結(jié)構(gòu)的響應(yīng)S(內(nèi)力、應(yīng)力、位移、固有頻率、振型以及失穩(wěn)波形)等都是隨設(shè)計(jì)變量X而變化的，是設(shè)計(jì)變量的函數(shù)[7].將S在X0處進(jìn)行多元泰勒展開(kāi)，得到型如式(10)形式的多元泰勒展開(kāi)式.對(duì)于梁、桿、板、殼這類單元，其應(yīng)力可由內(nèi)力與截面幾何性質(zhì)的關(guān)系式求出，而截面幾何性質(zhì)就是設(shè)計(jì)變量，因此只要求得內(nèi)力，就可以得到應(yīng)力約束的顯函數(shù)表達(dá)式.

4 算例驗(yàn)證

為了驗(yàn)證基于泰勒展開(kāi)的響應(yīng)面方法的求解精

度，這里主要分析第一節(jié)中介紹的算例.

4.1 二桿靜定桁架

我們選擇在設(shè)計(jì)變量的中心點(diǎn)處進(jìn)行泰勒展開(kāi)得到如式(10)形式的基于泰勒展開(kāi)的顯式化的響應(yīng)面方程;對(duì)構(gòu)建的響應(yīng)面方程進(jìn)行精度檢驗(yàn)時(shí)，選擇在各設(shè)計(jì)變量的相鄰節(jié)點(diǎn)附近進(jìn)行精度的驗(yàn)算，利用基于泰勒展開(kāi)的響應(yīng)面法求解的結(jié)果與通過(guò)結(jié)構(gòu)的解析解求得結(jié)果的相對(duì)誤差見(jiàn)表3.

表3 二桿桁架的泰勒展開(kāi)形式及檢驗(yàn)點(diǎn)處的相對(duì)誤差%

由表3可以看出利用線性形式的泰勒展開(kāi)求解的精度較差，增加了二階泰勒展開(kāi)的響應(yīng)面求解的精度有所提高，增加了二次交叉項(xiàng)的泰勒展開(kāi)的響應(yīng)面驗(yàn)算點(diǎn)處求解的精度變化不大，平均相對(duì)誤差可以控制在2%左右，驗(yàn)算點(diǎn)處的精度能夠滿足工程實(shí)際的需要，可認(rèn)為求解結(jié)果是準(zhǔn)確的，同時(shí)與表1對(duì)比分析可知基于泰勒展開(kāi)的響應(yīng)面方法求解的精度比利用傳統(tǒng)多項(xiàng)式基響應(yīng)面方法求解的精度有了明顯的提高，因此可以構(gòu)建二次形式的基于泰勒展開(kāi)的響應(yīng)面方程，實(shí)現(xiàn)該結(jié)構(gòu)的隱函數(shù)的顯示化.

4.2 六桿超靜定桁架

利用基于泰勒展開(kāi)的響應(yīng)面法求解驗(yàn)算點(diǎn)的結(jié)果與通過(guò)結(jié)構(gòu)的解析解求得結(jié)果的相對(duì)誤差見(jiàn)表4.

表4 六桿桁架的泰勒展開(kāi)形式及檢驗(yàn)點(diǎn)處的相對(duì)誤差%

由表4可以看出來(lái)，增加了二階泰勒展開(kāi)的二次型和完整二次型的響應(yīng)面方法，求解的精度比線性展開(kāi)的形式有所提高，基本能夠保證工程實(shí)際的需要；同時(shí)與表2對(duì)比分析可知基于泰勒展開(kāi)的響應(yīng)面方法求解的精度比利用傳統(tǒng)多項(xiàng)式基響應(yīng)面方法求解的精度有了明顯的提高.因此可以構(gòu)建含交叉項(xiàng)的完整二次形式的基于泰勒展開(kāi)的響應(yīng)面方程，實(shí)現(xiàn)該結(jié)構(gòu)的隱函數(shù)的顯示化.

傳統(tǒng)的響應(yīng)面方法本質(zhì)上是在零點(diǎn)作泰勒展開(kāi)，而文中提出的方法是在設(shè)計(jì)變量的中間值做泰勒展開(kāi)，因此擬合精度會(huì)比傳統(tǒng)方法高；同時(shí)通過(guò)數(shù)據(jù)的變化趨勢(shì)分析可知越遠(yuǎn)離展開(kāi)點(diǎn)相對(duì)誤差越大，本文的方法是在一個(gè)較小的區(qū)域內(nèi)才能獲得較好的精度，此時(shí)可采用文獻(xiàn)[4]中介紹的運(yùn)動(dòng)極限的計(jì)算方法給出運(yùn)動(dòng)極限的粗估，并限制泰勒展開(kāi)的步長(zhǎng)來(lái)保證解的有效性和準(zhǔn)確性.

5 結(jié)論

1)通過(guò)桁架算例分析表明傳統(tǒng)多項(xiàng)式基響應(yīng)面在處理桁架問(wèn)題時(shí)擬合精度較差.

2)基于傳統(tǒng)多項(xiàng)式基響應(yīng)面算法存在的問(wèn)題，提出了基于泰勒展開(kāi)的響應(yīng)面方法,并通過(guò)區(qū)域網(wǎng)格劃分和差分公式的推導(dǎo)給出了該算法中各階偏導(dǎo)數(shù)的求解公式.

3)通過(guò)桁架結(jié)構(gòu)的算例分析，驗(yàn)證了本文提出的基于泰勒展開(kāi)的響應(yīng)面方法，能夠有效地提高桁架結(jié)構(gòu)的位移求解精度問(wèn)題，并通過(guò)運(yùn)動(dòng)極限的合理設(shè)置能夠高質(zhì)量的實(shí)現(xiàn)桁架結(jié)構(gòu)優(yōu)化設(shè)計(jì)中隱函數(shù)的顯式化.

[1] 孫煥純,柴山,王躍方,等.離散變量結(jié)構(gòu)優(yōu)化設(shè)計(jì)[M].大連：大連理工大學(xué)出版社，2001.

[2] Box G E P,Wilson K B.On the experimental attainment of optimum conditions (with discussion)[J]. Journal of Royal Statistical Society,1951,13(1):1-45.

[3] Myers R H,Montgomery D C. Response surface methodology: process and product optimization using designed experiments, [M].New York:John Wiley&Sons,Inc,1995.

[4] 隋允康,宇慧平.響應(yīng)面方法的改進(jìn)及其對(duì)工程優(yōu)化的應(yīng)用[M].北京：科學(xué)出版社，2011.

[5] 隋允康,張立新,杜家政.基于響應(yīng)面方法的桁架截面敏度分析和優(yōu)化[J].力學(xué)季刊,2006,27(1):96-102.

[6] 李廉錕.結(jié)構(gòu)力學(xué)[M].北京：高等教育出版社.2004.

[7] 錢令希.工程結(jié)構(gòu)優(yōu)化設(shè)計(jì)[M].北京:水利電力出版社.1983.