• Systèmes complexes et transition énergétique › Statistical physics in mechanics

• Laboratory of Physics and Chemistry of Surfaces and Interfaces • Service de Physique et Chimie des Surfaces et des Interfaces

• Laboratory of Nano-Objects and Complex Systems (LNOSC) • GMT-MSIN : Modélisation des Surfaces Interfaces et Nanostructures

A Study of Irreversible Deformation in Amorphous Silica: Using Molecular Dynamics Simulations

Cindy Lynn ROUNTREE, Elisabeth BOUCHAUD

*Coll: Mehdi TALAMALI, Damien VANDEMBROUCQ, (UMR CNRS/Saint-Gobain, Aubervilliers)
*

In order to understand the origin and nature of irreversible deformations in silica glass, we performed MD simulations of amorphous silica systems which are subjected to (1) a cyclic loading and unloading of the hydrostatic pressure and (2) a shearing force.

These MD simulations employ reliable interatomic potentials including steric repulsion, charge transfer, and electronic polarizability of atoms through pair-wise interaction terms. Covalent effects in silica are included through bond-bending and bond-stretching three-body terms. Comparing MD results with various structural, dynamic, and mechanical measurements validates the potential. Over the past decade, the LSU group also developed highly efficient, scalable, and portable multiresolution algorithms to carry out largescale MD simulations (10^{7}-10^{9} atoms) on parallel architectures and algorithms that allow real-time immersive visualization of up to a billion atoms. The simulations described here employ this suite of algorithms.

Simulations presented herein involve 12 independent 5184 atoms systems of amorphous silica. Amorphous silica was generated by heating β-cristobalite to 3200K and then quenching the molten system to 0K. At this point a conjugate gradient was run to relax the system. Afterwards the system was heated up to ambient temperature.

Subsequently, the effects of cyclic densification in 12 independent systems 5184 atoms each are examined. The systems are densified in a stepwise process to a peak point. Afterwards, the systems are expanded back out to (1) their initial volume; or (2) until the pressure is approximately 0GPa.

(Top) Image depicts the average residual densification verses the imposed maximum densification. The system’s residual densification departs from ~0% around 20% maximum densification. Thus onset of permanent plastic deformation occurs after the system has been relaxed from 20% maximum densification which is in agreement with experiment. (Bottom) Image depicts the average residual densification verses the pressure at the imposed maximum densification. The system’s residual densification departs from ~0% around 12GPa. Thus onset of permanent plastic deformation occurs after the system has been relaxed from 10GPa which is in agreement with experiment.

The evolution of the average stress of 12 systems, composed of 5184 atoms each, when sheared: 1) from 0% to 60% (error bars in gray); 2) from +5.91% to 0% (error bars in blue); 3) from +11.91% to 0% (error bars in yellow); 4) from +17.91% to 0% (error bars in orange); 5) from +29.91% to -60% (error bars in light green); and 6) from +41.91% t0 -60% (error bars in light blue). Negative deformation implies that the shear direction is opposite to the initial shear direction.

This process is repeated continuously by increasing slightly each time the value of the peak densification point. After the system reaches ~21% maximum densification, it was found to have a residual densification of ~0.6% (i.e. the system was unable to return to its original density, see figure 1). Similarly the pressure corresponding to peak densification was examined. Once the system exceeds ~12GPa, the system is found to have a residual densification (see figure 1). This 21% maximum densification and the 12GPa pressure define the onset of plasticity, in agreement with experimental values.

The effects of cyclic shearing in 12 independent systems 5184 atoms each are examined. The systems were initially sheared from 0% shear strain to 60% shear strain (see figure 2, black line, error bars in grey). As one can see up to about 5% shear strain the system experiences a linear elastic behavior. Also it is important to note that the shear modulus can be calculated in this linear region and is in good agreement with experiments. Afterwards the system enters into the non-linear elastic regime. The stress reaches a maximum around 15% at which point the stress begins to decrease slightly and the systems begins to flow.

#783 - Màj : 16/04/2007