xcur in fld3d1.f

[Shu]

Dear All,

I am going through the large deformation code in FEAP. I have a question when I am looking at file fld3d1.f. When it is computing Cauchy stresses and spatial tangent tensor at t_n+1, there is a statement xcur(i)=xr(i,l)+ur(i,l). However, above this line, when compute coordinates, ur(i,l) = ur(i,l) + xr(i,l). I am wondering why xcur include xr twice here. Thanks for your attention.

Shu

[FEAP_Admin]

I think this a bit of incomplete (incorrect) code. This value of xcur( ) does not appear to be used within fld3d1.f or any of the routines that it calls.  At most it is used in the strain computation routines but it is re-computed (correctly) there.

Thanks for pointing out the errant line; after double checking again we will remove it for the next release of the code.

[JStorm]

I also observed this bug.

Following the manual, xcur can be used in user materials. That is why also FEAP standard elements need to provide that variable.
Seems to be that in fld3d2.f, fld3d1u, strnen.f, fld3d2u, strn3d.f and strn3m.f that bug is already fixed, but not in fld3d1.f, mech3f.f