陳立,曾孝云,黃文,張建飛
中國工程物理研究院 機械制造工藝研究所,綿陽 621900
隨著中國航空航天事業(yè)的大力發(fā)展,針對航空航天飛行器結(jié)構(gòu)的輕量化設(shè)計變得尤為重要。近年來點陣結(jié)構(gòu)因其具有超輕質(zhì)、高比剛度、高比強度等優(yōu)良的機械性能,以及減振降噪等功能屬性應用到航空航天工程領(lǐng)域[1]。航空航天器中存在大量支撐功能部件的支撐結(jié)構(gòu),這些功能部件在飛行過程中受到基礎(chǔ)加速度激勵工況,這就要求功能部件所處的位置不能發(fā)生過大的振動位移響應。在航空航天工程結(jié)構(gòu)中,往往是采用蒙皮和均勻點陣填充并配合力學性能校核進行設(shè)計[2]。然而,具有可變單胞密度的非均勻點陣結(jié)構(gòu)顯然能獲得更好的性能,因此面向位移響應最小化的非均勻點陣結(jié)構(gòu)優(yōu)化設(shè)計是一個值得研究的問題。
點陣結(jié)構(gòu)可以看成是由一系列周期性單胞微結(jié)構(gòu)組成。均勻化方法可以對單胞微結(jié)構(gòu)進行優(yōu)化設(shè)計,Sigmund[3]提出了基于均勻化理論的拓撲優(yōu)化方法可以對各向同性材料單胞構(gòu)型進行逆向設(shè)計,通常采用均勻化方法設(shè)計單胞結(jié)構(gòu)是從宏觀的結(jié)構(gòu)拓撲設(shè)計和微觀結(jié)構(gòu)設(shè)計兩個尺度進行。該方法需要宏觀結(jié)構(gòu)尺度遠大于微觀尺度。為了克服尺度問題,Wang 等[4-5]提出了一種參數(shù)化點陣結(jié)構(gòu),實現(xiàn)了宏觀拓撲結(jié)構(gòu)和微觀點陣結(jié)構(gòu)協(xié)同優(yōu)化。Wu 等[6]基于子結(jié)構(gòu)法提出了近似降階懲罰模型,將點陣單胞視為一個超單元,實現(xiàn)了點陣結(jié)構(gòu)與宏觀結(jié)構(gòu)協(xié)同設(shè)計。避免了均勻化理論帶來的尺度效應。Wu 和Li[7]采用擴展有限元方法實現(xiàn)了共形點陣結(jié)構(gòu)的尺寸優(yōu)化。Zhou 等[8]基于參數(shù)化點陣和多尺度有限元方法提出了實體-點陣的混合變量優(yōu)化模型,避免了均勻化理論中的尺度分離。上述工作均針對結(jié)構(gòu)的靜力學性能進行研究。在動力學研究領(lǐng)域,Xiao 等[9]同種單胞結(jié)構(gòu)不同密度的力學性能進行計算,實現(xiàn)在點陣夾心結(jié)構(gòu)的最大化基礎(chǔ)頻率的優(yōu)化。Vicente 等[10]提出了一種基 于BESO 的協(xié)同拓撲優(yōu)化方法,可以同時優(yōu)化周期微結(jié)構(gòu)和宏觀結(jié)構(gòu)兩個尺度的頻率響應問題,Zhang 等[11]通過構(gòu)造可連接性微結(jié)構(gòu)的方法實現(xiàn)了多尺度點陣結(jié)構(gòu)頻率響應最小的優(yōu)化設(shè)計。
點陣結(jié)構(gòu)中的桿件也可以看成是由一系列梁組成,采用梁單元對點陣結(jié)構(gòu)進建模,避免了均勻化方法帶來的尺度分離效應。點陣結(jié)構(gòu)的優(yōu)化設(shè)計就可以轉(zhuǎn)化為對每個梁截面尺寸的優(yōu)化。Tang 等[12]提出了一種基于格框算法的點陣生成方法,實現(xiàn)了點陣結(jié)構(gòu)的生成和最小柔順度的優(yōu)化設(shè)計,Chen 等[13]提出了一種基于網(wǎng)格的點陣結(jié)構(gòu)生成方法,實現(xiàn)了點陣結(jié)構(gòu)最小柔順度的優(yōu)化設(shè)計,Terriault 和Brailovski[14]提出了基于網(wǎng)格的共形點陣生成方法,并對其力學性能進行了預測。Zhou 和Li[15-16]采 用Michell 桁架方法對類桁架結(jié)構(gòu)的自然頻率進行了優(yōu)化設(shè)計。
上述對點陣結(jié)構(gòu)靜力學優(yōu)化設(shè)計、固有頻率優(yōu)化方面都已分別有了相應的研究。但點陣結(jié)構(gòu)在動力學激勵下的優(yōu)化方面僅有很少一部分研究內(nèi)容[11]。實際工程中的點陣結(jié)構(gòu)往往是三維復雜結(jié)構(gòu),有些結(jié)構(gòu)即包含蒙皮也包含內(nèi)部的點陣填充,這些結(jié)構(gòu)會受到外部各種周期性載荷和簡諧激勵,而這些周期性載荷可以通過傅里葉變化轉(zhuǎn)化為不同頻率下的簡諧激勵的疊加。因此研究點陣結(jié)構(gòu)在簡諧基礎(chǔ)加速度激勵下的優(yōu)化問題將變得很有意義。
受Zhang 等[1]文章的啟發(fā),結(jié)合航空航天工程對三維點陣結(jié)構(gòu)在關(guān)鍵位置的位移響應控制的迫切需求開展研究。提出通過優(yōu)化點陣結(jié)構(gòu)桿件的截面尺寸來降低結(jié)構(gòu)動響應的優(yōu)化方法,建立了簡諧加速度激勵下點陣結(jié)構(gòu)優(yōu)化的數(shù)學模型、推導了設(shè)計目標的靈敏度。實現(xiàn)了非均勻點陣結(jié)構(gòu)在基礎(chǔ)加速度激勵下的優(yōu)化設(shè)計。
首先給出點陣結(jié)構(gòu)的建模方法、簡諧加速度激勵下的點陣結(jié)構(gòu)優(yōu)化設(shè)計的基礎(chǔ)理論,包括大質(zhì)量法、結(jié)構(gòu)動響應分析方法等,然后給出點陣結(jié)構(gòu)在簡諧加速度激勵下的優(yōu)化模型。
點陣結(jié)構(gòu)是由三維空間中重復排列的一些列桿件組成。單胞是定義點陣結(jié)構(gòu)基本單元,一些典型的點陣單胞結(jié)構(gòu)如圖1 所示,單胞結(jié)構(gòu)可以看成是由桿系結(jié)構(gòu)組成,將單胞中桿件采用梁單元進行建模,桿件的每個交點都看成是單元節(jié)點,最終點陣結(jié)構(gòu)優(yōu)化問題轉(zhuǎn)為梁截面尺寸優(yōu)化問題[13]。面向大規(guī)模點陣結(jié)構(gòu)尺寸優(yōu)化問題,Zhou[17]和Stolpe[18]等都進行過相關(guān)研究。

圖1 一些典型的點陣單胞Fig.1 Some typical lattice unit cells
為實現(xiàn)大規(guī)模點陣結(jié)構(gòu)的建模,這里采用Tang 等[12]所提出的方法,將不同點陣單胞生成單胞庫,根據(jù)工程結(jié)構(gòu)的性能要求選擇單胞結(jié)構(gòu)進行填充。根據(jù)單胞的幾何尺寸對幾何體進行離散,然后生成大規(guī)模點陣結(jié)構(gòu),具體的優(yōu)化設(shè)計流程如圖2 所示。具體步驟為

圖2 點陣結(jié)構(gòu)優(yōu)化設(shè)計流程Fig.2 Optimization design process of lattice structure
1)確定優(yōu)化對象的優(yōu)化計算邊界條件。
2)對優(yōu)化對象的設(shè)計域進行劃分,并從點陣單胞庫中選擇單胞生成點陣結(jié)構(gòu)。
3)在簡諧加速度激勵下優(yōu)化點陣結(jié)構(gòu)的動力學響應。
4)完成優(yōu)化模型。
結(jié)構(gòu)在承受基礎(chǔ)激勵時的動力學平衡方程為
采用分塊矩陣形式,式(1)可以改寫為
在大質(zhì)量法中,通過在基礎(chǔ)處附著質(zhì)量很大的集中質(zhì)量ML,得到對應于Mbb的大質(zhì)量矩陣,此時對結(jié)構(gòu)響應起主要作用的是大質(zhì)量受到基礎(chǔ)傳遞來的加速度激勵而產(chǎn)生的慣性力。為保證基礎(chǔ)處原有的加速度場不變,需在大質(zhì)量點激勵自由度方向上施加相應的外力載荷。假設(shè)基礎(chǔ)節(jié)點激勵自由度方向受到的加速度激勵向量為則式(2)可以變換為
將式(3)左側(cè)第2 項展開為
通過大質(zhì)量法就可以將基礎(chǔ)加速度激勵等效施加到結(jié)構(gòu)與基礎(chǔ)交界面處的節(jié)點上。即基礎(chǔ)加速度為
式中:Γg為加速度激勵幅值向量;ω為加速度激勵頻率。此時等效的力激勵向量為
通過大質(zhì)量法將簡諧加速度激勵等效轉(zhuǎn)化為簡諧力激勵。并且在用模態(tài)位移法分析簡諧響應時,當不考慮結(jié)構(gòu)剛體模態(tài)的情況下,可獲得響應點的相對位移響應。
簡諧響應分析,即求在簡諧載荷作用下解線性結(jié)構(gòu)的穩(wěn)態(tài)響應。首先給出具有n自由度的結(jié)構(gòu)在簡諧力激勵下的動力學平衡方程:
式中:f(t)=Fejωt(j2=-1)為簡諧力向量;F為簡諧載荷振幅向量;ω為簡諧載荷激勵圓頻率。
在采用模態(tài)疊加法求解簡諧響應時,通過求解自由振動下的平衡方程,求解結(jié)構(gòu)的特征頻率和特征值:
則第i階固有頻率i和對應的第i階振型向量Φi可通過求解以下特征方程求得:
設(shè)質(zhì)量歸一化振型矩陣為Φ=[Φ1,Φ2,…,Φn],在經(jīng)典阻尼條件下求得:
式中:I為單位矩陣;ξi為第i階模態(tài)阻尼比。
假設(shè)系統(tǒng)的瑞利阻尼矩陣是質(zhì)量矩陣和剛度矩陣的線性組合[20],即
式中:α和β為瑞利阻尼系數(shù),取α=10-3、β=10-6。根據(jù)式(11)可以求解得到結(jié)構(gòu)模態(tài)阻尼比為
通過引入廣義坐標系變化公式,位移響應可以寫為
式中:y(t)為廣義坐標下的位移向量。將式(14)代入式(8)中,并左乘ΦiT,可以得到n個解耦的單自由度振動方程為
式(15)為標準的二階系統(tǒng)形式,可以求得:
式(16)右端的第1 部分即為第i個單自由度系統(tǒng)的頻率響應函數(shù),引入表達式:
根據(jù)Zhu 等[21]文章中的結(jié)論,在簡諧加速度激勵下,模態(tài)位移法(MDM)和模態(tài)加速度(MAM)有著相同的計算精度,通過模態(tài)位移法(MDM)可求得:
考慮到計算效率,求解中取結(jié)構(gòu)前q階模態(tài)進行計算[22],則:
在單頻簡諧基礎(chǔ)加速度激勵下的點陣結(jié)構(gòu)優(yōu)化設(shè)計中,以目標自由度相對位移響應幅值最小為優(yōu)化目標。以體積為約束的尺寸優(yōu)化模型,設(shè)計變量為單元的半徑rk。點陣結(jié)構(gòu)優(yōu)化的數(shù)學模型可表述為
式中:rk為在點陣結(jié)構(gòu)中第k個梁單元的半徑設(shè)計變量;ne為設(shè)計單元的數(shù)量;Vk為第k個梁單元的體積;V0為所有點陣結(jié)構(gòu)單元總和的體積上限值;rmin和rmax為點陣結(jié)構(gòu)桿件的最小半徑和最大半徑。
首先求解用模態(tài)位移法得到的簡諧響應關(guān)于設(shè)計變量的靈敏度,因為本文的優(yōu)化目標為某個自由度c 的位移響應,現(xiàn)假設(shè)a是一個在自由度c 上為1 其余項為0 的n維向量。則有:
自由度c 的位移響應幅值為
式中:Xc復位移幅值;xc(t)=Xcejωt,|xc(t)|關(guān) 于設(shè)計變量的靈敏度求解公式為
采用模態(tài)位移法求解簡諧響應xc(t)關(guān)于設(shè)計變量rk的靈敏度。對式(18)兩邊關(guān)于設(shè)計變量直接進行求導即可求得:
頻率靈敏度可以通過式(25)進行求解[23]:
振型的靈敏度可以通過不完備振型疊加法[24]進行求解,即采用前q階振型的線性組合表示第i階振型的靈敏度:
通過鏈式求導法則即可求得式(21)中簡諧位移響應的靈敏度。
對如圖3 所示的三維懸臂梁點陣結(jié)構(gòu)進行優(yōu)化設(shè)計,模型的幾何尺寸為210 mm×100 mm×40 mm。模型灰色部分為質(zhì)量塊,模擬結(jié)構(gòu)所受到的負載,灰色部分的中心為響應點。從點陣單胞庫中隨機選擇BC、Octa點陣單胞對設(shè)計與進行填充,單胞的尺寸為10 mm×10 mm×10 mm,填充后的模型如圖4 所示。材料的彈性模量、泊松比和密度分別為1.95 GPa、0.39、1 200 kg/m3,設(shè)點陣結(jié)構(gòu)中桿件半徑上下限分別為rmin=0.25 mm 和rmax=1.5 mm。設(shè)點陣結(jié)構(gòu)的體積上限V0=1.2×105mm3,此時BC 點陣結(jié)構(gòu)的一階共振頻率為48.9 Hz、Octa點陣結(jié)構(gòu)的一階共振頻率為126 Hz。

圖3 三維懸臂梁Fig.3 3D cantilever beam

圖4 懸臂梁點陣結(jié)構(gòu)Fig.4 Cantilever beam with lattice structure
給定模型左端支撐處沿Y方向施加激勵頻率ωp=40 Hz、幅值為5 m/s2的簡諧加速度激勵。此時激勵頻率低于結(jié)構(gòu)一階共振頻率。優(yōu)化目標為響應點沿Y方向的相對位移響應幅值最小。
優(yōu)化過程和最終優(yōu)化構(gòu)型如圖5 所示。結(jié)果顯示BC 點陣結(jié)構(gòu)的目標自由度的相對位移響應幅值由20.5 mm 下降到0.97 mm,Octa 點陣結(jié)構(gòu)的目標自由度的相對位移響應幅值由1.14 mm下降到0.30 mm,優(yōu)化后的結(jié)構(gòu)點陣密度分布清晰,優(yōu)化效果明顯。采用Octa 點陣單胞的結(jié)構(gòu)動響應更小,因為Octa 點陣結(jié)構(gòu)單胞桿件數(shù)量更多,優(yōu)化空間更大。

圖5 點陣結(jié)構(gòu)低頻激勵下的相對位移響應幅值迭代過程Fig.5 Iteration history of relative displacement response amplitudes under low-frequency excitation of lattice structure
當對BC 點陣結(jié)構(gòu)和Octa點陣結(jié)構(gòu),分別施加簡諧加速度激勵頻率50 Hz和130 Hz,且均高于一階共振頻率時。經(jīng)過優(yōu)化,結(jié)構(gòu)目標自由度的相對位移響應幅值分別由149 mm 下降到5.9 mm 和16.2 mm 下降到0.91 mm。但此時通過圖6 的優(yōu)化構(gòu)型可以看出,點陣結(jié)構(gòu)中材料主要集中在響應點附近,而支撐處材料變少,這一結(jié)果使得靜剛度很差,無法應用于工程當中。點陣結(jié)構(gòu)高頻激勵下的一階共振頻率迭代曲線如圖7所示。

圖6 點陣結(jié)構(gòu)高頻激勵下的相對位移響應幅值迭代過程Fig.6 Iterative history of relative displacement response amplitudes under high-frequency excitation of lattice structure

圖7 點陣結(jié)構(gòu)高頻激勵下的一階共振頻率迭代曲線Fig.7 Iterative curves of first resonance frequency under high-frequency excitation of lattice structure
針對高于共振頻率下簡諧激勵存在的問題,Olhoff 和Du[25]曾提出一種增加靜柔順度約束的方法用于解決無阻尼連續(xù)結(jié)構(gòu)在簡諧激勵下拓撲優(yōu)化問題。
在實際工程中,結(jié)構(gòu)通常受到某個頻段范圍內(nèi)的簡諧加速度激勵。針對這種載荷下的結(jié)構(gòu)優(yōu)化問題,Liu[26]和Zhu[21]等采用的以頻段內(nèi)相對位移響應幅值的積分最小為目標,很好的解決了這個問題。優(yōu)化數(shù)學模型如式(28)所示,具體詳見文獻[21]:
在圖4 所示的2 種單胞結(jié)構(gòu)的懸臂梁上分別施加頻段為[0,500]Hz 的簡諧加速度激勵,頻段內(nèi)相對位移響應幅值的積分迭代過程和優(yōu)化結(jié)構(gòu)如圖8 所示,可以看出,優(yōu)化過程中頻段內(nèi)相對位移響應幅值的積分數(shù)值均逐漸降低,并趨于收斂。兩種優(yōu)化結(jié)構(gòu)相較圖6 更符合靜剛度的要求,傳力路徑上的點陣結(jié)構(gòu)更粗,并且與頻點優(yōu)化的結(jié)構(gòu)有著明顯區(qū)別。

