ABAQUS二次開發在動態剪切擠壓 連接仿真中的應用

戎江 馬世博 閆華軍 張雙杰 邵海濤 張超

摘 要：針對在采用有限元軟件分析多參數變化對研究對象的影響時往往需要反復建模、重復定義、修改參數等問題，提出了改進的編程方法和圖形用戶界面建立方法。以中間坯動態剪切擠壓連接為例，基于Python語言對ABAQUS前處理進行了二次開發，編寫了動態剪切擠壓連接仿真模擬的內核腳本程序，借助RSG插件建立了動態剪切擠壓連接的圖形用戶界面，最后進行了仿真分析。結果表明，用戶可根據實際情況輸入不同工藝參數、剪具參數、坯料參數，自動完成前處理操作，在有效解決了模型裝配時繁瑣、易錯等問題的同時還提高了模型裝配和仿真結果的準確性，驗證了所提出的二次開發程序及圖形用戶界面的可行性。研究結果為后續研究不同壓下量、搭接量、刃口寬等參數對中間坯動態剪切連接質量的影響提供了有力工具，同時為ABAQUS二次開發在其他領域的應用提供了借鑒。

關鍵詞：特種加工工藝;運動仿真;Python語言;ABAQUS二次開發;剪切連接

中圖分類號：TG335.19 文獻標識碼：A

doi：10.7535/hbkd.2021yx03014

Applied research of secondary development of ABAQUS in the simulation of dynamic shear extrusion connection

RONG Jiang1， MA Shibo1，2， YAN Huajun1，2， ZHANG Shuangjie1，2， SHAO Haitao1， ZHANG Chao1

（1.School of Material Science and Engineering，Hebei University of Science and Technology，Shijiazhuang，Hebei 050018，China;2.Hebei Key Laboratory of Material Near-net Forming Technology，Shijiazhuang，Hebei 050018，China）

Abstract：Aiming at the problem that when the finite element software is used to analyze the influence of multi parameter changes on the research object，it is necessary to repeatedly model，define and modify parameters，an improved programming method and graphical user interface establishment method were provided.Taking the dynamic shear extrusion connection of intermediate billet as an example，the python language was used for secondary development of ABAQUS pretreatment，and the kernel script of dynamic shear extrusion connection simulation was programmed.Then the graphic user interface of dynamic shear extrusion connection was established by using RSG plug-in.Finally，simulation analysis was conducted.The results show that users can input different process parameters，die parameters and blank parameters according to the actual situation to complete the pre-processing operation automatically，which not only effectively solves the problems of complicated and error prone model assembly，but also improves the accuracy of model assembly and simulation results，so the feasibility of the secondary development program and graphical user interface is verified.The study provides a powerful tool for the follow-up study of the influence of different reduction，overlap，edge width and other parameters on the dynamic shear connection quality of intermediate billet，and provides a reference for the secondary development of ABAQUS in other fields.

Keywords：

special processing technology;motion simulation;Python language;ABAQUS secondary development;shear connection

ABAQUS作為一款大型通用有限元仿真分析軟件，在工程模擬中應用廣泛，可以分析大多數工程材料的線性、非線性行為[1-4]。為了滿足不同用戶的特定需求，ABAQUS軟件提供了基于Python語言的腳本接口[5-6]，以方便用戶開發自定義的前后處理工具，實現參數化建模、數據處理等功能。HUANG等[7]基于ABAQUS二次開發和統計分析軟件SPSS，利用Python腳本實現了空心樓板的自動建模、分析和數據處理。TANG[8]以Python作為編程語言對ABAQUS進行二次開發，并以彎輥為例，對其進行了參數分析。YI等[9]在研究分析齒輪箱減振降噪時，運用Python語言對ABAQUS的前后處理進行二次開發，提高了前后處理的分析和計算效率。

中間坯動態剪切擠壓連接技術是將兩塊搭接在一起的坯料，通過傳送帶輸送至動態剪切擠壓連接裝置的上、下剪具處，在上、下剪具擺動瞬間完成剪切擠壓連接[10-11]。在研究工藝及參數對中間坯連接質量的影響和對工藝進行優化時，往往需要反復建模、重復定義、修改參數等許多繁瑣工作[12-13]。基于數值模擬[14-17]及ABAQUS二次開發參數化建模思路，可以減小前處理過程中建模、裝配、定義載荷等工作量，提高有限元分析效率。因此，本文擬采用Python語言對ABAQUS進行二次開發，編寫動態剪切擠壓連接仿真模擬的內核腳本程序，借助RSG插件建立剪切擠壓連接的圖形用戶界面，實現前處理參數化，用戶可根據具體輸入的工藝參數、模具參數、坯料參數進行動態剪切擠壓連接仿真模擬。此外，通過與非參數化模擬進行對比，驗證該二次開發程序的可行性。

