Bugs in an element routine

[nguyenvinhhiep]

Dear FEAP Admin,

I am not sure but maybe there are some bugs in subroutine fld3d1.
In the coding line 113, ur(i,l) = ur(i,l) + xr(i,l) should leads to current coodinates of element which are assigned again to ur.
And In the line 150, xcur(i) = xr(i,l) + ur(i,l) doest not make sense due to the combination between reference and current coordinates.
Maybe I have confused somewhere or bugs really exist.

Thank you.

Best regards,
Hiep Nguyen

[Prof. R.L. Taylor]

I do not think that xref or xcur are used anymore for the 3-d element.  The update of ur is to give the current radius for computing angular velocity effects  r*omega**2.

Some clean up is clearly needed!

[Prof. R.L. Taylor]

I seem to recall now that these were added to allow users to have access to coordinates at the gauss point in case they were needed in their constitutive equations.  Currently, only some of the values are ever used.  So this is a bug and needs to be corrected to give correct values.  Thank you for pointing out the error,  It probably exists elsewhere too.