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?
A project by DCS Computing and CFDEMresearch
This website uses cookies to ensure you get the best experience on our website. By continuing to browse the site you are agreeing to our use of cookies. This website uses cookies. By continuing to browse the site you are agreeing to our use of cookies. More info
richti83 | Mon, 01/04/2016 - 18:22
git pull special LIGGGHTS(R)
git pull special LIGGGHTS(R) Version with bondspackage enabled:
https://github.com/richti83/LIGGGHTS-WITH-BONDS
make it and have a look at /examples/bondspackage/beam/
richti83 | Mon, 01/04/2016 - 18:25
git clone special LIGGGHTS(R)
git clone special LIGGGHTS(R) Version with bondspackage enabled:
https://github.com/richti83/LIGGGHTS-WITH-BONDS
make it and have a look at /examples/bondspackage/beam/
klchungae | Wed, 01/06/2016 - 15:01
Hi richti83 !
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.
richti83 | Thu, 01/07/2016 - 12:50
Yes we can
See:
http://www.cfdem.com/comment/15310#comment-15310
klchungae | Fri, 01/08/2016 - 12:27
YOU ARE GOOD
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 | Fri, 01/08/2016 - 12:59
as far as I remember the
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
klchungae | Sun, 01/31/2016 - 16:43
Cannot visualize bondings
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