基于偽密度法的汽車懸架擺臂拓?fù)鋬?yōu)化

孫旭

摘 要：通過建立汽車懸架的三維運(yùn)動(dòng)模型，運(yùn)用COSMOSmotion進(jìn)行動(dòng)力學(xué)仿真，獲取零件結(jié)構(gòu)優(yōu)化所需的載荷條件，以上擺臂為例，在HyperWorks中建立拓?fù)鋬?yōu)化模型，確定設(shè)計(jì)區(qū)域，并針對(duì)擺臂的不同制造工藝，采取相應(yīng)的制造約束條件，采用多工況單一目標(biāo)的方案，基于偽密度法對(duì)擺臂的結(jié)構(gòu)進(jìn)行拓?fù)鋬?yōu)化，實(shí)現(xiàn)了擺臂的概念設(shè)計(jì)目標(biāo)。

關(guān)鍵詞：懸架擺臂；偽密度法；制造約束；拓?fù)鋬?yōu)化

中圖分類號(hào)： TP202+.7 文獻(xiàn)標(biāo)志碼： A 文章編號(hào)：1005-2550（2014）02-0001-05

擺臂也稱控制臂，是汽車獨(dú)立懸架系統(tǒng)的重要安全件和功能件。擺臂在比較復(fù)雜的受力狀態(tài)下工作，要承受牽引力、制動(dòng)力、側(cè)向力和力矩等。在擺臂的設(shè)計(jì)中，其強(qiáng)度、固有頻率以及縱向和橫向的剛度要滿足指定要求。 擺臂采取的制造工藝一般為鍛造或鑄造。對(duì)擺臂的結(jié)構(gòu)設(shè)計(jì)，文獻(xiàn)[1][2]都通過拓?fù)鋬?yōu)化的設(shè)計(jì)方法進(jìn)行了初步的研究，文獻(xiàn)[1]通過ADAMS軟件仿真車輛行駛工況，獲得載荷條件，進(jìn)行了拓?fù)鋬?yōu)化，但在優(yōu)化中沒有考慮制造約束情況。文獻(xiàn)[2]的載荷條件直接由試驗(yàn)結(jié)果給定，僅考慮制造拔模約束。本文所進(jìn)行的汽車麥弗遜懸架上擺臂的概念設(shè)計(jì)是在HyperWorks平臺(tái)上，基于三維連續(xù)體偽密度法進(jìn)行的拓?fù)鋬?yōu)化，并充分考慮多種制造約束。

拓?fù)鋬?yōu)化是一種數(shù)學(xué)方法，能在給定的空間結(jié)構(gòu)中生成優(yōu)化的形狀與材料分布。結(jié)構(gòu)拓?fù)鋬?yōu)化的基本思想是將尋求結(jié)構(gòu)的最優(yōu)拓?fù)鋯栴}轉(zhuǎn)化為在給定的設(shè)計(jì)區(qū)域內(nèi)尋求最優(yōu)材料分布的問題。通過將區(qū)域離散成有限元網(wǎng)格，HyperWorks的優(yōu)化模塊（OptiStruct）為每個(gè)單元計(jì)算材料特性，在給定的約束條件下，利用近似與優(yōu)化算法更改材料的分布，以優(yōu)化用戶給定的設(shè)計(jì)目標(biāo)，當(dāng)目標(biāo)函數(shù)在任意連續(xù)三次迭代中的改變量低于給定公差時(shí)，即得到收斂結(jié)果。

在拓?fù)鋬?yōu)化中，每個(gè)單元的密度值應(yīng)取0或1，單元分別定義為空或?qū)嶓w。但大量離散變量的優(yōu)化是無法計(jì)算的，因此，必須用連續(xù)變量的方式表示材料的分配問題。偽密度法是結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)中的一種有效方法，通過人為引入一種假想材料，材料的密度是可變的?；谶@種方法，材料的剛度被假想成與密度成線性關(guān)系，每個(gè)單元的材料密度直接作為設(shè)計(jì)變量，在0-1之間連續(xù)變化，優(yōu)化的結(jié)果是材料的最優(yōu)化分布，材料的分布反映了結(jié)構(gòu)的最優(yōu)拓?fù)?。在OptiStruct中偽密度法也是實(shí)施制造約束的唯一方法[3] [4]。

本文就兩種不同的工藝采取相應(yīng)的約束條件，通過拓?fù)鋬?yōu)化功能對(duì)上擺臂進(jìn)行概念設(shè)計(jì)，優(yōu)化得到的結(jié)構(gòu)不僅要滿足所有載荷工況的約束要求且質(zhì)量更輕，為產(chǎn)品的具體結(jié)構(gòu)設(shè)計(jì)提供理論支撐。

1 上擺臂的拓?fù)鋬?yōu)化設(shè)計(jì)前準(zhǔn)備工作

1.1 模型基本參數(shù)及載荷的來源與確定

圖1 懸架運(yùn)動(dòng)模型

根據(jù)新款車型整車總布置，確定硬點(diǎn)位置，在SolidWorks建立麥弗遜懸架總體運(yùn)動(dòng)模型，如圖1所示。選取汽車典型行駛工況時(shí)[1][2]，運(yùn)用COSMOSmotion對(duì)其進(jìn)行動(dòng)力學(xué)分析，得到擺臂外端與轉(zhuǎn)向節(jié)連接處的作用力，以便獲取優(yōu)化設(shè)計(jì)有限元分析中所需的載荷條件[5]。

1.2 有限元模型建立

