1 前言
隨著經(jīng)濟(jì)建設(shè)的發(fā)展，工程設(shè)計(jì)人員需要有更系統(tǒng)、更科學(xué)的設(shè)計(jì)思想和方法，以達(dá)到提高產(chǎn)品開(kāi)發(fā)效率、節(jié)約原材料、降低成本及提高產(chǎn)品質(zhì)量的目的，結(jié)構(gòu)優(yōu)化設(shè)計(jì)則是實(shí)現(xiàn)這些目的較佳手段[1]。由于設(shè)計(jì)變量類(lèi)型的不同，結(jié)構(gòu)優(yōu)化設(shè)計(jì)可以分為由易到難的四個(gè)不同層次：尺寸優(yōu)化、形狀優(yōu)化、形貌優(yōu)化和拓?fù)洳季謨?yōu)化。由于拓?fù)鋬?yōu)化設(shè)計(jì)的難度較大，被公認(rèn)為是當(dāng)前結(jié)構(gòu)優(yōu)化領(lǐng)域內(nèi)最具有挑戰(zhàn)性的課題之一。但是在工程應(yīng)用中，拓?fù)鋬?yōu)化可以提供概念性設(shè)計(jì)方案，取得的經(jīng)濟(jì)效益比尺寸優(yōu)化、形狀優(yōu)化更大，因此，拓?fù)鋬?yōu)化技術(shù)對(duì)工程設(shè)計(jì)人員更具吸引力，已經(jīng)成為當(dāng)今結(jié)構(gòu)優(yōu)化設(shè)計(jì)研究的一個(gè)熱點(diǎn)。本文結(jié)合作者的研究工作，對(duì)結(jié)構(gòu)拓?fù)鋬?yōu)化的發(fā)展和研究現(xiàn)狀進(jìn)行了回顧，并介紹了拓?fù)鋬?yōu)化技術(shù)在發(fā)動(dòng)機(jī)零部件設(shè)計(jì)中的成功應(yīng)用案例。
2 拓?fù)鋬?yōu)化設(shè)計(jì)的研究概述
結(jié)構(gòu)拓?fù)鋬?yōu)化按其研究對(duì)象的不同，一般分為兩類(lèi)：桁架類(lèi)結(jié)構(gòu)拓?fù)鋬?yōu)化和連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化。而目前的研究成果主要是從桁架結(jié)構(gòu)拓?fù)鋬?yōu)化方面發(fā)展起來(lái)的。1904 年Michell 對(duì)桁架結(jié)構(gòu)拓?fù)鋬?yōu)化進(jìn)行研究，他用解析方法研究了一個(gè)載荷作用下應(yīng)力約束的結(jié)構(gòu)，得到重量最輕的最優(yōu)桁架所應(yīng)滿足的條件，后來(lái)稱(chēng)為Michell 準(zhǔn)則，并把滿足Michell 準(zhǔn)則的桁架稱(chēng)為Michell 桁架。這是結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)理論研究的里程碑。Michell 理論在近幾十年得到了重要發(fā)展。Cox證明了Michell 桁架就是最小柔度設(shè)計(jì)；Hegeminer等將Michell 準(zhǔn)則推廣到剛度、動(dòng)力參數(shù)優(yōu)化以及非線性彈性等情況；Hemp糾正了其中的一些錯(cuò)誤，求解了多種不同荷載形式下Michell 桁架的具體形式；Rozvany對(duì)Michell 桁架的唯一性以及桿件的正交性做了討論，對(duì)Michell 準(zhǔn)則做了進(jìn)一步修正，求解了多種不同邊界約束條件下Michell 桁架的具體形式；周克民等采用有限元方法計(jì)算出了Michell 桁架。這些研究對(duì)于整個(gè)結(jié)構(gòu)優(yōu)化學(xué)科的發(fā)展都有深刻的意義。
必須指出，對(duì)于Michell 桁架的研究，均是在上個(gè)世紀(jì)60 年代以后的事情，正是在那時(shí)由于有限元方法、數(shù)學(xué)規(guī)劃法和計(jì)算機(jī)的發(fā)展導(dǎo)致結(jié)構(gòu)優(yōu)化成為一門(mén)獨(dú)立學(xué)科，數(shù)值方法成為重要手段，截面與幾何層次的優(yōu)化發(fā)展為拓?fù)鋵哟蔚难芯刻峁┝嘶A(chǔ)。這個(gè)時(shí)期代表性的研究工作開(kāi)始主要是以1964 年Dorn，Gomory，Greenberg等提出的“基結(jié)構(gòu)法”為基礎(chǔ)，廣泛開(kāi)展連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化研究。1988 年，Bends.e 和Kikuchi從程耿東和Olhoff等在最小柔順性實(shí)心彈性薄板優(yōu)化設(shè)計(jì)的研究工作中得到啟迪，引入了具有空心的單胞微結(jié)構(gòu)，提出了基于均勻化理論的均勻化方法和連續(xù)體拓?fù)鋬?yōu)化的概念，開(kāi)創(chuàng)了結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)研究的新局面。隨后，通過(guò)Tenek 和Hagiwara、程耿東和王健、周克民和胡云昌 等人的研究，發(fā)展和推進(jìn)了變厚度法理論的應(yīng)用。
變密度法則是受均勻化方法的啟發(fā)而發(fā)展的，它不必構(gòu)造微結(jié)構(gòu)，而是人為地引入一種假想密度在0～1 之間可變的材料，以材料密度為拓?fù)湓O(shè)計(jì)變量，直接定義一個(gè)經(jīng)驗(yàn)公式來(lái)表達(dá)密度與彈性模量間假定的函數(shù)關(guān)系，這樣結(jié)構(gòu)的拓?fù)鋬?yōu)化問(wèn)題就被轉(zhuǎn)換為材料的最優(yōu)分布問(wèn)題。1999 年Bends.e 和Sigmund 證實(shí)了該方法物理意義的存在性。曾金玲和雷雨成等 利用變密度法將拓?fù)鋬?yōu)化技術(shù)成功應(yīng)用于汽車(chē)零部件的設(shè)計(jì)和改進(jìn)。吳中博和李書(shū)則在航天航空方面利用變密度法拓?fù)鋬?yōu)化技術(shù)解決工程問(wèn)題。變密度法避免了均勻化方法的微結(jié)構(gòu)構(gòu)造，但是代之以人造密度的假定，造成最終必須處理中間密度的困難，這就要求研究人員具有豐富的拓?fù)鋬?yōu)化設(shè)計(jì)經(jīng)驗(yàn)，避免優(yōu)化結(jié)果出現(xiàn)棋盤(pán)格、分析不收斂等問(wèn)題。本文案例則全部基于變密度法，所以下文闡述了變密度法的基本原理。
受結(jié)構(gòu)靜力拓?fù)鋬?yōu)化發(fā)展的推動(dòng)，近年來(lái)，在結(jié)構(gòu)動(dòng)力拓?fù)鋬?yōu)化研究方面也取得了一些令人鼓舞的進(jìn)展。Lim 等用均勻化方法對(duì)膜結(jié)構(gòu)進(jìn)行了多個(gè)固有頻率極大化的拓?fù)鋬?yōu)化，Kikuchi 等 用改進(jìn)的均勻化方法求解了特征值拓?fù)鋬?yōu)化問(wèn)題，Min 等利用均勻化和直接積分方法對(duì)沖擊載荷作用下的薄板結(jié)構(gòu)進(jìn)行了拓?fù)鋬?yōu)化設(shè)計(jì)。Xie 等用進(jìn)化結(jié)構(gòu)優(yōu)化方法對(duì)連續(xù)體結(jié)構(gòu)頻率約束問(wèn)題進(jìn)行了拓?fù)鋬?yōu)化，包括最大化或最小化指定階頻率，最大化指定兩頻率間間隔等問(wèn)題，隨后又對(duì)結(jié)構(gòu)隨機(jī)響應(yīng)約束下的連續(xù)體結(jié)構(gòu)進(jìn)行了拓?fù)鋬?yōu)化。
2 拓?fù)鋬?yōu)化案例
利用澳汰爾公司OptiStruct 商業(yè)優(yōu)化軟件，作者在多種發(fā)動(dòng)機(jī)產(chǎn)品開(kāi)發(fā)和技術(shù)改進(jìn)項(xiàng)目中，對(duì)發(fā)動(dòng)機(jī)零部件進(jìn)行拓?fù)鋬?yōu)化設(shè)計(jì)，改善結(jié)構(gòu)性能，實(shí)現(xiàn)了產(chǎn)品的自主創(chuàng)新設(shè)計(jì)。
2.1 發(fā)動(dòng)機(jī)主軸承座的拓?fù)鋬?yōu)化設(shè)計(jì)（位移作約束）
發(fā)動(dòng)機(jī)運(yùn)轉(zhuǎn)期間，主軸承座承受多種載荷，這些載荷包括：螺栓預(yù)緊載荷、軸瓦過(guò)盈載荷及曲軸動(dòng)載荷等。目前，主軸承座的主要評(píng)價(jià)指標(biāo)是結(jié)構(gòu)的強(qiáng)度、剛度是否滿足設(shè)計(jì)需求。在明確主軸承座承載情況和設(shè)計(jì)要求的前提下，作者對(duì)某大馬力發(fā)動(dòng)機(jī)原有主軸承座進(jìn)行了最大爆發(fā)壓力工況下的有限元分析，分析模型及主軸承座軸瓦徑向變形量見(jiàn)圖1（a）、圖1 （b）和圖1（c）。通過(guò)主軸承座的強(qiáng)度分析和動(dòng)態(tài)疲勞安全系數(shù)分析可以得知：主軸承座的動(dòng)態(tài)疲勞安全系數(shù)為1.843，遠(yuǎn)遠(yuǎn)大于安全系數(shù)閥值1，所以主軸承座的強(qiáng)度足以滿足設(shè)計(jì)需求。而從圖1（b）可以得知軸瓦在變形后水平方向徑向減小0.0739mm ，已經(jīng)接近曲軸、軸瓦徑向間隙最小值0.079mm，這容易導(dǎo)致曲軸與軸瓦間缺少油膜潤(rùn)滑，形成干摩擦，最終導(dǎo)致曲軸磨損加劇，發(fā)動(dòng)機(jī)動(dòng)載荷增加，甚至機(jī)毀人亡的悲劇；另外從圖1（c）可以得知軸瓦在變形后上下方向徑向增加0.0971mm ，小于軸瓦徑向變形許可值0.147mm 。所以，根據(jù)有限元分析結(jié)果可以判斷：主軸承座在水平方向的剛度不足夠，應(yīng)該改進(jìn)現(xiàn)有結(jié)構(gòu)，提高其剛度性能。
圖1 發(fā)動(dòng)機(jī)主軸承座有限元分析模型及軸瓦徑向變形圖
圖2 發(fā)動(dòng)機(jī)主軸承座拓?fù)鋬?yōu)化模型和優(yōu)化后的幾何模型圖
根據(jù)上述分析，作者針對(duì)主軸承座水平方向剛度小的問(wèn)題，建立了以下拓?fù)鋬?yōu)化模型：
設(shè)計(jì)變量為整個(gè)主軸承座的單元相對(duì)密度；約束為軸瓦水平方向半徑變形量小于0.030mm； 目標(biāo)為主軸承座質(zhì)量最小。通過(guò)18 次迭代，優(yōu)化后的主軸承座拓?fù)湫螤詈凸饣蟮膸缀涡螤钊鐖D2 所示。將所得優(yōu)化模型重新進(jìn)行主軸承座強(qiáng)度、剛度校核，通過(guò)計(jì)算可以得知，優(yōu)化后的主軸承座較原有主軸承座質(zhì)量輕了6.87%；動(dòng)態(tài)疲勞安全系數(shù)為1.833，強(qiáng)度指標(biāo)略微下降0.54%，但依然遠(yuǎn)遠(yuǎn)超過(guò)安全系數(shù)閥值1，所以?xún)?yōu)化后的主軸承座仍滿足設(shè)計(jì)強(qiáng)度需求；通過(guò)安裝優(yōu)化主軸承座后的軸瓦變形量分析，可以得知，軸瓦在變形后水平方向徑向減小0.0626mm ，較原有安裝情況下的軸瓦變形減小15.29%，同時(shí)遠(yuǎn)小于曲軸、軸瓦徑向間隙最小值0.079mm，所以?xún)?yōu)化后的主軸承座已滿足水平剛度要求；雖然安裝優(yōu)化主軸承座后的軸瓦在變形后上下方向徑向增加0.0974mm ，較原有安裝情況下的軸瓦變形增大0.27%，但依然滿足小于軸瓦徑向變形許可值0.147mm 的剛度要求。因此，優(yōu)化后的主軸承座取代了原有設(shè)計(jì)方案，并實(shí)現(xiàn)了結(jié)構(gòu)創(chuàng)新設(shè)計(jì)。
2.2 發(fā)動(dòng)機(jī)后懸置支架的拓?fù)鋬?yōu)化設(shè)計(jì)（應(yīng)力作約束）
重型車(chē)上的發(fā)動(dòng)機(jī)后懸置支架因?yàn)榻?jīng)常承受較大的動(dòng)載荷而斷裂，導(dǎo)致嚴(yán)重車(chē)禍，因此，該結(jié)構(gòu)應(yīng)該保證具有足夠的動(dòng)態(tài)安全儲(chǔ)備系數(shù)。為了降低產(chǎn)品開(kāi)發(fā)費(fèi)用，縮短生產(chǎn)周期，某重型車(chē)發(fā)動(dòng)機(jī)后懸置支架在開(kāi)發(fā)設(shè)計(jì)階段要求校核結(jié)構(gòu)強(qiáng)度。在進(jìn)行有限元分析時(shí)，考慮發(fā)動(dòng)機(jī)后懸置支架經(jīng)常承受較大范圍的動(dòng)載荷，因此計(jì)算時(shí)增加了相應(yīng)的安全系數(shù)，載荷變化范圍設(shè)在：-4Fz ～ + 6Fz 之間，F(xiàn)z為發(fā)動(dòng)機(jī)后懸置支架勻速行駛于平路上的載荷。通過(guò)有限元分析和動(dòng)態(tài)疲勞安全系數(shù)分析，從圖3 的疲勞安全系數(shù)分布云圖中可以得知，該支架的最小安全系數(shù)為0.87，小于安全系數(shù)閥值1，因此該支架設(shè)計(jì)方案不滿足強(qiáng)度設(shè)計(jì)要求。
根據(jù)設(shè)計(jì)要求，在不改變?cè)撝Ъ芑A(chǔ)拓?fù)湫螤畹幕A(chǔ)上，通過(guò)增加加強(qiáng)筋的方式提高發(fā)動(dòng)機(jī)后懸置支架強(qiáng)度。因此，作者以支架的單元相對(duì)密度為設(shè)計(jì)變量；以支架危險(xiǎn)處的應(yīng)力作為約束，用支架質(zhì)量最小為目標(biāo)，在優(yōu)化過(guò)程中還增加了制造約束條件，保證所得結(jié)果能夠制造。通過(guò)29 次迭代，得到了圖4 所示的拓?fù)鋬?yōu)化形狀(圖4（a）)和光滑處理后的幾何形狀(圖4（b）)。并將所得優(yōu)化支架進(jìn)行強(qiáng)度校核，優(yōu)化后的支架安全系數(shù)分布云圖如圖5 所示，支架最小安全系數(shù)轉(zhuǎn)移到新加筋的根部，值為1.03，大于安全系數(shù)閥值1，而原有支架危險(xiǎn)處的安全系數(shù)升高到1.30，所以，優(yōu)化后的發(fā)動(dòng)機(jī)后懸置支架已經(jīng)滿足強(qiáng)度設(shè)計(jì)要求。
圖3 發(fā)動(dòng)機(jī)原有后懸置支架疲勞安全系數(shù)云圖
圖4 發(fā)動(dòng)機(jī)后懸置支架拓?fù)鋬?yōu)化模型和優(yōu)化后的幾何模型圖
圖5 優(yōu)化后的發(fā)動(dòng)機(jī)后懸置支架疲勞安全系數(shù)云圖
2.3 發(fā)動(dòng)機(jī)排氣制動(dòng)閥支架的拓?fù)鋬?yōu)化設(shè)計(jì)（頻率座約束）
排氣制動(dòng)閥支架主要起支撐發(fā)動(dòng)機(jī)渦輪增壓器的作用。某發(fā)動(dòng)機(jī)排氣制動(dòng)閥支架在使用過(guò)程中出現(xiàn)斷裂，為了找出斷裂原因，對(duì)該支架進(jìn)行了強(qiáng)度校核和排氣系統(tǒng)模態(tài)分析，強(qiáng)度分析模型和系統(tǒng)模態(tài)分析模型分別見(jiàn)圖6(a)和圖6(b)。
圖6 支架有限元分析模型和系統(tǒng)模態(tài)分析模型
經(jīng)過(guò)支架的靜態(tài)強(qiáng)度校核和動(dòng)態(tài)疲勞安全系數(shù)分析，其危險(xiǎn)點(diǎn)疲勞安全系數(shù)云圖如7 所示，支架的最小安全系數(shù)為1.267，已經(jīng)大于人們默認(rèn)的結(jié)構(gòu)安全最小閥值1，所以支架在靜態(tài)和動(dòng)態(tài)載荷作用下已經(jīng)滿足結(jié)構(gòu)的強(qiáng)度需求，不會(huì)導(dǎo)致支架損壞。另一方面，通過(guò)模態(tài)分析發(fā)現(xiàn)系統(tǒng)沿發(fā)動(dòng)機(jī)橫向1 階振動(dòng)固有頻率為103.7 Hz，支架振型如圖8 所示，其中灰色結(jié)構(gòu)為沒(méi)有發(fā)生振動(dòng)的形態(tài)，深色結(jié)構(gòu)為系統(tǒng)沿發(fā)動(dòng)機(jī)橫向1 階振動(dòng)時(shí)支架的振型。而該發(fā)動(dòng)機(jī)在常運(yùn)轉(zhuǎn)轉(zhuǎn)速下要求支架系統(tǒng)固有頻率必須在110 Hz 以上，所以可以判斷支架損壞主要是由于共振造成的。根據(jù)上述結(jié)論，為避免支架損壞，只需要在滿足現(xiàn)有結(jié)構(gòu)強(qiáng)度的基礎(chǔ)上，提高支架系統(tǒng)的第一階固有頻率。
圖7 支架疲勞安全系數(shù)分布云圖
在進(jìn)行動(dòng)態(tài)優(yōu)化設(shè)計(jì)時(shí)還存在與靜態(tài)優(yōu)化設(shè)計(jì)完全不同的難題。首先，動(dòng)態(tài)的目標(biāo)函數(shù)不像靜態(tài)目標(biāo)函數(shù)那么穩(wěn)定。例如，在頻率拓?fù)鋬?yōu)化設(shè)計(jì)過(guò)程中，由于各階模態(tài)振型容易切換，從而導(dǎo)致以固有頻率為目標(biāo)的標(biāo)準(zhǔn)函數(shù)不再光滑。第二，動(dòng)態(tài)拓?fù)鋬?yōu)化時(shí)，目標(biāo)函數(shù)的靈敏度分析可能是非連續(xù)的。第三，如果動(dòng)態(tài)拓?fù)鋬?yōu)化設(shè)計(jì)成為一個(gè)非凸性求解問(wèn)題，適用于靜態(tài)拓?fù)鋬?yōu)化設(shè)計(jì)的迭代原理在此時(shí)就失去作用，導(dǎo)致迭代終止，優(yōu)化無(wú)法進(jìn)行。由于這些原因，動(dòng)態(tài)拓?fù)鋬?yōu)化設(shè)計(jì)在工程實(shí)際上的應(yīng)用受到了較大的限制。本文在進(jìn)行支架的頻率拓?fù)鋬?yōu)化設(shè)計(jì)時(shí)，綜合考慮上述三點(diǎn)，建立了以支架頻率為約束的優(yōu)化模型，并且采用模態(tài)追蹤控制，避免了模態(tài)振型之間的相互切換。經(jīng)過(guò)25 次迭代，拓?fù)鋬?yōu)化后的模型和光滑后的幾何模型分別如圖9(a) 和圖9(b)所示。
圖8 原支架系統(tǒng)一階固有頻率振型（103.7 Hz ）
圖9 支架頻率拓?fù)鋬?yōu)化模型和優(yōu)化后的幾何模型圖
優(yōu)化后的支架經(jīng)過(guò)與原支架在相同條件下進(jìn)行靜、動(dòng)態(tài)強(qiáng)度分析和模態(tài)分析，優(yōu)化后的支架動(dòng)態(tài)疲勞安全系數(shù)都大于8，整個(gè)支架受力均勻，沒(méi)有應(yīng)力集中部位，強(qiáng)度較原有支架大大提高；而且優(yōu)化后的支架系統(tǒng)沿發(fā)動(dòng)機(jī)橫向1 階振動(dòng)固有頻率達(dá)到111.5 Hz，該頻率已經(jīng)達(dá)到目標(biāo)要求；另外，通過(guò)質(zhì)量分析對(duì)比，該支架比原有支架減輕25.2% ；化后的支架還簡(jiǎn)化了原有支架的加工工藝，一方面降低了生產(chǎn)成本，另一方面減少了沖壓過(guò)程中形成彎角褶皺的風(fēng)險(xiǎn)。
3 結(jié)論
根據(jù)本文發(fā)動(dòng)機(jī)各零部件拓?fù)鋬?yōu)化的例子成功展示了拓?fù)鋬?yōu)化技術(shù)是一個(gè)非常有用的工具，讓設(shè)計(jì)師可以根據(jù)產(chǎn)品的性能要求，在指定的設(shè)計(jì)空間內(nèi)快速、準(zhǔn)確的實(shí)現(xiàn)產(chǎn)品設(shè)計(jì)，另一方面，可以通過(guò)拓?fù)鋬?yōu)化改善結(jié)構(gòu)性能、減輕產(chǎn)品質(zhì)量，并找到一種全新的設(shè)計(jì)方案。
（轉(zhuǎn)載）
