Introduction

Natural cellulosic fibres such as wood pulps are well known reinforcing agents in composites and the high strength and large aspect ratio of nano-scale cellulose fibrils make them suitable as reinforcement in composite materials (Berglund and Peijs 2010; Oksman et al. 2016). Recent developments in cellulose pulp processing have made nano-sized cellulose fibrils commercially available (Saito et al. 2007). Since cellulose is one of the most abundant organic polymer in the world, it has the potential to be an inexpensive reinforcing filler material in thermoplastics. A major drawback of using cellulose nanofibrils (CNF) as reinforcement in thermoplastics is however the difference in hydrophilic nature between cellulose and the usually hydrophobic matrix phase polymer. Moreover, it can be argued that dewatering and drying of cellulose nanomaterials promotes irreversible bonding between nanofibrils, as discussed in Fernandes Diniz et al. (2004). This incompatibility can affect the mixing of the cellulose in the matrix and can lead to aggregation of cellulose fibre/fibrils resulting in a non-uniform thermoplastic composite (Ariño and Boldizar 2012). Insufficient compatibility can also result in poor fibre–matrix adhesion, leading to poor stress transfer and inferior mechanical properties. In addition, the low thermal degradation temperature of cellulose limits the range of processing temperatures and the selection of applicable polymer matrices. Common matrices used in melt extrusion of CNF are polyesters such as poly(lactic acid) (PLA) (Jonoobi et al. 2010), due to its renewable origin and biodegradability. Other common matrices used are polyethylene and polypropylene (Boldizar et al. 1987), which have relatively low melting temperatures and are common bulk polymers.

To create well dispersed CNF composite materials, methods that focus mainly on cellulose surface chemistry have been used (Dufresne 2017). These include grafting hydrophobic compounds to the cellulose surface to make the fibrils more hydrophobic, often combined with a solvent exchange step (Siqueira et al. 2009; de Menezes et al. 2009). Surfactants or the surface adsorption of amphiphilic co-polymers have also been used to improve the dispersion and compatibility of the CNF in the polymer matrix (Volk et al. 2015). These methods have been effective, but the use of non-polar solvents and multiple reaction steps may limit their industrial relevance.

The method of solvent casting CNF together with a water-soluble polymer is known to yield a well dispersed composite with excellent mechanical properties (Sehaqui et al. 2011; Srithep et al. 2012). However, the use of a water-soluble matrix polymer makes the resulting composite moisture-sensitive. To avoid this, an insoluble aqueous-dispersed polymer can be used to solvent cast composites. This approach is not entirely new, early efforts include the use of cellulose nanocrystals in combination with natural rubber latex as a matrix phase (Favier et al. 1995). Water-assisted mixing of CNF and micrometer-sized PLA latex has been used to create fully biodegradable composites with improved tensile properties (Larsson et al. 2012). More recently, a aqueous dispersion of an elastomeric co-polymer using poly(ethylene-acrylic acid) as stabilizing surfactant was used to create CNF composite films with improved tensile stiffness (Maia et al. 2017).

In the present work, poly(ethylene-co-acrylic acid) (EAA) reinforced with unmodified CNF was mixed according to the water-assisted method, with the purpose to explore the possibility to disperse high contents of CNF in a polymer matrix. The EAA polymer was selected mainly due to it being insoluble in water and having the ability to form a stable aqueous dispersion. Also, this modified polyethylene had previously shown good adhesion and compatibility with regenerated cellulose (Saarikoski et al. 2012). Composites with up to 70 vol% dry content CNF could be made. To better understand the general capability of CNF to reinforce a polymer matrix, the thermal and mechanical properties were studied in comparison with those of plain cellulose pulp composites. The cellulose pulp selected for comparison consisted mainly of spruce sulphite and was of a highly beaten type with a similar surface structure as the CNF used.

Methods

Materials

