Shape Memory Polyurethane Materials Containing Ferromagnetic Iron Oxide and Graphene Nanoplatelets

Intelligent materials, such as memory shape polymers, have attracted considerable attention due to wide range of possible applications. Currently, intensive research is underway, in matters of obtaining memory shape materials that can be actuated via inductive methods, for example with help of magnetic field. In this work, an attempt was made to develop a new polymer composite—polyurethane modified with graphene nanoplates and ferromagnetic iron oxides—with improved mechanical properties and introduced magnetic and memory shape properties. Based on the conducted literature review, gathered data were compared to the results of similar materials. Obtained materials were tested for their thermal, rheological, mechanical and shape memory properties. Structure of both fillers and composites were also analyzed using various spectroscopic methods. The addition of fillers to the polyurethane matrix improved the mechanical and shape memory properties, without having a noticeable impact on thermal properties. As it was expected, the high content of fillers caused a significant change in viscosity of filled prepolymers (during the synthesis stage). Each of the studied composites showed better mechanical properties than the unmodified polyurethanes. The addition of magnetic particles introduced additional properties to the composite, which could significantly expand the functionality of the materials developed in this work.


Introduction
Memory shape materials are the center of growing attention as they offer a way to solve many problems in such a wide array of fields as aeronautics, medicine, transportation, electronics, design and various everyday products [1]. It is a group of intelligent materials that, when influenced by an external stimuli-heat, light, humidity, pH level, electric or magnetic field-can alter some of their physical parameters, usually shape, but also color or stiffness. Systems based on memory shape properties, unlike other intelligent systems, require no additional actuators. They react on their own to the change in their environment. Currently, the most researched group of materials in this category are polymers (SMP), as shape memory effect can be easily obtained-the only requirement is the presence of phase separation, such as in polyurethanes, acrylane copolymers and styrene based copolymers [1].
SMPs have many advantages, to name few: material properties can be fine-tuned to requirements via changes in materials, additives and synthesis method; SMP can be sensitive to different stimuli and with an adequate process the multiple shape memory effect can be achieved; high repeated deformation; low density; and possibility to modify the material by introducing fillers. Of course, they possess some disadvantages as well: their low stiffness that can influence the ability to maintain

Modification of Fe 3 O 4 Particles
Iron oxide particles were modified with KH-550 in accordance to the procedure described by Cai et al. [11]. In short, 15 g of Fe 3 O 4 were dispersed in 500 mL of water-ethanol solution (1:1). Afterwards 8 mL of the compounding agent was added to the mixture and mixed for 4 h. The solution was left overnight for the sedimentation to occur. Next, the percipate was filtrated under reduced pressure and washed with the water-ethanol solution until neutral pH was obtained. Final participate [modified Fe 3 O 4 (mFe 3 O 4 )] was dried under vacuum at 60 • C for 24 h.

Preparation of PU/Fe 3 O 4 /GNP Composites
For the prepolymer preparation, melted PTMG and MDI were mixed at 65 • C for 30 min and then secondly 85 • C for 2 h. Next, composites were synthesized by taking the prepolymer and mixing it with fillers for 60 s at high speed. Then, adequate amounts of BDO and DABCO were added. After mixing for 10 s, the mixture was poured out into hot forms which were consequently placed inside an oven for 24 h at 80 • C and then further incubated for 48 h at 60 • C. The overall ratio of PTMG/MDI/BDO was 1/4/1. Following Cai et al. [11], the Fe 3 O 4 content was constant at 30 wt %, while the GNP content varied between 0 wt % and 1.2 wt %. The exact composition of the materials and their marking can be found in Table 1. The average thickness of materials was ca. 4 mm.

Characterization
Scanning Electron Microscopy (SEM) measurements of the fillers were conducted using the Quanta FEG 250 microscope (FEI Company, Hillsboro, OR, USA) with the ETD detector with 500-50,000× magnification at energy of 10 keV HV. Fourier transformation infrared spectra (FT-IR) were obtained using the Nicolet 8700 Spectrophotometer (Thermo Scientific, Madison, WI, USA with the measurement range of 500-4500 cm −1 ). For wide-angle X-ray scattering (XRD), Bragg-Brentano X'PERT PHILIPS diffractometer (PANalitical, Almelo, The Netherlands) was used to gather the XRD diffractograms (40 kV, 30 mA, λ Cu Kα = 0.1542 nm). Samples were scanned in 2θ range 5-80 • . Differential Scanning Calorimetry (DSC) has been carried by DSC209F1 (NETZSCH, Selb, Germany) instrument in the temperature range from −85 • C to 250 • C. The samples weight was about 10 mg. Analysis was conducted in nitrogen atmosphere with flow 40 mL/min. The samples were pre-heated (1st-heating) to eliminate thermal history, then cooling and 2nd-heating at 20 • C/min. Dynamical Mechanical Analysis (DMA) has been carried by DMA Q800 analyzer (TA Instruments, New Castle, Materials 2017, 10, 1083 4 of 23 DE, USA. Temperature investigation mode has been performed in bending mode (single cantilever), in temperature range from −100 to 100 • C, 1 Hz frequency amplitude and heating rate 4 • C/min using liquid nitrogen as cooling medium. Thermogravimetry (TGA) of obtained materials was conducted by TG 209 F3 analyzer (NETZSCH, Selb, Germany). Degradation process has been performed in the temperature range from 35 • C to 600 • C, nitrogen atmosphere with flow rate 40 mL/min and heating rate 20 • C/min. The rheological properties of the filled prepolymers were investigated using rheometer R/S-CPS+ apparatus (AMETEK Brookfield, Middleboro, MA, USA) with the controlled shear rate (CSR) 100 or 300 s −1 at 70 • C. The mechanical properties were measured using the Static Testing Machine Z020 (Zwick/Roell, Ulm, Germany), with the stretching speed of 50 mm/min.

