Next Article in Journal
Correlation between Friction and Wear in Cylindrical Anchorages Simulated with Wear Machine and Analyzed with Scanning Probe and Electron Microscope
Next Article in Special Issue
Predicting Failure of Additively Manufactured Specimens with Holes
Previous Article in Journal
FRP-Reinforced/Strengthened Concrete: State-of-the-Art Review on Durability and Mechanical Effects
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Materials
Volume 16
Issue 5
10.3390/ma16051992
materials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interfacial Shear Strength of Single-Walled Carbon Nanotubes-Cement Composites from Molecular Dynamics and Finite Element Studies

by
Carlos Talayero
1,*,
Isabel Lado-Touriño
1,
Omar Aït-Salem
2,
Ismael Sánchez Ramos
2,
Alicia Páez-Pavón
1 and
Rosario G. Merodio-Perea
1
1
Engineering Department, School of Architecture, Engineering and Design, Universidad Europea de Madrid, 28670 Villaviciosa de Odón, Spain
2
Hexagon HMI, 28050 Madrid, Spain
*
Author to whom correspondence should be addressed.
Materials 2023, 16(5), 1992; https://doi.org/10.3390/ma16051992
Submission received: 6 February 2023 / Revised: 23 February 2023 / Accepted: 26 February 2023 / Published: 28 February 2023
(This article belongs to the Special Issue Damage Analysis for Composite Materials: Methods, Testing and Evaluation)
Browse Figures
Versions Notes

Abstract

:
Carbon nanotubes (CNTs) are nanometer-sized structures that can be used to reinforce cement matrices. The extent to which the mechanical properties are improved depends on the interfacial characteristics of the resulting materials, that is, on the interactions established between the CNTs and the cement. The experimental characterization of these interfaces is still impeded by technical limitations. The use of simulation methods has a great potential to give information about systems lacking experimental information. In this work, molecular dynamics (MD) and molecular mechanics (MM) were used in conjunction with finite element simulations to study the interfacial shear strength (ISS) of a structure formed by a pristine single-walled CNT (SWCNT) inserted in a tobermorite crystal. The results show that, for a constant SWCNT length, ISS values increase when the SWCNT radius increases, while for a constant SWCNT radius, shorter lengths enhance ISS values.

1. Introduction

Cement and cement composites are the most used materials in the building sector due to their high compressive strength, low prices and mature production technology. However, cement also has low tensile strength, good permeability to several substances and easily cracks. To improve these weaknesses, different types of reinforcements are added to cement. Carbon nanotubes (CNTs) are one of the most promising candidates due to their excellent mechanical properties. The number of research studies on cement-CNTs composites has steadily grown in recent years [1,2,3,4,5]. CNTs have elastic moduli greater than 1 TPa and tensile strengths as high as 63 GPa [6,7]. When CNTs are homogeneously dispersed within the cement matrix, density [8], permeability [9], mechanical properties [10] and durability [11] of the composites are improved. The effect of CNTs is particularly noticeable on mechanical properties, as they can delay cracks propagation.
The CNTs—cement interface is a key factor in mechanical properties as this is typically a weak area and determines the material fracture behavior [12]. Diverse simulation methods, such as coarse-grained and all-atoms molecular dynamics [13], ref. [14] show that when CNTs are incorporated into a cement matrix, the way of fracture is changed and mechanical properties are clearly enhanced.
Interfacial bond strength can be determined from Interfacial Shear Strength (ISS) measurements, such as fragmentation, pull-out, microdroplet, push-out and push-in tests [15,16,17,18,19,20] and most experimental studies have focused on composites made of CNTs and polymeric matrices due to the difficulties in studying CNT-cement composites. The lack of experimental information can be supplemented by the use of simulation methods such as molecular mechanics (MM) and molecular dynamics (MD). Both techniques provide atomic information about the interface, which cannot be obtained from experiments. ISS can be calculated using a pull-out model [21]. Pulling out a CNT from a matrix provide information about the interface strength. There are numerous molecular modelling studies of ISS of CNTs-reinforced polymer composites [22,23,24,25,26], while results on carbon-based reinforcements are scarce and most of them concern graphene [27,28,29,30].
The Finite Element Method (FEM) is other type of study that can be used to calculate mechanical properties of composite materials. It can be based on different formulations, such as beam models [31], representative volume element FEM analysis [3,32] or nonlinear analysis for dynamic events [33] or when crack analysis is needed [34].
In this work, SWCNT pull-out simulations from a tobermorite crystal were performed to calculate the ISS of several nanocomposites. MD, MM were applied to obtain interfacial energies, which were subsequently used as input in a FE calculation. Pristine SWCNTs with different radii (2.70 to 4.74 Å) and lengths varying from 9.84 Å to 19.68 Å, were inserted into a tobermorite were explained as a function of SWCNT diameter and length and related to non-bonded interactions between the SWCNT and cement. Our work aims to get more insight into the effect that SWCNT geometry plays on ISS, which is a key factor for enhancing mechanical properties of these materials. We think that the synergy achieved by using both atomistic and FE methods is a powerful way of attaining our goals.

