圓柱殼體阻尼材料布局拓?fù)鋬?yōu)化研究

李 超，李以農(nóng)，施 磊，鄭 玲

(重慶大學(xué) 機(jī)械傳動國家重點(diǎn)實(shí)驗(yàn)室，重慶 400044)

為了抑制板、殼等結(jié)構(gòu)的振動，可采用在結(jié)構(gòu)表面敷設(shè)自由阻尼或約束阻尼的方法［1-3］。通常阻尼材料被完全覆蓋于待控結(jié)構(gòu)表面，結(jié)構(gòu)設(shè)計(jì)的目的就是在盡量增大結(jié)構(gòu)振動能量耗散前提下，確定阻尼材料的類型、層數(shù)及各層厚度。由于阻尼材料布局的空間拓?fù)涫枪潭ǖ?，因此，要?shí)現(xiàn)阻尼結(jié)構(gòu)動力學(xué)特性的最優(yōu)化是很困難的。

Bendsoe 和 Kikuchi［4］基于均勻化理論，提出了均勻化拓?fù)鋬?yōu)化方法，成為連續(xù)體拓?fù)鋬?yōu)化的重要方法。為解決求解過程中數(shù)值不穩(wěn)定問題，Sigmund等［5］提出了實(shí)體罰函數(shù)法(SIMP)。連續(xù)體拓?fù)鋬?yōu)化方法還有變厚度法和變密度法。盡管這些方法數(shù)學(xué)上嚴(yán)謹(jǐn)，但實(shí)際應(yīng)用中卻較為困難。后來，Xie等［6］提出了“漸進(jìn)優(yōu)化算法”，它通過逐步刪除結(jié)構(gòu)中利用率低的材料而得到一個(gè)滿應(yīng)力設(shè)計(jì)。漸進(jìn)優(yōu)化算法的關(guān)鍵是確定一個(gè)合適的準(zhǔn)則來評價(jià)一個(gè)單元是“高效率單元”或是“低效率單元”，從而決定單元的刪除或保留。漸進(jìn)優(yōu)化算法概念盡管在數(shù)學(xué)上并不很嚴(yán)格，但算法比較簡單，計(jì)算效率較高，有較強(qiáng)的工程實(shí)用性。

楊德慶［7］針對阻尼結(jié)構(gòu)的耗能特點(diǎn)，提出了阻尼胞單元和阻尼拓?fù)涿舳鹊雀拍?，對自由阻尼板進(jìn)行了阻尼材料布局優(yōu)化。韋勇［8］研究了自由阻尼及約束阻尼板的拓?fù)鋬?yōu)化問題，通過計(jì)算能耗最大的阻尼材料粘貼位置，提出了阻尼結(jié)構(gòu)快速拓?fù)鋬?yōu)化方法。此外，郭中澤等［9-10］對阻尼材料的布局優(yōu)化進(jìn)行了研究，榮見華、杜海珍等［11-14］提出一種基于應(yīng)變能的結(jié)構(gòu)拓?fù)潆p方向漸進(jìn)優(yōu)化算法，劉寅東［15］提出了通過力傳遞的路徑來構(gòu)造拓?fù)涞耐負(fù)鋬?yōu)化思路求解連續(xù)結(jié)構(gòu)的拓?fù)鋬?yōu)化問題。李以農(nóng)等［16］利用漸進(jìn)優(yōu)化算法，以模態(tài)阻尼比為目標(biāo)函數(shù)對約束阻尼梁進(jìn)行拓?fù)鋬?yōu)化，取得了良好的效果。以上文獻(xiàn)針對的優(yōu)化模型都較為簡單，梁和平板較多，而本文對圓柱殼體這一常用結(jié)構(gòu)，做了自由阻尼和約束阻尼結(jié)構(gòu)的拓?fù)鋬?yōu)化分析，得到了兩種不同結(jié)構(gòu)下的最優(yōu)拓?fù)洳季郑⒎治隽似淠B(tài)對拓?fù)鋬?yōu)化的影響。

圓柱薄殼結(jié)構(gòu)在工程上有廣泛的應(yīng)用，如魚雷外殼、飛機(jī)機(jī)身、艦船船體等，研究圓柱薄殼結(jié)構(gòu)的阻尼減振及材料布局拓?fù)鋬?yōu)化具有重要的理論意義和工程應(yīng)用價(jià)值。本文采用漸進(jìn)優(yōu)化算法，研究圓柱殼體阻尼材料布局的拓?fù)鋬?yōu)化問題。建立圓柱殼體阻尼減振拓?fù)鋬?yōu)化模型，導(dǎo)出結(jié)構(gòu)模態(tài)損耗因子對阻尼胞單元的靈敏度表達(dá)式，根據(jù)拓?fù)鋬?yōu)化準(zhǔn)則，通過逐步刪除阻尼材料中靈敏度最低的單元，達(dá)到對圓柱殼體阻尼材料布局優(yōu)化的目的。

1 圓柱殼體阻尼結(jié)構(gòu)拓?fù)鋬?yōu)化模型

阻尼結(jié)構(gòu)振動能量消耗主要來源于阻尼材料剪切變形。因此，為實(shí)現(xiàn)結(jié)構(gòu)最大的振動能量耗散，通常以結(jié)構(gòu)模態(tài)損耗因子或結(jié)構(gòu)某關(guān)鍵點(diǎn)的位移為目標(biāo)函數(shù)，以阻尼材料刪除率作約束條件，建立阻尼結(jié)構(gòu)的拓?fù)鋬?yōu)化模型。本文以圓柱殼體的模態(tài)損耗因子為目標(biāo)函數(shù)，阻尼材料刪除率作約束條件，阻尼胞單元有或無為設(shè)計(jì)變量，建立拓?fù)鋬?yōu)化模型。設(shè)計(jì)變量為:

其中，n是阻尼胞單元總數(shù)。βi=1表明:殼體結(jié)構(gòu)表面第i個(gè)單元位置有阻尼材料，βi=0表明:殼體結(jié)構(gòu)表面第i個(gè)單元位置無阻尼材料。拓?fù)鋬?yōu)化模型描述為:

2 靈敏度分析

由模態(tài)應(yīng)變能方法可知，結(jié)構(gòu)的第k階模態(tài)損耗因子的計(jì)算公式為:

式中:Edk為結(jié)構(gòu)中阻尼材料的第k階模態(tài)應(yīng)變能;Esk為結(jié)構(gòu)的第k階模態(tài)應(yīng)變能;ηd為阻尼材料損耗因子。

當(dāng)?shù)趇個(gè)單元阻尼材料被刪除時(shí)，結(jié)構(gòu)的第k階模態(tài)損耗因子的變化量為:

在漸進(jìn)優(yōu)化(ESO)方法中，由于刪除單元數(shù)量在每次迭代中是有限的，因此，做如下近似:

其中:Edki為第k階模態(tài)中第i個(gè)阻尼材料單元的應(yīng)變能。

由(4)式～(6)式可得:

式(7)表示第i個(gè)阻尼材料單元?jiǎng)h除后，引起結(jié)構(gòu)第k階模態(tài)損耗因子的變化量，即為目標(biāo)函數(shù)對設(shè)計(jì)變量的靈敏度:

前M階模態(tài)損耗因子對設(shè)計(jì)變量的靈敏度則為:

這里，ωk為權(quán)重系數(shù)，可以根據(jù)模型對某階模態(tài)的敏感程度，有一定的變化，且滿足:

對阻尼材料進(jìn)行布局優(yōu)化，就是逐步刪除靈敏度絕對值最小的阻尼胞單元，使結(jié)構(gòu)模態(tài)損耗因子沿著最緩慢的路徑下降，使結(jié)構(gòu)在一定阻尼材料用量下，保持損耗因子最大化。

3 優(yōu)化計(jì)算流程

采用ANSYS參數(shù)化設(shè)計(jì)語言(APDL)，編制優(yōu)化設(shè)計(jì)程序。

