×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Jobs

Node co-ordinates and connectivity for Vumat.
2

Node co-ordinates and connectivity for Vumat.

Node co-ordinates and connectivity for Vumat.

(OP)
Hello all,

I'l looking for some suggestions on how to get the node co-ordinate data and element connectivity for use in a Vumat routine. I need to perform some integration with this data in my routine. The keyword list gives the coordinate of the material points and characteristic length of the element but this is not sufficient for my purpose. I would like to have lets say..the height/length of the element in z-dir.
Is it possible to get this through some keyword in Vumat? If anybody has an idea pls share.

Im using 3d brick elements with reduced integration.

RE: Node co-ordinates and connectivity for Vumat.

Seems as though we're both trying to do something very similar! I've been trying to do this for a good while now..best I can come up with is reading in the element and nodal connectivity from the INP the first time the VUMAT is called and storing the info in a common block. Its not a very nice approach but can't seem to find any other way of doing it. Would be very interested in hearing if someone else has a better idea...

RE: Node co-ordinates and connectivity for Vumat.

(OP)
Yeah..that is one way of doing it but as you said not ideal. Can you share how you are doing this but? I'm sorry I know only as much Fortran as just to make my routine work. I'm learning more about it on the fly.

RE: Node co-ordinates and connectivity for Vumat.

2
Yeah no problem. I'll attach a copy of the subroutine, its not so well commented so I'll give a brief run-through of what it does. The subroutine is only called once in the whole analysis (first integration point, first increment) so expensive calls like 'vgetinternal' aren't a big deal.

1) Open input file for reading
2) Skip down to *Element in input file
3) Read in list of nodes for each element
4) Loop through each element and store the number of nodes common to each element (e.g. if two elements share four nodes then they are beside each other...based on C3D8R element).
5) Store the resulting element map in a common block.

The code returns an element map for use in the main program. This means that I know which elements are beside each other from within my VUMAT. I'm using it at the moment for non-local damage but with a small bit of modification you could use it to store nodal co-ordinates too.

Using common blocks like this are a last resort for me (aside from looking into coding a UEL). They can result in some nasty bugs and the code doesn't work so well over multiple processors. Let me know of you want more details.

BTW I'm fairly new to Fortran coding too so if anyone sees anything that can be improved on in the attached routine let me know!

RE: Node co-ordinates and connectivity for Vumat.

(OP)
Many thanks for sharing the routine MechIrl. Really appreciate it!
I was away for a while so couldn't reply earlier. I'm looking at a different formulation insead of using the node details..so waiting to see if the reults are good enough.
So it might happen that I finally don't have to use the element and node connectivity for my calculations.

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members!


Resources