Poly(ethylene-co-acrylic acid) (EAA) dispersion with a 20 wt% solids content was obtained from BIM Kemi AB, Sweden. The EAA polymer had an acrylic acid content of 15%, a density of 0.994 g/cm3, a melting point of \(88\,{^{\circ }}\hbox {C}\) and a melt flow rate of 36 g/10 min (ISO 1133, \(190\,{^{\circ }}\hbox {C}\), 2.16 kg), according to the supplier. The dispersed EAA was neutralized to its ionomer form with NaOH during the dispersing process. The pH of the EAA dispersion was 9.7. CNF paste made from spruce sulphite pulp was kindly provided by Borregaard AS, Norway. The CNF had a solids content of 10 wt% in water. The CNF had a Klason lignin content of 3.2 wt%, a hemicellulose content (arabinoxylan and galactoglucomannan) of 2.8 wt% and a cellulose content of 94 wt%. These contents were estimated from monosaccharide carbohydrate analysis. The cellulose fibre pulp used was a highly beaten never-dried bleached softwood mixture of 80% spruce sulphite and 20% spruce sulphate pulp, kindly provided by Nordic Paper Seffle AB, Sweden. The pulp had a Klason lignin content of 1.9 wt%, a hemicellulose content of 17 wt% and a cellulose content of 80 wt%.

Rheological measurements

The rheological properties of the suspensions were measured using an Anton Paar MCR702, Twin-drive shear rheometer (Graz, Austria). The suspensions studied were the CNF paste, with a dry content of 10 wt% and additionally the CNF–EAA suspensions containing 0.7 wt% CNF and 1.8 wt% CNF, respectively. The linear viscoelastic range of these suspensions was characterized through oscillatory shear strain measurements in the range of 0.01–100% strain and a frequency of 1 Hz. The steady shear flow was investigated in the shear rate range of 0.1–\(100\,\hbox {s}^{-1}\), the measurements being performed at \(25\,{^{\circ }}\hbox {C}\) using a cone-plate fixture. The cone-plate geometry had a diameter of 25 mm, a cone angle of \(2^\circ \) and a gap distance of 0.106 mm. Temperature sweeps were performed on the 70 and 10 vol% CNF–EAA suspensions from 27 to \(90\,{^{\circ }}\hbox {C}\) at a frequency of 1 Hz and a strain amplitude of 0.2%. A parallel-plate fixture with a plate diameter of 15 mm and a gap distance of 0.5 mm was used for these measurements and silicone oil was used to reduce the evaporation of water.

Fibre analysis

The fibre length and width was determined according to Tappi standard T271 using a Kajaani FS300 fibre analyser (Metso Automation, Finland), where fines content is based on the centerline length. Fibre length and width are reported as mean values based on approximately 30,000 fibre pieces.

Dynamic light scattering

The hydrodynamic diameter of the EAA particles in suspension was measured using a Malvern Instruments Zetasizer ZEN3600. The data reported is the z-average from cumulant fit. The EAA dispersion was diluted to 0.01 wt% with deionized water and sonicated for 10 min. Measurements were conducted in a standard quartz cuvette at room temperature.

Composite manufacture

EAA dispersion was added to CNF paste or pulp to obtain the intended volume fractions of cellulose. The mixtures had a solids content of 2.5–4.5 wt% cellulose, and excess of water was added to achieve a total water content of \(96 \pm 1\)%. These suspensions were then mixed at room temperature in an L&W pulp disintegrator (Lorentzen & Wettre, Sweden) for 60,000 revolutions at 2900 rpm. Volume fractions were calculated assuming a cellulose density of 1.5 g/cm3 and an EAA density of 0.994 g/cm3. The mixed cellulose and EAA suspensions were dried at room temperature to about 0.1 mm thick sheets. Approximately 12 g of the dry sheets were stacked and compression moulded to composite plates with an area of \(100 \times 100\) mm2 and a thickness of about 1 mm using a Bucher–Guyer KHL 100, Switzerland. Stacks of dry sheets were placed in the open mould tempered at \(105\,{^{\circ }}\hbox {C}\). The mould was then closed and a pressure of 20 bar was applied until the distance between the mould halves had stabilized, which took approximately 5 min. The pressure was then raised to 500 bar and the mould was cooled to \(30\,{^{\circ }}\hbox {C}\), which took about 2 min. Finally, with the mould cooling turned off, the mould temperature was allowed to increase to \(40\,{^{\circ }}\hbox {C}\), which took about 1 min. The pressure was then released and the sample was removed from the mould.

