UMAT - coordinates of Gauss integration point

[mubeen]

Hi all,
is it possible in FEAP to get the Physical/Global coordinates of the Gauss integration point upon which the UMAT subroutine is being run ?

Reason:
The non-local material model includes the gradients of some scalar internal variables, and I don't intend to treat the internal variables as indepent variables (e.g. displacements) in a User element routine, as it increases the size of stiffness matrix (K) of the element 5-folds. So the other possibility is to inform the UMAT about the coordinates of  integration point (and keep those coordinates in a history array) upon which it is returning the material behavior, so that it can calculate the gradients on that intergation point.

[FEAP_Admin]

Gauss point coordinates are not available in umatlXX.f

However if you want, you can make the changes relatively easily.  The essential idea would be to add the array to the list of arguments.  The changes would have to be propagated up the call chain all the way up to the elements.   So for example make the changes in:

user/umatl1.f,
user/umodel.f,
elements/material/small/modlsd.f,
elements/solid2d/sdl2d1.f  [in this routine the coordinates will be in xx and yy]