$$V(r) = 4\varepsilon\left[\left(\frac{\sigma}{r}\right)^{12}- \left(\frac{\sigma}{r}\right)^{6}\right]$$ Force: $F = -dV/dr$
Verlet integration:
$$x_{n+1}= 2x_n - x_{n-1}+ \frac{F_n}{m}\Delta t^2$$
Visualize metallic crystal lattice deformation with a spring-mass model. Adjust strain, lattice type, and defects to explore elastic deformation, dislocation nucleation, and plastic flow in real time.
The simulator models interatomic forces using the Lennard-Jones potential, which describes how the potential energy between two atoms changes with distance. It has a repulsive term at very close range and an attractive term at moderate range, creating a stable equilibrium.
$$ V(r) = 4\varepsilon\left[\left(\frac{\sigma}{r}\right)^{12}- \left(\frac{\sigma}{r}\right)^{6}\right] $$Here, $V(r)$ is the potential energy, $r$ is the distance between atoms, $\varepsilon$ is the depth of the potential well (bond strength), and $\sigma$ is the finite distance where the potential is zero. The $r^{-12}$ term is the strong repulsion, and the $r^{-6}$ term is the attraction.
The motion of each atom (mass) under the force from this potential is calculated using the Verlet integration algorithm, a common method in molecular dynamics. It updates an atom's position based on its current and previous positions and the acceleration from the net force.
$$ x_{n+1}= 2x_n - x_{n-1}+ \frac{F_n}{m}\Delta t^2 $$Here, $x_{n+1}$ is the new position, $x_n$ is the current position, $x_{n-1}$ is the previous position, $F_n$ is the net force on the atom at step $n$, $m$ is its mass, and $\Delta t$ is the time step. This equation allows the simulator to animate the lattice's dynamic response to deformation.
Alloy Design for Stronger Materials: By understanding how lattice type and defects control deformation, materials scientists design new alloys. For instance, adding carbon to iron (creating steel) pins dislocations in the BCC lattice, making it much harder for them to move and thus strengthening the material dramatically.
Predicting Metal Forming Limits: In automotive manufacturing, sheet metal is stamped into car body panels. CAE simulations using crystal plasticity models help predict where the metal might tear or wrinkle during this plastic forming process, optimizing the design and process parameters.
Fatigue Life Analysis in Aerospace: Cyclic loading on aircraft components can cause dislocations to accumulate and form cracks. Modeling deformation at the crystal level helps engineers predict fatigue life and schedule maintenance, preventing catastrophic failures.
Development of Nanoscale Devices: In micro-electromechanical systems (MEMS), components are so small that their mechanical behavior is governed by the motion of just a few dislocations. Atomic-scale deformation models are essential for designing reliable micro-sensors and actuators.
First, understand that this simulator does not perfectly replicate "the atom itself." It is strictly a model that "approximates interatomic bonds as springs." For example, in real metals, increased temperature intensifies atomic thermal vibration, making dislocations easier to move (facilitating workability), but this tool does not directly account for temperature effects. Think of those spring constants and potentials as implicitly representing "behavior at a certain fixed temperature."
Next, when increasing the "defect density" in the parameter settings, you might tend to think strength will always increase. While dislocation motion is indeed hindered, too many defects can conversely make the material brittle. In real material design, excessive precipitates can become crack initiation sites. In the simulator as well, if you apply tension with an extremely high number of defects, you should be able to observe fracture suddenly initiating from another location without dislocation motion. It's a great example to experience the strength-toughness trade-off.
Finally, be mindful of the simulation's "scale." What you're observing here is the nano- to micron-scale world. Predicting the deformation of a single automobile hood panel presents the next major challenge: "how to scale up these results to macroscopic behavior." It's important to understand that this tool is the first step, intended for learning the "fundamental mechanisms of deformation."
FCC copper lattice (a₀=3.61 Å, E≈110 GPa) under 5% uniaxial strain: spring constants k=10 N/m per bond. At 5% strain, edge displacements reach ~0.18 Å; maximum bond stress approaches 2.8 GPa. Introducing a single edge dislocation (b=2.56 Å) concentrates stress locally to 4.2 GPa within 3 lattice spacings. Plastic relaxation via atom shuffle reduces peak stress 15%.