孔婕，邓璐璐，刘亿，闵峻英

轻合金精密成形专题

基于GISSMO断裂准则的6016铝合金断裂行为研究

孔婕1，邓璐璐2，刘亿1，闵峻英1

（1. 同济大学 机械与能源工程学院，上海 201804；2. 泛亚汽车技术中心有限公司，上海 201206）

研究零部件在成形与碰撞过程中，6016铝合金在不同应力状态下的断裂行为。通过准静态拉伸实验，获得了6016铝合金的基本力学性能。利用Nakajima成形极限实验，获得了6016铝合金材料的断裂成形极限曲线。设计了7种涵盖成形及碰撞过程中应力状态的断裂极限测试试样，采用数字图像相关技术（DIC）记录了试样在变形过程中的全场应变。利用实验-有限元反求方法标定了6016铝合金的GISSMO断裂准则的参数，并用帽形件三点弯曲实验验证了模型的合理性。相比于传统断裂成形极限图的预测结果，基于GISSMO断裂准则的仿真结果与实验具有更好的一致性。所建立的GISSMO模型可以用于预测6016铝合金在复杂应力状态下的断裂行为。

GISSMO；断裂极限；6016铝合金；失效准则

近些年来，中国提出在2030年前与2060年前分别实现“碳达峰”与“碳中和”的目标。在汽车车身制造中使用汽车轻量化材料是实现“双碳”目标的重要途径，铝合金因其塑性好、密度小、比强度高等优点被广泛应用于车身零部件的制造[1]。随着汽车零件制造工艺的不断发展，材料的断裂问题在成形过程中变得尤为突出，在碰撞过程中材料的断裂行为也更加复杂。断裂极限作为评价板材成形和碰撞性能的重要指标，是指导汽车零部件设计及优化的重要依据，因此对材料断裂极限的精准预测显得尤为重要。

成形极限图（Forming Limit Diagram，FLD）最早由Keeler[2]提出，经过M. G. Goodwin[3]进一步补充和完善，成为判断板料成形性能最简便和最直观的方法。FLD可以定量衡量材料的极限变形能力，有助于分析薄板成形及碰撞过程中的断裂失效问题。Nakajima和Marciniak实验是测定板材成形极限常用的实验方法，可以获得线性应变路径下的断裂极限数据，被工业和学术界广泛采用。Graf等[4]发现2008-T4铝合金板材在几种典型预应变加载条件下的成形极限曲线与无预应变时测量得到的成形极限曲线相比发生了很大的变化，材料的成形极限表现出应变路径依赖性。材料在成形及碰撞的过程中通常会经历复杂应变路径，因此使用成形极限图对断裂材料的断裂行为进行预测并不准确。金属材料断裂准则为预测材料的断裂行为提供了有力途径，根据材料在变形过程中损伤的累积是否与材料的塑性行为相耦合，可以将断裂准则分为耦合与非耦合2种。耦合断裂准则的2个重要分支是连续损伤力模型（Continuum Damage Mechanics，CDM）和以Gurson模型[5]为代表的细观力学损伤模型。非耦合断裂准则假设材料在变形过程中的损伤累积对材料的本构模型参数没有影响，且形式简单、参数标定容易，因此在金属塑性成形领域得到了广泛应用[6]。Neukamm等[7]在2009年提出了一种广义增量应力状态相关模型（Generalized Incremental Stress State Dependent Model，GISSMO），GISSMO非耦合韧性断裂准则考虑了不同应力状态下材料失效行为的差异及非线性损伤累积方式[8]，已被广泛应用于预测成形与碰撞过程中的材料损伤与失效行为[9-12]。

文中以6016铝合金为研究对象，通过单向拉伸和表征不同应力状态的韧性断裂实验，分别获得了材料基本力学性能参数及断裂极限。基于数字图像相关技术（Digital Image Correlation，DIC）获得的实验数据，采用实验-有限元参数反求方法标定了GISSMO韧性断裂准则的参数，通过帽形件三点弯曲实验进行了验证，并对比了传统成形极限方法与GISSMO韧性断裂准则的预测精度，验证了该韧性断裂准则的合理性，为金属板材失效行为的预测提供一定的指导。

1 材料与实验方法

1.1 基本力学性能实验

实验采用厚度为0.8 mm的6016铝合金板材，根据GB/T 228.1—2010，在板材上取与轧制方向成0°、45°、90°的试样进行准静态单向拉伸实验，拉伸试样尺寸如图1所示。拉伸试样采用线切割方式进行加工并使用砂纸打磨消除加工产生的毛刺。准静态拉伸实验在MTS E45.105-ATBC万能材料实验机上进行，拉伸速度为3 mm/min，名义应变速率约为0.001 s−1，拉伸直至试样断裂为止。

