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.

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.

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.

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.

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.