Bolted joints play an important role in aerospace, machinery manufacturing, weapons and other fields, and the contact pressure distribution at the connection interface seriously affects the service performance of the bolts. Contact pressure analysis is an essential basis and reference for structural design, calibration, inspection, and safety monitoring of bolted connection structures. Considering the preload force and frictional contact between the joint components, a block mapping hexahedron mesh generation and finite element modeling method for bolted connections is presented, and the finite element models of single and double bolted structures are constructed. Then, an experimental test platform for measuring contact pressure distribution is constructed. Comparison of finite element analysis results with experimental results for the contact pressure distribution in the single and double-bolted joints; the root mean squared error of contact pressure distribution is less than 5%, which verifies the effectiveness of the finite element models. After that, the models are used to study the contact pressure distribution in single-bolted joints and contact pressure coupling effect of double-bolted joints, and the normal stress distribution features within members and contact pressure distribution contour are revealed for single and double-bolted joint. Finally, the finite element models are adopted to investigate the effects of the clamping length, preload, material properties, hole clearance, and bolt size on contact pressure distribution in single-bolted joints. The coupling effect of contact pressure distribution in double-bolted joints under different preloads and geometric parameters are also examined.