2. Models and Methods

In this section, the modelling methods are described.

2.1. Atomistic Model Systems

The periodic simulation cell used for the pull-out atomistic calculations had dimensions of 29.54 × 67.46 × 500 Å3. The cell in the pull-out direction was large enough and included a vacuum layer (Figure 1).
The atomic models built to represent the cement matrix and the composite material are shown in Figure 2. Despite the fact that cementitious systems are not fully crystalline [35], tobermorite 11 Å was employed to model cement by many different authors [36,37,38]. As we are mainly interested in studying interfacial interactions, crystalline tobermoritecan be a good approximation to this interface.
To build the model of the composite material, a SWCNT was inserted into the tobermorite crystal (Figure 2). Different pristine SWCNTs ((CNT (4, 4) = 2.70 Å, CNT (5, 5) = 3.39 Å, CNT (6, 6) = 4.07 Å and CNT (7, 7) = 4.74 Å) with lengths of 9.64 Å, 19.68 Å and 24.60 were used to study the influence of diameter and length on ISS.

2.2. Molecular Mechanics and Molecular Dynamics

The first step in creating a good model for the pull-out process was to relax the simulation lattice by means of a 100 ps simulation in the NPT ensemble at 298 K and atmospheric pressure. The temperature and pressure were controlled by a Nose-Hoover thermostat [39] and a Berendsen barostat [40] respectively. A time step of 1 fs was used during the simulation. Three independent simulations were carried out to obtain average values. After relaxing the structures, they were used to simulate the pull-out process. The size of the cell in x direction was increased and a vacuum layer was included. Before the pull-out study, these new models were subjected to a process of optimization by molecular mechanics (MM). Then, pulling-out was done by displacing the SWCNT along the x axis with increments of 10 Å. Initial and final structures before and after the SWCNT pull-out process are shown in Figure 3. After each pull-out step, the energy of the system was calculated by MM. Materials Studio 7.0 software [41] was used to perform the simulations.
The condensed-phase optimized molecular potential for atomistic simulation studies forcefield COMPASSII [42] was employed to calculate the potential energy of all systems. This forcefield has been successfully applied in numerous MD studies of systems containing CNTs and cement derived materials [28,43,44,45,46]. Coulomb interactions were calculated by the Ewald summation method while an atom-based cutoff method was applied for van der Waal interactions. The cutoff distance for both interactions was 12.5 Å.
The interfacial energy (Einterfacial) was calculated as the energy difference between the fully embedded configuration of the SWCNT and the complete pulled-out configuration (Figure 3). This energy was then used to calculate the ISS by FEM.

2.3. Finite Element Method Model

A Finite Element Model based on Axisymmetric conditions has been built in MSC.Marc software by HEXAGON [47] to represent the composite material. Different FE models have been defined to cover the different combinations of radius and lengths. Same loads, boundary conditions, physical and material properties and analysis methodology are used for each geometrical combination.

2.3.1. Loads & Boundary Conditions

Loads and Boundary Conditions were defined using contact bodies and contact interactions. Five contact bodies were defined (Figure 4): Armchair (a flexible body corresponding to Pristine armchair SWCNTs material), tobermorite (a flexible body corresponding totobermoritematrix material, Base (a rigid body defined by a curve), Tirador (a rigid body defined by a curve) and Sym_Axy (a rigid body defined by a curve). The contact interactions defined between contact bodies represented the load and boundary conditions where “Tirador–Armchair” was a glue contact interaction. This interaction defined the load as Tirador rigid body and it was controlled by its position. Tirador body was moved 1.1 × 10−9 mm in X-Axis, so it pulled Armchair body due to glue contact interaction. The “Sym_Axy–Armchair” contact was defined as a touching contact interaction. This interaction represents the axisymmetric boundary condition. The “Base–Armchair” pair was defined as a glue contact interaction. This interaction represented a fixed boundary condition. Finally, the “Base–tobermorite” interaction was defined as a glue contact interaction. This interaction represented a fixed boundary condition.
Furthermore, a delamination condition was defined in the interface between Armchair and tobermorite elements. This delamination condition automatically created interface cohesive elements (CZM) between both components when stresses were higher than the stablished (10−9 MPa). Very low values (numerically zero) were set to allow the software Marc always automatically split the mesh between Armchair and tobermorite elements, and introduce interface cohesive elements.

2.3.2. Physical Properties

Axisymmetric properties were defined for all elements of the model, for both Armchair and tobermorite components. Marc Element types 10 (Arbitrary Quadrilateral Axisymmetric Ring) and 2 (Axisymmetric Triangular Ring) were used [48].

2.3.3. Material Properties

Armchair and tobermorite were defined using standard isotropic elastic material model. The Young’s Modulus (E) of the Armchair is set to 98,413 MPa and the Poisson Ratio (υ) to 0.13. The Young’s Modulus (E) of the tobermorite is set to 85,000 MPa and the Poisson Ratio (υ) to 0.3.
Interface Cohesive elements (CZM) were defined using Exponential formulation with the following properties derived from the where the values obtained with MM and MD methods in order to get a coherent traction vs opening displacement curve: Cohesive Energy (Gc) of 9.64 × 10−9 mJ and Critical Opening Displacement (vc) of 7.81 × 10−11 mm.
The effective traction (t) was introduced as a function of the effective opening displacement (v) and was characterized by an initial reversible response followed by an irreversible response as soon as a critical effective opening displacement (vc) has been reached. The irreversible part was characterized by increasing damage ranging from 0 (onset of delamination) to 1 (full delamination). The exponential formulation has the following Equation (1) [48].
t = G c   v v c 2 e − v v c
in which Gc is the energy release rate (cohesive energy).

2.3.4. Analysis Methodology

An Implicit Non-Linear analysis using MSC.Marc is performed using a constant load step (1% of total load is applied in each load step), a convergence method based on residual forces criterion with automatic switching to displacement criterion was used, a relative convergence force tolerance of 10% was defined, large strain formulation was set and “Node-to-Segment” contact algorithm was selected.

3. Results and Discussion

3.1. Interfacial Energy

The variation in total potential energy with pull-out displacement for the model containing a SWCNT(6, 6) is shown in Figure 5. As the SWCNT was pulled-out from the matrix, the total potential energy increased and reached a constant value beyond a certain pull-out distance. This indicated that the SWCNT and the matrix no longer interacted. The energy of the fully embedded configuration was lower than the energy of the complete pulled-out configuration. Many other authors found similar results [23,26,27,28]. The interfacial energy is the energy difference between the fully embedded configuration of the SWCNT and the complete pulled-out configuration (Figure 5).
The interfacial energy values are presented in Figure 6 for all systems. The standard deviation of each value is also shown in this figure. In general, interfacial energy increases when SWCNT diameter and length increase. This is probably due to the greater number of atoms in longer and bigger SWCNTs, which increases the contact points between the SWCNT and the tobermorite and, hence, the interfacial energy.

3.2. Interfacial Shear Strength

For each FE model, Contact Force in X-Axis of contact body Base can be plotted with respect to position in X-Axis of contact body Tirador. The maximum value of this force can be considered as the limit value from which delamination begins to occur. This force will be used to calculate ISS dividing it by contact area between Armchair and tobermorite contact bodies (Figure 7).
ISS values are improved when SWCNT radius increases, and SWCNT length shortens. This is given by a similar extraction force (Table 1) for each diameter that leads to lower stress values when the contact area increases (length). This trend was also found by Li et al. [24] for pristine SWCNTs/polyethylene matrices. They attributed the positive correlation between SWCNT dimensions and ISS values to the increase of interfacial atoms.
It can be seen from Table 1 that, depending upon the SWCNT length and diameter, calculated ISS values ranged from 8.09 to 27.70 MPa. The effect of CNTs length on mechanical properties is controversial, as its efficiency as reinforcing agents depends on several factors such as defects on the CNT surface, concentration or dispersion in the cement matrix, although, in general, the addition of CNTs to cement leads to high improvements in mechanical properties. Some authors found that long CNTs performed better, and, to improve mechanical properties, highest amounts of short CNTs were needed [49,50].
Other authors [51] found that composites reinforced with the long CNTs showed worse properties and they attributed the efficiency of the shortest CNTs to a better filling of the nanopores within the cement matrix.
In a recent work, it was observed that the CNT length had a minimal effect on the properties of cement [52], as it is obtained in our simulations. For a given radius, the contact force can consider as constant, resulting in higher ISS for smaller areas (longer CNTs).
Thus, to validate our approach, a comparison of our results with experimental data of single SWCNT pull-out from cement matrices would be extremely useful. To the best of our knowledge, none is currently available as the experimental characterization of these interfaces still faces numerous technical problems.
It is noticeable that, for a given radius of the SWCNT, the energy rises with the length (increase of interfacial atoms) but the ISS decreases (the contact area is higher, and the ISS is obtained by dividing the force by the area).

3.3. Damage Level

If a more detailed post-processing is done in the FE models, a Damage Level can be studied. Damage value indicates the amount of irreversible cohesive energy that has been lost. A value of 1.0 indicates complete loss of cohesive energy and an area where delamination will begin to occur. Damage results can be defined as follows [48]:
D = G − G elastic G c − G elastic ; 0   ≤ D   ≤ 1 ; v e   ≥   v c
For instance, in the case of Radius 2.71 Å and Length 9.84 Å, delamination starts when cohesive elements begin to be damaged at 24% of total load applied (Displacement × Tirador: 2.64 × 10−10 mm).
In Figure 8a it can be observed that crack initiates in the interface between Armchair and tobermorite when Cohesive Interface Elements begin to have irreversible behavior (2.238 % of Damage Level). It must be noticed that deformed shape is not in real scale but scaled with a factor in order to be able to see Damage Level in the Cohesive Interface Elements.
At the end of the analysis (Figure 8b)) (Displacement × Tirador: 1.1 × 10−9 mm), it can be observed how the delamination has grown through the interface between Armchair and tobermorite and how Damage Level of Cohesive Interface Elements have also grown until 100% of Damage Level (deformed shape is also scaled with a magnification factor).
The Damage level analysis permits to confirm that the pull-out test results are affected by the cohesion energy given by molecular models. On the other hand, the results of the test are limited by the weakest element, obtaining ISS values coherent with the cement-based materials limits [53].

4. Conclusions

Pristine SWCNTs with different geometrical characteristics were pulled-out from a tobermorite matrix to study the interfacial characteristics of the resulting composites by MM and MD simulations. Only non-bonded interactions between the SWCNT and the cement matrix were considered. The interfacial energy obtained from these atomistic simulations were subsequently used as input in FEM calculations to determine ISS values. ISS values increased with larger SWCNT radius, which could be attributed to a more extensive contact surface between the SWCNT and the tobermorite. However, ISS showed an opposite trend, being larger for shorter CNTs. This is due to the increased area for longer CNTs for an almost constant force for a given radius. The weakest zone of the material seems to be the tobermorite matrix. Delamination and crack growth is mainly observed at the tobermorite side of the interface. Thus, the introduction of a SWCNT of varying geometry modifies the mechanical behavior of the matrix as confirmed by the Damage Level analysis test. The pull-out results are affected by the cohesion energy, which, in turn, depends on the CNT geometry. We think the joint use of atomistic and FEM methods seems to be a good approach to describe the interfacial properties of this kind of material. It would be very interesting to be able to compare the results of this simulation with experiments, but, to our knowledge, none is currently available.

Author Contributions

C.T., O.A.-S. and I.S.R. worked on the FEM simulations and results interpretation. I.L.-T., A.P.-P. and R.G.M.-P. worked on the Molecular models. I.L.-T. and C.T. prepared the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by Universidad Europea de Madrid and MSC. Software Simulating Reality, S.A.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated and analyzed during this study are generated by commercial software and the conditions and properties are described in the article. For that reason, they are not available as supplementary information files but are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Reales, O.A.M.; Filho, R.D.T. A review on the chemical, mechanical and microstructural characterization of carbon nanotubes-cement based composites. Constr. Build. Mater. 2017, 154, 697–710. [Google Scholar] [CrossRef]
  2. Shi, T.; Li, Z.; Guo, J.; Gong, H.; Gu, C. Research progress on CNTs/CNFs-modified cement-based composites—A review. Constr. Build. Mater. 2019, 202, 290–307. [Google Scholar] [CrossRef]
  3. Rocha, V.V.; Ludvig, P. Influence of Carbon Nanotubes on the Mechanical Behavior and Porosity of Cement Pastes Prepared by A Dispersion on Cement Particles in Isopropanol Suspension. Materials 2020, 13, 3164. [Google Scholar] [CrossRef] [PubMed]
  4. Makul, N. Advanced smart concrete—A review of current progress, benefits and challenges. J. Clean. Prod. 2020, 274, 122899. [Google Scholar] [CrossRef]
  5. Metaxa, Z.; Tolkou, A.; Efstathiou, S.; Rahdar, A.; Favvas, E.; Mitropoulos, A.; Kyzas, G. Nanomaterials in Cementitious Composites: An Update. Molecules 2021, 26, 1430. [Google Scholar] [CrossRef] [PubMed]
  6. Li, G.Y.; Wang, P.M.; Zhao, X. Mechanical behavior and microstructure of cement composites incorporating surface-treated multi-walled carbon nanotubes. Carbon N. Y. 2005, 43, 1239–1245. [Google Scholar] [CrossRef]
  7. Li, G.Y.; Wang, P.M.; Zhao, X. Pressure-sensitive properties and microstructure of carbon nanotube reinforced cement composites. Cem. Concr. Compos. 2007, 29, 377–382. [Google Scholar] [CrossRef]
  8. Gao, F.; Tian, W.; Wang, Z.; Wang, F. Effect of diameter of multi-walled carbon nanotubes on mechanical properties and microstructure of the cement-based materials. Constr. Build. Mater. 2020, 260, 120452. [Google Scholar] [CrossRef]
  9. Gao, Y.; Jing, H.; Zhou, Z.; Shi, X.; Li, L.; Fu, G. Roles of carbon nanotubes in reinforcing the interfacial transition zone and impermeability of concrete under different water-to-cement ratios. Constr. Build. Mater. 2021, 272, 121664. [Google Scholar] [CrossRef]
  10. Adhikary, S.K.; Rudžionis, Ž.; Rajapriya, R. The Effect of Carbon Nanotubes on the Flowability, Mechanical, Microstructural and Durability Properties of Cementitious Composite: An Overview. Sustainability 2020, 12, 8362. [Google Scholar] [CrossRef]
  11. Du, Y.; Gao, P.; Yang, J.; Shi, F.; Shabaz, M. Experimental Analysis of Mechanical Properties and Durability of Cement-Based Composite with Carbon Nanotube. Adv. Mater. Sci. Eng. 2021, 2021, 8777613. Available online: https://go.gale.com/ps/i.do?p=AONE&sw=w&issn=16878434&v=2.1&it=r&id=GALE%7CA696851757&sid=googleScholar&linkaccess=fulltext (accessed on 6 October 2022). [CrossRef]
  12. Yu, Z.; Lau, D. Evaluation on mechanical enhancement and fire resistance of carbon nanotube (CNT) reinforced concrete. Coupled Syst. Mech. 2017, 6, 335–349. [Google Scholar] [CrossRef]
  13. Qin, R.; Zhou, A.; Yu, Z.; Wang, Q.; Lau, D. Role of carbon nanotube in reinforcing cementitious materials: An experimental and coarse-grained molecular dynamics study. Cem. Concr. Res. 2021, 147, 106517. [Google Scholar] [CrossRef]
  14. Sindu, B.S.; Sasmal, S. Molecular dynamics simulations for evaluation of surfactant compatibility and mechanical characteristics of carbon nanotubes incorporated cementitious composite. Constr. Build. Mater. 2020, 253, 119190. [Google Scholar] [CrossRef]
  15. Yilmaz, Y.I. Analyzing Single Fiber Fragmentation Test Data by Using Stress Transfer Model. J. Compos. Mater. 2002, 36, 537–551. [Google Scholar] [CrossRef]
  16. Zhao, Y.R.; Xing, Y.M.; Lei, Z.K.; Lang, F.C. Interfacial stress transfer behavior in a specially-shaped fiber/matrix pullout test. Acta Mech. Sin. 2009, 26, 113–119. [Google Scholar] [CrossRef]
  17. Zu, M.; Li, Q.; Zhu, Y.; Dey, M.; Wang, G.; Lu, W.; Deitzel, J.M.; Gillespie, J.W.; Byun, J.-H.; Chou, T.-W. The effective interfacial shear strength of carbon nanotube fibers in an epoxy matrix characterized by a microdroplet test. Carbon N. Y. 2012, 50, 1271–1279. [Google Scholar] [CrossRef]
  18. Battisti, A.; Ojos, D.E.-D.L.; Ghisleni, R.; Brunner, A.J. Single fiber push-out characterization of interfacial properties of hierarchical CNT-carbon fiber composites prepared by electrophoretic deposition. Compos. Sci. Technol. 2014, 95, 121–127. [Google Scholar] [CrossRef]
  19. Rodríguez, M.; Molina-Aldareguía, J.; González, C.; Llorca, J. A methodology to measure the interface shear strength by means of the fiber push-in test. Compos. Sci. Technol. 2012, 72, 1924–1932. [Google Scholar] [CrossRef] [Green Version]
  20. Jean, J.-H.; Yu, C.; Park, C.-G. Bonding Characteristics of Macrosynthetic Fiber in Latex-Modified Fiber-Reinforced Cement Composites as a Function of Carbon Nanotube Content. Int. J. Polym. Sci. 2016, 2016, 7324975. [Google Scholar] [CrossRef] [Green Version]
  21. Li, Y.; Wang, Q.; Wang, S. A review on enhancement of mechanical and tribological properties of polymer composites reinforced by carbon nanotubes and graphene sheet: Molecular dynamics simulations. Compos. Part B Eng. 2019, 160, 348–361. [Google Scholar] [CrossRef]
  22. Liao, K.; Li, S. Interfacial characteristics of a carbon nanotube–polystyrene composite system. Appl. Phys. Lett. 2001, 79, 4225–4227. [Google Scholar] [CrossRef]
  23. Chowdhury, S.; Okabe, T. Computer simulation of carbon nanotube pull-out from polymer by the molecular dynamics method. Compos. Part A Appl. Sci. Manuf. 2007, 38, 747–754. [Google Scholar] [CrossRef]
  24. Li, Y.; Liu, Y.; Peng, X.; Yan, C.; Liu, S.; Hu, N. Pull-out simulations on interfacial properties of carbon nanotube-reinforced polymer nanocomposites. Comput. Mater. Sci. 2011, 50, 1854–1860. [Google Scholar] [CrossRef]
  25. Chandra, Y.; Scarpa, F.; Adhikari, S.; Zhang, J.; Flores, E.S.; Peng, H.-X. Pullout strength of graphene and carbon nanotube/epoxy composites. Compos. Part B Eng. 2016, 102, 1–8. [Google Scholar] [CrossRef] [Green Version]
  26. Chawla, R.; Sharma, S. Molecular dynamics simulation of carbon nanotube pull-out from polyethylene matrix. Compos. Sci. Technol. 2017, 144, 169–177. [Google Scholar] [CrossRef]
  27. Fan, D.; Lue, L.; Yang, S. Molecular dynamics study of interfacial stress transfer in graphene-oxide cementitious composites. Comput. Mater. Sci. 2017, 139, 56–64. [Google Scholar] [CrossRef]
  28. Alkhateb, H.; Al-Ostaz, A.; Cheng, A.H.-D.; Li, X. Materials Genome for Graphene-Cement Nanocomposites. J. Nanomech. Micromech. 2013, 3, 67–77. [Google Scholar] [CrossRef]
  29. Eftekhari, M.; Mohammadi, S. Molecular dynamics simulation of the nonlinear behavior of the CNT-reinforced calcium silicate hydrate (C–S–H) composite. Compos. Part A Appl. Sci. Manuf. 2016, 82, 78–87. [Google Scholar] [CrossRef]
  30. Kai, M.; Zhang, L.; Liew, K. Graphene and graphene oxide in calcium silicate hydrates: Chemical reactions, mechanical behavior and interfacial sliding. Carbon N. Y. 2019, 146, 181–193. [Google Scholar] [CrossRef]
  31. Joshi, U.A.; Joshi, P.; Harsha, S.P.; Sharma, S.C. Evaluation of the Mechanical Properties of Carbon Nanotube Based Composites by Finite Element Analysis. Int. J. Eng. Sci. Technol. 2010, 2, 1098–1107. [Google Scholar]
  32. Talayero, C.; Aït-Salem, O.; Gallego, P.; Páez-Pavón, A.; Merodio-Perea, R.G.; Lado-Touriño, I. Computational Prediction and Experimental Values of Mechanical Properties of Carbon Nanotube Reinforced Cement. Nanomaterials 2021, 11, 2997. [Google Scholar] [CrossRef] [PubMed]
  33. RAnsari, R.; Hemmatnezhad, M. Nonlinear finite element analysis for vibrations of double-walled carbon nanotubes. Nonlinear Dyn. 2011, 67, 373–383. [Google Scholar] [CrossRef]
  34. Alhassan, M.; Betoush, N.; Al-Huthaifi, N.; Al Dalou, A. Estimation of the fracture parameters of macro fiber-reinforced concrete based on nonlinear elastic fracture mechanics simulations. Results Eng. 2022, 15, 100539. [Google Scholar] [CrossRef]
  35. Bauchy, M.; Qomi, M.J.A.; Ulm, F.J.; Pellenq, R.J.M. Order and disorder in calcium–silicate–hydrate. J. Chem. Phys. 2014, 140, 214503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Richardson, I.G. Tobermorite/jennite- and tobermorite/calcium hydroxide-based models for the structure of C-S-H: Applicability to hardened pastes of tricalcium silicate, β-dicalcium silicate, Portland cement, and blends of Portland cement with blast-furnace slag, metakaolin, or silica fume. Cem. Concr. Res. 2004, 34, 1733–1777. [Google Scholar] [CrossRef]
  37. Hou, D.; Lu, Z.; Li, X.; Ma, H.; Li, Z. Reactive molecular dynamics and experimental study of graphene-cement composites: Structure, dynamics and reinforcement mechanisms. Carbon N. Y. 2017, 115, 188–208. [Google Scholar] [CrossRef]
  38. Allen, A.J.; Thomas, J.J.; Jennings, H.M. Composition and density of nanoscale calcium-silicate-hydrate in cement. Nat. Mater. 2007, 6, 311–316. [Google Scholar] [CrossRef] [PubMed]
  39. Hoover, W.G. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A 1985, 31, 1695. [Google Scholar] [CrossRef] [Green Version]
  40. Berendsen, H.J.C.; Postma, J.P.M.; Van Gunsteren, W.F.; Dinola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1998, 81, 3684. [Google Scholar] [CrossRef] [Green Version]
  41. BIOVIA Materials Studio—BIOVIA—Dassault Systèmes®. Available online: https://www.3ds.com/products-services/biovia/products/molecular-modeling-simulation/biovia-materials-studio/ (accessed on 6 October 2022).
  42. Sun, H. The COMPASS force field: Parameterization and validation for phosphazenes. Comput. Theor. Polym. Sci. 1998, 8, 229–246. [Google Scholar] [CrossRef]
  43. Nano Mechanical Properties on the Mineralogical Array of Calcium Silicate Hydrates and Calcium Hydroxide through Molecular Dynamics-CSIR-SERC|Request PDF. Available online: https://www.researchgate.net/publication/272702301_Nano_mechanical_properties_on_the_mineralogical_array_of_calcium_silicate_hydrates_and_calcium_hydroxide_through_molecular_dynamics-CSIR-SERC (accessed on 6 October 2022).
  44. Du, J.; Bu, Y.; Shen, Z. Interfacial properties and nanostructural characteristics of epoxy resin in cement matrix. Constr. Build. Mater. 2018, 164, 103–112. [Google Scholar] [CrossRef]
  45. Al-Ostaz, A.; Wu, W.; Cheng, D.; Song, D.A.; Associate, C.R.; Al-Ostaz, A. A molecular dynamics and microporomechanics study on the mechanical properties of major constituents of hydrated cement. Compos. Part B Eng. 2010, 41, 543–549. [Google Scholar] [CrossRef]
  46. Hajilar, S.; Shafei, B. Nano-scale investigation of elastic properties of hydrated cement paste constituents using molecular dynamics simulations. Comput. Mater. Sci. 2015, 101, 216–226. [Google Scholar] [CrossRef]
  47. Marc 2021.3—Online Help (HTML). Available online: https://help.hexagonmi.com/bundle/marc_2021.3/page/marc_main.htm (accessed on 7 October 2022).
  48. MSC Marc-Mentat—Help Reference. Available online: https://simcompanion.hexagon.com/customers/s/article/msc-marc-mentat---help-reference-doc9250 (accessed on 7 October 2022).
  49. Konsta-Gdoutos, M.S.; Metaxa, Z.S.; Shah, S.P. Multi-scale mechanical and fracture characteristics and early-age strain capacity of high performance carbon nanotube/cement nanocomposites. Cem. Concr. Compos 2010, 32, 110–115. [Google Scholar] [CrossRef]
  50. Abu Al-Rub, R.K.; Ashour, A.I.; Tyson, B.M. On the aspect ratio effect of multi-walled carbon nanotube reinforcements on the mechanical properties of cementitious nanocomposites. Constr. Build. Mater. 2012, 35, 647–655. [Google Scholar] [CrossRef]
  51. Manzur, T.; Yazdani, N.; Emon, A.B. Effect of Carbon Nanotube Size on Compressive Strengths of Nanotube Reinforced Cementitious Composites. J. Mater. 2014, 2014, 960984. [Google Scholar] [CrossRef]
  52. Mohsen, M.O.; Alansari, M.; Taha, R.; Senouci, A.; Abutaqa, A. Impact of CNTs’ treatment, length and weight fraction on ordinary concrete mechanical properties. Constr. Build. Mater. 2020, 264, 120698. [Google Scholar] [CrossRef]
  53. Centro de Publicaciones; Ministerio de Fomento; Secretaría General Técnica. Instrucción de Hormigon Estructural. EHE-08, 5th ed.; Ministerio de Transportes, Movilidad y Agenda urbana: Madrid, Spain, 2011; p. 702. [Google Scholar]
Materials 16 01992 g001 550
Figure 1. Cell used for the atomistic pull-out simulations.
Figure 1. Cell used for the atomistic pull-out simulations.
Materials 16 01992 g001
Materials 16 01992 g002 550
Figure 2. Models used in the simulation; cement: tobermorite crystal (left) and composite: tobermorite/SWCNT.
Figure 2. Models used in the simulation; cement: tobermorite crystal (left) and composite: tobermorite/SWCNT.
Materials 16 01992 g002
Materials 16 01992 g003 550
Figure 3. Initial and final structures before (top) and after (bottom) the SWCNT pull-out process.
Figure 3. Initial and final structures before (top) and after (bottom) the SWCNT pull-out process.
Materials 16 01992 g003
Materials 16 01992 g004 550
Figure 4. Geometry and boundary conditions for the FEM model.
Figure 4. Geometry and boundary conditions for the FEM model.
Materials 16 01992 g004
Materials 16 01992 g005 550
Figure 5. Total potential energy versus pull−out displacement.
Figure 5. Total potential energy versus pull−out displacement.
Materials 16 01992 g005
Materials 16 01992 g006 550
Figure 6. Interfacial energy as a function of SWCNT diameter and length.
Figure 6. Interfacial energy as a function of SWCNT diameter and length.
Materials 16 01992 g006
Materials 16 01992 g007 550
Figure 7. ISS as a function of SWCNT diameter and length.
Figure 7. ISS as a function of SWCNT diameter and length.
Materials 16 01992 g007
Materials 16 01992 g008 550
Figure 8. Damage level for for Radius 2.71 Å and Length 9.84 Å (a) Damage level when delamination starts. (b) Damage level when damage reaches 100%.
Figure 8. Damage level for for Radius 2.71 Å and Length 9.84 Å (a) Damage level when delamination starts. (b) Damage level when damage reaches 100%.
Materials 16 01992 g008
Table
Table 1. Contact force and ISS of the different systems.
Table 1. Contact force and ISS of the different systems.
Radius (Å)Length (Å)Maximum Contact Force (N × 10−11)Contact Area (mm2 × 10−12)ISS (MPa)
2.719.843.481.6820.740
19.683.403.3510.158
24.63.394.198.091
3.399.845.022.1023.932
19.684.874.1911.620
24.64.835.249.210
4.079.846.512.5225.859
19.686.495.0312.896
24.66.446.2910.237
4.749.848.092.9327.602
19.688.275.8614.115
24.68.207.3311.194
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Talayero, C.; Lado-Touriño, I.; Aït-Salem, O.; Ramos, I.S.; Páez-Pavón, A.; G. Merodio-Perea, R. Interfacial Shear Strength of Single-Walled Carbon Nanotubes-Cement Composites from Molecular Dynamics and Finite Element Studies. Materials 2023, 16, 1992. https://doi.org/10.3390/ma16051992

AMA Style

Talayero C, Lado-Touriño I, Aït-Salem O, Ramos IS, Páez-Pavón A, G. Merodio-Perea R. Interfacial Shear Strength of Single-Walled Carbon Nanotubes-Cement Composites from Molecular Dynamics and Finite Element Studies. Materials. 2023; 16(5):1992. https://doi.org/10.3390/ma16051992

Chicago/Turabian Style

Talayero, Carlos, Isabel Lado-Touriño, Omar Aït-Salem, Ismael Sánchez Ramos, Alicia Páez-Pavón, and Rosario G. Merodio-Perea. 2023. "Interfacial Shear Strength of Single-Walled Carbon Nanotubes-Cement Composites from Molecular Dynamics and Finite Element Studies" Materials 16, no. 5: 1992. https://doi.org/10.3390/ma16051992

APA Style

Talayero, C., Lado-Touriño, I., Aït-Salem, O., Ramos, I. S., Páez-Pavón, A., & G. Merodio-Perea, R. (2023). Interfacial Shear Strength of Single-Walled Carbon Nanotubes-Cement Composites from Molecular Dynamics and Finite Element Studies. Materials, 16(5), 1992. https://doi.org/10.3390/ma16051992

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop