yade-users team mailing list archive

[Nasibeh Moradi <nasibeh.moradi@xxxxxxxxx>]

2010/4/30 Václav Šmilauer <eudoxos@xxxxxxxx>

>
> For your code, it looks normal;

now I don't know what I must do for solving it, what is your suggestion?

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.
>
you are right, thanks. :)

Nasibeh