1 動態剪切擠壓連接仿真模型的建立

1.1 剪切擠壓連接簡介

如圖1所示動態剪切擠壓連接初末狀態，實際生產中上、下剪具在剪切連接裝置的作用下，以近橢圓軌跡在豎直面上運動，止料刃與中間坯接觸后限制金屬沿坯料方向的塑性流動，刃口咬入坯料，在上、下剪具刃口作用下搭接坯料發生剪切-擠壓劇烈變形，使中間坯瞬間完成連接[18]，去除坯料接頭后進入精軋機。

1.2 剪切擠壓連接模型的建立與網格劃分

動態剪切擠壓連接作為一種新工藝，研究其工藝參數、剪具參數、坯料參數等對連接質量的影響尤為重要。仿真模擬時一旦有參數發生變化必將重新建模并裝配，這將花費大量的時間和精力，故提出前處理參數化。針對上、下剪具外形輪廓規則，為提高模擬效率將其設置為解析剛體并創建參考點和集合，坯料設置為三維可變形體。上、下剪具為解析剛體無需劃分網格，坯料網格采用六面體結構，坯料在剪切連接過程中連接區的網格會發生較大的扭曲，故單元類型選擇線性減數積分（C3D8RT）?？紤]到連接過程中只是連接區內的網格會發生扭曲，因此將坯料分區為普通區和細分區，

并對細分區進行網格細分，這樣可在保證仿真精度的同時提高模擬效率。坯料網格劃分后如圖2所示，細分區網格數量為24 000，是普通區的5倍。

1.3 材料屬性

材料屬性涉及熱傳導、比熱、彈塑性、膨脹、密度等，其中熱傳導和比熱通過大量試驗測得，

且任一參數的設置有誤都會導致中間坯連接失敗，經過反復模擬及實驗得到可實現連接的材料參數，部分材料屬性見表1。另外，坯料材料定義中合理的損傷準則可以得到準確的仿真模擬結果，動態剪切擠壓連接具有連接速度快、連接區應力應變大、坯料溫度高的特點，故采用Johnson-Cook[19] 損傷模型，見式（1）。

εf=[D1+D2exp（D3η） ]（1+D4ln ε*eq）（1+D5T* ），（1）

式中：D1，D2，D3是受應力三軸度影響的材料參數; D4是應變率敏感參數;D5是溫度影響參數;η是應力三軸度，η=Sm /Seq，其中Sm是平均應力（ 3個主應力的平均值） ，Seq是等效應力;ε*eq是無量綱化應變率，ε*eq=εeq/ε0，其中εeq和ε0分別為等效塑性應變率和參考應變率; T*為無量綱溫度，T* = （ T - Tr） /（ Tm - Tr） ，其中T，Tr和 Tm分別為材料的當前溫度、參考溫度和熔點。

1.4 載荷施加

動態剪切擠壓連接過程中無法通過簡單的加載力和速度實現上、下剪具的復雜運動，因此本文采用幅值加載軌跡，可更準確地模擬剪具真實運動狀態。為此，根據圖3 a）所示動態剪切擠壓連接的機構原理圖[20]設計出如圖3 b）所示的曲柄搖桿和曲柄擺桿機構，通過UG軟件對機構進行運動仿真分析得到如圖4所示上、下剪具的運動軌跡。上、下剪具的運動軌跡可分為擠壓連接行程、退刀行程、空行程，為了節約模擬時間，在ABAQUS中取軌跡的擠壓連接行程進行加載。

為了更真實地模擬實際生產線中的中間坯動態剪切擠壓連接過程，在有限元仿真軟件ABAQUS中將坯料設置為隨動，即坯料與上、下剪具在水平方向上的速度相同，豎直方向上可隨剪具壓力發生移動。上、下剪具除了設置水平、豎直方向的運動，還有繞Z軸的公轉和繞止料刃的自轉。圖5所示為剪切擠壓連接過程部分重要片段。初始時上下剪具前斷面與水平方向呈14°夾角，上、下剪具在水平、豎直方向運動的同時繞止料刃逆時針自轉。止料刃與坯料接觸后限制變形區金屬在水平運動方向上的塑性流動，剪具在豎直方向繼續剪切，刃口開始切除料頭，完成連接后剪具與坯料分離，剪具前斷面與水平方向夾角變為0°。

2 動態剪切擠壓連接仿真模型的參數化

