Soft Matter
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
PAPER
View Article Online View Journal
Cite this: DOI: 10.1039/c4sm01816h
Temperature memory effect in amorphous shape memory polymers† Kai Yu and H. Jerry Qi* Temperature memory effect (TME) refers to the ability of shape memory polymers (SMPs) to memorize the temperature at which pre-deformation was conducted. In the past few years, this TME was experimentally demonstrated by comparing the applied programming temperature (Td) with a characteristic recovery temperature (Tc), which corresponds to either the maximum recovery stress or free recovery speed. In these well-designed experiments, Tc was observed to be close to Td, which is consistent with the intuitive understanding of ‘memorization’. However, since the polymer recovery behavior has been proved to be strongly dependent on various programming and recovery conditions, a new question that whether Tc is always equal to Td in any thermo-temporal conditions remains to be addressed. In this paper, we answered this question by examining the free recovery profile of an acrylate based amorphous SMP. The recovery Tc, which is the temperature with the maximum recovery speed, versus the recovery temperature is shown to be strongly dependent on both programming and recovery conditions. Their detailed influence could be explained by using the reduced time. During a thermomechanical working cycle of SMPs, in addition to the Td, any other thermo-temporal conditions, such as the holding time (th),
Received 16th August 2014 Accepted 18th September 2014
cooling rate, recovery heating rate (q), etc., can affect the observed Tc by changing the reduced programming or recovery time. In this manner, the relationship between Tc and Td is not uniquely determined. Besides, the TME in SMPs can only be achieved within a given temperature range. Both
DOI: 10.1039/c4sm01816h
onset and offset of this temperature range are shown to be influenced by the programming history, but are independent of the recovery conditions.
www.rsc.org/softmatter
1. Introduction The conventional denition of shape memory polymers (SMPs) is “a class of stimuli-responsive smart materials that can be deformed at high temperatures and subsequently be frozen upon cooling to x the programmed temporary shape, which remains stable unless an environmental stimulus is applied, such as temperature,1–8 magnetic elds,9–14 light,15–20 moisture,21,22 etc.”. Depending on whether or not the external load still exists when the SMP is activated, the polymer will demonstrate the shape memory (SM) effect by providing recovery stresses or showing shape recovery, namely a constrained or free recovery process. Among the variety of external stimuli that might be utilized as the triggering mechanism, the most common one is direct heating to increase the temperature. Although the basic concept of SMPs has been known for decades, recent advances have led to the discovery of previously uncovered memory properties that challenge our understanding of this type of polymer material.23–31 For example, the
The George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA. E-mail: [email protected] † Electronic supplementary 10.1039/c4sm01816h
information
(ESI)
available.
This journal is © The Royal Society of Chemistry 2014
See
DOI:
programming condition has been shown to signicantly affect the shape memory behavior of SMPs even if the same thermotemporal recovery condition is used.27,29,30,32–34 A temperature memory effect (TME), in this occasion, refers to the capability of SMPs to memorize the applied programming temperatures (Td). Previous studies in this area demonstrated that both free recovery curves and maximum recovery stress under strain constraint would be shied to a higher recovery temperature if a higher programming temperature was used.4,35 Subsequently, Miaudet et al.28 tested the shape memory property in a polyvinyl alcohol (PVA) based nanocomposite, where the temperature memory is described such that the temperature corresponding to the maximum recovery stress is roughly identical to Td. Xie et al.23,25,36 also reported the free recovery based TME in Naon polymer. The temperature at which a maximum shape recovery rate was observed was also shown to be close to the applied Td. From these pioneer experimental investigations, one can gain recognition that the shape memory characteristic temperature (Tc), which either corresponds to the maximum recovery stress or shape recovery speed, can be quantitatively related to the programming temperature. However, the question remains whether Tc would always be equal to Td in any SM cycle, or their detailed relationship is also dependent on other experimental parameters during the
Soft Matter
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Soft Matter
programming and recovery steps. Many previous studies have shown that the free recovery curves of SMPs are strongly dependent on the packaging programming conditions, such as the loading rate, holding time (th), cooling rate, etc.27,30,37–40 Tobushi et al.41 also experimentally demonstrated the inuence of the recovery heating rate (q) on the free recovery behavior of SMPs. In the case of a higher q being applied, the amount of recovered strain is indeed lower when a given recovery temperature is reached, namely a delayed shape recovery. In this manner, Tc is also dependent on the recovery heating rate. More recently, Wang et al.42 also showed that the free recovery based TME could only be observed within a specic range of recovery temperature, but not in its entire range of glass transition. These observations indicate that the relationship between Tc and Td is not so straightforward as we expected. The detailed programming and recovery conditions during the SM cycle, as well as the polymer temperature dependent thermoviscoelastic properties, may be used to determine Tc during the free recovery of SMPs, and should be considered when we study the TME. In this paper, we studied the TME of thermally responsive amorphous SMPs under the aforementioned considerations, and explored the detailed relationship between Tc and Td during the polymer free shape recovery. Experimental observations rstly revealed the inuence of the programming temperature and heating rate on free shape recovery behavior of an acrylate based SMP material. A one dimensional (1D) multi-branch rheological model was then used to capture the experimental observations, and parametrically predict the free recovery Tc under various programming and recovery conditions. We conclude that both the temperature range to observe the TME in the SMPs, and the relationship between Tc and Td, are not unique. They both are dependent on the thermo-temporal history within a shape memory (SM) cycle. The TME in the SMPs is indeed a synergistic effect of both programming and recovery conditions on the material free recovery behavior.
Paper
2. 2.1
Experimental Materials
The SMP material used in this study is an acrylate-based network polymer. Its liquid monomer tert-butyl acrylate (tBA), liquid cross-linker poly(ethylene glycol) dimethacrylate (PEGDMA) and the photo-polymerization initiator (2,2-dimethoxy-2-phenylacetophenone) in the powder form were ordered from Sigma Aldrich (St Louis, MO, USA) and used in their asreceived conditions without any further purication. The synthesis procedure of the SMP followed Yakacki et al.43 by mixing the tBA monomer, cross-linker and photo-polymerization initiator with a weight ratio of 4000 : 1000 : 1, and the solution was then exposed to UV irradiation for 10 min before the nal heat treatment at 90 C for 1 hour to assist the complete conversion of monomers. All the experiments involved in this paper were conducted on a Dynamic Mechanical Analyzer (DMA, TA Instruments, Model Q800, New Castle, DE, USA). The SMP was cast into rectangular specimens with a uniform dimension of 10 mm 3 mm 0.2 mm. 2.2
DMA test
DMA tests were performed to characterize the glass transition behavior of the acrylate SMPs. The SMP specimen was rstly heated to 95 C on the DMA machine and stabilized for 20 minutes to reach thermal equilibrium. During the experiment, the strain was oscillated at a frequency of 1 Hz with a peak-topeak amplitude of 0.1% while the temperature was decreased from 95 C to 5 C at a rate of 1 C min1. 2.3
Thermomechanical working cycle of SMPs
Free recovery tests of acrylate SMP were conducted on the DMA machine by using strain controlled methods, and the detailed procedure followed the SM cycle shown in Fig. 1a. In the programming step, the SMP was stretched to a target strain emax (20%) with a constant loading rate ˙e (0.01 s1) at the
Fig. 1 A schematic illustration of (a) the thermomechanical history of the programming and free recovery process in a SM cycle, and (b) the characteristic free recovery temperature Tc.
Soft Matter
This journal is © The Royal Society of Chemistry 2014
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Paper
Soft Matter
programming temperature Td, followed by a specied holding time (th) at Td before being cooled to the shape xing temperature TL (20 C) at a rate of 2.5 C min1. Once TL was reached, the specimen was held for one hour then the tensile force was removed. In the free recovery step, the temperature was increased to the recovery temperature Tr at a given heating rate of q. The strain evolution during this step was used to calculate the time dependent shape recovery ratio. Throughout the paper, the shape recovery ratio Rr and the shape recovery speed Vr are respectively dened as25,44 Rr ¼ 1 e(t)/er
(1a)
and Vr ¼
vRr ; vTH
(1b)
where e(t) is the strain evolution during the free recovery step, er is the strain measured at the start of the free recovery process, TH(t) ¼ qt is the time dependent recovery temperature, and the characteristic recovery temperature Tc is taken to be the temperature corresponding to the maximum recovery speed (see Fig. 1b).
3.
Constitutive modeling of SMPs
Currently, theoretical studies of shape memory phenomena in polymers are mainly based on two modeling approaches: phase evolution45–50 and thermoviscoelastic modeling method.27,51–55 Here in this study, to capture the shape memory properties of the acrylate-based SMP, we use a multi-branched thermoviscoelastic constitutive model. The model resembles a generalized Maxwell element that consists of one equilibrium branch and several thermoviscoelastic nonequilibrium branches (number of n) to represent the multiple relaxation processes (for example, Rouse modes) of the polymer. Eeq and Ei are the elastic moduli in the equilibrium and nonequilibrium branches, respectively. si(T) denotes the temperature dependent relaxation time in the dashpots. A detailed description of the multibranch model for SMP application can be found in the study by
Yu et al.27 which is also shown in the ESI† of this paper. Its extension for a 3D nonlinear model can be found in Westbrook et al.54 and Diani et al.53 Here, it is assumed that the time–temperature shi for each branch follows the same rule. According to the well-established “thermorheological simplicity” principle56 under a nonisothermal condition, the relaxation times (or viscosity) of each nonequilibrium branch vary as the temperature changes: si ¼ aT(T)s0i ,
(2)
where aT(T) is a time–temperature superposition (TTSP) shiing factor, and s0i is the reference relaxation time at the temperature when aT(T) ¼ 1. Depending on whether the temperature is above or below a shiing temperature Ts,57 the shi factor is respectively calculated by the following two equations, namely the Williams– Landel–Ferry (WLF) equation58 and the Arrhenius-type equation:59 C1 ðT Tr Þ log aT T ¼ ; C2 þ ðT Tr Þ
(3a)
AFc 1 1 ln aT T ¼ : Tg kb T
(3b)
where Tr is the WLF reference temperature, Tg is the glass transition temperature of SMP, C1 and C2 are material constants that depend on the choice of Tr, A is a material constant, Fc is the congurational energy and kb is Boltzmann's constant. Ts is taken to be the crossing point of two curves representing eqn (3a) and (3b) on an aT(T) vs. T plot. The model parameters were determined by using tan d and storage modulus of the SMP obtained from the DMA tests. The identication procedure is briey introduced in the ESI.† A detailed model parameters identication method of the multibranch model can be found in previous studies.27,52,53 It should be noted that in this paper, the applied modeling method is based on linear thermoviscoelasticity. Caution should be taken if a SMP demonstrates highly nonlinear or inelastic behavior. The small strain DMA test will be insufficient to characterize the
(a) 1D estimation for the DMA result of a SMP by using NLREG software. The original experimental curves are shown in solid lines. (b) Shifting factors at different temperatures: above 27.5 C, shifting factors follow the WLF equation; below 27.5 C, shifting factors follow the Arrhenius-type behavior. Fig. 2
This journal is © The Royal Society of Chemistry 2014
Soft Matter
View Article Online
Soft Matter
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Table 1
Paper
Parameters of the multi-branched constitutive model
Parameters
Values
C1 AFckb1 Ts Eeq Ei
16.44 C2 51.6 C 40 000 K Tg 42.7 C 27.5 C TM 30 C 3.05 MPa (424.2, 424.2, 399.2, 276.6, 329.7, 245.3, 153.4, 64.4, 23.5, 7.51, 3.02, 1.50, 0.71, 0.43, 0.12, 0.16) MPa for i ˛ [1,16] (100.5, 100, 100.5, 101, 101.5, 102, 102.5, 103, 103.5, 104, 104.5, 105, 105.5, 106, 106.5, 107) s, for i ˛ [1,16]
s0i
Parameters
Values
thermoviscoelastic (or thermoviscoplastic) behavior and thus may not be able to provide an accurate estimation of the SM performance. In this case, a nonlinear modeling approach such as those presented by Westbrook et al.,54 Abrahamson et al.,60 Tobushi et al.61 and Srinivasa et al.62 can be used. As shown in Fig. 2a, the identied parameters could enable the multi-branch model to successfully capture the experimental storage modulus and tan d curves within the entire testing temperature range (5–95 C). The temperature dependent shi factors were then calculated by using eqn (3) and are plotted in Fig. 2b as a function of temperature. The nal set of model parameters is listed in Table 1. Compared with our previous work, where seven nonequilibrium branches are constructed in the multi-branch model, the amount of nonequilibrium branches in this study is extended to sixteen with rened relaxation time distribution, which will enable a smooth and accurate prediction on the recovery speed curves versus the recovery temperature.
4. Results and discussion 4.1
Experimental observations
We start our discussion by rstly showing some experimental observations on the free recovery ratio and speed under different programming and recovery conditions. As a baseline SM cycle, the programming temperature Td, the holding time th at Td and the recovery heating rate q are respectively set to be 35 C, 60 s and 1 C min1. The experimental setup is the same as
that described in Section 2.3. As a comparison, two more SM cycles were performed: one with an increased programming temperature (Td ¼ 40 C), while the other one with a decreased heating rate (q ¼ 0.1 C min1). The shape recovery ratio and speed in these three SM cycles are respectively plotted in Fig. 3a and b as a function of recovery temperature (TH). As shown in Fig. 3, both the programming temperature and the heating rate signicantly affect the material free recovery prole: on one hand, increasing Td from 35 C to 40 C decreases the shape recovery speed under the same heating rate, and elevates Tc from 34 C to 38 C. On the other hand, similar to the previous study by Tobushi et al.,41 decreasing q increases the recovery speed versus the recovery temperature, and consequently shis Tc to a small value. This is because a lower heating rate, which means a longer heating time before reaching a given recovery temperature, would allow the SMPs to evolve more towards the equilibrium state and hence give a larger recovery strain. When the heating rate is set to be 0.1 C min1, Tc is seen to decrease from 34 C to 30 C. In addition to the experimental observations, simulation results based on the multi-branch model and predetermined model parameters are also plotted in Fig. 3 as dashed lines. It can be seen that the numerical simulations adequately capture the temperature dependent shape recovery ratio and normalized recovery speed. In the following, parametric studies based on the multi-branch model will be performed to investigate the inuence of programming and recovery conditions on the TME in SMPs. 4.2
Inuence of the programming condition
The inuence of Td on the polymer free recovery prole is predicted and demonstrated in Fig. 4a and b, where Td is gradually increased from 25 C to 59 C (2 C temperature interval) during the simulation. The holding time th at Td and the heating rate during the free recovery are respectively set to be 60 s and 1 C min1 in each case. Fig. 4c further plots the Tc during the polymer free recovery as a function of applied Td. Two main features could be found from the numerical simulations in Fig. 4. Firstly, within a specic temperature range (30 C to 47.5 C), increasing the programming temperature will decrease the recovery speed versus the recovery
Fig. 3 Influence of the programming temperature and heating rate on the (a) shape recovery ratio and (b) normalized shape recovery speed. Note: the dashed lines denote simulation curves by using the 1D multi-branch model.
Soft Matter
This journal is © The Royal Society of Chemistry 2014
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Paper
Soft Matter
Fig. 4 Parametric study on the influence of the programming temperature. Td increases from 25 C to 59 C with a 2 C temperature interval. (a) Shape recovery ratio as a function of recovery temperature. (b) Shape recovery speed as a function of recovery temperature. (c) Tc as a function of Td.
temperature. Tc increases with Td in exhibiting a roughly linear relationship. This obeys the recognition of the traditional TME that Tc observed during the free recovery step is close to the Td applied in the programming step. Secondly, beyond this temperature range, changing Td does not obviously affect Tc. This prediction, that the TME could only be realized within a given temperature range, is consistent with the recent study of Wang et al.42 The two features shown in Fig. 4 could be explained by examining the stress distribution within each branch of the multi-branch model. Previously, we demonstrated that if the same thermo-temporal recovery condition is applied on SMPs, their free recovery proles are only dependent on the stress distribution aer unloading, or namely the stress distribution at the beginning of recovery.27 Fig. 5a shows this stress distribution at the end of seven different programming conditions with different Td values (27.5 C, 30 C, 40 C, 47.5 C and 60 C respectively). At a temperature higher than 47.5 C, Tc is observed to be independent of Td. This phenomenon can be explained as follows: when the SMPs are stretched at higher temperatures, the elastic strain development within each nonequilibrium branch is small due to the high mobility of dashpots. The stress relaxation in the subsequent holding and cooling steps will further release these nonequilibrium elastic strains. During the unloading at low temperatures, removing the external load will cause a re-distribution of strains (due to the rebalance of forces
This journal is © The Royal Society of Chemistry 2014
to maintain overall equilibrium) in all the branches,24,27 and an internal stress state aer unloading is that only the equilibrium spring remains in tension while the rest springs are being compressed (see Fig. 5a for the cases of Td ¼ 47.5 C and 60 C). In this manner, when the temperature is raised again during the recovery step, only the equilibrium branch provides driven force for viscous strain release (see Fig. 5b) and material shape recovery, while the rest branches keep on resisting until the contained dashpots are activated. Since the nonequilibrium elastic stress is already released essentially before unloading for the case of Td ¼ 47.5 C, further increasing the programming temperature to 60 C will not affect the rheological stress state, and therefore, the subsequent free recovery curve of SMPs is almost the same. This is the reason for Td ¼ 47.5 C being observed as the upper temperature limit for the TME in the acrylate SMP as seen in Fig. 4c. However, when the programming temperature is gradually decreased from 47.5 C to 30 C, mobility of the dashpots decreases, and considerable elastic strain could be developed when stretching the multi-branch model, especially in branches with larger relaxation times. These elastic strains cannot be totally released during the subsequent holding and cooling, and will be stable when the temperature is low. As seen in Fig. 5a, decreasing the programming temperature will change the rheological stress state aer unloading in two ways: rstly, the stress magnitude in each nonequilibrium branch is increased. The intensied stress state will proportionally increase the
Soft Matter
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Soft Matter
Paper
(a) Stress distribution after unloading in seven different programming conditions. (b) Viscous strain release in each nonequilibrium branch and strain evolution in the equilibrium branch for the case of Td ¼ 47.5 C and Td ¼ 30 C.
Fig. 5
releasing rate of viscous strain in the dashpots (see the bottom gure in Fig. 5b). Secondly, aer the stress re-distribution among all the branches during the unloading, more and more nonequilibrium springs will remain stretched in addition to the equilibrium spring. These branches will provide additional driven forces for the shape recovery of the multi-branch model, and consequently increase the recovery speed. In the rheological multi-branch model, this additional driven force, namely the tensile stress in the nonequilibrium branches at the beginning of recovery, essentially represents the internal energy stored in the real polymer system aer the programming step. When stretching the polymer network at lower temperatures, the mobility of polymer chains might be insufficient for them to adjust the self-conformation and fully comply with the external deformation. In this manner, in addition to the network entropy change, some internal energy will also be generated and stored within the SMPs. This internal energy originates from the alteration in the macromolecule bond length and angle during the deformation, and will assist the spontaneous shape recovery of the polymer network upon reheating. With regard to the inuence of the programming temperature as revealed in Fig. 5a, it seems that decreasing the programming temperature is a good choice to increase the shape recovery speed. But the material's shape xing ability will be compromised because the material tends to be unstable during unloading and a large portion of the programmed deformation is recovered right aer unloading,27 which is unfavorable in most engineering applications. For example, when Td ¼ 30 C, only 41% of the programmed strain is xed aer unloading. Further decreasing the programming temperature will continually decrease the shape xity of the polymer, but does not affect the stress state too much. In this manner, the recovery Tc will remain unchanged when Td is below 30 C.
Soft Matter
Now it is clear that the rheological stress state at the beginning of free recovery is the key to determine Tc. Previously, we have demonstrated that this stress state is only dependent on the reduced programming time tr.27 For a non-isothermal process, ðt ds ; (4) tr ¼ aðTÞ 0 where t is the physical time and T is the temperature evolution during the programming step. Therefore, in addition to Td, any other thermo-temporal programming conditions, such as holding time th aer loading, cooling rates, shape xing temperature TL, etc., can also affect the recovery speed and the Tc by changing the reduced programming time and then consequently the rheological stress state. Here, we select th as the parametric variable, and the predicted free recovery curves are shown in Fig. 6a. Td and the heating rate q are set to be 35 C and 1 C min1 respectively in each simulation case. As shown in the gure, increasing the holding time th will decrease the recovery speed because it plays the same role as the increasing Td in determining the reduced programming time. Besides, two boundary holding times are revealed from the gure: one represents the extreme condition that no holding time is applied (black curve). While for the second one, it indicates that a 6 104 s (17 hours) holding time at Td ¼ 35 C is long enough to essentially release the elastic stresses in nonequilibrium branches during the programming, and the corresponding recovery prole is therefore the same as that with th ¼ 6 108 s (19 years). Subsequently, Fig. 6b summarizes the recovery Tc as a function of both Td and the holding time th. Clearly, the relationship between Tc and Td is not unique under the inuence of th. A longer holding time will increase the Tc towards a saturated level at 47.5 C. Besides, both the onset and offset values of the
This journal is © The Royal Society of Chemistry 2014
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Paper
Soft Matter
Fig. 6 Influence of both Td and holding time on the free recovery profiles and Tc of SMPs. (a) The predicted free recovery curves with different holding times. Td ¼ 35 C. (b) Tc plotted as a function of both holding time and Td. (c) The constructed master curve of Tc. The reference holding time th0 ¼ 6 101 s. (d) The holding time dependent shift factors during the construction of the master curve.
temperature range with the TME are observed to decrease with increment of the holding time. For the curves in Fig. 6b, aer shiing each of them horizontally to a reference holding time th0, namely, by Tc(Td,th) ¼ Tc(athTd,th0),
(5a)
where ath is the holding time dependent shi factor, a master curve can be constructed as shown in Fig. 6c. The corresponding shi factors are plotted in Fig. 6d as a function of holding time. A linear relationship between ath and th is then revealed in the semi-log scale. By using the nonlinear regression method, we approximate the relationship as th ath ¼ 0:85 þ 0:083 log (5b) th0 Within the temperature range with a clear TME observed, eqn (5) quantitatively demonstrates the packaging effect of both the programming temperature and the holding time on Tc. That is, when a large holding time is applied during the programming step, the corresponding shi factor will be exponentially ramped and therefore a smaller Td should be used if an equivalent Tc is expected. Otherwise Tc will monotonically increase. The analysis on the effect of increasing Td and th can be extended to the inuence of other thermo-temporal
This journal is © The Royal Society of Chemistry 2014
programming conditions, such as decreasing loading rate, cooling rates, or increasing the stabilization time at TL, etc. 4.3
Inuence of recovery conditions
In addition to the programming condition, we will show in this section that the relationship between Tc and Td is also dependent on the free recovery condition. Considering the SM cycle shown in Fig. 1a, the only parametric variable under investigation is the recovery heating rate q. Fig. 7a shows the predicted free recovery curves with q increasing from 103 C min1 to 102 C min1. All the programming procedures are set to be identical in each simulation case with Td ¼ 35 C and a holding time equal to 60 s. As demonstrated in our previous study,27 if the same programming condition is used, the free recovery ratio of SMPs is only dependent on the reduced recovery time, and a longer reduced recovery time will lead to a higher recovery ratio, namely Rr ¼ Rr(tr)
(6a)
with tr ¼
ð TH TL
dT ; qaðTÞ
(6b)
where TH is a given recovery temperature. Eqn (6b) indicates that a larger heating rate will decrease the reduced recovery
Soft Matter
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Soft Matter
Paper
Fig. 7 Influence of both Td and heating rate q on the free recovery profiles and Tc of SMPs. (a) The predicted free recovery curves with different heating rates. Td ¼ 35 C and holding time equals to 6 s. (b) Tc plotted as a function of both Td and q. (c) The constructed master curve of Tc. The reference heating rate q0 ¼ 100 C min1. (d) The heating rate dependent shift factors during the construction of the master curve.
time and then consequently the recovery ratio when TH is reached, as conrmed in Fig. 7a. Fig. 7b summarizes the recovery Tc as a function of both Td and q. In addition to the revealed relationship between Tc and Td that is dependent on q, the temperature range to observe the TME is shown to be unchanged upon changing the heating rate. Subsequently, by shiing each curve in Fig. 7b vertically to a reference heating rate (q0 ¼ 100 C min1) according to Tc ðTd ; qÞ ¼
1 Tc ðTd ; q0 Þ aq
(7a)
where aq is a heating rate dependent shi factor, the curves are seen to be rejoined together to form a master curve, as shown in Fig. 7c. The corresponding shi factors are summarized in Fig. 7d. In the semi-log scale, the shi factor decreases linearly with increment of the heating rate, namely q aq ¼ 1 0:092 log (7b) q0 Eqn (7) quantitatively approximates the inuence of the recovery heating rate on Tc with two features revealed: rstly, when a larger heating rate q is applied during the polymer free recovery step, both the shi factor and the corresponding Tc will decrease. Secondly, different from the above-mentioned effect of programming conditions, there is no boundary heating rate condition that inuences the heating rate. Continuous increase of q will continually elevate Tc towards the target TH. An
Soft Matter
innitely large value represents an isothermal recovery condition, where Tc is always equal to TH regardless of the applied programming temperature. In contrast, if no heat is applied during the recovery step, namely q ¼ 0, a thermally simple SMP can still recover to its original shape, even though this process may take several years.8,27 Under this condition, Tc is equal to the shape xing temperature TL.
5.
Conclusion
In this paper, from both the experimental and theoretical point of view, we investigated the temperature memory effect (TME) in the acrylate based amorphous shape memory polymer (SMP). Different from the traditional denition of TME in the SMPs, our study indicates that the characteristic recovery temperature (Tc), which is the temperature with the maximum recovery speed, is not always equal to the programming temperature (Td). By revealing the underlying mechanism, we can gain a deeper understanding of how the programming and recovery conditions will affect the free recovery behavior of SMPs synergistically, which consequently provides a facile strategy to optimize the shape memory performance of SMPs by adjusting their thermomechanical working conditions, instead of developing new polymer systems. Such a strategy offers huge convenience for engineering SMPs in the industry, especially in the elds of medicine, aeronautics and astronautics, where a long period of performance verication and official approval is typically expected.
This journal is © The Royal Society of Chemistry 2014
View Article Online
Paper
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
References 1 A. Lendlein and S. Kelch, Shape-memory polymers, Angew. Chem., Int. Ed. Engl., 2002, 41(12), 2035–2057. 2 C. Liu, H. Qin and P. T. Mather, Review of progress in shapememory polymers, J. Mater. Chem., 2007, 17(16), 1543–1558. 3 T. Chung, A. Rorno-Uribe and P. T. Mather, Two-way reversible shape memory in a semicrystalline network, Macromolecules, 2008, 41(1), 184–192. 4 K. Gall, et al., Thermomechanics of the shape memory effect in polymers for biomedical applications, J. Biomed. Mater. Res., Part A, 2005, 73(3), 339–348. 5 P. T. Mather, X. F. Luo and I. A. Rousseau, Shape Memory Polymer Research, Annu. Rev. Mater. Res., 2009, 39, 445–471. 6 K. Yu, et al., Carbon nanotube chains in a shape memory polymer/carbon black composite: To signicantly reduce the electrical resistivity, Appl. Phys. Lett., 2011, 98(7), 0741021–0741023. 7 K. Yu, Y. J. Liu and J. S. Leng, Conductive Shape Memory Polymer Composite Incorporated with Hybrid Fillers: Electrical, Mechanical, and Shape Memory Properties, J. Intell. Mater. Syst. Struct., 2011, 22(4), 369–379. 8 Q. Ge, et al., Prediction of temperature-dependent free recovery behaviors of amorphous shape memory polymers, So Matter, 2012, 8(43), 11098–11105. 9 R. Mohr, et al., Initiation of shape-memory effect by inductive heating of magnetic nanoparticles in thermoplastic polymers, Proc. Natl. Acad. Sci. U. S. A., 2006, 103(10), 3540–3545. 10 P. R. Buckley, et al., Inductively heated shape memory polymer for the magnetic actuation of medical devices, IEEE Trans. Biomed. Eng., 2006, 53(10), 2075–2083. 11 A. M. Schmidt, Electromagnetic activation of shape memory polymer networks containing magnetic nanoparticles, Macromol. Rapid Commun., 2006, 27(14), 1168–1172. 12 Z. W. He, et al., Remote Controlled Multishape Polymer Nanocomposites with Selective Radiofrequency Actuations, Adv. Mater., 2011, 23(28), 3192–3196. 13 U. N. Kumar, et al., Non-contact actuation of triple-shape effect in multiphase polymer network nanocomposites in alternating magnetic eld, J. Mater. Chem., 2010, 20(17), 3404–3415. 14 K. Yu, et al., Design considerations for shape memory polymer composites with magnetic particles, Adv. Mater., 2013, 47(1), 51–63. 15 H. Y. Jiang, S. Kelch and A. Lendlein, Polymers move in response to light, Adv. Mater., 2006, 18(11), 1471–1475. 16 H. Koerner, et al., Remotely actuated polymer nanocomposites – stress-recovery of carbon-nanotube-lled thermoplastic elastomers, Nat. Mater., 2004, 3(2), 115–120. 17 A. Lendlein, et al., Light-induced shape-memory polymers, Nature, 2005, 434(7035), 879–882. 18 M. H. Li, et al., Light-driven side-on nematic elastomer actuators, Adv. Mater., 2003, 15(7–8), 569–572.
This journal is © The Royal Society of Chemistry 2014
Soft Matter
19 T. F. Scott, R. B. Draughon and C. N. Bowman, Actuation in crosslinked polymers via photoinduced stress relaxation, Adv. Mater., 2006, 18(16), 2128–2132. 20 T. F. Scott, et al., Photoinduced plasticity in cross-linked polymers, Science, 2005, 308(5728), 1615–1617. 21 W. M. Huang, et al., Water-driven programmable polyurethane shape memory polymer: Demonstration and mechanism, Appl. Phys. Lett., 2005, 86(11), 1141051– 1141053. 22 Y. C. Jung, H. H. So and J. W. Cho, Water-responsive shape memory polyurethane block copolymer modied with polyhedral oligomeric silsesquioxane, J. Macromol. Sci., Part B: Phys., 2006, 45(4), 453–461. 23 T. Xie, Tunable polymer multi-shape memory effect, Nature, 2010, 464(7286), 267–270. 24 K. Yu, et al., Mechanisms of multi-shape memory effects and associated energy release in shape memory polymers, So Matter, 2012, 8(20), 5687–5695. 25 T. Xie, K. A. Page and S. A. Eastman, Strain-Based Temperature Memory Effect for Naon and Its Molecular Origins, Adv. Funct. Mater., 2011, 21(11), 2057–2066. 26 T. Xie, X. C. Xiao and Y. T. Cheng, Revealing Triple-Shape Memory Effect by Polymer Bilayers, Macromol. Rapid Commun., 2009, 30(21), 1823–1827. 27 K. Yu, G. Qi and H. J. Qi, Reduced time as a unied parameter determining xity and free recovery of shape memory polymers, Nat. Commun., 2014, 5, 3066. 28 P. Miaudet, et al., Shape and temperature memory of nanocomposites with broadened glass transition, Science, 2007, 318(5854), 1294–1296. 29 X. Chen and T. D. Nguyen, Inuence of thermoviscoelastic properties and loading conditions on the recovery performance of shape memory polymers, Mech. Mater., 2011, 43, 127–138. 30 F. Castro, et al., Time and Temperature Dependent Recovery of Epoxy-Based Shape Memory Polymers, J. Eng. Mater. Technol., 2011, 133(2), 0210251–0210259. 31 L. Sun and W. M. Huang, Mechanisms of the multi-shape memory effect and temperature memory effect in shape memory polymers, So Matter, 2010, 6(18), 4403–4406. 32 J. L. Hu, F. L. Ji and Y. W. Wong, Dependency of the shape memory properties of a polyurethane upon thermomechanical cyclic conditions, Polym. Int., 2005, 54(3), 600–605. 33 J. Cui, K. Kratz and A. Lendlein, Adjusting shape-memory properties of amorphous polyether urethanes and radioopaque composites thereof by variation of physical parameters during programming, Smart Mater. Struct., 2010, 19(6), 0650191–0650199. 34 M. Bothe, F. Emmerling and T. Pretsch, Poly(ester urethane) with Varying Polyester Chain Length: Polymorphism and Shape-Memory Behavior, Macromol. Chem. Phys., 2013, 214(23), 2683–2693. 35 Y. Miyamoto, et al., Memory effect on the glass transition in vulcanized rubber, Phys. Rev. Lett., 2002, 88(22), 2255041– 2255044.
Soft Matter
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Soft Matter
36 J. J. Li and T. Xie, Signicant Impact of Thermo-Mechanical Conditions on Polymer Triple-Shape Memory Effect, Macromolecules, 2011, 44(1), 175–180. 37 X. Chen and T. D. Nguyen, Inuence of thermoviscoelastic properties and loading conditions on the recovery performance of shape memory polymers, Mech. Mater., 2011, 43(3), 127–138. 38 A. McClung, G. Tandon and J. Baur, Strain rate- and temperature-dependent tensile properties of an epoxybased, thermosetting, shape memory polymer (Veriex-E), Mech. Time-Depend. Mater., 2011, 16(2), 205–221. 39 A. McClung, G. Tandon and J. Baur, Deformation rate-, hold time-, and cycle-dependent shape-memory performance of Veriex-E resin, Mech. Time-Depend. Mater., 2013, 17(1), 39–52. 40 A. J. W. McClung and M. B. Ruggles-Wrenn, Strain Rate Dependence and Short-Term Relaxation Behavior of a Thermoset Polymer at Elevated Temperature: Experiment and Modeling, J. Pressure Vessel Technol., 2009, 131(3), 0314051–0314058. 41 H. Tobushi, et al., Thermomechanical constitutive model of shape memory polymer, Mech. Mater., 2001, 33(10), 545–554. 42 Y. R. Wang, et al., Relation between temperature memory effect and multiple-shape memory behaviors based on polymer networks, RSC Adv., 2014, 4(39), 20364–20370. 43 C. M. Yakacki, et al., Unconstrained recovery characterization of shape-memory polymer networks for cardiovascular applications, Biomaterials, 2007, 28(14), 2255–2263. 44 A. Lendlein and R. Langer, Biodegradable, elastic shapememory polymers for potential biomedical applications, Science, 2002, 296(5573), 1673–1676. 45 Y. P. Liu, et al., Thermomechanics of shape memory polymers: Uniaxial experiments and constitutive modeling, Int. J. Plast., 2006, 22(2), 279–313. 46 W. W. Muller and T. Pretsch, Hydrolytic aging of crystallizable shape memory poly(ester urethane): Effects on the thermo-mechanical properties and visco-elastic modeling, Eur. Polym. J., 2010, 46(8), 1745–1758. 47 V. Kaa, Shape memory polymers: A mesoscale model of the internal mechanism leading to the SM phenomena, Int. J. Plast., 2008, 24(9), 1533–1548. 48 Y. C. Chen and D. C. Lagoudas, A constitutive theory for shape memory polymers. Part I-Large deformations, J. Mech. Phys. Solids, 2008, 56(5), 1752–1765.
Soft Matter
Paper
49 K. N. Long, M. L. Dunn and H. J. Qi, Mechanics of so active materials with phase evolution, Int. J. Plast., 2010, 26(4), 603– 616. 50 H. J. Qi, et al., Finite deformation thermo-mechanical behavior of thermally induced shape memory polymers, J. Mech. Phys. Solids, 2008, 56(5), 1730–1751. 51 J. R. Lin and L. W. Chen, Shape-memorized crosslinked ester-type polyurethane and its mechanical viscoelastic model, J. Appl. Polym. Sci., 1999, 73(7), 1305–1319. 52 K. Yu, et al., A thermomechanical constitutive model for an epoxy based shape memory polymer and its parameter identications, Mech. Time-Depend. Mater., 2014, 18(2), 453–474. 53 J. Diani, et al., Predicting thermal shape memory of crosslinked polymer networks from linear viscoelasticity, Int. J. Solids Struct., 2012, 49(5), 793–799. 54 K. K. Westbrook, et al., A 3D Finite Deformation Constitutive Model for Amorphous Shape Memory Polymers: A MultiBranch Modeling Approach for Nonequilibrium Relaxation Processes, Mech. Mater., 2011, 43(12), 853–869. 55 T. D. Nguyen, et al., A thermoviscoelastic model for amorphous shape memory polymers: Incorporating structural and stress relaxation, J. Mech. Phys. Solids, 2008, 56(9), 2792–2814. 56 M. Rubinstein and R. H. Colby, Polymer Physics, Oxford University Press, Oxford, New York, 2003. 57 P. A. O'Connell and G. B. McKenna, Arrhenius-type temperature dependence of the segmental relaxation below T-g, J. Chem. Phys., 1999, 110(22), 11054–11060. 58 M. L. Williams, R. F. Landel and J. D. Ferry, Temperature Dependence of Relaxation Mechanisms in Amorphous Polymers and Other Glass-Forming Liquids, Phys. Rev., 1955, 98(5), 1549. 59 E. A. Di Marzio and A. J. M. Yang, Congurational entropy approach to the kinetics of glasses, J. Res. Natl. Inst. Stand. Technol., 1997, 102(2), 135–157. 60 E. R. Abrahamson, et al., Shape memory mechanics of an elastic memory composite resin, J. Intell. Mater. Syst. Struct., 2003, 14(10), 623–632. 61 H. Tobushi, et al., Thermomechanical constitutive modeling in shape memory polymer of polyurethane series, J. Intell. Mater. Syst. Struct., 1997, 8(8), 711–718. 62 P. Ghosh and A. R. Srinivasa, A Two-network thermomechanical model of a shape memory polymer, Int. J. Eng. Sci., 2011, 49(9), 823–838.
This journal is © The Royal Society of Chemistry 2014
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
PAPER
View Article Online View Journal
Cite this: DOI: 10.1039/c4sm01816h
Temperature memory effect in amorphous shape memory polymers† Kai Yu and H. Jerry Qi* Temperature memory effect (TME) refers to the ability of shape memory polymers (SMPs) to memorize the temperature at which pre-deformation was conducted. In the past few years, this TME was experimentally demonstrated by comparing the applied programming temperature (Td) with a characteristic recovery temperature (Tc), which corresponds to either the maximum recovery stress or free recovery speed. In these well-designed experiments, Tc was observed to be close to Td, which is consistent with the intuitive understanding of ‘memorization’. However, since the polymer recovery behavior has been proved to be strongly dependent on various programming and recovery conditions, a new question that whether Tc is always equal to Td in any thermo-temporal conditions remains to be addressed. In this paper, we answered this question by examining the free recovery profile of an acrylate based amorphous SMP. The recovery Tc, which is the temperature with the maximum recovery speed, versus the recovery temperature is shown to be strongly dependent on both programming and recovery conditions. Their detailed influence could be explained by using the reduced time. During a thermomechanical working cycle of SMPs, in addition to the Td, any other thermo-temporal conditions, such as the holding time (th),
Received 16th August 2014 Accepted 18th September 2014
cooling rate, recovery heating rate (q), etc., can affect the observed Tc by changing the reduced programming or recovery time. In this manner, the relationship between Tc and Td is not uniquely determined. Besides, the TME in SMPs can only be achieved within a given temperature range. Both
DOI: 10.1039/c4sm01816h
onset and offset of this temperature range are shown to be influenced by the programming history, but are independent of the recovery conditions.
www.rsc.org/softmatter
1. Introduction The conventional denition of shape memory polymers (SMPs) is “a class of stimuli-responsive smart materials that can be deformed at high temperatures and subsequently be frozen upon cooling to x the programmed temporary shape, which remains stable unless an environmental stimulus is applied, such as temperature,1–8 magnetic elds,9–14 light,15–20 moisture,21,22 etc.”. Depending on whether or not the external load still exists when the SMP is activated, the polymer will demonstrate the shape memory (SM) effect by providing recovery stresses or showing shape recovery, namely a constrained or free recovery process. Among the variety of external stimuli that might be utilized as the triggering mechanism, the most common one is direct heating to increase the temperature. Although the basic concept of SMPs has been known for decades, recent advances have led to the discovery of previously uncovered memory properties that challenge our understanding of this type of polymer material.23–31 For example, the
The George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA. E-mail: [email protected] † Electronic supplementary 10.1039/c4sm01816h
information
(ESI)
available.
This journal is © The Royal Society of Chemistry 2014
See
DOI:
programming condition has been shown to signicantly affect the shape memory behavior of SMPs even if the same thermotemporal recovery condition is used.27,29,30,32–34 A temperature memory effect (TME), in this occasion, refers to the capability of SMPs to memorize the applied programming temperatures (Td). Previous studies in this area demonstrated that both free recovery curves and maximum recovery stress under strain constraint would be shied to a higher recovery temperature if a higher programming temperature was used.4,35 Subsequently, Miaudet et al.28 tested the shape memory property in a polyvinyl alcohol (PVA) based nanocomposite, where the temperature memory is described such that the temperature corresponding to the maximum recovery stress is roughly identical to Td. Xie et al.23,25,36 also reported the free recovery based TME in Naon polymer. The temperature at which a maximum shape recovery rate was observed was also shown to be close to the applied Td. From these pioneer experimental investigations, one can gain recognition that the shape memory characteristic temperature (Tc), which either corresponds to the maximum recovery stress or shape recovery speed, can be quantitatively related to the programming temperature. However, the question remains whether Tc would always be equal to Td in any SM cycle, or their detailed relationship is also dependent on other experimental parameters during the
Soft Matter
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Soft Matter
programming and recovery steps. Many previous studies have shown that the free recovery curves of SMPs are strongly dependent on the packaging programming conditions, such as the loading rate, holding time (th), cooling rate, etc.27,30,37–40 Tobushi et al.41 also experimentally demonstrated the inuence of the recovery heating rate (q) on the free recovery behavior of SMPs. In the case of a higher q being applied, the amount of recovered strain is indeed lower when a given recovery temperature is reached, namely a delayed shape recovery. In this manner, Tc is also dependent on the recovery heating rate. More recently, Wang et al.42 also showed that the free recovery based TME could only be observed within a specic range of recovery temperature, but not in its entire range of glass transition. These observations indicate that the relationship between Tc and Td is not so straightforward as we expected. The detailed programming and recovery conditions during the SM cycle, as well as the polymer temperature dependent thermoviscoelastic properties, may be used to determine Tc during the free recovery of SMPs, and should be considered when we study the TME. In this paper, we studied the TME of thermally responsive amorphous SMPs under the aforementioned considerations, and explored the detailed relationship between Tc and Td during the polymer free shape recovery. Experimental observations rstly revealed the inuence of the programming temperature and heating rate on free shape recovery behavior of an acrylate based SMP material. A one dimensional (1D) multi-branch rheological model was then used to capture the experimental observations, and parametrically predict the free recovery Tc under various programming and recovery conditions. We conclude that both the temperature range to observe the TME in the SMPs, and the relationship between Tc and Td, are not unique. They both are dependent on the thermo-temporal history within a shape memory (SM) cycle. The TME in the SMPs is indeed a synergistic effect of both programming and recovery conditions on the material free recovery behavior.
Paper
2. 2.1
Experimental Materials
The SMP material used in this study is an acrylate-based network polymer. Its liquid monomer tert-butyl acrylate (tBA), liquid cross-linker poly(ethylene glycol) dimethacrylate (PEGDMA) and the photo-polymerization initiator (2,2-dimethoxy-2-phenylacetophenone) in the powder form were ordered from Sigma Aldrich (St Louis, MO, USA) and used in their asreceived conditions without any further purication. The synthesis procedure of the SMP followed Yakacki et al.43 by mixing the tBA monomer, cross-linker and photo-polymerization initiator with a weight ratio of 4000 : 1000 : 1, and the solution was then exposed to UV irradiation for 10 min before the nal heat treatment at 90 C for 1 hour to assist the complete conversion of monomers. All the experiments involved in this paper were conducted on a Dynamic Mechanical Analyzer (DMA, TA Instruments, Model Q800, New Castle, DE, USA). The SMP was cast into rectangular specimens with a uniform dimension of 10 mm 3 mm 0.2 mm. 2.2
DMA test
DMA tests were performed to characterize the glass transition behavior of the acrylate SMPs. The SMP specimen was rstly heated to 95 C on the DMA machine and stabilized for 20 minutes to reach thermal equilibrium. During the experiment, the strain was oscillated at a frequency of 1 Hz with a peak-topeak amplitude of 0.1% while the temperature was decreased from 95 C to 5 C at a rate of 1 C min1. 2.3
Thermomechanical working cycle of SMPs
Free recovery tests of acrylate SMP were conducted on the DMA machine by using strain controlled methods, and the detailed procedure followed the SM cycle shown in Fig. 1a. In the programming step, the SMP was stretched to a target strain emax (20%) with a constant loading rate ˙e (0.01 s1) at the
Fig. 1 A schematic illustration of (a) the thermomechanical history of the programming and free recovery process in a SM cycle, and (b) the characteristic free recovery temperature Tc.
Soft Matter
This journal is © The Royal Society of Chemistry 2014
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Paper
Soft Matter
programming temperature Td, followed by a specied holding time (th) at Td before being cooled to the shape xing temperature TL (20 C) at a rate of 2.5 C min1. Once TL was reached, the specimen was held for one hour then the tensile force was removed. In the free recovery step, the temperature was increased to the recovery temperature Tr at a given heating rate of q. The strain evolution during this step was used to calculate the time dependent shape recovery ratio. Throughout the paper, the shape recovery ratio Rr and the shape recovery speed Vr are respectively dened as25,44 Rr ¼ 1 e(t)/er
(1a)
and Vr ¼
vRr ; vTH
(1b)
where e(t) is the strain evolution during the free recovery step, er is the strain measured at the start of the free recovery process, TH(t) ¼ qt is the time dependent recovery temperature, and the characteristic recovery temperature Tc is taken to be the temperature corresponding to the maximum recovery speed (see Fig. 1b).
3.
Constitutive modeling of SMPs
Currently, theoretical studies of shape memory phenomena in polymers are mainly based on two modeling approaches: phase evolution45–50 and thermoviscoelastic modeling method.27,51–55 Here in this study, to capture the shape memory properties of the acrylate-based SMP, we use a multi-branched thermoviscoelastic constitutive model. The model resembles a generalized Maxwell element that consists of one equilibrium branch and several thermoviscoelastic nonequilibrium branches (number of n) to represent the multiple relaxation processes (for example, Rouse modes) of the polymer. Eeq and Ei are the elastic moduli in the equilibrium and nonequilibrium branches, respectively. si(T) denotes the temperature dependent relaxation time in the dashpots. A detailed description of the multibranch model for SMP application can be found in the study by
Yu et al.27 which is also shown in the ESI† of this paper. Its extension for a 3D nonlinear model can be found in Westbrook et al.54 and Diani et al.53 Here, it is assumed that the time–temperature shi for each branch follows the same rule. According to the well-established “thermorheological simplicity” principle56 under a nonisothermal condition, the relaxation times (or viscosity) of each nonequilibrium branch vary as the temperature changes: si ¼ aT(T)s0i ,
(2)
where aT(T) is a time–temperature superposition (TTSP) shiing factor, and s0i is the reference relaxation time at the temperature when aT(T) ¼ 1. Depending on whether the temperature is above or below a shiing temperature Ts,57 the shi factor is respectively calculated by the following two equations, namely the Williams– Landel–Ferry (WLF) equation58 and the Arrhenius-type equation:59 C1 ðT Tr Þ log aT T ¼ ; C2 þ ðT Tr Þ
(3a)
AFc 1 1 ln aT T ¼ : Tg kb T
(3b)
where Tr is the WLF reference temperature, Tg is the glass transition temperature of SMP, C1 and C2 are material constants that depend on the choice of Tr, A is a material constant, Fc is the congurational energy and kb is Boltzmann's constant. Ts is taken to be the crossing point of two curves representing eqn (3a) and (3b) on an aT(T) vs. T plot. The model parameters were determined by using tan d and storage modulus of the SMP obtained from the DMA tests. The identication procedure is briey introduced in the ESI.† A detailed model parameters identication method of the multibranch model can be found in previous studies.27,52,53 It should be noted that in this paper, the applied modeling method is based on linear thermoviscoelasticity. Caution should be taken if a SMP demonstrates highly nonlinear or inelastic behavior. The small strain DMA test will be insufficient to characterize the
(a) 1D estimation for the DMA result of a SMP by using NLREG software. The original experimental curves are shown in solid lines. (b) Shifting factors at different temperatures: above 27.5 C, shifting factors follow the WLF equation; below 27.5 C, shifting factors follow the Arrhenius-type behavior. Fig. 2
This journal is © The Royal Society of Chemistry 2014
Soft Matter
View Article Online
Soft Matter
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Table 1
Paper
Parameters of the multi-branched constitutive model
Parameters
Values
C1 AFckb1 Ts Eeq Ei
16.44 C2 51.6 C 40 000 K Tg 42.7 C 27.5 C TM 30 C 3.05 MPa (424.2, 424.2, 399.2, 276.6, 329.7, 245.3, 153.4, 64.4, 23.5, 7.51, 3.02, 1.50, 0.71, 0.43, 0.12, 0.16) MPa for i ˛ [1,16] (100.5, 100, 100.5, 101, 101.5, 102, 102.5, 103, 103.5, 104, 104.5, 105, 105.5, 106, 106.5, 107) s, for i ˛ [1,16]
s0i
Parameters
Values
thermoviscoelastic (or thermoviscoplastic) behavior and thus may not be able to provide an accurate estimation of the SM performance. In this case, a nonlinear modeling approach such as those presented by Westbrook et al.,54 Abrahamson et al.,60 Tobushi et al.61 and Srinivasa et al.62 can be used. As shown in Fig. 2a, the identied parameters could enable the multi-branch model to successfully capture the experimental storage modulus and tan d curves within the entire testing temperature range (5–95 C). The temperature dependent shi factors were then calculated by using eqn (3) and are plotted in Fig. 2b as a function of temperature. The nal set of model parameters is listed in Table 1. Compared with our previous work, where seven nonequilibrium branches are constructed in the multi-branch model, the amount of nonequilibrium branches in this study is extended to sixteen with rened relaxation time distribution, which will enable a smooth and accurate prediction on the recovery speed curves versus the recovery temperature.
4. Results and discussion 4.1
Experimental observations
We start our discussion by rstly showing some experimental observations on the free recovery ratio and speed under different programming and recovery conditions. As a baseline SM cycle, the programming temperature Td, the holding time th at Td and the recovery heating rate q are respectively set to be 35 C, 60 s and 1 C min1. The experimental setup is the same as
that described in Section 2.3. As a comparison, two more SM cycles were performed: one with an increased programming temperature (Td ¼ 40 C), while the other one with a decreased heating rate (q ¼ 0.1 C min1). The shape recovery ratio and speed in these three SM cycles are respectively plotted in Fig. 3a and b as a function of recovery temperature (TH). As shown in Fig. 3, both the programming temperature and the heating rate signicantly affect the material free recovery prole: on one hand, increasing Td from 35 C to 40 C decreases the shape recovery speed under the same heating rate, and elevates Tc from 34 C to 38 C. On the other hand, similar to the previous study by Tobushi et al.,41 decreasing q increases the recovery speed versus the recovery temperature, and consequently shis Tc to a small value. This is because a lower heating rate, which means a longer heating time before reaching a given recovery temperature, would allow the SMPs to evolve more towards the equilibrium state and hence give a larger recovery strain. When the heating rate is set to be 0.1 C min1, Tc is seen to decrease from 34 C to 30 C. In addition to the experimental observations, simulation results based on the multi-branch model and predetermined model parameters are also plotted in Fig. 3 as dashed lines. It can be seen that the numerical simulations adequately capture the temperature dependent shape recovery ratio and normalized recovery speed. In the following, parametric studies based on the multi-branch model will be performed to investigate the inuence of programming and recovery conditions on the TME in SMPs. 4.2
Inuence of the programming condition
The inuence of Td on the polymer free recovery prole is predicted and demonstrated in Fig. 4a and b, where Td is gradually increased from 25 C to 59 C (2 C temperature interval) during the simulation. The holding time th at Td and the heating rate during the free recovery are respectively set to be 60 s and 1 C min1 in each case. Fig. 4c further plots the Tc during the polymer free recovery as a function of applied Td. Two main features could be found from the numerical simulations in Fig. 4. Firstly, within a specic temperature range (30 C to 47.5 C), increasing the programming temperature will decrease the recovery speed versus the recovery
Fig. 3 Influence of the programming temperature and heating rate on the (a) shape recovery ratio and (b) normalized shape recovery speed. Note: the dashed lines denote simulation curves by using the 1D multi-branch model.
Soft Matter
This journal is © The Royal Society of Chemistry 2014
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Paper
Soft Matter
Fig. 4 Parametric study on the influence of the programming temperature. Td increases from 25 C to 59 C with a 2 C temperature interval. (a) Shape recovery ratio as a function of recovery temperature. (b) Shape recovery speed as a function of recovery temperature. (c) Tc as a function of Td.
temperature. Tc increases with Td in exhibiting a roughly linear relationship. This obeys the recognition of the traditional TME that Tc observed during the free recovery step is close to the Td applied in the programming step. Secondly, beyond this temperature range, changing Td does not obviously affect Tc. This prediction, that the TME could only be realized within a given temperature range, is consistent with the recent study of Wang et al.42 The two features shown in Fig. 4 could be explained by examining the stress distribution within each branch of the multi-branch model. Previously, we demonstrated that if the same thermo-temporal recovery condition is applied on SMPs, their free recovery proles are only dependent on the stress distribution aer unloading, or namely the stress distribution at the beginning of recovery.27 Fig. 5a shows this stress distribution at the end of seven different programming conditions with different Td values (27.5 C, 30 C, 40 C, 47.5 C and 60 C respectively). At a temperature higher than 47.5 C, Tc is observed to be independent of Td. This phenomenon can be explained as follows: when the SMPs are stretched at higher temperatures, the elastic strain development within each nonequilibrium branch is small due to the high mobility of dashpots. The stress relaxation in the subsequent holding and cooling steps will further release these nonequilibrium elastic strains. During the unloading at low temperatures, removing the external load will cause a re-distribution of strains (due to the rebalance of forces
This journal is © The Royal Society of Chemistry 2014
to maintain overall equilibrium) in all the branches,24,27 and an internal stress state aer unloading is that only the equilibrium spring remains in tension while the rest springs are being compressed (see Fig. 5a for the cases of Td ¼ 47.5 C and 60 C). In this manner, when the temperature is raised again during the recovery step, only the equilibrium branch provides driven force for viscous strain release (see Fig. 5b) and material shape recovery, while the rest branches keep on resisting until the contained dashpots are activated. Since the nonequilibrium elastic stress is already released essentially before unloading for the case of Td ¼ 47.5 C, further increasing the programming temperature to 60 C will not affect the rheological stress state, and therefore, the subsequent free recovery curve of SMPs is almost the same. This is the reason for Td ¼ 47.5 C being observed as the upper temperature limit for the TME in the acrylate SMP as seen in Fig. 4c. However, when the programming temperature is gradually decreased from 47.5 C to 30 C, mobility of the dashpots decreases, and considerable elastic strain could be developed when stretching the multi-branch model, especially in branches with larger relaxation times. These elastic strains cannot be totally released during the subsequent holding and cooling, and will be stable when the temperature is low. As seen in Fig. 5a, decreasing the programming temperature will change the rheological stress state aer unloading in two ways: rstly, the stress magnitude in each nonequilibrium branch is increased. The intensied stress state will proportionally increase the
Soft Matter
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Soft Matter
Paper
(a) Stress distribution after unloading in seven different programming conditions. (b) Viscous strain release in each nonequilibrium branch and strain evolution in the equilibrium branch for the case of Td ¼ 47.5 C and Td ¼ 30 C.
Fig. 5
releasing rate of viscous strain in the dashpots (see the bottom gure in Fig. 5b). Secondly, aer the stress re-distribution among all the branches during the unloading, more and more nonequilibrium springs will remain stretched in addition to the equilibrium spring. These branches will provide additional driven forces for the shape recovery of the multi-branch model, and consequently increase the recovery speed. In the rheological multi-branch model, this additional driven force, namely the tensile stress in the nonequilibrium branches at the beginning of recovery, essentially represents the internal energy stored in the real polymer system aer the programming step. When stretching the polymer network at lower temperatures, the mobility of polymer chains might be insufficient for them to adjust the self-conformation and fully comply with the external deformation. In this manner, in addition to the network entropy change, some internal energy will also be generated and stored within the SMPs. This internal energy originates from the alteration in the macromolecule bond length and angle during the deformation, and will assist the spontaneous shape recovery of the polymer network upon reheating. With regard to the inuence of the programming temperature as revealed in Fig. 5a, it seems that decreasing the programming temperature is a good choice to increase the shape recovery speed. But the material's shape xing ability will be compromised because the material tends to be unstable during unloading and a large portion of the programmed deformation is recovered right aer unloading,27 which is unfavorable in most engineering applications. For example, when Td ¼ 30 C, only 41% of the programmed strain is xed aer unloading. Further decreasing the programming temperature will continually decrease the shape xity of the polymer, but does not affect the stress state too much. In this manner, the recovery Tc will remain unchanged when Td is below 30 C.
Soft Matter
Now it is clear that the rheological stress state at the beginning of free recovery is the key to determine Tc. Previously, we have demonstrated that this stress state is only dependent on the reduced programming time tr.27 For a non-isothermal process, ðt ds ; (4) tr ¼ aðTÞ 0 where t is the physical time and T is the temperature evolution during the programming step. Therefore, in addition to Td, any other thermo-temporal programming conditions, such as holding time th aer loading, cooling rates, shape xing temperature TL, etc., can also affect the recovery speed and the Tc by changing the reduced programming time and then consequently the rheological stress state. Here, we select th as the parametric variable, and the predicted free recovery curves are shown in Fig. 6a. Td and the heating rate q are set to be 35 C and 1 C min1 respectively in each simulation case. As shown in the gure, increasing the holding time th will decrease the recovery speed because it plays the same role as the increasing Td in determining the reduced programming time. Besides, two boundary holding times are revealed from the gure: one represents the extreme condition that no holding time is applied (black curve). While for the second one, it indicates that a 6 104 s (17 hours) holding time at Td ¼ 35 C is long enough to essentially release the elastic stresses in nonequilibrium branches during the programming, and the corresponding recovery prole is therefore the same as that with th ¼ 6 108 s (19 years). Subsequently, Fig. 6b summarizes the recovery Tc as a function of both Td and the holding time th. Clearly, the relationship between Tc and Td is not unique under the inuence of th. A longer holding time will increase the Tc towards a saturated level at 47.5 C. Besides, both the onset and offset values of the
This journal is © The Royal Society of Chemistry 2014
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Paper
Soft Matter
Fig. 6 Influence of both Td and holding time on the free recovery profiles and Tc of SMPs. (a) The predicted free recovery curves with different holding times. Td ¼ 35 C. (b) Tc plotted as a function of both holding time and Td. (c) The constructed master curve of Tc. The reference holding time th0 ¼ 6 101 s. (d) The holding time dependent shift factors during the construction of the master curve.
temperature range with the TME are observed to decrease with increment of the holding time. For the curves in Fig. 6b, aer shiing each of them horizontally to a reference holding time th0, namely, by Tc(Td,th) ¼ Tc(athTd,th0),
(5a)
where ath is the holding time dependent shi factor, a master curve can be constructed as shown in Fig. 6c. The corresponding shi factors are plotted in Fig. 6d as a function of holding time. A linear relationship between ath and th is then revealed in the semi-log scale. By using the nonlinear regression method, we approximate the relationship as th ath ¼ 0:85 þ 0:083 log (5b) th0 Within the temperature range with a clear TME observed, eqn (5) quantitatively demonstrates the packaging effect of both the programming temperature and the holding time on Tc. That is, when a large holding time is applied during the programming step, the corresponding shi factor will be exponentially ramped and therefore a smaller Td should be used if an equivalent Tc is expected. Otherwise Tc will monotonically increase. The analysis on the effect of increasing Td and th can be extended to the inuence of other thermo-temporal
This journal is © The Royal Society of Chemistry 2014
programming conditions, such as decreasing loading rate, cooling rates, or increasing the stabilization time at TL, etc. 4.3
Inuence of recovery conditions
In addition to the programming condition, we will show in this section that the relationship between Tc and Td is also dependent on the free recovery condition. Considering the SM cycle shown in Fig. 1a, the only parametric variable under investigation is the recovery heating rate q. Fig. 7a shows the predicted free recovery curves with q increasing from 103 C min1 to 102 C min1. All the programming procedures are set to be identical in each simulation case with Td ¼ 35 C and a holding time equal to 60 s. As demonstrated in our previous study,27 if the same programming condition is used, the free recovery ratio of SMPs is only dependent on the reduced recovery time, and a longer reduced recovery time will lead to a higher recovery ratio, namely Rr ¼ Rr(tr)
(6a)
with tr ¼
ð TH TL
dT ; qaðTÞ
(6b)
where TH is a given recovery temperature. Eqn (6b) indicates that a larger heating rate will decrease the reduced recovery
Soft Matter
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Soft Matter
Paper
Fig. 7 Influence of both Td and heating rate q on the free recovery profiles and Tc of SMPs. (a) The predicted free recovery curves with different heating rates. Td ¼ 35 C and holding time equals to 6 s. (b) Tc plotted as a function of both Td and q. (c) The constructed master curve of Tc. The reference heating rate q0 ¼ 100 C min1. (d) The heating rate dependent shift factors during the construction of the master curve.
time and then consequently the recovery ratio when TH is reached, as conrmed in Fig. 7a. Fig. 7b summarizes the recovery Tc as a function of both Td and q. In addition to the revealed relationship between Tc and Td that is dependent on q, the temperature range to observe the TME is shown to be unchanged upon changing the heating rate. Subsequently, by shiing each curve in Fig. 7b vertically to a reference heating rate (q0 ¼ 100 C min1) according to Tc ðTd ; qÞ ¼
1 Tc ðTd ; q0 Þ aq
(7a)
where aq is a heating rate dependent shi factor, the curves are seen to be rejoined together to form a master curve, as shown in Fig. 7c. The corresponding shi factors are summarized in Fig. 7d. In the semi-log scale, the shi factor decreases linearly with increment of the heating rate, namely q aq ¼ 1 0:092 log (7b) q0 Eqn (7) quantitatively approximates the inuence of the recovery heating rate on Tc with two features revealed: rstly, when a larger heating rate q is applied during the polymer free recovery step, both the shi factor and the corresponding Tc will decrease. Secondly, different from the above-mentioned effect of programming conditions, there is no boundary heating rate condition that inuences the heating rate. Continuous increase of q will continually elevate Tc towards the target TH. An
Soft Matter
innitely large value represents an isothermal recovery condition, where Tc is always equal to TH regardless of the applied programming temperature. In contrast, if no heat is applied during the recovery step, namely q ¼ 0, a thermally simple SMP can still recover to its original shape, even though this process may take several years.8,27 Under this condition, Tc is equal to the shape xing temperature TL.
5.
Conclusion
In this paper, from both the experimental and theoretical point of view, we investigated the temperature memory effect (TME) in the acrylate based amorphous shape memory polymer (SMP). Different from the traditional denition of TME in the SMPs, our study indicates that the characteristic recovery temperature (Tc), which is the temperature with the maximum recovery speed, is not always equal to the programming temperature (Td). By revealing the underlying mechanism, we can gain a deeper understanding of how the programming and recovery conditions will affect the free recovery behavior of SMPs synergistically, which consequently provides a facile strategy to optimize the shape memory performance of SMPs by adjusting their thermomechanical working conditions, instead of developing new polymer systems. Such a strategy offers huge convenience for engineering SMPs in the industry, especially in the elds of medicine, aeronautics and astronautics, where a long period of performance verication and official approval is typically expected.
This journal is © The Royal Society of Chemistry 2014
View Article Online
Paper
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
References 1 A. Lendlein and S. Kelch, Shape-memory polymers, Angew. Chem., Int. Ed. Engl., 2002, 41(12), 2035–2057. 2 C. Liu, H. Qin and P. T. Mather, Review of progress in shapememory polymers, J. Mater. Chem., 2007, 17(16), 1543–1558. 3 T. Chung, A. Rorno-Uribe and P. T. Mather, Two-way reversible shape memory in a semicrystalline network, Macromolecules, 2008, 41(1), 184–192. 4 K. Gall, et al., Thermomechanics of the shape memory effect in polymers for biomedical applications, J. Biomed. Mater. Res., Part A, 2005, 73(3), 339–348. 5 P. T. Mather, X. F. Luo and I. A. Rousseau, Shape Memory Polymer Research, Annu. Rev. Mater. Res., 2009, 39, 445–471. 6 K. Yu, et al., Carbon nanotube chains in a shape memory polymer/carbon black composite: To signicantly reduce the electrical resistivity, Appl. Phys. Lett., 2011, 98(7), 0741021–0741023. 7 K. Yu, Y. J. Liu and J. S. Leng, Conductive Shape Memory Polymer Composite Incorporated with Hybrid Fillers: Electrical, Mechanical, and Shape Memory Properties, J. Intell. Mater. Syst. Struct., 2011, 22(4), 369–379. 8 Q. Ge, et al., Prediction of temperature-dependent free recovery behaviors of amorphous shape memory polymers, So Matter, 2012, 8(43), 11098–11105. 9 R. Mohr, et al., Initiation of shape-memory effect by inductive heating of magnetic nanoparticles in thermoplastic polymers, Proc. Natl. Acad. Sci. U. S. A., 2006, 103(10), 3540–3545. 10 P. R. Buckley, et al., Inductively heated shape memory polymer for the magnetic actuation of medical devices, IEEE Trans. Biomed. Eng., 2006, 53(10), 2075–2083. 11 A. M. Schmidt, Electromagnetic activation of shape memory polymer networks containing magnetic nanoparticles, Macromol. Rapid Commun., 2006, 27(14), 1168–1172. 12 Z. W. He, et al., Remote Controlled Multishape Polymer Nanocomposites with Selective Radiofrequency Actuations, Adv. Mater., 2011, 23(28), 3192–3196. 13 U. N. Kumar, et al., Non-contact actuation of triple-shape effect in multiphase polymer network nanocomposites in alternating magnetic eld, J. Mater. Chem., 2010, 20(17), 3404–3415. 14 K. Yu, et al., Design considerations for shape memory polymer composites with magnetic particles, Adv. Mater., 2013, 47(1), 51–63. 15 H. Y. Jiang, S. Kelch and A. Lendlein, Polymers move in response to light, Adv. Mater., 2006, 18(11), 1471–1475. 16 H. Koerner, et al., Remotely actuated polymer nanocomposites – stress-recovery of carbon-nanotube-lled thermoplastic elastomers, Nat. Mater., 2004, 3(2), 115–120. 17 A. Lendlein, et al., Light-induced shape-memory polymers, Nature, 2005, 434(7035), 879–882. 18 M. H. Li, et al., Light-driven side-on nematic elastomer actuators, Adv. Mater., 2003, 15(7–8), 569–572.
This journal is © The Royal Society of Chemistry 2014
Soft Matter
19 T. F. Scott, R. B. Draughon and C. N. Bowman, Actuation in crosslinked polymers via photoinduced stress relaxation, Adv. Mater., 2006, 18(16), 2128–2132. 20 T. F. Scott, et al., Photoinduced plasticity in cross-linked polymers, Science, 2005, 308(5728), 1615–1617. 21 W. M. Huang, et al., Water-driven programmable polyurethane shape memory polymer: Demonstration and mechanism, Appl. Phys. Lett., 2005, 86(11), 1141051– 1141053. 22 Y. C. Jung, H. H. So and J. W. Cho, Water-responsive shape memory polyurethane block copolymer modied with polyhedral oligomeric silsesquioxane, J. Macromol. Sci., Part B: Phys., 2006, 45(4), 453–461. 23 T. Xie, Tunable polymer multi-shape memory effect, Nature, 2010, 464(7286), 267–270. 24 K. Yu, et al., Mechanisms of multi-shape memory effects and associated energy release in shape memory polymers, So Matter, 2012, 8(20), 5687–5695. 25 T. Xie, K. A. Page and S. A. Eastman, Strain-Based Temperature Memory Effect for Naon and Its Molecular Origins, Adv. Funct. Mater., 2011, 21(11), 2057–2066. 26 T. Xie, X. C. Xiao and Y. T. Cheng, Revealing Triple-Shape Memory Effect by Polymer Bilayers, Macromol. Rapid Commun., 2009, 30(21), 1823–1827. 27 K. Yu, G. Qi and H. J. Qi, Reduced time as a unied parameter determining xity and free recovery of shape memory polymers, Nat. Commun., 2014, 5, 3066. 28 P. Miaudet, et al., Shape and temperature memory of nanocomposites with broadened glass transition, Science, 2007, 318(5854), 1294–1296. 29 X. Chen and T. D. Nguyen, Inuence of thermoviscoelastic properties and loading conditions on the recovery performance of shape memory polymers, Mech. Mater., 2011, 43, 127–138. 30 F. Castro, et al., Time and Temperature Dependent Recovery of Epoxy-Based Shape Memory Polymers, J. Eng. Mater. Technol., 2011, 133(2), 0210251–0210259. 31 L. Sun and W. M. Huang, Mechanisms of the multi-shape memory effect and temperature memory effect in shape memory polymers, So Matter, 2010, 6(18), 4403–4406. 32 J. L. Hu, F. L. Ji and Y. W. Wong, Dependency of the shape memory properties of a polyurethane upon thermomechanical cyclic conditions, Polym. Int., 2005, 54(3), 600–605. 33 J. Cui, K. Kratz and A. Lendlein, Adjusting shape-memory properties of amorphous polyether urethanes and radioopaque composites thereof by variation of physical parameters during programming, Smart Mater. Struct., 2010, 19(6), 0650191–0650199. 34 M. Bothe, F. Emmerling and T. Pretsch, Poly(ester urethane) with Varying Polyester Chain Length: Polymorphism and Shape-Memory Behavior, Macromol. Chem. Phys., 2013, 214(23), 2683–2693. 35 Y. Miyamoto, et al., Memory effect on the glass transition in vulcanized rubber, Phys. Rev. Lett., 2002, 88(22), 2255041– 2255044.
Soft Matter
View Article Online
Published on 18 September 2014. Downloaded by Northeastern University on 31/10/2014 13:31:49.
Soft Matter
36 J. J. Li and T. Xie, Signicant Impact of Thermo-Mechanical Conditions on Polymer Triple-Shape Memory Effect, Macromolecules, 2011, 44(1), 175–180. 37 X. Chen and T. D. Nguyen, Inuence of thermoviscoelastic properties and loading conditions on the recovery performance of shape memory polymers, Mech. Mater., 2011, 43(3), 127–138. 38 A. McClung, G. Tandon and J. Baur, Strain rate- and temperature-dependent tensile properties of an epoxybased, thermosetting, shape memory polymer (Veriex-E), Mech. Time-Depend. Mater., 2011, 16(2), 205–221. 39 A. McClung, G. Tandon and J. Baur, Deformation rate-, hold time-, and cycle-dependent shape-memory performance of Veriex-E resin, Mech. Time-Depend. Mater., 2013, 17(1), 39–52. 40 A. J. W. McClung and M. B. Ruggles-Wrenn, Strain Rate Dependence and Short-Term Relaxation Behavior of a Thermoset Polymer at Elevated Temperature: Experiment and Modeling, J. Pressure Vessel Technol., 2009, 131(3), 0314051–0314058. 41 H. Tobushi, et al., Thermomechanical constitutive model of shape memory polymer, Mech. Mater., 2001, 33(10), 545–554. 42 Y. R. Wang, et al., Relation between temperature memory effect and multiple-shape memory behaviors based on polymer networks, RSC Adv., 2014, 4(39), 20364–20370. 43 C. M. Yakacki, et al., Unconstrained recovery characterization of shape-memory polymer networks for cardiovascular applications, Biomaterials, 2007, 28(14), 2255–2263. 44 A. Lendlein and R. Langer, Biodegradable, elastic shapememory polymers for potential biomedical applications, Science, 2002, 296(5573), 1673–1676. 45 Y. P. Liu, et al., Thermomechanics of shape memory polymers: Uniaxial experiments and constitutive modeling, Int. J. Plast., 2006, 22(2), 279–313. 46 W. W. Muller and T. Pretsch, Hydrolytic aging of crystallizable shape memory poly(ester urethane): Effects on the thermo-mechanical properties and visco-elastic modeling, Eur. Polym. J., 2010, 46(8), 1745–1758. 47 V. Kaa, Shape memory polymers: A mesoscale model of the internal mechanism leading to the SM phenomena, Int. J. Plast., 2008, 24(9), 1533–1548. 48 Y. C. Chen and D. C. Lagoudas, A constitutive theory for shape memory polymers. Part I-Large deformations, J. Mech. Phys. Solids, 2008, 56(5), 1752–1765.
Soft Matter
Paper
49 K. N. Long, M. L. Dunn and H. J. Qi, Mechanics of so active materials with phase evolution, Int. J. Plast., 2010, 26(4), 603– 616. 50 H. J. Qi, et al., Finite deformation thermo-mechanical behavior of thermally induced shape memory polymers, J. Mech. Phys. Solids, 2008, 56(5), 1730–1751. 51 J. R. Lin and L. W. Chen, Shape-memorized crosslinked ester-type polyurethane and its mechanical viscoelastic model, J. Appl. Polym. Sci., 1999, 73(7), 1305–1319. 52 K. Yu, et al., A thermomechanical constitutive model for an epoxy based shape memory polymer and its parameter identications, Mech. Time-Depend. Mater., 2014, 18(2), 453–474. 53 J. Diani, et al., Predicting thermal shape memory of crosslinked polymer networks from linear viscoelasticity, Int. J. Solids Struct., 2012, 49(5), 793–799. 54 K. K. Westbrook, et al., A 3D Finite Deformation Constitutive Model for Amorphous Shape Memory Polymers: A MultiBranch Modeling Approach for Nonequilibrium Relaxation Processes, Mech. Mater., 2011, 43(12), 853–869. 55 T. D. Nguyen, et al., A thermoviscoelastic model for amorphous shape memory polymers: Incorporating structural and stress relaxation, J. Mech. Phys. Solids, 2008, 56(9), 2792–2814. 56 M. Rubinstein and R. H. Colby, Polymer Physics, Oxford University Press, Oxford, New York, 2003. 57 P. A. O'Connell and G. B. McKenna, Arrhenius-type temperature dependence of the segmental relaxation below T-g, J. Chem. Phys., 1999, 110(22), 11054–11060. 58 M. L. Williams, R. F. Landel and J. D. Ferry, Temperature Dependence of Relaxation Mechanisms in Amorphous Polymers and Other Glass-Forming Liquids, Phys. Rev., 1955, 98(5), 1549. 59 E. A. Di Marzio and A. J. M. Yang, Congurational entropy approach to the kinetics of glasses, J. Res. Natl. Inst. Stand. Technol., 1997, 102(2), 135–157. 60 E. R. Abrahamson, et al., Shape memory mechanics of an elastic memory composite resin, J. Intell. Mater. Syst. Struct., 2003, 14(10), 623–632. 61 H. Tobushi, et al., Thermomechanical constitutive modeling in shape memory polymer of polyurethane series, J. Intell. Mater. Syst. Struct., 1997, 8(8), 711–718. 62 P. Ghosh and A. R. Srinivasa, A Two-network thermomechanical model of a shape memory polymer, Int. J. Eng. Sci., 2011, 49(9), 823–838.
This journal is © The Royal Society of Chemistry 2014