基于以上運(yùn)動(dòng)模型，建立上擺臂設(shè)計(jì)邊界，在SolidWorks環(huán)境建立拓?fù)鋬?yōu)化三維實(shí)體簡(jiǎn)化模型，并轉(zhuǎn)化為iges格式文件導(dǎo)入HyperWorks，進(jìn)行網(wǎng)格劃分。根據(jù)計(jì)算機(jī)的運(yùn)算能力和求解精度，選擇網(wǎng)格尺寸（單元尺寸）為30 mm，單元類型選用四面體，劃分網(wǎng)格后，整個(gè)擺臂共有27 725個(gè)節(jié)點(diǎn)，11 122個(gè)單元，所建立的有限元模型如圖2所示。為了真實(shí)準(zhǔn)確的反映結(jié)構(gòu)本身的實(shí)際情況，在擺臂接觸孔周圍以及與車架連接處將所有的節(jié)點(diǎn)通過剛性單元（RBE2）連接在一起，這些RBE2是一些沒有質(zhì)量的單元，連接比較容易，且接近實(shí)際結(jié)構(gòu)受力特性。

圖2 優(yōu)化前控制臂有限元模型

1.3 材料選取

在Collectors面板下的create子面板中創(chuàng)建一個(gè)材料集合器，card image類型選擇為MAT1：MAT1是指所選材料為線形、不受溫度影響、各向同行材料。若上擺臂采取鍛造工藝，選取40Cr 材料，彈性模量為2e+05，泊松比為0.3。

1.4 邊界條件的建立與載荷施加

在模型上施加載荷之前，要先建立一個(gè)載荷集合器，以便把要施加的載荷經(jīng)處理后存入計(jì)算機(jī)，為以后載荷施加的調(diào)用提供支持。麥弗遜懸架控制臂分析時(shí)通常固定前襯套X、Y、Z三個(gè)方向平動(dòng)，后襯套點(diǎn)Y、Z方向平動(dòng)，外球銷點(diǎn)Z方向平動(dòng)。

根據(jù)以上由COSMOSmotion獲得的制動(dòng)、轉(zhuǎn)向、凹坑上跳工況時(shí)的縱向、側(cè)向力和垂向力，在上擺臂與轉(zhuǎn)向節(jié)連接處添加X、Y、Z方向的載荷，施加到相應(yīng)的單元（RBE2單元）和節(jié)點(diǎn)（3698）上。如表1及如圖3所示。

表1 載荷的設(shè)定

圖3 多工況下節(jié)點(diǎn)載荷施加情況

利用創(chuàng)建出的邊界條件及載荷分別組合起來創(chuàng)建出制動(dòng)、轉(zhuǎn)向和凹坑上跳三種工況，表2所示為模型優(yōu)化前處理信息。

表2 上擺臂模型優(yōu)化前處理信息

2 上擺臂的拓?fù)鋬?yōu)化

2.1 優(yōu)化設(shè)計(jì)數(shù)學(xué)模型建立

根據(jù)概念設(shè)計(jì)要求，明確拓?fù)鋬?yōu)化問題，建立數(shù)學(xué)模型如下：

Minimize：VOL（X）=VOL（x1，x2，...，xn）

Subject To VonMises≤760

Disp1（X）≤0.05

Disp2（X）≤0.02

Disp3（X）≤0.04

式中：X=x1，x2，..，xn 是設(shè)計(jì)變量，為材料單元密度；VOL（X）為目標(biāo)函數(shù)，體積的最小化；VonMises為應(yīng)力約束；Disp1（X）為施加載荷的節(jié)點(diǎn)在工況1下的位移約束，合位移小于0.05；Disp2（X）為施加載荷的節(jié)點(diǎn)在工況2下的位移約束，合位移小于0.02；Disp3（X）為施加載荷的節(jié)點(diǎn)在工況3下的位移約束，合位移小于0.04。endprint

2.2 OptiStruct中的拓?fù)鋬?yōu)化設(shè)計(jì)過程

拓?fù)鋬?yōu)化設(shè)計(jì)過程主要分成以下三個(gè)部分：

（1）定義設(shè)計(jì)需要的拓?fù)鋬?yōu)化區(qū)域與設(shè)計(jì)變量。

拓?fù)鋬?yōu)化就是要在給定的設(shè)計(jì)空間內(nèi)找到最優(yōu)的材料分布，因此在優(yōu)化前必須要確定設(shè)計(jì)空間。也即要?jiǎng)澐謨?yōu)化設(shè)計(jì)區(qū)域及非設(shè)計(jì)區(qū)域，因?yàn)橥負(fù)鋬?yōu)化區(qū)域是由屬性來識(shí)別的，所以通過多個(gè)屬性可區(qū)分設(shè)計(jì)及非設(shè)計(jì)區(qū)域。進(jìn)入Optimization 模塊Topology 菜單選擇設(shè)計(jì)區(qū)域?qū)傩裕愋瓦x擇PSOLID 創(chuàng)建優(yōu)化變量，定義優(yōu)化的設(shè)計(jì)空間[3][6]。本例中將上擺臂與車架、轉(zhuǎn)向節(jié)相連接的三處銷孔區(qū)域?yàn)榉窃O(shè)計(jì)區(qū)域（如圖2黃色網(wǎng)格區(qū)），其余部分為優(yōu)化設(shè)計(jì)區(qū)域（如圖2中間部分藍(lán)色網(wǎng)格區(qū)）。

（2）設(shè)置優(yōu)化參數(shù)。

①定義響應(yīng)

在這個(gè)優(yōu)化問題中，目標(biāo)是體積的最小化，而約束是受力的3698號(hào)節(jié)點(diǎn)的位移。將在Responses面板中創(chuàng)建兩個(gè)響應(yīng)：一個(gè)是定義目標(biāo)函數(shù)的體積響應(yīng)，另一個(gè)是位移約束。位移響應(yīng)是線性靜態(tài)分析的結(jié)果，體積響應(yīng)是全局響應(yīng)，它定義組件的整體結(jié)構(gòu)、個(gè)體特征和材料。

②定義目標(biāo)函數(shù)

本文的目的是結(jié)構(gòu)零部件的輕量化設(shè)計(jì)，因此在Optimization的子面板Objective中，將體積響應(yīng)定義為目標(biāo)函數(shù)，指定其為最小。

③定義設(shè)計(jì)約束

在Optimization的子面板Dconstraints中，對(duì)每個(gè)子工況，將對(duì)已定義的合位移響應(yīng)加一個(gè)上限約束；在Optimization的子面板下的Stress Constraint中定義應(yīng)力約束。

④給定迭代約束條件

Hyperworks軟件可以通過定義最大迭代次數(shù)控制迭代過程，也可以通過定義收斂公差，自動(dòng)進(jìn)行多次迭代，直到滿足優(yōu)化參數(shù)中所規(guī)定的收斂公差。本文中用的是定義收斂公差，讓其自動(dòng)進(jìn)行多次迭代的方法，通過指定最大迭代次數(shù)收斂公差來控制迭代過程。

（3）進(jìn)行拓?fù)鋬?yōu)化。

如果優(yōu)化計(jì)算收斂，就可以得到上擺臂的拓?fù)鋬?yōu)化結(jié)果。

2.3 拓?fù)鋬?yōu)化結(jié)果分析

Hypermesh的后處理功能為拓?fù)鋬?yōu)化的結(jié)果分析和結(jié)果轉(zhuǎn)化提供了多種方式，以下通過查看密度等值面圖進(jìn)行分析。

本文采用的偽密度法是在承認(rèn)實(shí)際材料密度不變的情況下，引入一種假想的密度可變材料，將連續(xù)結(jié)構(gòu)體離散為有限元模型后，首先指定結(jié)構(gòu)中的每個(gè)有限單元的密度相同，再以每個(gè)單元的相對(duì)密度為設(shè)計(jì)變量[2]，通過調(diào)整每個(gè)單元的相對(duì)密度值Xe來實(shí)現(xiàn)結(jié)構(gòu)的增刪，則有：

ρ=Xe*ρ0 （1）

式中：Xe為每個(gè)單元的相對(duì)密度，即單元密度閥值；ρ0為在設(shè)計(jì)域里的每個(gè)單元的固有密度；ρ為拓?fù)湓O(shè)計(jì)變量。

經(jīng)過18步的迭代，結(jié)果收斂。Xe作為單元密度閥值，它決定了結(jié)構(gòu)材料的除去量。經(jīng)過多次調(diào)整Xe的值，認(rèn)為Xe值為0.15時(shí)優(yōu)化結(jié)果比較符合設(shè)計(jì)。在HyperView中，將單元密度閥值設(shè)置為0.15，得到的可設(shè)計(jì)區(qū)域內(nèi)的最優(yōu)化布局的密度等值面圖，如圖4所示。

圖4 最優(yōu)化布局密度等值面圖

為將以上得到的優(yōu)化布局結(jié)構(gòu)圖生成CAD模型，需返回HyperMesh窗口，在后處理Post 面板將第18步結(jié)果讀入，并使用OSSmooth 命令將優(yōu)化后的單元格光順成IGES 曲面，單元密度閥值同樣取0.15。參考輸出的結(jié)果曲面，利用CAD 軟件對(duì)控制臂模型進(jìn)行參數(shù)化重建，重建的CAD 模型必須要經(jīng)過CAE 分析方可凍結(jié)，關(guān)于CAE 分析，本文不再討論。

基于以上拓?fù)鋬?yōu)化結(jié)果，可以在不知道上擺臂結(jié)構(gòu)拓?fù)湫螤畹那疤嵯?，根?jù)已知邊界條件和載荷條件確定較合理的結(jié)構(gòu)形式，它不涉及具體結(jié)構(gòu)尺寸設(shè)計(jì)，但可以為設(shè)計(jì)人員提供全新的設(shè)計(jì)和最優(yōu)的材料分布方案。

2.4 基于鑄造工藝的上擺臂拓?fù)鋬?yōu)化

對(duì)經(jīng)過以上拓?fù)鋬?yōu)化后的上擺臂結(jié)構(gòu)，如果采用鑄造方式制造，模具比較復(fù)雜，結(jié)構(gòu)也不對(duì)稱，難以實(shí)現(xiàn)。為了得到易于制造的結(jié)構(gòu)，本文對(duì)于采用鑄造工藝上擺臂，在拓?fù)鋬?yōu)化過程中添加對(duì)稱和脫模約束的條件來解決以上問題。

（1）考慮新增加的約束條件對(duì)模型的要求，使用的有限元模型與之前相比，單元?jiǎng)澐指?xì)，如圖5所示，但設(shè)計(jì)目標(biāo)、設(shè)計(jì)變量及已有的約束類型不變。

圖5 添加制造約束的上擺臂有限元模型

（2）為拓?fù)鋬?yōu)化設(shè)計(jì)變量定義對(duì)稱約束。

對(duì)稱約束也稱模式組（Pattern Grouping）約束，對(duì)設(shè)計(jì)空間施加對(duì)稱約束可以生產(chǎn)對(duì)稱設(shè)計(jì)。無論初始的網(wǎng)格、載荷和邊界條件如何，在拓?fù)鋬?yōu)化中加入對(duì)稱約束就可以得到實(shí)體模型的對(duì)稱結(jié)構(gòu)。在Optistruct中，對(duì)稱面有指定的Anchor和First節(jié)點(diǎn)定義，對(duì)稱面通過Anchor節(jié)點(diǎn)，并垂直于Anchor指向First節(jié)點(diǎn)的矢量。

