Equivalent thermo-mechanical parameters for perfect crystals

Thermo-elastic behavior of perfect single crystal is considered. The crystal is represented as a set of interacting particles (atoms). The approach for determination of equivalent continuum values for the discrete system is proposed. Averaging of equations of particles' motion and long wave approximation are used in order to make link between the discrete system and equivalent continuum. Basic balance equations for equivalent continuum are derived from microscopic equations. Macroscopic values such as Piola and Cauchy stress tensors and heat flux are represented via microscopic parameters. Connection between the heat flux and temperature is discussed. Equation of state in Mie-Gruneisen form connecting Cauchy stress tensor with deformation gradient and thermal energy is obtained from microscopic considerations.Comment: To be published in proceedings of IUTAM Simposium on "Vibration Analysis of Structures with Uncertainties", 2009; 14 pages

Lattice-switch Monte Carlo

We present a Monte Carlo method for the direct evaluation of the difference between the free energies of two crystal structures. The method is built on a lattice-switch transformation that maps a configuration of one structure onto a candidate configuration of the other by `switching' one set of lattice vectors for the other, while keeping the displacements with respect to the lattice sites constant. The sampling of the displacement configurations is biased, multicanonically, to favor paths leading to `gateway' arrangements for which the Monte Carlo switch to the candidate configuration will be accepted. The configurations of both structures can then be efficiently sampled in a single process, and the difference between their free energies evaluated from their measured probabilities. We explore and exploit the method in the context of extensive studies of systems of hard spheres. We show that the efficiency of the method is controlled by the extent to which the switch conserves correlated microstructure. We also show how, microscopically, the procedure works: the system finds gateway arrangements which fulfill the sampling bias intelligently. We establish, with high precision, the differences between the free energies of the two close packed structures (fcc and hcp) in both the constant density and the constant pressure ensembles.Comment: 34 pages, 9 figures, RevTeX. To appear in Phys. Rev.

Effect of stress-triaxiality on void growth in dynamic fracture of metals: a molecular dynamics study

The effect of stress-triaxiality on growth of a void in a three dimensional single-crystal face-centered-cubic (FCC) lattice has been studied. Molecular dynamics (MD) simulations using an embedded-atom (EAM) potential for copper have been performed at room temperature and using strain controlling with high strain rates ranging from 10^7/sec to 10^10/sec. Strain-rates of these magnitudes can be studied experimentally, e.g. using shock waves induced by laser ablation. Void growth has been simulated in three different conditions, namely uniaxial, biaxial, and triaxial expansion. The response of the system in the three cases have been compared in terms of the void growth rate, the detailed void shape evolution, and the stress-strain behavior including the development of plastic strain. Also macroscopic observables as plastic work and porosity have been computed from the atomistic level. The stress thresholds for void growth are found to be comparable with spall strength values determined by dynamic fracture experiments. The conventional macroscopic assumption that the mean plastic strain results from the growth of the void is validated. The evolution of the system in the uniaxial case is found to exhibit four different regimes: elastic expansion; plastic yielding, when the mean stress is nearly constant, but the stress-triaxiality increases rapidly together with exponential growth of the void; saturation of the stress-triaxiality; and finally the failure.Comment: 35 figures, which are small (and blurry) due to the space limitations; submitted (with original figures) to Physical Review B. Final versio