14.02.2013 Views

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

Thesis - Leigh Moody.pdf - Bad Request - Cranfield University

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Appendix I / Utilities / Quaternions<br />

_ _<br />

B<br />

A<br />

B<br />

B<br />

B B<br />

( t + ∆t<br />

) : = ( Q&<br />

⋅ ∆t<br />

) ⊗ Q : = ∆Q<br />

Q<br />

Q ⊗<br />

A<br />

Equation 22.10-20<br />

Over time period (∆t) the rotating frame moves through an angle (ωAB.∆t)<br />

about a vector defined by the direction cosines (ωA,B/ωA,B). If the angular<br />

rate remains constant over this period the quaternion representing the<br />

rotation,<br />

∆Q<br />

B<br />

A<br />

∆Q<br />

B<br />

A<br />

: =<br />

: =<br />

⎛<br />

⎜ ⎛ ∆t<br />

cos ⎜ ⋅ ω<br />

⎜<br />

⎝ ⎝ 2<br />

⎛<br />

⎜<br />

⎜<br />

⎜<br />

⎜<br />

⎜<br />

⎜<br />

⎝<br />

ω<br />

ω<br />

B<br />

A,<br />

B<br />

B<br />

A,<br />

B<br />

⋅<br />

⎛<br />

⎜<br />

⎝<br />

B<br />

A,<br />

B<br />

1<br />

1 −<br />

8<br />

1<br />

2<br />

−<br />

22.10-4<br />

⎞<br />

⎟<br />

⎠<br />

⋅<br />

,<br />

ω<br />

ω<br />

A<br />

B<br />

A,<br />

B<br />

B<br />

A,<br />

B<br />

⋅ sin<br />

B ( ω ⋅ ∆t<br />

)<br />

1<br />

48<br />

⋅<br />

A,<br />

B<br />

⎛<br />

⎜<br />

⎝<br />

+ K<br />

A<br />

∆t<br />

⋅ ω<br />

2<br />

A<br />

B<br />

A,<br />

B<br />

⎞<br />

⎟<br />

⎠<br />

⎞<br />

⎟<br />

⎟<br />

⎠<br />

Equation 22.10-21<br />

( ) ⎟ ⎟⎟⎟⎟⎟<br />

B 2 ⎞<br />

ω ⋅ ∆t<br />

+ K ⎟ ⋅ ∆t<br />

Equation 22.10-22<br />

The general expression for propagating a quaternion over time interval (∆t),<br />

B ( Q )<br />

A<br />

k + 1<br />

: = Q<br />

B<br />

A<br />

+<br />

N<br />

∑<br />

n : = 1<br />

⎛<br />

⎜ ⎡<br />

0<br />

n<br />

⎜ ∆t<br />

⎢<br />

⎜ ⋅ ⎢ L<br />

⎜<br />

n + 1 ⎢ B<br />

⎜ ⎢<br />

ωA,<br />

⎝ ⎣<br />

A,<br />

B<br />

B<br />

M<br />

M<br />

M<br />

2<br />

−<br />

−<br />

⎠<br />

B ( ω )<br />

A,<br />

B<br />

L<br />

B [ ω × ]<br />

A,<br />

B<br />

T<br />

n<br />

⎞<br />

⎠<br />

⎤<br />

⎥<br />

⎥<br />

⎥<br />

⎥⎦<br />

⋅ Q<br />

B<br />

A<br />

⎞<br />

⎟<br />

⎟<br />

⎟<br />

⎟<br />

⎠<br />

Equation 22.10-23<br />

Q_RATE takes a quaternion representing the transformation from frame (A)<br />

to frame (B), and the angular velocity vector of (B) with respect to (A). It<br />

returns the quaternion time rate of change representing the rate of rotation of<br />

(B) with respect to (A), using 1 st order of approximation.<br />

B<br />

A<br />

Q_RATE<br />

B B B<br />

( Q , ω , Q&<br />

)<br />

A<br />

A,<br />

B<br />

A<br />

Equation 22.10-24<br />

B B B<br />

B B B<br />

( − ω • q , q ( 0 ) ⋅ ω − ω q )<br />

2 ⋅ Q : =<br />

×<br />

&<br />

A,<br />

B<br />

A<br />

A<br />

A,<br />

B<br />

A,<br />

B<br />

A<br />

Equation 22.10-25

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!