1 Introduction

One of the most important aspects that concerns all the structural materials and, more in particular composite structures, is constituted by the behavioural uncertainties associated with impacts with foreign objects. Indeed, because of their weak resistance to through-the-thickness solicitations, these materials are susceptible to delamination damages under impact loads that can compromise their integrity and lead to catastrophic failures. Based on this premise, multifunctional composite systems play a fundamental role in the development of new generation impact resistant materials due to the possibility to modify their specific properties by adding or modifying specific phases within the laminate structure. Indeed, by following this approach, it is possible to enhance the impact properties of a traditional laminate by manufacturing a Shear Thickening Fluid (STF) that can subsequently be included within the laminate sequence as an active layer for impact damage suppression.

STFs are non-Newtonian fluids with remarkable properties that behave like a liquid when there are no forces applied and turn in a very stiff solid-like structure in presence of high shear rates [1]. Since this liquid-solid transition is the result of a rapid increase of viscosity, these highly nonlinear fluids have been exploited for the design of damping systems and shock absorbers.

Liquid body armour based on the inclusion of STF have been studied intensively during the last decade [25], and several commercial products (mostly reinforced soft protections based generally on aramid fibres) are available on the market. However, very few works have been done on the design of stiff composite materials with embedded STF (without the use of aramid woven fabric). In the work carried out by Fisher et al., an STF was sandwiched between polyvinyl chloride beams, showing good results in damping the vibrations of the entire structure [6], while Soutrenon et al. realised a structural damper by embedding the STF within an open-cell foam encapsulated in silicone [7]. Another multifunctional system patented by Picken [8] is based on the preparation of an STF using a supramolecular polymer (polymers made of monomers held together by noncovalent interactions [9]) that is then used to impregnate the fibres and produce a composite material. The main focus of this work is the development a engineered STF that is can be embedded directly as an active layer within the fibrous reinforcement of a composite laminate in order to manufacture stiff CFRP parts characterised by higher levels of energy absorption due to the thickening of the solution during impacts with foreign objects. The technique does not imply any modifications of the thermoset resin used for the composite manufacturing, therefore the procedure can be carried out using commercial thermoset/carbon fibres prepregs.

2 Theoretical Aspects

Generally, a non-Newtonian fluid can be defined as a material whose viscosity is dependent on the shear rate or shear rate history. However, included in the group there are some fluids in which it is possible to observe normal stress-differences or other non-Newtonian behaviour even if they are characterised by shear-independent viscosity.

Figure 1 shows the rheology of different kinds of Non Newtonian fluids, according to the different relationships between shear stress and shear rate.

Fig. 1
figure 1

Rheology of non-Newtonian fluids

A Shear Thickening Fluid (see red curve in the figure above) is a non-Newtonian fluid obtained by dispersing spherical small particles (dispersed phase) within a liquid phase (continuous phase). Increasing the particles’ concentration beyond dilute limit leads to the activation of particle-particle interactions and the creation of long range order, resulting in non-Newtonian effects such as shear thinning or shear thickening. These effects can be explained with transitions from an ordered to a disordered state that occur when high stresses are applied [10]. Indeed, when this fluid is subjected to a relatively low shear rate, the repulsive interactions between particles are able to create an ordered and layered structure that consents a more effective particles flow through the liquid phase, reducing the viscosity of the solution. At this point, if the shear rate is increased above a critical point (\( \dot{\gamma} \) c), the particle-particle interactions are overcome by shear forces and the particles are squeezed together, moving the fluid out from its equilibrium state and resulting in a disordered and more complex structure. In particular, shear thickening occurs when the repulsive Brownian forces are overcome by attractive hydrodynamic shear forces, leading to the formation and percolation of shear induced transient aggregates called “hydroclusters” that thickens the fluid [11].

As the hydroclusters are made of particles momentarily compressed together to form agglomerates, the particles’ ability to flow around each other is hindered and requires more energy, leading to an abrupt increase in viscosity [12, 13]. When the shear stress is removed, the hydroclusters spread apart and the particles are separated, forming again a stable suspension (Fig. 2).

Fig. 2
figure 2

Schematics of Shear Thinning and Shear Thickening Effect

