For robot-assisted pelvic fracture reduction, at least two bone needles need to be inserted into the ilium of the affected pelvis, and the robot clamping device is connected with the bone needles. The biomechanical properties of the pelvic musculoskeletal tissues are different with the different Spatial Position and Orientation (SPO) of the bone needles. In order to determine the optimal SPO of bone needle pairs, the constraints between the bone needles and the pelvis are analyzed, and the SPO vectors of 150 groups bone needles are obtained by the KNN-hierarchical clustering method; a batch modeling method of bone needles with different SPO is proposed. 150 finite element models of damaged pelvic musculoskeletal tissue with different SPO of bone needles are established and simulated. The stress and strain distribution homogenization of musculoskeletal tissue with bone needles as evaluation index, the simulation results of 150 models are evaluated. Results show that, the anterior superior iliac spine and the anterior inferior iliac spine are suitable regions to place bone needles in the pelvis, and the optimal distribution of the needle combination is found in this region. The overall stress and strain distribution of the damaged pelvic musculoskeletal tissue under the large reduction force is the best.Copyright © 2023. Published by Elsevier Ltd.