《Mechanics of Solid Polymers》5.3.8 Yeoh超弹模型

Jiayu Zhou 1, Zhaodong Xia 1, Dongfang Ma 2,*, Huanran Wang 1Materials (Basel). 2024 Apr 10;17(8):1728. doi: 10.3390/ma17081728摘要PC/ABS 复合材料是安全气囊盖板的常用材料，本文对 PC/ABS 复合材料进行了不同温度和应变速率下的单轴拉伸试验。结果表明，温度和加载速率对 PC/ABS 复合材料的力学性能有影响，随着温度的升高，屈服应力减小，断裂时刻的应变增大，但同一温度下应变速率对力学性能影响较小，与延性材料相似。将试验结果代入考虑热效应的 Abaqus 模型，采用逆方法计算得到准确的 Johnson-Cook 本构参数。基于本构模型和 DIC 获得的失效分析结果，模拟了室温和不同应变速率下的单轴拉伸试验，并与试验结果进行了比较，准确再现了试验过程。设计了 PC/ABS 复合材料靶板压入试验，建立了有限元模型模拟试验过程。研究结果与试验结果进行了对比，表明本构断裂应变和失效断裂应变计算结果正确。关键词： 失效分析、聚合物和塑料、Johnson-Cook 本构、数值模拟、反演方法 本文通过单轴拉伸试验及弹道冲击试验标定了PC/ABS复合材料的Johnson-Cook模型塑性及失效参数，取得较好的效果。本文获得的PC/ABS 材料参数如下： 密度 1120 kg/m3 比热容 机械-热转化β 0.9 摩擦系数 0.25 杨氏模量 1750MPa 泊松比 0.38 A 57.5MPa d1 0.825 B 120MPa d2 -0.043 C 0.032 d3 2.6 n 1.734 d4 / m 1.02 d5 / 一 、拉伸试验 实验采用型号为TC-45M的PC/ABS复合材料，试样设计如图1b所示，并在MTS-810动静万能材料试验机（美国MTS公司）上进行测试（图1a）。该设备应变率范围为10⁻⁴–10 s⁻¹，配备温控箱以实现不同温度下的中低应变率拉伸试验。高速拉伸时需考虑试样尺寸对力平衡误差的影响，因此选用15mm的标距长度，以减小尺寸对实验结果的影响，同时适用于准静态实验。 由于PC/ABS复合材料广泛应用于汽车制造（如安全气囊盖），其在高应变率下的失效行为是研究重点，因此额外选用ZwickRoell-5020高速液压拉伸试验机系统（图1c，应变率范围10–1000 s⁻¹）进行高速冲击测试。试样尺寸统一采用图1b所示规格（标距15 mm），具体实验条件如表1所列，涵盖应变率（0.01–1000 s⁻¹）和温度（238.15 K–293.15 K）的多组合工况。 Table 1.单轴拉伸试验条件Strains (s−1) Temperature (K) Equipment Model Loading Rate (mm/s) 0.01 238.15/273.15/293.15 MTS-810 0.15 1 15 10 150 100 ZwickRoell-5020 1500 1000 15,000 图1 273.15 K 和应变率 0.01 s−1.准静态测试三次的应力-应变曲线二、弹道冲击试验本研究设计了针对PC/ABS复合材料的弹道冲击试验，以评估其在高速冲击下的动态响应特性。试验设计靶板规格：采用100 mm×100 mm×3 mm的矩形PC/ABS复合材料靶板，四角设置固定孔以确保试验过程中靶板稳定。弹体参数：弹体由Cr12MoV工具钢制成，采用圆柱形设计（长度24 mm，直径12 mm的半球形头部），具体尺寸见图2a。试验装置与条件冲击系统：使用口径12 mm、枪管长度4 m的高速气 枪发射弹体，弹体平均撞击速度控制在34 m/s（高于高速公路车辆平均速度30 m/s[22]）。通过测速仪实时监测弹体速度，确保试验条件的一致性。数据采集系统：采用高速摄像系统（采样频率4×10⁴ Hz，分辨率384×288像素）记录弹体侵彻过程，拍摄间隔为25 μs。摄像方向垂直于弹体入射路径，并在靶板后方设置45°反射镜，以同步捕捉靶板背面的变形与断裂形貌（图2b）。三、材料参数标定冲击载荷下的粘弹热力学耦合模型数量较少，且存在参数测定过多等问题。由于 PC/ABS 在整车制备中应用广泛，研究其失效行为应关注应变率和温度的影响，尤其是高应变率下的塑性变形与失效，而非其蠕变或松弛行为。因此，本文选取 JC 本构模型。JC 本构模型是描述材料塑性硬化、应变率效应和热软化的唯象模型。这三个唯象公式在 JC 本构模型中以乘法关系连接，JC 本构模型主要适用于大变形、高应变率、高温下的材料，适用于大多数材料的数值模拟。方程（1）的形式为：(1) 在哪里 ——等效塑性应变； ——等效塑性应变率； — 参考应变率； — 参考温度； ——材料的熔点温度； — 测试温度。 是材料参数。以下是本研究确定参数的原理 ：在参考温度为 238.15 K、应变速率为 0.01 s−1时，由力-位移曲线转化而来的应力-应变曲线得出 和 238.15 K 下不同应变速率下的应力-应变曲线的值为 0.01 s−1时不同温度下的应力-应变曲线的值为 根据该原则，JC 本构模型的初始参数确定为： 通过上述反演方法确定了 JC 本构参数，如 Table 2 所示。数值结果与试验结果对比如 Figure 8 所示，需要指出的是，由于没有破坏准则，所以数值模拟曲线没有出现陡峭的下降，此时在中、低应变速率下变形过程基本可以再现。表 2.PC/ABS 复合材料的 JC 本构参数。A B C n m 57.5兆帕120兆帕0.032 1.734 1.02 Figure 8四、JC 本构参数验证 反演方法所得结果仅适用于中、低应变率，而气囊展开时覆盖材料经常受到高速冲击载荷。为了检验反演方法确定的 JC 本构模型在高应变下的有效性，采用高应变率拉伸试验对其进行了校准。使用 ZwickRoell-5020 高速拉伸试验机，在特定温度室的三个环境温度下，以 100 s −1 和 1000 s −1 的应变率进行单轴拉伸试验。将试验结果与反演方法获得的 JC 本构模型数值模拟数据进行比较。JC 本构模型的数值模拟结果与试验结果较为接近。 Figure 9 给出了对比结果，说明反演方法得到的 JC 本构模型适用于大应变率。需要特别指出的是，实验过程可视为高速撞击下的绝热过程[ 33 ]，实验中会产生大量热量，图案中会产生很大程度的温度变化。因此，在模拟中加入了机械热效应，以校正绝热增温的影响。Figure 9一维拉伸状态下的失效行为分析上面构建的 JC 本构参数不包括失效准则。研究人员关注的是工程中材料的失效行为。本章研究了这种复合材料的失效行为。高速摄像机记录了 PC/ABS 复合材料在室温下以不同应变速率的变形和断裂过程。请注意，其他温度是通过环境温度箱实现的。因此，其他温度下的过程无法用高速摄像机捕捉到。 Figure 10 显示典型结果。数字图像相关 (DIC) 可用于使用光学传感器实时测量绘图图案中的尺寸变化。 利用 DIC 获取了不同应变速率下试件拉伸变形信息及试件破坏时刻的局部变形情况，并计算了断裂应变 计算公式（2）显示计算结果，其中 表示材料的原始横截面积， 表示材料断裂的横截面积）[ 34 ]。失效区域的局部应变列于 Table 3 中。(2) Figure 10表 3.断裂应变和样品温度。应变速率（s −1 ）0.01 1 10 100 1000 断裂应变0.825 0.833 0.922 0.646 0.527 模拟样品温度（K）293 313 328 311 309 理论样品温度（K）293 320.6 326.5 313.9 312.1 本文基于反演方法构建的 JC 本构模型，结合断裂应变根据试件的变形信息计算得到 JC 本构方程。通过数值模拟分析了不同应变率下 PC/ABS 复合材料室温（293.15 K）单轴拉伸行为，模拟将该过程与高速摄像机的记录进行比较。 Figure 10 展示了典型的对比结果， Figure 11 将不同应变率下试验的力-位移曲线与数值模拟结果进行了比较。从图中可以看出，利用反演方法确定的 JC 本构和利用 DIC 获得的失效参数可以再现室温下不同应变率下 PC/ABS 复合材料单轴拉伸试验过程。JC 失效模型如公式 (3) 所示，该模型以解耦形式描述应力三轴度、应变速率和温度的影响，以便可以删除对研究不重要的因素。 ( 3) 其中 —静水应力； ——等效强度； — 样品温度； —参考温度。 需要注意的是，由于拉伸过程中产生大量的热，样品的温度会发生很大的变化，因此 指的是试样断裂前材料内部的温度，而不是环境温度。高分子材料对温度高度敏感，在拉伸过程中，试样内会产生大量的热量，导致材料温度迅速升高。为了准确分析失效行为，有限元模拟中的失效温度选择为试件的断裂温度。本文主要研究应变速率和温度对 PC/ABS 材料失效行为的影响，其中 JC 失效模型降级为公式（4）。(4) 其中 , , 。 拟合后可以确定 如图Table 3所示，其理论温度也可通过公式（5）计算得出，其中密度ρ=1120Kg/m3，比热c=1400J/(kg.K),机械热效应 ， 是升高的温度。 在 0.01 s−1的应变速率下，样品与外界处于热平衡状态， 可视为与环境温度相同。(5) 弹道冲击失效行为分析 对弹道冲击试验进行了数值模拟。输入是 通过DIC 计算的失效断裂应变和 JC 本构参数。使用 0.25 的摩擦系数来描述子弹与靶板之间的相互作用。Figure 12 显示了 PC/ABS 复合材料侵彻实验的有限元模拟结果，以及与侵彻靶板的实验过程的对比。时间（t0= 0 μs）定义为子弹开始接触靶板的时刻。然后，我们将模拟结果与以下各个时刻的靶板损伤形状和子弹位置进行了对比：t0= 0 μs、t1= 200 μs、t2= 400 μs、t3= 600 μs、t4= 800 μs 和 t5= 1000 μs。发现有限元模拟结果与实验结果高度一致。将有限元模拟的损伤形状与实验回收的靶板进行了对比，如Figure 13所示。考虑到实验过程的复杂性、靶板固定误差以及气压不稳定、子弹射出偏差等诸多不可抗拒的因素，有限元模拟结果是可以接受的。总之，本文利用试件的变形信息计算得到的失效断裂应变和利用反演方法建立的 PC/ABS 复合材料的 JC 本构参数都是准确的。 Figure 12 模拟与实验中靶板侵彻破坏过程对比 来源：ABAQUS仿真世界