3.2. Quasi-Static Nano-Indentation
In an initial experimental investigation aimed at surface characterization of nanocomposite materials, a nanoscale indentation following the Oliver-Pharr method [
43] is adopted to evaluate the hardness (H) and reduced elastic modulus (Er) of the samples as well as the respective contact depth and the coefficient of friction (COF) at scratch. In particular, the nanoindentation hardness is determined by dividing the peak load (P
max) by the contact area of the indenter (A), whereas the process for calculating the reduced modulus of elasticity entails fitting the unloading section of the load-displacement curve with a power-law function to ascertain stiffness. Subsequently, Er is estimated by establishing the mathematical relationship between unloading stiffness and the projected contact area under load. Notably, this calculation considers elastic deformations in the specimen and the indenter. Indeed, reduced modulus (E
r) in the nanoindentation test is determined using the Sneddon’s formula [
45]:
where E and ν represent the elasticity modulus and Poisson's ratio of the test specimen, while E
i and ν
i denote the elasticity modulus and Poisson's ratio of the indenter tip.
The experimental results of aforementioned mechanical properties are graphically reported in
Figure 8 a), b), c) and d) respectively, whereas all numerical results, including statistical information, are summarized in
Table S2 of section E.S.M.
The findings presented in this study reveal that the hardness of sample EP exhibits an approximately 36.7% increase when compared to samples EP-R-160-DBA and EP-R-160-T (see
Figure 8a). The unfunctionalized unfilled epoxy resin (EP sample) tends to be harder than rubber-functionalized epoxy due to its higher cross-linking density and the absence of a rubber phase, which reduces the matrix's rigidity. The addition of finely well-dispersed globular domains of rubber acting as a plasticizer reduces the overall hardness.
Rubber-functionalized epoxy with self-healing fillers shows a phase composition with nanodispersions domains (as highlighted by the SEM investigation) causing localized deformation and further reducing hardness.
On the other hand, the same observations just reported can justify the results regarding the reduced modulus (see
Figure 8 b) of the sample EP, which is about 7% higher compared with the other samples. As a consequence of this higher stiffness of the not functionalized epoxy matrix (EP sample) with respect to the remaining samples analyzed, and as evident from the
Figure 8 c), the contact depth for it is significantly lower (876.1 nm) in comparison with the other measured valued. Finally, regarding the coefficient of friction (COF) during scratch testing reported in
Figure 8 d), sample EP-R-160-T demonstrates the highest coefficient, showing an enhancement of approximately 11.4% and 8.2% compared to sample EP and EP-R-160, respectively.
This enhancement in scratch resistance can be attributed to the incorporation of rubber domains, which alter the material's mechanical properties and enhances its ability to withstand external forces in lateral displacements.
In particular, in the rubber-functionalized epoxy, the dispersed rubber phase acts as energy-dissipating regions upon scratching. These rubber domains can absorb and distribute the applied force, reducing the severity of scratches and preventing crack propagation. To complete this analysis,
Figure 9 shows the time evolution of the of friction during scratching (defined as the ratio between the lateral force LF and normal force NF) for these two mechanically opposing formulations considered in this study, i.e., the unmodified resin (EP) and the EP-R-160-T sample in a) and b), respectively.
It is worth noting that, in line with the previously mentioned considerations, the EP-R-160-T sample exhibited greater resistance over time to the strain caused by the indenter scratching movement than the EP sample.
Figure 10 displays average load-displacement curves (from 49 tests), depicting different levels of indenter penetration within the composite samples (see
Figure 10 b).
Applying a trapezoidal load function with a maximum force of 8000 µN (see
Figure 10a) led to a deeper displacement exceeding 700 nm within the material, thereby enabling a more comprehensive assessment of the combined effects on elastic modulus and hardness. The reduced maximum displacement observed in the curves corresponding to the EP sample indicates greater resistance to indenter movement within the material.
This suggests, once again, that this composite has higher mechanical properties, particularly in terms of stiffness, compared to those containing rubber domains (the remaining composites investigated).
It is interesting to note how the results from nanomechanical characterization align with the experimental findings of Dynamic Mechanical Analysis (DMA) conducted in our previous study [
29] and hereafter referenced for clarity and ease of comparison (see
Figure 11 which reports the Tan δ and storage modulus vs temperature for unmodified (EP) and rubber functionalized (EP-R-160) samples, in a) and b), respectively.
The sample of unmodified resin (EP), which exhibited higher values of hardness and reduced modulus compared to those functionalized with rubber, showed the highest peak in the mechanical spectrum, related to the glass transition (i.e., Tg), which is centered at 200°C, whereas for the rubber-functionalized sample (EP-R-160), it is appears around 132°C. Furthermore, it is noticeable that the introduction of rubber also affects the storage modulus. The primary decrease in the storage modulus occurs at lower temperature values for EP-R-160 as a result of incorporating the functionalized precursor into the resin. This is a very relevant result. In fact, a big challenge addressed in the formulation of these self-healing resins was to improve the resin's dynamic properties and reduce the matrix's rigidity, acting on its phase composition. Therefore, the design of a material with small domains of the polymer (rubber phase) at higher mobility, finely interpenetrated in the resin, has allowed higher mobility of the chains around the elastomeric domains, and a higher efficiency in the autorepair action. DMA analysis evidences a reduction of the rigidity of the matrix (at the macroscopic level), as evidenced before by the results of “Quasi-static nano-indentation” tests. In
Figure 11 it is clearly visible that the peak of Tan δ vs. temperature, for the sample containing rubber phase, opens around 60°C. This peak is shifted to a lower temperature range with respect to the sample without rubber phase. It involves a temperature range from 60 °C to 180 °C (with a max around 132 °C), where the sample without the rubber phase shows a transition from 140°C to 250°C.
3.4. In-Situ Scanning Probe Microscopy (SPM) Inspection
In-situ scanning probe microscopy (SPM) integrates nanomechanical testing with high-resolution imaging using a single probe, improving accuracy, repeatability, and efficiency. Unlike ex-situ methods, it avoids the need to relocate the sample, thereby saving time and enhancing positioning accuracy.
Concerning the foregoing,
Figure 14 (a-e) shows 2D SPM images (left part) and 3D SPM images (right part) of the XPM nanoindentation test trace made on the surface of all samples taken into consideration in our work.
Observing mainly
Figure 14 a) related to the not-modified epoxy resin (EP sample), it is possible to confirm previous considerations about the surface smoothness and its quality. Rubber-functionalized epoxy resin displays visible imperfections and irregularities and a not-uniform surface finish. These observations support the statements reported in the previous section concerning the XPM analysis, which has shown better uniform surface color maps for the not-modified epoxy, compared to the other samples.
Another peculiar aspect which provides information on the nano-scratch behavior of composite materials is the trace profile of the groove left by the indenter tip, as depicted in 2D view (up part) and 3D view (down part) of
Figure 15 for the sample EP, EP-R-160 DBA and EP-R-160 - T in a), b) and c), respectively.
The scratch appears noticeably more linear in the case of the unmodified resin (EP sample, see
Figure 15 a)). Conversely, for the other two formulations compared (EP-R-160 DBA, see
Figure 15 b) and EP-R-160-T, see
Figure 15 c)), the path of the indenter tip appears hindered by the rubbery domains, resulting in jagged and less defined scratch outlines There is also a perceptible difference in the width of the incision, which can be correlated with the deeper notch made by the probe.
These graphical results agree with the lower scratch resistance of these rubber functionalized-based samples, which was discussed in terms of the coefficient of friction in the previous sub-section 3.2. Quasi-static nano-indentation.
Furthermore, from the 3D views of
Figure 15 a), b) and c), it is possible to observe the depth of the tip’s incision, which in the case of EP sample is lower compared to the other two formulations (EP-R-160-DBA and EP-R-160-T) considered for this comparison via images. This result is consistent with the contact depth values already observed and discussed for these specimens.
3.5. Multhyphisics Simulation Study: Results
Finite Element (FE) simulation is a well-recognized technique for analyzing the material's mechanical behavior. Therefore, an experimental-simulation approach is adopted to compare simulation and experimental data. Once the FE model had been validated, it was used to investigate further mechanical properties.
The simulations are conducted focusing exclusively on two materials, EP and EP-R-160-T, which, based on the observed experimental results, are the two most mechanically distinct formulations. The numerical analysis begins with a preliminary investigation over the entire time interval in which the load is applied, aiming to identify the most significant time points for further numerical insights.
Figure 16 reports the z-axis displacement versus the entire time interval (loading, holding, unloading) for EP sample in a) and EP-R-160-T sample in b). The curves provide detailed insights into the material's mechanical behavior during the nanoindentation test. At the beginning of this characterization, as the indenter applies force, the penetration depth increases steadily in linear mode. This initial loading stage is pivotal for comprehending the material's reaction to mechanical strain. Subsequently, there is a holding period during which the load is maintained. This allows for the observation of any viscoelastic behavior or time-dependent deformation of the material. Finally, during the unloading phase, the applied load is removed, and the indenter retreats from the material. The depth rate profile in this phase reveals significant details about the material's elastic recovery and residual indentation depth.
Henceforth, based on the mechanical response over time observed in the previous figure, only some specific time points will be considered as references. More specifically, as indicated in
Figure 17 a), three instances (0.025 s, 0.05 s, and 0.075 s) were selected during the loading phase, while during the holding interval, the time instance t = 0.13 s was considered, at which the maximum indentation depth was recorded. The corresponding load values are explicitly reported in the figure. Obviously, for the time instances without load (t = 0 and t = 0.3 s), as shown in
Figure 17 b), the samples are at rest with the indenter tip just tangent to their upper surfaces.
Figure 18 reports the simulation findings on the z-axis displacement observed for the EP (left part) and EP-R-160-T (right part) for the three-time instances selected during the holding phase, i.e. 0.025s, 0.5s and 0.075 in a), b) and c), respectively. These views allow us to highlight that, consistent with the results of the previous Figure and regardless of the sample, the indentation tip, or equivalently the contact depth, progressively increases over time. In particular, this depth is greater (please refer to the color bar or to the quote lines drawn as guidelines for the eyes) for the EP-R-160-T sample due to its lower hardness and reduced modulus, which have experimentally determined and already been discussed as results.
The time instant t=0.13 s is considered in
Figure 19 a) and b) for the EP and EP-R-160-T samples, respectively. At this instant, the maximum contact depth is recorded, measuring -877 nm and -1010 nm for the two samples, which are close to the experimentally obtained results. This particular cross-sectional view allows for a better exploration of the nanoindentation profile. Specifically, the contour lines highlight the greater depth for the EP-R-160-T sample compared to the EP reference according to their hardness values. For the unmodified resin (EP), the z-displacement during nanoindentation is consistently lower. This is because it exhibits higher hardness and stiffness. As the indenter penetrates the surface, the material resists deformation, resulting in a smaller indentation depth. This behavior indicates that the material can withstand higher loads without undergoing significant plastic deformation. In contrast, the sample EP-R-160-T tends to exhibit a greater z-displacement under the same indentation load, reflecting its softer nature. This increased z-displacement is fully compatible with the lower hardness compared to the reference EP. These differences are essential for tailoring materials to specific applications where a balance between hardness and toughness is required.
The z-axis displacement (D) versus time (depth rate) describes the change in penetration depth of the indenter over time. Once again, this parameter is related to the material's mechanical properties and, in particular, to its hardness.
Figure 20 reports the results concerning the depth rate of the two considered samples evaluated as the slope of the fitting curve (linear fit, R
2 strictly close to 1) of the displacements of the indenter recording during the loading phase. The EP sample, characterized by a higher hardness compared to the EP-R-10-T, exhibits greater resistance to indenter penetration, and therefore its depth rate is notably lower: the calculated slopes result in 172.7 and 203.8, respectively, for the two samples.
In a nanomechanical test of a composite material, the z-axis displacement versus thickness can provide information on the deformation behavior of the material under an applied load. More in detail, as the thickness of the material varies, the z-axis displacement measures the vertical displacement of the material's surface in response to the applied force. This relationship is strictly close to the mechanical properties of materials such as its hardness, elasticity, and resistance to deformation.
Figure 21 shows, for the specific time instants considered in this numerical analysis, the z-axis displacement vs the thickness (alongside the segment highlighted in the inset of the same figure) for EP and EP-R-160-T samples, in a) and b), respectively.
As expected, EP, being a stiffer material, exhibits less displacement for a given applied load, resulting in a smaller z-axis displacement versus thickness curve. Conversely, the EP-R-160-T exhibits greater displacement, resulting in a larger curve since it is a softer material.
Figure 22 shows a 3D view of the indentation imprint at the maximum contact depth (
hi), i.e., 877nm and 1010 nm, for the EP sample and the EP-R-160-T sample in a) and b), respectively, whereas the corresponding 2D top views are shown in
Figure 22 c) and d). As expected and evident from the analysis of the figures, the projected areas are different according to the relationship A
proj=
f(h
i). With the greater contact depth reached by the tip in the EP-R-160-T sample, its area is significantly larger than that of the EP sample.
Using equation 4, the numerically estimated area values are Aproj (EP) =18.89 µm2 and Aproj (EP-R-160-T) = 25.05 µm2, respectively.
Once these values were determined and with the maximum applied load (8000 μN) known, the hardness (H) of the materials was calculated using Equation 5. The resulting values were 0.42 GPa for EP and 0.32 GPa for EP-R-160-T, closely matching the experimentally obtained values. These numerical results are summarized and compared with the experimental data in
Table 1, showing the percentage variation between the two data sets.
In nanoindentation investigation, when comparing different composite materials, the displacement magnitude provides further insights into their mechanical properties due to their hardness, elasticity, plasticity, and contact depth differences.
Figure 23 shows the deformation contours, detected at the instant t=0.13 s corresponding to the maximum load of 8000 μN, for the EP sample in a) and EP-R-160-T in b). The 3- dimensional cross-sectional views are chosen to provide more detail of the resulting indentation profile.
Regardless of the sample, the deformations are highly localized at the contact area between the material/indenter tip interface and along its periphery. It progressively extends into the bulk of the samples, primarily along the z-axis direction, where the load is applied. The maximum values recorded in both cases (993 nm for EP and 869 nm for EPR-160-T, respectively) are consistent with the previously presented and discussed contact depth data and hardness values. For the EP sample, being more rigid and thus more prone to crack propagation, it is evident that the deformations extend more significantly in other spatial directions compared to the EP-R-160-T sample. The latter, being relatively softer, better absorbs and dissipates the mechanical energy, reducing the extent of deformations.
Once the model is validated with experimental data and the reliability of the results is confirmed, it can be used to investigate additional mechanical properties that have not yet been here observed experimentally.
For example, in this study, the Von Mises stress is numerically investigated since it is a critical measure in engineering and materials science, particularly for predicting when a material will yield under complex loading conditions. It is especially relevant for composite materials, which combine two or more distinct phases to achieve superior mechanical properties compared to the individual components.
Figure 24 depicts the variation of Von Mises stress (average values) across the total time span of force application (loading, holding, and unloading phases) for both composites, EP and EP-R-160-T, analyzed over the entire domain and on their upper surfaces (see the inserts of the same figure), in a) and b), respectively.
From the analysis of these results, it is evident that, regardless of the exploration area, whether volumetric domain or surface, the Von Mises stress remains higher for the EP sample than the EP-R-160-T sample. The maximum stress peak for both samples is reached at t=0.13 s, corresponding to the maximum penetration depth. Therefore, this particular time instant will be considered for the next investigations. These maximum stress peaks are 49.9 [µN/µm²] and 41.8 [µN/µm²] over the domain for EP and EP-R-160-T, respectively, and 79.1 [µN/µm²] and 69.3 [µN/µm²] on their upper surfaces.
Figure 25, panels a) and b), respectively, presents a 3D view of the two samples properly sectioned to better visualize the spatial distribution of Von Mises stress.
The pictures highlight that the main stress values are highly localized at the interface indenter probe/contact surface and near its periphery to extend into the surrounding space progressively. In accordance with the results and considerations already reported in the previous
Figure 24 a), the higher stress intensity for the EP sample is visually noticeable compared to that of the EP-R-160-T.
Finally, the Von Mise stress profiles detected on the upper surfaces of the materials are further analyzed graphically, with 3D views and top views, in
Figure 26a) and b) correspondingly for the EP and EP-R-160-T samples. These additional graphical representations again highlight and quantify (see the respective color bars) how the higher Von Mises stress is concentrated in the indentation area.
The traced contour lines provide information on the Von Mises stress distribution on the material surfaces: the appreciable reduced extent exhibited by the EP-R-160-T sample confirms its facilitated energy dissipation and improved ability to absorb impacts, contributing to the overall toughness and resistance to mechanical stresses.