陸善彬， 呂 婕， 陳 偉， 馮兆玄， 戚桂悅， 葉 輝

(1. 吉林大學(xué)汽車工程學(xué)院，吉林 長春 130025；2. 南京工程學(xué)院汽車與軌道交通學(xué)院，江蘇 南京 211167)

基于HyperMesh二次開發(fā)的無鉚釘鉚接有限元快速建模

陸善彬1， 呂 婕1， 陳 偉2， 馮兆玄1， 戚桂悅1， 葉 輝1

(1. 吉林大學(xué)汽車工程學(xué)院，吉林 長春 130025；
2. 南京工程學(xué)院汽車與軌道交通學(xué)院，江蘇 南京 211167)

目前采用約束單元 SPR2建立無鉚釘鉚接的簡化有限元模型，建模簡單，仿真精度較高，但HyperMesh 12.0版本不支持SPR2建模，導(dǎo)致了整車建模效率較低，因此本文應(yīng)用Tcl/Tk 語言和 HyperWorks內(nèi)置函數(shù)開發(fā)了整車鉚接SPR2的有限元快速建模模塊。并以某轎車鉚接建模為例，驗(yàn)證了該模塊的可操作性。該模塊可根據(jù)用戶選擇的鉚接點(diǎn)自動識別連接板件，極大地提高了建模效率，同時可有效地避免軟件的升級更新，節(jié)約費(fèi)用，并為其他新型連接的有限元快速建模提供了一種參考方法。

無鉚釘鉚接；二次開發(fā)；快速建模；Tcl/Tk

無鉚釘自沖鉚接[1-2]以其異種金屬連接的優(yōu)勢而被廣泛應(yīng)用在車身結(jié)構(gòu)的連接中。無鉚釘鉚接是一種機(jī)械連接技術(shù)，不需要額外的材料連接，是通過凸模將上板料壓入下板料形成鎖扣進(jìn)行連接。雙層板鉚接斷面如圖1。在整車碰撞模擬中，為了減少仿真計算時間，需要建立無鉚釘鉚接的簡化有限元模型。無論在動態(tài)試驗(yàn)還是靜態(tài)試驗(yàn)中，在十字拉伸工況下，SPR2模型能夠較準(zhǔn)確預(yù)測出失效點(diǎn)，同時在搭接剪切工況下，SPR2模型能模擬出韌性失效模式。因此，SPR2模型作為鉚接的簡化模型較為理想[3]。

圖1 雙層板鉚接的斷面

HyperMesh是一款應(yīng)用廣泛，功能強(qiáng)大的有限元前處理軟件。目前我國大多數(shù)企業(yè)對于HyperMesh還只限于使用階段，而對HyperMesh二次開發(fā)的研究并不是很多[4]。曹文鋼和范超[5]運(yùn)用Tcl/Tk語言和HyperWorks內(nèi)置函數(shù)開發(fā)了CAE流程自動化系統(tǒng)，為產(chǎn)品設(shè)計、分析提供一種快捷方法。孫靜等[6]研究了 HyperWorks 流程自動化系統(tǒng)的開發(fā)過程與關(guān)鍵技術(shù)，利用系統(tǒng)對鋁合金車體進(jìn)行靜態(tài)強(qiáng)度分析。許宇能等[7]利用HyperMesh二次開發(fā)功能開發(fā)了行人保護(hù)模擬自動化前處理系統(tǒng)，實(shí)現(xiàn)自動劃分參考線、確定硬點(diǎn)和定位模型功能。目前HyperMesh 12.0版本可以建立點(diǎn)焊、膠粘、螺栓等多種連接方式的有限元模型，但不支持無鉚釘鉚接SPR2的建模，只能在所有建模工作結(jié)束后，通過文本編輯完成，導(dǎo)致SPR2建模繁瑣，尤其在整車建模環(huán)境下，建模效率極低。

針對SPR2建模繁瑣，整車鉚接效率低等問題，本文應(yīng)用 HyperMesh二次開發(fā)功能開發(fā)整車無鉚釘鉚接SPR2有限元快速建模模塊，從而大大提高了建模的效率，避免頻繁的軟件升級更新，節(jié)省費(fèi)用，并為其他新型連接方式的有限元快速建模提供參考。

1 基于Tcl/Tk的HyperMesh二次開發(fā)技術(shù)

HyperMesh 提供了Tcl接口，擴(kuò)展軟件功能，以實(shí)現(xiàn)用戶的需求。Tcl(工具命令語言)是一種腳本語言，其語法固定，易于閱讀、平臺獨(dú)立且不需要編譯,可以作為一個獨(dú)立的應(yīng)用程序或嵌入式應(yīng)用程序。而Tk是Tcl的圖形界面工具包，用于定義Tcl命令、創(chuàng)建和操縱用戶界面部件[8]。

HyperMesh提供了兩種常見的功能函數(shù)，一種為“*”命令函數(shù)，如*createvector用于向量的創(chuàng)建。另一種為“hm”命令函數(shù)，如hm_getdistance可以測得兩節(jié)點(diǎn)間距離。

2 整車鉚接快速有限元建模模塊的實(shí)現(xiàn)

整車鉚接快速有限元建模模塊包括3個部分：文件導(dǎo)入、SPR2參數(shù)和局部調(diào)整區(qū)域參數(shù)，模塊界面如圖2。

圖2 模塊的界面

2.2 模塊關(guān)鍵技術(shù)

2.2.1 自動識別連接板件

在進(jìn)行整車鉚接有限元建模時，用戶先在鉚接點(diǎn)位置建立節(jié)點(diǎn)。完成所有節(jié)點(diǎn)創(chuàng)建后，用戶定義上層板件，選擇對應(yīng)的鉚接節(jié)點(diǎn)，模塊自動識別下層板件comp的ID號。

利用命令hm_getclosetnode在除上板以外comp的單元中找到離鉚接節(jié)點(diǎn)最近的節(jié)點(diǎn)，該節(jié)點(diǎn)所在的comp即為下板。

2.2.2 局部調(diào)整區(qū)域的確定

為了讓模擬結(jié)果更精準(zhǔn)，需要調(diào)整鉚接點(diǎn)附近的殼單元屬性，一般鉚接成型后板材附近母材會發(fā)生硬化，屈服強(qiáng)度會增加，同時板材會減薄[3]。上層板料減薄和硬化較為明顯，因此將鉚接點(diǎn)上板周圍殼單元的屈服應(yīng)力增大，殼單元厚度減小，作為局部調(diào)整區(qū)域。

紅色文化作為一種資源，既包括物質(zhì)資源，也包括精神資源。紅色文化的物態(tài)類是指中國共產(chǎn)黨在革命斗爭與建設(shè)過程中所形成的革命精神、革命傳統(tǒng)、革命思想等文化的物質(zhì)載體與物態(tài)表現(xiàn)，其表現(xiàn)形式多樣，有在革命戰(zhàn)爭年代為革命事業(yè)浴血奮戰(zhàn)、做出了偉大貢獻(xiàn)，甚至犧牲生命的革命前輩與革命先烈，作為歷史人物，他們的革命事跡也已融入紅色文化；有革命前輩與革命先烈在革命與戰(zhàn)爭年代使用過的物品以及居住、工作過的舊址，這些物品與舊址因其主人也成為紅色文化的重要組成部分；此外，革命與建設(shè)年代發(fā)生的重大事件與重要活動等也屬于紅色的范疇。

根據(jù)用戶選擇鉚接節(jié)點(diǎn)模塊可以自動建立局部調(diào)整區(qū)域，流程如圖3。區(qū)域近似為正方形，區(qū)域大小和壓邊圈尺寸有關(guān)，由用戶定義。模塊將具有相同屬性和材料的局部調(diào)整區(qū)域單元放一個comp中，并新建與其對應(yīng)的屬性和材料，便于用戶編輯。

模塊先計算附近單元的尺寸，根據(jù)用戶定義的局部調(diào)整區(qū)域大小計算出區(qū)域邊長的單元個數(shù)。

當(dāng)區(qū)域邊長的單元數(shù)為偶數(shù)時，利用離鉚接點(diǎn)最近的節(jié)點(diǎn)找到與其相連的單元，命名為偶數(shù)中心單元。

圖3 建立局部調(diào)整區(qū)域流程圖

若單元數(shù)為2，則中心單元為加強(qiáng)區(qū)單元，若單元數(shù)為4，則再尋找與中心單元相連的單元，中心單元以及與其相連的單元則為局部調(diào)整區(qū)域單元。若單元數(shù)為6,8,10,···,以此類推。

當(dāng)區(qū)域邊長單元數(shù)為奇數(shù)時，模塊要求局部調(diào)整區(qū)域個數(shù)不能為 1。若個數(shù)為 1，鉚接點(diǎn)位置可能在局部調(diào)整區(qū)域的邊界，局部調(diào)整區(qū)域的作用不明顯，模塊會彈窗警告。若單元個數(shù)是3，在上板內(nèi)找到與鉚接節(jié)點(diǎn)距離最小的單元即為奇數(shù)中心單元，奇數(shù)中心單元以及與其相連的單元則為局部調(diào)整區(qū)域單元。若單元個數(shù)5，則個數(shù)為3的局部調(diào)整區(qū)域及其相連的單元則為局部調(diào)整區(qū)域單元。若單元數(shù)為7,9,11,···,以此類推。

在搜索奇數(shù)中心單元時，對鉚接點(diǎn)向上板翻邊平面進(jìn)行映射獲得映射節(jié)點(diǎn)，映射節(jié)點(diǎn)所在單元為奇數(shù)中心單元。流程如圖4所示。

圖4 尋找奇數(shù)單元流程圖

判斷映射點(diǎn)是否在單元內(nèi)，如圖5，創(chuàng)建映射點(diǎn)與單元各節(jié)點(diǎn)N1, N2, N3, N4之間向量Vy1, Vy2, Vy3, Vy4，創(chuàng)建節(jié)點(diǎn)N2和N1的向量V21，N2和N3的 向 量 V23， 計 算 V1=Vy1×Vy2, V2=Vy2×Vy3, V3=Vy3×Vy4, V4=Vy4×Vy1, VNt=V21×V23。若映射點(diǎn)在單元內(nèi)及邊界上，則 V1×VNt≥0，V2×VNt≥0，V3×VNt≥0，V4×VNt≥0，否則映射點(diǎn)在單元外。

圖5 判斷映射點(diǎn)在單元內(nèi)

2.2.3 SPR2鉚接模型建立

利用關(guān)鍵字*CONSTRAINT-SPR2建立鉚接模型，卡片信息如圖6所示。

圖6 CONSTRAINT-SPR2卡片信息

MID和SID：分別為鉚釘連接的主從板comp的 ID號，用戶選擇主板，從板由模塊自動搜索獲得。

NSID：鉚接節(jié)點(diǎn)集的ID號。用戶可以一次性選擇多個鉚接節(jié)點(diǎn)來建立SPR2模型，對于連接相同板件的鉚接節(jié)點(diǎn)建立一個節(jié)點(diǎn)set集，以節(jié)點(diǎn)集中第一個節(jié)點(diǎn)ID號命名。

THICK：主從板件厚度之和。根據(jù)主從板comp的ID號，利用命令hm_getentityvalue，可以分別獲得主從板的厚度，并計算兩板件厚度之和。

其余參數(shù)由試驗(yàn)獲得，用戶可以在模塊界面定義，如圖7(a)，雙擊任一參數(shù)，出現(xiàn)彈窗7(b)，用戶填入?yún)?shù)。點(diǎn)擊OK，完成參數(shù)輸入,如圖7(c)，點(diǎn)擊按鈕clear可以清除數(shù)據(jù)。

確定SPR2的所有參數(shù)后，模塊自動導(dǎo)出文件，覆蓋原文件，此時的.K文件(LS_Dyna的輸入文件)已經(jīng)包含局部調(diào)整區(qū)域信息。考慮到 HyperMesh 12.0版本沒有提供直接建立SPR2的模塊，模塊通過讀寫.K文件的方式創(chuàng)建SPR2鉚接模型。模塊利用Open命令打開文件，對文件進(jìn)行逐行讀取，當(dāng)讀取到“*end”時，文件的訪問位置為結(jié)尾處，將訪問位置設(shè)為最后一行首字符處，然后按照圖6卡片中格式將SPR2信息寫入.K文件，完成SPR2模型建立。

圖7 SPR2參數(shù)輸入界面

3 應(yīng)用實(shí)例

利用某轎車的有限元模型進(jìn)行整車無鉚釘鉚接有限元模型的建立，具體步驟如下：

(1) 在HyperMesh中調(diào)用.tcl文件，打開模塊界面，選擇模型所在目錄，導(dǎo)入模型。

(2) 在所有需要鉚接的部件之間建立鉚接節(jié)點(diǎn)，保證鉚接節(jié)點(diǎn)在兩層板的中心位置，整車及局部的鉚接節(jié)點(diǎn)如圖 8所示。選擇鉚接節(jié)點(diǎn)和上板comp的ID號，并在界面上填寫SPR2及局部調(diào)整區(qū)域的參數(shù)，具體參數(shù)如圖9所示。

圖8 整車及局部鉚接點(diǎn)圖

圖9 無鉚釘鉚接建模參數(shù)

(3) 點(diǎn)擊按鈕Connect，完成整車鉚接的SPR2模型建立，共創(chuàng)建730個鉚接點(diǎn)，整車及局部的鉚接建模如圖10，圖中黑色單元為局部調(diào)整區(qū)域。

圖10 整車及局部鉚接建模圖

需要說明的是，本文所選的鉚接點(diǎn)位置只是為了驗(yàn)證本模塊而所建，實(shí)際生產(chǎn)時，鉚接點(diǎn)位置應(yīng)以實(shí)際需要為準(zhǔn)。

4 結(jié) 論

本文系統(tǒng)地介紹了整車鉚接快速有限元建模的開發(fā)流程及關(guān)鍵技術(shù)，用戶選擇鉚接點(diǎn)，輸入SPR2和局部調(diào)整區(qū)域參數(shù)，模塊自動建立整車鉚接有限元模型，解決了HyperMesh無法實(shí)現(xiàn)SPR2建模的問題，提高了SPR2建模的效率。利用模塊實(shí)現(xiàn)某轎車的無鉚釘鉚接的SPR2建模，驗(yàn)證了模塊的有效性。

[1] Mucha J. The analysis of lock forming mechanism in the clinching joint [J]. Materials & Design, 2011, 32(10): 4943-4954.

[2] Jayasekara V, Min K H, Noh J H, Kim M T, Seo J M, Lee H Y, Hwang B B. Rigid-plastic and elastic-plastic finite element analysis on the clinching joint process of thin metal sheets [J]. Metals and Materials International, 2010, 16(2): 339-347.

[3] 周璐瑤, 呂 婕, 陸善彬, 葉 輝. 無鉚釘自沖鉚接有限元簡化模型建模方法研究[J]. 鍛壓技術(shù), 2014, 39(5): 126-131.

[4] 朱金光, 冷 峻, 劉安寧, 張 鵬. CAE分析在企業(yè)技術(shù)創(chuàng)新中的初步探索與實(shí)踐[J]. 農(nóng)業(yè)裝備與車輛工程, 2008(10): 35-37.

[5] 曹文剛, 范 超. 基于HyperWorks的CAE流程自動化系統(tǒng)設(shè)計開發(fā)[J]. 工程圖學(xué)學(xué)報, 2011, 32(1): 16-21.

[6] 孫 靜, 黃雪飛, 李慧萍, 米小珍, 王 楓. 基于HyperWorks的流程自動化系統(tǒng)開發(fā)與應(yīng)用[J]. 鐵路計算機(jī)應(yīng)用, 2012, 21(12): 30-33.

[7] 許宇能, 陸善彬, 葉 輝. 基于HyperMesh二次開發(fā)的汽車行人保護(hù)數(shù)值仿真自動化前處理系統(tǒng)[C]//Altair 2012 HyperWorks技術(shù)大會論文集, 上海, 2012: 1-9.

[8] Ousterhout J K, Jones K. Tcl/Tk入門經(jīng)典[M]. 2版. 張?jiān)? 譯. 北京: 清華大學(xué)出版社, 2010: 211.

Finite Element Fast Modeling of Clinching Joints Based on HyperMesh Secondary Development

Lu Shanbin1, Lv Jie1, Chen Wei2, Feng Zhaoxuan1, Qi Guiyue1, Ye Hui1
(1. Automotive Engineering College, Jilin University, Changchun Jilin 130025, China; 2. School of Automative&Rail Transit, Nanjing Institute of Technology, Nanjing Jiangsu 211167, China)

At present, the constraint element SPR2 is applied to establish a simplified finite element model of clinching joints for its simple way for modeling and accurate simulation. However, HyperMesh 12.0 version can not support SPR2 modeling, which leads to the low efficiency of vehicle modeling. Therefore the Tcl/Tk and the built-in function of HyperWorks are used to develop finite element fast modeling for the vehicle clinching joint of SPR2. A certain car riveting modeling is taken as an example to verify the availability of the module. The module can automatically identify the connected plates according to the nodes selected, which greatly improves the modeling efficiency, avoids the software updates and saves the cost, and it provides a reference method for other new joints fast modeling.

clinching joint; secondary development; fast modeling; Tcl/Tk

TP 311.52

A

2095-302X(2014)05-0804-05

2014-05-05；定稿日期：2014-06-13

國家自然科學(xué)基金資助項(xiàng)目(51105166)

陸善彬(1978–)，男，江蘇海門人，副教授，工學(xué)博士。主要研究方向?yàn)槠嚺鲎舶踩浴-mail：lusb@jlu.edu.cn