asdeasoft.net

You have 2 options:

1. Use the "localForce" result. In this case, you can get results only at the 2 end nodes because this is a result on the nodes of the element. In this case, you need the PythonAPI scripting approach.
2. Use the "section.force" result. This result is on the Gauss points, i.e. on each cross section. This is easier: you can visualize it using Gauss Plot, and you can select a Gauss point and extract chart data (of one or multiple locations).
   Also here, if you want, you can make the extraction programmatic using Python.

   

   GP.png (57.85 KiB) Viewed 665 times

Approach 2 (section.force) is more versatile, but as of now you cannot use it: you modeled your beam as an Euler Bernoulli beam (non-shear-deformable), so you cannot see the shear forces in the cross-section, but only at the element ends for equilibrium.
However, You can easily turn this element into a shear deformable one using a ForceBeamColumn element.
You need to activate the Optional parameter in the Elastic section (if you are using the Elastic section) or use a section Aggregator if you are using the Fiber section