Implementation of V model in LIGGGHTS
UNDER CONSTRUCTION
Soon: pictures, more detailed explanation, example run-file commands for liggghts.
„Implementation of the V Model in LIGGGHTS” is a joint project at SPbSTU with JKU.
Содержание
Supporters and Contributors
- Patrick Fodor, JKU
- Christoph Kloss, JKU
- Oleg Kovalev, SPbSTU
- Vitaly Kuzkin, SPbSTU
LIGGGHTS
LIGGGHTS is an Open Source Discrete Element Method Particle Simulation Software.
V-Model
The article is based on the paper V.A. Kuzkin and I.E. Asonov "Vector-based model of elastic bonds for simulation of granullar solids "
Implementing V-Model into LIGGGHTS
In order to implement V-Model into LIGGGHTS, necessary functions are split in 4 classes. One to create the bonds, one to calculate the forces and pass them to liggghts, one to propagate intact and remove broken bonds, and one to initialize an atomstyle containing all the properties necessary for each atom.
- Bond Creation: fix_bond_create_v.cpp/h
- connected to inputscript by “FixStyle(bond/create/v,FixBondCreateV)”
- Checking for Bonding conditions
- If Bonding conditions are fulfilled, bond will be created and initialized.
- Force Calculation: bond_v.cpp/h
- connected to inputscript by “BondStyle(vbond/BondV)”
- Loop through all bonds and calculate forces on particles.
- Bond Propagation: fix_bond_propagate_v.cpp/h
- connected to inputscript by “FixStyle(bond/propagate/v,FixBondPropagateV)”
- propagate bond contact history
- remove bonds
- Atom Properties: atom_vec_bond_v.cpp/h
- connected to inputscript by “AtomStyle(bond/v,AtomVecBondV)”
- contain atom properties (location, velocities,…)
Results
A new applicable Bondstyle based on the V Model has been succesfully implemented into liggghts framework and can be used via inputscripts just as other functions of liggghts.
Contacts
Patrick Fodor:
- Email: polyun.p@gmail.com
- Skype: polyun-
