Conventional analytical and numerical methods for the mechanical properties of helical threads are relied on many assumptions and approximations and thus hardly yield satisfied results. A parameterized 3D finite element model of bolted joints with real helical thread geometry is established and meshed with refined hexahedral elements. The Von Mises plasticity criterion, kinematic hardening rule of materials and interfacial contacts are employed to make it possible for the suggested model be able to approach real assembly conditions. Then, the mechanical properties of bolted joints with different thread pitches, thread numbers and modular ratios are investigated, including the contact pressure distribution at joint interfaces, the axial load distribution and stress concentration in screw threads during the loading and unloading process. Simulation results indicate that the load distribution in screw threads produced by the suggested model agrees well the results from CHEN’s photoelastic tests. In addition, an interesting phenomenon is found that tightening the bolt with a large preload first and then adjusting the clamping force by unloading can make the load distribution more uniform and reduce the maximum residual equivalent stress in thread roots by up to 40%. This research provides a simple and practical approach to constructing the 3D finite element model and predicting the mechanical properties of helical thread connection. Conventional analytical and numerical methods for the mechanical properties of helical threads are relied on many assumptions and approximations and hereby hardly yield satisfied results. A parameterized 3D finite element model of bolted joints with real helical thread geometry is established and meshed with refined hexahedral elements. The Von Mises plasticity criterion, kinematic hardening rule of materials and interfacial contacts are employed to make it possible for the suggested model be able to approach real assembly conditions. Then, the mechanical properties of bolted joints with different thread pitches, thread numbers and modular ratios are investigated, including the contact pressure distribution at screw threads during the loading and unloading process. The simulation results indicate that the load distribution in screw threads produced by the suggested model agrees well the results from CHEN’s photoelastic tests. In addition, an interesting phenomenon is found that tightening the bolt with a large preload first and then adjusting the clamping force by unloading can make the load distribution more uniform and reduce the maximum residual equivalent stress in the thread roots by up to 40%. This research provides a simple and practical approach to constructing the 3D finite element model and predicting the mechanical properties of helical thread connection.