首先，建立有限元分析模型，計(jì)算各階模態(tài)的應(yīng)變能，然后按照式(9)進(jìn)行靈敏度分析，設(shè)定要?jiǎng)h除的阻尼單元數(shù)量，去除靈敏度最小的那部分單元。刪除最小靈敏度單元后，留下的單元即為求得的每一輪的拓?fù)鋬?yōu)化結(jié)果，每一輪優(yōu)化結(jié)束后，都要檢查約束條件是否滿足，如果不滿足，迭代繼續(xù);如果滿足，停止迭代輸出結(jié)果。

此時(shí)，阻尼材料空間拓?fù)浞植技礊橥負(fù)鋬?yōu)化結(jié)果。因此，設(shè)計(jì)步驟可描述為:

(1)建立有限元分析模型;

(2)模態(tài)分析，計(jì)算結(jié)構(gòu)阻尼單元和非阻尼單元每階模態(tài)對應(yīng)的模態(tài)應(yīng)變能;

(3)由式(8)或式(9)計(jì)算目標(biāo)函數(shù)對各阻尼單元設(shè)計(jì)變量的靈敏度;

圖1 拓?fù)鋬?yōu)化程序流程圖Fig.1 Topological optimization program flow chart

(4)設(shè)定每次刪除的阻尼單元數(shù)量;

(5)刪除靈敏度最小的阻尼胞單元;

(6)判斷是否達(dá)到約束條件，如果不滿足，迭代繼續(xù)，如果滿足，迭代停止，輸出計(jì)算結(jié)果。

拓?fù)鋬?yōu)化流程圖見圖1所示。

4 數(shù)值算例

4.1 自由阻尼材料最優(yōu)布局

一圓柱形薄殼，兩端面受到徑向和軸向約束。長為100 mm，內(nèi)徑為250 mm，外徑260 mm，基層的楊氏模量為109 GPa，密度為4 440 kg/m3，泊松比為0.34;在圓柱形薄殼外表面覆蓋一層粘彈性材料，厚度為1 mm，楊氏模量57 MPa，密度為1 140 kg/m3，泊松比為0.3，損耗因子為1。采用solid45單元對圓柱基殼和粘彈性阻尼材料劃分網(wǎng)格。

在控制阻尼材料用量條件下，搜尋粘彈性阻尼材料的最優(yōu)拓?fù)浞植?，使?fù)合圓柱殼體的6、8、10階模態(tài)損耗因子平均值最大，即各階加權(quán)系數(shù)為1/3。

圖2～圖5分別是粘彈性阻尼材料刪除率為20%、40%、60%和80%條件下的阻尼材料最優(yōu)拓?fù)浞植?，靈敏度分別為 1.624 5 ×10-7、2.23 ×10-7、2.79 ×10-7和4.3×10-7。圖6是損耗因子平均值隨阻尼材料刪除率的變化。圖7為材料刪除率隨迭代次數(shù)增加的變化曲線。從圖2～圖5可知，由于兩端面存在約束，模態(tài)損耗因子對阻尼胞單元的靈敏度較小，先被刪除。此外，周向波節(jié)線處的靈敏度也較小，隨著刪除率增加，也將被刪除。從圖6可以看出:損耗因子平均值隨阻尼材料刪除率增加而減少，但下降速度并不是很快，表明:采用拓?fù)鋬?yōu)化技術(shù)，可充分利用阻尼材料的空間拓?fù)洳季?，使結(jié)構(gòu)保持較高的振動能量耗散特性。圖7可以看出每次迭代刪除的阻尼單元數(shù)大致不變，說明了迭代過程比較穩(wěn)定。

圖7 阻尼材料刪除率隨迭代次數(shù)變化Fig.7 Change of damping material delete rate with iteration number

4.2 約束阻尼材料最優(yōu)布局

