In this example, we will use Tutorial Lesson 3 (Tied back excavation) [link] for the Plaxis file. In this case, we want to determine the force in the top anchor that is connected on the left retaining wall. This anchor is activated in Phase_3, and is connected with the wall at X = 40 m and Y = 27 m.

In order to retrieve this anchor force, we will use Output’s Remote Scripting environment with the Python wrapper:

• in Output program open SciTE program under Expert>Python>Editor
  note that this will launch the Configure remote scripting server window.
• add the boilerplate code at the beginning of the script [link to detailed article]
• and we should make sure the Plaxis results are opened in the Output program
• then we will obtain the anchor force for each phase in which the anchor is active. Since the ground anchor is activated in Phase_3, we can use all phases starting from this phase ( g_out.Phases[3:] )
• for each phase we will obtain the anchor results: coordinates X and Y, and the anchor force N.
• from which we will filter out the node that belongs to the top anchor of the left retaining wall (at X= 40.0, Y = 27 m). We can identify an anchor by one of its endpoints. Because the order of the nodes is the same in each of these getresults-lists, we can use the X and Y coordinates to filter out the correct anchor force.

# geometric data for correct anchor: 
Xanchor = 40.0 
Yanchor = 27.0  

# start from Phase_3 since the anchors are
# not active before this phase 
for phase in g_o.Phases[3:]: 
    anchorF = g_o.getresults(phase, 
                             g_o.ResultTypes.NodeToNodeAnchor.AnchorForce2D, 
                             'node') 
    anchorX = g_o.getresults(phase, 
                              g_o.ResultTypes.NodeToNodeAnchor.X, 
                             'node') 
    anchorY = g_o.getresults(phase, 
                             g_o.ResultTypes.NodeToNodeAnchor.Y, 
                             'node') 
    
    for x,y,N in zip(anchorX,anchorY,anchorF): 
        # check for the correct anchor by location: 
        if abs(x - Xanchor) < 1E-5 and abs(y - Yanchor) < 1E-5: 
            print("Anchor force N in {}: {:.2f} kN/m".format(phase.Name,N))

The result in the Python console would look similar to this:

  Response in Python console
>>>  
Anchor force N in Phase_3: 500.00 kN/m 
Anchor force N in Phase_4: 544.18 kN/m 
Anchor force N in Phase_5: 468.12 kN/m 
Anchor force N in Phase_6: 460.47 kN/m 
>>>