Remote interaction between particles

Submitted by klchungae on Mon, 01/04/2016 - 05:18

Hi everyone,

I want to simulate a remote interaction between particles and a bond's breakage situation when the tensile is over a certain limit. How can I do that?

Forums:

klchungae | Wed, 01/06/2016 - 15:01

It is really helpful. Thank a lot!
I have a further question.
Can ParaView demonstrate bondings? I want to view the bond as an object.

klchungae | Fri, 01/08/2016 - 12:27

I have another question. I have formed a 50*50 particles grid mesh by bonds and I want to fix one end of the mesh.
However, the freeze command doesn't work. The mesh keeps falling under the action of gravity.

atom_style hybrid granular bond/gran n_bondtypes 1 bonds_per_atom 6
atom_modify map array
boundary f f f
newton off
communicate single vel yes
units si

region reg block 0 2 0 2 0 2 units box
create_box 1 reg

neighbor 0.1 bin
neigh_modify delay 0

#New pair style
pair_style gran model hertz tangential history
pair_coeff * *

#bond
bond_style gran

bond_coeff 1 0.011 1e5 1e5 1 1e5 1e5

#Material properties required for new pair styles
fix m1 all property/global youngsModulus peratomtype 1.e7
fix m2 all property/global poissonsRatio peratomtype 0.3
fix m3 all property/global coefficientRestitution peratomtypepair 1 0.2
fix m4 all property/global coefficientFriction peratomtypepair 1 0.3

mass 1 1.0

#Create atoms
region bc block 0 2 0 2 1 1.005 units box
lattice sc 0.04
create_atoms 1 region bc

group bonded region bc
set group bonded density 6300 diameter 0.011
run 1

group freezed id <= 50

run 1

fix f1 freezed freeze

fix bondcr all bond/create/gran 1 1 1 0.041 1 6

timestep 0.00001

fix gravi all gravity 9.81 vector 0.0 0.0 -1.0

#region and insertion
group nve_group region reg

#apply nve integration to all particles that are inserted as single particles
fix integr nve_group nve/sphere

#output settings
thermo_style custom step atoms numbond
thermo 1000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes

#insert the first particles so that dump is not empty
shell mkdir post
dump dmp bonded custom 100 post/dump*.liggghts id type x y z vx vy vz fx fy fz omegax omegay omegaz radius

#insert particles
run 1
fix_modify bondcr every 0

run 20000 upto

richti83's picture

richti83 | Fri, 01/08/2016 - 12:59

as far as I remember the order of the fixes is important, try fix integr nve_group nve/sphere in front of fix f1 freezed freeze

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact

klchungae | Sun, 01/31/2016 - 16:43

I have downloaded libliggghts_bondreader.so from:
https://github.com/richti83/ParaView_Reader_for_LIGGGHTS/tree/master/pre...
, and used manage plugins to load new libliggghts_bondreader.so.

But error comes out when I open the bonds*.bond files.
Warning: In /home/kitware/Dashboards/MyTests/NightlyMaster/ParaViewSuperbuild-Release/paraview/src/paraview/VTK/Rendering/Core/vtkRenderer.cxx, line 1029
vtkOpenGLRenderer (0x2f16330): Resetting view-up since view plane normal is parallel