一圓柱形薄殼，兩端面受到徑向和軸向約束。長100 mm，內(nèi)徑250 mm，外徑260 mm，基殼楊氏模量為109 GPa，密度為4 440 kg/m3，泊松比為0.34;在基殼上覆蓋一層粘彈性阻尼材料，厚度為1 mm，楊氏模量為57 MPa，密度為 1 140 kg/m3，泊松比為 0.3，損耗因子為1;在粘彈性阻尼材料上再覆蓋一層約束層，厚度為1 mm，楊氏模量為210 GPa，密度為7 800 kg/m3，泊松比為0.3?；鶜づc粘彈性夾心層，采用solid45單元?jiǎng)澐志W(wǎng)格，約束層采用shell63單元?jiǎng)澐志W(wǎng)格。

在控制約束阻尼材料用量條件下，搜尋約束阻尼材料的最優(yōu)拓?fù)浞植迹辜s束阻尼圓殼的6、8、10階模態(tài)損耗因子的平均值最大，即各階加權(quán)系數(shù)為1/3。

圖8～圖11為約束阻尼材料刪除率為20%、40%、60%和80%條件下的拓?fù)鋬?yōu)化布局，靈敏度分別為5.43 ×10-4、8.41 × 10-4、1.25 × 10-4和 1.36 × 10-3。從圖8～圖11可知，靈敏度較低的阻尼胞單元處于兩端面以及周波節(jié)線處，當(dāng)刪除率為60%時(shí)，這種效應(yīng)變得更加明顯。圖12為損耗因子平均值隨約束阻尼材料刪除率的變化。圖12表明:刪除率小于20%時(shí)，模態(tài)損耗因子比全覆蓋結(jié)構(gòu)還略有升高，刪除率達(dá)40%時(shí)，模態(tài)損耗因子與全覆蓋結(jié)構(gòu)相當(dāng)，繼續(xù)增大刪除率，則模態(tài)損耗因子開始下降，說明刪除的前40%阻尼材料對模態(tài)損耗因子的影響不大，但整個(gè)附加質(zhì)量卻減輕了。由此帶來的結(jié)果是阻尼材料用量減少，附加質(zhì)量減輕，而耗能作用不變，這就是采用拓?fù)鋬?yōu)化布局的效應(yīng)。圖13為材料刪除率隨迭代次數(shù)增加的變化曲線。從圖13可以看出:每次迭代刪除的阻尼單元數(shù)大致不變，說明了迭代過程比較穩(wěn)定。

通過對優(yōu)化模型的模態(tài)分析可知，圓柱殼體模態(tài)振動關(guān)于殼體中心軸對稱，是典型的空間對稱，而從自由阻尼結(jié)構(gòu)和約束阻尼結(jié)構(gòu)拓?fù)鋬?yōu)化布局來看，阻尼材料分布也是關(guān)于圓柱殼體中心軸對稱的，說明模態(tài)振型和模態(tài)損耗因子靈敏度有著密切的聯(lián)系。

圖8 約束阻尼材料最優(yōu)拓?fù)浞植?刪除率20%)Fig.8 Constraint damping topological placement(20%)

圖9 約束阻尼材料最優(yōu)拓?fù)浞植?刪除率40%)Fig.9 Constraint damping topological placement(40%)

圖10 約束阻尼材料最優(yōu)拓?fù)浞植?刪除率60%)Fig.10 Constraint damping topological placement(60%)

圖11 約束阻尼材料最優(yōu)拓?fù)浞植?刪除率80%)Fig.11 Constraint damping topological placement(80%)

圖12 約束阻尼模態(tài)損耗因子隨刪除率變化Fig.12 Change of constraint damping loss factor of modality with delete rate

圖13 約束阻尼材料刪除率隨迭代次數(shù)變化Fig.13 Change of constraint damping material delete rate with iteration number

5 結(jié)論

