劉瑞軍 胡 平 張 杰

1.北華大學,吉林,132013

2.大連理工大學工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室,大連,116024

3.鄭州宇通客車股份有限公司,鄭州,450016

0 引言

常規(guī)車門設(shè)計主要解決結(jié)構(gòu)尺寸和附件安裝等問題,基本不考慮剛度因素。事實上,車門剛度對車輛的安全性和舒適性有很大的影響[1-5]。隨著計算機技術(shù)的發(fā)展,現(xiàn)在已能夠利用計算機技術(shù)對所設(shè)計車門的剛度進行仿真,將仿真結(jié)果作為優(yōu)選設(shè)計方案的影響因素之一,就可以獲得滿足剛度要求且結(jié)構(gòu)合理、整體輕量、成本低、性能高的設(shè)計方案[6-8],因此,仿真結(jié)果的精度在很大程度上影響車門設(shè)計方案的優(yōu)劣。采用結(jié)構(gòu)分析軟件ABAQUS能夠?qū)崿F(xiàn)車門的常規(guī)剛度仿真,但仿真結(jié)果與試驗結(jié)果總是存在不同程度的差異,如何利用現(xiàn)有的CAE軟件提高車門剛度的仿真精度成了亟待解決的問題。常規(guī)CAE仿真采用的是等厚度、零應力和零應變的板材零件裝配成的車門模型,實際上車門模型的內(nèi)外板經(jīng)過沖壓成形后,上述狀態(tài)值均已發(fā)生了變化,因此,車門仿真模型不能反映車門真實狀態(tài)是導致仿真結(jié)果不精確的主要原因[9-13]。為提高車門仿真精度,本文提出將車門內(nèi)外板的成形結(jié)果(厚度、應力和應變)作為“初始場”引入車剛度仿真模型的精細仿真方法。仿真結(jié)果表明,精細仿真能夠真實地反映車門剛度值。

1 車門剛度

車門剛度是指車門在一定載荷作用下抵抗變形的能力,以載荷與載荷引起變形量的比值來表示,單位N/mm。車門剛度包括側(cè)向剛度、垂直剛度、腰線剛度和表面剛度等。

1.1 車門剛度對車輛性能的影響

表面剛度不足的車門在碰撞時產(chǎn)生的大變形會影響駕駛員和乘客的生存空間,扭轉(zhuǎn)剛度不足引起的車門變形會影響門鎖、升降器等附件的使用性能,嚴重時可導致門鎖和升降器不能正常工作,這些均會影響到乘員的安全性。

整體剛度不足的車門的變形會影響車門的密封性能,無論是漏風、漏水,還是因整體變形引起的噪聲、振動等都會影響到車輛的舒適性。

2 數(shù)據(jù)轉(zhuǎn)換

2.1 數(shù)據(jù)轉(zhuǎn)換規(guī)則

不同的CAD/CAE/CAM軟件采用不同的文件格式表示其模型,例如 *.IGES、*.DXF、*.SAT、*.STL等。為了實現(xiàn)不同 CAD/CAE/CAM 間的數(shù)據(jù)共享,不同的文件格式需要進行數(shù)據(jù)轉(zhuǎn)換。目前,CAD與CAE軟件間的數(shù)據(jù)都是基于零件的幾何模型進行轉(zhuǎn)換的,而基于幾何模型的數(shù)據(jù)轉(zhuǎn)換存在諸多不足:①文件結(jié)構(gòu)復雜、語法眾多,數(shù)據(jù)的讀、寫、存儲和轉(zhuǎn)換困難;②模型的參數(shù)化造型特征經(jīng)過轉(zhuǎn)換后將完全喪失,特征樹結(jié)構(gòu)將不存在;③實體模型轉(zhuǎn)換后,有時只存在線框或面模型,需要重新構(gòu)造面和實體;④復雜模型轉(zhuǎn)換后往往會產(chǎn)生數(shù)據(jù)丟失或冗余,面與面之間還可能產(chǎn)生縫隙;⑤轉(zhuǎn)換后的模型不能進行參數(shù)化修改。鑒于CAE軟件均采用單元號、節(jié)點號和節(jié)點坐標表示有限元網(wǎng)格模型,本文提出了基于CAE網(wǎng)格模型的數(shù)據(jù)轉(zhuǎn)換方法。網(wǎng)格模型是指劃分有限元網(wǎng)格后的幾何模型,網(wǎng)格模型特征通過節(jié)點信息(節(jié)點編號、三維坐標值和轉(zhuǎn)動分量)進行描述,根據(jù)靜力等效原則或虛功原理將網(wǎng)格模型中單元的位移、約束和載荷移置到節(jié)點上,所以,網(wǎng)格模型轉(zhuǎn)換規(guī)則的實質(zhì)是節(jié)點的轉(zhuǎn)換,節(jié)點的轉(zhuǎn)換可以實現(xiàn)網(wǎng)格模型的“零”失真轉(zhuǎn)換。

