LUNAR: Automated Input Generation and Analysis for Reactive LAMMPS Simulations.

Generating simulation-ready molecular models for the LAMMPS molecular dynamics (MD) simulation software package is a difficult task and impedes the more widespread and efficient use of MD in materials design and development. Fixed-bond force fields generally require manual assignment of atom types, bonded interactions, charges, and simulation domain sizes. A new LAMMPS pre- and postprocessing toolkit (LUNAR) is presented that efficiently builds molecular systems for LAMMPS. LUNAR automatically assigns atom types, generates bonded interactions, assigns charges, and provides initial configuration methods to generate large molecular systems. LUNAR can also incorporate chemical reactivity into simulations by facilitating the use of the REACTER protocol. Additionally, LUNAR provides postprocessing for free volume calculations, cure characterization calculations, and property predictions from LAMMPS thermodynamic outputs. LUNAR has been validated via building and simulation of pure epoxy and cyanate ester polymer systems with a comparison of the corresponding predicted structures and properties to benchmark values, including experimental results from the literature. LUNAR provides the tools for the computationally driven development of next-generation composite materials in the Integrated Computational Materials Engineering (ICME) and Materials Genome Initiative (MGI) frameworks. LUNAR is written in Python with the usage of NumPy and can be used via a graphical user interface, a command line interface, or an integrated design environment. LUNAR is freely available via GitHub.

Optimal Molecular Dynamics System Size for Increased Precision and Efficiency for Epoxy Materials.

Molecular dynamics (MD) simulation is an important tool for predicting thermo-mechanical properties of polymer resins at the nanometer length scale, which is particularly important for efficient computationally driven design of advanced composite materials and structures. Because of the statistical nature of modeling amorphous materials on the nanometer length scale, multiple MD models (replicates) are typically built and simulated for statistical sampling of predicted properties. Larger replicates generally provide higher precision in the predictions but result in higher simulation times. Unfortunately, there is insufficient information in the literature to establish guidelines between MD model size and the resulting precision in predicted thermo-mechanical properties. The objective of this study was to determine the optimal MD model size of epoxy resin to balance efficiency and precision. The results show that an MD model size of 15,000 atoms provides for the fastest simulations without sacrificing precision in the prediction of mass density, elastic properties, strength, and thermal properties of epoxy. The results of this study are important for efficient computational process modeling and integrated computational materials engineering (ICME) for the design of next-generation composite materials for demanding applications.

Multiscale Process Modeling of Semicrystalline PEEK for Tailored Thermomechanical Properties.

Polyether ether ketone (PEEK) is a semicrystalline thermoplastic that is used in high-performance composites for a wide range of applications. Because the crystalline phase has a higher mass density than that of the amorphous phase, the evolution of the crystalline phase during high-temperature annealing processing steps results in the formation of residual stresses and laminate deformations, which can adversely affect the composite laminate performance. Multiscale process modeling, utilizing molecular dynamics, micromechanics, and phenomenological PEEK crystal kinetic laws, is used to predict the evolution of volumetric shrinkage, elastic properties, and thermal properties, as a function of crystalline phase evolution, and thus annealing time, in the 306-328 °C temperature range. The results indicate that lower annealing temperatures in this range result in a faster evolution of thermomechanical properties and shrinkage toward the pure crystalline values. Therefore, from the perspective of composite processing, it may be more advantageous to choose the higher annealing rates in this range to slow the volumetric shrinkage and allow PEEK stress relaxation mechanisms more time to relax internal residual stresses in PEEK composite laminates and structures.

Evolution of Glassy Carbon Derived from Pyrolysis of Furan Resin.