（3）為拓?fù)鋬?yōu)化設(shè)計(jì)變量定義脫模約束。

對(duì)于鑄造件，必須考慮制造過程中的脫模，因此脫模方向上不能有材料的阻擋。脫模約束有單向脫模（Single）和沿給定方向分模（Split）兩種，施加脫模約束只需指定脫模方向既可。本文采取單方向脫模，脫模方向由作為Anchor和First的兩個(gè)節(jié)點(diǎn)確定，沿從Anchor Node指向First Node的矢量方向[3] [6]。

（4）提交Optistruct運(yùn)算。

（5）使用HyperView進(jìn)行結(jié)果后處理分析。

圖6 施加對(duì)稱和拔摸約束的拓?fù)鋬?yōu)化結(jié)果

經(jīng)過調(diào)整單元密度閥值為0.3，得到基于偽密度法的拓?fù)鋬?yōu)化的密度云圖。并使用OSSmooth 命令將優(yōu)化后的單元格光順成IGES 曲面，生成iges文件，得到優(yōu)化結(jié)果，如圖6所示。

3 結(jié)語

本文運(yùn)用COSMOSmotion進(jìn)行動(dòng)力學(xué)仿真，獲取擺臂結(jié)構(gòu)優(yōu)化所需的載荷條件。在HyperWorks中建立拓?fù)鋬?yōu)化模型，著重闡述針對(duì)擺臂的不同制造工藝，采取了相應(yīng)的制造約束條件，采用多工況單一目標(biāo)的方案，基于偽密度法對(duì)擺臂的結(jié)構(gòu)進(jìn)行拓?fù)鋬?yōu)化，實(shí)現(xiàn)了擺臂的概念設(shè)計(jì)目標(biāo)。

通過優(yōu)化使擺臂材料達(dá)到一個(gè)最優(yōu)的分布，在既定的條件下實(shí)現(xiàn)控制臂的輕量化的概念設(shè)計(jì)。通過拓?fù)鋬?yōu)化概念設(shè)計(jì)，開發(fā)人員可以全面了解產(chǎn)品的結(jié)構(gòu)和功能特征，可以有針對(duì)性地對(duì)總體結(jié)構(gòu)和具體結(jié)構(gòu)進(jìn)行設(shè)計(jì)。特別在產(chǎn)品設(shè)計(jì)初期，僅憑經(jīng)驗(yàn)和想象進(jìn)行零部件的設(shè)計(jì)是不夠的，在適當(dāng)?shù)募s束條件下，充分利用拓?fù)鋬?yōu)化進(jìn)行分析，并結(jié)合豐富的設(shè)計(jì)經(jīng)驗(yàn)，可設(shè)計(jì)出滿足最佳技術(shù)條件和工藝條件的產(chǎn)品。但需注意的是，如果想得到可用于生產(chǎn)實(shí)踐的更為準(zhǔn)確的模型，還需進(jìn)一步的細(xì)節(jié)優(yōu)化和可制造性設(shè)計(jì)。

參考文獻(xiàn)：

[1]呂寶剛.越野車獨(dú)立懸架關(guān)鍵零部件的輕量化設(shè)計(jì)[D].吉林大學(xué)，2007.

[2]上官文斌，蔣翠翠等.汽車懸架控制臂的拓?fù)鋬?yōu)化與性能計(jì)算[J].汽車工程，2008，30（8）：709～712.

[3]張勝蘭等.基于HyperWorks的結(jié)構(gòu)優(yōu)化設(shè)計(jì)技術(shù)[M].北京：機(jī)械工業(yè)出版社，2008.

[4]曾慶強(qiáng).結(jié)構(gòu)拓?fù)鋬?yōu)化中的若干問題[D].大連理工大學(xué)，2006.

[5]張晉西，郭學(xué)琴.SolidWorks及COSMOSMotion機(jī)械仿真設(shè)計(jì)[M].北京：清華大學(xué)出版社，2007.

[6]Diaz A.and Sigmund O. Checkerboard patterns in layout optimization.Structural Optimization，1995， 10：40～45.endprint

2.2 OptiStruct中的拓?fù)鋬?yōu)化設(shè)計(jì)過程

拓?fù)鋬?yōu)化設(shè)計(jì)過程主要分成以下三個(gè)部分：

（1）定義設(shè)計(jì)需要的拓?fù)鋬?yōu)化區(qū)域與設(shè)計(jì)變量。

拓?fù)鋬?yōu)化就是要在給定的設(shè)計(jì)空間內(nèi)找到最優(yōu)的材料分布，因此在優(yōu)化前必須要確定設(shè)計(jì)空間。也即要?jiǎng)澐謨?yōu)化設(shè)計(jì)區(qū)域及非設(shè)計(jì)區(qū)域，因?yàn)橥負(fù)鋬?yōu)化區(qū)域是由屬性來識(shí)別的，所以通過多個(gè)屬性可區(qū)分設(shè)計(jì)及非設(shè)計(jì)區(qū)域。進(jìn)入Optimization 模塊Topology 菜單選擇設(shè)計(jì)區(qū)域?qū)傩裕愋瓦x擇PSOLID 創(chuàng)建優(yōu)化變量，定義優(yōu)化的設(shè)計(jì)空間[3][6]。本例中將上擺臂與車架、轉(zhuǎn)向節(jié)相連接的三處銷孔區(qū)域?yàn)榉窃O(shè)計(jì)區(qū)域（如圖2黃色網(wǎng)格區(qū)），其余部分為優(yōu)化設(shè)計(jì)區(qū)域（如圖2中間部分藍(lán)色網(wǎng)格區(qū)）。

（2）設(shè)置優(yōu)化參數(shù)。

①定義響應(yīng)

在這個(gè)優(yōu)化問題中，目標(biāo)是體積的最小化，而約束是受力的3698號(hào)節(jié)點(diǎn)的位移。將在Responses面板中創(chuàng)建兩個(gè)響應(yīng)：一個(gè)是定義目標(biāo)函數(shù)的體積響應(yīng)，另一個(gè)是位移約束。位移響應(yīng)是線性靜態(tài)分析的結(jié)果，體積響應(yīng)是全局響應(yīng)，它定義組件的整體結(jié)構(gòu)、個(gè)體特征和材料。

②定義目標(biāo)函數(shù)

本文的目的是結(jié)構(gòu)零部件的輕量化設(shè)計(jì)，因此在Optimization的子面板Objective中，將體積響應(yīng)定義為目標(biāo)函數(shù)，指定其為最小。

③定義設(shè)計(jì)約束

在Optimization的子面板Dconstraints中，對(duì)每個(gè)子工況，將對(duì)已定義的合位移響應(yīng)加一個(gè)上限約束；在Optimization的子面板下的Stress Constraint中定義應(yīng)力約束。

④給定迭代約束條件

Hyperworks軟件可以通過定義最大迭代次數(shù)控制迭代過程，也可以通過定義收斂公差，自動(dòng)進(jìn)行多次迭代，直到滿足優(yōu)化參數(shù)中所規(guī)定的收斂公差。本文中用的是定義收斂公差，讓其自動(dòng)進(jìn)行多次迭代的方法，通過指定最大迭代次數(shù)收斂公差來控制迭代過程。

（3）進(jìn)行拓?fù)鋬?yōu)化。

如果優(yōu)化計(jì)算收斂，就可以得到上擺臂的拓?fù)鋬?yōu)化結(jié)果。

2.3 拓?fù)鋬?yōu)化結(jié)果分析

Hypermesh的后處理功能為拓?fù)鋬?yōu)化的結(jié)果分析和結(jié)果轉(zhuǎn)化提供了多種方式，以下通過查看密度等值面圖進(jìn)行分析。

本文采用的偽密度法是在承認(rèn)實(shí)際材料密度不變的情況下，引入一種假想的密度可變材料，將連續(xù)結(jié)構(gòu)體離散為有限元模型后，首先指定結(jié)構(gòu)中的每個(gè)有限單元的密度相同，再以每個(gè)單元的相對(duì)密度為設(shè)計(jì)變量[2]，通過調(diào)整每個(gè)單元的相對(duì)密度值Xe來實(shí)現(xiàn)結(jié)構(gòu)的增刪，則有：

ρ=Xe*ρ0 （1）

式中：Xe為每個(gè)單元的相對(duì)密度，即單元密度閥值；ρ0為在設(shè)計(jì)域里的每個(gè)單元的固有密度；ρ為拓?fù)湓O(shè)計(jì)變量。

經(jīng)過18步的迭代，結(jié)果收斂。Xe作為單元密度閥值，它決定了結(jié)構(gòu)材料的除去量。經(jīng)過多次調(diào)整Xe的值，認(rèn)為Xe值為0.15時(shí)優(yōu)化結(jié)果比較符合設(shè)計(jì)。在HyperView中，將單元密度閥值設(shè)置為0.15，得到的可設(shè)計(jì)區(qū)域內(nèi)的最優(yōu)化布局的密度等值面圖，如圖4所示。

圖4 最優(yōu)化布局密度等值面圖

為將以上得到的優(yōu)化布局結(jié)構(gòu)圖生成CAD模型，需返回HyperMesh窗口，在后處理Post 面板將第18步結(jié)果讀入，并使用OSSmooth 命令將優(yōu)化后的單元格光順成IGES 曲面，單元密度閥值同樣取0.15。參考輸出的結(jié)果曲面，利用CAD 軟件對(duì)控制臂模型進(jìn)行參數(shù)化重建，重建的CAD 模型必須要經(jīng)過CAE 分析方可凍結(jié)，關(guān)于CAE 分析，本文不再討論。

基于以上拓?fù)鋬?yōu)化結(jié)果，可以在不知道上擺臂結(jié)構(gòu)拓?fù)湫螤畹那疤嵯拢鶕?jù)已知邊界條件和載荷條件確定較合理的結(jié)構(gòu)形式，它不涉及具體結(jié)構(gòu)尺寸設(shè)計(jì)，但可以為設(shè)計(jì)人員提供全新的設(shè)計(jì)和最優(yōu)的材料分布方案。

2.4 基于鑄造工藝的上擺臂拓?fù)鋬?yōu)化

對(duì)經(jīng)過以上拓?fù)鋬?yōu)化后的上擺臂結(jié)構(gòu)，如果采用鑄造方式制造，模具比較復(fù)雜，結(jié)構(gòu)也不對(duì)稱，難以實(shí)現(xiàn)。為了得到易于制造的結(jié)構(gòu)，本文對(duì)于采用鑄造工藝上擺臂，在拓?fù)鋬?yōu)化過程中添加對(duì)稱和脫模約束的條件來解決以上問題。

（1）考慮新增加的約束條件對(duì)模型的要求，使用的有限元模型與之前相比，單元?jiǎng)澐指?xì)，如圖5所示，但設(shè)計(jì)目標(biāo)、設(shè)計(jì)變量及已有的約束類型不變。