由于動態剪切擠壓連接仿真的前處理建模及其他模塊的參數設置比較復雜，無論是錄制宏文件還是修改“.rpy” [21]文件，得到的程序繁瑣陳雜，不易修改。因此本文采用手動編程，圖6為動態剪切擠壓連接仿真界面的創建流程。通過定義RSG（really simple GUI）對話框構造器開發所需要的function函數，借助RSG中的GUI創建圖形用戶界面。在Kernel頁面加載該腳本程序，并使腳本中的參數與GUI窗體控件中的keywords保持一致，從而產生關聯，注冊生成動態剪切擠壓連接仿真插件程序。

2.1 模型參數化

實際生產中中間坯的厚度、寬度等尺寸往往是變化的，這就要求設計不同的剪具來適應生產，其中對連接質量起主要作用的模具參數是止料刃的高度和刃口寬。需要將中間坯的厚度（thickness—t）、寬度（width—w1）以及剪具止料刃的高度（height—h）和刃口寬（width—w2）設置為變量以供用戶選擇。表2為本文參數化模擬的相關物理量及含義。

2.2 工藝參數的參數化

在動態剪切擠壓連接過程中搭接量（lap amount—L）、壓下量（reduction—R）是影響連接質量的2個關鍵工藝參數，如圖7所示。本文研究搭接量和壓下量對中間坯連接質量的影響時，二者無法直接用參數替換。根據裝配關系分析得出二者由上、下剪具的初始位置所控制，因此將搭接量、壓下量與上、下剪具的初始位置關聯。根據裝配關系進行方程求解，將其轉化為可參數化的變量。

壓下量[10]（R）計算式為

R=st×100%=2t-yt×100%，（2）

式中：s為坯料下壓量，mm;t為中間坯厚度，mm;y為剪切完成后坯料連接區的垂直距離，mm。

Y上=52+52102y，（3）

Y下=50-50102y，（4）

式中：Y上為上剪具的初始縱坐標;Y下為下剪具的初始縱坐標;上剪具在剪切連接過程中豎直方向的行程是52 mm;下剪具在剪切連接過程中豎直方向的行程是50 mm。

搭接量（L）為刃口之間的重疊區域，搭接量越大參與變形的坯料越多。剪切連接過程中認為整體在水平方向上具有相同的速度，故上、下剪具在水平方向上的相對位移是0。如此便可通過裝配經驗先確定下剪具的初始橫坐標，再根據式（5）確定上剪具的初始橫坐標。

X上= X下-28×2+2L，（5）

式中：X上為上剪具的初始橫坐標;X下為下剪具的初始橫坐標;L為搭接量，mm;28是刃口到止料刃的距離，mm。

確定模型初始位置的部分腳本代碼如下：

#下剪具初始位置

a1.translate（instanceList=（′xiamo-1′，），vector=（-25.0，-50-0.98*t+0.0049*t*R，0.0））

#上剪具初始位置

a1.translate（instanceList=（′shangmo-1′，），vector=（-25-56+L，52+1.02*t-0.0051*t*R，0.0））

2.3 其他前處理參數化

軌跡加載的部分程序代碼如下：

#上剪具運動軌跡

DisplacementBC（name=′S-X′，createStepName=′Step-1′，region=region，u1=1.0，u2=UNSET，u3=UNSET，ur1=UNSET，ur2=UNSET，ur3=UNSET，amplitude=′X′，fixed=OFF，distributionType=UNIFORM，fieldName=′′，localCsys=None）

DisplacementBC（name=′S-Y′，createStepName=′Step-1′，

region=region，u1=UNSET，u2=1.0，u3=UNSET，ur1=UNSET，ur2=UNSET，ur3=UNSET，amplitude=′S-Y′，fixed=OFF，distributionType=UNIFORM，fieldName=′′，localCsys=None）

DisplacementBC（name=′S-jiaodu′，createStepName =′Step-1′，region=region，u1=UNSET，u2=UNSET，u3=UNSET，ur1=UNSET，ur2=UNSET，ur3=-0.017，amplitude=′jiaodu′，fixed=OFF，distributionType =UNIFORM，fieldName=′′，localCsys=None）

#材料屬性定義的部分程序代碼如下：

JohnsonCookDamageInitiation （ table=（（0.101，1.68，-2.353，-0.634，-0.255，1425.0，717.0，100.0），））

2.4 圖形用戶界面

通過RSG對話框構造器構建圖形用戶界面。在Kernel標簽頁中選擇剛剛保存的程序文件，function函數選擇程序文件中定義的函數。在GUI標簽頁中對關鍵字進行定義，并與腳本文件中的關鍵字相對應。為了便于用戶理解，還在對話框中添加了動態剪切擠壓連接模型的圖形文件，使整個對話框顯得直觀形象，生成如圖8所示的圖形用戶界面。

3 仿真結果對比

用上述開發的圖形用戶界面，輸入與非參數化模擬相同的參數，圖9、圖10所示分別為參數化模擬和非參數化模擬的應力、應變曲線，從中可以看出二者變形基本一致，測量應變云圖中二者連接后的坯料厚度，代入式（3）得出參數化模擬的結果更接近理論壓下量的值。經多組參數模擬對比得出表3所示結果。從表3可以看出參數化模擬有效降低了模擬時間，提高了仿真效率，同時提高了仿真結果的準確性。

4 結 語

1）基于ABAQUS有限元軟件前處理模塊進行二次開發，編寫了動態剪切擠壓連接仿真模擬的內核腳本程序，借助RSG插件建立了動態剪切擠壓連接的圖形用戶界面，使用戶可以根據實際參數快速建立仿真模型進行模擬，有效解決了模型裝配時繁瑣、易錯等問題，提高了模擬分析效率。

2）運用UG軟件建立了動態剪切機構的三維模型，并裝配進行了運動仿真。通過調整刀塊的裝配位置得到了剪切連接所需的合理軌跡，將其運用到ABAQUS中完成了中間坯動態剪切擠壓連接的模擬分析。

3）針對工藝參數中搭接量、壓下量無法直接用參數替換的問題，根據裝配關系將二者與上、下剪具的初始位置關聯，進行方程求解，轉化為可參數化的變量。通過與非參數化模擬結果進行對比，驗證了該方法的正確性，且模擬結果與靜態試驗中的變形規律具有一致性。

4）本文主要對ABAQUS的前處理進行二次開發，實現了中間坯動態剪切擠壓連接仿真模擬的參數化，后續將對坯料的組織模擬及后處理進行相應開發，以實現整個過程及結果的參數化。

參考文獻/References：

[1] WU Jian，HUANG Xiaochuan.Parametric modeling of CFRP steel tube concrete based on Abaqus secondary development[J].Journal of Progress in Civil Engineering，2020，2（7）：25-28.

[2] 吳澤兵，張帥，郭龍龍，等.ABAQUS二次開發在PDC單齒破巖仿真中的應用[J].西安石油大學學報（自然科學版），2020，35（1）：104-109.

WU Zebing，ZHANG Shuai，GUO Longlong，et al.Application of ABAQUS secondary development in rock breaking simulation of PDC cutter[J].Journal of Xi′an Shiyou University （Natural Science Edition），2020，35（1）：104-109.

[3] 程可慧.基于Python的ABAQUS在特征值屈曲分析上的應用[J].山西建筑，2020，46（19）：196-198.

CHENG Kehui.In the eigenvalue buckling analysis application of ABAQUS based on Python[J].Shanxi Architecture，2020，46（19）：196-198.

[4] 張春會，鄭曉明，田英輝，等.組合荷載作用下平板錨承載能力的數值預測[J].河北科技大學學報，2016，37（3）：283-287.

ZHANG Chunhui，ZHENG Xiaoming，TIAN Yinghui，et al.Numerical prediction of the bearing capacity of plate anchors subjected to combined loads[J].Journal of Hebei University of Science and Technology，2016，37（3）：283-287.

[5] SHIY L，LI H W，WANG W D，et al.A fiber model based on secondary development of ABAQUS for elastic-plastic analysis[J].International Journal of Steel Structures，2018，18（5）：1560-1576.

[6] 朱兆華，黃菊花，張庭芳，等.ABAQUS前、后處理模塊二次開發的應用[J].組合機床與自動化加工技術，2009（1）：30-33.

ZHU Zhaohua，HUANG Juhua，ZHANG Tingfang，et al.Application of Second-developed on ABAQUS pre-precess and post-process[J].Modular Machine Tool & Automatic Manufacturing Technique，2009（1）：30-33.

[7] HUANGC T，WANG Z J，QU Y Y.Parameters analysis with regard to the hollow floor based on the ABAQUS secondary development[J].Advanced Materials Research，2012，446/447/448/449：1337-1342.

[8] WEI T.Application of ABAQUS secondary development in finite element analysis of the bend roller[J].Advanced Materials Research，2011，187：609-613.

[9] YIG L，SUI Y K，DU J Z.Application of python-based Abaqus preprocess and postprocess technique in analysis of gearbox vibration and noise reduction[J].Frontiers of Mechanical Engineering，2011，6（2）：229-234.

[10]張雙杰，李云匣，馬世博，等.金屬快速剪切擠壓連接變形的關鍵參數優化[J].塑性工程學報，2017，24（3）：69-77.