2.2 車門模型轉(zhuǎn)換

為了達到將板料成形結(jié)果加載到結(jié)構(gòu)分析軟件的目的,需要將自主開發(fā)的板料成形CAE軟件KMAS/one—step產(chǎn)生的成形結(jié)果文件(*.nas),通過編寫的模型轉(zhuǎn)換接口轉(zhuǎn)換為結(jié)構(gòu)分析軟件ABAQUS能夠識別的文件(*.inp),轉(zhuǎn)換流程如圖1所示。

以長安汽車集團開發(fā)的某新車型右前門內(nèi)板CAD模型為例,根據(jù)網(wǎng)格模型轉(zhuǎn)換規(guī)則將KMAS/one—step生成的網(wǎng)格數(shù)據(jù)導入ABAQUS并顯示,顯示結(jié)果如圖 2所示。KMAS/one—step網(wǎng)格數(shù)據(jù)文件和轉(zhuǎn)換后的ABAQUS網(wǎng)格數(shù)據(jù)文件中的單元數(shù)量、單元形態(tài)和相應節(jié)點坐標均相同,說明按照網(wǎng)格模型轉(zhuǎn)換規(guī)則實現(xiàn)網(wǎng)格零件單元、節(jié)點數(shù)據(jù)轉(zhuǎn)換的方法可行、正確。

3 車門剛度精細仿真

為了提高車門剛度仿真精度和計算效率,車門內(nèi)外板的光順區(qū)域采用四邊形單元化分網(wǎng)格,邊、角等復雜區(qū)域采用三角形單元化分網(wǎng)格。另外,在劃分網(wǎng)格時,需對車門邊框、車窗部位和密封圈接觸面等狹長結(jié)構(gòu)進行局部網(wǎng)格加密處理。采用公用節(jié)點相連的方式完成車門內(nèi)外板的裝配,使車門成為由連續(xù)單元構(gòu)成的整體,便于進行結(jié)構(gòu)分析。

結(jié)合車門結(jié)構(gòu)特點及技術(shù)要求,根據(jù)車門剛度性能校核標準,確定車門剛度精細仿真分析的約束和加載。設(shè)UX、UY、UZ分別為車門在X 、Y 、Z軸方向的位移量,URX、URY、URZ分別為車門繞X、Y、Z軸的轉(zhuǎn)動量。車門剛度精細仿真流程如圖3所示。

4 車門剛度精細仿真實例

以長安汽車開發(fā)的整體式車門為例,采用精細仿真方法對車門側(cè)向剛度進行分析。車門側(cè)向剛度試驗約束和加載形式見表1。

表1 車門側(cè)向剛度試驗約束和加載形式

車門板采用四邊形和三角形殼單元,鉸鏈采用八節(jié)點六面體和六節(jié)點楔形體單元,車門劃分網(wǎng)格后共有29 845個單元、30 783個節(jié)點。車門內(nèi)外板材料為ST14,材料性能參數(shù)見表2。

表2 車門內(nèi)外板材料性能參數(shù)表

采用ABAQUS軟件對車門分別進行常規(guī)剛度仿真和精細剛度仿真,車門側(cè)向剛度的常規(guī)仿真和精細仿真結(jié)果如圖4、圖5所示。常規(guī)仿真、精細仿真與試驗的最大應力、最大位移和側(cè)向剛度見表3。分析表3數(shù)據(jù)可知,常規(guī)仿真的誤差為14.28%,精細仿真的誤差為3.45%,因此,精細仿真能夠較真實地反映車門的剛度。