Shape Memory Measurements
The shape memory properties of the composites were measured using the common fold-deploy method [20], however because of poor results several changes were made trying to find the best programming process for the materials. This procedure contains three steps: (1) rectangular samples (90 × 10 × 1 mm) were heated to programming temperature (T p ) in water bath and bent to a U-shape using a rod with 5 mm diameter to secure temporary shape; (2) samples were quickly placed into the ice-water bath for 5 min; and (3) samples were "released" and left for 5 min at rest temperature (T l )-the initial angle of the samples was recorded (θ f ). Next, the samples were placed in the water bath at temperature of 45 • C and the change of shape was observed with time (t) meanwhile final angle of the sample (θ r ) was being recorded. The procedure was repeated 3 times for each material, the temperature values were kept at constant value ±2 • C. The temporary shape retention rate (R f ) and the shape recovery rate (R r ) were calculated using the following formulas: Due to low melting temperature of soft segments, partial shape recovery occurred when the material was left at the room temperature. This prompted the placement of samples at the lower temperature (+4 • C) and measurement of the initial change of shape (∆ 0 ) while the material was "resting": Thus far, to modify SMP carbon nanotubes (CNT, MWCNT) [21], graphene oxide [22], its reduced form [23] and their modifications were used-with only few articles referring solely to GNP. This is why one of the objectives of the research was to understand the influence of the temperature on such materials. Initial results could hardly be called promising, therefore more modifications to the programming were introduced (see Table 2). However, the basic procedure was left unchanged. To further investigate the influence of temperature on the sample, the heating of the material was performed at three temperatures (T p ): 60, 80 and 90 • C. In other words, the shape memory investigation was performed under four different sets of conditions.

Results and Discussion
Fillers Morphology Figure 1 shows the SEM analysis results of the used fillers. The Fe 3 O 4 particles have average diameter of ca. 45 µm, but they can create lumps due to temporal auto-orientation of the magnetic domains. GNP are seen as thin flakes with diameters between 2 and 10 nm.

Spectroscopic Measurements
To analyze the structure of used fillers and obtained nanocomposites, FT-IR measurements were performed. Results are presented below on Figures 2 and 3. Lack of any signals on the GNP spectra points towards the lack of functional groups, which can be related to results obtained by Geng et al. [24]. Fe3O4 and mFe3O4 particles were dried in vacuum directly before the measurements, therefore spectra do not contain signal at 3429 cm −1 corresponding to stretching vibrations of Fe-O-H (water absorbed on the particles surface) [25].

Spectroscopic Measurements
To analyze the structure of used fillers and obtained nanocomposites, FT-IR measurements were performed. Results are presented below on Figures 2 and 3. Lack of any signals on the GNP spectra points towards the lack of functional groups, which can be related to results obtained by Geng et al. [24]. Fe 3 O 4 and mFe 3 O 4 particles were dried in vacuum directly before the measurements, therefore spectra do not contain signal at 3429 cm −1 corresponding to stretching vibrations of Fe-O-H (water absorbed on the particles surface) [25].