ZHANG Shuangjie，LI Yunxia，MA Shibo，et al.Optimization on key parameters for the metal deformation of rapid shear extrusion bonding[J].Journal of Plasticity Engineering，2017，24（3）：69-77.

[11]何陽.金屬快速剪切擠壓連接的關鍵參數研究[D].石家莊：河北科技大學，2015.

HE Yang.Research on the Key Parameters of Metal Fast Shear Extrusion Connections[D].Shijiazhuang：Hebei University of Science and Technology，2015.

[12]陳春城，胡亞輝，鄭清春，等.基于Python的ABAQUS二次開發在骨鉆削仿真中的應用[J].現代制造工程，2018（6）：54-57.

CHEN Chuncheng，HU Yahui，ZHENG Qingchun，et al.Secondary development of ABAQUS in the simulation of bone drilling based on Python[J].Modern Manufacturing Engineering，2018（6）：54-57.

[13]陳博偉，胡亞輝，楊翠蕾，等.ABAQUS二次開發在切削鈦合金仿真中的應用[J].現代制造工程，2017（7）：105-108.

CHEN Bowei，HU Yahui，YANG Cuilei，et al.Secondary development of ABAQUS in the simulation of cutting of titanium alloy[J].Modern Manufacturing Engineering，2017（7）：105-108.

[14]李爭，程立源，高夢海，等.鸚鵡螺等角螺線型風力機數值模擬及建筑擾流分析[J].河北科技大學學報，2020，41（6）：551-556.

LI Zheng，CHENG Liyuan，GAO Menghai，et al.Numerical simulation of nautilus equiangular spiral wind turbine and analysis of building turbulence[J].Journal of Hebei University of Science and Technology，2020，41（6）：551-556.

[15]張蘇英，郭寶樑，趙國花，等.智能疏散系統導向標志計算及疏散路徑仿真[J].河北科技大學學報，2020，41（1）：40-49.

ZHANG Suying，GUO Baoliang，ZHAO Guohua，et al.Intelligent evacuation system guidance sign calculation and evacuation path simulation[J].Journal of Hebei University of Science and Technology，2020，41（1）：40-49.

[16]劉東輝，王金豪，孫曉云，等.Z字形飛機建模與半實物仿真研究[J].河北科技大學學報，2019，40（6）：520-525.

LIU Donghui，WANG Jinhao，SUN Xiaoyun，et al.Modeling and hardware-in-the-loop simulation of Z-shaped aircraft[J].Journal of Hebei University of Science and Technology，2019，40（6）：520-525.

[17]吳復濤，王娟，李鑫，等.AlSi9Mg變速箱下殼低壓鑄造工藝模擬及優化[J/OL].河北科技大學學報.[2020-12-25].http：//xuebao.hebust.edu.cn/hbkjdx/ch/reader/create_pdf.aspx？file_no=b2020yx008&flag=1&journal_id=hbkjdx&year_id=.

WU Futao，WANG Juan，LI Xin，et al.Simulation and optimization of low pressure casting process for AlSi9Mg gearbox shell[J/OL].Journal of Hebei University of Science and Technology.[2020-12-25].http：//xuebao.hebust.edu.cn/hbkjdx/ch/reader/create_pdf.aspx？file_no=b2020yx008&flag=1&journal_id=hbkjdx&year_id=.

[18]馬世博，侯瑞東，邵明杰，等.固態金屬剪切連接組織及界面演變研究[J].塑性工程學報，2017，24（2）：204-211.

MA Shibo，HOU Ruidong，SHAO Mingjie，et al.Microstructure and interface evolution of solid metal shear bonding[J].Journal of Plas-ticity Engineering，2017，24（2）：204-211.

[19]許鑫.基于快速剪切擠壓變形Q345鋼本構關系研究[D].石家莊：河北科技大學，2015.

XU Xin.Based on the Rapid Shear Extrusion Deformation Q345 Steel Constitutive Relation Research[D].Shijiazhuang：Hebei University of Science and Technology，2015.

[20]姚云峰，李嶺沖，張雙杰，等.中間坯快速剪切擠壓連接機構特性分析與優化[J].機械設計，2016，33（12）：85-89.

YAO Yunfeng，LI Lingchong，ZHANG Shuangjie，et al.Characteristic analysis and optimization of fast shear-extrusion joining mechanism for intermediate billets[J].Journal of Machine Design，2016，33（12）：85-89.

[21]LUO L L，ZHAO M Y.The secondary development of ABAQUS by using python and the application of the advanced GA[J].Physics Procedia，2011，22：68-73.