This research advances hybrid soft-rigid robot simulations, achieving up to 1000 times faster computations through analytical derivatives in the GVS framework.