靈敏度分析在鑄件優(yōu)化設(shè)計中的應(yīng)用

段小勇 尹中保

摘? 要：在鑄件設(shè)計時一般都是先經(jīng)過拓?fù)鋬?yōu)化然后再進行結(jié)構(gòu)設(shè)計，但拓?fù)鋬?yōu)化無法得到確定結(jié)構(gòu)的最優(yōu)尺寸，僅能得到給定邊界條件下的基本形式，因此在拓?fù)鋬?yōu)化后如何進行詳細(xì)尺寸設(shè)計對設(shè)計者要求較高，有時會設(shè)計出很多方案進行強度驗證，且這些方案往往不是最優(yōu)方案。在拓?fù)鋬?yōu)化分析完畢后再加入靈敏度優(yōu)化技術(shù)，可以快速得到滿足設(shè)計要求且重量最輕結(jié)構(gòu)的詳細(xì)尺寸，大大提高了設(shè)計效率。

關(guān)鍵詞：靈敏度分析;鑄件;優(yōu)化設(shè)計

中圖分類號：U463.85? ? 文獻標(biāo)識碼：A? ? 文章編號：1005-2550（2020）01-0082-04

Application of sensitivity analysis in casting optimal design

DUAN Xiao-yong ， YIN Zhong-bao

（ Dongfeng Liu Zhou Motor Co.， Ltd， Liu Zhou Guang Xi 541007， China ）

Abstract： In casting design， topology optimization is usually followed structural design， but topology optimization cant determine the optimal size of the structure， can only obtain the basic form under the given boundary conditions， so how to design the detailed size after topology optimization is more demanding for the designer， sometimes many schemes are designed to verify the strength， and these schemes are often not the optimal scheme. With the addition of sensitivity optimization technology after topology optimization analysis， the detailed size of the structure satisfying the design requirements and the lightest weight can be obtained quickly， and the design efficiency is greatly improved.

前? ? 言

鑄件在工業(yè)生產(chǎn)中已有非常廣泛的應(yīng)用，在零件設(shè)計要求越來越輕量化、精細(xì)化的當(dāng)下，鑄件結(jié)構(gòu)也呈現(xiàn)出越來越復(fù)雜化的特點，尤其是重型卡車上的鑄件。由于重型卡車使用工況非常惡劣，且很多鑄件都是連接多個功能模塊的，因此重型卡車上的鑄件不僅結(jié)構(gòu)異常復(fù)雜且強度要求非常高。同時國內(nèi)對卡車超載限制越來越嚴(yán)格，客戶希望車輛自重越輕越好，這樣才能在整車總重不變的情況下多拉貨物。為了滿足以上要求，鑄件的設(shè)計就必須在滿足強度要求的情況下重量盡可能的輕。因此如何做好鑄件的輕量化設(shè)計顯得尤為重要。目前國內(nèi)鑄件設(shè)計的一般步驟如圖1所示：

通過拓?fù)鋬?yōu)化設(shè)計所生成的三維模型還存在以下幾個問題：

1. 通過拓?fù)鋬?yōu)化設(shè)計僅得到給定約束條件下的基本形式;

2. 無法確定最優(yōu)結(jié)構(gòu)的具體尺寸[1]。

因此通過拓?fù)鋬?yōu)化設(shè)計后得到的零部件三維模型還需進行強度、剛度等校核，尤其是在強度不滿足使用要求時，采用什么樣的加強方案既能保證滿足強度要求又能使零件重量盡可能輕，這對設(shè)計者的經(jīng)驗要求很高，有時甚至需要做出十多個設(shè)計方案進行強度驗證、對比，浪費了大量的人力、物力，且設(shè)計出的方案重量還不一定達到相對最輕。以上問題，隨著CAE優(yōu)化技術(shù)的完善，使得通過軟件尋找最優(yōu)尺寸方案成為可能，在滿足強度等設(shè)計要求的同時，還可以得到零件具體尺寸。

1? ? 靈敏度分析基本理論

靈敏度優(yōu)化分析就是設(shè)計響應(yīng)對優(yōu)化變量的偏導(dǎo)數(shù)，對于有限元方程：

[K]{U}={P}

兩邊對設(shè)計變量X求偏導(dǎo)數(shù)就可得到關(guān)于位移的靈敏度[2]：

則對位移向量U的偏導(dǎo)數(shù)為：

一般來說，設(shè)計響應(yīng)是位移向量U的函數(shù)：

g={Q}T{U}

所以設(shè)計響應(yīng)對設(shè)計變量的偏導(dǎo)數(shù)為：

