Crystalline hydrates of swelling clay minerals (smectites) exhibit a strong coupling between their ion exchange and hydration/dehydration reactions. The uptake or removal of water from smectite interlayers as a result of a change in the environmental conditions also leads to the partitioning of cations. Three factors, the solid ion composition, the solid basal spacing/water content, and the aqueous solution composition, are all implicated in controlling the thermodynamics of ion exchange. However, conventional approaches to measuring the exchange free energy cannot separate the influence of each of these individual factors. Here, we explore the energetics of the swelling and ion exchange reactions in montmorillonite using a potential of mean force approach and the thermodynamic integration method within molecular simulations. We investigate the influence of solution and clay composition on the spontaneity of the reactions, focusing on the 2 water-layer hydration state. The swelling simulations provide the equilibrium water content, interlayer water structure, and basal spacings, while thermodynamic integration of sodium–potassium exchange in the aqueous solution and solid phase are combined to calculate ion exchange free energies as a function of solution composition. Results confirm the tendency of the clay to collapse to lower hydration states as the concentration of the solution increases. Changes to the equilibrium water content, even at fixed hydration states, and the composition of the mixed electrolyte solution play a critical role in driving ion exchange and the selectivities of the clay to the exchanged cation, while the composition of the solid phase is shown to be insignificant. These findings underscore the extreme sensitivity of clay swelling and ion exchange thermodynamics to small (tenths of an Angstrom) deviations in layer spacing.

@en

The effectiveness of carbonaceous nanoparticles in arresting and delaying damage in nanocomposites has been attributed to multiscale toughening mechanisms. To explore their application in joined interfaces of composites, this study investigates the use of carbon nanotube (CNT) interleaved films for co-cured joining of composite parts and their consequent effects on the interfacial fracture toughness. Carbon nanotubes dispersed in a thermoset resin into thin films of two discrete thicknesses (200 μ and 500 μ) and three concentrations of CNT dispersion were chosen for this study (0.5% wt., 1% wt., and 2% wt.). The films were semi-cured in the oven before being incorporated as interleaves in the composite laminate interface. Fracture toughness of the interface in mode I loading conditions was determined through double cantilever beam (DCB). Micrographs of the fracture surfaces reveal a slip-and-stick based crack jump and arrest phenomena in mode I when nanoparticles are added to the interleaved interface. The thickness of the interleaves has a more significant effect on mode I toughening mechanisms than the concentration of the nanoparticles. @en

Bonded and co-cured composites are popular alternatives to structures joined with mechanical fasteners in aircraft but the complex and coupled damage mechanisms in the co-cured/bonded region are poorly understood, thus making the evaluation of their strength and durability difficult with current modelling strategies. This study explores the potential of interleaf inclusion in failure-prone, critical regions of co-cured composite specimens in improving the joint strength and interface fracture toughness and strives to advance the understanding of damage initiation in the co-cured region using an atomistic model. A two-pronged approach is pursued here with bench-scale experimental testing and molecular modelling in this study. Experiments are performed for mode I fracture toughness with double cantilever beam (DCB) on composite laminates with an epoxy interleaf layer. Two epoxy resins and three methods for interleaf inclusion are explored in this study; we supplement the results from DCB testing with insights from confocal microscopy on the crack tip and the interleaf layer pre- and post-testing. Molecular dynamic (MD) simulations capture the cohesive interactions at the threephase interface containing the carbon fiber, the prepreg epoxy, and the interleaf epoxy. Results highlight that an interleaf layer made from partially-cured and filmed epoxy, further consolidated in the composite lay-up is the most effective way to suppress void formation, improve dispersion, and maximize cohesive interactions at the interface of co-cured composites. @en

Montmorillonite (MMT) clay - a layered porous nanomaterial used as seals in engineered waste containment barriers for spent nuclear fuel - adopts discrete hydration/swelling states depending upon surrounding water and ion activities and confining pressure. The structure of nanoconfined water and charge-balancing counterions in the clay mineral interlayers dictate the swelling and mechanical behavior of MMT, so a molecular model for this clay with high structural fidelity is required to accurately predict the reliability of long-term nuclear waste storage. Here, we present a molecular model for MMT that is based on high resolution transmission electron microscopy of Wyoming-MMT single crystals. Imaging data unambiguously show a cis-vacant arrangement of structural hydroxyl groups in the octahedral sheet, whereas existing molecular models assume a centrosymmetric trans-vacant configuration for MMT. Using atomistic simulations, we find that the cis-vacant arrangement of structural hydroxyl groups significantly affects the structure of adsorbed water yielding a larger population of hydrogen bonds with bridging oxygens on the tetrahedral sheet and weak hydrogen bonding between the hydroxyl groups in the octahedral sheet and water in the clay mineral interlayers. As a result, water adsorbed in the interlayer is more "ice-like", with stronger ordering and lower density, although the diffusivity of the interlayer species is not significantly diminished. Our improved structural model for MMT provides insight into the energetics of water adsorption, which ultimately dictates its pore- to macro-scale swelling, transport, and fracture properties.

@en

An atomistically-informed multiscale modeling framework to investigate damage evolution and failure in radially-grown carbon nanotube (CNT) architecture is detailed in this paper. Molecular dynamics (MD) simulations are performed to investigate the effects of nano-reinforcements on the elastic-plastic characteristics of the constituent interphase in the radially-grown nanocomposite. An interphase damage model is developed using the continuum damage mechanics approach with damage evolution equations derived using atomistic simulations. The developed damage model is integrated with a high-fidelity micromechanical analysis and captures the underlying physical behavior that could be attributed to the enhancement of the out-of-plane properties at higher length scales. The mechanical properties of the nanocomposite obtained from micromechanical simulations are compared to experimental values reported in the literature, to validate the developed modeling framework. Conclusions are presented by comparing the material response of radially-grown CNT architectures with the traditional fiber reinforced polymer (FRP) with dispersed CNT architecture.

@en

This paper presents the implementation of an atomistic computational framework to investigate a fuzzy fiber nanocomposite architecture. A polymeric functional coating for the carbon fiber surface, which also serves as a substrate for the CNT growth, is explicitly modeled. Additionally, the carbon fiber surface is modeled through irregularly stacked graphene layers with voids. Radially grown CNTs and epoxy resin and hardener molecules are also included at the nanocomposite interface. The epoxy curing simulation is carried out first, followed by virtual deformation of the simulation volume. The tensile and transverse moduli of the fuzzy fiber nanocomposite interface is computed from the virtual deformation simulations. Results indicate that the out-of-plane interface modulus is significantly improved, compared to the traditional fiber/matrix interface by the addition of the polymer coating and radially-grown CNTs. The in-plane transvers modulus also shows improvement and the atomistic modeling framework is able to capture the physical mechanisms that lead to improvement/degradation of material properties.@en

An atomistic modeling framework to investigate the interface/interphase of composite architecture with carbon fibers containing radially-grown carbon nanotubes (often called fuzzy fibers) is detailed in this paper. A polymeric functional coating for the carbon fiber surface, which also serves as a substrate for the CNT growth, is explicitly modeled. The tensile and transverse moduli of the fuzzy fiber/epoxy interphase is computed from virtual deformation simulations and compared to experimental values reported in literature, in order to validate the nanoscale model. Furthermore, the effect of the polymer substrate is studied by modeling the local interphase mechanics. Various modes of virtual loading provide the cohesive behavior of the local substrate/epoxy interphase. Conclusions are presented by comparing the material response of the interphase with and without the polymeric substrate. The integration of results from the nanoscale to an atomistically-informed subcell-based continuum level model is also demonstrated in the paper.

@en