In order to observe shear-thickening behaviour in a suspension, there are three important factors to consider:

  • Increasing the concentration of particles within the fluid will lead to an enhancement of the thickening effect, therefore the weight fraction has to be evaluated according to the system viscosity requirement [14]. For our purposes the concentration required is around 20–25 % wt. as it gives good results in terms of impact resistance and workability for composite applications. Further increase of the particles concentration leads to a very high level of initial viscosity of the solution (which is also dependent on the molecular weight of the continuous phase), decreasing its workability during the embedding process within the laminate. This problem can be partially solved by increasing the temperature during the manufacturing process. Moreover, volume fraction is strictly connected with particles dimensions as it can change the electrostatic interactions, leading to different forces interaction in the liquid/particles system.

  • Suspension must not flocculate (destabilisation of the suspension): particles have to be neutral or repel one another by electrostatic, steric or entropic interactions. To achieve this, the continuous phase must induce steric or electrosteric stabilisation in order to enhance the colloid stability. Several carriers can be used for the preparation of the STF, such as ethylene glycol [15, 16], water [17] and other organic solvents, however the most used one for structural application is polyethylene glycol due to its tendency to form very stable suspension [4, 18] coupled with its biocompatibility.

  • Particle volume fraction at which shear thickening behaviour is first observed is affected by the particle ratio. Studies present in literature report thickening behaviour observed for particles with a dimension range from 1 nm to 5000 nm [19]. Another important effect is given from the shape of the dispersed phase, however, in case of irregular particles, it is more difficult to understand the energy dispersion mechanism due to the superimposition of the thickening effect and the reciprocal friction between the inclusions [20].

Based on these considerations, the STF realised for this study forms a complex system which is characterised by a double change in its microstructure with the increase of the applied shear rate. Indeed, at low shear rate, the particles are rearranged from the equilibrium and become more organised in the flow, decreasing the fluid viscosity (shear thinning). A further increase of the shear rate leads to the formation of the hydroclusters and, consequently, to the thickening effect.

As for the particles used, there are several examples in literature of the different materials that can be employed, such as CaCO3 [20], PMMA [21, 22], PVC [23], PS and SiO2 [6, 13, 20]. For our purposes, different solutions were produced using two types of commercially available SiO2 with average diameters of 7 and 14 nm (Sigma Aldrich). Good thickening effects have been obtained for particles with a nominal diameter of 200 nm.

3 Experimental: Solution Preparing

The components needed for manufacturing an STF are nanosized particles and a liquid carrier. For this work, several solutions were produced by mixing an increasing percentage of silica (SiO2) particles (from 10 to 25 % wt.) with PolyEthylene Glycol (HO-CH2-(CH2-O-CH2-)n-CH2-OH).

To ensure a good distribution of the silica within the carrier and to avoid the formation of agglomerates, an ultrasonic cavitation mixing method was used. Small quantities of silica particles were gradually added to the PEG and stirred in manually. Once a viscous gel was formed, the ultrasonic probe was immersed in gel and activated for 10 min. High amplitude compressive waves generate localised low-pressure zones in the carrier, ensuring homogenous mixing of silica and PEG. At the end of the process the solution were stored under vacuum for 12 h in order to remove air bubbles and further improve the dispersibility.

4 Rheological Behaviour of STF Solutions

By following the approach of Stickel et al. [24] it is possible to define the viscosity of a suspension (η) as a general function of multiple parameters which are based on the characteristics of the dispersed phase (the particles, indicated with the subscript dp) and the continuous phase (the suspending media, indicated with the subscript cp):

$$ \eta =f\left({r}_{dp},{\rho}_{dp},{C}_{dp},{\rho}_{cp},\eta, {k}_BT,\dot{\gamma},t\right) $$
(1)

In particular, r dp , ρ dp and C dp represent the hydrodynamic radius (which is equal to the particle radius in case of spherical particles), density and number concentration of the particles, while ρ cp  , η cp refer to density and viscosity of the suspending media. The additional terms that affects the viscosity are the thermal energy (k B T, where k B is the Boltzmann’s constant and T is the absolute temperature), the shear rate (\( \dot{\gamma} \)) and time (t). These variables can be reduced from eight to five by grouping the properties of the two phases in dimensionless quantities, leading to the definition of the following parameters:

$$ {\eta}_r=\frac{\eta }{\eta_{cp}} $$
$$ {\rho}_r\frac{\rho_{dp}}{\rho_{cp}} $$
$$ \phi =\left(\frac{3}{4}\pi {r}_{dp}^3\right){C}_{dp} $$
$$ {t}_r=\frac{t{k}_BT}{\eta_{cp}{r}_{dp}^3} $$

This reduces Eq. 1 to:

$$ {\eta}_r=f\left(\phi, {\rho}_r,{k}_BT,\dot{\gamma},{t}_r\right) $$
(2)

where ϕ represents the volume fractions of the particles within the solution.

In order to further reduce the variables that influence the viscosity it is possible to introduce the Reynolds number for a spherical particle immersed in a fluid:

$$ {Re}_{\dot{\gamma}}=\frac{\rho_{cp}{v}_{dp}{r}_{dp}}{\eta_{cp}}=\frac{\rho_{cp}{r}_{dp}^2\dot{\gamma}}{\eta_{cp}} $$

Considering that the particles dimensions are below one micron, it is possible to assume very low values of the Reynolds number for the solution under study. As a consequence, if we consider a single particle of the dispersed phase, it is possible to correlate its diffusion through the continuous phase with the thermal energy and its mobility μ (the ratio of the particle’s terminal drift velocity to an applied force) by following the Stokes-Einstein-Sutherland fluctuation-dissipation relationship, as discussed by Wagner et al. [10]:

$$ D=\mu {k}_BT=\frac{k_BT}{6\pi {\eta}_{cp}{r}_{dp}} $$

where the Stokes’ law was used to express the mobility as the inverse of the drag coefficient (\( \mu =\frac{1}{\zeta }=1/6\pi {\eta}_{cp}{r}_{dp}\Big) \).

From this expression, it is possible to estimate the time required for a particle to diffuse through a distance equal to its hydrodynamic radius as \( {r}_{dp}^2/ \) D 0, from which high and low shear rates \( \dot{\gamma} \)can be defined.

The diffusion rate D measures the interaction between the particles and the solvent molecules and can be related to the stress rate \( \dot{\tau} \) or the the shear rate \( \dot{\gamma} \) introducing the dimensionless Peclét number (Pe):

$$ P{e}_{\dot{\gamma}}=\frac{\dot{\tau}{r}_{dp}^36\pi {\eta}_{cp}}{k_BT}=\frac{\dot{\gamma}{r}_{dp}^3}{D_0} $$

Using the expressions of the Reynolds and Peclét numbers, it is possible to rewrite the expression of the viscosity as:

$$ {\eta}_r=f\left(\phi, {\rho}_r,R{e}_{\dot{\gamma}},P{e}_{\dot{\gamma}},{t}_r\right) $$

In addition, by considering that the particles of the dispersed phase are neutrally buoyant, it is possible to neglect both ρ r and t r . Leading to:

$$ {\eta}_r=f\left(\phi, R{e}_{\dot{\gamma}},P{e}_{\dot{\gamma}}\right) $$
(3)

Since it correlates the shear rate with the rheology of a solution, the Peclét number can be used to analyse the behaviour of Non-Newtonian fluids. Indeed, while for low shear rates (\( {Pe}_{\dot{\gamma}}<1 \)) the Brownian motion has the time to restore the fluid structure keeping it in an equilibrium state, when the shear rate increases, the fluid deformation operates in a shorter time scale than the Brownian motion, moving the solution out of its equilibrium state.

Following Eq. 3, it appears clear that in order to investigate the properties of the STF, it is important to control two different parameters: the volume fraction of particles (ϕ) and the shear rate (\( \dot{\gamma}\Big) \). Based on these considerations, several solutions were prepared with an increasing percentage of silica and the flow properties were measured using a Bohlin C-VOR Rheometer in parallel plate geometry (PP40) increasing the shear rate.

The results obtained from the rheometric analyses are shown in Fig. 3 and are quite revealing in several ways for all the tested solutions. First, it is clearly shown from the curves that both thinning and thickening effects are present for highly concentrated solutions, while for lower concentration the presence of secondary flows prevents accurate measurements of the weak thickening effect. Indeed, after an initial equilibrium state in which the fluid behaves almost like a Newtonian fluid, it is possible to observe a slight trend of decreasing viscosity for relatively low shear rates, while, when higher rates of shear are applied, the viscosity reaches a critical point in which the magnitudes of repulsive inter-particles forces are balanced by shear forces. When the shear stress exceeds this critical level, the shear forces that push the particles together become larger than the repulsive interactions so that the particles are moved from their equilibrium state, disrupting the ordered structure and leading to an abrupt increase of viscosity that thickens the fluid. Since the concentration of hydroclusters within the continuous phase is strictly related on the particles concentration [24], increasing the silica volume fraction within the STF makes the nonlinear effects more evident, enhancing thinning/thickening behaviours (see Fig. 4).

Fig. 3
figure 3

Rheological analysis for STF at different concentrations: a) STF/15; b) STF/20; c) STF/25; d) comparison between the different curves

Fig. 4
figure 4

viscosity increase with increasing silica concentration

Initial viscosity of the solution depends strongly also on the molecular weight of the carrier. In this case the carrier used was a PEG 200 g mol−1. Results in terms of viscosity change due to thickening effects (two orders of magnitude increase) are not dissimilar to the ones found in literature [25].

Another important result can be observed from Fig. 5, which illustrates the existence of a strong relationship between volume fraction and the critical shear rate. Indeed, for higher particle concentrations, \( \dot{\gamma} \) c moves to lower values due to the reduced intra-particles distance, anticipating the transition between thinning and thickening effects, as reported by Maranzano et al. [14].

Fig. 5
figure 5

Behaviour of critical shear rate with increasing particles concentration

5 Composite Mechanical Characterisation

The structural composite material was manufactured sandwiching a thin layer of STF between several layers of Carbon Fabric prepregs (see Fig. 6). The amount of fluid used in every sample has been kept constant for all the tests undertaken (10 g for 300x150mm plates), therefore we will refer to “STF percentage” as for the amount of silica within the fluid. Moreover, because of the small amount of STF involved in the hybridisation process, the total weight of the final laminate is not affected. Samples were prepared via Vacuum Assisted stratification followed by Autoclave consolidation in order to activate the curing reaction of the epoxy resin, and cut in 100x150mm plates for impact testing (Fig. 7).

Fig. 6
figure 6

Schematics of STF/CFRP composite preparation

Fig. 7
figure 7

manufacturing procedure of STF/CFRP composite samples

5.1 Damping Analysis

In order to fully understand how the material would behave when subjected to an impact, experiments were conducted to evaluate the energy absorption of the STF/CFRP composite and a comparison was made with a traditional laminate, obtained with the same prepreg and the same lay-up sequence. A first series of test was undertaken to analyse the behaviour of the samples when subjected to a forced vibration using n Macro Fibre Composite (MFC) piezoelectric actuator to input a chirp signal over a 0-10 kHz range of frequencies. Fast Fourier Transform (FFT) was performed in order to obtain the frequency content of the acquired vibration signals and evaluate the resonance frequencies of the samples.

The resonance frequencies of a structure are proportional to its stiffness, according to the relation \( {f}_i\propto \sqrt{k/m} \), therefore it is possible to evaluate this mechanical characteristic by comparing the relative shift of frequencies peak patterns of hybridised samples and traditional unreinforced laminates [26]. Moreover, since the relative amplitude of the resonance peaks is a measurement of the ability of the material to store energy when subjected to an external load, it is possible to have an estimation of the damping ability of the structure by analysing its variation between the hybridised samples and the control ones.

Results from the tests are illustrated in Fig. 8 which clearly shows a visible shift to higher frequencies between the control sample and the hybridised one, indicating an increase in stiffness for the STF sample. Similar results have been observed for CNT reinforced composite structures [27]. In addition, the shifting of the amplitude of the entire frequency spectrum to lower levels suggests an enhancement of the damping behaviour of the structure due to the dynamic change in viscosity of the embedded STF.

Fig. 8
figure 8

Forced vibration signals for STF and Neat samples

A more in depth analysis on the damping behaviour of the STF based composites was carried out with free vibration tests in order to quantify the dissipation effect of the hybrid phase. Samples were tested in single cantilever beam geometry, clamping them to an optical bench and vibrations signals were collected using piezoelectric sensors.

Logarithmic decrement δ was evaluated from the acquired signals, following Eq. 4:

$$ \delta =\frac{1}{n}\mathit{\ln}\frac{x(t)}{x\left(t+nT\right)} $$
(4)

Where x(t) is the amplitude at time t, and x(t + nT) is the amplitude calculated for a peak n periods away. Subsequently, results obtained can be used to evaluate the damping ratio, by using Eq. 5:

$$ \zeta =\frac{1}{\sqrt{1+{\left(\frac{2\pi }{\delta}\right)}^2}} $$
(5)

Results from this series of tests demonstrate the effectiveness of the STF within the CFRP structure as structural damper (see Fig. 9). The trends of both logarithmic decrement and damping ratio (+50 % for both δ and ζ), clearly indicate how the inclusion of STF changes the damping properties of the material. This difference can be explained by the dynamic microstructure modification of the non-Newtonian fluid in the composite structure. Indeed, when the composite is subjected to an impact, part of the stress is transferred to the embedded STF and it activates the transient formation of hydroclusters. The abrupt increase in viscosity enhances the energy dissipation rate for the entire structure.

Fig. 9
figure 9

Free vibration signals and results for STF and unreinforced samples

5.2 Impact Testing

In order to fully characterise the impact properties of the hybridised composite, several test were conducted on STF reinforced samples and a comparison was made with traditional laminates obtained with the same prepreg and the same layup procedure. Two different variables have been taken in account during these tests: silica particles dimensions (7 and 14 nm) and particles concentration within the STF (from 10 to 25 %).

Tests were conducted using a drop tower rig with an impactor mass of 12,684 kg and a semi-spherical head geometry (20 mm diameter). Impacts were undertaken at 6, 20 and 40 J with velocities of 0.8, 1.7 and 2.5 m/s.

Figure 10 shows the curves obtained with the same impact energy for samples with increasing 7 nm particles concentration within the STF (herein referred as STF7/XX where XX stands for the silica concentration),

Fig. 10
figure 10

Example of the force displacement curves from the impact test on 7 nm Silica STF reinforced CFRP: a) 6 J; b) 20 J: c) 40 J

Results are summarised in Fig. 11 in which it is possible to observe that for low velocity impacts (6 J) the peak force stays linear (or within the experimental error) showing no large variations between control and reinforced samples.

Fig. 11
figure 11

Summary of the results for impact on STF7 reinforce samples: a) Force Peak variation at different energies and different silica concentrations in STF; b) Energy absorption variation at different energies and different silica concentrations in STF

For higher velocity impacts the effects of the hybridisation process becomes more consistent and affects the value of both force peak and the energy absorption rate.

Indeed, for the 20 J impacts, the presence of the STF leads to a slight increase (~12 %) in the value of the force peaks, going from 4252 N (unreinforced samples) to 4724 N (STF25 samples), while there is no variation in the energy absorption rate.

The effects of the hybridisation are more evident for the 40 J impacts, where the curves show a dramatic increase in the force peaks, going from 5197 N for the unreinforced samples to a maximum of 6929 N for the STF based ones (~30 %). It is possible to observe a relationship between this increment and the concentration of silica particles within the solution, as the force peaks values are enhanced by 12 % for the STF15 and by an additional 20 % for STF25.

Same trend can be observed for the energy absorption rate during the impact. In this case the effect of the percentage of SiO2 within the STF is even more evident. Indeed as it is possible to observe from Fig. 11b, while for 6 and 20 J impacts the trend is basically linear, when the impact velocity is raised up to 2.5 m/s the presence of the STF increases the energy absorption rate by 38 % for STF7/15 and almost 50 % for STF7/25. These results confirm what observed in the previous paragraphs and they can be explained due to the dynamic behaviour of the STF, which activates the formation of the hydroclusters only when it is subjected to high shear stresses, thickening the fluid and enabling improved energy absorption mechanisms.

In order to understand the role of particles dimensions within the STF, a second series of test was undertaken using 14 nm silica particles for the STF reinforcement.

Figure 12 represents the Force-Displacement curves grouped in terms of energies and in terms of percentage of SiO2 particles (14 nm) within the STF. As it is possible to observe from the results (Fig. 13), in this case the presence of bigger particles leads to an additional enhancement of the hybridisation effects which is particularly visible especially from the energy absorption curves. Indeed, by analysing the energy absorbed during the different impacts it is possible to observe an increase of more than 50 % for the STF14/20, while for the 7 nm particles the percentage required in order to reach this value was 25 % (STF7/20). For lower impact velocities the energy absorption levels stay basically linear as for the previous samples.

