Differential action model for contact inverse dynamics in multibody systems. More...
#include <contact-invdyn.hpp>
Classes | |
class | ResidualModelActuation |
Actuation residual. More... | |
class | ResidualModelContact |
Contact-acceleration residual. More... | |
Public Types | |
typedef ActuationModelAbstractTpl< Scalar > | ActuationModelAbstract |
typedef DifferentialActionModelAbstractTpl< Scalar > | Base |
typedef ConstraintModelManagerTpl< Scalar > | ConstraintModelManager |
typedef ContactItemTpl< Scalar > | ContactItem |
typedef ContactModelMultipleTpl< Scalar > | ContactModelMultiple |
typedef CostModelSumTpl< Scalar > | CostModelSum |
typedef DifferentialActionDataContactInvDynamicsTpl< Scalar > | Data |
typedef DataCollectorAbstractTpl< Scalar > | DataCollectorAbstract |
typedef DifferentialActionDataAbstractTpl< Scalar > | DifferentialActionDataAbstract |
typedef MathBaseTpl< Scalar > | MathBase |
typedef MathBase::MatrixXs | MatrixXs |
typedef StateMultibodyTpl< Scalar > | StateMultibody |
typedef MathBase::VectorXs | VectorXs |
![]() | |
typedef DifferentialActionDataAbstractTpl< Scalar > | DifferentialActionDataAbstract |
typedef MathBaseTpl< Scalar > | MathBase |
typedef MathBase::MatrixXs | MatrixXs |
typedef ScalarSelector< Scalar >::type | ScalarType |
typedef StateAbstractTpl< Scalar > | StateAbstract |
typedef MathBase::VectorXs | VectorXs |
Public Member Functions | |
DifferentialActionModelContactInvDynamicsTpl (std::shared_ptr< StateMultibody > state, std::shared_ptr< ActuationModelAbstract > actuation, std::shared_ptr< ContactModelMultiple > contacts, std::shared_ptr< CostModelSum > costs) | |
Initialize the contact inverse-dynamics action model. More... | |
DifferentialActionModelContactInvDynamicsTpl (std::shared_ptr< StateMultibody > state, std::shared_ptr< ActuationModelAbstract > actuation, std::shared_ptr< ContactModelMultiple > contacts, std::shared_ptr< CostModelSum > costs, std::shared_ptr< ConstraintModelManager > constraints) | |
Initialize the contact inverse-dynamics action model. More... | |
virtual void | calc (const std::shared_ptr< DifferentialActionDataAbstract > &data, const Eigen::Ref< const VectorXs > &x) override |
virtual void | calc (const std::shared_ptr< DifferentialActionDataAbstract > &data, const Eigen::Ref< const VectorXs > &x, const Eigen::Ref< const VectorXs > &u) override |
Compute the system acceleration, cost value and constraint residuals. More... | |
virtual void | calcDiff (const std::shared_ptr< DifferentialActionDataAbstract > &data, const Eigen::Ref< const VectorXs > &x) override |
virtual void | calcDiff (const std::shared_ptr< DifferentialActionDataAbstract > &data, const Eigen::Ref< const VectorXs > &x, const Eigen::Ref< const VectorXs > &u) override |
Compute the derivatives of the dynamics, cost and constraint functions. More... | |
template<typename NewScalar > | |
DifferentialActionModelContactInvDynamicsTpl< NewScalar > | cast () const |
Cast the contact-invdyn model to a different scalar type. More... | |
virtual bool | checkData (const std::shared_ptr< DifferentialActionDataAbstract > &data) override |
Checks that a specific data belongs to the contact inverse-dynamics model. | |
virtual std::shared_ptr< DifferentialActionDataAbstract > | createData () override |
Create the contact inverse-dynamics data. More... | |
EIGEN_MAKE_ALIGNED_OPERATOR_NEW | CROCODDYL_DERIVED_CAST (DifferentialActionModelBase, DifferentialActionModelContactInvDynamicsTpl) typedef _Scalar Scalar |
const std::shared_ptr< ActuationModelAbstract > & | get_actuation () const |
Return the actuation model. | |
const std::shared_ptr< ConstraintModelManager > & | get_constraints () const |
Return the constraint model manager. | |
const std::shared_ptr< ContactModelMultiple > & | get_contacts () const |
Return the contact model. | |
const std::shared_ptr< CostModelSum > & | get_costs () const |
Return the cost model. | |
virtual const VectorXs & | get_g_lb () const override |
Return the lower bound of the inequality constraints. | |
virtual const VectorXs & | get_g_ub () const override |
Return the upper bound of the inequality constraints. | |
virtual std::size_t | get_ng () const override |
Return the number of inequality constraints. | |
virtual std::size_t | get_ng_T () const override |
Return the number of equality terminal constraints. | |
virtual std::size_t | get_nh () const override |
Return the number of equality constraints. | |
virtual std::size_t | get_nh_T () const override |
Return the number of equality terminal constraints. | |
pinocchio::ModelTpl< Scalar > & | get_pinocchio () const |
Return the Pinocchio model. | |
virtual void | print (std::ostream &os) const override |
Print relevant information of the contact inverse-dynamics model. More... | |
virtual void | quasiStatic (const std::shared_ptr< DifferentialActionDataAbstract > &data, Eigen::Ref< VectorXs > u, const Eigen::Ref< const VectorXs > &x, const std::size_t maxiter=100, const Scalar tol=Scalar(1e-9)) override |
Computes the quasic static commands. More... | |
![]() | |
DifferentialActionModelAbstractTpl (std::shared_ptr< StateAbstract > state, const std::size_t nu, const std::size_t nr=0, const std::size_t ng=0, const std::size_t nh=0, const std::size_t ng_T=0, const std::size_t nh_T=0) | |
Initialize the differential action model. More... | |
bool | get_has_control_limits () const |
Indicates if there are defined control limits. | |
std::size_t | get_nr () const |
Return the dimension of the cost-residual vector. | |
std::size_t | get_nu () const |
Return the dimension of the control input. | |
const std::shared_ptr< StateAbstract > & | get_state () const |
Return the state. | |
const VectorXs & | get_u_lb () const |
Return the control lower bound. | |
const VectorXs & | get_u_ub () const |
Return the control upper bound. | |
VectorXs | quasiStatic_x (const std::shared_ptr< DifferentialActionDataAbstract > &data, const VectorXs &x, const std::size_t maxiter=100, const Scalar tol=Scalar(1e-9)) |
void | set_g_lb (const VectorXs &g_lb) |
Modify the lower bound of the inequality constraints. | |
void | set_g_ub (const VectorXs &g_ub) |
Modify the upper bound of the inequality constraints. | |
void | set_u_lb (const VectorXs &u_lb) |
Modify the control lower bounds. | |
void | set_u_ub (const VectorXs &u_ub) |
Modify the control upper bounds. | |
Protected Attributes | |
VectorXs | g_lb_ |
Lower bound of the inequality constraints. | |
VectorXs | g_ub_ |
< Lower bound of the inequality constraints | |
std::size_t | ng_ |
< Upper bound of the inequality constraints | |
std::size_t | nh_ |
< Number of inequality constraints | |
std::size_t | nu_ |
< Number of equality constraints | |
std::shared_ptr< StateAbstract > | state_ |
< Control dimension | |
![]() | |
VectorXs | g_lb_ |
Lower bound of the inequality constraints. | |
VectorXs | g_ub_ |
Lower bound of the inequality constraints. | |
bool | has_control_limits_ |
std::size_t | ng_ |
Number of inequality constraints. | |
std::size_t | ng_T_ |
Number of inequality terminal constraints. | |
std::size_t | nh_ |
Number of equality constraints. | |
std::size_t | nh_T_ |
Number of equality terminal constraints. | |
std::size_t | nr_ |
Dimension of the cost residual. | |
std::size_t | nu_ |
Control dimension. | |
std::shared_ptr< StateAbstract > | state_ |
Model of the state. | |
VectorXs | u_lb_ |
Lower control limits. | |
VectorXs | u_ub_ |
Upper control limits. | |
VectorXs | unone_ |
Neutral state. | |
Additional Inherited Members | |
![]() | |
EIGEN_MAKE_ALIGNED_OPERATOR_NEW typedef _Scalar | Scalar |
![]() | |
void | update_has_control_limits () |
Update the status of the control limits (i.e. if there are defined limits) | |
Differential action model for contact inverse dynamics in multibody systems.
This class implements forward kinematic with holonomic contact constraints (defined at the acceleration level) and inverse-dynamics computation using the Recursive Newton Euler Algorithm (RNEA). The stack of cost and constraint functions are implemented in CostModelSumTpl
and ConstraintModelManagerTpl
, respectively. The acceleration and contact forces are decision variables defined as the control inputs, and the under-actuation and contact constraint are under the name tau
and its frame name, thus the user is not allowed to use it.
Additionally, it is important to note that calcDiff()
computes the derivatives using the latest stored values by calc()
. Thus, we need to first run calc()
.
Definition at line 44 of file contact-invdyn.hpp.
DifferentialActionModelContactInvDynamicsTpl | ( | std::shared_ptr< StateMultibody > | state, |
std::shared_ptr< ActuationModelAbstract > | actuation, | ||
std::shared_ptr< ContactModelMultiple > | contacts, | ||
std::shared_ptr< CostModelSum > | costs | ||
) |
Initialize the contact inverse-dynamics action model.
It describes the kinematic evolution of the multibody system with contacts, and computes the needed torques using inverse-dynamics.
[in] | state | State of the multibody system |
[in] | actuation | Actuation model |
[in] | contacts | Multiple contacts |
[in] | costs | Cost model |
DifferentialActionModelContactInvDynamicsTpl | ( | std::shared_ptr< StateMultibody > | state, |
std::shared_ptr< ActuationModelAbstract > | actuation, | ||
std::shared_ptr< ContactModelMultiple > | contacts, | ||
std::shared_ptr< CostModelSum > | costs, | ||
std::shared_ptr< ConstraintModelManager > | constraints | ||
) |
Initialize the contact inverse-dynamics action model.
[in] | state | State of the multibody system |
[in] | actuation | Actuation model |
[in] | contacts | Multiple contacts |
[in] | costs | Cost model |
[in] | constraints | Constraints model |
|
overridevirtual |
Compute the system acceleration, cost value and constraint residuals.
It extracts the acceleration value from control vector and also computes the cost and constraints.
[in] | data | Contact inverse-dynamics data |
[in] | x | State point \(\mathbf{x}\in\mathbb{R}^{ndx}\) |
[in] | u | Control input \(\mathbf{u}\in\mathbb{R}^{nu}\) |
Implements DifferentialActionModelAbstractTpl< _Scalar >.
|
overridevirtual |
Reimplemented from DifferentialActionModelAbstractTpl< _Scalar >.
|
overridevirtual |
Compute the derivatives of the dynamics, cost and constraint functions.
It computes the partial derivatives of the dynamical system and the cost and contraint functions. It assumes that calc()
has been run first. This function builds a quadratic approximation of the time-continuous action model (i.e., dynamical system, cost and constraint functions).
[in] | data | Contact inverse-dynamics data |
[in] | x | State point \(\mathbf{x}\in\mathbb{R}^{ndx}\) |
[in] | u | Control input \(\mathbf{u}\in\mathbb{R}^{nu}\) |
Implements DifferentialActionModelAbstractTpl< _Scalar >.
|
overridevirtual |
Reimplemented from DifferentialActionModelAbstractTpl< _Scalar >.
|
overridevirtual |
Create the contact inverse-dynamics data.
Reimplemented from DifferentialActionModelAbstractTpl< _Scalar >.
DifferentialActionModelContactInvDynamicsTpl<NewScalar> cast | ( | ) | const |
Cast the contact-invdyn model to a different scalar type.
It is useful for operations requiring different precision or scalar types.
NewScalar | The new scalar type to cast to. |
|
overridevirtual |
Computes the quasic static commands.
The quasic static commands are the ones produced for a reference posture as an equilibrium point with zero acceleration, i.e., for \(\mathbf{f^q_x}\delta\mathbf{q}+\mathbf{f_u}\delta\mathbf{u}=\mathbf{0}\)
[in] | data | Contact inverse-dynamics data |
[out] | u | Quasic-static commands |
[in] | x | State point (velocity has to be zero) |
[in] | maxiter | Maximum allowed number of iterations (default 100) |
[in] | tol | Tolerance (default 1e-9) |
Reimplemented from DifferentialActionModelAbstractTpl< _Scalar >.
|
overridevirtual |
Print relevant information of the contact inverse-dynamics model.
[out] | os | Output stream object |
Reimplemented from DifferentialActionModelAbstractTpl< _Scalar >.