第34卷第1O期 2018年10月 高分子材料科学与工程 POLYMER MATERlALS SCIENCE AND ENGINEERING Vo1.34.No.10 Oct.2018 纤维增强形状记忆聚合物复合材料的热变形行为 顾建平 ,一,孙慧玉2,张小朋 ,濮 琦3 (1.南京工程学院材料工程学院江苏省先进结构材料与应用技术重点实验室,江苏南京211167; 2.南京航空航天大学机械结构力学及控制国家重点实验室 江苏南京210016; 3.苏州混凝土水泥制品研究院江苏省高耐久混凝土工程技术研究中心,江苏苏州215000) 摘要:通过将热力学内变量理论与Adam—Gibbs结构松弛模型相结合的方法,构建了适用于描述非晶态形状记忆聚合物 (SMP)玻璃化转变过程中热变形行为的模型。在此基础上,利用复合材料细观力学分析方法进一步建立了纤维增强形 状记忆聚合物复合材料(SMPC)热膨胀系数模型。模型中考虑了温度变化对基体弹性模量和泊松比在玻璃化转变过程 中的影响。研究了一种含T300碳纤维增强相SMPC的热变形行为,结果表明,SMPC的各向热应变在升降温阶段几乎 一致,仅在玻璃化转变附近微小分离;较快降温率导致的低温态热应变略小,而较快升温率导致平衡态温度更高;x向热 应变随纤维倾角增大而增大,随纤维含量增大而减小。研究结果可对纤维增强SMPC的设计与应用提供一定理论参考。 关键词:形状记忆聚合物复合材料;纤维;热变形;玻璃化转变;结构松弛 中图分类号:O631.2 2 文献标识码:A 文章编号:1000—7555(2018)10—0059—07 非晶态形状记忆聚合物(SMP)是一种通过对聚合 物进行分子组合和改性得到的特殊高分子材料。研究 表明,SMP及其复合材料的应用前景广阔,遍及航空 航天、生物医疗和纺织工业等领域 J。当前的研究 力学内变量演化方程。相比传统用标量形式表示的结 构松弛模型,该方法获得的是一种三维模型,模型中的 热力学势取决于材料所受的机械载荷、温度历程和热 力学内变量,因此具备更明确的物理意义和更广泛的 适用性。实际上,将Adam—Gibbs结构松弛模型与上 在开发新型SMP和纤维增强形状记忆聚合物复合材 料(SMPC)方面进展较快,然而在构建材料与结构的 述方法结合能更合理、准确地描述其热力学内变量演 化规律。依据此方法,本文研究了非晶态SMP在形状 记忆温度循环中的热变形行为,并以此为基础系统地 研究了纤维增强SMPC的热变形行为,以期能对这类 材料的设计与应用提供理论参考。 1模型建立 1.1 SMP热变形模型 非晶态SMP处于热力学参考平衡态时的温度、应 模型方面进展相对滞后。由于热致SMP和SMPC会 经历热一力学耦合的形状记忆过程,因此在其热一力学 模型的建立过程中,热变形行为的研究较为重要。当 前对SMP热变形行为的研究方法主要分为两类:与相 变理论结合的方法和考虑结构松弛效应的方法[4--7 J。 研究表明,上述SMP热变形研究方法虽能较好地与实 验相符,但存在物理意义不够明确或者通用性不强的 缺点。 基于热力学理论,Adam和Gibbs的结构松弛模型 力张量、应变张量、比熵和内变量张量可依次表示为 ref,丁ref,E ref,s ef和。1 ref,OL2ref‘一,O/mref。在Lion 能较好地描述非晶态聚合物的结构松弛时间在玻璃化 转变过程中的温度相关性|8 J。Lion等 ,“l_在研究非 等_9.“’J结构松弛模型中,当材料所受温度( )和应力 晶态聚合物的结构松弛效应时,通过引入内变量描述 材料所受热一力过程对其当前分子构型影响,然后基于 热力学Clausius—Duhem不等式推导出平衡态附近的热 doi:10.16865/j.cnki.1000.7555.2018.10.O11 (7’)分别达到 和T f并保持不变时,其他变量也将 趋于平衡态时的参考值。该模型适用于描述在参考平 衡态附近时各变量的演变规律。当前状态中各变量取 收稿日期:2017.10.11 基金项目:国家自然科学基金资助项目(11572153);江苏省自然科学基金资助项目(BK20170759,BK20151467);南京工程学院校级科研基 金(ZKJ201603,YKJ201312);绿色建材国家重点实验室开放基金(YA-533);江苏高校优秀科技创新团队 通讯联系人:孙惹玉,主要从事智能材料及三维编织复合材料力学模型的研究,E.mail:hysun@nuaa.edu.cn 60 高分子材料科学与工程 值可定义为: =臼rri+0(t) 、(1 a) 1 =T r+ (,) (1b) E=E f+r( ) =S rct"+,7(t) 。k= k. (1 c) (1d) (1e) ,. f+ k(t) 式中: (f),=(t),r( ),77(z)和 k(t)——均为时 式中:7r _, 和u, k——均为材料参数,且满足: O O O间相关的摄动函数。由于方程(1)中的张量均为二阶 对称张量,下面采用Voigt标记简化表达。例如应力 张量(7、)通常表示为: ±√2 l硼tk— k i (10) A 为温度相关的材料参数矩阵,可表示为: A^=A kI6 O O (1】) r7、11 7’l2 T13] 0 T=f r,、2】 7、22 7、13 l L了、31 T32 7、 』 采用Voigt标记则可表示为: (2) 参数A k与松弛时问rk具备如下关系: k —A kd (12)0 O 0 k 因此,方程(4(c))可改写为: {,,i=[丁l 丁2 7、3 ,』、4 l、6] (3) + + 式中:7 1=7、l】,7、2=T12,7、3:7、33,了、4=7、 7、5= ] ) (14) T23,7、6=7、13。利用Clausius—I)l】}1em不等式可推导得 应变和比熵摄动函数的本构关系以及内变量 (£)的 演化方程如下: -兰 }+{足0}, 0 0 , 0 , { }=D0{ 采用Adam—Gibbs模型,松弛时间 k可表示为: rk= eXpf\ 0S 0 1 10∑u,k{ k} (4a) 式中:rk。、B 和S ——依次为参考松弛时间、激活参 数和构型熵。通过将式(14)中的构型熵变化为比熵可 ,7 一 {ek}.{ +吉 】}. k k 0 0 0 (4b) 得 (5) (6) xp(鲁) { }=一A k(dk{艿k}+{ k} +叫k {∑}) (4c) 0, , O , 0 式中,比熵(S)的摄动项可通过式(4(b))获得。式(4 式中:f。,0——材料的比热容;f。——材料密度;{尼()}和 {ek}——材料参数向量,可表示为: O 0 0 0 0 Z (a))描述的变形耦合了热变形和机械变形两部分。考 虑在三(t):0条件下,材料产生的变形即为热变形, 可表示为: {r}={。(】} 一』。∑Wk{ k} (16) {是o}:[是。是。是0 0 0 0] 0 0 0 0 o {e k}:[Pk k ek 0 0 0] ’ 式中:是[1——玻璃态热膨胀系数;ek——反映了内变量 与温度的耦合效应。D(,——材料在玻璃态时的柔度 1●● ●●●, ● ● ● 矩阵,可表示为: 0 o 0 0 0 考虑材料为各向同性材料,热变形不存在剪切变形部 分,所以相应内变量表示为: { k}={ k】 k2 艿k3 0 0 0} (17) 1.2 SMPC热膨胀系数模型 考虑增强纤维与SMP基体之间的过渡层很薄,因 此可将纤维增强SMI C视为两相复合材料。采用复 合材料细观力学分析方法可推导得到纤维纵向分布 SMPC的轴向热膨胀系数为: (7) —Q'Jr  ̄E【l1u『lI+。f】Ef1【 f 式中:E和 ——分别为弹性模量与泊松比。材料参 数矩阵 k反映了内变量摄动对Gibbs自由能的影响, 可表示为: dt=d kJn (8) (18) 式中:。l1 和af1——分别为基体和纤维轴向热膨胀系 数;E 和Ef ——分别为基体和纤维轴向弹性模量; 和-Uf-一——分别为基体和纤维的体积分数。在此基 式中: k一~材料参数;I ——6阶单位矩阵。材料参 数矩阵议,k反映了内变量与应力之间的耦合效应,可 表示为 础上,考虑横向变形与轴向变形关系并利用混合律可 得横向热膨胀系数: a2=口, l(1+ 。)um+af2(1+, f)uf一日l ft2l (1 9a) 第lO期 顾建平等:纤维增强形状记忆聚合物复合材料的热变形行为 61 tz21= u + fOf (19b) ay=aisin。 +a2COS (23b) 式中:。f2——纤维横向热膨胀系数; 和 f——基体 a y=2(al一口2)sin cos (23c) 和纤维的泊松比。通常来说,常见纤维增强SMPC中 的纤维增强相材料参数受温度影响较小,可以忽略。 式中: ——纤维倾斜,定义为纤维与 轴正向的夹 角。z方向的热膨胀系数由于不受纤维取向的影响, 因此始终与a2相等。 2结果与讨论 2.1模型和材料参数 利用上述模型,本文将从理论上对纤维增强 考虑到玻璃化转变区间的温度变化对SMP基体的材 料参数影响显著,为了讨论温度对a1和a:的影响,必 须首先对其进行研究。 由聚合物理论可知,随着温度的变化,非晶态聚合 物会经历?转变、玻璃化转变和黏性流动这3个阶段。 此处采用Mahieux和Reifsnider[川提出的通用模型来 SMPC的热变形行为进行预报。作为算例,可以采用 描述温度变化对SMP基体弹性模量影响: 文献[5]中所研究的丙烯酸酯基聚合物作为该SMPC E :(E 1一E 2) ̄exp[一f r]十(E 2一E 3). 的SMP基体。选用日本东丽公司的T300碳纤维作 为增强相(材料参数Ef, f, f1和af2见Tab.1),考虑 exp[一( )” ]+E'3"exp[一( )” ] (20) SMPC中纤维体积含量为20%且在基体中均匀分布。 式中:oB、o 和 f——依次分别为p转变温度、玻璃化 转变温度和黏性流动起始温度;E l、E 2和E 3——依 次为上述转变区间的模量初始值;参数 1、 z 2和 z ——材料的we.bul1系数。为了进一步研究温度对 SMP基体泊松比的影响,首先可利用唯象学的方法将 SMP视为玻璃相和橡胶相组成,而后利用Qi等[ ]的 相变模型结合混合律构建SMP泊松比在玻璃化转变 过程中的变化规律: =/zg + (1一.厂g) (21) Fig.2 Parameter fitting for thermal strain(experimental data ean 式中: 和 ——分别为玻璃相和橡胶相的泊松比; be obtained from Fig.2 of Reference[ 1) ——玻璃相的体积分数,可表示为: 1 1一 7 (22) 式中: ——相变参考温度,z——材料参数描述了转 变区间宽度。 Fig.3 Parameter fitting for storage modulus(experimental data can be obtained from Fig.2 of Reference E ]) 需要说明的是,为了便于计算和参数确定,文中 仅取1个热力学内变量 ,而不考虑参数屉的取值对 Fig.1 Relations of coordinate systems of a fiber reinforced 模型描述SMP玻璃化转变的影响。如Fig.2所示,通 composite 过将方程(16)描述的SMP轴向热应变与文献[5]的热 如Fig.1所示,考虑纤维取向对SMPC热变形的 膨胀系数实验数据拟合的方法,获取SMP热变形模型 影响,利用应变转轴公式,可推导得 r平面内任意方 参数(Tab.1中参数0 f,S ref…lD)。拟合过程中模型参 向的热膨胀系数如下: 数的初值确定可参考文献[10]中Tab.1。SMP模量 ax=盘1COS + ̄2 2sin。 (23a) 可分解为储能模量与损耗模量,由于储能模量通常远 62 高分子材料科学与工程 大于损耗模量,故可将SMP弹性模量视为与储能模量 相等。通过将方程(20)与DMA实验数据在Matlab 所示。不失一般性,SMP基体玻璃相 和橡胶相 可分别取为0.35和0.5。与Qi等 l 类似,令相变参 考温度0 :27.5℃,转变区间参数Z=7,从而可得 SMP泊松比随温度变化规律,如Fig.4所示。 中非线性拟合,获取SMP弹性模量随温度变化的模型 参数(Tab.1中参数0p,0g…nz3),拟合结果如Fig.3 Tab.1 Parameters of the modeI and the SMPC 线性变化且两者几乎一致,然而在0 附近两者出现较 大分离且均呈现非线性变化。Fig.6所示为此温度历 程中热力学内变量 的变化情况。从式(16)可知,正 是由于玻璃化转变过程中8的非线性变化才导致了上 述热变形在此区间的显著差异。 Fig·4 Poisson s ratio as a function of temperature 2.2热变形行为预报 下面将利用上述模型预报SMP基体和纤维增强 SMPC在典型的形状记忆温度循环中的热变形行为。 Fig.5所示为SMP基体以2.5℃/min的变温率在l0 ℃到60℃之间循环变化时的热应变演变规律。从图 中可知,SMP基体在高温和低温时的热应变基本符合 Fig.5 Thermal strain as a function of temperature in the thermal cycle 第lO期 顾建平等:纤维增强形状记忆聚合物复合材料的热变形行为 63 向热膨胀系数,并由式(23)可得到 v平面内任意方 向的等效热膨胀系数,进而获取该平面内任意方向的 热应变随温度变化规律。 Fig.7描述了在相同温度历程中,碳纤维体积分 数为2O%且沿轴向均布SMPC的轴向和横向热变形 规律。观察发现,该SMPC在此温度区间发生的轴向 变形非常小,降温和升温阶段的热应变几乎一致,仅在 附近发生微小分离。这主要是由于SMP基体的弹 性模量与碳纤维轴向弹性模量相差悬殊且碳纤维轴向 Fig.6 Internal state variable as a function Of temperature in tile thermal cycle 热膨胀系数为负值所致。该SMPC在此温度区间发 生的横向变形与SMP基体的热变形规律基本一致,也 满足高温和低温时的热应变随温度呈线性变化且几乎 重合,但在0 附近两者发生明显分离。由于SMPC横 向热变形被视为基体和增强相的横向热变形组合而 将SMP基体的等效热膨胀系数(a )定义为: 1一 : (24) 式中: 0——初始温度; ——当前温度;£ ——热应 变。将Fig.5中热应变数值代人式(24),得到此温度 得,该SMPC中基体的含量较高且碳纤维横向热膨胀 系数较小,导致其横向热应变与基体类似但总体略大 区间SMP的等效热膨胀系数。通过将a 代人式 (18)和式(19)获得纤维轴向均布的SMPC轴向和横 Fig.7 Longitudinal strain(a)and transverse strain(b)of fiber reinforced SMPC versus temperature in the thermal cycle 使得低温时的热应变略小于同温条件下较慢降温率对 应的热应变。较快的升温率导致升温过程中材料热变 形达到平衡态所需时间增长,从而使其需在更高的温 度下才能进入平衡态。 Fig.9所示为纤维倾角为30。,45。,60。SMPC的 方向热应变随温度变化情况。从图中可知,各类倾角 的纤维增强SMPC的 方向热应变变化规律均与纯 SMP基体较为类似。z方向热应变随纤维倾角的增 大而增大,也导致在降温和升温阶段中0 附近的热应 Fig.8 X-direction strain of the SMPC at a temperature changing 变分离越明显。 Fig.10所示为纤维倾角为30。,45。SMPC的xy 平面内切向热应变随温度变化曲线。可以发现,不同 倾角SMPC在高温和低温时的热应变基本符合线性 变化且两者几乎一致,但在 附近两者出现明显分离 且均呈现非线性变化。这也是由于该SMPC中基体 rate of 2.5℃/min and 7.5℃/min versus temperature in the thermal cycle 为了研究变温速率对纤维增强SMPC热变形的 影响,Fig.8以纤维体积分数为2O%,倾角45。SMPC 的。,方向热应变为例进行分析。研究表明,较快的降 温率导致材料发生热变形结构松弛的时间较短,从而 体积含量远超纤维含量,导致基体的热变形行为主导 高分子材料科学与工程 了该SMPC的热变形行为。 3 结论 本文构建了适用于非晶态SMP在玻璃化转变区 间的热变形模型并在此基础上建立了纤维增强SMI C 的热膨胀系数模型。文中以一种含T300碳纤维增强 相的SMI,C为例,全面研究了其在形状记忆温度循环 过程中的热变形行为,结论如下: (1)碳纤维轴向均布SMPC在变温区间的轴向变 形非常小,降温和升温阶段的热应变几乎一致,仅在 0 附近发生微小分离;横向变形与SMP基体的热变 形规律基本一致。 l:ig.9 -direction strain of the SMPC with fiber inclination angles of 3O。.45。and 60。versus temperature in the thermal cycle (2)较快的降温率使得碳纤维增强SMPC低温时 的热应变略小于同温条件下较慢降温率对应热应变; 较快的升温率导致材料需在更高的温度下才能进入平 衡态。 (3) 方向热应变随纤维倾角的增大而增大;随碳 纤维含量的增加而减小。 (4)不同纤维倾角SMI C的xy平面内切向热应 变在高低温阶段时基本重合且符合线性变化,仅在 附近出现明显分离且均呈现非线性变化。 参考文献: Fig.10 xy·direction strain of the SMPC with fiber inclination angles of 30。and 45。versus temperature in the thermal cycle [1]Leng J s,LanX,LiuY J,et a1.Shape.memory polymers andtheir composites:stimulus methods and applications[J].Ptog.Mater. Sci.,2011,56:1077一I135. 下面以碳纤维体积含量不同但倾角均为45。 SMPC为例研究纤维体积含量对材料z方向热应变 的影响,模型计算结果如Fig.11所示。可以发现,随 着碳纤维含量的增加,材料.27方向的热应变减小。实 际上,研究发现,其他方向的热应变也都符合这一规 [2]Liu Y J,Ln H B,Lan x,et a1.Review of electro-active shape- memory polymer composite[J].Compos.Sci.Techno1.,2009, 69:2064-2068. [3]LuHB,Yin J Y。YaoYT,et a1.Functionally graded carbon nanotube and nation/silica nanofibre for electrical actuation of carbon fibre reinforced shape memory polymer【j].Pigra.Resin Techno1., 2016,45(2):93.98. 律。这主要是由于碳纤维相对基体具有较小的热膨胀 系数和较大的弹性模量所致。 量0 [4]l iu Y P,Gall K,Dunn M L,et a1.Thermomechania of cshape memory polmerys:uniaxial experiments and constitutive modeling [J].Int.J.Plast.,2006,22:279.313. [5]Westbrook K K,Kao P H,Castro F,et a1.A 3D finite deformation constitutive model for amorphous shape memory polmers:A mulyti '茎 号一0 C branch modeling approach for nonequ ̄brium relaxation processes 兰t-0 一[J].Mech.Mater.,201l,43:853—869. [6]Xiao R,Guo j K,Nguyen T D.Modeling the multiple shape memory effect and temperature memory effcte in amorphous polmers y三 亏 兰一0 () 20 30 40 |;℃ 5O 60 [J].RSC Adv.,2015,5:416-423. [7]顾建平,孙慧玉,方建士,等.热致非晶态形状记忆聚合物的热粘 弹性参数及回复行为[J].高分子材料科学与工程,2016,32(6): 107.112. Fig.11 x—direction Strain of the SMPC with fiber voitmae fraction of 20%,40%and 60%versus temperature in the thermal cycle Gu J P,Sun H Y,Fang J S,“ .Thernmvieodastie parameters and recovery behaviors of thermally activated amorphous shape memory polymers【J].Polymer Materials Science&Engineering, 2016,32(6):107—112. [8]Adam G,Gibbs J H.On the temperature dependence of cooperative 第1O期 顾建平等:纤维增强形状记忆聚合物复合材料的热变形行为 65 relaxation properties in the glass-forming liquids[j].J.Chem. Phys.,1965,43:139-146. thermodynamics with internal state variables[J].Thermochim. Acta,2011,522:182-193. Mahieux C A,Reifsnider K L.Property modeling across transition [9]Lion A,LiebI C,Kolmeder S,et a,.Representation of the glass- transition in mechanical and thermal properties of glass-forming materials:a three-dimensional theory based on thermodynamics with temperatures in polymers:a robust stiffness-temperature model[J]. Polmer,2001,42:3281·3291.y internal state variables[j].j.Mech.Phys.Solids,2010,58: 1338·1360. QiH J,NguyenTD,CastroF。et a1.Finite deformationthermo- mechanical behavior of thermally induced shape memol'y polymers [10]Lion A,Peters J,Kolmeder S.Simulation of temperature history. dependent phenomena of glass—forming materials based on [J].J.Mech.Phys.Solids,2008,56:1730.1751. Thermal Deformation Behaviors of Fiber Reinforced Shape Memory Polymer Composites Jianping Gu ,Huiyu Sun ,Xiaopeng Zhang ,Qi Pu (1.Jiangsu Key Laboratory ofAdvanced Structural Materials and Application Technology,School of Materials Science and Engineering,Nanjing Institute of Technology,Nanjing 211167,China;2.State Key Laboratory f oMechanics and Control of Mechanical Structures,Na ing University of Aeronautics and Astronautics,Na ing 210016,China;3.Suzhou Concrete and Cement Products Research Institute Co.Ltd.,Jiangsu H Durability oncretCe Engineering Technology Research Center,Suzhou 215000,China) ABSTRACT:A therma1 deformation model incorporated with an internal state variable modeling approach based on thermodynamics and Adam—Gibbs structural relaxation model Was developed for amorphous shape memory polymers (SMPs)in the paper.The model can be used to capture the material thermal deformation behavior in the glass transition.Then the thermal expansion coefficient mode1 of fiber reinforced shape memory polymer composites (SMPCs)was proposed based on the mesoscopic mechanical analysis.In the model,the temperature dependence of the elastic modulus and Poisson’s ratio of the SMPs was taken into consideration.The thermal deformation of the SMPC which was reinforced with carbon fiber T300 was studied.It is found that therma1 deformations in different directions are almost the same during heating and cooling,except a slight difference in the glass transition.The thermal deformation in the glassy state decreases with the increase of cooling rate,and the equilibrium temperature increases with the decrease of heating rate.Besides,the therma1 deformation in x direction increases with the increase of the fiber orientation,while it decreases with the increase of the fiber volume content.Therefore,the guidance for design and application of fiber reinforced SMPCs can be provided by the research findings of the paper. Keywords:shape memory polymer composites;fiber;thermal deformation;glass transition;structural relaxation