中国学者一作兼通讯！复合材料粘弹性仿真，登上顶刊！

来源公众号：材料科研SCI

文献解读

单向复合材料层合板和可展开结构粘弹性行为的多尺度建模

后台私信回复“20230505”，即可获取PDF全文！

研究背景

航天器通常由先进的复合材料层合板和可展开结构制作而成。由于聚合物固有的粘弹性，纤维增强聚合物复合材料表现出或多或少的模量松弛行为。了解它们的粘弹性行为对于实现复合材料的优化设计至关重要，特别是对于航空航天工业中存储应变能可展开结构的设计。

内容简介

针对上述情况，四川大学安宁教授团队提出了一种基于有限元(FE)的双尺度计算策略，能够评估单向复合材料层合板和一般壳结构的粘弹性行为。在Abaqus/CAE中开发了RVE分析插件工具，通过将纤维和基体的已知性质以及微观结构几何形状作为输入，快速估计单向复合材料的有效正交各向异性粘弹性。解决了一些基准问题，验证了所提插件工具的准确性和效率。

成果以题为《Multiscale modeling of viscoelastic behavior of unidirectional composite laminates and deployable structures》发表于《Materials & Design》上。

点击末尾“阅读原文”即可链接至文章页面，文章来源如下：

https://doi.org/10.1016/j.matdes.2022.110754

结果讨论

1、弹性RVE有效性能的计算

利用粘弹性RVE计算器将纤维体积分数Vf设置为0.1 ~ 0.6，并在此范围内变化，生成了正方形、菱形和六边形阵列三种不同纤维排列方式的弹性RVE，计算了其有效材料性能，并与图1中随机纤维分布的弹性有限元RVE进行了对比。

结果表明，在每个子图中，不同纤维排列方式下的三种模型计算的预测结果一致，但在特定情况下仍然存在一些偏差，这反映了纤维排列的影响。

2、粘弹性RVE有效松弛特性的计算

为了进一步验证粘弹性复合材料有效松弛模量预测工具的功能，作者使用粘弹性RVE计算器插件计算了粘弹性RVE的正交异性应力松弛行为，然后将有效粘弹性特性实现到用户自定义的材料子程序UMAT中进行均匀化。

粘弹性RVE计算器插件和UMAT预测的松弛计算结果如图2所示。

最近，Rodríguez-Ramos等人利用渐近均匀化方法(AHM)建立了一种计算纤维粘弹性复合材料松弛有效模量的封闭理论，并将AHM预测的理论结果绘制在这里进行比较。可以看出，所有的结果都非常吻合。

最后，可以得出结论，由UMAT定义的单材料均匀模型可以从插件中由双材料定义的RVE获得的松弛数据中获得等效的均匀松弛响应。

3、纤维排列和环境温度的影响

在解决了基准问题并验证了所提出的分析工具的准确性之后，作者继续应用该工具来分析RVE微观结构几何形状(即纤维体积分数和纤维排列)以及环境温度对粘弹性复合材料平均正交异性弛豫性能的影响。

图3给出了在环境温度恒定且等于参考温度40℃时，不同纤维排列的RVEs在Vf=0.2、0.4、0.64三个体积分数下的松弛结果。

结果表明，具有六边形纤维阵列的rve的松弛曲线位于三组数据的中间，并且在菱形阵列中可以观察到较高的值，而在正方形阵列中可以观察到较低的值，并且发现偏差随着纤维体积分数的增加而增加。

作者在图4中报道了相同纤维体积分数为0.64的RVEs在30℃、40℃、50℃等不同环境温度下的松弛结果。

从图中可以看出，实验温度越高，复合材料的松弛过程越快，而实验温度越低，松弛过程越慢。

4、复合胶带弹簧铰链的应用

在本节中，作者重点研究可展开复合胶带弹簧铰链，并使用提出的分析工具揭示粘弹性松弛对其展开动力学的影响。图5显示了在这项工作中考虑的复合磁带弹簧铰链的几何形状。

图6为考虑边界条件和载荷的复合带簧铰链有限元网格图。复合材料的密度设置为1580.76 kg/m3，在整个分析过程中，整个结构在负Z方向上定义重力载荷；粘度系数通过反复试验过程确定为1.5*10-8。

图7(a)给出了展开角随时间的演变，其中数值结果与Fernandes等人报道的实验数据(方形标记)吻合较好。需要指出的是，在模拟中随机出现了两组不同的展开角-时间响应曲线。这是由于折叠过程中形成的铰链结构的可变性，如图7(b)所示。

图8给出了两种配置下的模拟部署快照，与参考文献中的实验观察结果很好地吻合。

图9为复合材料带弹簧铰链在50°C下不同松弛期的展开动力学变化。

仿真结果表明，无论铰链以完全折叠构型还是半折叠构型α或β存放，装载后铰链的展开速度都比初始不放松情况慢几分之一秒左右，并且存储时间越长，展开旋转速度越慢，这与Fernandes等人的研究结果一致。

结论

在微观尺度上，本文开发了一种数值分析工具，并通过一个名为粘弹性RVE计算器的图形用户界面插件与Abaqus/CAE相连接，该工具可以通过将组成材料的微观结构几何形状和已知性质作为输入，计算单向复合材料的正交异性粘弹性。

1、六角形纤维阵列复合材料的松弛曲线始终处于正方形和菱形阵列的松弛曲线之间；较高的实验温度会加速松弛过程，较低的实验温度会延迟松弛过程。

2、插件输出的随时间变化的层状特性随后可以使用Abaqus内置函数在宏观尺度上对复合材料层状板的粘弹性行为进行建模，以定义堆叠顺序并相应地更新材料特性。

3、建立了一个三维壳体有限元模型，并量化了长期装载对复合材料带弹簧铰链展开动力学的影响。结果表明，与以前的实验结果相一致，与未储存的情况相比，长期储存将导致铰链的展开率显著降低。

招人写稿啦！！！

“材料科研SCI"微信公众号自创办以来得到了广大科研工作者和研究生的广泛关注和支持。为更好地服务材料学科的研究与应用，本公众号现需长期供稿作者！

您可发送简历到小编邮箱:

likeapoem@163.com

邮件主题注明“姓名+供稿作者”。

非常欢迎材料相关专业的科研爱好者加入“材料科研SCI"团队！