Microscopy and tomography

Dried CNF paste was studied with a LEO Ultra 55 FEG Scanning Electron Microscopy (SEM). Samples were goldsputtered in vacuum for 80 s at 10 mA, to give a gold coating approximately 10 nm thick. The atomic force microscope (AFM) used was a Nanoscope IIIa with a type G scanner (Digital Instruments Inc.), with a Micro Masch silicon cantilever NSC 15. The measurements were performed in air and in the tapping mode. The 3D data of the internal structure of composites were scanned with a Zeiss Xradia XRM520 X-ray tomograph. The scanned volume were reconstructed to a size of 1 \(\upmu {\text {m}}^3\)/voxel. Graphs of voxel gray scale were generated with the ImageJ software. 3D images of composite microstructure is based on surface reconstructions made in the software 3D slicer.

Differential scanning calorimetry

Thermograms were prepared by Differential Scanning Calorimetry (DSC), using a Perkin Elmer DSC7 from 10 to \(250\,{^{\circ }}\hbox {C}\). The scan rate was \(10\,{^{\circ }}\hbox {C}\)/min and the sample mass was approximately 10 mg. Peak deconvolution was done using a three-term Gaussian function. The crystallinity of the EAA matrix was calculated using the Eq. 1.

$$\begin{aligned} X_c=\frac{\varDelta {\text {H}}_c}{w_{ EAA} \varDelta {\text {H}}_0} \end{aligned}$$
(1)

Here, \(\varDelta \hbox {H}_c\) is the specific heat of fusion of the composite, \(\hbox {w}_{ EAA}\) is the mass fraction EAA of the sample and \(\varDelta \hbox {H}_0\) is the specific heat of fusion for a polyethylene crystal (277.1 J/g) (Brandrup et al. 1999).

Dynamic mechanical thermal analysis

The viscoelastic properties of composites were evaluated using dynamic mechanical thermal analysis (DMTA). The measurements were performed with a Rheometrics Solids Analyzer RSA II at a frequency of 1 Hz, a heating rate of \(3\,{^{\circ }}\hbox {C}\)/min and a strain amplitude between 0.1 and 0.15%. The temperature range used was \(-\,80\,{^{\circ }}\hbox {C}\) to \(110\,{^{\circ }}\hbox {C}\). The sample thickness varied between 0.8 and 1 mm and the width between 4 and 7 mm. The samples were measured in tension within the linear viscoelastic range. Prior to the temperature sweep, a strain sweep was performed for each sample series at room temperature.

Tensile testing

The Young’s modulus (E), stress at break (\(\sigma _{{\text {b}}}\)) and elongation at break (\(\varepsilon _{{\text {b}}}\)), were evaluated in accordance with ISO 527-3. Samples were cut into dumb-bell shapes having a width of 4 mm and a thickness of approximately 1 mm and conditioned for at least 2 days at \(23 \pm 2\,{^{\circ }}\hbox {C}\) and \(50 \pm 5\)% relative humidity (RH) prior to testing, The tensile measurements were also made at \(23 \pm 2\,{^{\circ }}\hbox {C}\) and \(50 \pm 5\)% RH, with a Zwick Z1/Roell with a grip separation of 40 mm and the cross-head speed was 6 mm/min. A Zwick Eye UI 1540M video extensometer was used, with an initial gauge length of 20 mm. A load of 0.1 N was applied before measurement. The tensile stiffness measured was compared to the Cox–Krenchel model according to Eq. 2.

$$\begin{aligned} E_{\text {c}}=\eta _{d}E_{\text {f}}\left( 1-\frac{\tanh \left( \frac{\beta l}{2}\right) }{\frac{\beta l}{2}}\right) \phi +\left( 1-\phi \right) E_{\text {m}} \end{aligned}$$
(2)

where

$$\begin{aligned} \beta ^2=\frac{E_{\text {m}}/\left( 1+\nu _{\text {m}}\right) }{r^2E_{\text {f}}\ln \left( \frac{R}{r}\right) } \qquad {\text {and}} \qquad R=\sqrt{\frac{\pi r^2}{4\phi }} \end{aligned}$$

\(\hbox {E}_{\text {m}}\), \(\hbox {E}_{\text {f}}\) and \(\hbox {E}_{\text {c}}\) are the the elastic moduli of the EAA matrix, the cellulose fibre and the composite, respectively (Cox 1952; Thomason and Vlug 1996). \(\hbox {E}_{\text {m}}\) was measured to be 0.29 GPa and \(\hbox {E}_{\text {f}}\) was set to 21 GPa for cellulose fibres (Ehrnrooth and Kolseth 1984) and 32 GPa for the CNF (Tanpichai et al. 2012). \(\phi \) is the volume fraction of cellulose and \(\eta _{\text {d}}\) is a orientation factor which was assumed to be 3/8 corresponding to a random in-plane fibre orientation (Thomason and Vlug 1996). The Poisson ratio, \(\nu _{\text {m}}\) for the EAA matrix was assumed to be 0.3 and r is the fibre radius. Fibre dimensions were taken from the results of the fibre analysis. For CNF composites, a diameter of \(18.4\,\upmu \hbox {m}\) and a length of 0.34 mm were used as the fibre analysis neglects both the pulp fines and the cellulose fibrils. For pulp composites, a diameter of \(27\,\upmu \hbox {m}\) and a length of 1.3 mm were used.

The effective stiffness of CNF (\(\hbox {E}_{\text {f,l}}\)) was calculated using the method of Ansari et al. (2014), through the Halpin–Tsai model.

$$\begin{aligned} E_{\text {f,l}}=\frac{\frac{3}{8} \left( E_{\text {c}} - \frac{5}{8}E_{\text {T}}\right) -E_{\text {m}}(1-\phi )}{\phi } \end{aligned}$$
(3)

where

$$\begin{aligned} E_{\text {T}}=E_{\text {m}}\left( \frac{1+2\eta \phi }{1-\eta \phi }\right) \qquad {\text {and}} \qquad \eta =\frac{E_{{\text {f,t}}}-E_{\text {m}}}{E_{\text {f,t}}+2E_{\text {m}}} \end{aligned}$$

The transverse fibre modulus in Eq. 3 \(\hbox {E}_{\text {f,t}}\) was set to 15 GPa (Diddens et al. 2008).

Results and discussion

Rheological properties

The shear storage modulus, (\({\text {G}}^\prime \)) at 1 Hz is shown as a function of the shear strain in Fig. 1a. The graph clearly shows a difference in \({\text {G}}^\prime \) between the CNF paste and the CNF–EAA suspension with a water content of \(96 \pm 1\) wt%. The lower \({\text {G}}^\prime \) of the CNF–EAA suspensions was expected, due to the reduction in CNF concentration by the addition of EAA dispersion and the excess of water. The transition from the linear viscoelastic region to the non-linear region increased with decreasing CNF concentration, the addition of EAA and the excess of water as expected (Quennouz et al. 2015). The critical strain value, indicated by the drop in storage modulus, increased from 1.1% for the CNF paste to 2.3% for the CNF–EAA suspension with 0.7 wt% CNF dry content. These critical strains measured were estimated to corresponds to critical stress values of around 250, 7.6 and 2.4 Pa for the CNF paste, the 1.8 and 0.7 wt% suspensions respectively, and was coupled to the yield stress of the samples. Interestingly, the CNF content in the suspensions did not appear to affect the slope in the region of shear thinning. This observation was in contrast to previous work by Naderi et al. (2014), who showed an increase in the slope after the linear viscoelastic range with increasing CNF content of the suspension. Furthermore, it was here seen that the loss modulus (\({\text {G}}^{\prime \prime }\)) was lower than the storage modulus \({\text {G}}^\prime \) up to about 6 % strain indicating the more prominent elastic character of the CNF–EAA suspensions. Figure 1b shows the steady shear viscosity as a function of shear rate. A decreasing shear viscosity with increasing shear rate was observed, in agreement with earlier work (Moberg and Rigdahl 2012). The shear viscosity deceased with decreasing CNF content as expected. The viscosity could not be measured on the CNF paste at shear rates above 5 \(\hbox {s}^{-1}\) because of slippage between the sample and the cone-plate fixture. Figure 2 shows the \({\text {G}}^\prime \) as a function of temperature for the CNF–EAA suspensions with 0.7 and 1.8 wt% CNF dry content, respectively. The high water content (\(96 \pm 1\)%) in both the suspensions limited the temperature range of the experiments due to the evaporation of water, as indicated by the increase in \({\text {G}}^\prime \) with temperature. The CNF–EAA suspension with 0.7 wt% CNF dry content and 3.95 wt% EAA showed a peak in \({\text {G}}^\prime \) at \(59\,{^{\circ }}\hbox {C}\), which is close to the first melting point of EAA. This peak suggests that the increase in \({\text {G}}^\prime \) could be related to morphological changes in the EAA dispersion such as flocculation or coalescence, rather than just evaporation, at these temperatures. The peak was observed only in the sample with 0.7 wt% CNF dry content due to the high EAA concentration in the sample.

Fig. 1
figure 1

a The dynamic shear storage modulus (\({\text {G}}^{\prime }\), solid line) and the shear loss modulus, (\({\text {G}}^{\prime \prime }\), dashed line) as a function of shear strain and b the steady shear viscosity as a function of shear rate for CNF paste and CNF–EAA suspensions with 1.8 and 0.7 wt% CNF dry contents

Fig. 2
figure 2

The variation of dynamic shear storage modulus (\({\text {G}}^\prime \)) as a function of temperature for CNF–EAA suspensions with 1.8 and 0.7 wt% CNF dry contents

Fig. 3
figure 3

Scanning electron micrograph of a dried CNF and b AFM image of dried CNF

Cellulose and EAA morphology

Fibre analysis showed that fibres or large fibre fragments were still present in the CNF starting material. These fibres were however shorter than those in the pulp (Table 1). The large fibre fragments present in the CNF were probably the result of an incomplete fibrillation process. The fibrillated nature of the CNF was reflected in the fines content, see Table 1, which was much higher for the CNF than for the highly beated pulp.

Table 1 Fibre analysis of CNF and pulp
Fig. 4
figure 4

EAA–CNF-20 a X-ray microtomography data showing three density phases for the EAA–CNF-20 material, b 3D rendering of the (blue) high density fibre phase and the (green) low density phase. (Color figure online)

Fig. 5
figure 5

EAA–pulp-20 a X-ray microtomography data showing three density phases for the EAA–Pulp-20 material, b 3D rendering of the (blue) high density fibre phase and the (green) low density phase

The microstructure of the CNF used in this work is shown in Fig. 3a. Scanning electron microscopy revealed that the cellulose fibrils of the studied CNF had diameters and lengths largely on the micrometer scale. The smallest fibrils present in the CNF had diameters on the nanoscale (< 100 nm) as shown by AFM (Fig. 3b). The CNF appears to be a mixture of fibres, micro- and nanofibrils. The relative contents of these fibres, micro- and nanofibrils could not be determined, since none of the analysis techniques covered the whole size range of the CNF material. The EAA dispersion used to make the composites was milky in appearance and the mean particle size was 50 nm, determined by dynamic light scattering. This particle size was considerably smaller than other micrometer sized latex dispersions that have been used to produce CNF composites (Larsson et al. 2012).

