一种针对Mullins效应参数的球压痕表征方法与流程

本发明属于压痕实验技术领域，具体涉及一种高分子聚合物材料的mullins效应表征方法。

背景技术：

压痕是一种深度相关的测试技术，源于18世纪中叶由矿物学家发展的硬度划痕测试，随后被用来测量金属的硬度。压痕实验利用压头压入被测材料，通过分析压痕载荷–压入深度曲线获得相关力学性质，被广泛应用于宏观及微纳观尺度下材料的线弹性、塑性及粘弹性性质力学表征。

1881年，hertz发展了接触理论，hertz解也被广泛应用于球压头压痕实验中的弹性接触问题。lee和radok将hertz接触理论扩展至粘弹性材料，通过蠕变函数预测粘弹性材料对外加载荷的响应。johnson等将hertz解扩展至粘附接触的情况。高玉臣等发展了新的应变能密度函数，研究了橡胶类材料与刚性锥之间的大变形接触问题。

在连续介质力学框架下，高分子聚合物材料可以用线弹性、超弹性、线性粘弹性、非线性粘弹性、多孔弹性或其它双相、多相本构关系来描述，需要根据材料表现出的力学响应选择合适的本构模型。mullins效应，即在典型橡胶材料的应力-应变循环试验中，当橡胶经历“加载-卸载-再加载循环”时，卸载应力和重加载应力远远小于初始加载时的应力。在重加载时，随着应变的增加，应力-应变曲线首先沿着卸载路径变化，然后随着应变的进一步增加，应力-应变曲线与主曲线重合。mullins效应是高分子聚合物材料的一个显著特性，对减振器等承载元件的设计效果影响很大，值得重点关注。

mullins和tobin分别于1953年和1957年对填充硫化橡胶的应力软化现象进行了试验和模型预测对照。shen于2001年对聚氨酯泡沫体受压缩时的mullins效应所做的试验，获得了与ogden模型相吻合的应力应变曲线。cheng和chen于2003年对三元乙丙橡胶(epdm)的mullins效应做了预测，并与模型曲线进行了对比。johnson和beatty对丁纳橡胶和氯丁橡胶的mullins效应做了拉伸试验。lu、以及hao分别对粘弹性和mullins效应进行了本构模型参数表征。tang等研究了纳米复合材料水凝胶的mullins效应，考虑了残余变形和非理想循环加载的影响。zhang等用mullins效应和内聚力模型对含裂纹水凝胶等弹性体进行了断裂分析。此外，还有众多研究均对mullins效应进行了表征。上述表征多为拉伸、压缩等常规方法，在试验中需要制备特定尺寸和形状的样件，且在样件机加工过程中很可能对材料表观应力软化效应产生影响。相比而言，压痕方法的测量更为简便，更易实现原位、多尺度和局部表征测量，但是目前利用压痕方法进行mullins效应表征的工作较为少见。mullins效应虽早已有本构模型嵌入至abaqus，但目前鲜有专门利用压痕方法表征其本构参数的报道。

技术实现要素：

有鉴于此，本发明提供了一种针对mullins效应参数的球压痕表征方法，能够实现原位、多尺度和局部表征高分子聚合物材料的mullins效应特性。

一种针对mullins效应参数的球压痕表征方法，该方法实现的步骤如下：

步骤一：通过在压痕实验中引入无量纲函数来建立压痕实验中各物理量之间的关系；

步骤二：根据含mullins效应的有限元模型计算获得上述无量纲函数的显示表达；

步骤三：进行球压头压痕实验获得实际材料的加卸载载荷-位移曲线；

步骤四：根据无量纲函数的显示表达以及球压头压痕实验加载载荷-位移曲线反演得到材料的初始剪切模量和超弹性力学参数；

步骤五：根据无量纲函数的显示表达以及球压头压痕实验卸载载荷-位移曲线反演得到材料的mullins效应参数。

进一步地，所述步骤一的内容包括：将π定理应用于neo-hookean、mooney–rivlin、fung以及arruda–boyce四种超弹性本构模型并分别对应ogden-roxburghmullins效应模型描述的球压头压痕加载实验，得到球压头压痕实验中物理量之间的关系：

