夏天,劉瑩,喬鑫
?
基于python程序設計的發(fā)動機艙蓋多學科聯(lián)合優(yōu)化
夏天,劉瑩*,喬鑫
(華晨汽車工程研究院,遼寧 沈陽 110141)
以某SUV車型為研究對象,對發(fā)動機艙蓋進行優(yōu)化設計。采用基于python語言的二次開發(fā)技術(shù),建立了發(fā)動機艙蓋的參數(shù)化模型。選取發(fā)動機艙蓋內(nèi)板加強梁拓撲結(jié)構(gòu)形式的變化、內(nèi)板橫梁的移動、縱梁的移動作為設計變量,以發(fā)動機艙蓋模態(tài)、剛度為約束條件,以重量及行人頭部保護性能為優(yōu)化目標,建立各項性能指標的徑向基神經(jīng)網(wǎng)絡近似模型,采用非支配遺傳算法對發(fā)動機艙蓋進行多學科聯(lián)合優(yōu)化,在保證各項性能滿足要求的前提下,獲得了發(fā)動機艙蓋的最優(yōu)結(jié)構(gòu)。
發(fā)動機艙蓋;python;多學科優(yōu)化;近似模型;非支配遺傳算法
發(fā)動機艙蓋(簡稱發(fā)蓋,以下同)是整車的重要組成部件,起到空氣導流、保護發(fā)動機及管線、保護行人等作用,發(fā)蓋的結(jié)構(gòu)應滿足模態(tài)、剛度、行人保護以及輕量化的要求。其中,發(fā)蓋內(nèi)板加強梁的結(jié)構(gòu)形式對發(fā)蓋的各項性能都有重要影響,內(nèi)板加強梁的設計是發(fā)蓋結(jié)構(gòu)設計的重要問題。
發(fā)蓋的結(jié)構(gòu)設計涵蓋多個學科領域,各個響應之間相互關(guān)聯(lián)、相互影響,因此發(fā)蓋的設計需要多學科聯(lián)合優(yōu)化設計方法。近年來汽車結(jié)構(gòu)的優(yōu)化研究主要集中在車身[1-4]、車門[5-8]等部件中,對發(fā)蓋的優(yōu)化研究比較少。本文采用python語言二次開發(fā)技術(shù),以某處于研發(fā)階段SUV車型發(fā)動機艙蓋為對象,以內(nèi)板加強梁的拓撲結(jié)構(gòu)以及加強梁的位置為變量,綜合考慮發(fā)蓋的模態(tài)、剛度、重量及行人頭部保護性能,對發(fā)蓋進行多學科聯(lián)合優(yōu)化,獲得了發(fā)蓋的優(yōu)化結(jié)構(gòu)。
發(fā)動機艙蓋剛度分析與模態(tài)分析采用Nastran求解,行人保護分析采用LS-DYNA求解。
取發(fā)蓋鎖鉤中心點、左右鉸鏈安裝點的中心點及左側(cè)鉸鏈安裝點,建立局部坐標系,'方向為鉸鏈軸向,'為以上三點所確定的平面的法向。約束鉸鏈安裝孔的6個自由度,約束前側(cè)緩沖塊的局部坐標系的'自由度;在鎖鉤中心點施加局部坐標系的'方向300N的載荷。發(fā)蓋彎曲剛度分析有限元模型如圖1所示。
圖1 發(fā)蓋彎曲剛度模型
約束鉸鏈安裝孔的6個自由度,約束鎖鉤中心點局部坐標系的'方向;前側(cè)左右緩沖塊在局部坐標系方向施加240N?m的扭矩。發(fā)蓋扭轉(zhuǎn)剛度分析有限元模型如圖2所示。
圖2 發(fā)蓋扭轉(zhuǎn)剛度模型
采用Lanczos方法計算發(fā)蓋的前五階自由模態(tài),并設置最小頻率為1Hz,以節(jié)約計算時間。分析結(jié)束后提取一階扭轉(zhuǎn)模態(tài)頻率和一階彎曲模態(tài)頻率。
在行人保護分析模型中,截取整車A柱以前的模型。約束模型的截斷位置和減震器支座的全部自由度;使用兒童和成人頭部碰撞器模型:兒童頭部碰撞器沖擊速度為40km/h,與水平面呈50°角,成人頭部碰撞器沖擊速度為40km/h,與水平面呈65°角;前端半載地面線為Z=-354mm;計算時間為30ms。行人保護分析有限元模型如圖3所示。
圖3 行人保護分析模型
本文通過python語言的二次開發(fā)實現(xiàn)對發(fā)蓋模型的參數(shù)化建模。在優(yōu)化流程中,通過變量來控制內(nèi)板加強梁的拓撲結(jié)構(gòu)的變化、橫梁位置的移動、縱梁位置的移動,從而實現(xiàn)內(nèi)板加強梁結(jié)構(gòu)的改變。
根據(jù)發(fā)蓋的結(jié)構(gòu)特點及設計經(jīng)驗,本文建立了三種內(nèi)板加強梁的拓撲結(jié)構(gòu),如圖4所示。其中,加強梁拓撲結(jié)構(gòu)的變化設置為變量1,橫梁的縱向位置移動設置為變量2(如圖4豎直白色箭頭),縱梁的橫向位置移動設置為變量3(如圖4水平白色箭頭,另一側(cè)縱梁的移動量為-3,不設置為獨立變量)。
圖4 x2、x3為0時各加強梁結(jié)構(gòu)形式
為了實現(xiàn)發(fā)蓋的參數(shù)化,首先需要分別建立三種拓撲結(jié)構(gòu)的內(nèi)板加強梁模型以及發(fā)蓋的基體模型,如圖5、6所示。
圖5 加強梁結(jié)構(gòu)模型
圖6 發(fā)蓋基體模型
其次,編寫python程序。python的實現(xiàn)流程如下:1在發(fā)蓋基體模型中,根據(jù)輸入的1的值,python程序調(diào)用圖5中相應的的內(nèi)板加強梁模型;2根據(jù)輸入的2、3分別在縱向、橫向上移動橫梁、縱梁以及相應的接頭;3創(chuàng)建加強梁與發(fā)蓋外板間的膨脹膠;4在加強梁接頭與基體模型的接口間創(chuàng)建單元,連接發(fā)蓋基體與加強梁;5補充創(chuàng)建新生成的內(nèi)板加強梁與外板間的膨脹膠;6合并基體內(nèi)板與加強梁的單元節(jié)點;7重置加強梁的屬性信息。圖7為python程序設計實現(xiàn)發(fā)蓋參數(shù)化控制的過程。
圖7 python程序設計參數(shù)化實現(xiàn)過程(x1=3)
通過python程序設計建立了發(fā)蓋參數(shù)化模型,利用Isight開放平臺調(diào)用發(fā)蓋基體的CAE模型文件,再在CAE模型文件中調(diào)用python二次開發(fā)程序,進而輸出不同發(fā)蓋結(jié)構(gòu)的求解文件,最后分別調(diào)用Nastran、DYNA求解器計算樣本點的響應值。為了提升優(yōu)化效率,采用“試驗設計-近似模型-優(yōu)化”的優(yōu)化策略,優(yōu)化流程如圖8所示。
圖8 優(yōu)化設計流程圖
定義加強梁的拓撲結(jié)構(gòu)的變化及橫梁、縱梁的位置移動為設計變量,變量取值范圍如表1所示,各個梁的位置變化方向請參見圖4;以發(fā)蓋彎曲剛度、扭轉(zhuǎn)剛度、發(fā)蓋的重量、一階彎曲模態(tài)、一階扭轉(zhuǎn)模態(tài)頻率及行人保護頭部得分為響應;采用優(yōu)化的拉丁超立方方法生成100個樣本。基于Isight的DOE流程如圖9所示。
圖9 基于Isight的DOE流程
表1 優(yōu)化變量信息
采用徑向基RBF近似模型替代仿真模型來建立設計變量與各個響應之間的關(guān)系,其擬合精度可通過決定系數(shù)2來驗證,2值越接近1,表明近似模型有更高的精度。本文中各響應的2值如表2所示。近似模型的誤差如圖10所示??梢娝薪颇P途哂休^高的精度,滿足優(yōu)化要求。
表2 優(yōu)化響應信息
以行人保護頭部得分和發(fā)蓋質(zhì)量為優(yōu)化目標,二者設置相同的權(quán)重;以彎曲剛度、扭轉(zhuǎn)剛度、一階彎曲模態(tài)頻率、一階扭轉(zhuǎn)模態(tài)頻率為約束,對發(fā)蓋進行多目標優(yōu)化。采用非支配遺傳算法NSGA-Ⅱ搜索優(yōu)化解,種群大小設置為100,迭代代數(shù)設置為100,優(yōu)化得到的pareto前沿如圖11所示。
圖11 Pareto最優(yōu)前沿
從pareto前沿解集中選擇一個優(yōu)化解,將變量值帶入有限元模型中,提交到求解器中計算,得到的結(jié)果與優(yōu)化解對比,誤差均在3%以內(nèi),說明優(yōu)化的結(jié)果是合理的。所選擇的的優(yōu)化優(yōu)化方案變量及響應結(jié)果如表3所示,在滿足剛度與模態(tài)的要求下,發(fā)蓋的質(zhì)量為20.809kg,行人保護頭部得分為8.29。
表3 優(yōu)化結(jié)果及驗證
所選取的優(yōu)化的解梁的拓撲結(jié)構(gòu)為第一種形式;橫梁向車前方移動了43.75mm;縱梁向中間方向移動了64.68mm,如圖12所示。
通過python的程序設計,實現(xiàn)了發(fā)蓋的參數(shù)化建模,并進而完成了多學科聯(lián)合優(yōu)化。在滿足發(fā)蓋剛度、模態(tài)的條件下,以內(nèi)板加強梁拓撲結(jié)構(gòu)、梁的位置為優(yōu)化變量,以發(fā)蓋質(zhì)量與行人保護得分值為優(yōu)化目標,對發(fā)蓋進行了多目標優(yōu)化,獲得了pareto前沿解集,并驗證了優(yōu)化解的有效性。通過python的二次開發(fā)實現(xiàn)的參數(shù)化建模,很容易實現(xiàn)程序的擴展,而且變量的變化均是在CAE模型中進行的,不會引入建模誤差。本文的優(yōu)化思路具有重要的工程應用價值。
[1] TORSTENFELT B, KLARBRINGA A. Conceptual optimal design of modular car product families using simultaneous size, shape and topology optimization[J].Finite Elements in Analysis & Design, 2007,43(14):1050-1061.
[2] 王登峰,蔡珂芳,馬明輝,等.基于隱式參數(shù)化模型的白車身輕量化設計[J].汽車工程,2018,40(5):610-616.
[3] 王震虎,周巧英,劉開勇,等.基于響應面模型的白車身多目標輕量化設計[J].中國機械工程,2018,29(1):75-81.
[4] Duddeck F, Zimmer H. Modular car body design and optimization by an implicit parameterization technique via SFE concept[J]. Lecture Notes in Electrical Engineering. 2012,195:413-424.
[5] 高云凱,申振宇,馮兆玄,等.多目標優(yōu)化在車門輕量化設計中的應用[J].吉林大學學報(自然科學版),2017,45(2):275-308.
[6] 朱茂桃,朱彩帆,郭佳歡,等.基于6σ穩(wěn)健性的軋制差厚板車門優(yōu)化設計[J].中國機械工程,2017,28(8):995-1001.
[7] 周定智,王霄,劉會霞,等.基于穩(wěn)健性和多目標優(yōu)化的車頂結(jié)構(gòu)輕量化設計研究[J].機電工程,2015,32(9):1170-1175.
[8] 袁廷輝,倪士林,王霄,等.基于穩(wěn)健性多目標優(yōu)化的微客尾門輕量化設計[J].機電工程,2017,34(8):860-864.
Multidisciplinary Design Optimization of Engine Hood Based on Python Program Design
Xia Tian, Liu Ying*, Qiao Xin
(Brilliance Automotive Engineering Research Institute, Liaoning Shenyang 110141 )
This paper takes a SUV as the study object, and the Multidisciplinary design optimization of the engine hood is conducted.A parametric model for the engine hood is established by using python program design.The topological structure changes of bonnet reinforcement and the movement of cross beam and longitudinal beam are selected as the design variables; and the constraints include the performance of mode, stiffness; the mass and pedestrian head protection performance are taken as the objectives.Then the Radial basis function approximate models are established for every performance.Finally, optimization is performed using NSGA-Ⅱ Algorithm. The optimum structure of the engine hood is obtained while beeping its performance.
engine hood;python; multidisciplinary optimization;approximate model;NSGA-ⅡAlgorithm
U464
A
1671-7988(2019)08-185-04
U464
A
1671-7988(2019)08-185-04
夏天(1985-),男,工程師,碩士,華晨汽車工程研究院,從事汽車結(jié)構(gòu)輕量化研究工作;
劉瑩(1982-),男,高級工程師,碩士,就職于華晨汽車工程研究院,從事結(jié)構(gòu)耐久輕量化研究工作。
10.16638/j.cnki.1671-7988.2019.08.057