Derivado de matriz de transformación para dinámica 6DoF

Estoy implementando un planificador para un robot submarino 6DOF y estoy usando la dinámica derivada en el capítulo 7.5 del Handbook of Marine Craft de Fossen. Estoy usando las ecuaciones de movimiento expresadas en NED usando posiciones y ángulos de Euler para usar el control de planitud diferencial. Véase la ecuación 7.190.

Parte de esto es la transformación de velocidades (lineales y angulares: [x_dot, y_dot, z_dot, p, q, r]) entre el marco del mundo (NED) y el marco del cuerpo del robot. Esto se describe en la ecuación 7.191 utilizando la matriz j , que transforma las velocidades lineales y angulares entre el marco del mundo fijo y el marco del cuerpo:

j Θ ( η ) = [ R b norte ( Θ norte b ) 0 3 × 3 0 3 × 3 T Θ ( Θ norte b ) ]
dónde:
η ˙ = j Θ ( η ) v
T Θ ( Θ norte b ) = [ 1 pecado ϕ broncearse θ porque ϕ broncearse θ 0 porque ϕ pecado ϕ 0 pecado ϕ / porque θ porque ϕ / porque θ ]
η es la posición/orientación en el marco del mundo fijo: [ X , y , z , ϕ , θ , ψ ]

η ˙ son las velocidades en el marco mundial: [ X ˙ , y ˙ , z ˙ , pag , q , r ]

v son las velocidades en la estructura del cuerpo.

Mi problema es que para encontrar la aceleración en el marco del mundo, necesito saber j ˙ , que parece que no puedo encontrar una definición en el libro de texto de Fossen. Véase la ecuación 7.192: C y η ¨ ambos dependen de j ˙ . Soy consciente de la derivada temporal de la matriz de rotación. R usando la matriz simétrica sesgada, pero no estoy seguro de cómo encontrar la derivada del todo j matriz. ¿Alguien sabe qué debo hacer o dónde buscar más información?

Ejemplo de uso de j ˙ :

η ¨ = j Θ ( η ) v ˙ + j ˙ Θ ( η ) v

¿por qué pones esta pregunta?

Respuestas (2)

En el ejemplo, creo que J podría resolverse usando un enfoque de ecuación diferencial lineal. ¿Quizás un método iterativo como el de Euler o el de Runge-Kuta?

Esto, por supuesto, supone que v y su derivada temporal no dependen de n(eta) [perdón por formatear, estoy en el móvil].

Una vez que se haya resuelto J, me imagino que las condiciones de contorno correctas lo llevarán a una respuesta útil.

Vale la pena intentarlo al menos.

empiezas con:

(1) R = F ( q )

dónde q son las coordenadas generalizadas

derivada del tiempo:

R ˙ = d d t ( F ( q ) ) = F q q ˙ = j q ˙

de este modo

R ˙ = R ˙ ( q , q ˙ ) = gramo ( q , q ˙ )

(2) R ¨ = gramo q ˙ j q ¨ + gramo q q ˙ = j q ¨ + j q ˙ q j ˙ q ˙

Observaciones: las ecuaciones de movimiento

Ecuaciones de NEWTON:

METRO R ¨ = F A

con ecuación (2)

(A) METRO ( j q ¨ + j ˙ q ˙ ) = F A + F C

dónde F A son las fuerzas aplicadas y F C son las fuerzas de restricción.

para eliminar las fuerzas de restricción, multiplique la ecuación (A) desde la izquierda con j T

j T METRO j q ¨ = j T F A j T METRO j ˙ q ˙