Hybrid simulation (granular & peri)

Hi all,

I am trying to perform a hybrid simulation with the pair_style granular and peri on LIGGGHTS.
For the moment, I did not receive reply from the forum of LIGGGHTS.

However, if anyone has already performed a hybrid simulation with the pair_style peri
Please, can he send me his script in order to correct mine ?

Thanks in advance
JF

My script:

units si
pair_style hybrid granular peri

pair_style hybrid /overlay gran/hooke/history cohesion sjkr lperi/pmb

atom_modify map array
boundary s s s
newton off

Read data

read_data hybrid.sys

small target

lattice sc 0.00005

Some grouping (useful for fixing)

group A type 1
group B type 2

communicate single vel yes

neighbor 0.001 bin
neigh_modify delay 0

timestep 0.00001

fix m1 all property/global youngsModulus peratomtype 1.e6 1.e6
fix m2 all property/global poissonsRatio peratomtype 0.3 0.3
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.7 0.7 0.7 0.7
fix m4 all property/global coefficientFriction peratomtypepair 2 0.85 0.85 0.85 0.85
fix m5 all property/global characteristicVelocity scalar 2.
fix m6 all property/global cohesionEnergyDensity peratomtypepair 2 500 500 500 500

#======= PERI ===============
pair_coeff * * 1.6863e22 0.0015001 0.0005 0.25
set group A density 2500
set group A volume 1.25e-10
velocity A set 0.0 0.0 0.0 sum no units box
fix integr_peri all nve
compute 1 A damage/atom

Neither your pair_style or pair_coeff commands are correct.
Why do you specify pair_style twice? The pair_coeff
commands require you to specify the sub-style. Both sub-styles
are not referenced. Please read the pair hybrid doc page.

Steve

Hi,

Thanks for your reply !!!

It works better like this:

pair_style hybrid/overlay gran/hertz/history cohesion sjkr peri/pmb

pair_coeff 1 1 peri/pmb 1.6863e22 0.0015001 0.0005 0.25

pair_coeff 2 2 gran/hertz/history

pair_coeff 1 2 gran/hertz/history

Regards
JF