Spectroscopic Measurements
To analyze the structure of used fillers and obtained nanocomposites, FT-IR measurements were performed. Results are presented below on Figures 2 and 3. Lack of any signals on the GNP spectra points towards the lack of functional groups, which can be related to results obtained by Geng et al. [24]. Fe3O4 and mFe3O4 particles were dried in vacuum directly before the measurements, therefore spectra do not contain signal at 3429 cm −1 corresponding to stretching vibrations of Fe-O-H (water absorbed on the particles surface) [25].   Obtained infrared absorption spectra of manufactured composites are very similar to each other-the synthesis did not include introduction of strong functional groups into the system. Characteristic signals for C-H stretching vibrations are visible at: 2828 and 2930 cm −1 (in the -CH2group), and 1390 cm −1 (in the -CH3 group). At 1100 cm −1 , the signal for stretching vibrations of C-O of ether group is found.
The presence of the urethane groups is the cause of the absorption signal at ca. 1700 cm −1 , which corresponds to the stretching vibrations of C=O bond-this signal is also called the 1st amide signal, while the 2nd amide signal is located at 1530 cm −1 , which is characteristic for deformative vibrations of N-H.
At 3300 cm −1 , signals from secondary amides, i.e., N-H stretching vibrations, can be found. Due to the nature of the MU5 polyurethane, these signals can cover/envelope vibrations from KH 550 used to modify previously mentioned iron oxides-interactions with the filler could also cause slight shifts of this band towards lower wavenumbers.
In the MU5 composite, a small signal at 563 cm −1 is visible corresponding to the O=Fe-O bond in mFe3O4 particles. The bond is clearly visible in the filler spectra. The spectral range started at 500 cm −1 , which is why the signals at low wave number are not very clear.
Analyzed FT-IR spectra clearly show that neither GNP nor mFe3O4 form new chemical bonds in polymer matrix-all changes in the composite properties are due to physical interactions between fillers and polymer chains [26,27].
The crystal structure and size distribution was analyzed using the XRD method, spectra for both fillers are presented on Figure 4. The diffraction image of GNP contains the characteristic signal at shift 2θ = 26.4°, which corresponds to d002 plane, while distance between layers is ca. 3.4 Å [28]. Diffractogram of Fe3O4 shows signals for multiple planes (220), (311), (400), (511) and (440) for shifts 30°, 35°, 42°, 57° and 62° accordingly [16]. Obtained infrared absorption spectra of manufactured composites are very similar to each other-the synthesis did not include introduction of strong functional groups into the system. Characteristic signals for C-H stretching vibrations are visible at: 2828 and 2930 cm −1 (in the -CH 2 -group), and 1390 cm −1 (in the -CH 3 group). At 1100 cm −1 , the signal for stretching vibrations of C-O of ether group is found.
The presence of the urethane groups is the cause of the absorption signal at ca. 1700 cm −1 , which corresponds to the stretching vibrations of C=O bond-this signal is also called the 1st amide signal, while the 2nd amide signal is located at 1530 cm −1 , which is characteristic for deformative vibrations of N-H.
At 3300 cm −1 , signals from secondary amides, i.e., N-H stretching vibrations, can be found. Due to the nature of the MU5 polyurethane, these signals can cover/envelope vibrations from KH 550 used to modify previously mentioned iron oxides-interactions with the filler could also cause slight shifts of this band towards lower wavenumbers.
In the MU5 composite, a small signal at 563 cm −1 is visible corresponding to the O=Fe-O bond in mFe 3 O 4 particles. The bond is clearly visible in the filler spectra. The spectral range started at 500 cm −1 , which is why the signals at low wave number are not very clear.
Analyzed FT-IR spectra clearly show that neither GNP nor mFe 3 O 4 form new chemical bonds in polymer matrix-all changes in the composite properties are due to physical interactions between fillers and polymer chains [26,27].
The crystal structure and size distribution was analyzed using the XRD method, spectra for both fillers are presented on Figure 4. The diffraction image of GNP contains the characteristic signal at shift 2θ = 26.4 • , which corresponds to d 002 plane, while distance between layers is ca. 3.4 Å [28]. Diffractogram of Fe 3 O 4 shows signals for multiple planes (220), (311), (400), (511) and (440) for shifts 30 • , 35 • , 42 • , 57 • and 62 • accordingly [16].  XRD spectra of all investigated composites ( Figure 5) shows wide bands at angle shift 2θ = 20°, which can mean that polyurethane has mostly amorphous structure. At higher shifts, several diffraction maxima can be seen that can be attributed to mFe3O4 particles. Lack of significant bands for GNP can be assigned to its low content. Distinctive bands at shift 2θ = 26.6° and 54.7° for composites MU4 and MU5 can be attributed to GNP [26], which, in those two composites, is the most abundant. Lower intensity of the maxima (in comparison to GNP spectra) is a result of exfoliation of polymer matrix exfoliation. Shifts between layers in Fe3O4 and GNP particles are additionally specified on Figure 6, where XRD spectra for composite MU5 is presented  XRD spectra of all investigated composites ( Figure 5) shows wide bands at angle shift 2θ = 20 • , which can mean that polyurethane has mostly amorphous structure. At higher shifts, several diffraction maxima can be seen that can be attributed to mFe 3 O 4 particles. Lack of significant bands for GNP can be assigned to its low content. Distinctive bands at shift 2θ = 26.6 • and 54.7 • for composites MU4 and MU5 can be attributed to GNP [26], which, in those two composites, is the most abundant. Lower intensity of the maxima (in comparison to GNP spectra) is a result of exfoliation of polymer matrix exfoliation. Shifts between layers in Fe 3 O 4 and GNP particles are additionally specified on Figure 6, where XRD spectra for composite MU5 is presented  XRD spectra of all investigated composites ( Figure 5) shows wide bands at angle shift 2θ = 20°, which can mean that polyurethane has mostly amorphous structure. At higher shifts, several diffraction maxima can be seen that can be attributed to mFe3O4 particles. Lack of significant bands for GNP can be assigned to its low content. Distinctive bands at shift 2θ = 26.6° and 54.7° for composites MU4 and MU5 can be attributed to GNP [26], which, in those two composites, is the most abundant. Lower intensity of the maxima (in comparison to GNP spectra) is a result of exfoliation of polymer matrix exfoliation. Shifts between layers in Fe3O4 and GNP particles are additionally specified on Figure 6, where XRD spectra for composite MU5 is presented

Thermal Investigation
The XRD analysis showed low degree of crystallinity in tested materials which is also visible on DSC thermograms-there are no maxima for changes in crystalline structure. Thus, it is not possible to use DSC as an accurate method verifying glass-transition temperature (it can be established from DMA measurements). Maxima visible on Figure 7 correspond to melting temperature of soft (SS) and hard segments (HS) in the polyurethane-temperature values are presented in Table 3. Additional analysis of surface area, below the maxima for hard segments melting range, allows estimating melt solid-liquid transition enthalpy. Melting temperature of soft segments is relatively low (from 13 to 20 °C), which can be used to obtain the shape memory effect for composites. In order to fix temporary shape, the material has to be cooled down (and kept) below 10 °C. During storage, the temperature should not exceed Tmss, even if the programmed recovery temperature is higher.

