861 |
|
\varphi _{b_m h}^2 \circ \varphi _{a_m h}^1 \circ \varphi _{b_{m - |
862 |
|
1} h}^2 \circ \ldots \circ \varphi _{a_1 h}^1 . |
863 |
|
\end{equation} |
864 |
< |
Careful choice of coefficient $a_1 \ldot b_m$ will lead to higher |
864 |
> |
Careful choice of coefficient $a_1 \ldots b_m$ will lead to higher |
865 |
|
order method. Yoshida proposed an elegant way to compose higher |
866 |
|
order methods based on symmetric splitting\cite{Yoshida1990}. Given |
867 |
|
a symmetric second order base method $ \varphi _h^{(2)} $, a |
1256 |
|
Using Equation (\ref{introEquation:motionHamiltonianCoordinate}, |
1257 |
|
\ref{introEquation:motionHamiltonianMomentum}), one can write down |
1258 |
|
the equations of motion, |
1259 |
< |
\[ |
1260 |
< |
\begin{array}{c} |
1261 |
< |
\frac{{dq}}{{dt}} = \frac{p}{m} \label{introEquation:RBMotionPosition}\\ |
1262 |
< |
\frac{{dp}}{{dt}} = - \nabla _q V(q,Q) \label{introEquation:RBMotionMomentum}\\ |
1263 |
< |
\frac{{dQ}}{{dt}} = PJ^{ - 1} \label{introEquation:RBMotionRotation}\\ |
1264 |
< |
\frac{{dP}}{{dt}} = - \nabla _Q V(q,Q) - 2Q\Lambda . \label{introEquation:RBMotionP}\\ |
1265 |
< |
\end{array} |
1266 |
< |
\] |
1259 |
> |
|
1260 |
> |
\begin{eqnarray} |
1261 |
> |
\frac{{dq}}{{dt}} & = & \frac{p}{m} \label{introEquation:RBMotionPosition}\\ |
1262 |
> |
\frac{{dp}}{{dt}} & = & - \nabla _q V(q,Q) \label{introEquation:RBMotionMomentum}\\ |
1263 |
> |
\frac{{dQ}}{{dt}} & = & PJ^{ - 1} \label{introEquation:RBMotionRotation}\\ |
1264 |
> |
\frac{{dP}}{{dt}} & = & - \nabla _Q V(q,Q) - 2Q\Lambda . \label{introEquation:RBMotionP} |
1265 |
> |
\end{eqnarray} |
1266 |
|
|
1267 |
|
In general, there are two ways to satisfy the holonomic constraints. |
1268 |
|
We can use constraint force provided by lagrange multiplier on the |
1343 |
|
\[ |
1344 |
|
\hat vu = v \times u |
1345 |
|
\] |
1347 |
– |
|
1346 |
|
Using \ref{introEqaution:RBMotionPI}, one can construct a skew |
1347 |
|
matrix, |
1348 |
|
\begin{equation} |
1349 |
< |
(\mathop \Pi \limits^ \bullet - \mathop \Pi \limits^ \bullet ^T |
1349 |
> |
(\mathop \Pi \limits^ \bullet - \mathop \Pi \limits^ {\bullet ^T} |
1350 |
|
){\rm{ }} = {\rm{ }}(\Pi - \Pi ^T ){\rm{ }}(J^{ - 1} \Pi + \Pi J^{ |
1351 |
|
- 1} ) + \sum\limits_i {[Q^T F_i (r,Q)X_i^T - X_i F_i (r,Q)^T Q]} - |
1352 |
|
(\Lambda - \Lambda ^T ) . \label{introEquation:skewMatrixPI} |
2011 |
|
\begin{array}{l} |
2012 |
|
\Xi _P^{tt} = \Xi _O^{tt} \\ |
2013 |
|
\Xi _P^{tr} = \Xi _P^{rt} = \Xi _O^{tr} - U_{OP} \Xi _O^{tt} \\ |
2014 |
< |
\Xi _P^{rr} = \Xi _O^{rr} - U_{OP} \Xi _O^{tt} U_{OP} + \Xi _O^{tr} U_{OP} - U_{OP} \Xi _O^{tr} ^{^T } \\ |
2014 |
> |
\Xi _P^{rr} = \Xi _O^{rr} - U_{OP} \Xi _O^{tt} U_{OP} + \Xi _O^{tr} U_{OP} - U_{OP} \Xi _O^{{tr} ^{^T }} \\ |
2015 |
|
\end{array} |
2016 |
|
\label{introEquation:resistanceTensorTransformation} |
2017 |
|
\end{equation} |