Learning Forum

[Giannis123]

Hello,

I'm trying to write a subroutine for a cohesive UMAT in dyn21b (FORTRAN 77)

Does anyone know how to get the nodal coordinates for each element in the global system?

I found some things online but it doesn't seem to work...

Thanks

Giannis

[Andreas Koutras]

You first need to get the element node numbers. These are available in the common block aux33loc. You can find an example in dyn21umat.f. The ix1 - ix8 are the internal node numbers.
C_TASKCOMMON (aux33loc)
common/aux33loc/
&ix1(nlq),ix2(nlq),ix3(nlq),ix4(nlq) &ix5(nlq),ix6(nlq),ix7(nlq),ix8(nlq)
The nodal coordinates can be accessed with the pointer dm_x.
If you address the dm_x array as:
real x(3,*)
then the coordinates of each node are in x(1:3,ix1) etc.
I hope this helps.

[Giannis123]

Thanks, it worked!