圖5 添加制造約束的上擺臂有限元模型

（2）為拓?fù)鋬?yōu)化設(shè)計(jì)變量定義對(duì)稱約束。

對(duì)稱約束也稱模式組（Pattern Grouping）約束，對(duì)設(shè)計(jì)空間施加對(duì)稱約束可以生產(chǎn)對(duì)稱設(shè)計(jì)。無論初始的網(wǎng)格、載荷和邊界條件如何，在拓?fù)鋬?yōu)化中加入對(duì)稱約束就可以得到實(shí)體模型的對(duì)稱結(jié)構(gòu)。在Optistruct中，對(duì)稱面有指定的Anchor和First節(jié)點(diǎn)定義，對(duì)稱面通過Anchor節(jié)點(diǎn)，并垂直于Anchor指向First節(jié)點(diǎn)的矢量。

（3）為拓?fù)鋬?yōu)化設(shè)計(jì)變量定義脫模約束。

對(duì)于鑄造件，必須考慮制造過程中的脫模，因此脫模方向上不能有材料的阻擋。脫模約束有單向脫模（Single）和沿給定方向分模（Split）兩種，施加脫模約束只需指定脫模方向既可。本文采取單方向脫模，脫模方向由作為Anchor和First的兩個(gè)節(jié)點(diǎn)確定，沿從Anchor Node指向First Node的矢量方向[3] [6]。

（4）提交Optistruct運(yùn)算。

（5）使用HyperView進(jìn)行結(jié)果后處理分析。

圖6 施加對(duì)稱和拔摸約束的拓?fù)鋬?yōu)化結(jié)果

經(jīng)過調(diào)整單元密度閥值為0.3，得到基于偽密度法的拓?fù)鋬?yōu)化的密度云圖。并使用OSSmooth 命令將優(yōu)化后的單元格光順成IGES 曲面，生成iges文件，得到優(yōu)化結(jié)果，如圖6所示。

3 結(jié)語

本文運(yùn)用COSMOSmotion進(jìn)行動(dòng)力學(xué)仿真，獲取擺臂結(jié)構(gòu)優(yōu)化所需的載荷條件。在HyperWorks中建立拓?fù)鋬?yōu)化模型，著重闡述針對(duì)擺臂的不同制造工藝，采取了相應(yīng)的制造約束條件，采用多工況單一目標(biāo)的方案，基于偽密度法對(duì)擺臂的結(jié)構(gòu)進(jìn)行拓?fù)鋬?yōu)化，實(shí)現(xiàn)了擺臂的概念設(shè)計(jì)目標(biāo)。

通過優(yōu)化使擺臂材料達(dá)到一個(gè)最優(yōu)的分布，在既定的條件下實(shí)現(xiàn)控制臂的輕量化的概念設(shè)計(jì)。通過拓?fù)鋬?yōu)化概念設(shè)計(jì)，開發(fā)人員可以全面了解產(chǎn)品的結(jié)構(gòu)和功能特征，可以有針對(duì)性地對(duì)總體結(jié)構(gòu)和具體結(jié)構(gòu)進(jìn)行設(shè)計(jì)。特別在產(chǎn)品設(shè)計(jì)初期，僅憑經(jīng)驗(yàn)和想象進(jìn)行零部件的設(shè)計(jì)是不夠的，在適當(dāng)?shù)募s束條件下，充分利用拓?fù)鋬?yōu)化進(jìn)行分析，并結(jié)合豐富的設(shè)計(jì)經(jīng)驗(yàn)，可設(shè)計(jì)出滿足最佳技術(shù)條件和工藝條件的產(chǎn)品。但需注意的是，如果想得到可用于生產(chǎn)實(shí)踐的更為準(zhǔn)確的模型，還需進(jìn)一步的細(xì)節(jié)優(yōu)化和可制造性設(shè)計(jì)。

參考文獻(xiàn)：

[1]呂寶剛.越野車獨(dú)立懸架關(guān)鍵零部件的輕量化設(shè)計(jì)[D].吉林大學(xué)，2007.

[2]上官文斌，蔣翠翠等.汽車懸架控制臂的拓?fù)鋬?yōu)化與性能計(jì)算[J].汽車工程，2008，30（8）：709～712.

[3]張勝蘭等.基于HyperWorks的結(jié)構(gòu)優(yōu)化設(shè)計(jì)技術(shù)[M].北京：機(jī)械工業(yè)出版社，2008.

[4]曾慶強(qiáng).結(jié)構(gòu)拓?fù)鋬?yōu)化中的若干問題[D].大連理工大學(xué)，2006.

[5]張晉西，郭學(xué)琴.SolidWorks及COSMOSMotion機(jī)械仿真設(shè)計(jì)[M].北京：清華大學(xué)出版社，2007.

[6]Diaz A.and Sigmund O. Checkerboard patterns in layout optimization.Structural Optimization，1995， 10：40～45.endprint

2.2 OptiStruct中的拓?fù)鋬?yōu)化設(shè)計(jì)過程

拓?fù)鋬?yōu)化設(shè)計(jì)過程主要分成以下三個(gè)部分：

（1）定義設(shè)計(jì)需要的拓?fù)鋬?yōu)化區(qū)域與設(shè)計(jì)變量。

