The dynamic balance equation of the system is set up by means of the Hamilton variational principles with the constraint conditions of the clamp points of spacers in the subconductors by the penalty function. The numerical simulation for galloping of twin bundle conductors is achieved with the time integral by Newmark method and the NewtonRaphson iteration solution for nonlinear equation. The validity is approved with simulation for the galloping of horizontal and vertical icing twin bundle conductor. The difference of aerodynamic loads, which is due to the impacts of wake flow of windward subconductor on lee subconductor, on the individual subconductors can be considered with the proposed approach. It is appropriate for further study in the field of bundle conductor galloping.