宋 凱 成艾國 陳少偉 王國春 鐘志華
湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國家重點(diǎn)實(shí)驗(yàn)室,長沙,410082
汽車車身各個(gè)關(guān)鍵截面的特性對整車的安全性、NVH特性和強(qiáng)度等性能至關(guān)重要。傳統(tǒng)的截面設(shè)計(jì)中,工程師先根據(jù)經(jīng)驗(yàn)在概念階段設(shè)計(jì)出汽車截面,只有在詳細(xì)設(shè)計(jì)階段才能通過CAE仿真評估乃至通過樣車試驗(yàn)才能知道截面的特性是否能滿足整車的性能要求。大量的工程實(shí)例證明,汽車設(shè)計(jì)進(jìn)入詳細(xì)設(shè)計(jì)階段以后再來修改車身關(guān)鍵截面是十分困難和昂貴的。
如何在汽車的概念設(shè)計(jì)階段就能找出設(shè)計(jì)中的缺陷、優(yōu)化汽車結(jié)構(gòu),對于成功的整車設(shè)計(jì)來說,是相當(dāng)重要的問題。近年來國內(nèi)外許多學(xué)者在汽車概念設(shè)計(jì)階段進(jìn)行了大量的分析與優(yōu)化,并取得了較好的成果[1]。Banichuk[2]提出了最大扭曲剛度的一種優(yōu)化方法,Egner等[3]提出了扭曲、彎曲、剪切相結(jié)合的一種優(yōu)化方法,Banichuk等[4]提出了具有正多邊形截面形狀的Bar單元,Kim[5]應(yīng)用拓?fù)鋬?yōu)化方法進(jìn)行截面設(shè)計(jì),Masataka等[6]應(yīng)用遺傳算法對截面形狀進(jìn)行了多目標(biāo)優(yōu)化。
在實(shí)際工程應(yīng)用中,工程師設(shè)計(jì)車身截面時(shí)不僅要考慮車身造型、內(nèi)部空間以及總體布置等方面的尺寸約束,還要考慮鈑金件沖壓成形等方面的形狀約束,更要在輕量化前提下設(shè)計(jì)出滿足車身性能要求或是具有指定截面特性的截面。顯然,采用傳統(tǒng)的經(jīng)驗(yàn)設(shè)計(jì)法較難快速設(shè)計(jì)出最優(yōu)的截面形狀。若能在考慮截面尺寸和形狀約束條件的前提下,建立一種能夠生成具有指定截面特性且滿足輕量化要求的截面形狀生成方法,將給結(jié)構(gòu)設(shè)計(jì)人員帶來較大的指導(dǎo)作用。
汽車車身的特性主要由車身關(guān)鍵截面特性、車身接頭特性以及車身的框架結(jié)構(gòu)特點(diǎn)決定。典型的車身關(guān)鍵截面一般是由若干層沖壓鈑金件通過焊接邊上的若干焊點(diǎn)連接而成的封閉截面。圖1所示為某車型門檻梁截面,本文以該截面為例來說明車身關(guān)鍵截面蟻群優(yōu)化算法的過程。
圖1 典型車身關(guān)鍵截面
汽車車身關(guān)鍵截面特性對整車的性能至關(guān)重要,它直接影響車身的剛度、振動(dòng)噪聲及舒適性、被動(dòng)安全性和強(qiáng)度等性能。截面特性可以用一些截面特性參數(shù)來表示,本文重點(diǎn)考慮如下幾個(gè)常用的截面特性參數(shù)。
(1)截面面積S:
式中,A為鈑金材料的面積區(qū)域;S為截面中鈑金的材料面積,較小的S有利于車身的輕量化。
(2)截面慣性矩Iy、Iz:
式中,y、z為局部坐標(biāo)系下截面微元的坐標(biāo)。
(3)截面主慣性軸和坐標(biāo)軸的夾角φ(逆時(shí)針方向?yàn)檎?/p>
由于汽車車身截面大多是由沖壓鈑金構(gòu)成的封閉截面,所以在進(jìn)行優(yōu)化時(shí)一定要滿足制造方面的形狀約束。首先,各層板不能出現(xiàn)交叉現(xiàn)象;其次,每層板不能出現(xiàn)沖壓負(fù)角,而且還要有一定的拔模角(本文暫不考慮有整形工序的沖壓件)。圖2所示為一些不滿足制造條件的典型截面形狀,這些截面形狀都不允許出現(xiàn)在截面優(yōu)化結(jié)果之中。
汽車車身關(guān)鍵截面的設(shè)計(jì)除了要考慮制造的可行性,還要考慮車身造型、內(nèi)部空間以及總體布置等方面的尺寸約束條件。
圖2 不滿足制造條件的車身截面形狀
圖3為某車型門檻梁在整車坐標(biāo)系中的示意圖,工程師在設(shè)計(jì)該車型門檻梁截面時(shí),要考慮諸多的外部設(shè)計(jì)約束。例如,確定的車門框開口尺寸決定了門檻梁截面門框膠條的位置;車輛內(nèi)部空間大小決定了門檻內(nèi)板和地板焊接處的位置;側(cè)門的密封膠條位置決定了門檻外板相應(yīng)密封面的位置;側(cè)門的包邊位置以及側(cè)門和門檻外板之間的間隙決定了門檻外板相應(yīng)部位的尺寸;整車的外部造型決定了門檻外板相應(yīng)部位的形狀;車身的最小離地間隙決定了門檻梁另外一個(gè)焊接邊的z向尺寸。
圖3 車身截面設(shè)計(jì)約束條件示意圖
一般來說,對于A柱、B柱、C柱和門檻梁之類的車身類梁型結(jié)構(gòu)而言,它們在整車坐標(biāo)系中的空間位置和截面本身特性是影響車身性能的重要因素。截面一般是指垂直于類梁型結(jié)構(gòu)形心軸的斷面,所以截面在車身坐標(biāo)系中的位置完全由類梁型結(jié)構(gòu)確定。在截面所在的平面內(nèi),不同局部坐標(biāo)系下的截面特性參數(shù)是不同的,所以在建立優(yōu)化模型之前,首先要確定截面的局部坐標(biāo)系。
坐標(biāo)系的選擇對帶有約束條件的蟻群優(yōu)化算法十分重要,它需要滿足兩個(gè)條件:首先,在該坐標(biāo)系下,截面的優(yōu)化結(jié)果不能出現(xiàn)沖壓負(fù)角或無拔模角等不滿足制造條件的情況;其次,在該坐標(biāo)系下,所得結(jié)果應(yīng)盡可能地包含所有滿足制造條件的截面形狀。
一般情況下,在確定拉伸件的沖壓方向時(shí)不把翻邊考慮在內(nèi)[7],翻邊和沖壓方向的角度可在后續(xù)工藝設(shè)計(jì)中實(shí)現(xiàn)。對于不包含需要整形工序零部件的典型截面而言,截面各零部件的沖壓方向大多近似平行,本文假設(shè)截面各零部件的沖壓方向是平行的,且截面的兩個(gè)焊接邊垂直于沖壓方向。
綜上所述,本文截面局部坐標(biāo)系的z軸平行于截面各零部件的沖壓方向,y軸垂直于z軸,取其中一個(gè)焊點(diǎn)為坐標(biāo)原點(diǎn)。圖4為某車型門檻梁截面局部坐標(biāo)示意圖,本文后續(xù)的論述均在該截面局部坐標(biāo)系下完成。
圖4 截面優(yōu)化坐標(biāo)系
如前所述,汽車車身關(guān)鍵截面的設(shè)計(jì)要受到車身造型、內(nèi)部空間以及總體布置等方面的尺寸約束,那么在截面設(shè)計(jì)與優(yōu)化之前會(huì)存在一個(gè)外部輸入的尺寸約束邊界,也就是截面形狀變化的可行域邊界。本文采用網(wǎng)格法的思想,將截面的可行域打網(wǎng)格,將每個(gè)網(wǎng)格節(jié)點(diǎn)作為截面形狀的控制點(diǎn)。
截面可行域邊界分為兩種:普通的可行域邊界表示截面形狀變化的最大允許邊界;固定的可行域邊界表示截面形狀必須采用該輪廓。例如汽車的外部造型曲線就是固定的可行域邊界,在設(shè)計(jì)截面時(shí),該造型曲線必須采用。一般情況下,固定的可行域邊界多為曲線,這里可以將其近似為若干段直線段,這樣只需給出若干個(gè)點(diǎn)的坐標(biāo)就能以若干直線段來近似代替復(fù)雜的曲線。
在截面局部坐標(biāo)系中,可行域邊界由上邊界、下邊界、左邊界和右邊界4個(gè)函數(shù)表示,其中,左邊界曲線為通過坐標(biāo)原點(diǎn)平行于z軸的直線,右邊界曲線為通過截面另一焊點(diǎn)平行于z軸的直線。截面的上下邊界函數(shù)曲線是由點(diǎn)和直線段組成的分段函數(shù)曲線,這些點(diǎn)由人工輸入,稱為邊界特征點(diǎn)。邊界特征點(diǎn)也分為兩種,其中,普通的邊界特征點(diǎn)表示截面邊界的變化情況,固定的邊界特征點(diǎn)表示的是截面形狀必須經(jīng)過該點(diǎn)。
截面上下邊界中的直線段是由邊界特征點(diǎn)確定的,在確定相鄰邊界特征點(diǎn)中的直線段時(shí)應(yīng)遵循以下兩個(gè)原則:①若該直線段為固定可行域邊界(例如外部造型邊界),則該處邊界曲線為連接相鄰固定邊界特征點(diǎn)的直線段;②若該直線段為普通可行域邊界,則該處邊界曲線為通過某點(diǎn)的平行線段,該點(diǎn)為相鄰邊界特征點(diǎn)中絕對值較大的一個(gè)。這樣可以保證產(chǎn)生的所有截面形狀均在截面可行域之內(nèi)。
綜上所述,設(shè)計(jì)人員只需用提供普通可行域邊界和固定可行域邊界的特征點(diǎn)的坐標(biāo),根據(jù)邊界曲線確定規(guī)則就可以快速確定邊界曲線。圖5為某車型門檻梁截面可行域邊界示意圖,對于該截面而言,u1和u2是由側(cè)門密封面確定的邊界特征點(diǎn),u3、u4和u5是門檻外部造型曲線簡化后的直線段端點(diǎn),u6是由外部空間確定的邊界特征點(diǎn),d1、d2和d3是由內(nèi)部空間確定的邊界特征點(diǎn)。該截面的可行域邊界函數(shù)如下:式中,yL(z)為截面可行域左邊界的函數(shù);yR(z)為截面可行域右邊界的函數(shù);yU(z)為截面可行域上邊界的函數(shù);yD(z)為截面可行域下邊界的函數(shù);uj為截面可行域上邊界的特征點(diǎn),j=1,2,…,J,yu1=0,yuJ=B,本文中J=6;B為截面可行域y方向的最大寬度;dk為截面可行域下邊界的特征點(diǎn),k=1,2,…,K,yd1=0,ydK=B,本文中K=3。
圖5 典型截面可行域邊界示意圖
確定了截面可行域邊界函數(shù)后,需要對截面可行域進(jìn)行離散化。從實(shí)際工程應(yīng)用的角度出發(fā),截面離散化時(shí)遵循兩個(gè)原則:首先,保證在每個(gè)截面邊界特征點(diǎn)處產(chǎn)生網(wǎng)格控制節(jié)點(diǎn);其次,根據(jù)截面的大小調(diào)整網(wǎng)格的密度,相鄰網(wǎng)格控制節(jié)點(diǎn)之間的最大距離不能過大,以保證截面形狀必要的精度。截面離散的具體步驟如下:
(1)在截面上邊界特征點(diǎn)uj處增加J′=J-1個(gè)縱斷面,則通過特征點(diǎn)uj的縱斷面的y坐標(biāo)為yuj。
比較截面下邊界特征點(diǎn)dk和上邊界特征點(diǎn)uj的y坐標(biāo)關(guān)系,假設(shè)滿足
的下邊界特征點(diǎn)有K′個(gè),在相應(yīng)的dk處增加縱斷面,則通過特征點(diǎn)dk的縱斷面的y坐標(biāo)為ydk。
(2)將截面可行域沿y方向平均分成L-1等份,則平分面al(l=1,2,…,L)處的y坐標(biāo)可表示為
比較yal和截面特征點(diǎn)uj、dk的坐標(biāo)關(guān)系,假設(shè)滿足
的平分面有L′個(gè),在相應(yīng)的平分面al處增加縱斷面。
(3)按照y坐標(biāo)由小到大的順序?qū)λ械目v斷面編號(hào)排序,則第m(m=1,2,…,M)個(gè)縱斷面的y坐標(biāo)ym是已知的,且縱斷面總數(shù)M可用下式求出:
(4)將每個(gè)縱斷面在可行域內(nèi)沿著z方向平均分成N-1份,則每個(gè)縱斷面上有N個(gè)節(jié)點(diǎn),截面可行域最后被離散為M×N個(gè)控制節(jié)點(diǎn)。第m個(gè)縱斷面上的第n(n=1,2,…,N)個(gè)節(jié)點(diǎn)坐標(biāo)為
圖6為典型截面離散化后的示意圖,在縱斷面和固定的可行域邊界交匯處會(huì)產(chǎn)生新的固定邊界特征點(diǎn),我們設(shè)fixm表示第m個(gè)縱斷面上的固定特征點(diǎn)。
圖6 典型截面離散化示意圖
將截面優(yōu)化的可行域離散為M×N個(gè)控制節(jié)點(diǎn)后,還需要確定截面每層鈑金的厚度。汽車框架結(jié)構(gòu)大多由鋼板構(gòu)成,雖然鋼板的厚度在理論上是個(gè)連續(xù)變量,但是工程師一般只能從0.7~2.5mm中間選取一些標(biāo)準(zhǔn)的厚度,因而車身關(guān)鍵截面的鈑金厚度實(shí)際上是一個(gè)離散變量。
若截面的鈑金層數(shù)為I,則增加I個(gè)虛擬的縱斷面。假設(shè)截面的第i(i=1,2,…,I)層鈑金可選厚度為Ti種,則可用虛擬縱斷面i上的第ti(ti=1,2,…,Ti)個(gè)節(jié)點(diǎn)來對應(yīng)截面第i層鈑金的第ti種厚度。
在每個(gè)縱斷面和虛擬縱斷面上都放置1只螞蟻,共有M+I(xiàn)只螞蟻,每只螞蟻具有以下特征:
(1)每只螞蟻只能在該縱斷面或虛擬縱斷面的各個(gè)節(jié)點(diǎn)間移動(dòng),它根據(jù)概率函數(shù)(以每個(gè)節(jié)點(diǎn)吸引強(qiáng)度為變量)選擇該縱斷面的下一個(gè)節(jié)點(diǎn)(設(shè)τmn為第m個(gè)縱斷面上第n個(gè)節(jié)點(diǎn)的吸引強(qiáng)度,τiti為第i個(gè)虛擬縱斷面上第ti個(gè)節(jié)點(diǎn)的吸引強(qiáng)度)。
(2)每只螞蟻具有記憶功能,除非一個(gè)循環(huán)完成,否則它不能移動(dòng)到該縱斷面或虛擬縱斷面上已選擇過的節(jié)點(diǎn)。設(shè)tabum為第m個(gè)縱斷面上螞蟻的禁忌表(tabum(n)表示禁忌表中對應(yīng)縱斷面m上第n個(gè)節(jié)點(diǎn)的元素),它用以記錄螞蟻m在一個(gè)循環(huán)中已選擇過的節(jié)點(diǎn),tabum在每個(gè)循環(huán)中隨著優(yōu)化過程做動(dòng)態(tài)調(diào)整。
(3)螞蟻每完成一個(gè)循環(huán)后,在訪問過的節(jié)點(diǎn)上留下信息量。
節(jié)點(diǎn)吸引強(qiáng)度的規(guī)則定義如下:
蟻群算法初始化時(shí),若每個(gè)節(jié)點(diǎn)的吸引強(qiáng)度都相同,則螞蟻很難在短時(shí)間內(nèi)找到較好的解,嚴(yán)重影響收斂速度。解決的方法是在初始化時(shí)就給出啟發(fā)性的信息量[8]:先產(chǎn)生大量的截面形狀(例如50個(gè)截面形狀),從中選擇比較優(yōu)的(例如5個(gè)截面形狀),使這些節(jié)點(diǎn)上留下相應(yīng)的信息量,以此來引導(dǎo)螞蟻進(jìn)行節(jié)點(diǎn)選擇。
第m個(gè)縱斷面上的螞蟻選擇該縱斷面第n個(gè)節(jié)點(diǎn)的移動(dòng)概率Pmn為
式中,allowedm表示螞蟻m下一步可以選擇的節(jié)點(diǎn)表,allowedm= {0,1,…,N}-tabum。
第i個(gè)虛擬縱斷面上的螞蟻選擇該虛擬縱斷面第ti個(gè)節(jié)點(diǎn)的移動(dòng)概率Piti為
縱斷面M上的螞蟻移動(dòng)一次為一個(gè)循環(huán),其他縱斷面上的螞蟻移動(dòng)I次為一個(gè)循環(huán);每個(gè)虛擬縱斷面上的螞蟻移動(dòng)一次為一個(gè)循環(huán)。
當(dāng)所有的螞蟻完成一次循環(huán)后,各節(jié)點(diǎn)的吸引強(qiáng)度根據(jù)下式更新:
式中,ρ為吸引強(qiáng)度的持久性系數(shù),一般0.5≤ρ≤0.9;Δτmn、Δτiti分別為螞蟻在本次循環(huán)中留在節(jié)點(diǎn)上的信息量;Q為一個(gè)正常數(shù);F為目標(biāo)函數(shù)值。
利用節(jié)點(diǎn)生成截面形狀的規(guī)則定義如下:當(dāng)所有的螞蟻完成一次循環(huán)后,所有縱斷面的禁忌表元素組成了截面的全部形狀控制點(diǎn)。坐標(biāo)原點(diǎn)作為截面的一個(gè)焊接點(diǎn),縱斷面M的禁忌表元素tabuM(n)作為截面的另外一個(gè)焊接點(diǎn);剩下的每個(gè)縱斷面禁忌表中的元素tabum(n)按照n的大小排序,n最大(即z坐標(biāo)最大)的節(jié)點(diǎn)為截面的最上層鈑金控制點(diǎn),n最小的節(jié)點(diǎn)為截面的最下層鈑金控制點(diǎn),依此類推確定截面各層鈑金的控制點(diǎn),這樣就保證了截面上層鈑金的每個(gè)控制節(jié)點(diǎn)都在下層鈑金控制點(diǎn)之上,避免了出現(xiàn)各層鈑金交叉的現(xiàn)象;對于每層鈑金上的所有控制點(diǎn),按照縱斷面m的大小順序(即y坐標(biāo)大小順序)依次連接起來就形成了該層鈑金的形狀,該方法保證了每層鈑金不會(huì)出現(xiàn)沖壓負(fù)角和沒有拔模角的現(xiàn)象。
當(dāng)所有的螞蟻完成一次循環(huán)后,虛擬縱斷面i上被選擇的節(jié)點(diǎn)所代表的厚度表示了截面i的厚度。
本文截面優(yōu)化的蟻群算法基本步驟如下:
(1)nc←0(nc為循環(huán)次數(shù));給τnm、τiti賦予相同的值,產(chǎn)生大量的截面形狀(例如50個(gè)截面形狀),從中選擇比較優(yōu)的(例如5個(gè)截面形狀),使這些節(jié)點(diǎn)留下相應(yīng)的信息量,給出ρ、Q的值。
(2)將M+I(xiàn)只螞蟻分別置于所有縱斷面和所有虛擬縱斷面的頂部節(jié)點(diǎn)N 和Ti上;清空tabum。
(3)螞蟻m按照轉(zhuǎn)移概率Pmn選擇該縱斷面中的下一個(gè)節(jié)點(diǎn)n,將節(jié)點(diǎn)n置于tabum中,每只螞蟻選擇I次完成一次循環(huán);螞蟻i按照轉(zhuǎn)移概率Piti選擇該虛擬縱斷面中的下一個(gè)節(jié)點(diǎn)ti,每只螞蟻選擇一次完成一次循環(huán)。
(4)計(jì)算截面的目標(biāo)函數(shù)F,按照更新方程修改節(jié)點(diǎn)吸引強(qiáng)度τmn和τiti;nc←nc+1;
(5)若nc大于規(guī)定的循環(huán)次數(shù),停止運(yùn)算,選擇縱斷面節(jié)點(diǎn)中τmn最大的前I個(gè)節(jié)點(diǎn)為該縱斷面的最終控制節(jié)點(diǎn),選擇虛擬縱斷面節(jié)點(diǎn)中τiti最大的節(jié)點(diǎn)為第i層鈑金的厚度;否則轉(zhuǎn)到步驟(2)。
現(xiàn)以筆者參與開發(fā)的國內(nèi)某新車型的門檻梁截面的設(shè)計(jì)為例來檢驗(yàn)程序的可靠性與實(shí)用性。該新車型門檻梁的基本特征和上文中提到的某車型門檻梁截面較為相似,故這里不再詳述該新車型門檻梁截面具體的離散化過程和算法。該門檻梁截面基本結(jié)構(gòu)由門檻外板、門檻加強(qiáng)板和門檻內(nèi)板組成,每層鋼板可供選擇的厚度級別有6種,厚度為0.7~2.0mm。表1列舉了該門檻梁截面外部輸入的設(shè)計(jì)約束邊界特征點(diǎn)坐標(biāo)。本文采用EXCEL軟件進(jìn)行了截面離散化的二次程序開發(fā),表2列舉了截面離散化過程中的參數(shù)取值以及部分離散化結(jié)果。
表1 門檻梁截面邊界特征點(diǎn)坐標(biāo)
表2 截面離散化的參數(shù)取值及結(jié)果
蟻群算法的各參數(shù)如下:ρ=0.95;Q=30;初始的τnm=τiti=1,首先產(chǎn)生50個(gè)初始截面形狀,從中選取5個(gè)較優(yōu)的形狀,在相應(yīng)的節(jié)點(diǎn)處賦予吸引強(qiáng)度50;規(guī)定的最大循環(huán)次數(shù)為1000。
截面的目標(biāo)函數(shù)F可用下式表示:
式中,wp(p=1,2,3,4)為截面各個(gè)特性的加權(quán)系數(shù),為相應(yīng)截面特性的指定值。
F越小表明優(yōu)化結(jié)果和指定目標(biāo)差距越小,且截面的重量越輕。在汽車的概念設(shè)計(jì)階段,建立由真實(shí)接頭模型、梁單元和大尺寸殼單元組成的車身概念有限元模型,以車身整體剛度和接頭部位強(qiáng)度為約束條件,對車身進(jìn)行輕量化優(yōu)化,可以得到車身各關(guān)鍵截面特性參數(shù)的合理組合,即車身各個(gè)截面特性的指定值,這里不再詳述,具體數(shù)值見表3。
表3 門檻梁截面指定截面特性和優(yōu)化結(jié)果
圖7給出了目標(biāo)值函數(shù)隨迭代過程的變化趨勢,圖8給出了最后得到的截面優(yōu)化形狀,每層鋼板的厚度優(yōu)化結(jié)果見表4,對應(yīng)的截面特性見表3。如圖8所示,截面中的某些鈑金出現(xiàn)了毛刺現(xiàn)象(目前針對毛刺的處理辦法主要靠工程師的經(jīng)驗(yàn)來修正,從算法上來消除毛刺現(xiàn)象還需要做進(jìn)一步的研究),結(jié)構(gòu)工程師根據(jù)工程經(jīng)驗(yàn)去除掉鈑金上的毛刺,參考生成的截面形狀變化趨勢,并增加相應(yīng)的導(dǎo)角,本著簡化零部件形狀的原則,最終設(shè)計(jì)的該門檻梁截面形狀如圖9所示。圖9中,實(shí)線為門檻梁截面各鈑金的輪廓線,虛線為該斷面處其他零部件的輪廓線,門檻梁各層鋼板的厚度同時(shí)也采用了優(yōu)化結(jié)果。采用該算法生成的車身關(guān)鍵截面形狀對后續(xù)的設(shè)計(jì)工作起到了較好的指導(dǎo)作用。
圖7 目標(biāo)值函數(shù)值隨迭代過程的變化趨勢
圖8 生成的門檻梁截面形狀示意圖
表4 門檻梁截面各層鋼板厚度優(yōu)化結(jié)果
圖9 門檻梁截面的最終設(shè)計(jì)結(jié)果示意圖
(1)提出了一種全新的截面形狀生成思路和方法,采用網(wǎng)格法的思想應(yīng)用蟻群算法生成具有指定截面特性且節(jié)省材料的截面形狀,同時(shí)還可得出每層鋼板的厚度。
(2)充分考慮到車身造型、內(nèi)部空間以及總體布置等方面的尺寸約束,同時(shí)滿足了截面鈑金件沖壓成形等方面的形狀約束要求,更加符合工程應(yīng)用的需要。
(3)生成的截面形狀包含了兩個(gè)焊接邊連線和沖壓方向不平行的情況,丟失的截面形狀解較少。
本文算法未考慮截面上孔洞的存在以及零件的導(dǎo)角,對存在沖壓負(fù)角(需要整形工序)的復(fù)雜截面形狀仍不適用。目前截面的變量不是連續(xù)變量。上述問題都是后續(xù)工作的重點(diǎn),這里不再詳述。
[1] Kojima Y.First Order Analysis as CAE for Design Engineers[J].R&D Review of Toyota CRDL,2001,37(1):1-3.
[2] Banichuk N V.Optimization of Elastic Bars in Torsion[J].International Journal of Solids and Structures,1976,12(4):275-286.
[3] Egner W,Zyczkowski M.Optimal Plastic Design of a Bar under Combined Torsion,Bending and Shear[J].Structural and Multidisciplinary Optimization,2001,22(5):394-406.
[4] Banichuk N V,Ragnedda F,Serra M.Optimum Shape of Bar Cross-sections[J].Structural and Multidisciplinary Optimization,2002,23(3):222-232.
[5] Kim Y Y,Kim T S.Topology Optimization of Beam Cross Sections[J].Int.J.Solids Struct.,2000,37(3):477-493.
[6] Masataka Yoshimura,Shinji Nishiwaki,Kazuhiro Izui.A Multiple Cross-sectional Shape Optimization Method for Automotive Body Frames[J].Journal of Mechanical Design,2005,127:49-57.
[7] 覃建周,吳伯杰.汽車覆蓋件沖壓方向的優(yōu)化[J].模具工業(yè),2002,2:5-8.
[8] 高尚,楊靜宇.群智能算法及其應(yīng)用[M].北京:中國水利水電出版社,2006.