首页/文章/ 详情

Abaqus 中的超弹性材料 – 它们是什么?它们如何工作?

7小时前浏览2

Abaqus 允许我们使用多种不同类型的稳健材料模型进行建模,具体取决于实际的物理特性。在本文中,我们将重点讨论超弹性材料,这种材料通常用于橡胶聚合物、生物材料和固体推进剂。本文旨在概述超弹性材料,并介绍如何对其进行建模和模型校准。

超弹性简介

自然界中存在各种类型的材料,它们在承受载荷时会表现出不同的机械行为。例如,钢、铝和铜等延展性金属在屈服前表现出线性弹性行为,然后在断裂前经历相当大的塑性(不可恢复)变形。陶瓷、玻璃和石墨等脆性材料在断裂前经历线性弹性变形,几乎没有塑性变形。并非所有材料都表现出线性弹性。事实上,有一类材料在施加载荷后,在应变超过 100% 时会发生巨大的弹性变形,并且在失效时没有塑性变形。这些材料被归类为超弹性材料。超弹性材料的应力-应变响应是高度非线性的,并且单调增加直至失效,如下图所示。



Hyperelastic Material In Abaqus



超弹性材料模型用于捕捉橡胶、填充弹性体、生物组织、固体推进剂等材料的弹性大应变行为。刚度随着样品的变形而变化,并且这些材料在弹性状态下储存的内部能量在卸载时完全恢复。



Hyperelastic Material In Abaqus



超弹性材料的压缩性

大多数超弹性材料(例如橡胶和生物组织)几乎或完全不可压缩。但其他材料(例如海绵和泡沫)是可压缩的,需要一类特殊的超弹性材料模型(例如 Ogden 泡沫)来捕捉其变形行为。不可压缩性项根据材料的初始体积模量 () 定义如下:


Hyperelastic Material In Abaqus


当 d 值减小时,材料趋于不可压缩。它没有范围,通常为正值。对于非约束材料,假设其具有不可压缩行为通常就足够了,但当材料受到约束时,体积行为起着至关重要的作用。

应变能势

当试样承受恒定应变时,材料中储存的应变能(等于应力-应变曲线下的面积)保持不变。因此,超弹性材料是以应变能势为模型的。



Hyperelastic Material In Abaqus



能量势 (Ψ) 定义为材料中应变张量 (ε) 的函数,如下所示:

Ψ = f (e)

应变能势可分为偏应变能和体积应变:

Ψ = Ψdev + Ψvol

材料的偏量行为取决于变形模式,即单轴拉伸、单轴压缩和剪切。Abaqus 中提供了多种应变能势来模拟超弹性材料,这些应变能势可分为物理模型和现象学模型。物理模型(例如 Arruda-Boyce 和 Van der Waals)考虑材料的微观结构来评估其响应。简单的现象学模型(例如 Neo-Hookean 和 Yeoh)效果良好,可以基于连续介质力学评估材料响应。下表显示了不同超弹性模型及其材料参数的应变能势。


Hyperelastic Material In Abaqus


在上述应变能势中,材料常数为“D”的分量表示体积项,Jel 表示弹性体积比。偏差分量可使用主伸长率(λ1、λ2、λ3)或应变不变量(ɪ1、ɪ2、ɪ3)计算。主伸长率是材料在主方向上的当前长度与原始长度之比。在 Ogden 模型和多项式模型中,N 表示函数的阶数。模型的非线性和材料常数的数量随函数阶数的增加而增加。

实验超弹性数据

需要实验数据来根据应变能势计算材料常数。为了与实验数据获得良好的相关性,我们需要考虑所有三种变形模式——拉伸、压缩和剪切。对于遵循胡克定律的金属,我们可以仅使用拉伸试验数据来准确预测材料的偏变形行为。可以根据金属的拉伸数据计算出独特的压缩和剪切响应。但是,对于橡胶和弹性体等软材料,它们在三种变形模式下都表现出非线性响应。因此,这三种模式下的刚度各不相同,仅使用拉伸试验数据无法计算出唯一的解。因此,我们需要均匀的单轴拉伸、压缩和简单剪切试验数据来捕捉超弹性材料的偏变形行为。此外,我们还需要体积试验数据来确定材料的体积响应。在不可压缩材料假设下,可以忽略体积响应。

应变能函数关于应变不变量的导数用于建立不可压缩材料的应力-应变关系。以主拉伸率 (λ1 , λ2 , λ3) 表示的变形梯度表示为:


Hyperelastic Material In Abaqus


主拉伸量(λi)与主公称应变(ϵi)的关系如下:


Hyperelastic Material In Abaqus


下图显示了 Abaqus 中用于校准超弹性材料模型的可接受实验测试数据以及每个测试中的主拉伸值。通过这些主拉伸值,可以计算出偏应变不变量和名义应变。



Hyperelastic Material In Abaqus



利用实验数据校准 Abaqus 材料常数

选择应变能势函数后,我们需要通过求解优化问题来计算材料系数。幸运的是,Abaqus 提供了一个非常好的工具,它直接内置于 Abaqus CAE 中。优化是通过最小化实验数据 (sE) 和模型预测数据 (sM) 之间的相对累积误差来实现的。这种方法最小化的是相对误差而不是绝对误差,以便在较低应变下获得更好的拟合效果。



Hyperelastic Material In Abaqus



相对最小二乘误差函数如下:


Hyperelastic Material In Abaqus


根据应变能势函数的性质和材料常数的数量,可以采用线性或非线性回归方法进行最小化。例如,对于应变能势与材料常数呈线性关系的多项式模型,可以使用线性回归。对于非线性模型,例如 Ogden 和 Van der Waals,则需要使用非线性回归来校准材料常数。误差函数的最小化值称为残差,用于确定不同材料模型之间的拟合优度。

最后的想法

Abaqus 提供了多种材料模型来捕捉超弹性材料(例如橡胶、弹性体和软组织)的非线性弹性行为。希望本文能帮助您对不同类型的应变能势、校准材料系数所需的测试数据以及校准过程提供一些基本的了解。

来源:ABAQUS仿真世界
Abaqus断裂非线性材料试验
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-09-10
最近编辑:7小时前
yunduan082
硕士 | 仿真主任工程... Abaqus仿真世界
获赞 178粉丝 321文章 403课程 0
点赞
收藏
作者推荐

HDPE热力学与动态力学性能的本构模型

https://doi.org/10.1016/j.polymertesting.2025.108935核心发现•标定并验证了三种本构模型,用于预测高密度聚乙烯(HDPE)材料随温度与应变率变化的力学行为•引入失效模型以模拟板材损伤并验证响应预测•证实三网络(TN)模型在预测HDPE变形响应方面具有更高精度•基于TN模型研究了不同标准尺寸比(SDR)HDPE管道的抗冲击性能摘要准确认知高密度聚乙烯(HDPE)的温度与应变率相关力学行为,对保障其工程应用(尤其是海底管道抗意外载荷设计)的长期完整性至关重要。鉴于HDPE对时间与温度的高度敏感性,需建立强健且复杂的本构与损伤模型。本研究对标定的三种HDPE本构模型(三网模型TN、三网粘塑性模型TNV及Johnson-Cook模型JC)进行了系统校准,校准基于多应变率与高温条件下的拉伸试验,并通过多速度落锤冲击试验验证。数值模拟结果表明:相较于JC与TNV模型,TN模型在预测HDPE力学响应与损伤行为方面更具优势,其冲击试验数据与数值预测在所有工况下均呈现显著一致性。该发现证实TN材料模型在HDPE构件低速冲击响应设计与评估中具有优异可靠性。研究还验证了校准模型在预测海底HDPE管道结构冲击载荷下损伤行为及穿孔极限方面的有效性。本构标定技术路线3.1Johnson-Cook本构模型Johnson-Cook(JC)模型[47]因其参数少、形式简洁的特点,被广泛应用于金属[48,49]和聚合物[43,50]等多种材料的力学表征。该模型以乘积形式综合考虑了应变硬化、应变率硬化及温度软化效应,其流动应力方程如式(1)所示[51]:式中:σ为等效流动应力ε_p为塑性应变ε̇_p与ε̇_0分别为当前应变率与参考应变率T、T_0、T_m分别表示瞬时温度、参考温度与熔点温度A为初始屈服强度B为应变硬化系数n为硬化指数C为应变率系数m为温度软化指数本研究中参考应变率ε̇_0与参考温度T_0分别取0.00333s^-1和23°C。JC失效模型与线性损伤累积模型分别由式(2)-(3)表示[51]:式中:σ_vm、σ_m、ε_f^p分别为vonMises等效应力、平均应力与塑性失效应变D_1~D_5分别代表初始断裂应变、指数因子、应力三轴度因子、应变率因子和温度因子,共同表征材料损伤状态需特别注意式(2)中D_3为三轴度乘子而非应力三轴度本身ω为损伤应力状态,D为累积损伤变量(当D=1时材料失效)[52]考虑到材料强度随损伤累积而退化,应力与损伤的线性演化关系由式(4)描述[52]:为模拟HDPE材料在落锤冲击试验中的动态响应,本研究对JC塑性模型与JC损伤模型进行了联合标定。通过图2所示的23°C、60°C和95°C真应力-应变曲线预测,获得弹性参数与模型参数(见表2)。基于数字图像相关(DIC)技术测得的轴向/横向应变场(图1),确定了PE100级HDPE在室温参考应变率下的弹性模量为1.2GPa,泊松比为0.419,该结果与文献[53,54]报道的HDPE参数范围一致。图1(b)中轴向应变场在试样下端呈现峰值,这是由于该区域相较于上端承受更大的纵向拉伸;而横向应变场则表现为垂直于加载方向的横向收缩,由此形成图1(b)所示的应变云图特征。Table2.HDPEJC本构标定参数EmptyCellVariableUnitValueElasticPropertiesEMPa1200–0.419PlasticityParametersAMPa8.86BMPa30.56n–0.384C–0.046m–0.762DamageConstantsD1–0.098D2–26.88D3–8.81D4–0.005D5–0.01Disp.DamageEvolutionm0.00133.2三网模型TN模型通过三个相互关联的分子网络(如图4示意结构所示)以并联形式描述HDPE的非线性粘塑性行为。其本构关系具有以下特征:应力叠加原理总柯西应力(σ_total)为各网络应力分量的叠加和:σ_total=σ_A+σ_B+σ_C变形协调条件所有网络承受相同的总变形梯度(F_total),即:F_total=F_A=F_B=F_C3.3三网粘塑性模型TNV模型在TN模型基础上进行了扩展,通过引入损伤演化机制和压力依赖特性增强了模型描述能力(如图6所示)。该模型具有以下核心特征:网络架构特性保留三网并联结构,但各分子支链的流动单元数量与TN模型存在差异采用超弹性本构方程描述网络力学行为支持根据应用需求灵活配置弹簧/阻尼单元(可选择性包含或排除特定支链元件)多物理场耦合可集成应力/应变失效模型,实现热-力-损伤耦合分析总变形梯度(F)采用乘法分解:F=Fe·Fvp,其中Fe为弹性分量,Fvp为粘塑性分量Table4.TNVmodelcalibratedandfailureparametersat23°C.HDPE板材实际变形形貌与TN模型有限元位移(Uy)云图对比:(a)2m/s冲击下的压痕损伤;(b)3m/s及(c)4m/s冲击下的贯穿损伤来源:ABAQUS仿真世界

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