1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
|
/*
* Copyright (c) 2011-2021, The DART development contributors
* All rights reserved.
*
* The list of contributors can be found at:
* https://github.com/dartsim/dart/blob/master/LICENSE
*
* This file is provided under the following "BSD-style" License:
* Redistribution and use in source and binary forms, with or
* without modification, are permitted provided that the following
* conditions are met:
* * Redistributions of source code must retain the above copyright
* notice, this list of conditions and the following disclaimer.
* * Redistributions in binary form must reproduce the above
* copyright notice, this list of conditions and the following
* disclaimer in the documentation and/or other materials provided
* with the distribution.
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND
* CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES,
* INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF
* MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
* DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR
* CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
* LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF
* USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED
* AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
* POSSIBILITY OF SUCH DAMAGE.
*/
#include "Controller.hpp"
using namespace dart;
using namespace dynamics;
using namespace math;
Controller::Controller(
const dynamics::SkeletonPtr& _skel,
constraint::ConstraintSolver* _collisionSolver,
double _t)
{
mSkel = _skel;
mCollisionHandle = _collisionSolver;
mTimestep = _t;
mFrame = 0;
int nDof = mSkel->getNumDofs();
mKp = Eigen::MatrixXd::Identity(nDof, nDof);
mKd = Eigen::MatrixXd::Identity(nDof, nDof);
mConstrForces = Eigen::VectorXd::Zero(nDof);
mTorques.resize(nDof);
mDesiredDofs.resize(nDof);
for (int i = 0; i < nDof; i++)
{
mTorques[i] = 0.0;
mDesiredDofs[i] = mSkel->getPosition(i);
}
// using SPD results in simple Kp coefficients
for (int i = 0; i < 6; i++)
{
mKp(i, i) = 0.0;
mKd(i, i) = 0.0;
}
for (int i = 6; i < 22; i++)
mKp(i, i) = 200.0; // lower body + lower back
for (int i = 22; i < nDof; i++)
mKp(i, i) = 20.0;
for (int i = 6; i < 22; i++)
mKd(i, i) = 100.0;
for (int i = 22; i < nDof; i++)
mKd(i, i) = 10.0;
mPreOffset = 0.0;
}
void Controller::computeTorques(
const Eigen::VectorXd& _dof, const Eigen::VectorXd& _dofVel)
{
// SPD tracking
// std::size_t nDof = mSkel->getNumDofs();
Eigen::MatrixXd invM = (mSkel->getMassMatrix() + mKd * mTimestep).inverse();
Eigen::VectorXd p = -mKp * (_dof + _dofVel * mTimestep - mDesiredDofs);
Eigen::VectorXd d = -mKd * _dofVel;
Eigen::VectorXd qddot
= invM * (-mSkel->getCoriolisAndGravityForces() + p + d + mConstrForces);
mTorques = p + d - mKd * qddot * mTimestep;
// ankle strategy for sagital plane
Eigen::Vector3d com = mSkel->getCOM();
Eigen::Vector3d cop = mSkel->getBodyNode("h_heel_left")->getTransform()
* Eigen::Vector3d(0.05, 0, 0);
Eigen::Vector2d diff(com[0] - cop[0], com[2] - cop[2]);
if (diff[0] < 0.1)
{
double offset = com[0] - cop[0];
double k1 = 20.0;
double k2 = 10.0;
double kd = 100.0;
mTorques[17] += -k1 * offset + kd * (mPreOffset - offset);
mTorques[25] += -k2 * offset + kd * (mPreOffset - offset);
mTorques[19] += -k1 * offset + kd * (mPreOffset - offset);
mTorques[26] += -k2 * offset + kd * (mPreOffset - offset);
mPreOffset = offset;
}
// Just to make sure no illegal torque is used
for (int i = 0; i < 6; i++)
mTorques[i] = 0.0;
mFrame++;
}
|
