SIH: add transport rate acceleration to local acceleration

This commit is contained in:
bresch
2024-11-26 11:45:46 +01:00
committed by Mathieu Bresciani
parent 674aa474e7
commit cd18138b1c
+5 -1
View File
@@ -445,7 +445,11 @@ void Sih::equations_of_motion(const float dt)
_specific_force_E = acceleration_E - gravity_acceleration_E;
_v_E_dot = acceleration_E + coriolis_acceleration_E;
_v_N_dot = _R_N2E.transpose() * _v_E_dot; //TODO: add transport rate
// add fictitious transport rate acceleration as the local navigation frame rotates
// to stay tangent to the ellipsoid
const Vector3f transport_rate = -_lla.computeAngularRateNavFrame(_v_N).cross(_v_N);
_v_N_dot = _R_N2E.transpose() * _v_E_dot + transport_rate;
// forward Euler velocity intergation
Vector3f v_E_prev = _v_E;