图1 准静态单向拉伸试样尺寸

1.2 成形极限实验

Nakajima实验方法[13]常用来评价板材的成形性能，被广泛应用于工业和学术界。实验试样通常为哑铃形（见图2），通过改变试样尺寸可以得到不同应力状态下材料的成形极限。根据ISO标准，外圆直径为180 mm的试样尺寸如表1所示。为了表征6016铝合金的成形性能，对哑铃形的试样（尺寸见表1）进行Nakajima成形极限实验。试样使用激光切割方式沿板材轧制方向切下，利用YAW4605-K电液伺服成形极限实验机对6016铝合金进行成形极限实验，实验冲头半径为50.8 mm，冲头上升速率为20 mm/min，每组实验进行3次保证重复性，采用DIC技术记录试样变形过程中的全场应变，实验装备示意图如图3所示。

图2 用于成形极限实验的哑铃形试样

表1 哑铃形试样尺寸参数

Tab.1 Geometry parameters of dumbbell specimens

图3 成形极限应变测试系统

1.3 韧性断裂测试

1.3.1 应力状态的表征

Bao等[14]利用不同几何尺寸试样进行实验，量化了应力状态对铝合金2024-T351断裂等效应变的影响，在韧性断裂的表征中，应力三轴度通常在断裂模型中用以表征应力状态。应力三轴度的定义是试样静水应力与等效应力的比值：

1.3.2 试样设计及韧性断裂实验

为了测定6016铝合金在不同应力状态下的断裂极限，文中设计了7种不同受力状态的试样进行准静态韧性断裂测试，其中6种试样的关键位置尺寸如图4所示[15]，而等双向拉伸实验采用直径为180 mm的圆盘形试样，实验在MTS电液伺服成形极限机上进行，实验冲头速度为20 mm/min，其余试样则在MTS E45.105-ATBC万能材料实验机上进行准静态拉伸，拉伸速度为1 mm/min。使用DIC技术记录试样变形过程中的全场应变，每种实验进行3次以保证重复性原则。

图4 不同应力状态韧性断裂试样关键尺寸

2 实验结果

2.1 单向拉伸实验结果与硬化模型标定

经单向拉伸实验数据处理后，6016铝合金的基本力学性能如表2所示。

材料3个方向的真实应力-应变曲线如图5所示，由于材料成形时的应变远大于单向拉伸实验的应变范围，因此需要根据实验得到的真实应力-应变曲线外推更大应变范围内的真实应力-应变曲线。

表2 6016铝合金基本力学性能

Tab.2 Basic mechanical property of 6016 aluminum alloy

图5 6016铝合金在不同方向的真实应力-应变曲线

Barnwal等[16]使用混合硬化模型准确描述了TRIP1180与DP980材料的硬化行为，结果显示，混合硬化模型综合了饱和与非饱和型硬化模型的优点，能更加灵活地描述材料的硬化行为。因此，文中采用混合Swift-Hockett_Sherby（SHS）硬化模型，外推6016铝合金在0°、45°、90°这3个方向的真实应力-应变曲线。Swift硬化模型见式（4），Hockett_Sherby硬化模型见式（5），混合SHS硬化模型见式（6）。

式中：为Swift模型在混合SHS模型中的权重；p为剔除弹性段数据后的塑性应变值；y为材料的屈服强度；0、、、、s、和为模型参数，其标定结果分别为0.002、0.9、440.29、0.31、329.41、9.38、0.94。

2.2 成形极限实验结果

实验结束后使用Vic-3D 8软件对试样的应变分布做后处理，提取试样断裂前一帧的面内主、次应变值绘制材料的断裂成形极限图。8种试样的等效断裂塑性应变值依次为0.436、0.428、0.408、0.393、0.413、0.493、0.538、0.584，成形极限图如图6所示。

图6 6016铝合金断裂成形极限图

2.3 韧性断裂测试实验结果

不同应力状态下韧性断裂实验得到的载荷-位移曲线如图7所示，每组实验的重复性较好，除等双向拉伸外，各组实验的载荷在断裂前均出现了“软化”现象，表明铝合金的断裂模式为韧性断裂。

图7 不同试样实验载荷-位移曲线

3 GISSMO韧性断裂准则的标定

