NodeRigidBodyRotVecLG

A 3D rigid body node based on rotation vector and Lie group methods for rigid bodies; the node has 3 displacement coordinates and three rotation coordinates and can be used in combination with explicit Lie Group time integration methods.

Authors: Gerstmayr Johannes, Holzinger Stefan

Additional information for NodeRigidBodyRotVecLG:

  • This Node has/provides the following types = Position, Orientation, RigidBody, RotationRotationVector
  • Short name for Python = RigidRotVecLG
  • Short name for Python visualization object = VRigidRotVecLG

The item NodeRigidBodyRotVecLG with type = ‘RigidBodyRotVecLG’ has the following parameters:

  • name [type = String, default = ‘’]:
    node’s unique name
  • referenceCoordinates [qref=[q0,q1,q2,ν0,ν1,ν2]Tref=[pTref,νTref]T, type = Vector6D, size = 3, default = [0.,0.,0., 0.,0.,0.]]:
    reference coordinates (position and rotation vector ν) of node ==> e.g. ref. coordinates for finite elements or reference position of rigid body (e.g. for definition of joints)
  • initialCoordinates [qini=[q0,q1,q2,ν0,ν1,ν2]Tini=[uTini,νTini]T, type = Vector6D, size = 3, default = [0.,0.,0., 0.,0.,0.]]:
    initial displacement coordinates u and rotation vector ν relative to reference coordinates
  • initialVelocities [q˙ini=[q˙0,q˙1,q˙2,ν˙0,ν˙1,ν˙2]Tini=[u˙Tini,ν˙Tini]T, type = Vector6D, size = 3, default = [0.,0.,0., 0.,0.,0.]]:
    initial velocity coordinate: time derivatives of displacement and angular velocity vector
  • visualization [type = VNodeRigidBodyRotVecLG]:
    parameters for visualization of item

The item VNodeRigidBodyRotVecLG has the following parameters:

  • show [type = Bool, default = True]:
    set true, if item is shown in visualization and false if it is not shown
  • drawSize [type = float, default = -1.]:
    drawing size (diameter, dimensions of underlying cube, etc.) for item; size == -1.f means that default size is used
  • color [type = Float4, size = 4, default = [-1.,-1.,-1.,-1.]]:
    Default RGBA color for nodes; 4th value is alpha-transparency; R=-1.f means, that default color is used

DESCRIPTION of NodeRigidBodyRotVecLG

The following output variables are available as OutputVariableType in sensors, Get…Output() and other functions:

  • Position: p0config=[p0,p1,p2]0configT=u0config+p0ref
    global 3D position vector of node; uref=0
  • Displacement: u0config=[q0,q1,q2]configT
    global 3D displacement vector of node
  • Velocity: v0config=[q˙0,q˙1,q˙2]configT
    global 3D velocity vector of node
  • Acceleration: a0config=[q¨0,q¨1,q¨2]configT
    global 3D acceleration vector of node
  • Coordinates: cconfig=[q0,q1,q2,ν0,ν1,ν2]Tconfig
    coordinate vector of node, having 3 displacement coordinates and 3 Euler angles
  • Coordinates\_t: c˙config=[q˙0,q˙1,q˙2,ν˙0,ν˙1,ν˙2]Tconfig
    velocity coordinates vector of node
  • RotationMatrix: [A00,A01,A02,A10,,A21,A22]configT
    vector with 9 components of the rotation matrix A0bconfig in row-major format, in any configuration; the rotation matrix transforms local (b) to global (0) coordinates
  • Rotation: [φ0,φ1,φ2]Tconfig
    vector with 3 components of the Euler/Tait-Bryan angles in xyz-sequence (A0bconfig=:A0(φ0)A1(φ1)A2(φ2)), recomputed from rotation matrix
  • AngularVelocity: ω0config=[ω0,ω1,ω2]0configT
    global 3D angular velocity vector of node
  • AngularVelocityLocal: ωbconfig=[ω0,ω1,ω2]bconfigT
    local (body-fixed) 3D angular velocity vector of node

Detailed information: For a detailed description on the rigid body dynamics formulation using this node, see Holzinger and Gerstmayr .

The node has 3 displacement coordinates [q0,q1,q2]T and three rotation coordinates, which is the rotation vector

ν=φn=νconfig+νref,

with the rotation angle φ and the rotation axis n. All coordinates cconfig lead to second order differential equations, However the rotation vector cannot be used as a conventional parameterization. It must be computed within a nonlinear update, using appropriate Lie group methods. The first 3 equations are residuals of translational forces in global coordinates, while the last 3 equations are residual of local (body-fixed) torques, compare the equations of motion of the rigid body.

The rotation matrix A(ν)0bconfig transforms a local (body-fixed) 3D position bv=[b0,b1,b2]bT to global 3D positions,

b0config=A(ν)0bconfigbb

Note that A(ν) is defined in function `` RotationVector2RotationMatrix``, see Section Function: RotationVector2RotationMatrix.

A Lie group integrator must be used with this node, which is why the is used, the rotation parameter velocities are identical to the local angular velocity ωb and thus the matrix Gb becomes the identity matrix.

Note, that the node automatically switches to Lie group integration of its rotational coordinates, both in explicit integration as well as for implicit time integration. This node avoids typical singularities of rotations and is therefore perfectly suited for arbitrary motion. Furthermore, nonlinearities are reduced, which may improve implicit time integration performance.

For creating a NodeRigidBodyRotVecLG together with a rigid body, there is a rigidBodyUtilities function AddRigidBody, see Section Function: AddRigidBody, which simplifies the setup of a rigid body significantely!

Relevant Examples and TestModels with weblink:

The web version may not be complete. For details, consider also the Exudyn PDF documentation : theDoc.pdf