其中，p为球压痕加载载荷，μ0为初始剪切模量，h为压入深度，r为球压头半径；πnh、πmr、πf和πab为无量纲函数，α、β和λm分别为对应本构模型中的超弹性力学参数，r，m，b为对应mullins效应本构模型中的参数；

将π定理应用于neo-hookean、mooney–rivlin、fung以及arruda–boyce四种超弹性本构模型并分别对应ogden-roxburghmullins效应模型描述的球压头压痕卸载实验，得到球压头压痕实验中物理量之间的关系：

其中，p^j为球压痕卸载载荷，μ0为初始剪切模量，h^j为压入深度，r为球压头半径；π^jnh、π^jmr、π^jf和π^jab为无量纲函数，α、β和λm分别为对应本构模型中的超弹性力学参数，r，m，b为对应mullins效应本构模型中的参数，下标j表示第j次卸载，每j次加载最大深度依次递增。

进一步地，所述步骤二的内容包括：利用商用软件开展大规模非线性有限元模拟，得到无量纲函数πnh、πmr、πf和πab的的曲线，从而确定无量纲函数的显示表达：

πmr(α，r，m，b)＝c1+c2α+c3α²；

πf(β，r，m，b)＝d1+d2β+d3β²；

其中，对于neo-hookean本构模型，无量纲函数πnh的系数κ1和κ2为常数；对于mooney–rivlin和fung本构模型，无量纲函数πmr和πf的系数ci和di(i＝1,2,3)依赖于h/r；对于arruda–boyce本构模型，无量纲函数πab的系数a0、ai和bi(i＝1,2,3)依赖于h/r；

同样地，利用商用软件开展大规模非线性有限元模拟，得到无量纲函数π^jnh、π^jmr、π^jf和π^jab的曲线，从而确定无量纲函数的显示表达：

其中，无量纲函数π^jnh、π^jmr、π^jf和π^jab的系数ei、fi、gi和hi(i＝1,2,3)依赖于r，m，b，超弹性力学参数α、β和λm以及卸载次数j。

进一步地，所述步骤三的内容包括：选取j个最大压入深度开展球压头压痕实验，记录不同压入深度的加载和卸载载荷-位移曲线。

进一步地，所述步骤四的内容包括：综合利用所述无量纲函数πnh、πmr、πf和πab的显示表达，可知当时，所有本构模型中球压痕加载载荷p和压入深度h呈简单函数关系，而且基本与mullins效应本构模型参数r，m，b无关：

根据球压头压痕实验加载载荷-位移曲线可以反演得到材料的初始剪切模量μ0；

在确定初始剪切模量之后，利用时的球压头压痕实验加载载荷-位移曲线确定不同本构模型的超弹性力学参数：

进一步地，所述步骤五的内容包括：根据无量纲函数π^jnh、π^jmr、π^jf和π^jab的显示表达以及球压头压痕实验中的卸载载荷p^j和压入深度h^j等信息反演得到材料的mullins效应参数r，m，b：

有益效果：

1、本发明的mullins效应参数的球压痕表征方法通过在球压头压痕实验中引入无量纲函数来建立压痕实验中各物理量之间的关系，然后确定该无量纲函数的显示表达，最后根据该无量纲函数的显示表达反演得到高分子聚合物材料的初始剪切模量、超弹性力学参数以及mullins效应参数。

2、本发明的mullins效应参数的球压痕表征方法是基于压痕实验，因此，可以直接作用于高分子聚合物材料实现原位表征；同时，可以根据高分子聚合物材料的不同尺寸调整球压头尺寸的大小，从而实现多尺度和局部表征高分子聚合物材料的力学性质；创新地建立了更完备的含mullins效应高分子聚合物材料的力学性质表征方法。

附图说明

图1为本发明表征方法的步骤流程图；

图2为球压头压痕实验的示意图；

图3为含mullins效应球压痕实验的加卸载载荷-位移曲线示意图；

图4为本发明中球压头压痕实验有限元模拟及网格划分的示意图；

图5为本构模型mooney–rivlin无量纲函数的曲线图；

图6为本构模型arruda–boyce无量纲函数的曲线图；

图7为本构模型fung无量纲函数的曲线图。

具体实施方式

下面结合附图并举实施例，对本发明进行详细描述。

如图1所示，本发明提供的mullins效应参数的球压痕表征方法，其实现的步骤如下：

步骤1：通过在压痕实验中引入无量纲函数来建立压痕实验中各物理量之间的关系；

步骤2：根据含mullins效应的有限元模型计算获得上述无量纲函数的显示表达；

步骤3：开展球压头压痕实验获得实际材料的加卸载载荷-位移曲线；

步骤4：根据无量纲函数的显示表达以及球压头压痕实验加载载荷-位移曲线反演得到材料的初始剪切模量和超弹性力学参数。

步骤5：根据无量纲函数的显示表达以及球压头压痕实验卸载载荷-位移曲线反演得到材料的mullins效应参数。

本发明提供的mullins效应参数的球压痕表征方法，通过在球压头压痕实验中引入无量纲函数来建立压痕实验中各物理量之间的关系，然后确定该无量纲函数的显示表达，最后根据该无量纲函数的显示表达反演得到高分子聚合物材料的初始剪切模量、超弹性力学参数以及mullins效应参数。

可知，本发明提供的mullins效应参数的球压痕表征方法，由于其基于压痕实验，因此，可以直接作用于高分子聚合物材料实现原位表征；同时，可以根据高分子聚合物材料的不同尺寸调整球压头尺寸的大小，从而实现多尺度和局部表征高分子聚合物材料的力学性质；创新地建立了更完备的含mullins效应高分子聚合物材料的力学性质表征方法。

在本发明的一种优选的实施例中，上述步骤1包括：

所述步骤1包括：将π定理应用于neo-hookean、mooney–rivlin、fung以及arruda–boyce四种超弹性本构模型以及ogden-roxburghmullins效应模型描述的球压头压痕加载实验，得到球压头压痕实验中物理量之间的关系。

其中，p为球压痕加载载荷，μ0为初始剪切模量，h为压入深度，r为球压头半径；πnh、πmr、πf和πab为无量纲函数，α、β和λm分别为对应本构模型中的超弹性力学参数，r，m，b为对应mullins效应本构模型中的参数。

将π定理应用于neo-hookean、mooney–rivlin、fung以及arruda–boyce四种超弹性本构模型以及ogden-roxburghmullins效应模型描述的球压头压痕卸载实验，得到球压头压痕实验中物理量之间的关系。

其中，p^j为球压痕卸载载荷，μ0为初始剪切模量，h^j为压入深度，r为球压头半径；π^jnh、π^jmr、π^jf和π^jab为无量纲函数，α、β和λm分别为对应本构模型中的超弹性力学参数，r，m，b为对应mullins效应本构模型中的参数。下标j表示第j次卸载，每j次加载最大深度依次递增。

上述物理量之间的关系的推理过程如下：

对于球压头压痕实验，当球压头半径r远小于待测材料尺寸时，可看作球压头作用于半无限大空间。这一物理过程中只包含少量的几何和物理参数，结合量纲分析和大规模有限元计算可得描述球压头压痕实验响应和材料参数间关系的显式表达式。在这里，我们针对mooney–rivlin、neo-hookean、fung以及arruda–boyce超弹性本构模型，以及mullins效应模型进行研究。

对于不可压缩neo-hookean材料的球压头压痕加载实验，加载压痕载荷pnh与下列独立参数相关

pnh＝f(μ0,h,r，r，m，b)；

式中μ0是材料在基态下的剪切模量，h为压入深度，r为球压头半径。应用π定理可得：

对于不可压缩neo-hookean材料的球压头压痕卸载实验，卸载压痕载荷pnh与下列独立参数相关：

p^jnh＝g(μ0,h,r，r，m，b)；

式中h^j为压入深度。应用π定理可得：

类似地，将π定理应用于mooney–rivlin、fung以及arruda–boyce三种不可压缩超弹性材料的压痕加载及卸载实验，可以得到：

在本发明的一种优选的实施例中，上述步骤2包括：

利用商用软件开展大规模非线性有限元模拟，得到无量纲函数πnh、πmr、πf和πab的的曲线，从而确定无量纲函数的显示表达：

πmr(α，r，m，b)＝c1+c2α+c3α²；

πf(β，r，m，b)＝d1+d2β+d3β²；

