为了表征超弹性材料，需要进行各种测试获取实验数据，包括承受单轴拉伸和压缩、双轴拉伸和压缩以及扭转测试。今天，我们向大家介绍如何使用通过单轴和双轴测试获得的拉伸和压缩测试数据，模拟由弹性泡沫材料制成的球体的压缩。通过案例模型，我们演示了如何使用可压缩的Storakers超弹性材料模型进行仿真计算，以及如何通过单轴和等双轴测试数据计算力-伸长率的关系。

使用测试数据进行压缩分析下图是进行单轴和双轴拉伸和压缩测试的示意图。在这里我们不讨论进行此类测试所需的设置，而是假设已获取实验数据。在之前关于获取材料数据用于结构力学仿真的文章中（详见文章底部“阅读原文”），我们说过要正确表征材料变形，获取单轴和等双轴测试的数据至关重要。如果仅具有单轴测试数据，那么对涉及多维应力变形状态模拟的结果将不准确。下图是泡沫材料的单轴和双轴测量数据：力-伸长率

。

值大于1代表张力数据，

值小于1代表压缩数据。由单轴和双轴测试得到的力-伸长率数据。Storakers超弹性材料模型Storakers材料模型通常用于模拟高度可压缩的泡沫材料。Storakers材料模型的应变能密度函数为(1)其中

，

和

是主伸长率；

是弹性体积比；

，

和

是Storakers材料参数。在超弹性材料仿真中，需要材料参数作为输入数据。我们通过将实验测试数据与代表应力或力对伸长率的解析表达式进行拟合，来计算这些参数。在通过各种超弹性材料模型对测量数据进行拟合的文章中，我们已经详细讨论了此一般过程。与之前文章中讨论的不可压缩材料模型不同，本篇文章研究的是可压缩材料的解析表达式。下面，让我们看一下Storakers材料模型的单轴和双轴测试中力与伸长率的关系。单轴测试参考上面单轴测试的示意图，假设用于材料表征的泡沫材料在主方向2上承受单轴载荷。假设仅允许弹性变形，那么各向同性超弹性材料的单轴变形的主伸长率由下式给出：

对于单轴测试，主方向1和3上的主柯西应力由下式给出

。对于可压缩的各向同性超弹性材料，主柯西应力与主伸长率之间的关系由下式给出：

(2)上述关系来自热力学第二定律和应力在主伸长率的频谱表示的影响之一。根据热力学第二定律，我们可知

，其中，

是第二类皮奥拉-基尔霍夫应力张量，

是弹性柯西-格林张量。频谱形式为，其中，

代表主要参考方向。已知柯西应力

，以频谱表示的变形梯度张量

可表示为

代表主要的空间方向），我们可以得到，其中

是主柯西应力（点击此处，了解有关不同应力测量及其方向的更多信息）。将等式1中的

代入

，我们得到然而，众所周知

，其中

和

是变形梯度张量。这种关系来自张量演算。当按照主伸长率重写公式时，我们可以得到

。将此关系代入上式并使

等于0，得到(3)

和

对于每个

，上面等式均适用。单轴力由下式计算(4)

其中，

是沿主方向的变形尺寸，

是沿主方向的原始尺寸。使用等式2，3，和4，我们可以将单轴力表示为

(5)等双轴测试假设泡沫材料在主方向1和2上受到等双轴力，并且仅允许弹性变形。对于各向同性超弹性材料的等双轴变形，主伸长率为

对于等双轴测试，主柯西应力

。由等式2，我们得到(6)

对于每个

，上面等式均适用。

等双轴力

由下式给出(7)

使用等式2，6，和7，我们得到

(8)使用优化接口计算Storakers材料参数根据等式5和8在N=2时得出的解析表达式，通过COMSOLMultiphysics?软件中的优化接口，对力-伸长率测试数据进行曲线拟合。下图显示了测得数据与拟合数据的对比，单轴和等双轴最小二乘拟合的权重相等。使用Storakers材料模型拟合材料参数。根据以上拟合，我们计算出用于测试的泡沫材料的Storakers材料参数：模拟球形泡沫材料的压缩上文中，我们已经计算出了泡沫材料的Storakers材料参数值，因此可以直接在球形泡沫压缩模拟中使用这些值。下图显示了模拟的设置，其中有一个中空的球形泡沫，该球形泡沫被一个较小的球形压头压在圆柱形刚性板上。右图为描述三维问题的二维轴对称模型。

泡沫材料压缩的三维示意图（左）和COMSOLMultiphysics中的二维轴对称设置（右）。

下图显示了用于模拟超弹性泡沫材料模型的设置。压头和刚性板被建模为线弹性材料。压头、泡沫和刚性板处于接触状态，因此我们在它们的边界之间定义了一个接触对。使用指定位移条件，给压头施加11毫米的向下移动，使得泡沫内部边界与其自身接触，称为自接触。通过定义第二个接触对建立自接触模型，我们可以将泡沫的内部边界同时定义为源和目标。假设我们不知道内表面的哪一部分将与相同表面的哪一部分接触，将底部刚性板设置为完全固定。超弹性材料模型的设置。结果与分析下面的三维图显示了旋转几何体部分中的vonMises应力图。（3D演示动画可点击阅读原文查看）通过压头压缩的球形泡沫中的VonMises应力。色标的最大值上限为10kPa。下图显示了在模拟的最终时间点时，泡沫边界上的接触压力。模拟的最终时间点时的接触压力图。当然，我们不仅限于模拟球形泡沫材料。现在，我们已经知道了给定伸长率范围内的Storakers材料参数，可以在任何其他情况下对相同的泡沫进行建模，只要其变形与用于计算参数的测试数据一致即可。但是，我们还应该记住，尽管通过曲线拟合实验测试数据获得了材料参数，但是如果材料模型违反了Drucker的稳定性标准，则材料模型的稳定性仍然会出现问题。

推荐阅读

COMSOL免费开放15个专题培训视频资料那些高抗疲劳和抗腐蚀的复合材料，你知道如何仿真吗？多孔材料的数字化表征特殊多孔弹性超材料的仿真分析材料科学的重点议题：永不忘本的形状记忆合金

本文内容来自COMSOL博客，点击“阅读原文”，阅读更多延伸文章。

预览时标签不可点收录于合集#个上一篇下一篇
转载请注明:http://www.aideyishus.com/lktp/208.html