Multi-Scale Modelling of Plastic Deformation, Damage and Relaxation in Epoxy Resins.

Epoxy resin plasticity and damage was studied from molecular dynamic simulations and interpreted by the help of constitutive modelling. For the latter, we suggested a physically motivated approach that aims at interpolating two well-defined limiting cases; namely, pulling at the vanishing strain rate and very rapid deformation; here, taken as 50% of the speed of sound of the material. In turn, to consider 0.1-10-m/s-scale deformation rates, we employed a simple relaxation model featuring exponential stress decay with a relaxation time of 1.5 ns. As benchmarks, deformation and strain reversal runs were performed by molecular dynamic simulations using two different strain rates. Our analyses show the importance of molecular rearrangements within the epoxy network loops for rationalizing the strain-rate dependence of plasticity and residual stress upon strain reversal. To this end, our constitutive model reasonably reproduced experimental data of elastic and visco-elastic epoxy deformation, along with the maximum stress experienced before fracturing. Moreover, we show the importance of introducing damage elements for mimicking the mechanical behavior of epoxy resins.

Addressing Surface Effects at the Particle-Continuum Interface in a Molecular Dynamics and Finite Elements Coupled Multiscale Simulation Technique.

Atomistic-to-continuum coupling methods are used to unravel molecular mechanisms of polymers and polymer composites. These multiscale techniques advantageously combine the computational efficiency of continuum approaches while keeping the accuracy of particle-based methods. The Capriccio method [Pfaller et al. Comput. Methods Appl. Mech. Eng. 2013, 260, 109-129.] is a well-proven multiscale technique, which connects finite elements (FE) with molecular dynamics (MD) in a partitioned-domain approach. A vital aspect of these multiscale methods is to provide physically sound boundary conditions to the particle domain suppressing any interface effects at the domain boundary occurring due to the coupling. These interfacial coupling artifacts still pose a significant problem, especially for amorphous polymers due to their highly irregular microstructure. We solve this problem by extending the particle-continuum interface by a layer of passive atoms which move with the outer continuum, thereby providing the missing interactions with a surrounding polymer bulk to the inner particle region. This solution allows us to successfully reproduce structural and mechanical properties obtained under conventional periodic boundary conditions, like density, stress, Young's modulus, and Poisson's ratio. Furthermore, we demonstrate the application of a nonaffine deformation by means of a simple bending test. In general, our revised method provides a framework to apply complex deformations for molecular scientists, while it allows the engineering community to examine challenging phenomena such as fracture behavior at a molecular level.

Nonperiodic stochastic boundary conditions for molecular dynamics simulations of materials embedded into a continuum mechanics domain.

A scheme is described for performing molecular dynamics simulations on polymers under nonperiodic, stochastic boundary conditions. It has been designed to allow later the embedding of a particle domain treated by molecular dynamics into a continuum environment treated by finite elements. It combines, in the boundary region, harmonically restrained particles to confine the system with dissipative particle dynamics to dissipate energy and to thermostat the simulation. The equilibrium position of the tethered particles, the so-called anchor points, are well suited for transmitting deformations, forces and force derivatives between the particle and continuum domains. In the present work the particle scheme is tested by comparing results for coarse-grained polystyrene melts under nonperiodic and regular periodic boundary conditions. Excellent agreement is found for thermodynamic, structural, and dynamic properties.

Extensive CGMD Simulations of Atactic PS Providing Pseudo Experimental Data to Calibrate Nonlinear Inelastic Continuum Mechanical Constitutive Laws.

In this contribution, we present a characterization methodology to obtain pseudo experimental deformation data from CG MD simulations of polymers as an inevitable prerequisite to choose and calibrate continuum mechanical constitutive laws. Without restriction of generality, we employ a well established CG model of atactic polystyrene as exemplary model system and simulate its mechanical behavior under various uniaxial tension and compression load cases. To demonstrate the applicability of the obtained data, we exemplarily calibrate a viscoelastic continuum mechanical constitutive law. We conclude our contribution by a thorough discussion of the findings obtained in the numerical pseudo experiments and give an outline of subsequent research activities. Thus, this work contributes to the field of multiscale simulation methods and adds a specific application to the body of knowledge of CG MD simulations.

Investigation of interphase effects in silica-polystyrene nanocomposites based on a hybrid molecular-dynamics-finite-element simulation framework.

A recently developed hybrid method is employed to study the mechanical behavior of silica-polystyrene nanocomposites (NCs) under uniaxial elongation. The hybrid method couples a particle domain to a continuum domain. The region of physical interest, i.e., the interphase around a nanoparticle (NP), is treated at molecular resolution, while the surrounding elastic continuum is handled with a finite-element approach. In the present paper we analyze the polymer behavior in the neighborhood of one or two nanoparticle(s) at molecular resolution. The coarse-grained hybrid method allows us to simulate a large polymer matrix region surrounding the nanoparticles. We consider NCs with dilute concentration of NPs embedded in an atactic polystyrene matrix formed by 300 chains with 200 monomer beads. The overall orientation of polymer segments relative to the deformation direction is determined in the neighborhood of the nanoparticle to investigate the polymer response to this perturbation. Calculations of strainlike quantities give insight into the deformation behavior of a system with two NPs and show that the applied strain and the nanoparticle distance have significant influence on the deformation behavior. Finally, we investigate to what extent a continuum-based description may account for the specific effects occurring in the interphase between the polymer matrix and the NPs.