794c8ec064
add sphere2.urdf move btSpatialAlgebra into LinearMath remove some warnings, introduce BT_ZERO, BT_ONE, BT_HALF as defines for 0.f/0., 1.f/1., 0.5f/0.5 respectively
230 lines
8.1 KiB
C++
230 lines
8.1 KiB
C++
/*
|
|
Bullet Continuous Collision Detection and Physics Library
|
|
Copyright (c) 2013 Erwin Coumans http://bulletphysics.org
|
|
|
|
This software is provided 'as-is', without any express or implied warranty.
|
|
In no event will the authors be held liable for any damages arising from the use of this software.
|
|
Permission is granted to anyone to use this software for any purpose,
|
|
including commercial applications, and to alter it and redistribute it freely,
|
|
subject to the following restrictions:
|
|
|
|
1. The origin of this software must not be misrepresented; you must not claim that you wrote the original software. If you use this software in a product, an acknowledgment in the product documentation would be appreciated but is not required.
|
|
2. Altered source versions must be plainly marked as such, and must not be misrepresented as being the original software.
|
|
3. This notice may not be removed or altered from any source distribution.
|
|
*/
|
|
|
|
#ifndef BT_MULTIBODY_LINK_H
|
|
#define BT_MULTIBODY_LINK_H
|
|
|
|
#include "LinearMath/btQuaternion.h"
|
|
#include "LinearMath/btVector3.h"
|
|
#include "BulletCollision/CollisionDispatch/btCollisionObject.h"
|
|
|
|
enum btMultiBodyLinkFlags
|
|
{
|
|
BT_MULTIBODYLINKFLAGS_DISABLE_PARENT_COLLISION = 1
|
|
};
|
|
|
|
#define BT_MULTIBODYLINK_INCLUDE_PLANAR_JOINTS
|
|
#define TEST_SPATIAL_ALGEBRA_LAYER
|
|
|
|
//
|
|
// Various spatial helper functions
|
|
//
|
|
|
|
//namespace {
|
|
|
|
#ifdef TEST_SPATIAL_ALGEBRA_LAYER
|
|
|
|
#include "LinearMath/btSpatialAlgebra.h"
|
|
|
|
#endif
|
|
//}
|
|
|
|
//
|
|
// Link struct
|
|
//
|
|
|
|
struct btMultibodyLink
|
|
{
|
|
|
|
BT_DECLARE_ALIGNED_ALLOCATOR();
|
|
|
|
btScalar m_mass; // mass of link
|
|
btVector3 m_inertiaLocal; // inertia of link (local frame; diagonal)
|
|
|
|
int m_parent; // index of the parent link (assumed to be < index of this link), or -1 if parent is the base link.
|
|
|
|
btQuaternion m_zeroRotParentToThis; // rotates vectors in parent-frame to vectors in local-frame (when q=0). constant.
|
|
|
|
btVector3 m_dVector; // vector from the inboard joint pos to this link's COM. (local frame.) constant. set for revolute joints only.
|
|
|
|
// m_eVector is constant, but depends on the joint type
|
|
// prismatic: vector from COM of parent to COM of this link, WHEN Q = 0. (local frame.)
|
|
// revolute: vector from parent's COM to the pivot point, in PARENT's frame.
|
|
btVector3 m_eVector;
|
|
|
|
#ifdef TEST_SPATIAL_ALGEBRA_LAYER
|
|
btSpatialMotionVector m_absFrameTotVelocity, m_absFrameLocVelocity;
|
|
#endif
|
|
|
|
enum eFeatherstoneJointType
|
|
{
|
|
eRevolute = 0,
|
|
ePrismatic = 1,
|
|
eSpherical = 2,
|
|
#ifdef BT_MULTIBODYLINK_INCLUDE_PLANAR_JOINTS
|
|
ePlanar = 3,
|
|
#endif
|
|
eFixed = 4,
|
|
eInvalid
|
|
};
|
|
|
|
|
|
|
|
// "axis" = spatial joint axis (Mirtich Defn 9 p104). (expressed in local frame.) constant.
|
|
// for prismatic: m_axesTop[0] = zero;
|
|
// m_axesBottom[0] = unit vector along the joint axis.
|
|
// for revolute: m_axesTop[0] = unit vector along the rotation axis (u);
|
|
// m_axesBottom[0] = u cross m_dVector (i.e. COM linear motion due to the rotation at the joint)
|
|
//
|
|
// for spherical: m_axesTop[0][1][2] (u1,u2,u3) form a 3x3 identity matrix (3 rotation axes)
|
|
// m_axesBottom[0][1][2] cross u1,u2,u3 (i.e. COM linear motion due to the rotation at the joint)
|
|
//
|
|
// for planar: m_axesTop[0] = unit vector along the rotation axis (u); defines the plane of motion
|
|
// m_axesTop[1][2] = zero
|
|
// m_axesBottom[0] = zero
|
|
// m_axesBottom[1][2] = unit vectors along the translational axes on that plane
|
|
#ifndef TEST_SPATIAL_ALGEBRA_LAYER
|
|
btVector3 m_axesTop[6];
|
|
btVector3 m_axesBottom[6];
|
|
void setAxisTop(int dof, const btVector3 &axis) { m_axesTop[dof] = axis; }
|
|
void setAxisBottom(int dof, const btVector3 &axis) { m_axesBottom[dof] = axis; }
|
|
void setAxisTop(int dof, const btScalar &x, const btScalar &y, const btScalar &z) { m_axesTop[dof].setValue(x, y, z); }
|
|
void setAxisBottom(int dof, const btScalar &x, const btScalar &y, const btScalar &z) { m_axesBottom[dof].setValue(x, y, z); }
|
|
const btVector3 & getAxisTop(int dof) const { return m_axesTop[dof]; }
|
|
const btVector3 & getAxisBottom(int dof) const { return m_axesBottom[dof]; }
|
|
#else
|
|
btSpatialMotionVector m_axes[6];
|
|
void setAxisTop(int dof, const btVector3 &axis) { m_axes[dof].m_topVec = axis; }
|
|
void setAxisBottom(int dof, const btVector3 &axis) { m_axes[dof].m_bottomVec = axis; }
|
|
void setAxisTop(int dof, const btScalar &x, const btScalar &y, const btScalar &z) { m_axes[dof].m_topVec.setValue(x, y, z); }
|
|
void setAxisBottom(int dof, const btScalar &x, const btScalar &y, const btScalar &z) { m_axes[dof].m_bottomVec.setValue(x, y, z); }
|
|
const btVector3 & getAxisTop(int dof) const { return m_axes[dof].m_topVec; }
|
|
const btVector3 & getAxisBottom(int dof) const { return m_axes[dof].m_bottomVec; }
|
|
#endif
|
|
|
|
int m_dofOffset, m_cfgOffset;
|
|
|
|
btQuaternion m_cachedRotParentToThis; // rotates vectors in parent frame to vectors in local frame
|
|
btVector3 m_cachedRVector; // vector from COM of parent to COM of this link, in local frame.
|
|
|
|
btVector3 m_appliedForce; // In WORLD frame
|
|
btVector3 m_appliedTorque; // In WORLD frame
|
|
|
|
btScalar m_jointPos[7];
|
|
btScalar m_jointTorque[6]; //TODO
|
|
|
|
class btMultiBodyLinkCollider* m_collider;
|
|
int m_flags;
|
|
|
|
|
|
int m_dofCount, m_posVarCount; //redundant but handy
|
|
eFeatherstoneJointType m_jointType;
|
|
|
|
// ctor: set some sensible defaults
|
|
btMultibodyLink()
|
|
: m_mass(1),
|
|
m_parent(-1),
|
|
m_zeroRotParentToThis(0, 0, 0, 1),
|
|
m_cachedRotParentToThis(0, 0, 0, 1),
|
|
m_collider(0),
|
|
m_flags(0),
|
|
m_dofCount(0),
|
|
m_posVarCount(0),
|
|
m_jointType(btMultibodyLink::eInvalid)
|
|
{
|
|
m_inertiaLocal.setValue(1, 1, 1);
|
|
setAxisTop(0, 0., 0., 0.);
|
|
setAxisBottom(0, 1., 0., 0.);
|
|
m_dVector.setValue(0, 0, 0);
|
|
m_eVector.setValue(0, 0, 0);
|
|
m_cachedRVector.setValue(0, 0, 0);
|
|
m_appliedForce.setValue( 0, 0, 0);
|
|
m_appliedTorque.setValue(0, 0, 0);
|
|
//
|
|
m_jointPos[0] = m_jointPos[1] = m_jointPos[2] = m_jointPos[4] = m_jointPos[5] = m_jointPos[6] = 0.f;
|
|
m_jointPos[3] = 1.f; //"quat.w"
|
|
m_jointTorque[0] = m_jointTorque[1] = m_jointTorque[2] = m_jointTorque[3] = m_jointTorque[4] = m_jointTorque[5] = 0.f;
|
|
}
|
|
|
|
// routine to update m_cachedRotParentToThis and m_cachedRVector
|
|
void updateCache()
|
|
{
|
|
//multidof
|
|
if (m_jointType == eRevolute)
|
|
{
|
|
m_cachedRotParentToThis = btQuaternion(getAxisTop(0),-m_jointPos[0]) * m_zeroRotParentToThis;
|
|
m_cachedRVector = m_dVector + quatRotate(m_cachedRotParentToThis,m_eVector);
|
|
} else
|
|
{
|
|
// m_cachedRotParentToThis never changes, so no need to update
|
|
m_cachedRVector = m_eVector + m_jointPos[0] * getAxisBottom(0);
|
|
}
|
|
}
|
|
|
|
void updateCacheMultiDof(btScalar *pq = 0)
|
|
{
|
|
btScalar *pJointPos = (pq ? pq : &m_jointPos[0]);
|
|
|
|
switch(m_jointType)
|
|
{
|
|
case eRevolute:
|
|
{
|
|
m_cachedRotParentToThis = btQuaternion(getAxisTop(0),-pJointPos[0]) * m_zeroRotParentToThis;
|
|
m_cachedRVector = m_dVector + quatRotate(m_cachedRotParentToThis,m_eVector);
|
|
|
|
break;
|
|
}
|
|
case ePrismatic:
|
|
{
|
|
// m_cachedRotParentToThis never changes, so no need to update
|
|
m_cachedRVector = m_dVector + quatRotate(m_cachedRotParentToThis,m_eVector) + pJointPos[0] * getAxisBottom(0);
|
|
|
|
break;
|
|
}
|
|
case eSpherical:
|
|
{
|
|
m_cachedRotParentToThis = btQuaternion(pJointPos[0], pJointPos[1], pJointPos[2], -pJointPos[3]) * m_zeroRotParentToThis;
|
|
m_cachedRVector = m_dVector + quatRotate(m_cachedRotParentToThis,m_eVector);
|
|
|
|
break;
|
|
}
|
|
#ifdef BT_MULTIBODYLINK_INCLUDE_PLANAR_JOINTS
|
|
case ePlanar:
|
|
{
|
|
m_cachedRotParentToThis = btQuaternion(getAxisTop(0),-pJointPos[0]) * m_zeroRotParentToThis;
|
|
m_cachedRVector = quatRotate(btQuaternion(getAxisTop(0),-pJointPos[0]), pJointPos[1] * getAxisBottom(1) + pJointPos[2] * getAxisBottom(2)) + quatRotate(m_cachedRotParentToThis,m_eVector);
|
|
|
|
break;
|
|
}
|
|
#endif
|
|
case eFixed:
|
|
{
|
|
m_cachedRotParentToThis = m_zeroRotParentToThis;
|
|
m_cachedRVector = m_dVector + quatRotate(m_cachedRotParentToThis,m_eVector);
|
|
|
|
break;
|
|
}
|
|
default:
|
|
{
|
|
//invalid type
|
|
btAssert(0);
|
|
}
|
|
}
|
|
}
|
|
};
|
|
|
|
|
|
#endif //BT_MULTIBODY_LINK_H
|