Help on Distribution of Bolt Load in 3D

I am kind of new in computational simulation except for CFD. I have a cylinder which has a flange around it and a lid on top. Lid and the flange are connected by bolts and sealed with an o ring. Is it possible to simulate the load distribution around the flange due to tightening the bolts? Which package is most suitable for this comsol or Abaqus? I want to do simple simulation just to check whether the odd number of bolts are able to distribute load uniformly around the flange.

Thanks in advance.