The solution of the dynamic equations of the six-axis accelerometer is a prerequisite for sensor calibration, structural optimization, and practical application. However, the forward dynamic equations (FDEs) and inverse dynamic equations (IDEs) of this type of system have not been completely solved due to the strongly nonlinear coupling relationship between the inputs and outputs. This article presents a comprehensive study of the FDEs and IDEs of the six-axis accelerometer based on a parallel mechanism. Firstly, two sets of dynamic equations of the sensor are constructed based on the Newton–Euler method in the configuration space. Secondly, based on the analytical solution of the sensor branch chain length, the coordination equation between the output signals of the branch chain is constructed. The FDEs of the sensor are established by combining the coordination equations and two sets of dynamic equations. Furthermore, by introducing generalized momentum and Hamiltonian function and using Legendre transformation, the vibration differential equations (VDEs) of the sensor are derived. The VDEs and Newton–Euler equations constitute the IDEs of the system. Finally, the explicit recursive algorithm for solving the quaternion in the equation is given in the phase space. Then the IDEs are solved by substituting the quaternion into the dynamic equations in the configuration space. The predicted numerical results of the established FDEs and IDEs are verified by comparing with virtual and actual experimental data. The actual experiment shows that the relative errors of the FDEs and the IDEs constructed in this article are 2.21% and 7.65%, respectively. This research provides a new strategy for further improving the practicability of the six-axis accelerometer.