> how to compute twist and bending (i.e. relative rotation) increment between > two spheres?
Maybe using an equation similar to the total formulation (ScGeom.cpp:94) but with spin instead of rotations? It would give : relativeSpin = angularVelocity2 - angularVelocity1 twistVel = relativeSpin.dot(normal) bendVel = relativeSpin - twistVel and then you multiply by dt. Alternatively, you could do this at each timestep : ScGeom6D::precomputeRotations() ScGeom6D::initRotations() then it would use orientations of the previous timestep instead of initial orientations, and precomputed angles would be increments. It would really give increments, not vel*dt. Bruno _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : yade-dev@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp