Deforming glassy polystyrene: Influence of pressure, thermal history, and deformation mode on yielding and hardening

Confinement effects on the mechanical heterogeneity of polymer fiber glasses

Both polymer fiber glasses and bulk polymer glasses exhibit nonlinear mechanical responses under uniaxial deformation. In polymer fibers, however, polymer chains are confined strongly and the surface area is relatively large compared to their volume. The confinement and the surface may lead to the spatially heterogeneous relaxation of chains in polymer fibers. In this work we perform molecular dynamics simulations and investigate the relation between the heterogeneous dynamics and the nonlinear mechanical responses at a molecular level. Our molecular simulations capture successfully not only the nonlinear mechanical response but also the dependence of mechanical properties on the strain rate of typical polymer glasses as in experiments. We find that the local elastic modulus and the nonaffine displacement are spatially heterogeneous in the pre-yield regime, which results in a lower elastic modulus for polymer fibers than bulk polymer glasses. In the post-yield regime, those mechanical properties become relatively homogeneous. Monomers with large nonaffine displacement are localized mainly at the interfacial region in the pre-yield regime while highly nonaffine monomers are distributed throughout the fibers in the post-yield regime. We show that the nonaffine displacement during deformation relates closely to the mechanical response of the polymer fibers. We also find that in the strain-hardening regime there is a significant difference in the energetic contribution to the stress between polymer fibers and bulk polymers, for which the modulus of the strain-hardening regime of the polymer fibers is smaller than that of bulk polymers.

Mechanisms of reinforcement in polymer nanocomposites

Qualitatively different stress–strain responses of polymer nanocomposites are shown to result from the dynamical evolution of three principal molecular structural motifs in the polymer–filler network.

Molecular dynamics simulations of uniaxial deformation of thermoplastic polyimides

