[lammps-users] Metal cutting simulation-reg

Sir,
I'm using LAMMPS software to simulate metal cutting process and was
getting some kind of errors.

my basic idea is to move a steel/copper tool over a copper work piece
and perform the operation of mill/ turn/ shear at a particular
velocity and observe the chip formation and compute the various forces
and stresses acting on each atom and on the whole system.

I'm totally new to LAMMPS and molecular dynamics. i have written the
code after reading the user manual, but i know there may be some
errors in it. kindly help me in this regard.

apart from the above mentioned, i'm not able to understand the Axis
system used in lammps {whether it is XYZ (wherein XY plane being the
horizontal and YZ being vertical) or ZXY (wherein ZX plane being the
horizontal and XY being vertical) ? }

also, other than that i had problem setting up the tool with two
opposite faces tilted towards each other at an angle and the angle of
inclination is a variable for my study. (right now its just set up as
a prism with one face inclined 5 degrees.) please refer the attached
schematics for better understanding

I'm herewith attaching the code that i have written and the schematic
of what i want to simulate also the error message that i got when i
ran the simulation for your reference.

Vishnu Vardhan Chandrasekaran
Graduate Teaching and Research Assistant
Department of Mechanical Engineering
Auburn University
Auburn, AL. USA.

22.o342864 (8.3 KB)

orthogonal Machining.docx (23.6 KB)

Sorry, but I don't typically have time to study people's input
scripts. You'll need to ask a more specific question about
how a particular command works to do what you want.

But I will say that there is no concept of "vertical" in LAMMPS.
It is just the usual x,y,z coordinate system. You can place
your atoms and geometric objects within the simulation
box at whatever orientation and position you want. If you are
using periodic boundaries, then you should think about what
restrictions that imposes. E.g. you can't have a tilted face of a tool
and expect it to cross periodic boundaries.

Steve