圖8 點陣結(jié)構(gòu)頻段內(nèi)相對位移響應幅值的積分迭代過程Fig.8 Iterative history of relative displacement response amplitude integration within band of lattice structure
對初始結(jié)構(gòu)和優(yōu)化后的結(jié)構(gòu)進行簡諧響應分析,結(jié)果如圖9 所示。相比于初始結(jié)構(gòu),優(yōu)化后的結(jié)構(gòu)在低頻處的相對位移響應幅值降低,一階共振頻率增大,在高頻處的響應得到很好的抑制。

圖9 點陣結(jié)構(gòu)頻段內(nèi)初始和優(yōu)化結(jié)構(gòu)的相對位移響應幅值Fig.9 Relative displacement response amplitudes of initial and optimized structures within band of lattice structure
某動量輪支架的設(shè)計構(gòu)型如圖10 所示,該結(jié)構(gòu)是基于模態(tài)加速度法提出的動力學響應拓撲優(yōu)化方法得到的優(yōu)化構(gòu)型進行重構(gòu)設(shè)計的結(jié)果。結(jié)構(gòu)通過底面4 個點進行固定,結(jié)構(gòu)頂面安裝有動量輪,采用BC 單胞對結(jié)構(gòu)進行填充,其單胞尺寸為10 mm×10 mm×10 mm,每個單元的初始半徑為0.3 mm,設(shè)點陣結(jié)構(gòu)中桿件半徑上下限分別為rmin=0.15 mm 和rmax=0.6 mm。給定點陣結(jié)構(gòu)的體積為8.8×105mm3,結(jié)構(gòu)蒙皮厚度為1 mm。結(jié)構(gòu)材料屬性為:彈性模量72 GPa,泊松比0.33,密度2 850 kg/m3。

圖10 動量輪支架組件Fig.10 Momentum wheel bracket assembly
結(jié)構(gòu)底部分別受到沿X軸和Y軸方向的簡諧加速度激勵,其中X方向的加速度激勵幅值為5 m/s2,Y方向的加速度激勵幅值為10 m/s2。在支架支撐的動量輪質(zhì)心處施加動量輪的質(zhì)量為5 kg,優(yōu)化目標設(shè)置為響應點在激勵頻段[0,1 500]Hz 內(nèi)沿X方向和Y方向的相對位移響應幅值曲線積分之和最小。
從圖11 中可以看出頻段內(nèi)相對位移響應幅值的積分完成收斂。從圖12 可以看出優(yōu)化結(jié)構(gòu)點陣密度分布清晰,在結(jié)構(gòu)4 個主承力柱上的點陣結(jié)構(gòu)更粗,說明4 個主承力柱對結(jié)構(gòu)的響應貢獻較大。

