Commit 0673f6f2 authored by Dan S. Bolintineanu's avatar Dan S. Bolintineanu
Browse files

Fixes a bug in pair granular that caused JKR cutoffs to be set incorrectly... 

Fixes a bug in pair granular that caused JKR cutoffs to be set incorrectly when multiple types with multiple cohesion values are used
parent ebee2b58

src/GRANULAR/pair_granular.cpp

+1 −1
Original line number Diff line number Diff line
@@ -1726,7 +1726,7 @@ double PairGranular::pulloff_distance(double radi, double radj, 
  double E, coh, a, Reff;

  Reff = radi*radj/(radi+radj);

  if (Reff <= 0) return 0;

  coh = normal_coeffs[itype][itype][3];

  coh = normal_coeffs[itype][jtype][3];

  E = normal_coeffs[itype][jtype][0]*THREEQUARTERS;

  a = cbrt(9*MY_PI*coh*Reff/(4*E));

  return a*a/Reff - 2*sqrt(MY_PI*coh*a/E);

CRA Git | Maintained and supported by SUSTech CRA and CCSE