Oh, you are right. It is a problem. I'll think about it.
You wrote the "throw" in the function in Dem3DofGeom; I assume you know
what it does... Unless you did cut&paste of the code, which you know
what I think about.
"throw" means "throw", and implies it has to be implemented in derived
class.
It IS implemented actually, in the class I was testing, I forgot
somebody else would use this code for facets too.
Is it really the case by the way? Who is using
Law2_ScGeom_FrictPhys_Basic with facets? Is it only working (I remember
it doesn't work for triaxial tests)?
(btw, a general question : why computing distances for energy? can't we
simply use cross products between force and disp. vectors?).
distance*|F|=work (|F| is already known, it is radius of plasticzone for
given sigmaN)? If you do cross product, then you have to compute the
norm of it anyway.
I wrote cross product?! Sorry, I meant dot product : u.f=work, there is
no norm().
Bruno
_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-dev
More help : https://help.launchpad.net/ListHelp