as title.
First, make sure those values are saved to the result file. Then either use result.normal_contact_force() operator, or retrieve the values through SMISC.
There is another option that is actually the default of Mechanical force probes for contacts. You can retrieve the element nodal forces (ENFO) for the nodes under the contact elements. Assuming the contact is the only entity transferring force, the ENFO will equal the contact forces.
There is an operator for element nodal forces and these must be written to the results file. You then sum the forces in XYZ to get the total force through the contact in each direction.
one benefit of this method is it doesn’t rely on the contact element output as it uses the underlying elements for the force data.