Materials and Interaction Parameters

Theo_Score's picture
Submitted by Theo_Score on Thu, 02/25/2016 - 05:27

Hi All,

I just want to have clarity on this one.

We have 2 materials say m1 & m2 (peaches and oranges) and they are conveyed on conveyor belt which is boundary material, b1. We have 5 possible interactions here i.e (m1-m1), , (m2-m2), (m1-m2=m2-m1), (b1-m1) and (b1-m2).

The create_box command should have 2 atom-types? and how are the 5 interactions specified say coefficient of restitution?, I see the Fix property/global has specified 2 vector length and 4 material properties, what are the 4 material properties defining in relation to the 5 possible interactions I pointed to above.

Thank you.

richti83's picture

richti83 | Thu, 02/25/2016 - 17:50

You need 3 types in create box (m1,m2,b1), 3 global peratomtype, and 3*3=9 global material peratomtypepair values²

property/global coefficientFriction peratomtypepair 3 &
m1m1 m1m2 m1b1 &
m1m2 m2m2 m2b1 &
m1b1 m2b1 b1b1

where b1b1 is not used bc. there is no wall-wall interaction in liggghts (the value can't be left because of the way how the matrix is read in c++)

²=As you can see the matrix is symetrec bc. m1b1=b1m1 and so on, so there are 6 =n*(n-1) independent values.

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