考慮結(jié)構(gòu)穩(wěn)定性的變密度拓?fù)鋬?yōu)化方法

考慮結(jié)構(gòu)穩(wěn)定性的變密度拓?fù)鋬?yōu)化方法

在結(jié)構(gòu)優(yōu)化設(shè)計(jì)中，強(qiáng)度、剛度和穩(wěn)定性是三個(gè)主要問題?，F(xiàn)有的結(jié)構(gòu)拓?fù)鋬?yōu)化方法多以考慮結(jié)構(gòu)剛度為主,如體積約束下的最大剛度和剛度約束下的最小體積等變密度法1優(yōu)化問題數(shù)學(xué)模型基于經(jīng)典變密度法,建立以單元相對(duì)密度為設(shè)計(jì)變量,結(jié)構(gòu)最小柔度為目標(biāo)函數(shù),體積和失穩(wěn)載荷因子為約束條件的優(yōu)化問題數(shù)學(xué)模型,如式(1)所示:式中:ρ為單元相對(duì)密度矢量;C為結(jié)構(gòu)柔度;n為結(jié)構(gòu)中單元的數(shù)量;U為結(jié)構(gòu)節(jié)點(diǎn)位移矢量;K為結(jié)構(gòu)總剛度矩陣;u2單元相對(duì)密度使用變密度法進(jìn)行連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化時(shí),需要對(duì)中間密度值進(jìn)行懲罰,從而使單元相對(duì)密度趨向0或1,避免灰度單元的出現(xiàn)。本文使用目前最流行的SIMP(solidisotropicmicrostructureswithpenalization,固體各向同性懲罰微結(jié)構(gòu))法式中:E3計(jì)算的靈敏度對(duì)于式(1)中的結(jié)構(gòu)柔度C(ρ)、體積V和第1階失穩(wěn)載荷因子λ3.1材料插值模型及剛度矩陣根據(jù)式(1)中結(jié)構(gòu)柔度的表達(dá)式,可得柔度C的靈敏度表達(dá)式:式中U=K由材料插值模型及平面四邊形單元的剛度矩陣,可得k將U=K式中c3.2體積敏度體積靈敏度如式(6)所示:3.3失穩(wěn)載荷因子由結(jié)構(gòu)穩(wěn)定性線性特征值理論,可得經(jīng)典的穩(wěn)定性分析特征方程式中:λ式(7)兩邊同乘φ式(8)對(duì)設(shè)計(jì)變量ρ由于K,K根據(jù)式(10)、式(9)可簡(jiǎn)化為:將式(7)代入式(11),得:式(12)經(jīng)變換得到失穩(wěn)載荷因子λ式中:k4考慮穩(wěn)定性的變密度法利用上述優(yōu)化問題數(shù)學(xué)模型及計(jì)算得到的靈敏度,采用拉格朗日乘子法推導(dǎo)優(yōu)化準(zhǔn)則。首先構(gòu)建數(shù)學(xué)模型的拉格朗日函數(shù):式中:β當(dāng)目標(biāo)函數(shù)C取得極值時(shí),拉格朗日函數(shù)(14)需滿足以下庫恩-塔克(Kuhn-Tucker)條件:針對(duì)設(shè)計(jì)變量ρ考慮上式中μ將式(17)展開并整理,得到:其中,β整理得到:式中:E式中:ρ要計(jì)算得到f同理展開第1階失穩(wěn)載荷因子λ在優(yōu)化過程中,考慮到V5幾何剛度矩陣顯式推導(dǎo)在式(19)中,還需使用e式中:G是與單元形函數(shù)有關(guān)的微分矩陣,S是與單元應(yīng)力有關(guān)的矩陣,V為單元體積。對(duì)于平面四邊形單元,G和S的表達(dá)式如下式中:N將G和S代入式(24),可得平面四邊形單元的幾何剛度矩陣的顯式表示:式中:式中:r=m,q,u,v;s=m,q,u,v。對(duì)于平面四邊形單元,有:其中ξ,η為局部坐標(biāo),它們與整體坐標(biāo)x,y的轉(zhuǎn)換關(guān)系為:其中2a,2b為四邊形單元的邊長(zhǎng)。因此,可求得形函數(shù)對(duì)坐標(biāo)x,y的微分為:將式(32)代入式(29),得:6低密度區(qū)域的出現(xiàn)變密度拓?fù)鋬?yōu)化方法由于引入了中間密度單元,會(huì)在優(yōu)化的過程中出現(xiàn)低密度區(qū)域,導(dǎo)致虛假模態(tài)的出現(xiàn),影響優(yōu)化進(jìn)程。本文采用文獻(xiàn)[9]提出的虛假模態(tài)處理方法,以識(shí)別虛假模態(tài)。將低密度區(qū)域模態(tài)應(yīng)變能貢獻(xiàn)比值作為判斷指標(biāo),設(shè)定一個(gè)密度閥值ρ式中:E7數(shù)學(xué)模型的建立及初始化考慮結(jié)構(gòu)穩(wěn)定性的變密度拓?fù)鋬?yōu)化方法的具體流程如圖1所示。首先建立優(yōu)化問題的數(shù)學(xué)模型并初始化基本參數(shù);然后進(jìn)行靜力學(xué)分析,并提取單元應(yīng)力(用于構(gòu)建應(yīng)力矩陣S)及應(yīng)變能(用于求解式(20)中的應(yīng)變能e8第14次迭代前后的優(yōu)化過程及優(yōu)化結(jié)果圖2所示為受壓柱模型,設(shè)計(jì)域?qū)挾萕=1m,高度H=2m,結(jié)構(gòu)的材料參數(shù)及優(yōu)化設(shè)計(jì)參數(shù)如表1所示。在設(shè)計(jì)域頂端中部施加均布載荷P=10圖3為柔度C、體積V和失穩(wěn)載荷因子λ圖4為優(yōu)化過程中結(jié)構(gòu)的拓?fù)渥兓^程、失穩(wěn)模態(tài)云圖以及相應(yīng)的優(yōu)化過程數(shù)據(jù)。由圖可見結(jié)構(gòu)下部?jī)蓚?cè)的材料被保留以保證穩(wěn)定性和一定的剛度,同時(shí)在上部?jī)蓚?cè)及中間區(qū)域去除材料以減小體積。與圖3相對(duì)應(yīng)來看,在第14次迭代之前,結(jié)構(gòu)主要在頂部?jī)蓚?cè)及底部中央去除材料,優(yōu)化進(jìn)程較緩。在第14次迭代以后,結(jié)構(gòu)內(nèi)部開始出現(xiàn)孔洞,優(yōu)化進(jìn)程加快。最終,失穩(wěn)載荷因子約束條件破壞,優(yōu)化終止。為實(shí)現(xiàn)對(duì)比說明,針對(duì)不考慮穩(wěn)定性的拓?fù)鋬?yōu)化問題,利用僅考慮體積約束下以結(jié)構(gòu)柔度最小為目標(biāo)的傳統(tǒng)變密度法對(duì)上述模型進(jìn)行拓?fù)鋬?yōu)化。在上述考慮穩(wěn)定性約束的優(yōu)化問題中,優(yōu)化后模型的體積收斂為初始模型體積的0.58,因此使用V對(duì)優(yōu)化過程進(jìn)一步研究后發(fā)現(xiàn),當(dāng)改變網(wǎng)格單元的大小后,優(yōu)化過程會(huì)受到一定的影響。當(dāng)減小網(wǎng)格單元的大小時(shí),即使用42×84個(gè)四節(jié)點(diǎn)四邊形平面應(yīng)變單元對(duì)設(shè)計(jì)域進(jìn)行網(wǎng)格離散時(shí),保持其余條件不變,迭代16次后,失穩(wěn)載荷因子為74.87,體積與原始體積比值為0.59,拓?fù)鋬?yōu)化后得到的優(yōu)化結(jié)構(gòu)如圖6所示,與使用40×80個(gè)單元迭代16次得到的優(yōu)化結(jié)構(gòu)不同。通過分析發(fā)現(xiàn),在優(yōu)化過程中網(wǎng)格單元大小的改變會(huì)導(dǎo)致局部結(jié)構(gòu)失穩(wěn)現(xiàn)象的出現(xiàn),從而影響最終的優(yōu)化結(jié)果。9優(yōu)化問題的求解通過引入失穩(wěn)載荷因子約束,提出了一種考慮結(jié)構(gòu)穩(wěn)定性的變密度拓?fù)鋬?yōu)化方法。其中,通過推導(dǎo)單元的幾何剛度矩陣,計(jì)算幾何應(yīng)變能,從而得到失穩(wěn)載荷因子的靈敏度是求解優(yōu)化問題的關(guān)鍵。優(yōu)化算例表明,該考慮穩(wěn)定性約束的變密度拓?fù)鋬?yōu)化方法能顯著提高結(jié)構(gòu)的穩(wěn)定性,需要注意的是,雖然該方法提高了結(jié)構(gòu)的穩(wěn)定性,但在柔度方面會(huì)有一定妥協(xié)。即在相同的約束條件下,結(jié)構(gòu)的剛度優(yōu)化和穩(wěn)定性優(yōu)化將會(huì)在一定程度上相互約