Based on the higher order shear deformation theory, the globallocal 1, 23 mode is derived in this paper. Given the conditions of the continuity of inplane deformation and shear stress, as well as the free conditions of the upper and lower surfaces, the number of the independent variables for each node of laminated plate is reduced to 13, which is independent of the number of laminate layers. The displacement mode is then developed through rigidflexible coupling dynamics for laminated plate mutibody system. Based on the principle of virtual work, the dynamic equations for mutibody system are established through the hybrid coordinate formulation, where the continuity of the interlaminar stress is taken into consideration for laminated plate. The globallocal 1,23 displacement mode is obtained and verified against the exact solution in statics example. Finally, the comparison of the results obtained by the proposed method and the traditional method for laminated plate mutibody system illustrates the necessity of considering the continuity of interlaminar shear stress in simulation of laminated plate mutibody system.