Thermal Investigation
The XRD analysis showed low degree of crystallinity in tested materials which is also visible on DSC thermograms-there are no maxima for changes in crystalline structure. Thus, it is not possible to use DSC as an accurate method verifying glass-transition temperature (it can be established from DMA measurements). Maxima visible on Figure 7 correspond to melting temperature of soft (SS) and hard segments (HS) in the polyurethane-temperature values are presented in Table 3. Additional analysis of surface area, below the maxima for hard segments melting range, allows estimating melt solid-liquid transition enthalpy.

Thermal Investigation
The XRD analysis showed low degree of crystallinity in tested materials which is also visible on DSC thermograms-there are no maxima for changes in crystalline structure. Thus, it is not possible to use DSC as an accurate method verifying glass-transition temperature (it can be established from DMA measurements). Maxima visible on Figure 7 correspond to melting temperature of soft (SS) and hard segments (HS) in the polyurethane-temperature values are presented in Table 3. Additional analysis of surface area, below the maxima for hard segments melting range, allows estimating melt solid-liquid transition enthalpy. Melting temperature of soft segments is relatively low (from 13 to 20 °C), which can be used to obtain the shape memory effect for composites. In order to fix temporary shape, the material has to be cooled down (and kept) below 10 °C. During storage, the temperature should not exceed Tmss, even if the programmed recovery temperature is higher. Melting temperature of soft segments is relatively low (from 13 to 20 • C), which can be used to obtain the shape memory effect for composites. In order to fix temporary shape, the material has to be cooled down (and kept) below 10 • C. During storage, the temperature should not exceed Tmss, even if the programmed recovery temperature is higher. DMA results are compatible with the mechanical measurements, as described below (in Section 12). Results presented on Figure 8 show relations between storage modulus and temperature-the addition of filler improved mechanical properties of almost all of the composites. Drop in mechanical properties of MU5 can be attributed to the "oversaturation" with the filler. MU4 sample once again shows the best mechanical properties, where the storage modulus is over four times higher than for the unfilled polymer material.  DMA results are compatible with the mechanical measurements, as described below (in Section 12). Results presented on Figure 8 show relations between storage modulus and temperature-the addition of filler improved mechanical properties of almost all of the composites. Drop in mechanical properties of MU5 can be attributed to the "oversaturation" with the filler. MU4 sample once again shows the best mechanical properties, where the storage modulus is over four times higher than for the unfilled polymer material.   Maxima of phase angle relations (tan δ) and temperature ( Figure 9) allowed the establishment of glass transition temperature of subsequent materials (results are summarized in Table 3). All of the tested materials have T g value in a narrow range of −50.7 to −48.8 • C. Similar shape of all curves suggests close if not the same pattern of energy loss. At room temperature, materials are rather elastic, and their stiffness is imposed by the materials thickness.  DMA results are compatible with the mechanical measurements, as described below (in Section 12). Results presented on Figure 8 show relations between storage modulus and temperature-the addition of filler improved mechanical properties of almost all of the composites. Drop in mechanical properties of MU5 can be attributed to the "oversaturation" with the filler. MU4 sample once again shows the best mechanical properties, where the storage modulus is over four times higher than for the unfilled polymer material.   Thermogravimetric analysis results are presented on Figures 10 and 11. The weight loss from temperature shows that the addition of fillers only slightly improves thermal stability of obtained composites, as the temperature values remain on similar level. The target content of fillers is at least 30%w, but weight loss in composites MU2-MU5 are higher, which points towards partial oxidation of the iron oxides. Maxima of phase angle relations (tan δ) and temperature ( Figure 9) allowed the establishment of glass transition temperature of subsequent materials (results are summarized in Table 3). All of the tested materials have Tg value in a narrow range of −50.7 to −48.8 °C. Similar shape of all curves suggests close if not the same pattern of energy loss. At room temperature, materials are rather elastic, and their stiffness is imposed by the materials thickness.
Thermogravimetric analysis results are presented on Figures 10 and 11. The weight loss from temperature shows that the addition of fillers only slightly improves thermal stability of obtained composites, as the temperature values remain on similar level. The target content of fillers is at least 30%w, but weight loss in composites MU2-MU5 are higher, which points towards partial oxidation of the iron oxides.  Differential thermogravimetric analysis (DTG, Figure 11) shows that the decomposition of materials started around 320 °C. The double maxima on DTG curve were caused by decomposition of hard and elastic segments in polyurethanes. Minima around 425 °C could originate from two sources: degradation of aromatic isocyanates [29] or degradation of iron oxide particles in the Maxima of phase angle relations (tan δ) and temperature ( Figure 9) allowed the establishment of glass transition temperature of subsequent materials (results are summarized in Table 3). All of the tested materials have Tg value in a narrow range of −50.7 to −48.8 °C. Similar shape of all curves suggests close if not the same pattern of energy loss. At room temperature, materials are rather elastic, and their stiffness is imposed by the materials thickness.
Thermogravimetric analysis results are presented on Figures 10 and 11. The weight loss from temperature shows that the addition of fillers only slightly improves thermal stability of obtained composites, as the temperature values remain on similar level. The target content of fillers is at least 30%w, but weight loss in composites MU2-MU5 are higher, which points towards partial oxidation of the iron oxides.  Differential thermogravimetric analysis (DTG, Figure 11) shows that the decomposition of materials started around 320 °C. The double maxima on DTG curve were caused by decomposition of hard and elastic segments in polyurethanes. Minima around 425 °C could originate from two sources: degradation of aromatic isocyanates [29] or degradation of iron oxide particles in the Differential thermogravimetric analysis (DTG, Figure 11) shows that the decomposition of materials started around 320 • C. The double maxima on DTG curve were caused by decomposition of hard and elastic segments in polyurethanes. Minima around 425 • C could originate from two sources: degradation of aromatic isocyanates [29] or degradation of iron oxide particles in the coupling agent shell [30]. Mentioned phenomena are present in all of the samples, which indicates the source being the isocyanates or products of the first degradation. The temperature values of successive degradation stages (DTG curve minima) are summarized in Table 4. Thermal properties of all examined materials do not show any significant differences between each other ( Table 3), regardless of filler content. However, thermomechanical properties change slightly with the addition of GNP in comparison to iron oxides which do not have much influence in that matter.

Rheological Properties
Investigated materials contain over 30%w of the fillers in relation to overall prepolymer mass. However, added fillers have entirely different particle size-micro mFe 3 O 4 and nano GNP. Therefore, the weight content can be deceiving in matters of discussing divergent bulk densities. That is why volume content could better illustrate the amount of fillers introduced to the system and the differences between each of them. In order to estimate the volume content (vt %) of fillers, the bulk density of all components had to be established. The bulk densities of prepolymer, mFe 3 O 4 and GNP were 1.04, 2.00 and 0.06 g/cm 3 , respectively. Those values are however estimated and applicable only in this discussion and their purpose is to better illustrate fillers content. Both ways of showing different filler content in the material are presented in Figure 12 and Table 5. coupling agent shell [30]. Mentioned phenomena are present in all of the samples, which indicates the source being the isocyanates or products of the first degradation. The temperature values of successive degradation stages (DTG curve minima) are summarized in Table 4. Thermal properties of all examined materials do not show any significant differences between each other ( Table 3), regardless of filler content. However, thermomechanical properties change slightly with the addition of GNP in comparison to iron oxides which do not have much influence in that matter.

Rheological Properties
Investigated materials contain over 30%w of the fillers in relation to overall prepolymer mass. However, added fillers have entirely different particle size-micro mFe3O4 and nano GNP. Therefore, the weight content can be deceiving in matters of discussing divergent bulk densities. That is why volume content could better illustrate the amount of fillers introduced to the system and the differences between each of them. In order to estimate the volume content (vt %) of fillers, the bulk density of all components had to be established. The bulk densities of prepolymer, mFe3O4 and GNP were 1.04, 2.00 and 0.06 g/cm 3 , respectively. Those values are however estimated and applicable only in this discussion and their purpose is to better illustrate fillers content. Both ways of showing different filler content in the material are presented in Figure 12 and Table 5.    (Figures 14 and 16, respectively). The higher the shear rate, the bigger the differences between analyzed prepolymers, but MU1 and MU3 are still very close to each other (Figures 17  and 18 (Figures 14 and 16, respectively). The higher the shear rate, the bigger the differences between analyzed prepolymers, but MU1 and MU3 are still very close to each other (Figures 17 and 18, respectively).     (Figures 14 and 16, respectively). The higher the shear rate, the bigger the differences between analyzed prepolymers, but MU1 and MU3 are still very close to each other (Figures 17 and 18, respectively).             With the increase in filler content, the average viscosity of the material increases, which is shown on Figure 19. It is not an entirely linear increase, furthermore the unfilled prepolymer has identical viscosity as the one with 20 vt %-they are both lower than the 13%v prepolymer. That difference was also seen while the materials were pressed during the shape memory tests. The MU2 composite was more problematic for thermal processing and the only one that required higher temperature and longer pressing, in order to obtain smooth surface (as seen in Figure 20 where both materials were pressed in the same conditions, however the surface appearance was significantly different). One can try to explain that in the following manner: the addition of iron oxide in MU2 impaired the polyurethanes flowability, but the consequent addition of GNP (MU3-MU5) improved thermal conductivity, which, in turn, enhanced the flowability of the system in higher temperatures (as in pressing). Nonetheless, this theory requires further investigation on that matter, for example if the effects can be seen as soon as temperature of 70 °C is reached.
High viscosity of filled prepolymers can prove quite an obstacle during composite production process, where the mixing stage is usually performed at higher shear rates. Again, that matter also needs further analysis-a good example of this could be change in mixing temperature.  With the increase in filler content, the average viscosity of the material increases, which is shown on Figure 19. It is not an entirely linear increase, furthermore the unfilled prepolymer has identical viscosity as the one with 20 vt %-they are both lower than the 13%v prepolymer. That difference was also seen while the materials were pressed during the shape memory tests. The MU2 composite was more problematic for thermal processing and the only one that required higher temperature and longer pressing, in order to obtain smooth surface (as seen in Figure 20 where both materials were pressed in the same conditions, however the surface appearance was significantly different). One can try to explain that in the following manner: the addition of iron oxide in MU2 impaired the polyurethanes flowability, but the consequent addition of GNP (MU3-MU5) improved thermal conductivity, which, in turn, enhanced the flowability of the system in higher temperatures (as in pressing). Nonetheless, this theory requires further investigation on that matter, for example if the effects can be seen as soon as temperature of 70 • C is reached. With the increase in filler content, the average viscosity of the material increases, which is shown on Figure 19. It is not an entirely linear increase, furthermore the unfilled prepolymer has identical viscosity as the one with 20 vt %-they are both lower than the 13%v prepolymer. That difference was also seen while the materials were pressed during the shape memory tests. The MU2 composite was more problematic for thermal processing and the only one that required higher temperature and longer pressing, in order to obtain smooth surface (as seen in Figure 20 where both materials were pressed in the same conditions, however the surface appearance was significantly different). One can try to explain that in the following manner: the addition of iron oxide in MU2 impaired the polyurethanes flowability, but the consequent addition of GNP (MU3-MU5) improved thermal conductivity, which, in turn, enhanced the flowability of the system in higher temperatures (as in pressing). Nonetheless, this theory requires further investigation on that matter, for example if the effects can be seen as soon as temperature of 70 °C is reached.
High viscosity of filled prepolymers can prove quite an obstacle during composite production process, where the mixing stage is usually performed at higher shear rates. Again, that matter also needs further analysis-a good example of this could be change in mixing temperature. Figure 19. Influence of filler content on average prepolymer viscosity at CSR = 100 s −1 . Figure 19. Influence of filler content on average prepolymer viscosity at CSR = 100 s −1 .

Mechanical Properties
When modifying polymer composites with powder fillers, the oversaturation effect of the material can sometimes be noticed-the amount added is so high that the properties of the material are heavily impaired. GNP is often used in polyurethane composites to improve the mechanical properties, while the main objective of mFe3O4 addition is to introduce new, unique properties-such as magnetic response from the material. Surprisingly, incorporating solely iron oxides enlarged material elongation at break parameter ( Figure 21).
The best mechanical properties were achieved for the MU4 composite: breaking strain was 11 MPa and elongation 870 mm, which is several times higher than results for unmodified material-2.6 MPa and 117 mm, respectively. Despite a seemingly small difference between MU4 and MU5 (0.2%w/3.7%v in amount of fillers), signs of oversaturation in the material can be seen-breaking strain was 9.4 N/mm 2 and elongation at break 840 mm. Young's modulus of materials is low (in comparison to "standard" polyurethane results), however it was also improved by 80% for MU3 material. A summary of all parameters describing mechanical properties of the materials can be found in Table 6.
When analyzing results of this particular investigation it should be remembered that the materials were foamy. This could be the reason for low mechanical endurance, yet when analyzing breaking profiles, it was not determined that the break originated from air bubbles trapped in the polymer matrix.  High viscosity of filled prepolymers can prove quite an obstacle during composite production process, where the mixing stage is usually performed at higher shear rates. Again, that matter also needs further analysis-a good example of this could be change in mixing temperature.

Mechanical Properties
When modifying polymer composites with powder fillers, the oversaturation effect of the material can sometimes be noticed-the amount added is so high that the properties of the material are heavily impaired. GNP is often used in polyurethane composites to improve the mechanical properties, while the main objective of mFe 3 O 4 addition is to introduce new, unique properties-such as magnetic response from the material. Surprisingly, incorporating solely iron oxides enlarged material elongation at break parameter ( Figure 21).
The best mechanical properties were achieved for the MU4 composite: breaking strain was 11 MPa and elongation 870 mm, which is several times higher than results for unmodified material-2.6 MPa and 117 mm, respectively. Despite a seemingly small difference between MU4 and MU5 (0.2%w/3.7%v in amount of fillers), signs of oversaturation in the material can be seen-breaking strain was 9.4 N/mm 2 and elongation at break 840 mm. Young's modulus of materials is low (in comparison to "standard" polyurethane results), however it was also improved by 80% for MU3 material. A summary of all parameters describing mechanical properties of the materials can be found in Table 6.
When analyzing results of this particular investigation it should be remembered that the materials were foamy. This could be the reason for low mechanical endurance, yet when analyzing breaking profiles, it was not determined that the break originated from air bubbles trapped in the polymer matrix.

Shape-Memory Properties
The shape memory properties were investigated under four different sets of conditions, as stated in the Table 2.
The objective of the rest period is to check the fixation of the temporary shape. It is crucial not only to determine the storage conditions of such materials, but also the time while the temporary shape is still lasting. The change of shape during rest (Δ0) illustrates the fixation in time increments of 5 minutes. Even if a satisfactory deformation is obtained at the initial stage, it has to remain on the same (or at least similar) level for a given period of time. Differences in the initial shape between various composites is shown on Figure 22. Results of this investigation are demonstrated in Tables 7-10

Shape-Memory Properties
The shape memory properties were investigated under four different sets of conditions, as stated in the Table 2.
The objective of the rest period is to check the fixation of the temporary shape. It is crucial not only to determine the storage conditions of such materials, but also the time while the temporary shape is still lasting. The change of shape during rest (∆ 0 ) illustrates the fixation in time increments of 5 min. Even if a satisfactory deformation is obtained at the initial stage, it has to remain on the same (or at least similar) level for a given period of time. Differences in the initial shape between various composites is shown on Figure 22. Results of this investigation are demonstrated in Tables 7-10

Shape-Memory Properties
The shape memory properties were investigated under four different sets of conditions, as stated in the Table 2.
The objective of the rest period is to check the fixation of the temporary shape. It is crucial not only to determine the storage conditions of such materials, but also the time while the temporary shape is still lasting. The change of shape during rest (Δ0) illustrates the fixation in time increments of 5 minutes. Even if a satisfactory deformation is obtained at the initial stage, it has to remain on the same (or at least similar) level for a given period of time. Differences in the initial shape between various composites is shown on Figure 22. Results of this investigation are demonstrated in Tables 7-10 Table 7. The percentage value of shape change of the investigated materials (Δ0) during the rest period. I  II  III  IV  MU1  23%  10%  5%  6%  MU2  27%  20%  11%  7%  MU3  49%  22%  11%  7%  MU4  44%  20%  9%  2%  MU5 39% 13% 10% 2%  Lower rest temperature has major influence on the shape retention during that time period, as the melt temperature of the soft segments is below the room temperature. Materials resting at 4 • C had lower ∆ 0 than materials programmed at the same temperature, but resting at 20 • C-the change was even three times lower, as seen in the MU5 material between the I and II set. Higher programming temperatures led to decrease of ∆ 0 , which is more prominent in materials with higher filler content. Lower rest temperature has major influence on the shape retention during that time period, as the melt temperature of the soft segments is below the room temperature. Materials resting at 4 °C had lower Δ0 than materials programmed at the same temperature, but resting at 20 °C-the change was even three times lower, as seen in the MU5 material between the I and II set. Higher programming temperatures led to decrease of Δ0, which is more prominent in materials with higher filler content. The materials were bent in half, yet the highest obtained deformation was 50% (90°) at Tp = 90 °C. The programming effectiveness seems to be dependent on temperature, as the highest chain mobility is achieved. The heating medium during the research was water, therefore the highest achievable temperature was 100 °C. It is possible that even higher programming temperature would lead to better deformation of the materials, however that hinders materials manipulation possibility. Table 8 and Figures 23 and 24 present the shape retention rate after the rest period.  The materials were bent in half, yet the highest obtained deformation was 50% (90 • ) at T p = 90 • C. The programming effectiveness seems to be dependent on temperature, as the highest chain mobility is achieved. The heating medium during the research was water, therefore the highest achievable temperature was 100 • C. It is possible that even higher programming temperature would lead to better deformation of the materials, however that hinders materials manipulation possibility. Table 8  Lower rest temperature has major influence on the shape retention during that time period, as the melt temperature of the soft segments is below the room temperature. Materials resting at 4 °C had lower Δ0 than materials programmed at the same temperature, but resting at 20 °C-the change was even three times lower, as seen in the MU5 material between the I and II set. Higher programming temperatures led to decrease of Δ0, which is more prominent in materials with higher filler content. The materials were bent in half, yet the highest obtained deformation was 50% (90°) at Tp = 90 °C. The programming effectiveness seems to be dependent on temperature, as the highest chain mobility is achieved. The heating medium during the research was water, therefore the highest achievable temperature was 100 °C. It is possible that even higher programming temperature would lead to better deformation of the materials, however that hinders materials manipulation possibility. Table 8 and Figures 23 and 24 present the shape retention rate after the rest period.   Materials presented a rather low shape retention rate (Rr). Furthermore, they spontaneously tried to regain their original shape (∆ 0 ) during the rest period. The process was faster when the rest temperature of the material was higher than the melting temperature of the soft segments. The speed of shape recovery process increased exponentially after heating up the material. Low Rr is a result of high filler content, which disturbs the creation of segment structure, that allows creating the temporary shape of the material. Cai et al. [31] embodied MWCNT into the polymer chains (instead of a classic chain extender), and therefore they did not disturb the 3D structure of the chains, as they were incorporated. Although obtained results are not encouraging regarding material applications, it is possible that fine-tuning of the procedures could lead to better results. R r and ∆ 0 are dependent on three factors: the chemical composition, programming temperature and the rest temperature. It is linked to materials thermal properties-GNP addition increases thermal conductivity. With connection with high programming temperatures, sufficient chain mobility can be reached and the chains can be frozen in the temporary shape.

Nanocomposite
An important parameter while discussing SMPs is the time required for the material to recover its permanent shape after reaching T trans (in this case 45 • C). Results of time measurements results are presented in Table 9 and illustrated on Figure 25. Quite a surprise is the influence of the rest temperature on recovery time for materials programmed at the same temperature (sets I and II) which is between 20% and 60%. The addition of mFe 3 O 4 had no noticeable influence on shape recovery, while 1.2% addition of GNP can reduce the time value by half. The programing temperature does not seem to have strong influence on the recovery time. However, the connection with T p and R f cannot be ignored, despite high R f values for the first cycle, recovery time was noticeably higher. In addition, the temperature of shape recovery process can also influence the time value-when the recovery process was conducted at 50 • C, sample MU5 needed only 2.3 s for a full shape recovery. Example photos of the recovery process are presented on Figure 26.  Some important criteria when judging the possibility of materials application, is the ability to regain the original shape. After the recovery process, all materials regained the 180° angle, the lowest value was 178° at one attempt. The filler content seems to have no strong influence on the recovery ability of the material. Performing multiple cycles (15 repetitions, for set III) on one sample did not show any change in the shape memory behavior of the sample).   Some important criteria when judging the possibility of materials application, is the ability to regain the original shape. After the recovery process, all materials regained the 180° angle, the lowest value was 178° at one attempt. The filler content seems to have no strong influence on the recovery ability of the material. Performing multiple cycles (15 repetitions, for set III) on one sample did not show any change in the shape memory behavior of the sample).  Some important criteria when judging the possibility of materials application, is the ability to regain the original shape. After the recovery process, all materials regained the 180 • angle, the lowest value was 178 • at one attempt. The filler content seems to have no strong influence on the recovery ability of the material. Performing multiple cycles (15 repetitions, for set III) on one sample did not show any change in the shape memory behavior of the sample). It is possible that such high recovery rate is a result of initial programming when samples were initially pressed. To check this theory further, investigation should be conducted by evaluating shape memory properties on a material that was not thermally processed. There is, however, a restriction regarding uniform (non-foamy), thick materials which could undergo previously mentioned programming procedure. Solvent casting could be used to obtain a thin film, however it requires the use of organic solvents.
In previous research, it is uncommon to find reference to SMPs that were initially pressed. However, references were found to PLA/GNP [28] and PCL/Fe 3 O 4 [8] composites where the shape recovery rate was 83% and 80%, respectively. As the influence of pressing is nowhere evaluated, the differences cannot be attributed to the chosen fillers, as in all materials analyzed in this work (including the unfilled polyurethane matrix) have very similar, high level of R r . Another thermal processing method used for SMPs was the injection method, yet that reference source did not contain the R f and R r values [32].
The next step would be to proceed with investigations regarding shape memory properties of composites placed in a magnetic field. As GNP shows no such properties, the magnetic response of these materials could resemble response of PCL/MWCNT/mFe 3 O 4 composites, as seen on Figure 27. It is possible that such high recovery rate is a result of initial programming when samples were initially pressed. To check this theory further, investigation should be conducted by evaluating shape memory properties on a material that was not thermally processed. There is, however, a restriction regarding uniform (non-foamy), thick materials which could undergo previously mentioned programming procedure. Solvent casting could be used to obtain a thin film, however it requires the use of organic solvents.
In previous research, it is uncommon to find reference to SMPs that were initially pressed. However, references were found to PLA/GNP [28] and PCL/Fe3O4 [8] composites where the shape recovery rate was 83% and 80%, respectively. As the influence of pressing is nowhere evaluated, the differences cannot be attributed to the chosen fillers, as in all materials analyzed in this work (including the unfilled polyurethane matrix) have very similar, high level of Rr. Another thermal processing method used for SMPs was the injection method, yet that reference source did not contain the Rf and Rr values [32].
The next step would be to proceed with investigations regarding shape memory properties of composites placed in a magnetic field. As GNP shows no such properties, the magnetic response of these materials could resemble response of PCL/MWCNT/mFe3O4 composites, as seen on Figure 27. Fe3O4 particles vibrations in the magnetic field are the source of thermal energy that leads to internal temperature increase up to Ttrans value (45 °C). Therefore, the magnetic shape recovery investigation should include the response time for the internal temperature to increase up to the required value. Carbon nanofillers have good thermal conductivity, which provides quick internal heating of the material in its entire volume which speeds up the shape recovery process. Multiwall carbon nanotubes and graphene nanoplates have the same thermal conductivity-3000 W/mK [33], therefore the influence of both fillers on the composites heating rate should be the same. Table 11 contains comparison of materials shape memory properties (time and Rr) programmed in the same conditions (Tp = 60 °C), unfortunately the source articles do not contain the data required here, i.e., respective Rf values. Polyurethane composites obtained in this research have higher shape recovery rate and lower response time. Fe 3 O 4 particles vibrations in the magnetic field are the source of thermal energy that leads to internal temperature increase up to T trans value (45 • C). Therefore, the magnetic shape recovery investigation should include the response time for the internal temperature to increase up to the required value. Carbon nanofillers have good thermal conductivity, which provides quick internal heating of the material in its entire volume which speeds up the shape recovery process. Multiwall carbon nanotubes and graphene nanoplates have the same thermal conductivity-3000 W/mK [33], therefore the influence of both fillers on the composites heating rate should be the same. Table 11 contains comparison of materials shape memory properties (time and R r ) programmed in the same conditions (T p = 60 • C), unfortunately the source articles do not contain the data required here, i.e., respective R f values. Polyurethane composites obtained in this research have higher shape recovery rate and lower response time. Table 11. Memory shape properties at temperature 45 • C in water bath of materials obtained by Cai et al. [11,31] and materials obtained in this work according to the cycle II.

