李佳霖,李勁松,李美
(海南大學(xué)機(jī)電工程學(xué)院,海南 ???570228)
在傳統(tǒng)產(chǎn)品結(jié)構(gòu)設(shè)計(jì)和優(yōu)化中,在獲得最終方案之前,通常需要根據(jù)方案參數(shù)的變化進(jìn)行多次模擬分析和實(shí)驗(yàn)驗(yàn)證。這是一個(gè)典型的循環(huán),名為“設(shè)計(jì)評(píng)估重新設(shè)計(jì)”的過(guò)程[1]。這個(gè)周期一直持續(xù)到它得到最佳方案或達(dá)到預(yù)定目標(biāo)為止。此外,在每個(gè)周期中,變更方案通常是基于經(jīng)驗(yàn)的,這在一定程度上影響了效率,其結(jié)果的可靠性也受到質(zhì)疑。為了消除人工干預(yù)的瓶頸,必須將仿真、分析和優(yōu)化過(guò)程相結(jié)合。計(jì)算機(jī)科學(xué)的蓬勃發(fā)展和計(jì)算機(jī)軟件技術(shù)的不斷完善,也為產(chǎn)品結(jié)構(gòu)設(shè)計(jì)和優(yōu)化中提供了強(qiáng)大的計(jì)算機(jī)平臺(tái)。例如,Isight平臺(tái)是一個(gè)數(shù)據(jù)集成的軟件框架,由四個(gè)應(yīng)用數(shù)學(xué)算法(測(cè)試設(shè)計(jì)、近似建模、勘探優(yōu)化和質(zhì)量設(shè)計(jì))有機(jī)地組成,是一個(gè)計(jì)算機(jī)自動(dòng)、智能化、重復(fù)性的設(shè)計(jì)過(guò)程,以實(shí)現(xiàn)更好、更快的生產(chǎn)目標(biāo)。
Isight平臺(tái)由于其強(qiáng)大的數(shù)據(jù)集成和數(shù)據(jù)處理功能[2],主要用于基于正交數(shù)組、拉丁超立方體等方法進(jìn)行DOE(實(shí)驗(yàn)設(shè)計(jì)),分析結(jié)構(gòu)參數(shù)變化對(duì)目標(biāo)函數(shù)的影響或使用實(shí)現(xiàn)產(chǎn)品結(jié)構(gòu)的多目標(biāo)優(yōu)化設(shè)計(jì),對(duì)優(yōu)化算法進(jìn)行了探討,目前,在Isight平臺(tái)的幫助下,陳曉平等人對(duì)影響曲軸疲勞強(qiáng)度的關(guān)鍵結(jié)構(gòu)參數(shù)進(jìn)行了敏感性分析、主要效應(yīng)分析和相互作用分析;程成等人完成了基于DOE方法的螺旋槳敞水性能設(shè)計(jì)與優(yōu)化研究,提高了螺旋槳的效率和最小壓力系數(shù);尹波等人采用兩種遺傳算法對(duì)基于MDO方法的三維機(jī)翼氣動(dòng)結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計(jì)。
然而,在有關(guān)Isight平臺(tái)應(yīng)用的文獻(xiàn)中,缺乏對(duì)其基本功能和相關(guān)概念的描述,也缺乏對(duì)其強(qiáng)大的數(shù)據(jù)集成和數(shù)據(jù)處理模塊的介紹。此外,由于與Isight平臺(tái)相關(guān)的出版書(shū)籍較少,用戶很難掌握自己的理論知識(shí)和應(yīng)用方法。因此,本文重點(diǎn)介紹了Isight平臺(tái)的數(shù)據(jù)集成和數(shù)據(jù)處理功能,并給出了一個(gè)工業(yè)實(shí)例。
本文的其余部分包括以下幾個(gè)部分:第一章節(jié)概述了Isight平臺(tái)的數(shù)據(jù)集成;第二章節(jié)介紹了平臺(tái)中的數(shù)據(jù)處理方法;第三章節(jié)介紹了產(chǎn)品結(jié)構(gòu)設(shè)計(jì)和優(yōu)化的整個(gè)過(guò)程。最后,第四章節(jié)得出結(jié)論并提出建議。
一般來(lái)說(shuō),在產(chǎn)品結(jié)構(gòu)設(shè)計(jì)和優(yōu)化過(guò)程中,通常需要使用一些軟件。因此,對(duì)于每組數(shù)據(jù),計(jì)算或分析的過(guò)程非常繁瑣。如果數(shù)據(jù)量很大,這個(gè)過(guò)程將持續(xù)很長(zhǎng)時(shí)間,這將花費(fèi)大量的人力和時(shí)間。然而,Isight平臺(tái)為其他軟件和自編程提供通用集成接口(如simcode、高級(jí)解析)和專業(yè)集成接口(如matlab、abaqus、ansys…)。然后,通過(guò)命令調(diào)用集成軟件或自編程,實(shí)現(xiàn)計(jì)算或分析的自動(dòng)操作[3]。
Isight中的數(shù)據(jù)集成策略是使仿真程序能夠輕松耦合。在Isight平臺(tái)中,如圖1所示,有兩種不同的集成方法供用戶選擇。第一種是通用集成接口,第二種是專業(yè)集成接口。由于第一種方法的靈活性和實(shí)用性,目前得到了廣泛地應(yīng)用。對(duì)于第一種方法,如圖2所示,積分方法主要有3個(gè)步驟[4-6]:
圖1 Isight平臺(tái)兩種集成方法
(1)對(duì)于simcode中的輸入或輸出模塊,應(yīng)集成來(lái)自軟件或自編程的輸入或輸出計(jì)算文件;
(2)對(duì)于simcode中的命令模塊,應(yīng)集成能夠?qū)崿F(xiàn)集成軟件獨(dú)立運(yùn)行的命令;
(3)通過(guò)圖形文件解析器從輸入和輸出文件映射特定參數(shù),如圖2所示的紅色標(biāo)記是輸入?yún)?shù),所示的綠色標(biāo)記是輸出參數(shù)。
圖2 文件分析器用戶界面
與第一種方法相比,第二種積分方法基本相同,但沒(méi)有第二步。
通過(guò)文件集成接口實(shí)現(xiàn)仿真程序的輸入和輸出映射到Isight后,由Isight驅(qū)動(dòng)產(chǎn)品設(shè)計(jì)和優(yōu)化工作,實(shí)現(xiàn)自動(dòng)化探索。
基于 Isight平臺(tái)的產(chǎn)品結(jié)構(gòu)設(shè)計(jì)與優(yōu)化,其強(qiáng)大的數(shù)據(jù)處理功能不容忽視。Isight的設(shè)計(jì)與優(yōu)化探索工具主要包括DOE技術(shù)、優(yōu)化算法、近似方法等。其中,DOE技術(shù)和優(yōu)化算法使用最廣泛。
DOE是指在明確測(cè)試目標(biāo)的前提下,對(duì)測(cè)試因子的具體設(shè)計(jì)和安排、各因子的水平和實(shí)驗(yàn)次數(shù),可以用來(lái)確定一個(gè)或多個(gè)因素對(duì)研究問(wèn)題特征的影響,并對(duì)變化效果進(jìn)行統(tǒng)計(jì)評(píng)價(jià)。DOE中有三個(gè)步驟:實(shí)驗(yàn)設(shè)計(jì)、執(zhí)行實(shí)驗(yàn)和結(jié)果分析。
對(duì)于實(shí)驗(yàn)設(shè)計(jì),Isight平臺(tái)提供了七種供用戶選擇的方法[1,7],如圖3所示。在這七種方法中,每種方法都有其優(yōu)缺點(diǎn)。針對(duì)具體問(wèn)題,用戶應(yīng)根據(jù)效率、準(zhǔn)確性和靈活性的原則選擇方法。
在結(jié)果分析方面,基于集成軟件的實(shí)驗(yàn)分析數(shù)據(jù),DOE后處理主要提供主效應(yīng)和交互效應(yīng)計(jì)算,并對(duì)數(shù)據(jù)進(jìn)行帕累托分析。因此,用戶可以得到主效應(yīng)圖、交互圖和帕累托圖,如圖4所示。此外,DOE后處理還為用戶提供了系數(shù)表、相關(guān)圖和工程數(shù)據(jù)挖掘。
圖3 操作界面和實(shí)驗(yàn)方法介紹
“主效應(yīng)”只考慮一個(gè)自變量對(duì)因變量的影響,而不考慮其他自變量。“交互效應(yīng)”提供了幾個(gè)自變量對(duì)因變量的協(xié)同效應(yīng)?!芭晾弁蟹治觥盵8]表示影響結(jié)果的自變量百分比(貢獻(xiàn)率);在pareto圖中,藍(lán)色條表示相應(yīng)的效果為正,紅色表示負(fù)。
圖4 后處理的主要圖表
圖5 16種優(yōu)化算法
在Isight平臺(tái)中,由于其強(qiáng)大的數(shù)據(jù)集成和處理能力,最廣泛使用的功能除了DOE就是MDO,其主要是對(duì)多個(gè)子目標(biāo)實(shí)現(xiàn)綜合性優(yōu)化,稱為多目標(biāo)優(yōu)化。然而,目標(biāo)可能相互沖突,因此優(yōu)化過(guò)程必須受到影響。因此,Isight平臺(tái)中有16種可用的優(yōu)化算法,供用戶選擇更好的算法和更合適的算法,以便通過(guò)標(biāo)量或非標(biāo)量方法實(shí)現(xiàn)產(chǎn)品結(jié)構(gòu)的優(yōu)化。如圖5所示,為16種優(yōu)化算法[9-10]。
本節(jié)介紹基于Isight平臺(tái)的數(shù)據(jù)集成和數(shù)據(jù)處理實(shí)現(xiàn)汽車橫向穩(wěn)定桿優(yōu)化的工業(yè)實(shí)例。
如圖6所示,汽車水平穩(wěn)定的具體參數(shù)如下:lt為760mm;l1為278mm;d為30mm;d為14mm;r為20mm;結(jié)構(gòu)截面為空心圓柱,D為截面外徑,d為截面內(nèi)徑,l1為杠桿臂長(zhǎng)[11]。
圖6 汽車橫向穩(wěn)定桿
根據(jù)汽車水平穩(wěn)定結(jié)構(gòu)輕量化和應(yīng)力最小化的優(yōu)化原理,本例中有三個(gè)輸入?yún)?shù)和三個(gè)輸出參數(shù),如表1所示。對(duì)于輸出參數(shù),Smax和Umax是結(jié)構(gòu)在工作狀態(tài)下的最大應(yīng)力和最大位移,v是汽車水平穩(wěn)定的體積。
表1 輸入和輸出參數(shù)
圖7 集成流程圖
本案例利用Isight平臺(tái)的通用集成接口,將CATIA軟件和ABAQUS軟件集成在一起,實(shí)現(xiàn)了汽車橫向穩(wěn)定的DOE和MDO。綜合流程圖如圖7所示。首先在CATIA中建立實(shí)體模型,對(duì)輸入?yún)?shù)進(jìn)行參數(shù)化后,得到.macro文件。然后將實(shí)體模型導(dǎo)入到ABAQUS軟件中,進(jìn)行有限元分析,得到模型工作狀態(tài)下的相關(guān)輸出參數(shù)以及.py文件。最后,基于平臺(tái)操作界面進(jìn)行了文件集成操作[12]。
3.4.1 DOE數(shù)據(jù)處理
本案例采用最優(yōu)拉丁超立方體設(shè)計(jì)方法來(lái)實(shí)現(xiàn)DOE。在設(shè)置輸入?yún)?shù)范圍和測(cè)試方案后,如圖8所示,DOE自動(dòng)計(jì)算,并根據(jù)結(jié)果生成帕累托圖[13]。其中,本案例的測(cè)試方案自動(dòng)生成20組,滿足測(cè)試需求的最少數(shù)量。
圖8 DOE的預(yù)處理界面
根據(jù)圖9中S、U和V的帕累托圖,可以得知輸入?yún)?shù)對(duì)輸出參數(shù)的貢獻(xiàn)率。例如,輸入?yún)?shù)d對(duì)v有正作用,對(duì)s和u有負(fù)作用,借助于帕累托圖,它們也可以指導(dǎo)我們進(jìn)行產(chǎn)品優(yōu)化研究。
圖9 DOE試驗(yàn)的帕累托圖
3.4.2 MDO數(shù)據(jù)處理
本案例使用NCGA的方法來(lái)執(zhí)行MDO[19,20]。在設(shè)置輸入?yún)?shù)、輸出參數(shù)和測(cè)試方案的范圍后,如圖10所示,MDO自動(dòng)計(jì)算,并根據(jù)結(jié)果生成Pareto集。在遺傳算法的設(shè)置中,初始種群數(shù)為8,代數(shù)為15。結(jié)果表明,MDO有120種不同的方案。
圖10 MDO預(yù)處理界面
從MDO的帕累托解集中,如圖11所示,可以得到更好的優(yōu)化,或者用戶可以根據(jù)實(shí)際情況選擇最佳的方案。在這種情況下,第三十九個(gè)方案是最終的優(yōu)化方案,優(yōu)化前后汽車水平穩(wěn)定性的數(shù)據(jù)比較如表2所示。
圖11 MDO帕累托解集
表2 優(yōu)化前后數(shù)據(jù)對(duì)比
相比之下,該方案的smax、umax和v值都有不同程度地降低,優(yōu)化后的方案達(dá)到了預(yù)期的目標(biāo)。優(yōu)化前后的位移云圖和應(yīng)力云圖如圖12和圖13所示。
圖12 位移云圖
圖13 應(yīng)力云圖
本文彌補(bǔ)了Isight軟件基礎(chǔ)知識(shí)在學(xué)術(shù)領(lǐng)域的空白,為那些想利用Isight軟件進(jìn)行產(chǎn)品設(shè)計(jì)和優(yōu)化研究的研究者提供了最基本的概念和方法,也為同類研究提供了一種思考和參考的方法。筆者相信,基于Isight平臺(tái)強(qiáng)大的功能,它將得到越來(lái)越廣泛的應(yīng)用。