劉玉敏,陳桂華
(漯河職業(yè)技術(shù)學(xué)院,河南 漯河 462000)
汽車車身包含大量的焊點(diǎn),在白車身結(jié)構(gòu)有限元分析過(guò)程中,焊點(diǎn)的模擬精度關(guān)系著白車身及其零部件強(qiáng)度分析和疲勞分析的精度。 受計(jì)算機(jī)軟硬件條件的制約,目前在白車身結(jié)構(gòu)疲勞耐久仿真分析過(guò)程中,仍無(wú)法建出精細(xì)的焊點(diǎn)模型[1],Liang Jianyong 等人的研究表明,ACM2 焊點(diǎn)模型在白車身模態(tài)分析中,有限元分析結(jié)果的誤差可控制在8%以內(nèi)[2]。 本文針對(duì)ACM(RBE3-HEXARBE3)這種焊點(diǎn)類型,對(duì)簡(jiǎn)單的薄壁梁結(jié)構(gòu)進(jìn)行彎曲剛度和扭轉(zhuǎn)剛度分析計(jì)算,找出焊點(diǎn)的最優(yōu)模擬方式,為結(jié)構(gòu)有限元仿真分析中焊點(diǎn)的模擬方法提供參考。
圖1 RBE3-HEXA-RBE3 焊點(diǎn)單元的表征示意圖
ACM 焊點(diǎn)即 Area contact model 焊點(diǎn)類型,如圖 1 所示,此種焊點(diǎn)單元類型會(huì)在兩個(gè)焊接面之間生成一個(gè)六面體實(shí)體單元,然后該六面體實(shí)體單元的8 個(gè)節(jié)點(diǎn)通過(guò)REB2 單元與焊接面進(jìn)行連接。 載荷通過(guò)RBE3 單元分配到被焊接單元的節(jié)點(diǎn)上,分配的大小由RBE2 單元的權(quán)重系數(shù)決定。 如圖1 所示,對(duì)于六面體單元的每一個(gè)節(jié)點(diǎn)均位于殼單元平面內(nèi),節(jié)點(diǎn)的載荷和位移相應(yīng)的由殼單元通過(guò)RBE3 單元加權(quán)平均獲得,從而實(shí)現(xiàn)兩層殼單元之間的載荷和位移的傳遞[3]。
為了對(duì)比及計(jì)算方便,采用如圖2 所示的簡(jiǎn)單梁結(jié)構(gòu)作為基礎(chǔ)模型,對(duì)不同的焊點(diǎn)類型進(jìn)行模擬驗(yàn)證[4]。 該梁?jiǎn)卧蓛蓚€(gè)厚度均為1mm 的鈑金焊接而成,將該梁?jiǎn)卧W(wǎng)格大小平均設(shè)置為8mm,在其兩側(cè)分別各均布10 個(gè)焊點(diǎn)。 該梁長(zhǎng)為380mm、寬為87mm、高為35.5mm,焊點(diǎn)間距均為35mm,均勻分布,共20 個(gè)焊點(diǎn)單元,最終建好的模型共含有1 426 個(gè)網(wǎng)格單元和1 426 個(gè)節(jié)點(diǎn)。
圖2 梁?jiǎn)卧P?/p>
RBE2 單元是剛性單元,具有極高的剛度,為了對(duì)比不同焊點(diǎn)單元類型間的剛度,以RBE2 單元為基準(zhǔn),將(ACM)剛度仿真分析值比RBE2 單元?jiǎng)偠确抡娣治鲋?,以此結(jié)果來(lái)評(píng)判這種焊點(diǎn)類型的剛度值。
具體評(píng)價(jià)標(biāo)準(zhǔn): 評(píng)價(jià)標(biāo)準(zhǔn)=(ACM 剛度仿真分析值/RBE2 單元?jiǎng)偠确抡娣治鲋担?。
2.2.1 扭轉(zhuǎn)工況的確定
如圖3 所示,約束端面A 的1-6 個(gè)方向的全部自由度,將端面B 處的節(jié)點(diǎn)用RBE2 抓取,在抓取的MPC 主節(jié)點(diǎn)處施加大小為2 000N·m 的扭矩,中間截面C 不做任何約束,讀取該梁?jiǎn)卧淖畲笮D(zhuǎn)角度。
圖3 工況示意圖
2.2.2 彎曲工況的確定
約束端面A 和端面B 的1-6 個(gè)方向的全部自由度,將中間截面C 處的節(jié)點(diǎn)用RBE2 單元抓取,在RBE2 主節(jié)點(diǎn)處施加-Z 向100N 的力,讀取RBE2 主節(jié)點(diǎn)的Z 向位移。
2.3.1 扭轉(zhuǎn)工況分析結(jié)果
按照扭轉(zhuǎn)分析工況的設(shè)置,在ABAQUS 中對(duì)不同工況下的模型進(jìn)行計(jì)算,最終得出分析結(jié)果。 根據(jù)圖4 的分析結(jié)果,得出在ACM 焊點(diǎn)及RBE2 下該梁結(jié)構(gòu)一端的扭轉(zhuǎn)角度分別為 0.011 9°、0.010 6°。 由此計(jì)算出不同焊點(diǎn)的扭轉(zhuǎn)剛度具體如表1 所示。
圖4 扭轉(zhuǎn)工況下分析結(jié)果
表1 不同焊點(diǎn)扭轉(zhuǎn)剛度分析結(jié)果
由表1 可知,ACM 焊點(diǎn)下梁模型的扭轉(zhuǎn)剛度為RBE2單元類型下梁?jiǎn)卧まD(zhuǎn)剛度的89.08%。
2.3.2 彎曲工況分析結(jié)果
彎曲工況下分析結(jié)果如圖5 所示。 對(duì)ACM 焊點(diǎn)及RBE2 單元下梁結(jié)構(gòu)模型進(jìn)行有限元分析,提取中間截面C 處RBE2 主節(jié)點(diǎn)在Z 方向上的位移,位移大小分別為0.047 7mm 和0.0443mm。 因此,計(jì)算出梁的彎曲剛度如表2 所示。
圖5 彎曲工況下分析結(jié)果
表2 不同焊點(diǎn)彎曲剛度分析結(jié)果
根據(jù)以上分析結(jié)果,ACM 焊點(diǎn)下梁模型的彎曲剛度約為RBE2 單元下剛度的92.87%。
1) 彎曲剛度和扭轉(zhuǎn)剛度方面,ACM 焊點(diǎn)的剛度大于CBAR 焊點(diǎn)的剛度,CBAR 焊點(diǎn)的剛度大于CWELD 的剛度。
2)有限元仿真分析過(guò)程中,不同焊點(diǎn)類型對(duì)結(jié)構(gòu)的扭轉(zhuǎn)剛度及彎曲剛度影響差異較大,應(yīng)該綜合考慮焊點(diǎn)的建模形式、建模時(shí)間以及焊點(diǎn)的模擬精度等因素,最終結(jié)合相應(yīng)的試驗(yàn)結(jié)果來(lái)決定焊點(diǎn)的模擬形式。