yade-users team mailing list archive

[Václav Šmilauer <eudoxos@xxxxxxxx>]

> In these Kn and Ks, E and v are macroscopic properties.
What ever paper you use, macroscopic and contact-level properties are
two different things, although they might be equal in some special
cases. 

For your code, it looks normal; comment out all non-elastic things
(including plasticity!!) and run the test again, i.e. you will have just
something like


 Real displN=geom->displacementN();
 Real Fn=phys->kn*displN;
 phys->normalForce=Fn*geom->normal;
 Vector3r Fs=phys->ks*geom->displacementT();
 applyForceAtContactPoint(....)

BTW it is quite confusing to call Fs "trialFs", since trial force (or
trial stress) is coloquially one before you do stress return in case of
plastic flow (you "try" if it is outside or inside, hence "trial"):

applyForceAtContactPoint(phys->normalForce+trialFs,geom->contactPoint,contact->getId1(),geom->se31.position,contact->getId2(),geom->se32.position,scene);


BTW2, at this line:

	Real crossSection=sqrt(2)*pow(geom->refR1,2);

order of id1 and id2 is arbitrary, and here you use just refR1 asymetrically.

HTH, v