robotoc
robotoc - efficient ROBOT Optimal Control solvers
Loading...
Searching...
No Matches
mpc_pace.hpp
Go to the documentation of this file.
1#ifndef ROBOTOC_MPC_PACE_HPP_
2#define ROBOTOC_MPC_PACE_HPP_
3
4#include <memory>
5
6#include "Eigen/Core"
7#include "Eigen/Geometry"
8
10#include "robotoc/ocp/ocp.hpp"
32
33
34namespace robotoc {
35
40class MPCPace {
41public:
48 MPCPace(const Robot& quadruped_robot, const double T, const int N);
49
54
59
63 MPCPace(const MPCPace&) = default;
64
68 MPCPace& operator=(const MPCPace&) = default;
69
73 MPCPace(MPCPace&&) noexcept = default;
74
78 MPCPace& operator=(MPCPace&&) noexcept = default;
79
88 void setGaitPattern(const std::shared_ptr<ContactPlannerBase>& foot_step_planner,
89 const double swing_height, const double swing_time,
90 const double stance_time, const double swing_start_time);
91
101 void init(const double t, const Eigen::VectorXd& q, const Eigen::VectorXd& v,
102 const SolverOptions& solver_options);
103
108 void reset();
109
114 void reset(const Eigen::VectorXd& q, const Eigen::VectorXd& v);
115
120 void setSolverOptions(const SolverOptions& solver_options);
121
131 void updateSolution(const double t, const double dt, const Eigen::VectorXd& q,
132 const Eigen::VectorXd& v);
133
138 const Eigen::VectorXd& getInitialControlInput() const;
139
144 const Solution& getSolution() const;
145
151
157 ControlPolicy getControlPolicy(const double t) const {
158 return ControlPolicy(ocp_solver_, t);
159 }
160
169 double KKTError(const double t, const Eigen::VectorXd& q,
170 const Eigen::VectorXd& v);
171
177 double KKTError() const;
178
183 std::shared_ptr<CostFunction> getCostHandle();
184
189 std::shared_ptr<ConfigurationSpaceCost> getConfigCostHandle();
190
195 std::shared_ptr<ConfigurationSpaceCost> getBaseRotationCostHandle();
196
201 std::vector<std::shared_ptr<TaskSpace3DCost>> getSwingFootCostHandle();
202
207 std::shared_ptr<CoMCost> getCoMCostHandle();
208
213 std::shared_ptr<Constraints> getConstraintsHandle();
214
219 std::shared_ptr<FrictionCone> getFrictionConeHandle();
220
225 const OCPSolver& getSolver() const { return ocp_solver_; }
226
231 const std::shared_ptr<ContactSequence>& getContactSequence() const {
232 return contact_sequence_;
233 }
234
240 void setRobotProperties(const RobotProperties& properties);
241
242 EIGEN_MAKE_ALIGNED_OPERATOR_NEW
243
244private:
245 std::shared_ptr<ContactPlannerBase> foot_step_planner_;
246 std::shared_ptr<ContactSequence> contact_sequence_;
247 std::shared_ptr<CostFunction> cost_;
248 std::shared_ptr<Constraints> constraints_;
249 OCPSolver ocp_solver_;
250 SolverOptions solver_options_;
251 ContactStatus cs_standing_, cs_left_swing_, cs_right_swing_;
253 double swing_height_, swing_time_, stance_time_, swing_start_time_,
254 T_, dt_, dtm_, ts_last_, eps_;
255 int N_, current_step_, predict_step_;
256 bool enable_stance_phase_;
257
258 std::shared_ptr<ConfigurationSpaceCost> config_cost_;
259 std::shared_ptr<ConfigurationSpaceCost> base_rot_cost_;
260 std::shared_ptr<TaskSpace3DCost> LF_foot_cost_, LH_foot_cost_,
261 RF_foot_cost_, RH_foot_cost_;
262 std::shared_ptr<CoMCost> com_cost_;
263 std::shared_ptr<MPCPeriodicConfigurationRef> base_rot_ref_;
264 std::shared_ptr<MPCPeriodicSwingFootRef> LF_foot_ref_, LH_foot_ref_,
265 RF_foot_ref_, RH_foot_ref_;
266 std::shared_ptr<MPCPeriodicCoMRef> com_ref_;
267 std::shared_ptr<FrictionCone> friction_cone_;
268
269 bool addStep(const double t);
270
271 void resetContactPlacements(const double t, const Eigen::VectorXd& q,
272 const Eigen::VectorXd& v);
273
274};
275
276} // namespace robotoc
277
278#endif // ROBOTOC_MPC_PACE_HPP_
Base interface of contact planners.
Definition: contact_planner_base.hpp:21
Contact status of robot model.
Definition: contact_status.hpp:32
The state feedback and feedforward policy of LQR subproblem at a time stage.
Definition: lqr_policy.hpp:16
MPC solver for the pace gait of quadrupeds.
Definition: mpc_pace.hpp:40
~MPCPace()
Destructor.
const Eigen::VectorXd & getInitialControlInput() const
Get the initial control input.
const std::shared_ptr< ContactSequence > & getContactSequence() const
Gets the const handle of the contact sequence.
Definition: mpc_pace.hpp:231
const aligned_vector< LQRPolicy > & getLQRPolicy() const
Gets of the local LQR policies over the horizon.
std::shared_ptr< ConfigurationSpaceCost > getBaseRotationCostHandle()
Gets the base rotation cost handle.
MPCPace()
Default constructor.
MPCPace & operator=(const MPCPace &)=default
Default copy assign operator.
std::shared_ptr< Constraints > getConstraintsHandle()
Gets the constraints handle.
void setGaitPattern(const std::shared_ptr< ContactPlannerBase > &foot_step_planner, const double swing_height, const double swing_time, const double stance_time, const double swing_start_time)
Sets the gait pattern.
MPCPace(MPCPace &&) noexcept=default
Default move constructor.
MPCPace(const MPCPace &)=default
Default copy constructor.
void reset()
Resets the optimal control problem solover via the solution computed by init().
void setSolverOptions(const SolverOptions &solver_options)
Sets the solver options.
std::shared_ptr< CoMCost > getCoMCostHandle()
Gets the com cost handle.
void setRobotProperties(const RobotProperties &properties)
Sets a collection of the properties for robot model in this MPC.
std::shared_ptr< FrictionCone > getFrictionConeHandle()
Gets the friction cone constraints handle.
double KKTError() const
Returns the l2-norm of the KKT residuals. MPCPace::updateSolution() must be computed.
ControlPolicy getControlPolicy(const double t) const
Gets the control policy at the specified time.
Definition: mpc_pace.hpp:157
std::shared_ptr< CostFunction > getCostHandle()
Gets the cost function handle.
MPCPace(const Robot &quadruped_robot, const double T, const int N)
Construct MPC solver.
double KKTError(const double t, const Eigen::VectorXd &q, const Eigen::VectorXd &v)
Computes the KKT residual of the optimal control problem.
void init(const double t, const Eigen::VectorXd &q, const Eigen::VectorXd &v, const SolverOptions &solver_options)
Initializes the optimal control problem solover.
const Solution & getSolution() const
Get the solution.
void updateSolution(const double t, const double dt, const Eigen::VectorXd &q, const Eigen::VectorXd &v)
Updates the solution by iterationg the Newton-type method.
const OCPSolver & getSolver() const
Gets the const handle of the MPC solver.
Definition: mpc_pace.hpp:225
std::shared_ptr< ConfigurationSpaceCost > getConfigCostHandle()
Gets the configuration space cost handle.
std::vector< std::shared_ptr< TaskSpace3DCost > > getSwingFootCostHandle()
Gets the swing foot task space costs (LF, LH, RF, RH feet) handle.
Optimal control problem solver by Riccati recursion.
Definition: ocp_solver.hpp:41
Dynamics and kinematics model of robots. Wraps pinocchio::Model and pinocchio::Data....
Definition: robot.hpp:32
Definition: constraint_component_base.hpp:17
aligned_vector< SplitSolution > Solution
Solution to the optimal control problem.
Definition: solution.hpp:16
std::vector< T, Eigen::aligned_allocator< T > > aligned_vector
std vector with Eigen::aligned_allocator.
Definition: aligned_vector.hpp:14
Control pocily constructed for the MPC solution.
Definition: control_policy.hpp:17
Collection of the robot properties, which can change after constructing robot models.
Definition: robot_properties.hpp:30
Options of optimal control solvers.
Definition: solver_options.hpp:17