CONSTRAINT VIOLATION STABILIZATION OF EULER-LAGRANGE EQUATIONS WITH NON-HOLONOMIC CONSTRAINTS

Two constraint violation stabilization methods are presented to solve the Euler Lagrange equations of motion of a multibody system with nonholonomic constraints. Compared to the previous works, the newly devised methods can deal with more complicated problems such as those with nonholonomic constraints or redundant constraints, and save the computation time. Finally a numerical simulation of a multibody system is conducted by using the methods given in this paper.

AN AUTOMATIC CONSTRAINT VIOLATION STABI- LIZATION METHOD FOR DIFFERENTIAL/ ALGEBRAIC EQUATIONS OF MOTION IN MULTIBODY SYSTEM DYNAMICS

A new automatic constraint violation stabilization method for numerical integration of Euler_Lagrange equations of motion in dynamics of multibody systems is presented. The parameters α,β used in the traditional constraint violation stabilization method are determined according to the integration time step size and Taylor expansion method automatically. The direct integration method, the traditional constraint violation stabilization method and the new method presented in this paper are compared finally.

Projected Runge-Kutta methods for constrained Hamiltonian systems

Projected Runge-Kutta （R-K） methods for constrained Hamiltonian systems are proposed. Dynamic equations of the systems, which are index-3 differential-algebraic equations （DAEs） in the Heisenberg form, are established under the framework of Lagrangian multipliers. R-K methods combined with the technique of projections are then used to solve the DAEs. The basic idea of projections is to eliminate the constraint violations at the position, velocity, and acceleration levels, and to preserve the total energy of constrained Hamiltonian systems by correcting variables of the position, velocity, acceleration, and energy. Numerical results confirm the validity and show the high precision of the proposed method in preserving three levels of constraints and total energy compared with results reported in the literature.

Constraining Lorentz invariance violation from the continuous spectra of short gamma-ray bursts

In some quantum gravity theories, a foamy structure of space-time may lead to Lorentz invariance violation（LIV）. As the most energetic explosions in the Universe, gamma-ray bursts（GRBs） provide an effect way to probe quantum gravity effects. In this paper, we use the continuous spectra of 20 short GRBs detected by the Swift satellite to give a conservative lower limit of quantum gravity energy scale MQG. Due to the LIV effect, photons with different energy have different velocities. This will lead to the delayed arrival of high energy photons relative to low energy ones. Based on the fact that the LIV-induced time delay cannot be longer than the duration of a GRB,we present the most conservative estimate of the quantum gravity energy scales from 20 short GRBs. The strictest constraint, M_（QG） 〉 5.05 × 10^（14） GeV in the linearly corrected case, is from GRB 140622 A. Our constraint on MQG,although not as tight as previous results, is the safest and most reliable so far.