A methodology that accurately simulates the brittle behavior of epoxy polymers initiating at the molecular level due to bond elongation and subsequent bond dissociation is presented in this paper. The system investigated in this study comprises a combination of crystalline carbon nanotubes (CNTs) dispersed in epoxy polymer molecules. Molecular dynamics (MD) simulations are performed with an appropriate bond order-based force field to capture deformation-induced bond dissociation between atoms within the simulation volume. During deformation, the thermal vibration of molecules causes the elongated bonds to re-equilibrate; thus, the effect of mechanical deformation on bond elongation and scission cannot be captured effectively. This issue is overcome by deforming the simulation volume at zero temperature—a technique adopted from the concept of quasi-continuum and demonstrated successfully in the authors’ previous work. Results showed that a combination of MD deformation tests with ultra-high strain rates at near-zero temperatures provides a computationally efficient alternative for the study of bond dissociation phenomenon in amorphous epoxy polymer. In this paper, the ultra-high strain rate deformation approach is extended to the CNT-epoxy system at various CNT weight fractions and the corresponding bond disassociation energy extracted from the simulation volume is used as input to a low-fidelity continuum damage mechanics (CDM) model to demonstrate the bridging of length scales and to study matrix failure at the microscale. The material parameters for the classical CDM model are directly obtained from physics-based atomistic simulations, thus improving the accuracy of the multiscale approach.

@en

A micromechanical modeling framework that integrates an interphase damage model is used to study damage evolution, and failure in radially-grown carbon nanotube (CNT) reinforced nanocomposites. The atomistically-informed damage model is developed using the continuum damage mechanics approach with damage evolution equations derived using molecular dynamics simulations. The developed damage model is integrated with a high-fidelity micromechanical analysis and captures the underlying physical behavior that could be attributed to the enhancement of the out-of-plane properties at the higher length scales. To illustrate the effects of nanoarchitecture, the response of the radially-grown CNT nanocomposite is compared with that of the traditional carbon fiber reinforced polymer (CFRP) composite.@en

A multiscale-modeling framework is presented to understand damage and failure response in carbon nanotube reinforced nanocomposites. A damage model is developed using the framework of continuum damage mechanics with a physical damage evolution equation inspired by molecular dynamics simulations. This damage formulation is applied to randomly dispersed carbon nanotube reinforced nanocomposite unit cells with periodic boundary conditions to investigate preferred sites and the tendency towards damage. The continuum model is seen as successfully capturing much of the unique nonlinear trends observed in the molecular dynamics simulations in a volume 1000 times greater than the molecular dynamics unit cell. Additionally, application of the damage model to the continuum unit cell revealed insights into the failure of carbon nanotube reinforced nanocomposites at the sub-microscale.@en

The carbon fiber/polymer matrix interphase region plays an important role in failure initiation and accurate modeling techniques are integral to study the effects of this complex region on the composite response. In composites infused with nanoparticles such as carbon nanotubes (CNT), the interphase region is more complex due to the presence of multiple constituents and their interactions with each other. An atomistic methodology to simulate the constituent interphases in carbon fiber reinforced CNT/epoxy nanocomposites is presented in this paper. The interphase model consisting of voids in multiple graphene layers enable the simulation of physical entanglement between the polymer matrix and the irregular carbon fiber surface. The voids in graphene layers are generated by removing carbon atoms and hydrogenating the end carbon atoms, which better represent the roughness of the carbon fiber surface. The epoxy curing studies and the response of fiber/matrix interphase under mechanical loading are investigated through molecular dynamic (MD) simulations with appropriate classical/harmonic and bond order-based force fields. Furthermore, the atomistic force-displacement behavior is also extracted to formulate a traction-separation law for interface cohesive zone models. The cohesive behavior determined from molecular models is parameterized in equations that can be integrated with an atomistically informed multiscale modeling framework.

@en

This paper examines microscale and sub-microscale damage mechanisms in carbon nanotube (CNT) reinforced nanocomposites. A multiscale modeling framework with a damage model developed from molecular dynamics simulation, is employed to study the physical mechanisms of damage initiation and propagation in CNT nanocomposites at the sub-microscale. Two CNT arrangements, randomly dispersed and entangled agglomerates, are examined. This investigation offers insights into damage properties of particular configurations of CNTs in a polymer matrix, in addition to specific understanding related to damage concentration effects around the filler material at the sub-microscale. High spatial CNT concentration differential is observed to affect damage initiation and rate of damage. It is further shown to result in sub-microscale crack initiation at low global strains, a phenomenon that is also observed at agglomeration boundaries, which results in CNT agglomerations to behave as crack initiation sites.

@en

The carbon fiber/polymer matrix interphase region plays an important role in the behavior and failure initiation of polymer matrix composites and accurate modeling techniques are needed to study the effects of this complex region on the composite response. This paper presents a high fidelity multiscale modeling framework integrating a novel molecular interphase model for the analysis of polymer matrix composites. The interphase model, consisting of voids in multiple graphene layers, enables the physical entanglement between the polymer matrix and the carbon fiber surface. The voids in the graphene layers are generated by intentionally removing carbon atoms, which better represents the irregularity of the carbon fiber surface. The molecular dynamics method calculates the interphase mechanical properties at the nanoscale, which are integrated within a high fidelity micromechanics theory. Additionally, progressive damage and failure theories are used at different scales in the modeling framework to capture scale-dependent failure of the composite. Comparisons between the current molecular interphase model and existing interphase models and experiments demonstrate that the current model captures larger stress gradients across the material interphase. These large stress gradients increase the viscoplasticity and damage effects at the interphase which are necessary for improved prediction of the nonlinear response and multiscale damage in composite materials.

@en

A novel molecular dynamics (MD) simulation methodology to capture brittle fracture in epoxy-based thermoset polymer under mechanical loading is presented. The ductile behavior of amorphous polymers has been captured through traditional MD simulation methods by estimating the stress-strain response beyond the yield point; however, brittle fracture in highly crosslinked polymer materials such as epoxy thermoset has not been addressed appropriately and is the primary objective of this work. In this study, a numerically cured epoxy system comprising molecules of epoxy resin and hardener is generated. During the virtual deformation test, it is observed that the inherent molecular vibration due to temperature re-equilibrates the elongated covalent bonds and this molecular vibration impedes further stretching the covalent bond leading to scission. In order to overcome the influence of thermal vibration, an approach that employs deformation tests at absolute zero temperature condition - a concept borrowed from the quasi-continuum method, is developed. Bond dissociation energy is measured to quantify the extent of failure in the system by calculating the bond potentials during the deformation tests. Applying zero temperature condition to the deformation test, however, requires a large amount of computational time due to intermediate energy minimization processes. To improve the computational efficiency, an ultra-high strain rate (UHSR ≈ 10¹³ s^-1) approach is developed by which the thermal vibration is decoupled from the deformation test using a strain rate higher than molecular vibration frequency. Note that the deformation test performed by traditional MD simulation methods used strain rates ranging 10⁸-10¹⁰ s^-1. Simulation results show that the UHSR approach successfully captures brittle fracture in epoxy polymer due to covalent bond dissociation with high computational efficiency.

@en

The piezoresistivity of carbon nanotube (CNT) reinforced nanocomposites is modeled using a multiscale damage modeling technique. Two phenomena of piezoresistivity are studied, the inherent piezoresistivity of the CNTs and the electrical tunneling effect. The damage model is developed under the framework of continuum damage mechanics (CDM) with a physical damage evolution equation inspired by Molecular Dynamics (MD) simulations. This damage model is applied to a nanocomposite unit cells with randomly dispersed CNTs. Orders of magnitude change in piezoresistivity is observed as the nanocomposite changes from non-damaged state to damaged state. This study provides insights into the prevailing mechanisms associated with piezoresistivity in the damaged and undamaged state of the CNT reinforced nanocomposites at the sub micro scale.

@en

This article presents a novel approach to model the mechanical response of smart polymeric materials. A cyclobutane-based mechanophore, named "smart particle" in this article, is embedded in an epoxy polymer matrix to form the self-sensing smart material. A spring-bead model is developed based on the results from molecular dynamics simulation at the nanoscale to represent bond clusters of a smart polymer. The spring-bead network model is developed through parametric studies and mechanical equivalence optimization to represent the microstructure of the material. A statistical network model is introduced, which is capable of bridging the high-accuracy molecular dynamics model at the nanoscale and the computationally efficient finite element model at the macroscale. A comparison between experimental and simulation results shows that the multiscale model can capture global mechanical response and local material properties.

@en

An atomistic methodology to simulate the constituent interphases in carbon fiber reinforced CNT/epoxy nanocomposites is presented in this paper. Two critical interphase regions are considered in the study: CNT/polymer interphase and the fiber/matrix interphase. The elastic and inelastic responses of the interphases are investigated through molecular dynamic (MD) simulations with appropriate force fields, and integrated with an atomistically informed multiscale modeling framework. The elastic behavior is studied at the molecular level using harmonic force fields whereas bond elongation and subsequent bond dissociation in epoxy polymer chains, and the fiber/matrix interphase are investigated using a bond order based force field. An MD simulation approach adopted from the concept of quasi-continuum (QC) is employed to calculate the bond dissociation energy during virtual deformation tests. The variation of bond dissociation energy density during the deformation tests is integrated into a continuum damage mechanics model to characterize microscale damage in the epoxy matrix. Furthermore, the atomistic forcedisplacement behavior is also extracted to formulate a tractionseparation law for the microscale cohesive zone models for the fiber/matrix interface.@en

In this paper a multiscale-modeling framework is presented wherein fundamental damage information at the atomic level is coupled with a sectional micromechanics model for the nonlinear and damage analysis of carbon fiber reinforced polymer (CFRP) composites. Damage information in the polymer matrix originating from the atomic scale, as investigated using molecular dynamics (MD) simulations, is transferred to the continuum length scale using a continuum damage mechanics approach with a physical damage evolution equation. Such a framework is shown to be computationally efficient for the linear and damage analysis of CFRP composites and reasonably accurate when compared to experimental data and verified models in literature. Furthermore, material uncertainty, such as curing degree variation in polymers, can be computationally calculated leading to a computational stochastic analysis of the CFRP composites. Thus, such a framework can be used to investigate the damage mechanics of ply level CFRP components at the nano, micro and macro length scales.@en

This paper presents the development of a novel methodology for modeling the interphase between carbon fiber and polymer matrix using atomistic scale simulations. The model is integrated within a multiscale framework for the analysis of polymer matrix composites. The interphase region that exists between carbon fiber and polymer matrix plays an important role in damage initiation and it is critical for modeling damage propagation across length scales. The developed interphase model in this work consists of multiple pseudo-damaged graphene layers, and is capable of capturing the physical entanglement between the polymer matrix and the carbon fiber through the use of molecular dynamics simulations. The pseudo-damage model simulates voids in the graphene layer by removing carbon atoms, while multiple pseudo-damaged graphene layers represent the amorphous structure of the carbon fiber surface. In order to accurately model the chemical structure of the polymer, a numerical curing process is employed that captures the stochastic interactions between resin and hardener molecules for crosslinking (carbon-nitrogen covalent bond generation). The molecular dynamics model also enables extraction of the material properties of the interphase at the nanoscale which is then integrated within a high fidelity generalized method of cells micromechanics theory. A modified Bodner-Partom viscoplastic theory is applied to the polymer matrix subcells. Additionally, progressive damage and failure theories are used at different scales in the modeling framework to accurately capture scale-dependent failure of the composite material. Comparative studies are performed by applying different interphase properties that are obtained from the literature to the interphase subcells. The results indicate that the responses for each type of interphase model are similar, but the transverse tensile strength, obtained using the MD simulated interphase properties, is approximately 7% lower than the other interphase types and also has a maximum difference of 10% in failure strain.@en

A multiscale methodology that accurately simulates the inelastic behavior of epoxy polymers initiating at the molecular level due to bond elongation and subsequent bond dissociation is presented in this paper. The system investigated in this study comprises a combination of crystalline carbon nanotubes (CNTs) dispersed in amorphous epoxy polymer molecules. Molecular dynamics (MD) simulations are performed with an appropriate bond order based force field to capture deformation-induced bond dissociation between atoms within the simulation volume. In order to overcome the influence of thermal vibrations of bonds on bond dissociation energy (BDE), a quasi-continuum (QC) approach, excluding the effects of temperature, is explored. Results indicate that a QC approach can simulate bond breakage leading to brittle behavior in amorphous epoxy polymer. The novel combination of MD deformation tests with high strain rates at near-zero temperatures, however, is seen to provide a more computationally efficient alternative for the study of bond dissociation phenomenon in amorphous epoxy polymer. The corresponding BDE extracted from the simulation volume is used as input to the continuum damage mechanics (CDM) model to study matrix failure at the microscale. The material parameters for the CDM model are directly obtained from physics-based atomistic simulations, thus, significantly reducing the propagation of errors in the multiscale framework.@en