Composite microstructure

The microstructure of the cellulose-EAA composites was determined using X-ray microtomography with a reconstructed voxel size of 1 \(\upmu \hbox {m}^3\)/voxel. Figure 4a shows the number of volume elements against the voxel gray scale for the EAA–CNF-20 composite material containing 20 vol% cellulose. As, the gray scale corresponds to the density of the material, three distinct phases could be distinguished. It was indicated that the phase with the highest voxel density corresponded to cellulose fibres and fibrils, since fibrous micro-structures were clearly seen in Fig. 4b. The volume fraction of the high density phase was however smaller than the total cellulose fraction of the composite, which suggests that the medium density phase was a mixture of EAA and cellulose fibrils smaller than the voxel size, 1 \(\upmu \hbox {m}^3\)/voxel. The third distinguishable phase had a low density, suggesting purely amorphous EAA or voids. In Fig. 4b, the high density fibre/fibril phase and the low density phase are shown together. Interestingly, the low density phase was situated on fibre/fibril surfaces and inside the lumen of fibres. These findings indicate a low density inter-phase near the cellulose surface. In Figure 5b it can be seen that a low density phase also existed in the pulp–EAA composites. In the pulp–EAA composite, this phase was also situated close to the fibre surface and inside the lumen of fibres. The volume fraction of 13% for the high density fibre phase in the pulp–EAA composite with 20 vol% pulp content indicates that there were small fibril elements in the pulp that fall below the detection limit of X-ray microtomography. These undetected cellulose fibrils were probably included in the medium density phase.

Thermal properties

The semi-crystalline nature of EAA can be seen in the thermograms in Fig. 6. Multiple endothermic melting peaks were detected in the temperature interval 50–\(90\,{^{\circ }}\hbox {C}\). A second heating after cooling at \(10\,{^{\circ }}\hbox {C}\)/min generated a more distinct melting peak at about \(88\,{^{\circ }}\hbox {C}\), consistent with the melting temperature of non-isothermally crystallized EAA (Zhang et al. 2009). For ionomers such as neutralized EAA, it is understood that the low temperature endotherm arises due to the gradual formation of thin crystals at room temperature after primary crystallization (Loo et al. 2005). The multiple melting peaks in all the samples after compression moulding indicated that thin crystals had formed. The melting peaks could be described by a three-term Gaussian function (Fig. 6a). The crystallinity of these two melting regions are shown in Table 2. The secondary crystallinity of EAA (\(\hbox {X}_{c1}\) in Table 2) formed gradually after compression moulding and the accuracy of this crystallinity was therefore uncertain since the time between compression moulding and DSC measurements was not the same for all samples. However, addition of cellulose to EAA did not significantly influence the primary crystallinity of EAA (\(\hbox {X}_{c2}\) in Table 2), which varied between 9 and 15% for the composite materials with no apparent trend. Previous works on cellulose composites have shown that cellulose reinforcement in composites can both induce and reduce crystallinity in the matrix phase (Samir et al. 2004; Yao et al. 2008; Maia et al. 2017). The crystallinity of EAA does not significantly change with the addition of cellulose, apart from a small reduction in crystallinity at higher cellulose contents, which suggests that any improvement in mechanical properties of the composites would probably be the result of the cellulose reinforcement and not due to the increased crystallinity of the EAA. The measured dynamic storage modulus and loss factor of the CNF–EAA composites over a broad temperature range are shown in Fig. 7. Three transitions can be discerned, the first at a temperature of \(-\,39 \pm 4\,{^{\circ }}\hbox {C}\), probably the glass transition temperature of the EAA matrix and a second transition at \(55 \pm 5\,{^{\circ }}\hbox {C}\), corresponding to the first melting temperature \(\hbox {T}_{\text {m1}}\) of the EAA matrix. The second melting peak, previously denoted as \(\hbox {T}_{\text {m2}}\) observed in the DSC analysis, is clearly seen by the loss factor, \(\tan \delta \), for samples containing 0–30 vol% CNF, see Fig. 7b. Addition of CNF not only increased the magnitude of \({\text {E}}^\prime \) up to the melting at about \(55\,{^{\circ }}\hbox {C}\) but also reduced the drop in \({\text {E}}^\prime \) after \(55\,{^{\circ }}\hbox {C}\), suggesting a general stiffening effect. The transitions decreased in magnitude with increasing CNF content, becoming insignificant above 40 vol% CNF.

