Class NewtonEuler1DR

Defined in Program listing for file kernel/src/modelingTools/NewtonEuler1DR.hpp

class NewtonEuler1DR : public NewtonEulerR

NewtonEuler1DR.

This class is an interface for a relation with impact. It implements the computation of the jacoboian of h from the points of contacts and the normal. Use this class consists in overloading the method computeh, by setting the member pc1, pc2, nc and y. The matrix jachq is used both for the building of the OSNSP (with T) and for the predictor of activation of deactivation of the Interaction.

Subclassed by NewtonEuler3DR, NewtonEuler5DR

Public Functions

NewtonEuler1DR()

constructorx

virtual ~NewtonEuler1DR()

destructor

virtual void computeh(double time, BlockVector &q0, SiconosVector &y)

to compute y = h(q,v,t) using plug-in mechanism

Parameters
  • time: current time

  • q0: the position

  • y: the output

virtual void computeJachq(double time, Interaction &inter, SP::BlockVector q0)

compute the jacobian of h w.r.t.

q

Parameters
  • time: current time

  • inter: the interaction using this relation

  • q0: the container of the block vector to the dynamical system

double distance() const

Return the distance between pc1 and pc, with sign according to normal.

virtual void initialize(Interaction &inter)

initialize components specific to derived classes.

Parameters

void setRelNc(SP::SiconosVector nnc)

Set the coordinates of inside normal vector at the contact point in ds2 frame.

It will be used to compute _Nc during computeh().

Parameters
  • nnc: new coordinates

void setRelPc1(SP::SiconosVector npc)

Set the coordinates of first contact point in ds1 frame.

It will be used to compute _Pc1 during computeh().

Parameters
  • npc: new coordinates

void setRelPc2(SP::SiconosVector npc)

Set the coordinates of second contact point in ds2 frame It will be used to compute _Pc2 during computeh().

Parameters
  • npc: new coordinates

Public Members

bool _isOnContact

V.A.

boolean _isOnCOntact ?? Why is it public members ? seems parametrize the projection algorithm the projection is done on the surface \(y=0\) or on \(y \geq 0\)