Glassy carbon (GC) material derived from pyrolyzed furan resin was modeled by using reactive molecular dynamics (MD) simulations. The MD polymerization simulation protocols to cure the furan resin precursor material are validated via comparison of the predicted density and Young's modulus with experimental values. The MD pyrolysis simulations protocols to pyrolyze the furan resin precursor is validated by comparison of calculated density, Young's modulus, carbon content, sp2 carbon content, the in-plane crystallite size, out-of-plane crystallite stacking height, and interplanar crystallite spacing with experimental results from the literature for furan resin derived GC. The modeling methodology established in this work can provide a powerful tool for the modeling-driven design of next-generation carbon-carbon composite precursor chemistries for thermal protection systems and other high-temperature applications.

Simple and convenient mapping of molecular dynamics mechanical property predictions of bisphenol-F epoxy for strain rate, temperature, and degree of cure.

It is well-known that all-atom molecular dynamics (MD) predictions of mechanical properties of thermoset resins suffer from multiple accuracy issues associated with their viscoelastic nature. The nanosecond simulation times of MD simulations do not allow for the direct simulation of the molecular conformational relaxations that occur under laboratory time scales. This adversely affects the prediction of mechanical properties at realistic strain rates, intermediate degrees of cure, and elevated temperatures. While some recent studies have utilized a time-temperature superposition approach to relate MD predictions to expected laboratory observations, such an approach becomes prohibitively difficult when simulating thermosets with a combination of strain rates, intermediate degrees of cure, and temperatures. In this study, a phenomenological approach is developed to map the predictions of Young's modulus and Poisson's ratio for a DGEBF/DETDA epoxy system to the corresponding laboratory-based properties for intermediate degrees of cure and temperatures above and below the glass transition temperature. The approach uses characterization data from dynamical mechanical analysis temperature sweep experiments. The mathematical formulation and experimental characterization of the mapping is described, and the resulting mapping of computationally-predicted to experimentally-observed elastic properties for various degrees of cure and temperatures are demonstrated and validated. This mapping is particularly important to mitigate the strain-rate effect associated with MD predictions, as well as to accurately predict mechanical properties at elevated temperatures and intermediate degrees of cure to facilitate accurate and efficient composite material process modeling.

Scalable High Tensile Modulus Composite Laminates Using Continuous Carbon Nanotube Yarns for Aerospace Applications.

An approach is established for fabricating high-strength and high-stiffness composite laminates with continuous carbon nanotube (CNT) yarns for scaled-up mechanical tests and potential aerospace structure applications. Continuous CNT yarns with up to 80% degree of nanotube alignment and a unique self-assembled graphitic CNT packing result in their specific tensile strengths of 1.77 ± 0.07 N/tex and an apparent specific modulus of 92.6 ± 3.2 N/tex. Unidirectional CNT yarn reinforced composite laminates with a CNT concentration of greater than 80 wt % and minimal microscale voids are fabricated using filament winding and aerospace-grade resin matrices. A specific tensile strength of up to 1.71 GPa/(g cm-3) and specific modulus of 256 GPa/(g cm-3) are realized; the specific modulus exceeds current state-of-the-art unidirectional carbon fiber composite laminates. The specific modulus of the laminates is 2.76 times greater than the specific modulus of the constituent CNT yarns, a phenomenon not observed in carbon fiber reinforced composites. The results demonstrate an effective approach for fabricating high-strength CNT yarns into composites for applications that require specific tensile modulus properties that are significantly beyond state-of-the-art carbon fiber composites and potentially open an unexplored performance region in the Ashby chart for composite material applications.

Probing the Influence of Surface Chemical Functionalization on Graphene Nanoplatelets-Epoxy Interfacial Shear Strength Using Molecular Dynamics.

In this work, a characterization study of the interfacial interaction between different types of graphene nanoplatelets and an epoxy matrix is computationally performed. To quantify the discrete mutual graphene-epoxy "interfacial interaction energy" (IIE) within the nanocomposite, molecular dynamics simulations with a reactive force field are performed on a localized model of the suggested nanocomposite. Pull-out molecular dynamics simulations are also performed to predict the interfacial shear strength between the two constituents. The results indicate a significant increase in interfacial adhesion of functionalized nanoplatelets with the hosting epoxy matrix relative to virgin graphene nanoplatelets. The obtained results also demonstrate a dramatic increase in the interfacial interaction energy (IIE) (up to 570.0%) of the functionalized graphene/epoxy nanocomposites relative to the unmodified graphene/epoxy nanocomposites. In the same context, the surface functionalization of graphene nanoplatelets with the polymer matrix leads to a significant increase in the interfacial shear strength (ISS) (up to 750 times). The reported findings in this paper are essential and critical to producing the next generation of lightweight and ultra-strong polymer-based nanocomposite structural materials.

A Micromechanical Study of Interactions of Cyanate Ester Monomer with Graphene or Boron Nitride Monolayer.

Polymer composites, hailed for their ultra-strength and lightweight attributes, stand out as promising materials for the upcoming era of space vehicles. The selection of the polymer matrix plays a pivotal role in material design, given its significant impact on bulk-level properties through the reinforcement/polymer interface. To aid in the systematic design of such composite systems, molecular-level calculations are employed to establish the relationship between interfacial characteristics and mechanical response, specifically stiffness. This study focuses on the interaction of fluorinated and non-fluorinated cyanate ester monomers with graphene or a BN monolayer, representing non-polymerized ester composites. Utilizing micromechanics and the density functional theory method to analyze interaction energy, charge density, and stiffness, our findings reveal that the fluorinated cyanate-ester monomer demonstrates lower interaction energy, reduced pull-apart force, and a higher separation point compared to the non-fluorinated counterpart. This behavior is attributed to the steric hindrance caused by fluorine atoms. Furthermore, the BN monolayer exhibits enhanced transverse stiffness due to increased interfacial strength, stemming from the polar nature of B-N bonds on the surface, as opposed to the C-C bonds of graphene. These molecular-level results are intended to inform the design of next-generation composites incorporating cyanate esters, specifically for structural applications.

Computational study of effect of radiation induced crosslinking on the properties of flattened carbon nanotubes.

Flattened carbon nanotubes (flCNTs) are a primary component of many carbon nanotube (CNT) yarn and sheet materials, which are promising reinforcements for the next generation of ultra-strong composites for aerospace applications. Significant improvements in the performance of CNT materials can be realized with improvements in the load transfer between flCNTs, which are generally oriented at different angles with respect to each other. An intriguing approach to improving the load transfer is via irradiation-induced chemical crosslinking between adjacent flCNTs. The objective of this research is to use molecular dynamics (MD) simulations to predict the behavior of flCNT junctions with 0- and 90-degree orientations and varying levels of crosslinking. The results indicate that crosslinking improves the flCNT interfacial load transfer for both orientations, but degrades the flCNT tensile response. The primary toughening mechanism at the flCNT/flCNT interface is the formation of carbon chains that provide load transfer up to the point of total rupture. Based on these results, it is clear that irradiation-induced crosslinking is beneficial in CNT-based composite systems in which interfacial load transfer between flCNTs is of primary importance, even though individual flCNTs may lose some mechanical integrity with crosslinking.

Accurate predictions of thermoset resin glass transition temperatures from all-atom molecular dynamics simulation.

To enable the design and development of the next generation of high-performance composite materials, there is a need to establish improved computational simulation protocols for accurate and efficient prediction of physical, mechanical, and thermal properties of thermoset resins. This is especially true for the prediction of glass transition temperature (Tg), as there are many discrepancies in the literature regarding simulation protocols and the use of cooling rate correction factors for predicting values using molecular dynamics (MD) simulation. The objectives of this study are to demonstrate accurate prediction the Tg with MD without the use of cooling rate correction factors and to establish the influence of simulated conformational state and heating/cooling cycles on physical, mechanical, and thermal properties predicted with MD. The experimentally-validated MD results indicate that accurate predictions of Tg, elastic modulus, strength, and coefficient of thermal expansion are highly reliant upon establishing MD models with mass densities that match experiment within 2%. The results also indicate the cooling rate correction factors, model building within different conformational states, and the choice of heating/cooling simulation runs do not provide statistically significant differences in the accurate prediction of Tg values, given the typical scatter observed in MD predictions of amorphous polymer properties.

Molecular Dynamics Modeling of Interfacial Interactions between Flattened Carbon Nanotubes and Amorphous Carbon: Implications for Ultra-Lightweight Composites.

Flattened carbon nanotubes (flCNTs) naturally form in many carbon nanotube-based materials and can exhibit mechanical properties similar to round carbon nanotubes but with tighter packing and alignment. To facilitate the design, fabrication, and testing of flCNT-based composites for aerospace structures, computational modeling can be used to efficiently and accurately predict their performance as a function of processing parameters, such as reinforcement/matrix cross-linking. In this study, molecular dynamics modeling is used to predict the load transfer characteristics of the interface region between the flat region of flCNTs (i.e., bi-layer graphene) and amorphous carbon (AC) with various levels and locations of covalent bond cross-linking and AC mass density. The results of this study show that increasing the mass density of AC at the interface improves the load transfer capability of the interface. However, a much larger improvement is observed when cross-linking is added both to the flCNT-AC interface and between the flCNT sheets. With both types of cross-linking, substantial improvements in interfacial shear strength, transverse tension strength, and transverse tension toughness are predicted. The results of this study are important for optimizing the processing of flCNT/AC composites for demanding engineering applications.

Computational Modeling of Hybrid Carbon Fiber/Epoxy Composites Reinforced with Functionalized and Non-Functionalized Graphene Nanoplatelets.

The mechanical properties of aerospace carbon fiber/graphene nanoplatelet/epoxy hybrid composites reinforced with pristine graphene nanoplatelets (GNP), highly concentrated graphene oxide (GO), and Functionalized Graphene Oxide (FGO) are investigated in this study. By utilizing molecular dynamics data from the literature, the bulk-level mechanical properties of hybrid composites are predicted using micromechanics techniques for different graphene nanoplatelet types, nanoplatelet volume fractions, nanoplatelet aspect ratios, carbon fiber volume fractions, and laminate lay-ups (unidirectional, cross-ply, and angle-ply). For the unidirectional hybrid composites, the results indicate that the shear and transverse properties are significantly affected by the nanoplatelet type, loading and aspect ratio. For the cross-ply and angle ply hybrid laminates, the effect of the nanoplate's parameters on the mechanical properties is minimal when using volume fractions and aspect ratios that are typically used experimentally. The results of this study can be used in the design of hybrid composites to tailor specific laminate properties by adjusting nanoplatelet parameters.

Prediction of the Interfacial Properties of High-Performance Polymers and Flattened CNT-Reinforced Composites Using Molecular Dynamics.

The next generation of ultrahigh-strength composites for structural components of vehicles for manned missions to deep space will likely incorporate flattened carbon nanotubes (flCNTs). With a wide range of high-performance polymers to choose from as the matrix component, efficient and accurate computational modeling can be used to efficiently downselect compatible resins and provide critical physical insight into the flCNT/polymer interface. In this study, molecular dynamics simulation is used to predict the interaction energy, frictional sliding resistance, and mechanical binding of flCNT/polymer interfaces for epoxy, bismaleimide (BMI), and benzoxazine high-performance resins. The results indicate that BMI has a stronger interfacial interaction and transverse tension binding with flCNT interfaces, while benzoxazine demonstrates the strongest levels of interfacial friction resistance.

Mechanical Properties and Characterization of Epoxy Composites Containing Highly Entangled As-Received and Acid Treated Carbon Nanotubes.

