Thermomechanical Response of Fullerene-Reinforced Polymers by Coupling MD and FEM

The aim of the present study is to provide a computationally efficient and reliable hybrid numerical formulation capable of characterizing the thermomechanical behavior of nanocomposites, which is based on the combination of molecular dynamics (MD) and the finite element method (FEM). A polymeric material is selected as the matrix-specifically, the poly(methyl methacrylate) (PMMA) commonly known as Plexiglas due to its expanded applications. On the other hand, the fullerene C₂₄₀ is adopted as a reinforcement because of its high symmetry and suitable size. The numerical approach is performed at two scales. First, an analysis is conducted at the nanoscale by utilizing an appropriate nanocomposite unit cell containing the C₂₄₀ at a high mass fraction. A MD-only method is applied to accurately capture all the internal interfacial effects and accordingly its thermoelastic response. Then, a micromechanical, temperature-dependent finite element analysis takes place using a representative volume element (RVE), which incorporates the first-stage MD output, to study nanocomposites with small mass fractions, whose atomistic-only simulation would require a substantial computational effort. To demonstrate the effectiveness of the proposed scheme, numerous numerical results are presented while the investigation is performed in a temperature range that includes the PMMA glass transition temperature, T_g.