Table 2 Thermal properties of cellulose-EAA composites
Fig. 6
figure 6

Thermograms and fitted Gauss-peaks of the first and second heating scans for a neat compression moulded EAA, b EAA–CNF-50 composite

Fig. 7
figure 7

a Storage modulus (\({\text {E}}^\prime \)) and b loss factor (\(\tan \delta \)) of the CNF–EAA nanocomposites studied in the temperature range from \(-\,80\) to \(110\,{^{\circ }}\hbox {C}\) measured at a frequency of 1 Hz and a heating rate of \(3\,{^{\circ }}\hbox {C}\)/min

Tensile properties

Typical stress-strain curves of the composites are shown in Fig. 8. Neat EAA was a highly ductile polymer with an elongation at break (\(\varepsilon _{\text {b}}\)) of about 500% and a low elastic modulus (E), see Table 3. The high \(\varepsilon _{\text {b}}\) of the EAA was significantly reduced by the addition of cellulose to the composites. As expected, the distinctive yield point seen in neat EAA at 13% strain disappeared when cellulose was added (see inset plot in Fig. 8). The cellulose composites instead showed a monotonous increase in stress until fracture. This behavior was observed for both CNF–EAA and pulp–EAA composites. Such deformation behavior has been attributed to deformation-induced rearrangements in the fibril/microfibril network (Ansari et al. 2014; Henriksson et al. 2008). In Fig. 9a, the Young’s modulus of the CNF–EAA composites is shown as a function of cellulose volume fraction. The modulus was observed to increase with increasing CNF content from 0.3 GPa for neat EAA to 6 GPa at 70 vol% CNF content, which gives a stiffening factor of 21. The Cox–Krenchel model in Fig. 9a was accurate at cellulose contents below 20 vol% and in the interval 40 to 60 vol%. The model underestimated the stiffness increase in the interval 20–40 vol%. The average fibril aspect ratio assumed in the model was however based on fibre analysis data which neglects both pulp fines and cellulose fibrils. Accurate modeling of the tensile properties would require more detailed knowledge about the fibril size and fibril content of the CNF. An increase in tensile properties of the CNF composites compared to those of the pulp composites could be expected since cellulose fibrils are known to have both a larger aspect ratio and a higher elastic modulus than pulp fibres (Tanpichai et al. 2012). The Cox–Krenchel model in Fig. 9b was accurate for the pulp–EAA composites, which was expected since this cellulose consisted mostly of pulp fibres and fibre analysis showed a lower fines content than the CNF.

Fig. 8
figure 8

Typical stress-strain curves of CNF–EAA and pulp–EAA composites. The inset shows the full stress–strain behavior for EAA

Fig. 9
figure 9

Elastic modulus along with the Cox–Krenchel prediction for a CNF–EAA composites and b pulp–EAA composites, tensile properties for all cellulose composites c stress at break and d strain at break

Table 3 Tensile properties of cellulose-EAA composites

