橡胶材料的超弹性本构模型建模方法
【专利摘要】本发明提供一种橡胶材料的超弹性本构模型建模方法，该方法首先提出聚合物分子链体积的影响模型；其次，改进八链模型，提出拓扑约束所引起的概率密度函数；然后，建立分子链的微观变形与宏观变形的关系，据此推导橡胶材料的应变能密度函数，得到非仿射超弹性本构模型。所建立的橡胶材料超弹性本构模型同时考虑了聚合物分子链体积、非仿射变形和拓扑约束的影响，建立了宏观的模型参数与微观的分子结构之间的物理联系，能够更深入地反应聚合物大分子力学变形的本质，更准确模拟多功能结构中橡胶材料超弹性特性。
【专利说明】
橡胶材料的超弹性本构模型建模方法
技术领域
[0001] 本发明设及结构力学建模技术领域，具体的设及一种橡胶材料的超弹性本构模型 建模方法。
【背景技术】
[0002] 橡胶材料的力学特性非常复杂，处于橡胶态的聚合物材料柔软而具有弹性，受到 很小的力就可W发生很大的变形，外力除去后拉伸的聚合物分子链又恢复到原来的騰曲状 态，无永久变形，运种特性也称为橡胶材料的高弹性或超弹性。
[0003] 超弹性状态橡胶材料的应力随着应变变化而变化，具有明显的非线性特征。目前， 对于橡胶材料的超弹性特性表征主要有两种途径:一种是基于唯像学理论而提出的超弹性 模型，另一种是基于分子链网络理论而提出的超弹性模型。虽然唯像学模型能够较好地描 述橡胶材料的超弹性特性，但是其模型参数往往缺少明确的物理含义，只能宏观的进行经 验上的数学逼近，不能与材料的微观物理结构联系起来。因此，也有学者基于分子网络理 论，提出了基于分子网络理论的超弹性本构模型。该本构模型将宏观模型参数与材料的微 观分子结构联系起来，W得到具有具体物理含义的结果。但是现有的基于分子网络理论的 超弹性本构模型无法同时考虑分子体积、拓扑约束和非仿射变形的影响。
【发明内容】
[0004] 为解决上述现有技术问题，本发明提供了一种橡胶材料的超弹性本构模型建模方 法。
[0005] 本发明提供一种橡胶材料的超弹性本构模型建模方法，包括W下步骤：
[0006] 1):获取橡胶材料的六个参数:剪切模量μ、非仿射变形系数K、橡胶材料中任一条 聚合物分子链上链节数目Ν、描述分子链体积的对橡胶材料影响程度参数V、聚合物网络拓 扑约束控制参数α与h;
[0007] 2)构建橡胶材料的超弹性本构模型，该步骤包括W下子步骤：
[000引步骤S1:将无变形时聚合物分子链两端的距离ro表示为：
[0009] r〇=n (1)
[0010] 其中，N为一条聚合物分子链共有个链节数目，1为每个链节长度，V描述分子链体 积的影响程度；
[0011] 步骤S2:在其中考虑了聚合物网络的拓扑约束，拓扑约束所引起的概率密度函数 如式(4)所示：
[0012]
(4)
[001引其中，Cc。。为归一化常数，h和α为聚合物网络拓扑约束的控制参数，人1为；个主方 向的宏观拉伸比；
[0014]步骤S3:分子链拉伸比λ。与宏观拉伸比λ"的关系式为：
[0015] λ〇 = Κ(λ"-1)+1 (9)
[0016] 其中，Κ为非仿射变形系数，当Κ = 1时，λ。与λη为传统的仿射变形关系，当Κ声1时，λ。 与λη为非仿射变形关系；
[0017]步骤S4:依据聚合物固体力学理论，根据S1、S2、S3中得到的^)、9[。。(、）、、，推导得 到NA模型的应变能密度函数，并根据应变能密度函数得到橡胶材料超弹性本构模型为：
[001 引
（14)
[0019]其中，μ为剪切模量，Gc〇n=化/8。
[0020] 进一步地，V作为变量且设v〉0.5。
[0021] 进一步地，橡胶材料为乙締基封端聚二甲基硅氧烷。
[0022] 本发明的技术效果：
[0023] 本发明提供的一种橡胶材料的超弹性本构模型建模方法，所得模型既建立了宏观 模型参数与微观分子结构的联系，又同时考虑了分子体积、拓扑约束和非仿射变形的影响， 其建模基础与实际物理情况更为接近，能够准备描述橡胶材料的超弹性特性。
[0024] 具体请参考根据本发明的橡胶材料超弹性本构模型建模方法提出的各种实施例 的如下描述，将使得本发明的上述和其他方面显而易见。
【附图说明】
[0025] 图1是本发明提供的一种橡胶材料超弹性本构模型建模方法流程示意图；
[0026] 图2是本发明优选实施例中改进的八链模型的几何构形结构示意图，其中a)为未 变形时的八链模型，b)为变形后的八链模型；
[0027] 图3是本发明优选实施例中NA模型简单拉伸试验验证结果示意图；
[002引图4是本发明优选实施例中NA模型纯剪切试验验证结果示意图；
[0029] 图5是本发明优选实施例中NA模型双轴拉伸试验验证结果示意图。
【具体实施方式】
[0030] 构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实 施例及其说明用于解释本发明，并不构成对本发明的不当限定。
[0031] 本发明提供的橡胶材料超弹性本构模型建模方法同时考虑分子体积、拓扑约束和 非仿射变形的橡胶材料超弹性本构关系建模方法，建立橡胶材料的非仿射超弹性本构模 型，后文简称NA(Non-Affine)模型，此处的橡胶材料可W为各类橡胶材料，尤其适于处理乙 締基封端聚二甲基硅氧烷材料，该材料为弹性高分子娃橡胶材料，具有可加工性强，可W根 据具体的应用需要，通过模具加注的方式制作成任意形状的特性。其化学分子结构式为：
[0032]
[0033] 参见图1，本发明提供的一种橡胶材料的超弹性本构模型建模方法，包括W下步 骤：
[0034] 1):获取橡胶材料的六个参数:剪切模量μ、非仿射变形系数Κ、橡胶材料中任一条 聚合物分子链上链节数目Ν、描述分子链体积的对橡胶材料影响程度参数V、聚合物网络拓 扑约束控制参数α与h;
[0035] 2)构建橡胶材料的超弹性本构模型，该步骤包括W下子步骤：
[0036] 步骤S1:将橡胶材料无变形时聚合物分子链两端的距离ro表示为：
[0037] r〇 = N"l (1)
[0038] 其中，N为一条聚合物分子链共有个链节数目，1为每个链节长度，V描述分子链体 积的影响程度。
[0039] 在传统当随机游走理论中，V = 0.5。考虑聚合物分子链体积的影响，将V作为变量 且设v〉0.5。按此设置无变形聚合物分子链两端的距离ro所得公式（1)能充分考虑分子体积 对橡胶材料的影响。
[0040] 步骤S2:考虑了聚合物网络的拓扑约束，拓扑约束所引起的概率密度函数如式(4) 所示：
[0041]
(4)
[0042] 其中，Cc。。为归一化常数，h和α为聚合物网络拓扑约束的控制参数，λι为Ξ个主方 向的宏观拉伸比；采用公式(4)表达实现了八链模型的改进，如图2所示。此处的Ξ个主方向 是指所形成坐标系中的Χ、Υ、ΖΞ个方向。
[0043] 步骤S3:分子链拉伸比λ。与宏观拉伸比λ。的关系式为：
[0044] λ〇 = Κ(λ"-1)+1 (9)
[0045] 其中，Κ为非仿射变形系数，当K=1时，λ。与λ。为传统的仿射变形关系，当Κ刮时，λ。 与λ。为非仿射变形关系。采用公式(9)所示的关系表达式，能将微观变形与宏观变形联系起 来。提高解算结果中微观变形与宏观变形的关联度。
[004引步骤S4:依据聚合物固体力学理论，根据S1、S2、S3中得到的；1"日、口(；。。(>0人，推导得 到NA模型的应变能密度函数，并根据应变能密度函数得到橡胶材料应力与应变关系（即本 构模型)为：
[0047]
(14)
[004引其中，μ为剪切模量，Gccn =化/8。相对于传统模型，公式（14)对应的本构模型考虑 了分子体积、拓扑约束和非仿射变形的影响，其建模基础与实际物理情况更为接近，能够准 备描述橡胶材料的超弹性特性。
[0049] 具体的，该方法包括W下步骤：
[0050] S1:考虑分子链体积影响，提出分子链拉伸比λ。的模型。
[0051] 根据聚合物分子链随机游走理论，对于一条无变形的聚合物分子链，其两端的距 离为r。= Vi兩。然而，随机游走理论忽略了分子链体积的影响。考虑聚合物分子链体积的影 响，将无变形聚合物分子链两端的距离ro表示为：
[0052] r〇 = N"l (1)
[0053] 其中，V描述了分子链体积的影响程度。当v = 0.5时，式（1)退化为随机游走理论的 结果;当考虑分子链体积的影响时，有v〉0.5，即，I;, > Vi万,/。
[0054] 利用式（1)得NA模型中的分子链拉伸比、为：
[0055]
口)
[0056] 将(2)式代入采用逆朗之万概率分布来描述的非高斯网络概率密度函数PnG，得PnG 为：
[0化7]
^3)
[0化引其中，CnG为归一化常数，β为朗之万函数。
[0059] S2:根据八链模型，考虑聚合物网络的拓扑约束，提出拓扑约束所引起的概率密度 函数。
[0060] 对于八链模型，当考虑聚合物网络的拓扑约束时，拓扑约束所引起的概率密度函 数采用一个经验模型来表达，如下式所示：
[0061]
(4)
[0062] 其中，Ccnn为归一化常数，h和α为聚合物网络拓扑约束控制参数，λι为Ξ个主方向 的宏观拉伸比。综合考虑非高斯网络概率密度函数PnG和拓扑约束所引起的概率密度函数， 提出NA模型概率密度函数：
[0063]
(5)
[0064] 由Bo 11 zmann方程，计算得到NA模型的构象赌为：
[0065]
(看）
[0066] 其中，kB为Boltzmann常数。根据Helmholtz自由能的表达式％4 庄化其中，T 为绝对溫度。将式(3)、（4)、（5)和(6)代入得：
[0067]
(7)
[006引 WNA表示单个分子链的应变能密度，其中，W0 = -SkBTlnCnG-化BTlnCcon。将WNA乘W分 子链密度η得到单位体积的应变能密度Wna，即：
[0069]
(85
[0070] 其中，y 二址bT为奧切板里，Gcon二化/8 , Wo 二-山nCriDn-Gaussian-0.75山nCcon。
[0071] 步骤S3:式(8)中的λ。为分子链拉伸比，反应的是微观变形，为了将微观变形与宏 观变形联系起来，提出将λ。与宏观拉伸比λ。为非仿射变形关系，即：
[0072] λ〇 = Κ(λ"-1)+1 (9)
[0073] 其中，Κ声1为非仿射变形系数，对于八链模型几何构形有：
[0074]
(1化)
[0075] S3:根据应变能密度函数Wna推导化uchy应力张量
[0076] 在得到NA模型的应变能密度函数Wna之后，根据连续力学理论[参见文献《Treloar L R G.The Physics of Rut)ber Elasticity.Oxford University Press,2005》]，得不可 压缩材料Ξ个主化uchy应力：
[0084] 所建立的ΝΑ模型包含六个参数:μ、K、N、v、Gc。n、α。NA模型能够表征复杂的橡胶材料 超弹性特性，如应力上翅趋势等。采用经典的化eloaH式验数据对提出的ΝΑ模型进行试验验 证，其采用的自然橡胶力学特性能够代表绝大部分橡胶材料的复杂超弹性特性。
[0085] W下面结合附图对一种多功能结构用的橡胶材料的超弹性本构关系建模方法作 详细说明。
[0086] 首先，利用单轴拉伸试验数据对八链模型和NA模型的模型参数分别进行拟合，得： y8 = 0.27MPa、N8 = 25.12、y = 0.2445MPa、K=1.023、N=27.9、v = 0.502、Gc〇n = 0.04]\ffa、a = 0.8602。然后，采用拟合所得参数对纯剪切和双轴拉伸行为的材料变形情况进行预测。采用 TreloaH式验数据时，NA模型和八链模型在简单拉伸变形、纯剪切变形、双轴拉伸变形模式 下的试验验证结果和相对误差分别如图3、图4、图5。从图3和图4可知，最大相对误差都在 20% W内，说明本发明提出的本构模型，能很好的描述八链模型和NA模型下的简单拉伸和 纯剪切行为。从图5可知，八链模型的最大相对误差达到了26%，而NA模型的最大相对误差 只有17%，最小相对误差达到了3%，说明NA模型对于双轴拉伸行为的预测准确性明显优于 八链模型。验证了本发明提出的针对橡胶材料的超弹性本构关系建模方法所构建的具体模 型准确性较高。
[0087] 本领域技术人员将清楚本发明的范围不限制于W上讨论的示例，有可能对其进行 若干改变和修改，而不脱离所附权利要求书限定的本发明的范围。尽管己经在附图和说明 书中详细图示和描述了本发明，但运样的说明和描述仅是说明或示意性的，而非限制性的。 本发明并不限于所公开的实施例。
[0088] 通过对附图，说明书和权利要求书的研究，在实施本发明时本领域技术人员可W 理解和实现所公开的实施例的变形。在权利要求书中，术语"包括"不排除其他步骤或元素， 而不定冠词"一个"或"一种"不排除多个。在彼此不同的从属权利要求中引用的某些措施的 事实不意味着运些措施的组合不能被有利地使用。权利要求书中的任何参考标记不构成对 本发明的范围的限制。
【主权项】
1. 一种橡胶材料的超弹性本构模型建模方法，其特征在于，包括W下步骤： 1) :获取橡胶材料的六个参数:剪切模量y、非仿射变形系数K、所述橡胶材料中任一条 聚合物分子链上链节数目N、描述分子链体积的对所述橡胶材料影响程度参数V、聚合物网 络拓扑约束控制参数a与]1; 2) 构建所述橡胶材料的超弹性本构模型，该步骤包括W下子步骤： 步骤Sl:将无变形时聚合物分子链两端的距离ro表示为： r〇 = N'"l (1) 其中，N为一条聚合物分子链共有个链节数目，1为每个链节长度，V描述分子链体积的 影响程度； 步骤S2:在其中考虑了聚合物网络的拓扑约束，拓扑约束所引起的概率密度函数如式 (4)所示：^4) 其中，Ccnn为归一化常数，h和a为聚合物网络拓扑约束的控制参数，、为S个主方向的宏 观拉伸比； 步骤S3 :分子链拉伸比、与宏观拉伸比Am的关系式为： 入C = K(入 m-l) + l (9) 其中，K为非仿射变形系数，当K= 1时，、与Am为传统的仿射变形关系，当K刮时，、与Am 为非仿射变形关系； 步骤S4:依据聚合物固体力学理论，根据S1、S2、S3中得到的^)、口[。。(、）、、，推导得到魁 模型的应变能密度函数，并根据所述应变能密度函数得到橡胶材料超弹性本构模型为：(14) 其中，y为剪切模量，Gcon =化/8。2. 根据权利要求1所述的橡胶材料的超弹性本构模型建模方法，其特征在于，所述V作 为变量且设V〉0.5。3. 根据权利要求1所述的橡胶材料的超弹性本构模型建模方法，其特征在于，所述橡胶 材料为乙締基封端聚二甲基硅氧烷。
【文档编号】G06F19/00GK105989244SQ201610303626
【公开日】2016年10月5日
【申请日】2016年5月10日
【发明人】廖寰, 廖一寰, 郝东, 李东旭, 刘望, 蒋建平
【申请人】中国人民解放军国防科学技术大学