表3 車門的最大應力、最大位移和側(cè)向剛度

圖4與圖5的顯示結(jié)果表明精細仿真的最大應力值大于常規(guī)仿真的最大應力值,說明精細仿真可以直觀地反映出應力集中部位。精細仿真的最大位移小于常規(guī)仿真的最大位移,說明板材的沖壓成形結(jié)果具有抵抗車門整體變形的能力。

為了觀察關(guān)鍵區(qū)域(應力集中區(qū))單元的位移變化情況,按照總載荷10%的幅值遞增加載,以更加精確地觀察成形結(jié)果對局部變形的影響程度,關(guān)鍵單元(單元號為25 053)的位置如圖6所示。

外部載荷以20N為基本單位遞增,最終載荷為200N,記錄關(guān)鍵單元的位移情況。單元位移隨外部載荷遞增時的變化情況如圖7所示。

由圖7可知,隨著載荷的增大,關(guān)鍵單元在精細仿真中的位移明顯大于常規(guī)仿真的位移。結(jié)果表明車門板沖壓成形后的應變時效硬化特性有利于提高車門整體的側(cè)向剛度,但板材變薄也會導致出現(xiàn)局部大變形現(xiàn)象,因此,應根據(jù)精細仿真結(jié)果對車門大應變區(qū)進行相應的強化設(shè)計。

5 結(jié)束語

本文提出了在CAE分析前處理過程中將板材成形結(jié)果作為“初始場”替代常規(guī)“零場”的精細仿真方法,該方法克服了常規(guī)車門剛度仿真模型不能反映實際狀態(tài)的弊端。采用網(wǎng)格轉(zhuǎn)換規(guī)則實現(xiàn)車門內(nèi)外板成形結(jié)果加載的精細仿真方法,明顯提高了車門剛度的仿真精度。

[1]Yuichi K,Chinmoy P.Evaluation of Vehicle Body Stiffness and Strength for Car to Car Compatibility[C]//International Body Engineering Conference&Exposition.Detroit,2001:2001-01-3098.

[2]Spencer C,Vartanian K L,White M A,et al.Light T ruck Frame Joint Stiffness Study[C]//SAE 2003 World Congress&Exhibition.Detroit,2003:2003-01-0241.

[3]Baskar S,Kamat A,Lange R.Analytical Robust DoorHingeSystem Design TaguchiApproach[C]//International Body Engineering Conference&Exposition.Detroit,1998:982306.

[4]Petniunas A,Otto N C,Simpson R.Door System Design for Improved Closure Sound Quality[C]//Noise&Vibration Conference&Exposition.T raverse City,MI,USA,1999:1999-01-1681.

[5]Staif J.Dimensional Management Process Applied to AutomotiveDoorSystem[C]//International Body Engineering Conference&Exposition.Detroit,1998:982305.

[6]施法中.板料沖壓成形過程的計算機仿真[J].中國機械工程,1998,9(11):70-71.

[7]莫健華,韓飛.金屬板材數(shù)字化漸進成形技術(shù)研究現(xiàn)狀[J].中國機械工程,2008,19(4):491-497.

[8]李峰,曹穎,申國哲,等.衛(wèi)星天線沖壓成形仿真分析[J].中國機械工程,2004,15(16):1487-1450.

[9]張君媛,陳建濤,唐建春,等.基于CATIA 平臺的車門附件參數(shù)化布置的研究[J].計算級輔助設(shè)計與圖形學學報,2003,24(6):765-768.

[10]高云凱,楊欣,金哲峰,等.轎車車身剛度優(yōu)化方法研究[J].同濟大學學報,2005,33(8):1094-1097.

[11]梁亮,韓立強,邵揚,等.材料性能及其結(jié)構(gòu)分析對車門設(shè)計的影響[J].機械設(shè)計與制造,2007(1):70-71.

[12]陳國華.基于有限元的汽車車門靜態(tài)強度剛度計算與分析[J].機械制造與自動化,2006,35(6):21-24.

[13]趙立紅,楊玉英,孫振忠.汽車覆蓋件剛度評價指標的研究[J].材料科學與工藝,2004,12(4):425-428.