Following the previous post, I found some errors in the Python code:

• Point loads applied at the left hand end generated incorrect results.
• The effect of beam deflections on shear forces was not included.

These problems have now been fixed, and updated code and a revised spreadsheet have been added to the download file at:

py_GeomK.zip

The screenshots below show results with end point loads, compared with results from Strand7 and the 3DFrame spreadsheet: