Bézier單元?jiǎng)偠扔成湎碌母咝Ф嘀鼐W(wǎng)格等幾何拓?fù)鋬?yōu)化方法

丁延冬 羅年猛 楊?yuàn)W迪 王書(shū)亭 朱浩然 謝賢達(dá)

華中科技大學(xué)機(jī)械科學(xué)與工程學(xué)院，武漢，430074

0 引言

拓?fù)鋬?yōu)化是一種在給定條件約束的設(shè)計(jì)域進(jìn)行結(jié)構(gòu)設(shè)計(jì)的有效工具，被廣泛應(yīng)用于散熱結(jié)構(gòu)[1]、柔性機(jī)構(gòu)[2]與彈性結(jié)構(gòu)[3]。經(jīng)典拓?fù)鋬?yōu)化方法主要有各向同性懲罰微結(jié)構(gòu)法(simple isotropic material with penalization，SIMP)[4]、進(jìn)化法[5]、水平集法[6-7]、移動(dòng)可變形組件/孔洞法[8-9]等。有限元法是傳統(tǒng)拓?fù)鋬?yōu)化的求解基礎(chǔ)，其求解精度對(duì)拓?fù)鋬?yōu)化過(guò)程的魯棒性至關(guān)重要，而有限元物理場(chǎng)離散形函數(shù)的C0連續(xù)無(wú)法滿足應(yīng)力約束等高階拓?fù)鋬?yōu)化問(wèn)題的求解要求。

無(wú)論是傳統(tǒng)有限元拓?fù)鋬?yōu)化程序還是等幾何拓?fù)鋬?yōu)化程序，更精確的優(yōu)化結(jié)果均需要規(guī)模更大的有限元求解，因此優(yōu)化過(guò)程需消耗大量的計(jì)算成本[17]。通過(guò)改進(jìn)有限元求解算法來(lái)提高拓?fù)鋬?yōu)化計(jì)算效率一直是拓?fù)鋬?yōu)化的研究熱點(diǎn)，主要方法有自適應(yīng)網(wǎng)格法[18-19]和多重網(wǎng)格法[20-22]。WANG等[23]將等幾何分析求解精度優(yōu)勢(shì)與多重網(wǎng)格求解效率的優(yōu)勢(shì)相結(jié)合，提出了多重網(wǎng)格等幾何優(yōu)化模型。然而，由于B樣條基函數(shù)空間的非單元一致性，上述多重網(wǎng)格等幾何拓?fù)鋬?yōu)化方法在每一層的拓?fù)鋬?yōu)化過(guò)程中需要進(jìn)行所有B樣條單元?jiǎng)偠鹊挠?jì)算和存儲(chǔ)，進(jìn)而實(shí)現(xiàn)優(yōu)化過(guò)程中等幾何單元?jiǎng)偠染仃嚨目焖俑拢瑢?dǎo)致多重等幾何拓?fù)鋬?yōu)化存在數(shù)據(jù)存儲(chǔ)負(fù)擔(dān)過(guò)重和預(yù)處理時(shí)間過(guò)長(zhǎng)等問(wèn)題。針對(duì)上述問(wèn)題，可采用Bézier提取算子對(duì)B樣條單元數(shù)據(jù)進(jìn)行預(yù)處理，實(shí)現(xiàn)任一層級(jí)網(wǎng)格B樣條基函數(shù)通過(guò)C0連續(xù)且單元空間保持的Bernstein基函數(shù)來(lái)等效表達(dá)。其中，Bézier提取算子作為一種基本的轉(zhuǎn)換方式，被用于保證B樣條基函數(shù)的高階連續(xù)性，已被廣泛應(yīng)用于等幾何分析框架中的T樣條[24]、截?cái)鄬哟蜝樣條[25-27]與極樣條[28]。

本文提出基于Bézier單元?jiǎng)偠扔成涞亩嘀鼐W(wǎng)格等幾何拓?fù)鋬?yōu)化模型，以實(shí)現(xiàn)等幾何拓?fù)鋬?yōu)化數(shù)據(jù)的存儲(chǔ)結(jié)構(gòu)和預(yù)處理過(guò)程效率的優(yōu)化。通過(guò)Bézier提取算子將各層級(jí)的所有B樣條單元?jiǎng)偠染仃嚨拇笠?guī)模存儲(chǔ)轉(zhuǎn)換為標(biāo)準(zhǔn)剛度矩陣與各層級(jí)Bézier提取矩陣的存儲(chǔ)，并將任一層級(jí)所有B樣條單元?jiǎng)偠染仃嚨挠?jì)算過(guò)程轉(zhuǎn)換為任一層級(jí)若干個(gè)(2或3)參數(shù)方向的Bézier提取矩陣的計(jì)算過(guò)程，進(jìn)而建立任一層級(jí)任一B樣條單元?jiǎng)偠染仃嚺c標(biāo)準(zhǔn)Bézier單元?jiǎng)偠染仃嚨挠成潢P(guān)系。利用B樣條單元?jiǎng)偠染仃嚺c標(biāo)準(zhǔn)Bézier單元?jiǎng)偠染仃囍g的映射法則，顯著降低多重網(wǎng)格等幾何拓?fù)鋬?yōu)化的剛度矩陣內(nèi)存占用并減少剛度矩陣預(yù)處理流程的時(shí)間消耗。最后，通過(guò)二維與三維數(shù)值算例的優(yōu)化結(jié)果展示本文模型在內(nèi)存空間和預(yù)處理時(shí)間兩個(gè)方面的優(yōu)勢(shì)。

1 理論基礎(chǔ)

首先介紹Bernstein多項(xiàng)式，接著對(duì)開(kāi)節(jié)點(diǎn)向量生成B樣條的遞歸定義進(jìn)行回顧，最后基于Bézier分解算法引出Bézier提取算子。其中，Bézier提取算子是實(shí)現(xiàn)Bernstein多項(xiàng)式與B樣條之間映射的理論基礎(chǔ)。

1.1 Bernstein多項(xiàng)式

對(duì)于任意位于[0，1]內(nèi)的參數(shù)坐標(biāo)ζ，相應(yīng)的p階Bernstein多項(xiàng)式由以下遞推式獲得：

Bi，p(ζ)=

(1)

其中，Bi，0(ζ)=1。階次p=1，2，3的Bernstein基函數(shù)如圖1所示。

(a)p=1 (b)p=2 (c)p=3

Bézier曲線可以由Bernstein多項(xiàng)式作為基函數(shù)并與控制點(diǎn)Pi的線性組合來(lái)表達(dá)，表達(dá)式為

(2)

1.2 B樣條

對(duì)于給定的開(kāi)節(jié)點(diǎn)向量Ξ=(ζ1，ζ2，…，ζn+p+1)，B樣條的基函數(shù)由Cox-de Boor遞推公式定義：

Nj，p(ζ)=

(3)

其中，p、n分別表示B樣條基函數(shù)的階次和相應(yīng)控制點(diǎn)個(gè)數(shù)。通過(guò)組合B樣條基函數(shù)與控制點(diǎn)可生成B樣條曲線，B樣條曲線的表達(dá)式如下：

(4)

多維參數(shù)空間下，B樣條基函數(shù)由張量積結(jié)構(gòu)生成，因此，二維B樣條基函數(shù)由以下公式定義：

N(j1，j2)，(p1，p2)(ζ，η)=Nj1，p1(ζ)Nj2，p2(η)

(5)

其中，(ζ，η)表示二維參數(shù)空間下的參數(shù)坐標(biāo)，jd(d=1，2)、pd(d=1，2)分別表示d維參數(shù)方向基函數(shù)的索引和階次。開(kāi)節(jié)點(diǎn)向量(0,0,0,1,2,3,3,3)的B樣條基函數(shù)空間和各節(jié)點(diǎn)向量單元的局部B樣條基函數(shù)空間如圖2所示，圖3為圖2所示B樣條基函數(shù)生成的B樣條曲線。

圖2 開(kāi)節(jié)點(diǎn)向量(0，0，0，1，2，3，3，3)的2階B樣條基函數(shù)

圖3 圖2中B樣條基函數(shù)所生成的B樣條曲線

1.3 Bézier提取算子

(6)

其中，αj的推導(dǎo)公式為

(7)

C(k)=

(8)

C(k)∈R(n+k-1)×(n+k)

整體的Bézier提取矩陣由C(k)(k=1，2，…，m)組成，其表達(dá)式為

Cext=C(1)C(2)…C(m)

(9)

N(ζ)=CextB(ζ)

(10)

Bézier曲線的控制點(diǎn)PB也可以通過(guò)Bézier提取算子從相應(yīng)B樣條控制點(diǎn)P獲取，如下式：

(11)

其中，PB的維度為(n+m)×d；P的維度為n×d。

借助于B樣條的張量積結(jié)構(gòu)，多變量的Bézier提取算子可方便地從單變量Bézier提取算子生成。雙變量情況下，Bézier提取算子的表達(dá)式如下：

C(ζ，η)=Cext(η)?Cext(ζ)

(12)

階次p=2時(shí)，圖2中任意單元的局部B樣條基函數(shù)均可映射為圖1b中具有空間保持特點(diǎn)的Bernstein基函數(shù)。圖4所示為圖2中單元3的B樣條基函數(shù)空間與Bernstein基函數(shù)空間之間的映射過(guò)程。圖5所示為B樣條基函數(shù)空間與相應(yīng)的Bernstein基函數(shù)空間，具體的映射關(guān)系表達(dá)如下：

圖4 圖2中單元3的基函數(shù)與圖1b的Bernstein基函數(shù)的映射關(guān)系

圖5 圖2中B樣條基函數(shù)與Bernstein基函數(shù)映射關(guān)系

(13)

2 多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型

首先介紹B樣條空間下等幾何拓?fù)鋬?yōu)化模型，在此基礎(chǔ)上對(duì)傳統(tǒng)多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型[23]進(jìn)行回顧；接著，基于Bézier提取算子對(duì)傳統(tǒng)等幾何拓?fù)鋬?yōu)化模型進(jìn)行重構(gòu)，獲得改進(jìn)多重網(wǎng)格模型；最后，通過(guò)多重網(wǎng)格各層設(shè)計(jì)變量之間的映射關(guān)系，完成改進(jìn)多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型的靈敏度分析。

2.1 傳統(tǒng)多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型

給定圖6所示的二維矩形設(shè)計(jì)域，等幾何拓?fù)鋬?yōu)化的設(shè)計(jì)目標(biāo)是在約束條件下尋找最優(yōu)的材料分布，即尋找具有最小應(yīng)變能的最優(yōu)結(jié)構(gòu)，目標(biāo)函數(shù)表示為

圖6 B樣條空間下的6×2網(wǎng)格的二維懸臂梁設(shè)計(jì)域

c(x)=fTu(x)

(14)

式中，x為所有B樣條單元的相對(duì)材料密度分布向量；c(x)為在材料分布x下的模型結(jié)構(gòu)柔度；f為外部力向量；u(x)為各個(gè)控制點(diǎn)上的位移向量。

綜合SIMP方法，B樣條空間下等幾何拓?fù)鋬?yōu)化的數(shù)學(xué)模型可以表示為

(15)

式中，ne為設(shè)計(jì)域離散成B樣條網(wǎng)格的單元個(gè)數(shù)；xe、ue、Ke(xe)分別為第e個(gè)單元的相對(duì)密度、位移向量和剛度矩陣；K(x)為全局剛度矩陣；V(x)為固體材料的體積；V0為設(shè)計(jì)域的總體積；vf為模型限制的材料體積分?jǐn)?shù)；為x的可容許空間。

式(15)中的Ke(xe)可由下式獲得：

(16)

Ee(xe)=Emin+(xe)p(E0-Emin)xe∈[0，1]

(17)

式中，E0為實(shí)體材料的彈性模量；Emin為空白單元的彈性模量，通常取0

基于多重網(wǎng)格的等幾何拓?fù)鋬?yōu)化模型是一種平衡優(yōu)化效率與優(yōu)化精度的有效手段，具體過(guò)程為：首先，等幾何拓?fù)鋬?yōu)化在較粗的網(wǎng)格下進(jìn)行優(yōu)化，此時(shí)迭代的速度較快、計(jì)算成本較低；然后，當(dāng)設(shè)計(jì)變量變化較小時(shí)，將粗網(wǎng)格優(yōu)化結(jié)果映射為更細(xì)網(wǎng)格的初始設(shè)計(jì)，并在細(xì)網(wǎng)格上進(jìn)行等幾何分析和優(yōu)化；最后，重復(fù)以上操作若干次，由最精細(xì)的網(wǎng)格生成等幾何拓?fù)鋬?yōu)化結(jié)果。

多重網(wǎng)格優(yōu)化框架如圖7所示，其中，等幾何拓?fù)鋬?yōu)化部分是基礎(chǔ)部分，是迭代過(guò)程中的迭代體；設(shè)計(jì)變量映射繼承部分是算法實(shí)現(xiàn)的核心，是迭代過(guò)程中完成代際繼承的關(guān)鍵。本文采用的設(shè)計(jì)變量映射方法為多維度擴(kuò)展方法，具體實(shí)現(xiàn)過(guò)程為：通過(guò)Kronecker乘積方法將粗網(wǎng)格模型的優(yōu)化設(shè)計(jì)變量以1對(duì)4(二維)或8(三維)的方式等值映射至細(xì)網(wǎng)格模型，進(jìn)而實(shí)現(xiàn)細(xì)網(wǎng)格模型相應(yīng)設(shè)計(jì)變量的初始化，如圖8所示。通過(guò)這樣的分層迭代過(guò)程與代際映射繼承的方式，在保證優(yōu)化精度的前提下顯著提高優(yōu)化求解的效率。

圖7 多重網(wǎng)格等幾何拓?fù)鋬?yōu)化的算法流程圖

(a)二維多重網(wǎng)格模型映射示意圖

將多重網(wǎng)格的框架與式(15)所示的數(shù)學(xué)模型相結(jié)合，得到多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型：

(18)

式中，l為多重網(wǎng)格模型的網(wǎng)格層級(jí)，l=1，2，…，L；xl的元素取值區(qū)間為[0，1]。

2.2 基于Bézier提取算子的多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型

首先，基于Bézier提取算子對(duì)傳統(tǒng)等幾何拓?fù)鋬?yōu)化模型進(jìn)行重構(gòu)；接著，將改進(jìn)等幾何拓?fù)鋬?yōu)化模型與圖7所示的多重網(wǎng)格的框架結(jié)合，提出B樣條單元?jiǎng)偠染仃嚱y(tǒng)一表達(dá)的多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型。

(19)

(20)

式(19)中，Be，Ber是關(guān)于Bernstein基函數(shù)的偏導(dǎo)數(shù)矩陣。根據(jù)圖5可知，任一B樣條單元的基函數(shù)均可使用相同的Bernstein基函數(shù)空間進(jìn)行等效表達(dá)，即各個(gè)單元的Be，Ber均相等。

根據(jù)上述特性及式(20)，可推導(dǎo)得到圖9所示的不同Bézier單元?jiǎng)偠染仃囍g的關(guān)系：

圖9 Bernstein空間下6×2網(wǎng)格的二維懸臂梁設(shè)計(jì)域

(21)

通過(guò)Bézier提取矩陣，圖9所示的B樣條單元?jiǎng)偠染仃嚤磉_(dá)式如下：

(22)

根據(jù)等幾何分析的等參特性，除了幾何空間可使用Bézier提取算子進(jìn)行映射之外，B樣條和Bézier控制點(diǎn)的位移向量也可使用Bézier提取算子實(shí)現(xiàn)二者之間的轉(zhuǎn)換。B樣條控制點(diǎn)位移u到Bézier控制點(diǎn)的位移ub的轉(zhuǎn)換方程如下：

(23)

基于以上的分析推導(dǎo)，對(duì)式(15)多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型進(jìn)行重構(gòu)，可得

(24)

相對(duì)于傳統(tǒng)多重網(wǎng)格模型，本文提出的多重網(wǎng)格等幾何模型的單元數(shù)據(jù)存儲(chǔ)結(jié)構(gòu)優(yōu)化如圖10所示，其中，L、n分別表示多重網(wǎng)格模型的層數(shù)與等幾何分析域參數(shù)空間的維度，KE0為各層統(tǒng)一的標(biāo)準(zhǔn)單元?jiǎng)偠染仃嚒鹘y(tǒng)的多重網(wǎng)格模型每一層都需要進(jìn)行所有B樣條單元?jiǎng)偠染仃嚨挠?jì)算和存儲(chǔ)，而本文多重網(wǎng)格模型采用標(biāo)準(zhǔn)Bézier單元?jiǎng)偠染仃嚺c相應(yīng)的Bézier提取算子實(shí)現(xiàn)任意層級(jí)的任一B樣條單元?jiǎng)偠染仃嚨慕y(tǒng)一表達(dá)，大幅減少了多重網(wǎng)格等幾何拓?fù)鋬?yōu)化的預(yù)處理時(shí)耗與單元數(shù)據(jù)內(nèi)存的消耗。

圖10 傳統(tǒng)與本文多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型的單元?jiǎng)偠染仃噧?nèi)存結(jié)構(gòu)對(duì)比

2.3 靈敏度分析

根據(jù)式(24)的數(shù)學(xué)模型以及圖7所示的算法流程，目標(biāo)函數(shù)相對(duì)于第l層的任意一個(gè)設(shè)計(jì)變量xe,l的靈敏度表達(dá)如下：

(25)

為避免數(shù)值不穩(wěn)定性現(xiàn)象，本文采用靈敏度過(guò)濾器修正目標(biāo)函數(shù)的靈敏度公式：

(26)

(27)

式中，wa,l為第l層的第e個(gè)單元與第a個(gè)單元之間的權(quán)重因子；dist(e，a)表示第a個(gè)單元與第e個(gè)單元的質(zhì)心之間的歐氏距離；rmin為由用戶定義的過(guò)濾器半徑。

3 數(shù)值算例

本文算例中的程序均運(yùn)行于Windows 10 系統(tǒng)的MATLAB R2021a，所用計(jì)算機(jī)的配置為Intel(R)Xeon(R)Gold 5102@2.20 GHz的CPU與64 GB內(nèi)存。圖11所示為激光雷達(dá)測(cè)試平臺(tái)的支撐結(jié)構(gòu)。對(duì)于激光雷達(dá)的主支撐結(jié)構(gòu)，可使用拓?fù)鋬?yōu)化方法進(jìn)行優(yōu)化設(shè)計(jì)，主支撐厚度較小時(shí)可用二維拓?fù)鋬?yōu)化以簡(jiǎn)化設(shè)計(jì)與制造流程，主支撐厚度較大時(shí)可使用三維拓?fù)鋬?yōu)化。二維和三維的激光雷達(dá)主支撐結(jié)構(gòu)設(shè)計(jì)模型如圖12所示。

圖11 激光雷達(dá)測(cè)試平臺(tái)

圖12a中二維主支撐結(jié)構(gòu)受到右側(cè)上端向下的外力負(fù)載，分別在該二維主支撐結(jié)構(gòu)的計(jì)算條件下應(yīng)用傳統(tǒng)多重網(wǎng)格模型[23]與基于Bézier提取的改進(jìn)多重網(wǎng)格模型。其中，等幾何網(wǎng)格的初始尺寸為40×20，多重網(wǎng)格的最大層數(shù)為4。二維主支撐結(jié)構(gòu)的多重網(wǎng)格等幾何拓?fù)鋬?yōu)化結(jié)果和優(yōu)化收斂曲線分別如圖13和圖14所示，其中，PDOF為自由度，k為迭代次數(shù)。從圖中可知，傳統(tǒng)多重網(wǎng)格等幾何模型[23]和所提的改進(jìn)多重網(wǎng)格等幾何模型具有相同的優(yōu)化收斂過(guò)程和優(yōu)化結(jié)果，驗(yàn)證了提出的多重網(wǎng)格B樣條單元?jiǎng)偠染仃嚱y(tǒng)一映射法則在二維優(yōu)化問(wèn)題中的正確性。

(a)二維算例模型

圖13 二維主支撐結(jié)構(gòu)的傳統(tǒng)[23]與改進(jìn)多重網(wǎng)格模型優(yōu)化結(jié)果對(duì)比

圖14 二維主支撐結(jié)構(gòu)的傳統(tǒng)[23]與改進(jìn)多重網(wǎng)格模型收斂曲線對(duì)比

此外，上述兩種多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型的各層剛度矩陣計(jì)算內(nèi)存消耗與預(yù)處理時(shí)間對(duì)比結(jié)果如圖15所示。由圖15可知：在各層剛度矩陣計(jì)算的內(nèi)存消耗上，改進(jìn)模型較傳統(tǒng)模型減少99.76%～99.96%；在各層的剛度矩陣計(jì)算的預(yù)處理時(shí)間上，改進(jìn)模型較傳統(tǒng)模型縮短95.89%～99.62%。因此，對(duì)于二維拓?fù)鋬?yōu)化問(wèn)題，所提出的多重網(wǎng)格剛度映射方法可大幅減少多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型各層等幾何單元?jiǎng)偠染仃嚨膬?nèi)存消耗與預(yù)處理時(shí)耗。

圖15 二維主支撐結(jié)構(gòu)各層等幾何剛度矩陣所需的內(nèi)存空間與預(yù)處理時(shí)間對(duì)比

當(dāng)多重網(wǎng)格的最大層數(shù)均為4且最大網(wǎng)格尺寸變化時(shí)，傳統(tǒng)模型與改進(jìn)模型進(jìn)行二維主支撐結(jié)構(gòu)優(yōu)化所需的單元?jiǎng)偠染仃囉?jì)算總內(nèi)存消耗與總預(yù)處理時(shí)間對(duì)比如圖16所示。在不同網(wǎng)格尺寸下，改進(jìn)模型剛度矩陣數(shù)據(jù)的總體內(nèi)存消耗較傳統(tǒng)模型可減少98.95%～99.48%，總體預(yù)處理時(shí)間可縮短99.93%～99.96%。相比于傳統(tǒng)模型，改進(jìn)模型在多種多重網(wǎng)格情況下進(jìn)行二維優(yōu)化問(wèn)題求解所需的等幾何單元?jiǎng)偠染仃嚨目傮w內(nèi)存消耗與總預(yù)處理時(shí)耗均可大幅減少。

圖16 二維主支撐結(jié)構(gòu)在不同最大網(wǎng)格尺寸下剛度矩陣計(jì)算的總體內(nèi)存消耗與總預(yù)處理時(shí)間變化對(duì)比

三維主支撐結(jié)構(gòu)算例模型如圖12b所示，結(jié)構(gòu)的右上部分受到負(fù)載的作用。其中，等幾何網(wǎng)格的初始尺寸為40×20×10，多重網(wǎng)格的最大層數(shù)為2。三維主支撐結(jié)構(gòu)等幾何拓?fù)鋬?yōu)化的優(yōu)化結(jié)果和優(yōu)化收斂曲線分別如圖17和圖18所示。從圖中可知，傳統(tǒng)模型和改進(jìn)模型具有相同的優(yōu)化收斂過(guò)程和優(yōu)化結(jié)果，驗(yàn)證了提出的多重網(wǎng)格B樣條單元?jiǎng)偠染仃嚱y(tǒng)一映射法則在三維優(yōu)化問(wèn)題中的正確性。

圖17 三維主支撐結(jié)構(gòu)的傳統(tǒng)[23]與改進(jìn)多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型優(yōu)化結(jié)果對(duì)比

圖18 三維主支撐結(jié)構(gòu)的傳統(tǒng)[23]與改進(jìn)多重網(wǎng)格模型收斂曲線對(duì)比圖

上述兩種多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型的各層剛度矩陣計(jì)算內(nèi)存消耗與預(yù)處理時(shí)間對(duì)比結(jié)果如圖19所示。從圖19中可知：在各層的剛度矩陣計(jì)算的內(nèi)存消耗上，改進(jìn)模型較傳統(tǒng)模型減少99.97%～99.99%；在各層的剛度矩陣計(jì)算的預(yù)處理時(shí)間上，改進(jìn)模型較傳統(tǒng)模型縮短94.96%～94.99%。因此，對(duì)于三維拓?fù)鋬?yōu)化問(wèn)題，本文提出的多重網(wǎng)格剛度映射方法可大幅減少多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型的各層等幾何單元?jiǎng)偠染仃嚨膬?nèi)存消耗與預(yù)處理時(shí)耗。

圖19 三維主支撐結(jié)構(gòu)各層等幾何剛度矩陣所需的內(nèi)存空間與預(yù)處理時(shí)間對(duì)比

當(dāng)多重網(wǎng)格模型的最大層數(shù)為2且最大網(wǎng)格尺寸變化時(shí)，傳統(tǒng)模型與改進(jìn)模型進(jìn)行三維主支撐結(jié)構(gòu)優(yōu)化所需的單元?jiǎng)偠染仃嚳們?nèi)存消耗與總預(yù)處理時(shí)間對(duì)比結(jié)果如圖20所示。在不同網(wǎng)格尺寸下，改進(jìn)模型單元?jiǎng)偠人璧目們?nèi)存消耗較傳統(tǒng)模型減少99.98%～99.99%，所需的總體預(yù)處理時(shí)間縮短94.97%～95.52%。因此，本文改進(jìn)模型在不同的多重網(wǎng)格情況下進(jìn)行三維優(yōu)化問(wèn)題求解所需的等幾何單元?jiǎng)偠染仃嚨目傮w內(nèi)存消耗與總預(yù)處理時(shí)間上都較傳統(tǒng)模型有明顯優(yōu)勢(shì)。

圖20 三維主支撐結(jié)構(gòu)在不同網(wǎng)格尺寸下剛度矩陣計(jì)算的總體內(nèi)存消耗與總預(yù)處理時(shí)間變化對(duì)比

綜合優(yōu)化模型設(shè)計(jì)分析與實(shí)驗(yàn)結(jié)果分析，本文提出的改進(jìn)多重網(wǎng)格模型具有正確性以及內(nèi)存和預(yù)處理時(shí)間方面的高效性，具體表現(xiàn)為：本文方案在優(yōu)化過(guò)程與優(yōu)化結(jié)果上與傳統(tǒng)方案保持一致，證明了本文方案是傳統(tǒng)方案正確有效的替代方案；二維與三維模型的拓?fù)鋬?yōu)化剛度計(jì)算中，改進(jìn)模型的各層內(nèi)存消耗與各層預(yù)處理時(shí)間都較傳統(tǒng)模型具有優(yōu)勢(shì)；二維與三維模型的拓?fù)鋬?yōu)化剛度計(jì)算中，改進(jìn)模型的總內(nèi)存消耗與總預(yù)處理時(shí)間都較傳統(tǒng)模型具有優(yōu)勢(shì)。

4 結(jié)語(yǔ)

本文提出一種基于Bézier單元?jiǎng)偠扔成涞母倪M(jìn)多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型，解決了傳統(tǒng)的多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型剛度矩陣計(jì)算過(guò)程存在的單元數(shù)據(jù)存儲(chǔ)內(nèi)存消耗過(guò)大和預(yù)處理時(shí)間過(guò)長(zhǎng)的問(wèn)題。通過(guò)Bézier提取算子將B樣條基函數(shù)等效表達(dá)為具有空間保持特點(diǎn)的Bernstein基函數(shù)，建立了標(biāo)準(zhǔn)Bézier單元?jiǎng)偠染仃嚺c各層Bézier提取算子進(jìn)行各層級(jí)B樣條單元?jiǎng)偠染仃嚨慕y(tǒng)一表達(dá)，實(shí)現(xiàn)了多重網(wǎng)格等幾何拓?fù)鋬?yōu)化內(nèi)存結(jié)構(gòu)的優(yōu)化與B樣條單元數(shù)據(jù)預(yù)處理過(guò)程的簡(jiǎn)化。在二維和三維拓?fù)鋬?yōu)化問(wèn)題中，本文提出的模型在數(shù)據(jù)存儲(chǔ)負(fù)擔(dān)與預(yù)處理效率上均得到大幅度的改善，進(jìn)而驗(yàn)證了本文提出的改進(jìn)多重網(wǎng)格等幾何拓?fù)鋬?yōu)化模型的有效性。