There are multiple approaches that probably vary widely depending on your preference and objective. Perhaps consider the following.
Put a circle approximately the size of a large washer (or slightly larger) on the interface surface. Put an anchor node at the circle centerpoint so you'll get a node exactly at the bolt centerline. "Shell-coat" the circle with shell elements, assigning a very high modulus of elasticity to these shell elements to simulate the much higher stiffness in the bolt grip region. Mesh the solid elements, which will follow the shell-coat mesh on the interface; or preferably extrude shell elements and a copy of the shell-coat elements, if possible, to create better solids. If bolting to ground, put a constraint (clamp) on the node at the circle centerpoint; release constraint torsional degree of freedom (dof). Read constraint reactions to get bolt forces. If bolting to another part, create a matching shell-coat circle on other part and run a very stiff, very short beam between the two centerpoint nodes, releasing beam torsional dof (and omit the above clamp).
Toggle modulus of elasticity of shell-coat elements (and beam element, if any) to a very high value (but not so high as to cause singularity and solution failure). When shell-coat elements (and beam, if any) are stiff enough, you'll see (in displacement display) near parallelism between bolted surfaces, locally in shell-coat circle, thus roughly simulating (approximating) the average heel-toe prying resistance of the clamped plates (regardless of direction of prying moment). If interface plates are nonparallel at shell-coat circle in displacement display, it means your shell coat (and/or beam, if any) is not yet stiff enough.
For any decent bolt pattern footprint (i.e., more than one line of bolts resisting moment in that direction), prying moment on individual bolts is often negligible and bolt loads can be read directly, ignoring moment reactions or beam moments. Good luck.