In this paper a highly scalable parallel formulation of the primal-dual technique is presented for index-3 constrained flexible multi-body dynamics system. The key features of the primal-dual approach are constraint preservation, preserving the original order of accuracy of time integration operators that are employed, and faster convergence rates of nonlinear iterations for the solution of flexible multi-body dynamical systems. In addition, this technique not only preserves the underlying properties of time integration operators for ordinary differential equations, but also eliminates the need for index reduction, constraint stabilization and regularization approaches. The key features of the parallel formulation of rigid and flexible modeling and simulation technology are capabilities such as adaptive high/low fidelity modeling that is useful from the initial design concept stage to the intermediate and to the final design stages in a single seamless simulation environment. The examples considered illustrate the capabilities and scalability of the proposed high performance computing (HPC) approach for large-scale simulations
