首页/文章/ 详情

《Mechanics of Solid Polymers》5.3.5 自由连接链模型

6天前浏览36

5.3.5 自由连接链模型


        本章迄今为止的内容主要集中于推导适用于给定变形历史的应力响应的连续力学表达式。具体而言,已经证明如果已知亥姆霍兹自由能的表达式,就可以直接计算应力。亥姆霍兹自由能可以来自现象学的(通常是参数化的)Ψ(·)表达式,或基于微观结构变形行为的微观机制表示。
        本小节介绍一个简单的微观机制启发分子链模型,称为Freely Jointed Chain自由连接链(FJC)模型。文献中还开发了其他微观力学链模型,如Kratky-Porod链模型[17]和蠕虫状链模型[18, 19]。本文不讨论这些模型。
        在FJC模型中,聚合物微观结构中的大分子被表示为n个"刚性链节"(称为Kuhn段),每个长度为l,总链轮廓长度为L = nl,见图5.11。这个模型也被称为随机游走模型,因为假设分子链不相互作用,且其分布对应于随机游走。
        这个分子链模型的亥姆霍兹自由能由各个分子的熵之和给出。这可以通过对可能的分子构型进行统计力学调查获得。

图5.11 随机生成的自由连接链,n = 100,l = 1

        当大分子被拉伸时,能够适应延伸的链端到端距离的可用构型数量减少,从而导致构型熵降低Δη。假设分子是自由连接的,且有固定的键长,可以从统计力学关系确定熵:

其中Ω(r)是分子链端到端距离的概率分布,N是每个参考单位体积中的链数,kB = 1.38 × 10^-23 J/K是玻尔兹曼常数。Flory [21]证明在这些条件下,概率分布可以写为:

对于n ≫ 1且r ≈ n的情况,积分(5.66)的一个良好近似是归因于Kuhn和Grün [22]的朗之万表达式:

其中β = L^-1(r/(nl)),L(x) = coth(x) - 1/x是朗之万函数。逆朗之万函数L^-1(x)无法用初等函数表示,但可以通过以下表达式[23]近似,最大相对误差为6.4 × 10^-4:

逆朗之万函数的图像见图5.12。

注意,在极限n → ∞时,方程(5.67)变为高斯分布:

图5.12 逆朗之万函数L^-1(x)。该函数仅在x值介于-1和+1之间定义


因此,朗之万表达式的链长熵变可以写为:

而对应的高斯链表达式变为:

其中N是每个参考单位体积中的链数,kB是玻尔兹曼常数,r是分子链的端到端距离。

        每单位参考体积的亥姆霍兹自由能定义为:

由于分子链模型中每个链节都是刚性的,因此没有内能储存(e0 ≡ 0),这使我们能够直接得到将分子链拉伸到给定端到端距离所需的力:

        这种微观力学启发的模型可用于更好地理解聚合物的力学响应行为。这种方法也直接应用于一些常用的超弹性材料模型中,如Arruda-Boyce Eight Chain模型(见5.3.10节)。

        接下来的几个小节将讨论一些在有限元程序中常用的超弹性模型。这个综述的目的是说明这些模型之间的一些相似之处,它们如何工作,以及在什么条件下它们能够准确预测真实聚合物材料的行为。

来源:ABAQUS仿真世界
材料
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-04-12
最近编辑:6天前
yunduan082
硕士 | 仿真主任工程... Abaqus仿真世界
获赞 163粉丝 290文章 363课程 0
点赞
收藏
作者推荐

《Mechanics of Solid Polymers》5.3.3 纯剪切与平面拉伸的相似性

5.3.3 纯剪切与平面拉伸的相似性 在特定条件下,平面拉伸试验的结果可等效转换为的纯剪切数据。当需要剪切数据但缺乏专门的剪切试验设备时,这种方法很有用。为说明纯剪切与平面拉伸的相似性,我们从不可压缩材料在纯剪切下的变形梯度开始:其中α为剪切位移。F11项取1/(1-α^2)是为了保持体积不变,即det(F)=1。图5.4展示了纯剪切下,正方形在不同α值时的变形情况。图5.4 展示了纯剪切下正方形在不同α值时的变形:注意,正方形变形成等边菱形。这与简单剪切下正方形变为平行四边形的情况不同现考虑不可压缩平面应变拉伸的变形梯度:其中α与公式(5.59)中的位移值相同。在此变形梯度下,2方向发生拉伸,而1方向无应变。如果我们将这个变形状态绕1方向旋转45°,则旋转矩阵为:旋转后的变形梯度变为:对于小应变,此方程可简化为:其中Δ是一个很小的值,用于保证不可压缩性。O(α^3)是大O记号,表示该项的阶为α^3。如果α≪1,则α^2及更高阶项可忽略,得到:这与纯剪切的变形梯度相同。这表明对不可压缩材料,在小应变下纯剪切和平面拉伸行为是一致的。 随着变形的增加,纯剪切和平面拉伸的差异会变大。图5.5展示了理想化不可压缩加载情况下,真实应力随变形α的变化。图中σ23表示纯剪切应力,σ22表示平面拉伸应力。应力计算采用剪切模量为1MPa的Neo-Hookean(NH)材料模型。图表明,当变形α>0.5时,两种变形模式开始出现明显差异。图5.5 纯剪切和对应平面拉伸在不同应用变形下的预测应力比较来源:ABAQUS仿真世界

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