Abstract: For the development of new concepts of advanced high performance aircraft and undersea vehicles, it is great significant to investigate the unsteady movement of creatures in nature.Based on the immersed boundary method(IB) and lattice Boltzmann method (LBM), an large-scale 3D IB-LBM parallel solver was developed. In order to investigate the influence of aspect ratio on the thrust performance of bionic motion wing under low Reynolds number, the numerical simulation of the NACA0012 bionic flapping wings with four different aspect ration were carried out in Guangzhou Chinese Tianhe 2 supercomputer. The detail of the vortex structure and evolution procedure were captured, and the unsteady force of NACA0012 wing were given. The numerical results show that the two-dimensional effect of the vortex structure of the wing becomes more apparent as the increases of the aspect ratio. However, the trust lift is increased first and then decreased with the increase of the the aspect ratio. It indicates that for the bionic vehicle design the aspect ration of the bionic motion wing must be selected and evaluated carefully.