圖11 動量輪支架頻段內(nèi)相對位移響應幅值的積分迭代曲線Fig.11 Iterative curves of relative displacement response amplitude integration within band of momentum wheel bracket

圖12 動量輪支架初始和優(yōu)化結(jié)構(gòu)Fig.12 Initial and optimized structures of momentum wheel bracket
對初始結(jié)構(gòu)和優(yōu)化結(jié)構(gòu)分別進行頻率響應分析,獲得其相對位移響應幅值曲線。圖13 為X方向相對位移響應幅值曲線,圖14 為Y方向相對位移響應幅值曲線,可以看出結(jié)構(gòu)各階共振頻率有所提高,給定頻段內(nèi)的相對位移響應的峰值均有所降低,相比于原構(gòu)型動響應降低25%。最終通過增材制造技術(shù)將其加工成型,如圖15 所示。

圖13 動量輪支架初始結(jié)構(gòu)和優(yōu)化結(jié)構(gòu)相對位移響應幅值(X 方向)Fig.13 Relative displacement response amplitudes of initial and optimized structures of momentum wheel bracket(X direction)

圖14 動量輪支架初始結(jié)構(gòu)和優(yōu)化結(jié)構(gòu)相對位移響應幅值(Y 方向)Fig.14 Relative displacement response amplitudes of initial and optimized structures of momentum wheel bracket(Y direction)

圖15 動量輪支架Fig.15 Momentum wheel bracket
為了驗證上述設(shè)計的正確性,選定基于頻段簡諧加速度激勵優(yōu)化的三維懸臂梁點陣結(jié)構(gòu)進行實驗。為了對實驗過程進行對比,將初始結(jié)構(gòu)和優(yōu)化結(jié)構(gòu)的懸臂梁模型均采用SLA 工藝進行增材制造,所得到的樣件如圖16 所示,實驗采用東菱科技有限公司的ES-3-150 型電動振動測試系統(tǒng)進行振動測試。在振動臺側(cè)面和零件的側(cè)面分別粘貼加速度傳感器,通過測量振動臺面和被測試件的加速度以及被測試件的相位角對零件的相對位移響應幅值進行計算。


圖16 振動試驗Fig.16 Vibration test
設(shè)置簡諧加速度激勵的掃描頻率范圍為20~1 000 Hz,掃描速率為2 Oct/min,激勵加速度幅值為5 m/s2,通過測得的加速度值計算的相對位移響應幅值曲線如圖17 所示,可以看出,優(yōu)化結(jié)構(gòu)較初始結(jié)構(gòu)的位移響應幅值更小,一階共振頻率也更高,說明了優(yōu)化方法的有效性。

圖17 懸臂梁點陣結(jié)構(gòu)相對位移響應幅值Fig.17 Relative displacement response amplitudes of cantilever beam with lattice structure
1)本文建立了簡諧加速度激勵下點陣結(jié)構(gòu)優(yōu)化問題的數(shù)學模型,并對目標自由度相對位移響應關(guān)于設(shè)計變量的靈敏度進行了推導。
2)當簡諧加速度激勵的頻率低于初始結(jié)構(gòu)的1 階共振頻率時,經(jīng)過優(yōu)化可以獲得符合工程要求的結(jié)構(gòu),優(yōu)化后結(jié)構(gòu)不僅一階共振頻率提高,而且相對位移響應峰值也有所下降;當激勵頻率在高于一階共振頻率或在某一頻率范圍內(nèi)時,通過引入頻段優(yōu)化方法,經(jīng)過優(yōu)化同樣可以獲得符合工程要求的結(jié)構(gòu)并提高結(jié)構(gòu)的靜力學性能。
3)采用該優(yōu)化方法對動量輪支架組件進行優(yōu)化,實現(xiàn)動響應降低25%。
4)通過對三維點陣懸臂梁結(jié)構(gòu)進行振動實驗對優(yōu)化結(jié)果進行了驗證,證明了優(yōu)化的有效性。