Numerical simulation plays an important role for the prediction of optical trapping based on plasmonic nano-optical tweezers. However, complicated structures and drastic local field enhancement of plasmonic effects bring great challenges to traditional numerical methods. In this article, an accurate and efficient numerical simulation method based on a dual-primal finite element tearing and interconnecting (FETI-DP) and Maxwell stress tensor is proposed, to calculate the optical force and potential for trapping nanoparticles. A low-rank sparsification approach is introduced to further improve the FETI-DP simulation performance. The proposed method can decompose a large-scale and complex problem into small-scale and simple problems by using non-overlapping domain division and flexible mesh discretization, which exhibits high efficiency and parallelizability. Numerical results show the effectiveness of the proposed method for the prediction and analysis of optical trapping at nanoscale.