The results of atomistic molecular-dynamics simulations of mechanical properties of heterocyclic polymer subjected to uniaxial deformation are reported. A new amorphous thermoplastic polyimide R-BAPO with a repeat unit consisting of dianhydride 1,3-bis-(3',4,-dicarboxyphenoxy)diphenyl (dianhydride R) and diamine 4,4'-bis-(4''-aminophenoxy)diphenyloxide (diamine BAPO) was chosen for the simulations. Our primary goal was to establish the impact of various factors (sample preparation method, molecular mass, and cooling and deformation rates) on the elasticity modulus. In particular, we found that the elasticity modulus was only slightly affected by the degree of equilibration, the molecular mass and the size of the simulation box. This is most likely due to the fact that the main contribution to the elasticity modulus is from processes on scales smaller than the entanglement length. Essentially, our simulations reproduce the logarithmic dependence of the elasticity modulus on cooling and deformation rates, which is normally observed in experiments. With the use of the temperature dependence analysis of the elasticity modulus we determined the flow temperature of R-BAPO to be 580 K in line with the experimental data available. Furthermore, we found that the application of high external pressure to the polymer sample during uniaxial deformation can improve the mechanical properties of the polyimide. Overall, the results of our simulations clearly demonstrate that atomistic molecular-dynamics simulations represent a powerful and accurate tool for studying the mechanical properties of heterocyclic polymers and can therefore be useful for the virtual design of new materials, thereby supporting cost-effective synthesis and experimental research.

Polystyrene Glasses under Compression: Ductile and Brittle Responses

Polystyrene of different molecular weights and their binary mixtures are studied in terms of their various mechanical responses to uniaxial compression at different temperatures. PS of M_w = 25 kg/mol is completely brittle until it is above its glass transition temperature T_g. In contrast, upon incorporation of a high molecular weight component, PS mixtures turn from barely ductile a few degrees below its T_g to ductile over 40° below T_g. In the upper limit, a PS of M_w = 319 kg/mol yields and undergoes plastic flow, even at T = -70 °C. The observed dependence of mechanical responses on molecular weight and molecular weight distribution can be adequately rationalized by the idea that yielding and plastic compression are caused by chain networking.

Strain Hardening During Uniaxial Compression of Polymer Glasses

The origin of high mechanical stresses in large deformation of polymer glasses has been elusive because both plasticity and elasticity take place. In this work on the nature of the mechanical responses, we carry out uniaxial compression experiments to make simultaneous mechanical and thermal measurements of polycarbonate. Our results confirm that two factors contribute to the growing mechanical stress in the post-yield regime, which is known as "strain hardening". Besides plastic deformation that is intersegmental in origin, chain tension as an intrasegmental component contributes considerably to the measured stress in post-yield. Such a conclusion modifies the previous consensus regarding the nature of strain hardening in mechanical deformation of polymer glasses.

Exploring the ability of a multiscale coarse-grained potential to describe the stress-strain response of glassy polystyrene

A new particle-based bottom-up method to develop coarse-grained models of polymers is presented and applied to polystyrene. The multiscale coarse-graining (MS-CG) technique of Izvekov et al. [J. Chem. Phys. 120, 10896 (2004)] is applied to a polymer system to calculate nonbonded interactions. The inverse Boltzmann inversion method was used to parametrize the bonded and bond-angle bending interactions. Molecular dynamics simulations were performed, and the CG model exhibited a significantly lower modulus compared to the atomistic model at low temperature and high strain rate. In an attempt to improve the CG model performance, several other parametrization schemes were used to build other models from this base model. The first of these models included standard frictional forces through use of the constant-temperature dissipative particle dynamics method that improved the modulus, yet was not transferrable to higher temperatures and lower strain rates. Other models were built by increasing the attraction between CG beads through direct manipulation of the nonbonded potential, where an improvement of the stress response was found. For these models, two parametrization protocols that shifted the force to more attractive values were explored. The first protocol involved a uniform shift, while the other protocol shifted the force in a more localized region. The uniformly shifted potential greatly affected the structure of the equilibrium model as compared to the locally shifted potential, yet was more transferrable to different temperatures and strain rates. Further improvements in the coarse-graining protocol to generate models that more satisfactorily capture mechanical properties are suggested.

Elastic yielding in cold drawn polymer glasses well below the glass transition temperature

This Letter reports elastic-driven internal yielding in strained ductile polymer glasses. After cold drawing of two different polymer glasses to neck at room temperature, we show that the samples display considerable retractive stress when warmed up above the storage temperature but still considerably below their glass transition temperatures. We conclude that the elastic yielding arises from the distortion of backbones leading to intra-segmental tension in the chain network.

Atomistic picture of isothermal volume relaxation behavior of atactic polystyrene glass provided by a molecular dynamics simulation

Using an atomistic molecular dynamics (MD) simulation, we study the volume relaxation behavior of atactic polystyrene and attempt to correlate this macroscopic behavior with certain microscopic aspects. To this end, the gyration radius, the dimensionless relative shape anisotropy, the mean-squared displacement, and the non-Gaussian parameter are examined simultaneously. Our result shows that the structures characterized at different length scales change in a self-similar way and these changes are intimately correlated to the translational mobility of atoms. The initial incubation of structural changes at the different scales originates from the restriction of mobility due to the cage effect. The applicability of the MD simulation to the investigation of the bulk properties is discussed.

Glassy dynamics and mechanical response in dense fluids of soft repulsive spheres. II. Shear modulus, relaxation-elasticity connections, and rheology

We apply the quiescent and mechanically driven versions of nonlinear Langevin equation theory to study how particle softness influences the shear modulus, the connection between shear elasticity and activated relaxation, and nonlinear rheology of the repulsive Hertzian contact model of dense soft sphere fluids. Below the soft jamming threshold, the shear modulus follows a power law dependence on volume fraction over a narrow interval with an apparent exponent that grows with particle stiffness. To a first approximation, the elastic modulus and transient localization length are controlled by a single coupling constant determined by local fluid structure. In contrast to the behavior of hard spheres, an approximately linear relation between the shear modulus and activation barrier is predicted. This connection has recently been observed for microgel suspensions and provides a microscopic realization of the elastic shoving model. Yielding, shear and stress thinning of the alpha relaxation time and viscosity, and flow curves are also studied. Yield strains are relatively weakly dependent on volume fraction and particle stiffness. Shear thinning commences at values of the effective Peclet number far less than unity, a signature of stress-assisted activated relaxation when barriers are high. Apparent power law reduction of the viscosity with shear rate is predicted with a thinning exponent less than unity. In the vicinity of the soft jamming threshold, a power law flow curve occurs over an intermediate reduced shear rate range with an apparent exponent that decreases as fluid volume fraction and/or repulsion strength increase.

Evolution of entanglements during the response to a uniaxial deformation of lamellar triblock copolymers and polymer glasses

Using coarse-grained molecular-dynamics simulations, a generic styrene-block-butadiene-block-styrene triblock copolymer under lamellar conformation is used in order to investigate the mutual entanglement evolution when a structure of alternating glassy (S)/rubbery (B) layers is submitted to an imposed deformation. By varying the amount of loop chains between each phase, i.e., noncrossing chains, it is possible to generate different types of S/B interface definitions. A specific boundary driven tensile strain protocol has been developed in order to mimic "real" experiments and measure the stress-strain curve. The same protocol is also applied to a reference state consisting in a directed glassy homopolymers, as well as to an isotropic glassy polymer. The evolution of initial mutual entanglements from the undeformed samples during the whole deformation process is monitored. It is shown for all considered systems that initial entanglements mostly participate to the preyield regime of the stress-strain curve and that this network is debonded during the strain-hardening regime. For triblocks with a non-null amount of crossing chains, the lower the amount is, the longer the memory effect of the initial entanglement network in the postyield regime is. On the fly distributions of entanglements, which depart from the postyield regime, depict memory effects and long-time correlations during the strain-hardening regime. For triblocks, loop chains reinforce these effects.

From macroscopic yield criteria to atomic stresses in polymer glasses

The relationship between macroscopic shear yield criteria and local stress distributions in deformed polymer glasses is investigated via molecular dynamics simulations on different scales of coarse-graining. Macroscopic shear stresses at the yield point obey a pressure-modified von Mises (pmvM) criterion for many different loading conditions and strain rates. Average local stresses in small volume elements obey the same yield criterion for volumes containing approx. 100 atoms or more. Qualitatively different behavior is observed on smaller scales: the average octahedral atomic shear stress has a simple linear relationship to hydrostatic pressure regardless of macroscopic stress state and failure mode. Local plastic events are identified through a threshold in the mean-squared nonaffine displacement and compared to the local stress state. We find that the pmvM criterion only predicts local yield events when stress and displacements are averaged over at least 100 atoms. By contrast, macroscopic shear yield criteria appear to lose their ability to predict plastic activity on the atomic scale.

Segmental dynamics in polymers: from cold melts to ageing and stressed glasses

Recent progress in developing statistical mechanical theories of supercooled polymer melts and glasses is reviewed. The focus is on those approaches that are either explicitly formulated for polymers, or are applications of more generic theories to interpret polymeric phenomena. These include two configurational entropy theories, a percolated free volume distribution model, and the activated barrier hopping nonlinear Langevin theory. Both chemically-specific and universal aspects are discussed. After a brief summary of classic phenomenological approaches, a discussion of the relevant length scales and key experimental phenomena in both the supercooled liquid and glassy solid state is presented including ageing and nonlinear mechanical response. The central concepts that underlie the theories in the molten state are then summarized and key predictions discussed, including the glass transition in oriented polymer liquids and deformed rubber networks. Physical ageing occurs in the nonequilibrium glass, and theories for its consequences on the alpha relaxation are discussed. Very recent progress in developing a segment scale theory for the dramatic effects of external stress on polymer glasses, including acceleration of relaxation, yielding, plastic flow and strain hardening, is summarized. The article concludes with a discussion of outstanding theoretical challenges.

Fracture in glassy polymers: a molecular modeling perspective

Over the past 25 years, molecular modeling and simulations have provided important insights into the physics of deformation and fracture of glassy polymers. This review presents an overview of key results discussed in the context of experimentally observed polymer behavior. Both atomistic and coarse-grained polymer models have been used in different deformation protocols to study elastic properties, shear yielding, creep, physical aging, strain hardening and crazing. Simulations reproduce most of the macroscopic features of plasticity in polymer glasses such as stress-strain relations and creep response, and reveal information about the underlying atomistic processes. Trends of the shear yield stress with loading conditions, temperature and strain rate, and the atomistic dynamics under load have been systematically explored. Most polymers undergo physical aging, which leads to a history-dependent mechanical response. Simulations of strain hardening and crazing demonstrate the nature of polymer entanglements in the glassy state and the role of local plasticity and provide insight into the origin of fracture toughness of amorphous polymers.