I am trying to use lammps to solve problems in contact mechanics. My system consists hemispherical asperity placed on top of a flat plate and this asperity will be deformed by a rigid plane. But first I want to apply the following heat treatment to this system: The system is to be maintained at constant temperature of 1350K for 0.25ns, then temperature drops to 1000K and is held constante for another 0.25ns; temperature drops again to 700K and is held constant for 0.25ns; new drop to 300K and another 0.25ns.
With the input file (file attached) I am using the simulation starts OK at 1350K but when the temperature drops to 1000K I get the error:
" Bad principal moments"
Can someone please help?
in.anneal (1.6 KB)