lecture 24: controlling the uncontrollable building the unbuildable
DESCRIPTION
Lecture 24: Controlling the Uncontrollable Building the Unbuildable. How do we deal with this uncontrollable system?. Can we extend the method of Zs to make it work for the bicycle?. We won’t finish the latter, and I’m not sure I know the answer although I think the answer is yes. - PowerPoint PPT PresentationTRANSCRIPT
1
Lecture 24: Controlling the UncontrollableBuilding the Unbuildable
How do we deal with this uncontrollable system?
Can we extend the method of Zs to make it work for the bicycle?
We won’t finish the latter, and I’m not sure I know the answeralthough I think the answer is yes.
2
We already know that the system is not controllable, but we don’t know what happens
It is a single input system, and, were it controllable, we’d write
€
τ 2 = −G ⋅x
and the closed reduced problem would become
€
˙ x = A − B ⋅G( ) ⋅x
We would write a characteristic polynomial incorporating the gainsand match the coefficients to the coefficients of a desired polynomial
3
Here the characteristic polynomial we obtain has no constant term
No matter what the gains, there is always a zero eigenvalue!
4
We can move the other eigenvalues using the gains;choose a simple set of eigenvalues, all equal to -1
I find g[2] through g[6] in terms of g[1], and then set g[1] to zero
(This choice fails us later on.)
5
Here is a plot of the gains vs. w0
6
We can now ask the question: what is the consequence of the zero eigenvalue?
But first, let’s get a set of eigenvalues that will be useful
the poles (as conjugate pairs)
7
The gains associated with these look not that different from the previous ones
the glitches are gains swapping roles
8
I would love to substitute this into a simulationbut we don’t have one, and are not likely to have one
I can look at the linear solutions and see what the zero eigenvalue does there.
The solution to the linear problem can be written
€
x = c iVi exp λ it( )i=1
6
∑
eigenvalues
eigenvectors
9
The constants are determined by the initial conditions
Recall that the state is
€
x =
q4
q5
q16
q17
u1
u2
⎧
⎨
⎪ ⎪ ⎪
⎩
⎪ ⎪ ⎪
⎫
⎬
⎪ ⎪ ⎪
⎭
⎪ ⎪ ⎪
=
φ1
θ1
φ3
θ3
˙ θ 1˙ ψ 2
⎧
⎨
⎪ ⎪ ⎪
⎩
⎪ ⎪ ⎪
⎫
⎬
⎪ ⎪ ⎪
⎭
⎪ ⎪ ⎪
10
The eigenvector associated with the zero eigenvalue has no 5, 6 components
An initial condition that involves u, motion, will not involve this eigenvector-eigenvalue
A general initial condition will lead to a long term nonzero result.
11
12
13
14
We can look at the Mathematica for this, and then come back here
15
The Spherical Bicycle: Extended Method of Zs
Let’s review the method of Zs as described in Lecture 19
16
The game is to isolate the us, writing the derivatives as coefficients times us
I will denote the coefficients in the momentum (Hamilton and reduced Hamilton) equations by Zs
This is going to seem a wee bit weirdBear with me. It’ll make sense eventually
We’ll have to do a lot more for the bicycle problem, but let’s start by trying to understand the process once more
17
We note that pi is linear in uj, which we can write as
€
pi = Z ijuj
We can obtain the Zs from the usual momentum expression
€
pi = M ij ˙ q j = M ik ˙ q k = M ikS jku j ⇒ Z ij = M ikS j
k
or
€
Zij = ∂pi
∂u j
The latter is often easier
18
€
Z ij = M ikS jk
It is clear from the expression
that Z does not depend on u
19
We already have the equations for the evolution of q:
€
˙ q i = S ji u j
Hamilton’s equations become
€
˙ p i = Z ij ˙ u j + ˙ Z ijuj = ∂L
∂qi + Qi
where any explicit
€
˙ q i = S ji u j
We need to replace
€
˙ Z ij =∂Z ij
∂qk ˙ q k =∂Z ij
∂qk Smk um
20
So that we have Hamilton’s equations in the form
€
Zij ˙ u j +∂Z ij
∂qk Smk umu j = ∂L
∂qi + Qi
Remember that these are not actually correct because I have not written the Lagrange multipliers
€
Z ij ˙ u j +∂Z ij
∂qk Smk umu j = ∂L
∂qi + Qi + λ kCik
is the correct form
21
We have learned to multiply by S to obtain the correct formthe reduced Hamilton’s equations
€
Zij ˙ u jSni +
∂Zij
∂qk Smk umu jSn
i = ∂L∂qi Sn
i + QiSni
We need to solve these simultaneously with the q equations
€
˙ q i = S ji u j
22
€
∂L∂qi is, formally, a terribly complicated expression
€
∂L∂qi = 1
2˙ q p
∂M pq
∂qi ˙ q q − ∂V∂qi
€
∂L∂qi = 1
2Sm
p um ∂M pq
∂qi Snqun − ∂V
∂qi
However, if we relegate the connectivity constraints to the pseudononholonomic worldmost of the derivatives of M are zero, and it is not at all bad
23
€
Z ij ˙ u jSni +
∂Z ij
∂qk Smk umu jSn
i = 12
Smp um ∂M pq
∂qi Snqun − ∂V
∂qi
⎛ ⎝ ⎜
⎞ ⎠ ⎟Sn
i + QiSni
The full formal expression for the reduced Hamilton’s equations
This is a set of nonlinear first order ordinary differential equations in uthe coefficients of which are functions of q, often quite complicated functions
This would be a snap to integrate numerically were the coefficients to be constant
The method of Zs pretends that they are
24
I have adapted the method of Zs from Kane and Levinson’s 1983 paper, cited in the text
Their work used Kane’s method, which I explore in the text, but not in class
My method of Zs is not as complete as it might be, because I allow some of the coefficients to remain as explicit functions of q
This relies on
€
∂L∂qi and Si
j actually being pretty simple functions of q
(and also the generalized forces)
This is where we shall have to extend our thought processes
25
€
ZijSni ˙ u j +
∂Zij
∂qk Smk Sn
i umu j = ∂L∂qi Sn
i + QiSni
Write Hamilton’s equations as follows
We replace the Z variables by constants
€
Z ij → Tij , ∂Z ij
∂qk Smk → Tijm
€
TijSni ˙ u j + TijmSn
i umu j = ∂L∂qi Sn
i + QiSni
This choice of substitutions aligns with Hamilton’s equations, not the reduced equations
26
Of course, they are not constants, so we have to add algebraic equationsto our system to allow us to calculate them as we go forward
€
Tij = Zij qk( ), Tijm =∂Z ij qp( )
∂qk Smk qq( )
Let me try to summarize this
27
the simple holonomic constraints
generalized coordinates
Lagrangian
nonholonomic and pseudononholonomic constraints
Hamilton’s equations
Start from square one
€
˙ q i = S ji u j
introduce the Zs
reduced Hamilton’s equations with Zs
generalized forces
28
Review the development, starting after application of the simple holonomic constraints
€
L = 12
M ij ˙ q i ˙ q j −V qk( )
€
pi = ∂L∂˙ q i
= M ij ˙ q j
€
C ji ˙ q j = 0
€
˙ q j = Skjuk
€
SkjCr
k = 0[ ]r
j
€
pi = M ijSkjuk
29
€
˙ p i = ∂L∂qi + Qi
€
Qi = ∂ ˙ W ∂˙ q i
€
∂L∂qi = 1
2∂Mmn
∂qi SrmSs
numun − ∂V∂qi
€
˙ p i = M ijSkj ˙ u k + M ij
˙ S kj ˙ u k + ˙ M ijSk
j ˙ u k
€
˙ p i = M ijSkj ˙ u k + M ij
∂Skj
∂qm Srmuruk +
∂M ij
∂qm SkjSr
muruk
30
€
M ijSkj ˙ u k + M ij
∂Skj
∂qm Srmuruk +
∂M ij
∂qm SkjSr
muruk = 12
∂Mmn
∂qi SrmSk
nuruk − ∂V∂qi + Qi
€
M ijSkj ˙ u k = −M ij
∂Skj
∂qm −∂M ij
∂qm Skj + 1
2∂Mmn
∂qi Skn
⎛ ⎝ ⎜
⎞ ⎠ ⎟Sr
muruk − ∂V∂qi + Qi
and the question arises, how many of these things do you want to consolidate?
One option is to keep them all symbolic, for which I will want symbolic
S, M, gradS, gradM
€
M ijSkjSp
i ˙ u k = −M ij∂Sk
j
∂qm −∂M ij
∂qm Skj + 1
2∂Mmn
∂qi Skn
⎛ ⎝ ⎜
⎞ ⎠ ⎟Sr
mSpi uruk − ∂V
∂qi Spi + QiSp
i
Finally the (correct) reduced Hamilton equations
31
Rather than introduce Zs at the beginning, let’s extend the idea to all the equations
How many of these things do you want to consolidate?Can we work with just M and S and their gradients?
S, M, gradS, gradM
€
M ijSkjSp
i ˙ u k = −M ij∂Sk
j
∂qm −∂M ij
∂qm Skj + 1
2∂Mmn
∂qi Skn
⎛ ⎝ ⎜
⎞ ⎠ ⎟Sr
mSpi uruk − ∂V
∂qi Spi + QiSp
i
32
But first, note that we’ll want to do this for qdot as wellbecause the qdot equations are also pretty messy
We have S and we want its symbolic equivalent
33
So we’ll have symbolic qdot equations
34
We have 22 simple odesaccompanied by 52 not very simple algebraic equations
The reduced Hamilton equations are more complicated.
35
€
M ijSkjSp
i ˙ u k = −M ij∂Sk
j
∂qm −∂M ij
∂qm Skj + 1
2∂Mmn
∂qi Skn
⎛ ⎝ ⎜
⎞ ⎠ ⎟Sr
mSpi uruk − ∂V
∂qi Spi + QiSp
i
One option is to build all the pieces of this from the four basic arrays
€
Spi , ∂Sk
j
∂qm , M ij , ∂Mmn
∂qi
36
We know S and M and can calculate the gradients directly
The symbolic version (both of them in one block)
37
We can do the same with S
And what we need to do is cobble together the equations of motion,and the accompanying algebraic equations
and I think I’ll look at this on the blackboard
38
€
M ijSkjSp
i ˙ u k = −M ij∂Sk
j
∂qm −∂M ij
∂qm Skj + 1
2∂Mmn
∂qi Skn
⎛ ⎝ ⎜
⎞ ⎠ ⎟Sr
mSpi uruk − ∂V
∂qi Spi + QiSp
i
Here’s the equation we are building