Accurate replication of the torque-tightening effect is important when simulating the mechanical behavior of composite joints using the finite element (FE) approach. In this study, a parametric mesh construction strategy was developed to automatically build a high-fidelity 3D FE model with lift angle to simulate bolt thread interactions. An experimental approach was proposed to calibrate the FE model. The relation between tightening torque and the axial force was studied by the FE method, experiment, and engineering algorithm. Results showed that traditional empirical formula overestimate the clamping force. By using the proposed model to calculate the axial force of the bolt, the accuracy of the axial force can be significantly improved when analyzing the composite bolted joints. Finally, parametric simulations were performed to investigate the influence of coefficients of friction at thread interfaces and the nut/washer interfaces. The pre-tightening force was significantly influenced by the frictional coefficients at these interfaces. Results indicated that increasing the coefficient of friction (CoF) at the nut/washer interface from 0.1 to 0.2 can lead up to 28% reduction in the clamping force, while this increase was 23% for the threads interfaces.