Notifications
Clear all

Understanding VUMAT code written in Mechanics of Solid polymers Textbook.


mzohebnawaz
Posts: 2
Topic starter
(@mzohebnawaz)
New Member
Joined: 2 months ago

Hello Everyone,

I am new to Hyperelastic materials, I am trying to learn through the Mechanics of Solid Polymers by Jorgen's. I am have done a uniaxial test on rubber. From the test data, I have generated Yeoh's coefficient through Abaqus inbuilt function. But when I tried to validate the tensile test, the results are not matching.  So I am trying to implement the Ne0-Hookean VUMAT model in the book.

Before Implementing I want to understand the code. Can anyone help with understanding the code, how it is written?

Along with that, I am unaware of implementing the code in ABAQUS. Can anyone please guide me with it?

 

Thank you

 

Regards,

 

 

5 Replies
jbRektor
Posts: 17
Admin
(@jbrektor)
Member
Joined: 1 year ago

I'm happy to help. The basic steps are as follows:

  1. Type in the code from my book into a text editor.
  2. Install a Fortran compiler that is supported by Abaqus.
  3. Run: abaqus make library=[name of fortran file]
  4. Run your abaqus test case using the new shared library file created by abaqus. Note that you may need to update your abaqus_v6.env file for abaqus to find the file.

/Jorgen

Reply
mzohebnawaz
Posts: 2
Topic starter
(@mzohebnawaz)
New Member
Joined: 2 months ago

I am using ABAQUS 6.13. I am unable to get Intel fortran 2013 version. Can you please help me with it.

Reply
jbRektor
Posts: 17
Admin
(@jbrektor)
Member
Joined: 1 year ago

In that case you may be able to use gfortran (which is free). Note that I cannot help with how to use gfortran with Abaqus, since I only use the Intel compilers.

/Jorgen

Reply
Chhorn buntheng
Posts: 1
(@chhorn-buntheng)
New Member
Joined: 2 months ago

Hi everyone,

I am simulating flexural RC beam strengthened using FRP plate. FRP plate was modeled as a solid element, so UMAT subroutine was used to assign material properties of FRP plate. After running, I got an error (300 element have been defined with zero hour glass stiffness. You may use *hourglass stiffness or change the element types. The element have been identified in element set ErrElemZeroHourGlassStiffness)

Please everyone help me to solve this problem.     

Thank you so much for your help in advance.

Sincerely

Buntheng

Reply
twhohen
Posts: 17
(@twhohen)
Active Member
Joined: 3 years ago

You can find an updated VUMAT posted on this website here:

https://polymerfem.com/community/user-material-subroutines/vumat-for-hyperelastic-materials-update-to-bergstroms-example/

I can share a paper with you that discusses VUMAT implementation in Abaqus. Search "Travis Hohenberger" in Researchgate and message me there. I'll send you the paper. Or if you have library services, you may have access. The article is in Rubber Chemistry & Technology under my name.

The oldest version I've set up subroutines on is 6.14. I used Visual Studio 2013, but I believe 2015 is compatible with it too. Are you aware that you don't necessarily need the VS or Parallel Studio versions to match Abaqus' year? You may want to check again that VS 2013 is not available. I was able to locate it in 2018, though it is a bit convoluted tracking down the older versions.

Regards,
Travis

Reply
Share: