yade-dev team mailing list archive

[Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>]

> All I would need is the relativeVelocity to compute incidentVn and 
> incidentVs, having the normal. If I inherit from ScGeom, 
> relativeVelocity is not a member of it. So how would you access it?
class ScGeomWithVel: public ScGeom: {
...
YADE_CLASS_BASE_DOC_ATTRS_INIT_CTOR_PY(ScGeomWithVel,ScGeom,....
       ((Vector3r,incidentV,"incident velocity. |ycomp|"))
}

Vector3r ScGeomWithVel:::rotateAndGetShear(){//I'd have to make it 
virtual in ScGeom too, not a problem
   ...
   Vector3r relativeVelocity = (rbp2->vel+rbp2->angVel.cross(c2x)) - 
(rbp1->vel+rbp1->angVel.cross(c1x));
   relativeVelocity+=shiftVel;
   incidentV = relativeVelocity; 
<------------------------------------------ updating your variable here
   //keep the shear part only
   relativeVelocity = relativeVelocity-normal.dot(relativeVelocity)*normal;
   ...
}

Later in the constitutive law, you just access ig->incidentV.
Other suggestions?

Bruno