Hello Mike,

Thank you for your response. You are right ! I checked the ds.dat file and there is:

/com,*********** Elements for Body 1 “SYS-35\Beam (Extracted Profile1)” ***********

et,1,188

keyo,1,3,2

eblock,19,solid,,1

This is interesting, because I tried also with the APDL commands under the Geometry \ SYS-35/Beam \ to force using linear shape functions along the length:

SECTYPE,1,BEAM,RECT

SECDATA,0.050,0.1,2.,2.

secoffset,cent

KEYOPT,1,1,0 !Six degrees of freedom per node, unrestrained warping (default)

KEYOPT,1,3,0 !Shape functions along the length is Linear.

With this APDL command the stiffness matrices doesn’t still match the hand calculation's stiffness matrix, but still the displacement results were equivalent (in the first post).🤔

From the file.full the dense stiffness matrix:

And the local stiffness matrix from the file.emat:

Best regards,

Pirkka Ruotsalainen