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

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

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

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

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

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

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

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

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

圖1 雙層板鉚接的斷面

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

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

1 基于Tcl/Tk的HyperMesh二次開發技術

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

HyperMesh提供了兩種常見的功能函數，一種為“*”命令函數，如*createvector用于向量的創建。另一種為“hm”命令函數，如hm_getdistance可以測得兩節點間距離。

2 整車鉚接快速有限元建模模塊的實現

整車鉚接快速有限元建模模塊包括3個部分：文件導入、SPR2參數和局部調整區域參數，模塊界面如圖2。

圖2 模塊的界面

2.2 模塊關鍵技術

2.2.1 自動識別連接板件

在進行整車鉚接有限元建模時，用戶先在鉚接點位置建立節點。完成所有節點創建后，用戶定義上層板件，選擇對應的鉚接節點，模塊自動識別下層板件comp的ID號。

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

2.2.2 局部調整區域的確定

為了讓模擬結果更精準，需要調整鉚接點附近的殼單元屬性，一般鉚接成型后板材附近母材會發生硬化，屈服強度會增加，同時板材會減薄[3]。上層板料減薄和硬化較為明顯，因此將鉚接點上板周圍殼單元的屈服應力增大，殼單元厚度減小，作為局部調整區域。

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

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

模塊先計算附近單元的尺寸，根據用戶定義的局部調整區域大小計算出區域邊長的單元個數。

當區域邊長的單元數為偶數時，利用離鉚接點最近的節點找到與其相連的單元，命名為偶數中心單元。

圖3 建立局部調整區域流程圖

若單元數為2，則中心單元為加強區單元，若單元數為4，則再尋找與中心單元相連的單元，中心單元以及與其相連的單元則為局部調整區域單元。若單元數為6,8,10,···,以此類推。

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

在搜索奇數中心單元時，對鉚接點向上板翻邊平面進行映射獲得映射節點，映射節點所在單元為奇數中心單元。流程如圖4所示。

圖4 尋找奇數單元流程圖

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

圖5 判斷映射點在單元內

2.2.3 SPR2鉚接模型建立

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

圖6 CONSTRAINT-SPR2卡片信息

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

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

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

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

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

圖7 SPR2參數輸入界面

3 應用實例

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

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

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

圖8 整車及局部鉚接點圖

圖9 無鉚釘鉚接建模參數

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

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

需要說明的是，本文所選的鉚接點位置只是為了驗證本模塊而所建，實際生產時，鉚接點位置應以實際需要為準。

4 結 論

本文系統地介紹了整車鉚接快速有限元建模的開發流程及關鍵技術，用戶選擇鉚接點，輸入SPR2和局部調整區域參數，模塊自動建立整車鉚接有限元模型，解決了HyperMesh無法實現SPR2建模的問題，提高了SPR2建模的效率。利用模塊實現某轎車的無鉚釘鉚接的SPR2建模，驗證了模塊的有效性。

[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]. 鍛壓技術, 2014, 39(5): 126-131.

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

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

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

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

[8] Ousterhout J K, Jones K. Tcl/Tk入門經典[M]. 2版. 張元章, 譯. 北京: 清華大學出版社, 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

國家自然科學基金資助項目(51105166)

陸善彬(1978–)，男，江蘇海門人，副教授，工學博士。主要研究方向為汽車碰撞安全性。E-mail：lusb@jlu.edu.cn