1 / 18

Molecular dynamic simulation of thermodynamic and mechanical properties and behavior of materials when dynamic loading

Molecular dynamic simulation of thermodynamic and mechanical properties and behavior of materials when dynamic loading . V . V . Dremov, A . V . Karavaev, F . A . Sapozhnikov, M . A . Vorobyova, RFNC-VNIITF , RUSSIA, E-mail: V.V.Dryomov@vniitf.ru L. Soulard

audra
Download Presentation

Molecular dynamic simulation of thermodynamic and mechanical properties and behavior of materials when dynamic loading

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Molecular dynamic simulation of thermodynamic and mechanical properties and behavior of materials when dynamic loading V.V. Dremov,A.V. Karavaev, F.A. Sapozhnikov, M.A. Vorobyova, RFNC-VNIITF, RUSSIA,E-mail:V.V.Dryomov@vniitf.ru L. Soulard CEA/DAM, France, Laurent.Soulard@cea.fr

  2. INTRODUCTION • Classical MD approach has been applied to modeling Be properties and behavior when dynamic loading. • Special attention has been paid to calculation of melting curve and physical properties when melting. Hugoniostat MD technique was applied to obtain Hugoniot of beryllium taking melting into account. • The results of direct MD simulation of shock wave loading of nano-polycrystalline beryllium (hcp grains, average grain size ~10nm) and the data on dynamic yield stress as depended on shock stress were obtained. • So as the length of Be samples used was about 0.2 μm only ultra-fast stage (time-scale ~20 ps) of the relaxation process behind shock front has been investigated.

  3. I. Parameters of Modified Embedded Atom Model [1] Dremov V.V., Karavaev A.V., Kutepov A. L.,Soulard L., AIP Conf. Proc., 955 (2008), 305-308. [2] Baskes M. I. and Johnson R. A., Modelling Simul.Mater. Sci. Eng., 2 (1994), 147-163.

  4. Thermodynamic and Mechanical Properties of Be when Static Loading Isotherm, 300K Russian Federal Nuclear Centre – Institute of Technical Physics

  5. Thermodynamic and Mechanical Properties of Be when Static Loading Isobar, P=0 Russian Federal Nuclear Centre – Institute of Technical Physics

  6. STATIC MD CALCULATIONS OF MELTING CURVE Figure 1. Melting curve of Be

  7. STATIC MD CALCULATIONS OF HUGONIOTS Figure 2. Melting curves and Hugoniots in (P,T) plane

  8. Thermodynamic and Mechanic Properties of Be when Static Loading. Bulk Modulus B(), T=300K Russian Federal Nuclear Centre – Institute of Technical Physics

  9. Thermodynamic and Mechanic Properties of Be when Static Loading. Bulk Modulus B(,T) Russian Federal Nuclear Centre – Institute of Technical Physics

  10. STATIC MD CALCULATIONS Figure 3. Bulk modulus versus temperature at P = 0. Figure 4. Bulk modulus vs temperature for different pressures. This work MD data

  11. Thermodynamic and Mechanic Properties of Be when Static Loading. Shear Modulus (), T=300K Russian Federal Nuclear Centre – Institute of Technical Physics

  12. Thermodynamic and Mechanic Properties of Be when Static Loading. Shear Modulus (,T) Russian Federal Nuclear Centre – Institute of Technical Physics

  13. STATIC MD CALCULATIONS Figure 5. Shear modulus versus temperature at P = 0. Figure 6. Shear modulus vs temperature for different pressures. This work MD data.

  14. STATIC MD CALCULATIONS Figure 7. Young modulus versus temperature at P = 0. Figure 8. Young modulus vs temperature for different pressures. This work MD data.

  15. STATIC MD CALCULATIONS Figure 9. Sound speed on the Hugoniot: experimental points from [13] are black boxes (CL) and rhombs (CB); experimental points from [14] are triangles; ab inito MD results on CL [8] are shown by the solid line 1; longitudinal and bulk sound speeds from MD calculations are shown by dashed and solid (2) lines respectively.

  16. DIRECT MD MODELING OF BERYLLIUM RESPONSE TO SHOCK Plane wave shock loading was simulated with the use of Be mono-crystalline and nano-polycrystalline samples having length ~0.18 μm and crossection 7070 rectangular unit cells (rectangular HCP u.c. contains 4 atoms). Figure 10. Evolution of shear stress profile with time when monocrystalline loading (Up=2000 m/s) in direction. Curves 1-5 corresponds to different moments of time (8,10,12,14,16 ps). Arrows indicate to shock and rarefaction wave fronts and to direction of propagation of the waves

  17. DIRECT MD MODELING OF BERYLLIUM RESPONSE TO SHOCK

  18. CONCLUSION • Obtained results suggest that • the new potential gives a temperature of melting at P=0, which is close to the experimental 1550K; • the temperature of melting commencement on the Hugoniot is about 1500K lower than that from quantum MD (possible effect of the system size?); • all Hugoniot data agree well with each other in P-T coordinates; • bulk and longitudinal sound speeds along the Hugoniot were obtained in a wide pressure range with the use of (P-T) tabulated elastic moduli and the Hugoniot. The data agree well with the available experimental data; • MD investigation into the elastic-plastic properties of Be suggests that at times typical for the MD simulation and the level of loading about 10-40 GPa, the yield stress for both monocrystalline and polycrystalline samples is much higher than in experiments. More properly, we should say about the effective yield stress whose value is defined by relatively slow kinetics of elastic-plastic deformation.

More Related