Peridynamics with LAMMPS

Hello!

I am working on the PMB Peridynamic model. I’have read some papers, but I don’t understand the physical meaning of the constant alpha which seems to have a link with breaks.
I have made some tests with different values for alpha but it changes nothing.
Can you enlight me about this?

Hugo

Again, Rezwan may advise.

Steve

alpha basically controls upped bound of the critical stretch s0. as its in the source code:

if (stretch > MIN(s0[i],s0[j])) partner[i][jj] = 0;

if (first)
s0_new[i] = s00[itype][jtype] - (alpha[itype][jtype] * stretch);
else
s0_new[i] = MAX(s0_new[i],s00[itype][jtype] - (alpha[itype][jtype] * stretch));

s0_new[i] is updated based on s00, alpha and bond stretch. with appropriate choice of alpha you should be able to see the effect. in the attached files stress (bar) at different timesteps are shown with fixed s00 and varying alpha: 1, 0.1, 0.001.

Again, Rezwan may advise.

Steve

alpha_1.png

alpha_0.1.png

alpha_0.001.png