Thank you for reply.
I think my hand calculation considered perimeter of the hole, it because the area that has stress is based on the cross sectional tube area that is close to the hole. the final simulation is very close to my hand calculation. I thought the max values of top and bottom is the displacements take place the most at the center, relative to the constraints (since the constraints will represent the initial position of the plates) Therefore I added two max values to find out the total displacement. 0.136+0.1389 = 0.2749.(nastran 2019) but this value turned out to be different when i used nastran 2020. Nastran 2020 showed almost half of 0.2749.
Cross sectional area for tube
Ac = π * (OD2 – ID2) / 4 = π * (32.0042 – 38.12) / 4 = 335.643 mm2
Member Stiffness
K = E * Ac / L
L=Washer thickness * 2 + Flange thickness * 2 + Tube spacer height
L = 3.4*2 + 38.1*2 +38.725 = 121.725
E = 1.90*10^5
Ac= 335.643
K=523872.5 N/mm
Displacement = K/ Bolt Preload = 523872.5/144094 = 0.275057 mm

Do you think applying constraints on the eight edges better?
I attached the files.
Thank you for your help!