拓?fù)鋬?yōu)化就是要在給定的設(shè)計(jì)空間內(nèi)找到最優(yōu)的材料分布，因此在優(yōu)化前必須要確定設(shè)計(jì)空間。也即要?jiǎng)澐謨?yōu)化設(shè)計(jì)區(qū)域及非設(shè)計(jì)區(qū)域，因?yàn)橥負(fù)鋬?yōu)化區(qū)域是由屬性來識(shí)別的，所以通過多個(gè)屬性可區(qū)分設(shè)計(jì)及非設(shè)計(jì)區(qū)域。進(jìn)入Optimization 模塊Topology 菜單選擇設(shè)計(jì)區(qū)域?qū)傩裕愋瓦x擇PSOLID 創(chuàng)建優(yōu)化變量，定義優(yōu)化的設(shè)計(jì)空間[3][6]。本例中將上擺臂與車架、轉(zhuǎn)向節(jié)相連接的三處銷孔區(qū)域?yàn)榉窃O(shè)計(jì)區(qū)域（如圖2黃色網(wǎng)格區(qū)），其余部分為優(yōu)化設(shè)計(jì)區(qū)域（如圖2中間部分藍(lán)色網(wǎng)格區(qū)）。

（2）設(shè)置優(yōu)化參數(shù)。

①定義響應(yīng)

在這個(gè)優(yōu)化問題中，目標(biāo)是體積的最小化，而約束是受力的3698號(hào)節(jié)點(diǎn)的位移。將在Responses面板中創(chuàng)建兩個(gè)響應(yīng)：一個(gè)是定義目標(biāo)函數(shù)的體積響應(yīng)，另一個(gè)是位移約束。位移響應(yīng)是線性靜態(tài)分析的結(jié)果，體積響應(yīng)是全局響應(yīng)，它定義組件的整體結(jié)構(gòu)、個(gè)體特征和材料。

②定義目標(biāo)函數(shù)

本文的目的是結(jié)構(gòu)零部件的輕量化設(shè)計(jì)，因此在Optimization的子面板Objective中，將體積響應(yīng)定義為目標(biāo)函數(shù)，指定其為最小。

③定義設(shè)計(jì)約束

在Optimization的子面板Dconstraints中，對(duì)每個(gè)子工況，將對(duì)已定義的合位移響應(yīng)加一個(gè)上限約束；在Optimization的子面板下的Stress Constraint中定義應(yīng)力約束。

④給定迭代約束條件

Hyperworks軟件可以通過定義最大迭代次數(shù)控制迭代過程，也可以通過定義收斂公差，自動(dòng)進(jìn)行多次迭代，直到滿足優(yōu)化參數(shù)中所規(guī)定的收斂公差。本文中用的是定義收斂公差，讓其自動(dòng)進(jìn)行多次迭代的方法，通過指定最大迭代次數(shù)收斂公差來控制迭代過程。

（3）進(jìn)行拓?fù)鋬?yōu)化。

如果優(yōu)化計(jì)算收斂，就可以得到上擺臂的拓?fù)鋬?yōu)化結(jié)果。

2.3 拓?fù)鋬?yōu)化結(jié)果分析

Hypermesh的后處理功能為拓?fù)鋬?yōu)化的結(jié)果分析和結(jié)果轉(zhuǎn)化提供了多種方式，以下通過查看密度等值面圖進(jìn)行分析。

本文采用的偽密度法是在承認(rèn)實(shí)際材料密度不變的情況下，引入一種假想的密度可變材料，將連續(xù)結(jié)構(gòu)體離散為有限元模型后，首先指定結(jié)構(gòu)中的每個(gè)有限單元的密度相同，再以每個(gè)單元的相對(duì)密度為設(shè)計(jì)變量[2]，通過調(diào)整每個(gè)單元的相對(duì)密度值Xe來實(shí)現(xiàn)結(jié)構(gòu)的增刪，則有：

ρ=Xe*ρ0 （1）

式中：Xe為每個(gè)單元的相對(duì)密度，即單元密度閥值；ρ0為在設(shè)計(jì)域里的每個(gè)單元的固有密度；ρ為拓?fù)湓O(shè)計(jì)變量。

經(jīng)過18步的迭代，結(jié)果收斂。Xe作為單元密度閥值，它決定了結(jié)構(gòu)材料的除去量。經(jīng)過多次調(diào)整Xe的值，認(rèn)為Xe值為0.15時(shí)優(yōu)化結(jié)果比較符合設(shè)計(jì)。在HyperView中，將單元密度閥值設(shè)置為0.15，得到的可設(shè)計(jì)區(qū)域內(nèi)的最優(yōu)化布局的密度等值面圖，如圖4所示。

圖4 最優(yōu)化布局密度等值面圖

為將以上得到的優(yōu)化布局結(jié)構(gòu)圖生成CAD模型，需返回HyperMesh窗口，在后處理Post 面板將第18步結(jié)果讀入，并使用OSSmooth 命令將優(yōu)化后的單元格光順成IGES 曲面，單元密度閥值同樣取0.15。參考輸出的結(jié)果曲面，利用CAD 軟件對(duì)控制臂模型進(jìn)行參數(shù)化重建，重建的CAD 模型必須要經(jīng)過CAE 分析方可凍結(jié)，關(guān)于CAE 分析，本文不再討論。

基于以上拓?fù)鋬?yōu)化結(jié)果，可以在不知道上擺臂結(jié)構(gòu)拓?fù)湫螤畹那疤嵯?，根?jù)已知邊界條件和載荷條件確定較合理的結(jié)構(gòu)形式，它不涉及具體結(jié)構(gòu)尺寸設(shè)計(jì)，但可以為設(shè)計(jì)人員提供全新的設(shè)計(jì)和最優(yōu)的材料分布方案。

2.4 基于鑄造工藝的上擺臂拓?fù)鋬?yōu)化

對(duì)經(jīng)過以上拓?fù)鋬?yōu)化后的上擺臂結(jié)構(gòu)，如果采用鑄造方式制造，模具比較復(fù)雜，結(jié)構(gòu)也不對(duì)稱，難以實(shí)現(xiàn)。為了得到易于制造的結(jié)構(gòu)，本文對(duì)于采用鑄造工藝上擺臂，在拓?fù)鋬?yōu)化過程中添加對(duì)稱和脫模約束的條件來解決以上問題。

（1）考慮新增加的約束條件對(duì)模型的要求，使用的有限元模型與之前相比，單元?jiǎng)澐指?xì)，如圖5所示，但設(shè)計(jì)目標(biāo)、設(shè)計(jì)變量及已有的約束類型不變。

圖5 添加制造約束的上擺臂有限元模型

（2）為拓?fù)鋬?yōu)化設(shè)計(jì)變量定義對(duì)稱約束。

對(duì)稱約束也稱模式組（Pattern Grouping）約束，對(duì)設(shè)計(jì)空間施加對(duì)稱約束可以生產(chǎn)對(duì)稱設(shè)計(jì)。無論初始的網(wǎng)格、載荷和邊界條件如何，在拓?fù)鋬?yōu)化中加入對(duì)稱約束就可以得到實(shí)體模型的對(duì)稱結(jié)構(gòu)。在Optistruct中，對(duì)稱面有指定的Anchor和First節(jié)點(diǎn)定義，對(duì)稱面通過Anchor節(jié)點(diǎn)，并垂直于Anchor指向First節(jié)點(diǎn)的矢量。

（3）為拓?fù)鋬?yōu)化設(shè)計(jì)變量定義脫模約束。

對(duì)于鑄造件，必須考慮制造過程中的脫模，因此脫模方向上不能有材料的阻擋。脫模約束有單向脫模（Single）和沿給定方向分模（Split）兩種，施加脫模約束只需指定脫模方向既可。本文采取單方向脫模，脫模方向由作為Anchor和First的兩個(gè)節(jié)點(diǎn)確定，沿從Anchor Node指向First Node的矢量方向[3] [6]。

（4）提交Optistruct運(yùn)算。

（5）使用HyperView進(jìn)行結(jié)果后處理分析。

圖6 施加對(duì)稱和拔摸約束的拓?fù)鋬?yōu)化結(jié)果

經(jīng)過調(diào)整單元密度閥值為0.3，得到基于偽密度法的拓?fù)鋬?yōu)化的密度云圖。并使用OSSmooth 命令將優(yōu)化后的單元格光順成IGES 曲面，生成iges文件，得到優(yōu)化結(jié)果，如圖6所示。

3 結(jié)語

本文運(yùn)用COSMOSmotion進(jìn)行動(dòng)力學(xué)仿真，獲取擺臂結(jié)構(gòu)優(yōu)化所需的載荷條件。在HyperWorks中建立拓?fù)鋬?yōu)化模型，著重闡述針對(duì)擺臂的不同制造工藝，采取了相應(yīng)的制造約束條件，采用多工況單一目標(biāo)的方案，基于偽密度法對(duì)擺臂的結(jié)構(gòu)進(jìn)行拓?fù)鋬?yōu)化，實(shí)現(xiàn)了擺臂的概念設(shè)計(jì)目標(biāo)。

通過優(yōu)化使擺臂材料達(dá)到一個(gè)最優(yōu)的分布，在既定的條件下實(shí)現(xiàn)控制臂的輕量化的概念設(shè)計(jì)。通過拓?fù)鋬?yōu)化概念設(shè)計(jì)，開發(fā)人員可以全面了解產(chǎn)品的結(jié)構(gòu)和功能特征，可以有針對(duì)性地對(duì)總體結(jié)構(gòu)和具體結(jié)構(gòu)進(jìn)行設(shè)計(jì)。特別在產(chǎn)品設(shè)計(jì)初期，僅憑經(jīng)驗(yàn)和想象進(jìn)行零部件的設(shè)計(jì)是不夠的，在適當(dāng)?shù)募s束條件下，充分利用拓?fù)鋬?yōu)化進(jìn)行分析，并結(jié)合豐富的設(shè)計(jì)經(jīng)驗(yàn)，可設(shè)計(jì)出滿足最佳技術(shù)條件和工藝條件的產(chǎn)品。但需注意的是，如果想得到可用于生產(chǎn)實(shí)踐的更為準(zhǔn)確的模型，還需進(jìn)一步的細(xì)節(jié)優(yōu)化和可制造性設(shè)計(jì)。

參考文獻(xiàn)：

[1]呂寶剛.越野車獨(dú)立懸架關(guān)鍵零部件的輕量化設(shè)計(jì)[D].吉林大學(xué)，2007.

[2]上官文斌，蔣翠翠等.汽車懸架控制臂的拓?fù)鋬?yōu)化與性能計(jì)算[J].汽車工程，2008，30（8）：709～712.

[3]張勝蘭等.基于HyperWorks的結(jié)構(gòu)優(yōu)化設(shè)計(jì)技術(shù)[M].北京：機(jī)械工業(yè)出版社，2008.

[4]曾慶強(qiáng).結(jié)構(gòu)拓?fù)鋬?yōu)化中的若干問題[D].大連理工大學(xué)，2006.

[5]張晉西，郭學(xué)琴.SolidWorks及COSMOSMotion機(jī)械仿真設(shè)計(jì)[M].北京：清華大學(xué)出版社，2007.

[6]Diaz A.and Sigmund O. Checkerboard patterns in layout optimization.Structural Optimization，1995， 10：40～45.endprint