其中，对于neo-hookean本构模型，无量纲函数πnh的系数κ1和κ2为常数；对于mooney–rivlin和fung本构模型，无量纲函数πmr和πf的系数ci和di(i＝1,2,3)依赖于h/r；对于arruda–boyce本构模型，无量纲函数πab的系数a0、ai和bi(i＝1,2,3)依赖于h/r。

同样地，利用商用软件开展大规模非线性有限元模拟，得到无量纲函数π^jnh、π^jmr、π^jf和π^jab的曲线，从而确定无量纲函数的显示表达：

其中，无量纲函数π^jnh、π^jmr、π^jf和π^jab的系数ei、fi、gi和hi(i＝1,2,3)依赖于r，m，b，超弹性力学参数α、β和λm以及卸载次数j。

具体地，可以利用商用有限元软件abaqus开展大规模非线性有限元模拟，以确定无量纲函数的显示表达。软件中已经有neo-hookean、mooney–rivlin和arruda–boyce三种超弹性本构模型，可以直接使用。利用用户材料子程序uhyper将fung本构模型引入abaqus软件。

mooney–rivlin本构模型中的参数α选取了可取的全部范围，即[0，1]。当α＝1时，neo-hookean本构模型是mooney–rivlin本构模型的特殊情况。arruda–boyce和fung本构模型中的参数选取了比较大的范围，以包括更广泛的生物软材料。确定无量纲函数时初始剪切模量可任意取值，这里μ0＝1mpa。计算中，变化范围为[0，1]。

图4为本发明中球压头压痕实验有限元模拟及网格划分的示意图。

如图4所示，在球压头压痕实验有限元模拟中，建立了轴对称模型，使用了41121个四节点双线性轴对称减缩积分单元(cax4rh)对基底进行离散。在压头与基底接触处适当加密网格以保证计算收敛。基底两侧边界节点无约束，下边界节点固定。在计算中，假设压头为刚性，球压头半径为2.5mm。以下结果中摩擦因数为0，分析验证了摩擦因数不为0时的影响，通过加密网格(单元数为71319个)验证了计算的收敛性。

请参考图5、图6和图7，为本构模型mooney–rivlin、fung、arruda–boyce的加载曲线对应的无量纲函数πmr、πf和πab。类似地，本构模型neo-hookean、mooney–rivlin、fung、arruda–boyce的卸载曲线对应的无量纲函数π^jnh、π^jmr、π^jf和π^jab也有类似形式，不再展示。

步骤3包括：

选取j个最大压入深度，开展球压头压痕实验，记录不同压入深度的加载和卸载载荷-位移曲线。

步骤4包括：

综合利用所述无量纲函数πnh、πmr、πf和πab的显示表达，可知当时，所有本构模型中球压痕加载载荷p和压入深度h呈简单函数关系，而且基本与mullins效应本构模型参数r，m，b无关：

根据球压头压痕实验加载载荷-位移曲线可以反演得到材料的初始剪切模量μ0。

在确定初始剪切模量之后，利用时的球压头压痕实验加载载荷-位移曲线确定不同本构模型的超弹性力学参数α、β和λm：

步骤5包括：

利用步骤4所得的材料初始剪切模量μ0、超弹性力学参数α、β和λm以及无量纲函数π^jnh、π^jmr、π^jf和π^jab的显示表达、球压头压痕实验中的卸载载荷p^j和压入深度h^j等信息反演得到材料的mullins效应参数组合r，m，b：

值得指出的是，每次实验卸载曲线均可用于mullins效应参数组合r，m，b的反演求解，可通过最小二乘法或者其他优化算法获得可与实验卸载曲线吻合的参数组合。

本发明提供的mullins效应参数的球压痕表征方法，确定了在压入深度与球压头半径比小于0.3时，压入深度与载荷间呈简单函数关系，该关系可用于求解超弹性材料的初始剪切模量。针对mooney–rivlin、fung及arruda–boyce本构模型，建立了利用球压头压痕实验表征超弹性参数的方法。结果表明，大压入高度高度和材料的锁定伸长会影响压痕实验表征其他超弹性材料参数的效果。利用球压痕实验卸载曲线，建立了反演求解mullins效应参数组合r，m，b的方法。

综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。