GISSMO韧性断裂准则与工程中常用的Johnson- Cook断裂模型的不同点在于，该模型考虑了材料从受损、非线性损伤累积到材料断裂失效的全过程，并且能预测材料在不同应力状态下的失效行为，适于分析与预测金属板材碰撞与成形过程中的断裂问题。

3.1 GISSMO韧性断裂准则简介

GISSMO韧性断裂准则主要使用非线性损伤累积方式定义材料的失效行为，如式（7）所示。在唯象断裂模型中，通常定义不与本构模型耦合的外部变量来描述材料的损伤行为。当仿真中单元的值达到1时，在仿真中就会被删除。

GISSMO模型除了能够表征不同应力状态下材料的断裂行为，还将韧性断裂中材料的失稳现象也考虑其中。材料的软化行为由变量定义，如式（8）所示：

当单元的值达到1时，其应力会发生软化现象，软化后的单元应力由式（9）计算得到。

式中：*为发生软化后的应力；为软化前应力；crit为=1时变量的临界值；为衰减系数。

3.2 GISSMO韧性断裂准则参数确定

图8 LS-OPT参数反求流程

在对6016铝合金进行参数反求的过程中，考虑到材料的各向异性，为提高仿真精度，选择LS_DYNA材料库中的MAT_36_3-PARAMETER_BARLAT_NLP材料模型，输入0°、45°、90°这3个方向上外推后的真实应力-应变曲线。为了保证仿真精度，目标区域网格大小划分为0.5 mm。

3.3 有限元仿真与实验结果对比

经LS-OPT软件优化后的GISSMO断裂模型参数如表3所示，对应得到的断裂及失稳曲线如图9所示。图10为使用优化后的6016铝合金GISSMO韧性断裂模型参数进行实验后，各组实验与仿真的载荷-位移曲线对比，仿真峰值力数值与断裂位移数据均与实验结果相吻合，且仿真断裂位移误差与实验值的误差均小于5%，满足工程使用要求，验证了6016铝合金GISSMO模型的有效性。

3.4 网格正则化过程

网格划分情况同样会影响有限元仿真的结果，在整车碰撞与汽车零件成形仿真中，为了兼顾仿真准确性与仿真效率，GISSMO韧性断裂准则通过引入网格依赖性因子来考虑尺寸效应的影响，使同一应力状态、不同网格尺寸的仿真结果相同。

为了校准不同网格尺寸下的网格依赖性因子值，建立了0.5、1、2.5、5、10 mm网格尺寸的单向拉伸实验仿真模型（见图11）。在未对不同网格大小的仿真模型进行归一化处理前，仿真结果如图12a所示，试样断裂位移随网格尺寸的增加而加大，导致仿真累积误差增加。通过参数反求的方式获取相应的因子值（见图13），使不同网格尺寸的模型均能得到与实验一致的载荷-位移曲线，如图12b所示。

表3 6016铝合金GISSMO韧性断裂模型参数

Tab.3 Parameters of GISSMO ductile fracture model for 6016 aluminum alloy

图9 输入GISSMO模型的6016铝合金断裂及失稳曲线

图10 试样仿真与实验载荷-位移曲线结果

图11 不同网格大小的单向拉伸仿真模型

图12 不同网格尺寸下单向拉伸工况载荷-位移曲线

图13 网格依赖性因子曲线

4 GISSMO韧性断裂准则的验证

4.1 帽形件三点弯曲实验

为了验证文中标定的6016铝合金GISSMO韧性断裂准则卡片的准确性，设计了如图14a所示的帽形件用于三点弯曲工况实验验证。使用折弯工艺加工零件，尽量消除初始损伤对零件的影响。加工好的帽形件与垫板焊接在一起后置于三点弯曲实验工装上（见图14b），实验工装的支撑辊子间距为110 mm，冲头与支撑辊子的直径为20 mm，实验时冲头下压速度为4 mm/min。

建立与实验工况相同的帽形件有限元模型对三点弯曲实验结果进行模拟，使用前文标定好的材料模型和断裂模型描述帽形件与垫板的变形行为，支撑辊子与冲头辊子设为刚性体。仿真模型中帽形件网格尺寸为2.5 mm，对变形较大的圆角区域进行网格细化确保仿真精度。

图14 三点弯曲实验零件尺寸与实验装备

4.2 仿真与实验结果对比

图15a为帽形件经三点弯曲实验后的变形情况，由于6016铝合金塑性较好，实验结束并未出现明显的开裂情况。但零件变形后期，其载荷出现了“软化”现象，即零件承载能力降低，图15b为冲头下压与实验相同距离时的仿真结果，与实验较为吻合。图16展示了三点弯曲工况下实验与仿真的载荷-位移对比情况，发现GISSMO模型能够准确预测材料的变形行为，预测结果尤其与材料变形后期时的软化行为较为吻合。