Huntsman-Merrimack MIRALON® carbon nanotubes (CNTs) are a novel, highly entangled, commercially available, and scalable format of nanotubes. As-received and acid-treated CNTs were added to aerospace grade epoxy (CYCOM® 977-3), and the composites were characterized. The epoxy resin is expected to infiltrate the network of the CNTs and could improve mechanical properties. Epoxy composites were tested for flexural and viscoelastic properties and the as-received and acid treated CNTs were characterized using Field-Emission Scanning and Transmission Electron Microscopy, X-Ray Photoelectron Spectroscopy, and Thermogravimetric Analysis. Composites containing 0.4 wt% as-received CNTs showed an increase in flexural strength, from 136.9 MPa for neat epoxy to 147.5 MPa. In addition, the flexural modulus increased from 3.88 GPa for the neat epoxy to 4.24 GPa and 4.49 GPa for the 2.0 wt% and 3.0 wt% as-received CNT/epoxy composites, respectively. FE-SEM micrographs indicated good dispersion of the CNTs in the as-received CNT/epoxy composites and the 10 M nitric acid 6 h treatment at 120 °C CNT/epoxy composites. CNTs treated with 10 M nitric acid for 6 h at 120 °C added oxygen containing functional groups (C-O, C=O, and O=C-O) and removed iron catalyst present on the as-received CNTs, but the flexural properties were not improved compared to the as-received CNT/epoxy composites.

Multiscale Modeling of Epoxy-Based Nanocomposites Reinforced with Functionalized and Non-Functionalized Graphene Nanoplatelets.

The impact on the mechanical properties of an epoxy resin reinforced with pristine graphene nanoplatelets (GNP), highly concentrated graphene oxide (GO), and functionalized graphene oxide (FGO) has been investigated in this study. Molecular dynamics (MD) using a reactive force field (ReaxFF) has been employed in predicting the effective mechanical properties of the interphase region of the three nanocomposite materials at the nanoscale level. A systematic computational approach to simulate the reinforcing nanoplatelets and probe their influence on the mechanical properties of the epoxy matrix is established. The modeling results indicate a significant degradation of the in-plane elastic Young's (decreased by ~89%) and shear (decreased by ~72.5%) moduli of the nanocomposite when introducing large amounts of oxygen and functional groups to the robust sp2 structure of the GNP. However, the wrinkled morphology of GO and FGO improves the nanoplatelet-matrix interlocking mechanism, which produces a significant improvement in the out-of-plane shear modulus (increased by 2 orders of magnitudes). The influence of the nanoplatelet content and aspect ratio on the mechanical response of the nanocomposites has also been determined in this study. Generally, the predicted mechanical response of the bulk nanocomposite materials demonstrates an improvement with increasing nanoplatelet content and aspect ratio. The results show good agreement with experimental data available from the literature.

Wetting Simulations of High-Performance Polymer Resins on Carbon Surfaces as a Function of Temperature Using Molecular Dynamics.

Resin/reinforcement wetting is a key parameter in the manufacturing of carbon nanotube (CNT)-based composite materials. Determining the contact angle between combinations of liquid resin and reinforcement surfaces is a common method for quantifying wettability. As experimental measurement of contact angle can be difficult when screening multiple high-performance resins with CNT materials such as CNT bundles or yarns, computational approaches are necessary to facilitate CNT composite material design. A molecular dynamics simulation method is developed to predict the contact angle of high-performance polymer resins on CNT surfaces dominated by aromatic carbon, aliphatic carbon, or a mixture thereof (amorphous carbon). Several resin systems are simulated and compared. The results indicate that the monomer chain length, chemical groups on the monomer, and simulation temperature have a significant impact on the predicted contact angle values on the CNT surface. Difunctional epoxy and cyanate ester resins show the overall highest levels of wettability, regardless of the aromatic/aliphatic nature of the CNT material surface. Tetrafunctional epoxy demonstrates excellent wettability on aliphatic-dominated surfaces at elevated temperatures. Bismaleimide and benzoxazine resins show intermediate levels of wetting, while typical molecular weights of polyether ether ketone demonstrate poor wetting on the CNT surfaces.

How to characterize interfacial load transfer in spiral carbon-based nanostructure-reinforced nanocomposites: is this a geometry-dependent process?

There is a great deal of attention given to spiral carbon-based nanostructures (SCBNs) because of their unique mechanical, thermal and electrical properties along with fascinating morphology. Dispersing SCBNs inside a polymer matrix leads to extraordinary properties of nanocomposites in diverse fields. However, the role of the interfacial mechanical properties of these nanocomposites remains unknown. Here, using molecular dynamics simulations, the characteristics of interfacial load transfer of SCBN-polyethylene nanocomposites are explored. Considering the geometric characteristics of SCBNs, new insight into the separation behavior of nanoparticles in normal and sliding modes is addressed. Interestingly, the results show that the maximum force and the separation energy of the SCBNs are much larger than those of graphene because of interlocking of the coils and polymer. The heavy influence of changes in the geometric characteristics of SCBNs on the separation behavior is observed. Pullout tests reveal that the influence of parameters such as the length and number of polyethylene chains, temperature, and functionalization of the SCBNs on the interfacial mechanical properties is also significant. This study sheds new light in understanding the crucial effect of the interaction of SCBNs with polymer chains on the interfacial mechanical properties, which can lead to better performance of nanocomposites.

Modeling Skeletal Muscle Stress and Intramuscular Pressure: A Whole Muscle Active-Passive Approach.

Clinical treatments of skeletal muscle weakness are hindered by a lack of an approach to evaluate individual muscle force. Intramuscular pressure (IMP) has shown a correlation to muscle force in vivo, but patient to patient and muscle to muscle variability results in difficulty of utilizing IMP to estimate muscle force. The goal of this work was to develop a finite element model of whole skeletal muscle that can predict IMP under passive and active conditions to further investigate the mechanisms of IMP variability. A previously validated hypervisco-poroelastic constitutive approach was modified to incorporate muscle activation through an inhomogeneous geometry. Model parameters were optimized to fit model stress to experimental data, and the resulting model fluid pressurization data were utilized for validation. Model fitting was excellent (root-mean-square error or RMSE <1.5 kPa for passive and active conditions), and IMP predictive capability was strong for both passive (RMSE 3.5 mmHg) and active (RMSE 10 mmHg at in vivo lengths) conditions. Additionally, model fluid pressure was affected by length under isometric conditions, as increases in stretch yielded decreases in fluid pressurization following a contraction, resulting from counteracting Poisson effects. Model pressure also varied spatially, with the highest gradients located near aponeuroses. These findings may explain variability of in vivo IMP measurements in the clinic, and thus help reduce this variability in future studies. Further development of this model to include isotonic contractions and muscle weakness would greatly benefit this work.

Simulating the effects of carbon nanotube continuity and interfacial bonding on composite strength and stiffness.

Molecular dynamics simulations of carbon nanotube (CNT) composites, in which the CNTs are continuous across the periodic boundary, overestimate the experimentally measured mechanical properties of CNT composites along the fiber direction. Since the CNTs in these composites are much shorter than the composite dimensions, load must be transferred either directly between CNTs or through the matrix, a mechanism that is absent in simulations of effectively continuous CNTs. In this study, the elastic and fracture properties of high volume fraction discontinuous carbon nanotube/amorphous carbon composite systems were compared to those of otherwise equivalent continuous CNT composites using ReaxFF reactive molecular dynamics simulations. These simulations were used to show how the number of nanotube-matrix interfacial covalent bonds affect composite mechanical properties. Furthermore, the mechanical impact of interfacial bonding was decomposed to reveal its effect on the properties of the CNTs, the interfacial layer of matrix, and the bulk matrix. For the composites with continuous reinforcement, it was found that any degree of interfacial bonding has a negative impact on axial tensile strength and stiffness. This is due to disruption of the structure of the CNTs and interfacial matrix layer by the interfacial bonds. For the discontinuous composites, the modulus was maximized between 4%-7% interfacial bonding and the strength continues to increase up to the highest levels of interfacial bonding studied. Areas of low stress and voids were observed in the simulated discontinuous composites at the ends of the tubes, from which fracture was observed to initiate. Experimental carbon nanotube yarn composites were fabricated and tested. The results were used to illustrate knockdown factors relative to the mechanical performance of the tubes themselves.