for the linearlized differential algebraic equation of the nonlinear constrained system, exact initial values of the accelerations are needed to solve itself. It may be very troublesome to perform the inverse operation for obtaining the incremental qu...
for the linearlized differential algebraic equation of the nonlinear constrained system, exact initial values of the accelerations are needed to solve itself. It may be very troublesome to perform the inverse operation for obtaining the incremental quantities since the mass matrix contains zero element in the diagonal. This fact makes the mass matrix impossible to be positive definite. To overcome this singularity phenomenon the mass matrix needs to be modified to allow the feasible application of predictor and corrector in the iterative computation. In this paper the proposed numerical algorithm based on the modified mass matrix combines the conventional implicit algorithm, Newton-Raphson method and Newmark method. The numerical example presents reliabilities for the proposed algorithm via comparisions of the 4th order Runge-Kutta method. The proposed algorithm seems to be satisfactory even though the acceleration, Lagrange multiplier, and energy show unstable behavior. Correpondingly, it provides one important clue to another algorithm for the enhancement of the numerical results.