Fig. 12
figure 12

Force Displacement curves from the impact test on 14 nm Silica STF reinforced CFRP: a 6 J; b 20 J: c 40 J

Fig. 13
figure 13

Summary of the results for impact on STF14 reinforce samples: a) Force Peak variation at different energies and different silica concentration in STF; b) Energy absorption variation at different energies and different silica concentrations in STF

It is possible to have a clear idea of the hybridisation effects given by the presence of the STF within the CFRP by analysing Fig. 14, in which the nanoreinforced sample is compared with an unreinforced one after the impact tests. As it is possible to see from Fig. 14a and b, the solicitations generated by the 40 J impacts were high enough to reach critical values on the bottom lamina of the traditional laminate, resulting in a total failure of the fibres on the top surface and a complete delamination on the bottom face. On the contrary, the STF based composite kept its structural integrity after the impact showing no superficial failure on the top surface and only a small internal delamination on the bottom one (see Fig. 14a and d).

Fig. 14
figure 14

Post impact image of STF reinforced sample in comparison with a traditional CFRP subjected to the same impact; a Control sample, top surface; b Control sample, bottom surface; c STF reinforced CFRP top surface; d STF reinforced CFRP bottom surface

5.2.1 C-Scan Analysis

The internal delamination extents of the impacted samples were evaluated analysing some of the CFRP/STF samples using a C-Scan system.

C-Scan is a non-destructive technique that uses a piezoelectric probe to send ultrasonic waves through a sample and measures the reflected signal obtained when a barrier (like a defect or the back of the sample) is reached. Analysing the difference between the signal sent and the one collected by the detector it is possible to assess the delamination region within a laminate. Figure 15 represents several ultrasonic images of the control samples in comparison with hybridised laminates. As it is possible to observe, the presence of the nanofluid reduces the extent of the delamination by more than 30 %. Increasing the concentration of silica within the STF leads to an additional 10 % reduction of the extension of the internal damage, confirming the results obtained with the impact tests.

Fig. 15
figure 15

C-Scan analysis of Impacted STF samples and traditional CFRP impacted at 40 J

6 Conclusions

The aim of this paper was to investigate the possibility to embed silica based Shear Thickening Fluids as shock absorber media within structural CFRP parts, in order to increase mechanical and impact properties. Several solutions were prepared using silica particles with different dimensions and mixed with Polyethylene Glycol at different concentrations. The viscosities of the solutions were analysed with rheological measurements and the results showed an evident increase of the viscosity once a critical shear rate \( {\dot{\gamma}}_c \) is reached. Further tests found out that increasing the concentration of particles within the solutions enhances the thickening effect, leading to higher level of final viscosity. Critical shear rate was found inversely proportional to the silica concentration within the STF.

Composite samples were manufactured sandwiching an STF layer between traditional carbon prepregs layers and were tested in free and forced vibration in order to estimate the damping properties. Forced vibration analyses showed that the amplitudes of the natural frequencies of the hybridised samples are decreased and shifted to higher frequencies in comparison with a traditional CFRP, suggesting an increase of stiffness and damping ratio for the hybrid material. Free vibration tests confirmed the previous results showing an increase in both damping ratio and logarithmic decrement of more than 50 % for the STF reinforced samples in comparison with the unreinforced laminates.

In order to understand the effect of the hybridisation on the out-of-plane properties, impact tests were conducted on samples with STF characterised by an increasing concentration of silica and with different dimensions. The results showed that for high velocity impacts (40 J) the presence of the embedded fluid within the CFRP structure is able to enhance the energy absorption rate by more than 50 %. This effect cannot be seen for lower velocity impacts (6 and 20 J) because the hydroclusters activation mechanism requires relatively higher shear solicitations. Raising the concentration of silica or increasing the dimension of the particles within the STF leads to an increase of the energy absorption.

These results have been confirmed by C-Scan analyses, which showed how the inclusion of the STF is able to reduce by more than 30 % the extent of the internal delamination for samples subjected to high velocity impacts. Increasing the concentration of silica within the STF enhances this effect.