A large reduction in effective stiffness \(\hbox {E}_{\text {f},l}\) with increasing CNF content has been linked to fibril aggregation in CNF composites (Ansari et al. 2014). The \(\hbox {E}_{\text {f},l}\) in Table 3, had a maximum value of 22.3 GPa for 30 vol% CNF, accompanied by a minor decrease to 18.6 GPa for 70 vol% CNF. This small decrease in effective stiffness suggests a small variation in dispersion of CNF in the EAA matrix in the range of 30–70 vol% CNF. Interestingly, the pulp composite were stiffer than the CNF composite at 20 vol% cellulose. The CNF composites had a stiffening threshold in the interval 20–30 vol% where the stiffness nearly doubled. This threshold concentration has been reported earlier in similar composites, Larsson et al. reported a significant stiffness increase in CNF/PLA composites at a cellulose concentration of 25 wt% (Larsson et al. 2012). Moreover, Maia et al. reported a doubling of the stiffness to 0.67 GPa in the interval 20–30 wt% for CNF and elastomer composites (Maia et al. 2017). The strength \(\sigma _{\text {b}}\) of the CNF composites was also affected by this concentration threshold (Fig. 9c). The fibril network in the composites seem to have a percolation point where it can form entanglements and enhance both stiffness and strength greatly. The CNF composites produced were measured to have a significantly higher strength than the pulp composites at a cellulose content above 20 vol%. Such an increase in strength was expected from previous work on similar composites mixed in the dry state (Ariño and Boldizar 2012). The present results indicate that the water-assisted mixing method has the potential to be better than dry mixing. Both CNF and pulp mixed in the wet state led to an increased strength by a factor of 1.5 at 20 vol% cellulose (Table 3). The composite with the highest loading of 70 vol% CNF had a strength of 83.3 MPa which was an increase by a factor of 3.5 compared to the EAA matrix.

The strain at break (\(\varepsilon _{\text {b}}\)) decreased with the addition of CNF, the lowest value being about 5% for the 60 vol% composite and the highest being 32% for the 10 vol% composite. The pulp–EAA composites generally had somewhat lower \(\varepsilon _{\text {b}}\) than the CNF–EAA composites (Fig. 9d). This ability to withstand deformation was generally larger for CNF-containing than for pulp-containing composites, implying higher mobility of the fibrils than of the fibres under stress and a better matrix-to-cellulose stress transfer. The morhpology of both CNF and pulp composites suggested potential voids near the cellulose surface (Figs. 4, 5). Such voids would have a negative impact on both the stiffness and strength of the EAA composites since matrix-to-cellulose stress transfer would not be optimal. The removal of such near-surface voids would most likely increase the mechanical properties of the CNF and pulp composites even further.

Conclusions

The water assisted mixing method was found to be useful for preparing thermoplastic composites consisting of EAA copolymer reinforced with up to 70 vol% cellulose nanofibrils, departing from a 10% suspension of cellulose nanofibrils and a water-dispersed EAA copolymer. After drying and compression moulding, the cellulose fibrils were well dispersed in the EAA matrix as seen with X-ray computer tomography. Some voids were, however, found in the interface between cellulose and the EAA matrix, suggesting that the adhesion between the matrix and cellulose surfaces could be improved.

The addition of cellulose nanofibrils improved the mechanical properties of the composites substantially. The tensile stiffness increased with the cellulose nanofibril content, from 0.3 GPa for the neat EAA to 6 GPa for the composite with 70 vol% cellulose nanofibrils. A maximum in effective stiffness was found at 30 vol% nanofibrils content. The stiffness of the composite reinforced with a cellulose pulp was 5 GPa for 70 vol% pulp in EAA. The tensile strength of the composites reinforced with cellulose nanofibrils increased from 24 MPa for the neat EAA to 83 MPa for composite with 70% cellulose nanofibrils. The strength of composites reinforced with cellulose pulp was 55 MPa for the 70 vol% pulp composite, which was substantially lower than for the composites with cellulose nanofibrils. The elongation at break of the composites with 70 vol% cellulose nanofibrils was about 5%, being twice as high as for the corresponding cellulose pulp composite.