本文由 Mohammad Shojaeifard、Roshanak Dolatabadi、Sara Sheikhi 和 Mostafa Baghani 联合完成,研究机构主要包括 University of Tehran 机械工程学院以及 Tehran University of Medical Sciences 药学院。文章发表于 Journal of Intelligent Material Systems and Structures,并标注为 2020 年的研究成果。该研究探讨了温度响应性水凝胶材料在受拉伸-扭转复合热-机械加载条件下的耦合作用。其 DOI 为 10.1177/1045389x20951273。
研究背景定位于智能材料及其在工程、生物医学等多个领域中的应用。温度响应性水凝胶是一种关键的材料,广泛用于药物输送、组织工程、传感器与执行器以及软机器人领域。该材料主要依赖其在温度变化下的体积膨胀与收缩特性,可实现快速且可逆变形。这种材料的外部刺激包括温度变化、pH 值、电场、磁场等,但温度在各种刺激中具有最广泛的应用,因为它是一种标量、便于控制的参数。
现有研究多集中于均匀温度条件下的水凝胶行为,忽视了热传导方程的求解。然而,现实应用中温度分布往往复杂且非均匀。因此,本研究的目的在于:
1. 探究在温度梯度条件下,温度响应性水凝胶空心圆柱体受拉伸-扭转复合加载时的膨胀行为;
2. 引入 Mazaheri 等人提出的构成模型,通过解析方法和有限元方法展开研究;
3. 通过建立温度分布与应力、形变场的综合分析,深化对温度敏感性水凝胶的力学行为理解。
本研究由以下组成部分构成:
采用 Mazaheri 等人在 2016 年提出的温度敏感性水凝胶的构成模型,该模型基于自由能密度公式,包含拉伸能与混合能的贡献。模型引入了温度 (T) 和变形梯度 (F),主张通过 Helmholtz 自由能的构造描述水凝胶材料的膨胀行为:
[ W = W{\text{stretch}}(T, F) + W{\text{mixing}}(J, T) ]
其中,拉伸能使用 Neo-Hookean 模型描述;混合能则基于扩展后的 Flory-Huggins 模型。此外,模型还考虑了温度和体积比对混合参数的影响,并因此提升了构成模型在邻近相转变温度情况下的稳定性。
研究对象为一个温度敏感性水凝胶空心圆柱,其初态为无应力状态。通过引入径向方向上的变形函数 (r®) 以及拉伸率 (g) 和扭转角度 (τ),定义变形梯度张量 (F)。
进一步引入 Left Cauchy Green 张量 (B) 获取其不变量 (I_1)、(I_2) 和 (J^2)。这为膨胀行为的进一步数学建模提供了关键输入。
第一种情况:等温假设
假定整个结构的温度在加载过程中均匀变化,即所有材料节点温度随时间一致下降。解析方案使用温度作为全局输入,依次求解赫尔姆霍兹方程和应力分布。
第二种情况:非均匀温度分布
通过傅里叶热传导方程求解内部的温度分布: [ \frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial T}{\partial r} \right) = 0 ]
边界条件设定为内半径 (T = 288K),外半径 (T = 300K)。得到的稳态温度分布函数 ((T®)) 会被代入以更新应力与形变场。这种方法允许对多种边界热条件的研究。
在模型验证阶段,通过有限元软件 Abaqus 及其自定义子程序 UHYPER 实现材料定义与模拟计算。针对多种测试用例,包括无约束膨胀和单方向受限膨胀实验,验证理论模型的有效性。
文中开发的分析方法能够进一步扩展到研究功能梯度材料及生物组织的非均匀分布应力行为,并结合新的实验验证算法的通用性。这为智能材料设计,尤其是在医疗设备、柔性机械臂等领域的应用铺平了道路。