以上方法適合約束很多而設(shè)計變量較少的優(yōu)化問題，稱作直接法[3]。靈敏度即求導(dǎo)信息，是一種評價因設(shè)計變量或者參數(shù)改變而引起結(jié)構(gòu)響應(yīng)特性變化率的方法，結(jié)構(gòu)靈敏度的研究是一個很特別的領(lǐng)域，是當(dāng)前結(jié)構(gòu)力學(xué)和結(jié)構(gòu)工程領(lǐng)域的主要研究方向之一。目前應(yīng)用較多的是尺寸優(yōu)化靈敏度分析，一般用在駕駛室、車架這種鈑金件的優(yōu)化設(shè)計上，而鑄件的靈敏度優(yōu)化分析方法目前國內(nèi)應(yīng)用的較少。

2? ? 零件優(yōu)化設(shè)計

本文以某重型卡車發(fā)動機后懸缸體支座為例。該零件連接發(fā)動機及發(fā)動機后懸軟墊，其許用應(yīng)力要求為不大于225MPa。經(jīng)過拓?fù)鋬?yōu)化設(shè)計得到的三維模型應(yīng)力云圖如圖2所示：

從應(yīng)力云圖可知該零件最大應(yīng)力為296MPa，遠(yuǎn)大于許用應(yīng)力，因此需要對該鑄件設(shè)計方案進行加強。我司相關(guān)系統(tǒng)設(shè)計人員對該結(jié)構(gòu)進行了十多個加強方案的強度驗證，浪費了大量時間，零件設(shè)計完成節(jié)點未能保證。在滿足強度要求的前提下設(shè)計出零件的最輕重量為4.31kg。對該方案進行靈敏度優(yōu)化分析如下：

2.1? ?目標(biāo)函數(shù)

優(yōu)化的目標(biāo)是在保證滿足設(shè)計條件情況下結(jié)構(gòu)達到最輕。

2.2? ?設(shè)計變量

發(fā)動機后懸缸體支座優(yōu)化分析由Hyper works軟件的優(yōu)化模塊來完成，軟件在生成計算變量時，需要把優(yōu)化的結(jié)構(gòu)分成若干的單元，設(shè)計單元的形狀由一組控制節(jié)點來描述，移動控制節(jié)點的位置可以改變單元的形狀，從而生成一組形狀設(shè)計變量。由于該發(fā)動機后懸缸體支座為接近三角形的懸臂結(jié)構(gòu)，主要承重發(fā)動機在X、Y、Z三個方向的沖擊，因此主要考慮增加該支座上下兩個面的尺寸，減小側(cè)面尺寸。選取發(fā)動機后懸缸體支座7個主要尺寸作為設(shè)計變量T1-T7，并根據(jù)實際安裝及加工要求，確定設(shè)計變量取值范圍如表1所示：

7個變量設(shè)計的設(shè)置情況如圖3所示：

靈敏度優(yōu)化設(shè)計變量選取應(yīng)注意以下事項：

1）應(yīng)盡量細(xì)化設(shè)計變量，靈敏度分析出得零件尺寸也才能更加細(xì)化、準(zhǔn)確，這樣得到的零件的結(jié)構(gòu)尺寸才最優(yōu)化;

2）應(yīng)充分考慮零件鑄造工藝要求和現(xiàn)場裝配空間要求;

3）由于螺栓孔周圍一般是采用剛性連接，應(yīng)力也不準(zhǔn)確，因此設(shè)計變量應(yīng)盡量排除螺栓孔，這樣才能使分析結(jié)果收斂。

2.3? ?約束條件

零件約束條件的設(shè)置是否合理至關(guān)重要，既要保證零件的強度滿足使用要求，又要盡量使零件的重量、成本達到最低，同時還要充分考慮到鑄件在鑄造過程中可能存在的缺陷問題。該零件材料初步選定為QT450，其材料抗拉強度為450MPa。根據(jù)以往設(shè)計經(jīng)驗以及與同類競品發(fā)動機后懸缸體支座的強度對比，并考慮到該零件分析工況為極限沖擊工況，安全系數(shù)大于2.0即可。因此本零件的許用應(yīng)力為225MPa，即優(yōu)化的約束條件設(shè)置為225MPa.

2.4? ?分析結(jié)果

經(jīng)過hyperworks 優(yōu)化模塊16次迭代優(yōu)化，得到各變量對重量及最大應(yīng)力單元靈敏度如表2：

從靈敏度分析結(jié)果可知，T1、T3變量增大，對于關(guān)注位置的應(yīng)力是較靈敏的，而T2、T4、T5、T6、T7對體積靈敏度較高，而對應(yīng)力靈敏度較低。

采用靈敏度優(yōu)化分析的優(yōu)勢在于可以很方便的得到各變量對目標(biāo)函數(shù)及約束條件的貢獻量，然后對貢獻量較高的設(shè)計變量進行相應(yīng)調(diào)整，快速使分析結(jié)果達到收斂，并得到滿足約束條件時重量最輕變量值組合，因此在設(shè)計變量較多、設(shè)計限制范圍較大時，可以先給出一個較小的變量范圍分析出靈敏度，然后根據(jù)設(shè)計變量對體積、應(yīng)力等的靈敏度對變量上、下限值進行相應(yīng)調(diào)整，并把對零件體積、應(yīng)力靈敏度均很低的設(shè)計變量進行剔除，從而加快靈敏度優(yōu)化分析的速度。本次靈敏度優(yōu)化分析得到的各設(shè)計變量最終優(yōu)化結(jié)果如表3所示：

靈敏度優(yōu)化分析后零件結(jié)構(gòu)尺寸變動如圖4所示：

為了驗證靈敏度優(yōu)化分析結(jié)果是否有效及零件強度是否滿足材料許用應(yīng)力要求，按照分析得到的各設(shè)計變量具體尺寸，在三維軟件中對發(fā)動機后懸缸體支座數(shù)模進行相應(yīng)修改，然后再導(dǎo)入有限元分析軟件里建立分析模型并求解，得到靈敏度優(yōu)化分析后零件應(yīng)力云圖如圖5所示：

通過對鑄件進行靈敏度優(yōu)化分析，最大應(yīng)力為222MPa，與約束條件225MPa接近，且單件重量為3.80kg，比經(jīng)過十多次方案修改得到的最優(yōu)方案還輕了0.51kg。對比發(fā)動機后懸缸體支座靈敏度優(yōu)化前后應(yīng)力云圖，發(fā)現(xiàn)優(yōu)化后應(yīng)力云圖分布更加均勻合理。

經(jīng)過靈敏度優(yōu)化后的發(fā)動機后懸支座經(jīng)過試制后裝在某重型牽引車上，如圖6所示。該車在帶掛并滿載總重55T的狀態(tài)下在襄陽試驗場進行了可靠性試驗驗證，最終優(yōu)化設(shè)計后的發(fā)動機后懸支座順利通過了試驗大綱規(guī)定的全部試驗，未出現(xiàn)失效現(xiàn)象。

3? ? 結(jié)論

1）通過拓?fù)鋬?yōu)化，僅能得到給定條件下結(jié)構(gòu)的基本形式，無法得到結(jié)構(gòu)的最優(yōu)尺寸，而靈敏度優(yōu)化分析能很好解決這一問題，因此鑄件類的優(yōu)化設(shè)計都應(yīng)先通過拓?fù)鋬?yōu)化得到最優(yōu)結(jié)構(gòu)形式，然后通過靈敏度分析確定該結(jié)構(gòu)形式的最優(yōu)尺寸，只有通過二者的結(jié)合，才能得到最優(yōu)結(jié)構(gòu);

2）本文分析中通過引入靈敏度優(yōu)化設(shè)計方法，經(jīng)過一次優(yōu)化分析就得到了滿足發(fā)動機后懸缸體支座強度要求的最優(yōu)結(jié)構(gòu)尺寸，且通過靈敏度優(yōu)化后的零件比經(jīng)過多次方案修改后的零件重量還輕了0.51kg，并通過了實車可靠性試驗驗證，大大提高了設(shè)計者的工作效率。

參考文獻：

[1]潘晉，陳昌亞等. 衛(wèi)星結(jié)構(gòu)的拓?fù)鋬?yōu)化和靈敏度分析[J]. 機械設(shè)計，2012 23（6）.

[2]M.Zhou，N. Pagaldipti， H.L.Thomas. 2000; An Integrated Approach For Topology， sizing And shape Optimization， AIAA-2000-4946.

[3]洪清泉，趙康等. OptiStruct﹠HyperSdudy理論基礎(chǔ)與工程應(yīng)用析[M]. 機械工業(yè)出版社，2012，12（1）：4-5.

[4]馮志華. 基于靈敏度分析的機械結(jié)構(gòu)動力學(xué)修改[J].蘇州絲綢工學(xué)院學(xué)報，2000，20（4）：46-50.

[5]王鈺棟，金磊，洪清泉等. Hyper mesh Hyper view應(yīng)用技巧與高級實例[M]. 機械工業(yè)出版社，2012.