采用漸進(jìn)優(yōu)化算法(ESO)，以阻尼結(jié)構(gòu)模態(tài)損耗因子最大化為目標(biāo)，阻尼材料刪除率為約束條件，阻尼胞單元為設(shè)計(jì)變量，對自由阻尼以及約束阻尼材料布局進(jìn)行了拓?fù)鋬?yōu)化。研究了阻尼結(jié)構(gòu)模態(tài)損耗因子對阻尼胞單元位置的靈敏度，導(dǎo)出靈敏度計(jì)算表達(dá)式。研究結(jié)果表明:約束阻尼結(jié)構(gòu)的阻尼效果優(yōu)于自由阻尼結(jié)構(gòu)，且圓柱殼體模態(tài)振型和模態(tài)損耗因子靈敏度服從相同的空間分布規(guī)律，對自由阻尼或約束阻尼進(jìn)行布局優(yōu)化，能使圓柱殼體以較小的附加質(zhì)量，獲得較大的振動能量耗散，顯著降低圓柱殼體的振動響應(yīng)。采用ANSYS參數(shù)化設(shè)計(jì)語言(APDL)，編寫基于漸進(jìn)優(yōu)化方法(ESO)的拓?fù)鋬?yōu)化程序，進(jìn)行阻尼材料布局優(yōu)化，在工程上是可行的，具有很好的工程實(shí)用性。

［1］林洪斌，趙德有，鄭云龍.敷設(shè)粘彈性阻尼的加筋板振動和阻尼分析［J］.中國造船，2002，43(3):40-47.

［2］Rao D K.Frequency and loss factor of sandwich beams under various boundary conditions［J］.Journal of Mechanical Engineering Science，1978，20(5):271-282.

［3］李恩奇，雷勇軍，唐國金，等.基于傳遞函數(shù)方法的約束層阻尼梁動力學(xué)分析［J］.振動與沖擊，2007，26(2):75-78.

［4］ Bendsoe M P，Kikuchi N.Generating optimal topologies in structural design using a homogenization method［J］.Comp.Mech.in Appl.Mech ＆ Engrg.，1988，77:197-224.

［5］Sigmund O.Design of material structure using topology optimization［D］.Denmark:Department of Solid Mechanics，Technical University of Denmark，1994.

［6］ Xie Y M，Steven G P.Evolutionary structural optimization［M］.Berlin，Heidelberg，New York:Springer，1997.

［7］楊德慶，柳擁軍.自由阻尼層結(jié)構(gòu)阻尼材料配置優(yōu)化的拓?fù)涿舳确ǎ跩］.振動工程學(xué)報(bào)，2003，16(4):420-425.

［8］韋 勇，陳國平.一般阻尼結(jié)構(gòu)的模態(tài)阻尼比優(yōu)化設(shè)計(jì)［J］.振動工程學(xué)報(bào)，2006，19(4):433-437.

［9］郭中澤，陳裕澤，侯 強(qiáng)，等.阻尼材料布局優(yōu)化研究［J］.兵工學(xué)報(bào)，2007，28(5):638-640.

［10］郭中澤，陳裕澤.基于準(zhǔn)則法的阻尼結(jié)構(gòu)拓?fù)鋬?yōu)化［J］.宇航學(xué)報(bào)，2009，30(6):2378-2391.

［11］Rong J H，Xie Y M，Yang X Y，et al.Topology optimization of structures under dynamic response constriants［J］.Journal of Sound and Vibration，2000，234(2):177-189.

［12］榮見華，姜節(jié)勝，胡德文，等.基于應(yīng)力及其靈敏度的結(jié)構(gòu)拓?fù)錆u進(jìn)優(yōu)化方法［J］.力學(xué)學(xué)報(bào)，2003，35(5):584-590.

［13］榮見華，唐國金，楊振興.一種三維結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)方法［J］.固體力學(xué)學(xué)報(bào)，2005，26(3):289-296.

［14］杜海珍，榮見華，傅建林，等.基于應(yīng)變能的雙方向結(jié)構(gòu)漸進(jìn)優(yōu)化方法［J］.機(jī)械強(qiáng)度，2005，27(1):72-77.

［15］劉寅東，卞 鋼.基于ANSYS的結(jié)構(gòu)拓?fù)鋬?yōu)化及其二次開發(fā)［J］.船舶力學(xué)，2006，10(2):120-125.

［16］李以農(nóng)，謝熔爐，王 宜，等.約束阻尼結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)的進(jìn)化算法［J］.重慶大學(xué)學(xué)報(bào)，2010，33(8):1-6.