Conclusions
After analyzing composite properties, we can notice a positive influence of added graphene nanoplates and ferromagnetic iron oxide.
Thermal properties of all composites remained on a similar level, however the mechanical properties increased significantly. Comparing to unfilled material (MU1), the breaking strain of composites MU4 and MU5 increased by 328% and 263%, respectively, while elongation at break increased by 641% and 616%.
With respect to memory shape properties, temperature trends described above cannot be ignored. At lower temperatures, the MU5 material had better retention rate (+17% compared to MU1), while, at higher temperatures, it increased by 43% and 40% for MU4 and MU5 materials, respectively. A seemingly small difference in their composition (0.2%w GNP) had high impact on time necessary for the shape recovery, it was 30%, or even 65% faster for MU4 and MU5 than for MU1 at the same temperature. However, discussing the "time" parameter, one has to bear in mind the actual difference between the temporary and permanent shape, as well as the desired application of such element. Presence of high strains during the shape recovery process might lead to micro cracks appearing in the material which could lower overall mechanical endurance.
In conclusion, a new, memory shape composite was successfully obtained. Fillers introduced to the polymer matrix lead to increase in some of the material properties, however, due to their powder nature, they might have disturbed overall segment structure inside the matrix.
Further investigation of this topic might include measurements in the magnetic field, while additional modifications of the composition and the programming process could lead to better results being obtained in near future.