Hi,
I am trying to understand why MD simulation crashes if larger time step (~100 fs) is applied to integrate the equation of motion. I know that there are two major things play role in solving F = ma and obtaining equilibrium solution. One is the accuracy, and the other is numerical stability of the integration algorithm.
Head responses subjected to impact loading are studied using the finite element method. The dynamic responses of the stress, strain, strain energy density and the intracranial pressure govern the intracranial tissues and skull material failures, and therefore, the traumatic injuries.
Welcome to join our research team in the Extreme Loading and Design, School of Mechanical, Aerospace and Civil Engineering, University of Manchester, Manchester, England !
Upon hearing the painful news of VPI tragedy, I am deeply saddened that we lost our young students and fellow faculty, in particular, Liviu Librescu, G. V. Loganathan and Kevin Granata. I would like to share the sorrow with the family and friends of the victims who lost their lives in the senseless crime. I am deeply compelled to write this message, since I am a faculty teaching young students applied mechanics/engineering like professors Librescu, Loganathan and Granata, and belong to a small community of Korean Americans. I feel helpless and guilty that I could not mentor our young people and hear such incomprehensible tragic news.
I would like to introduce my recent paper on field projection method with Prof. K.-S. Kim at Brown University. When atomistic simulations or atomic resolution experiments are carried out to analyze deformation energetics of atomic lattices near a crack tip, the results of the atomic positions and the total energy of the system have not been systematically converted to field quantities such as stresses or tractions near the crack tip. In this study, we introduce a mathematically consistent way of defining and measuring the cohesive tractions, separations and surface stresses in an atomic decohesion process zone using the atomic displacement data at some distance away from the crack tip. The method is called “nano-scale planar field projection method.” This is a generalization of the method developed for isotropic homogeneous solids by Hong and Kim (2003). The formulation introduced here is applicable to interface cracks between anisotropic solids as well as to cracks in homogeneous solids whether they are anisotropic or isotropic. The field projection method is then applied to a crack tip field in gold, simulated atomistically. The atomistic simulation is made with an embedded atom method (EAM) potential for a crystal decohesion along [112] direction in a (111) plane. Then, the details of energy partition in various modes of nano-scale separation processes are analyzed with the field projection method.
During the shootings occurred on the campus of Virginia Tech, the department of engineering science and mechanics lost two fine professors. Dr. Liviu Librescu was teaching solid mechanics class when the gunman initiated a two-minute massacre next door. Most of the students in his class escaped from the windows while Dr. Librescu held the door and did not survive when the gunman broke in. Another professor is Dr. Kevin Granata.
With a great many ties to VT, I have been following the tragic shootings there closely. While my former student (Nakhiah Goulbourne) and my former department head (Dick Benson) are both safe, I am afraid we may have lost one of our own this morning. From an AP report, a student stated that the instructor teaching a 9:05am mechanics class in 204 Norris Hall was killed. According to the timetable of classes at VT, this instructor would have been Liviu Librescu. I certainly hope the AP's report is wrong. As someone very interested in shell theories, I have spoken with Prof. Librescu many times and followed his work. His loss would be a great loss to our community.
