| |
Abstract:
Based
on the equidistant nodes and non-equidistant nodes such
as Chebyshev nodes and Legendre nodes, two methods of
block stabilization method were established in time
interval for the regular differential algebraic equation
of the Hamiltonian system in multi-body system dynamics.
The nonlinear algebraic equation of each node in a
single step interval was obtained. Newton iteration was
used in the solution process. The method is L stable.
Taking the crank slider system as an example, the
numerical simulation was carried out by using the method
in this paper. By changing the step size, simulation
time and the number of nodes, the regular differential
algebraic equation of the index-1,2,3 was verified. The
numerical results showed that our method has good
stability, high accuracy and high efficiency, and could
keep the Hamilton energy error and constraint error.
Especially for long-time simulation, Hamilton energy
error will not drift, so it is suitable for long-time
multi-body system dynamics simulation. |
|
|