图15 帽形件三点弯曲实验变形情况

图16 基于GISSMO模型的三点弯曲实验载荷-位移实验与仿真对比结果

成形极限曲线同样可以输入到MAT_36材料模型中用以预测材料的失效行为，由于6016铝合金的颈缩现象并不明显，因此输入Nakajima实验得到的断裂极限曲线用以预测三点弯曲实验结果。当仿真模型中单元的历史变量为1，即成形性能指数（Formability Index，FI）数值达到1时，材料即达到成形极限发生断裂。从图17可以看出，与冲头接触的圆角位置单元FI值在下压距离为25.6 mm时已经达到1，即材料已经发生了断裂现象。实验时冲头下压距离为35 mm时，材料尚未发生断裂，仅出现了承载能力的下降。因此，当使用传统成形极限曲线对材料的断裂行为进行预测时，会低估材料的断裂应变值，导致在仿真中发生提前开裂的情况。基于成形极限预测的结果进行零件设计时也会导致材料的浪费，不论是预测精度还是预测效率，GISSMO模型相较于成形极限图有较高的实用性。

图17 基于成形极限曲线的帽形件成形指数分布云图

5 结论

对6016铝合金进行了准静态单向拉伸实验与成形极限实验，得到了6016铝合金基本力学性能参数并对其成形性能做出了评价。设计了7种涵盖成形与碰撞过程中应力状态的韧性断裂试样，使用有限元-实验混合方法反求了GISSMO韧性断裂参数并设计了三点弯曲工况进行了验证，得到的结论如下。

1）提出的GISSMO断裂模型标定方法能够准确描述6016铝合金在不同应变路径下的断裂极限，仿真预测的断裂应变与实验误差在5%以内。

2）基于GISSMO的断裂模型对6016铝合金的断裂行为进行预测相较于采用传统断裂成形极限曲线进行预测具有更高的精度，能够为预测汽车铝合金构件失效行为提供一定的技术参考。

[1] 汤蓉蓉, 唐学帮. 铝合金汽车轻量化与焊接技术探究[J]. 机电工程技术, 2021, 50(12): 105-107.

TANG Rong-rong, TANG Xue-bang. Research on Aluminum Alloy Automobile Lightweight and Welding Technology[J]. Mechanical & Electrical Engineering Technology, 2021, 50(12): 105-107.

[2] KEELER S P. Plastic Instability and Fracture in Sheet Stretched over Rigid Punches[J]. ASM Transactions Quarterly, 1963, 56(11): 25-48.

[3] GOODWIN G M. Application of Strain Analysis to Sheet Metal Forming Problems in the Press Shop[J]. Metallurgical Transactions A, 1993, 24: 380-387.

[4] GRAF A, HOSFORD W. Effect of Changing Strain Paths on Forming Limit Diagrams of AI 2008-T4[J]. Metallurgical and Materials Transactions A, 1993, 24(11): 2503-2512.

[5] GURSON A L. Continuum Theory of Ductile Rupture by Void Nucleation and Growth: Part I-Yield Criteria and Flow Rules for Porous Ductile Media[J]. Journal of Engineering Materials and Technology, 1977, 99(1): 2-15.

[6] 杨婷, 熊自柳, 孙力, 等. 汽车用先进高强钢韧性断裂模型的研究与应用进展[J]. 锻压技术, 2021, 46(1): 10-16.

YANG Ting, XIONG Zi-liu, SUN Li, et al. Research and Application Progress on Ductile Fracture Model of Advanced High-Strength Steel (AHSS) for Automotive[J]. Forging & Stamping Technology, 2021, 46(1): 10-16.

[7] NEUKAMM F, FEUCHT M, HAUFE A. Considering Damage History in Crashworthiness Simulations[C]// Proceedings of the 7th European LS-DYNA Conference, Salzburg, 2009.

[8] 梁宾, 赵岩, 赵清江, 等. 基于Gissmo失效模型的6016铝合金板材断裂行为研究及应用[J]. 机械工程学报, 2019, 55(18): 53-62.

LIANG Bin, ZHAO Yan, ZHAO Qing-jiang, et al. On the Prediction of Failure in 6016 Aluminum Alloy Sheet by Gissmo Damage Model[J]. Journal of Mechanical Engineering, 2019, 55(18): 53-62.

[9] 张骥超, 连昌伟, 韩非. 第三代超高强钢QP1180硬化与失效行为研究[J]. 机械工程学报, 2021, 57: 1-9.

ZHANG Ji-chao, LIAN Chang-wei, HAN Fei. Study on Hardening and Failure Behavior of the 3rdGeneration Ultra-High Strength Steel QP1180[J]. Journal of Mechanical Engineering, 2021, 57: 1-9.

[10] 徐晨阳, 张骥超, 连昌伟. 基于GISSMO损伤模型的DH590高强钢断裂失效行为研究[J]. 塑性工程学报, 2021, 28(6): 68-74.

XU Chen-yang, ZHANG Ji-chao, LIAN Chang-wei. Study on Fracture Failure Behavior of DH590 High Strength Steel Based on GISSMO Damage Model[J]. Journal of Plasticity Engineering, 2021, 28(6): 68-74.

[11] 徐宽, 夏明生, 刘靖宝, 等. 基于Gissmo材料失效准则的圆盘剪剪刃应力模拟[J]. 塑性工程学报, 2021, 28(9): 110-116.

XU Kuan, XIA Ming-sheng, LIU Jing-bao, et al. Simulation of Stress on Shear Blade of Disc Shear Based on Gissmo Material Failure Criterion[J]. Journal of Plasticity Engineering, 2021, 28(9): 110-116.

[12] OTROSHI M, ROSSEL M, MESCHUT G. Stress State Dependent Damage Modeling of Self-Pierce Riveting Process Simulation Using GISSMO Damage Model[J]. Journal of Advanced Joining Processes, 2020, 1: 100015.

[13] NAKAZIMA K, KIKUMA T, HASUKA K. Study on the Formability of Steel Sheets[J]. Yamata Technical Report, 1968, 264: 8517-8530.

[14] BAO Y, WIERZBICKI T. On Fracture Locus in the Equivalent Strain and Stress Triaxiality Space[J]. International Journal of Mechanical Sciences, 2004, 46(1): 81-98.

[15] GORJI M B, MOHR D. Micro-Tension and Micro-Shear Experiments to Characterize Stress-State Dependent Ductile Fracture[J]. Acta Materialia, 2017, 131: 65-76.

[16] BARNWAL V K, LEE S, KIM J, et al. Failure Characteristics of Advanced High Strength Steels at Macro and Micro Scales[J]. Materials Science and Engineering: A, 2019, 754: 411-427.

[17] MOHR D, HENN S. Calibration of Stress-Triaxiality Dependent Crack Formation Criteria: A New Hybrid Experimental-Numerical Method[J]. Experimental Mechanics, 2007, 47(6): 805-820.

[18] ZHANG Zhi-qiang, CUI Yue-jie, YU Ge. Damaged and Failure Characterization of 7075-T6 Al Alloy Based on GISSMO Model[J]. Journal of Mechanical Science and Technology, 2021, 35(3): 1209-1214.

Fracture Behavior of 6016 Aluminum Alloy Based on GISSMO Criterion

KONG Jie1, DENG Lu-lu2, LIU Yi1, MIN Jun-ying1

(1. School of Mechanical Engineering, Tongji University, Shanghai 201804, China; 2. Pan Asia Technical Automotive Center Co., Ltd., Shanghai 201206, China)

The work aims to study the fracture behavior of 6016 aluminum alloy in forming and crashing simulation.The mechanical property of 6016 aluminum alloy was obtained by quasi-static uniaxial tension experiment and the fracture forming limit curves (FFLC) of 6016 aluminum alloy were obtained by conducting Nakajima experiment. Seven types of specimens were designed to obtain the fracture limits under different stress states which commonly occur during forming and crashing. The strain field of the specimens was recorded by digital image correlation (DIC) method. The hybrid experimental–numerical method was used to calculate the parameters of GISSMO criterion. The reasonability of the model was verified by the three-point bending experiment of hat-shaped parts. The results showed that the prediction results of GISSMO model were in better agreement with the experimental results compared with that of the conventional FFLC. The GISSMO model established in this paper is proved to be suitable for predicting the fracture behavior of 6016 aluminum alloy under complex stress states.

GISSMO; fracture limit; 6016 aluminum alloy; fracture criterion

10.3969/j.issn.1674-6457.2022.04.001

TG146.21

A

1674-6457(2022)04-0001-10

2022-01-04

国家自然科学基金（51805375）

孔婕（1998—），女，硕士生，主要研究方向为金属材料断裂表征。

闵峻英（1986—），男，博士，教授，主要研究方向为汽车轻量化与先进成形制造技术。

责任编辑：蒋红晨