User login

You are here

Sharing ABAQUS UMAT and VUMAT subroutines

azadpoor's picture

I have noticed that a lot of people are using UMAT and VUMAT material model
subroutines for their researches. However, most of these subroutines do not
find their way to free code market and remain in the hands of the developers.
Therefore, a lot of wheel re-inventing has to be done by other people. Since
development and testing of a material subroutine takes a damn long time, I
think it is a good idea if each of us can share the subroutines she has
developed with the others. I know that a couple of subroutines are publicly
available and many other people also know that. Please comment about less-known
freely available subroutines.

To be more specific, I think it would be useful if people can comment about
following highly demanded UMAT subroutines:

1.Cohesive zone elements models

2.Polycrystalline plasticity

3.Gurson-type ductile fracture models

 Note added by admin:  This topic has also been discussed in another thread:  node/2375


Yanfei Gao's picture


I agree with you that there are many UMAT, UEL, and other user subroutines hanging around the research community. In the past several years, I have been advocating my subroutines and provide them free to the public, as can be found from 

 It'll be worthwhile to compile a list of these for a broader public. And perhaps we can get ABAQUS to endorse us. As for your three questions:

(1) cohesive zone element: I have one in my webpage, which can be easily genearlized into other versions such as Camacho-Ortiz.

(2) Polycrystalline palsticity. This I don't know. I only know the famous code by Yonggang Huang (modified by Jeff Kysar) is available at Jeff Kysar's webpage (Columbia University).

(3) Gurson-type model: I don't know if someone provides this free. The infinestimal strain version is not difficult to implement.  

Mr. Gao

I am a new comer of imechanic,  I got your webpage though the search engine Google, but I could not enter it, I dont know why.

I knew you have a user subrotine about the cohesive zone and you are willing to share it with others, resently I am studing the CZM and what to simulate the surface problems with the FEM, but I am stumpt by the subrotine. I hope you can help me and send the subrotine to my email( Thank you very much.

Excuse for my poor english, I am a university student from HEU,China. Expect your reply.

Hi Gao

I tried visiting your page however Ierror "404 not found". Can you please suggest how can I go to your page or suggest other sources I can use to get a user subroutine. I am looking for a subroutine which can perform 3D composite damage in tension or compression. Thank you





Dear Yanfei 


the link you have redirected us , does not open , shall you please check and correct ? thanks a lot

I agree. Some 15 years ago, after I had written my first UEL, I was wondering if ABAQUS would be interested in establishing a clearinghouse or repository for various user-defined subroutines. This would be similar to accessible libraries - such as NIST - for mathematic/scientific application software such as equation solvers, eigenvalue extraction routines, etc. Over time I have written UMAT and UEL routines for progressive failure analysis, higher-order plate finite elements, specialized adhesive elements, and a score of decohesion finite element incorporating various coupled cohesive zone models to joint 1-D beam, and 2-D and 3-D continuum elements. 


These subroutines are just simply sitting in folders on my PC collecting electronic 'dust'. It would be nice to upload them to a site with a short README or document file for others to use.



Zhigang Suo's picture

Dear Erik:  Would you be willing to post your subroutines on iMechanica?

Perhaps we could adopt an open-source approach.  I'm not sure exactly how these projects get organized, but there are plenty of great examples out there on the web.  That way, someone could choose to contribute what they already have made, and the community could improve on it, provide independent validation, etc.   Could iMechanica host space for such a repository? 

Arun Prakash's picture

This is a very nice suggestion. It would be nice to have a repository of user material routines in imechanica. An alternative could be sourceforge, which would not just be a repository, but also a development hub, where different people could help develop/customize the UMAT as required.


It seems that quite a few people are interested in sharing subroutines. It would be nice if people would post their UMATs which they would like to share (or perhaps provide a link for the same) .   


I realise this is an old post. 

I am in the process of developing my first user element in Abaqus, and finding it rather difficult. 

The application requires modelling muscle behaviour. Hence a spring model (linear at this stage) with 3 elements: 2 in series with each other, and these 2, parallel to the third was implemented.

My problem involves static condensation of the stiffness matrix. (similar to mixed-methods)

Any advice/assistance would be great. 





attash3099's picture


i need uel subrotine(dynamic beam code)

please send for me uel tutorial about UEL


nice and good thought

Arun Prakash's picture

It seems that there are two different threads with the same topic and both have interesting comments. It would be nice if the moderators could merge these two.

Zhigang Suo's picture

We have not found a way to merge two threads.  However, a user can always write a comment to link the two threads, as you have just done.  More broadly, this method is probably more effective than merging threads, because the method can relate even loosely related threads.

thanks, interst. Spore

uubird's picture



can anyone upload ABAQUS umat for MTS model? I wrote one but I am afraid is not correct because it doesn't match the experiment result quite well.



In my graduation project, the instructor wanted from me to work the UMAT Subroutines in ABAQUS. I am looking for someone who can help to me about this issue. I really got an urge for completin this project and really need of somebody's help.

From now am thanking to the ones who will help to me.

I am lettin my msn adress for the ones who want to contact with me...:

I can't agree more! Furthermore, this problem doesnot exist on ABAQUS uniquely. This same idea can be also helpful to ANSYS, Nastran and so on!

Sprunger's picture

Hello, Re: Sharing Abaqus UMAT and VUMAT subroutines.SIMULIA has an “Abaqus Process Automation Portal” portal provides a mechanism to the SIMULIA (Abaqus) community for collaboration through discussion forums for GUI customization and scripting techniques. Users can upload, download, and share solutions and process automation components.

The Abaqus Process Automation Portal is approaching its 5 year anniversary and many engineers found it to be a useful resource. Almost 3000 registered members currently exist and this number continues to increase steadily.

Regards, Rene 


llguang's picture


Hello,everyone, i just have a question, if i have the UEL,can i use it to simulate the dynamic problem?

I want to model a reinforced concrete cylindrical shell with abaqus using concrete damaged plasticity model. I want to model meridional and circumferential rebars in cylindrical shell in abaqus cae. I want to know which rebar geometry (angular or constant) I should use for modeling of both meridional and cirunferential rebars.Can I plot crack pattern for reinforced concrete shell in abaqus?please guide me 

Ahmad Rafsanjani's picture

You can find a good tutorial for UMAT in the following post

Hi everyone.

Im a phd student who is learning ABAQUS for a short time, now I am using UMAT subroutine to define my own material, but when I finished coding fortran file with .f extendings, I dont know how to connect my subroutine file and input file(.inp) with ABAQUS.

When I write C:\Temp> abaqus job=test user=...

Input File                        : ...

Abaqus error : the following file(s) could not be located:...inp

I could not find any information in abaqus tutorial.


it seems that the input file (test.inp) is not in the same path with that that you are running ABAQUS in.

hi guys, im trying to make a UMAT code for Plasticaly Graded Materials, is there anyone have experience with the same topic?



i am new in abaqus umat. before i worked in windows xp but i have a windows 7 and os 64 bit

i installed abaqus 6.10, intel visual fortran 11.1.60 and vs2010 but i can not configure them, can anyone help and give the [procedure to configure abaqus in fortran

thank`s a lot 


In my simulation with vumat, the temperature will be higher than its melting temperature on several integral points, how sholud i update the stress and strain in that situation? Give them zero value? I think pressure shold be exist, but how to calculate it?

anyone can help me with this problem?


Gurudutt Chandrashekar's picture

Hello everyone,

I am a student of Msc in Computational Engineering at Ruhr University Bochum in Germany. I am implementing a new piezoelectric element as a part of my Master Thesis. I am supposed to write a UEL subroutine and implement it in abaqus.I would like to know whether anyone has written subroutines in Fortran 90 or further versions and implemented them in abaqus. All the Help documents in Abaqus have subroutines written in F77.  If it  possible in F90, then it will be very helpful if someone could point out the necessary modifications that one has to make to the F90 file to  run it in Abaqus. I think many ppl have this question .I myself have been trying to figure out an answer from past 1 month but could not find a convincing solution anywhere.

I have never done this, but I don't think there should be a problem in general. Make sure you've configured Abaqus to use a Fortran 90-aware compiler and if you're lucky it will Just Work.

I have used Abaqus user subroutines written in C++. In that case I just provide Abaqus the object file instead of the source code and it links it without caring what generated it. When you provide user=something, Abaqus automatically checks to see if it's source code or precompiled.

Gurudutt Chandrashekar's picture

Hi Graham,

Thanks for ur suggestion. I have made all my subroutines in F90 . I got a doubt  when i read in one of the abaqus documents on "Writing User Subrotines in Abaqus"  that we have to specifically use F77 conventions. I should now figure out a way of configuring my Abaqus to use my F90 subroutines.

A hack might be to rename your F90 compiler to your F77 compiler.


I am an MSc student and my dissertation topic requires me to model a split hopkinson pressure bar.I am modelling a SHPB assembly together with the specimen, and I am having trouble with placing the loads,creating the boundary constrain between the incident, transmitter bar and the specimen. I created the interaction between the assembly with a suppressed TIE, and after loading and submitting the job I kept on getting errors. The load did not transmit to up to a quarter of the incident bar.
I neeed suggestions on how to apply the attachment between the separate components that make up the assembly, how to time step properly and the loading sequence.
I made use of ABAQUS Explicit

badrinath's picture

Badri Veluri

Hi Erik,

Can you please share your knowledge as I am working on building the UEL for a cohesive zone model under fatigue loading. So that if I refer your knowledge base it would help me in getting some good point.


 Thanks in advance



Hi freinds,

I 'm preparing PHD Tthesis in mechanics feild, i arrived  to modelise densification of powders on Abaqus, but i can't modelise the thermal debinding behaviour (loss of polymer weight vs temperature). Anyone have an idea about this.

I have comsol and Abauqs

Thanks for your answers !


Hi dear friends
After much effort and time and confusion with VS, Intel Fortran Compiler and Abaqus (installaion,version,paths,...), I finally could run an analysis with UMAT subroutine successfully. This run and also some others were all extracted from Abaqus Documentation.
Now i wanna implement an analysis by my own. So i read the documentation and wrote a UMAT file by making my own changes in documentation file. From algorithm point of view, the file is correct. But i wanna know how i could link the Fortran file and Abaqus file. I mean how Abaqus realizes the Fortran file is the compatible subroutine for the model? Logically, a subroutine file shouldn't be compatible with another abaqus file. (I don't know if i could explain my problem !)

By the way, i've figured out that one of the links is about Property module. I should define a USER MATERIAL in "Edit material", and i also should define variables in there compatible with PROPS ( I'm mixed up right now!!!).

PLEASE guide me.

Hi friends,

I need umat subrotine(Lemaitre damage model)

please send for me if u have it.Thank u. 

My email:

Over time I have written UMAT and UEL routines for progressive failure analysis, higher-order plate finite elements, specialized adhesive elements, and a score of decohesion finite element incorporating various coupled cohesive zone models to joint 1-D beam, and 2-D and 3-D continuum elements. kidney stones

shahram's picture


I am doing approximately the  same thing and I have some problems, could you please send me an example of what you have done?

I have no idea where should I start actually to write a failure analysis code. :(

Thank you 

kawaminaii's picture


I have compiled a fortran code and used it as a subroutine but I have resieved this errore:

Error in job Base_subro: Problem during linking - Abaqus/Standard User Subroutines.   This error may be due to a mismatch in the Abaqus user subroutine arguments.   These arguments sometimes change from release to release, so user subroutines   used with a previous release of Abaqus may need to be adjusted.

Wht should I do with it?


    I have met this problem before.Maybe there is some error  in your Fortran file,you can check it out especially for the program grammar.

   Good luck!

Hi,Does any one worked with Vumat subrotine for anisotropic hardening ?If you have any sample code for anisotropic hardening just let me know.

Hi freinds

days ago i write a subroutine for contilever beam whit various young is true for 2D beam but for 3D beam write error.

ERROR:Too many attempts made for this increment

Abaqus/Standard Analysis exited with an error - Please see the  message file for possible error messages if the file exists.

i change the load and step increment freuqently but its not true yet.i dont khnow where is Incorrec.

please help me.............


mesbah.unm's picture


I have recently purchased Intel Visual Fortran 11.1.072 and Visual Studio 2010. I am currently using ABAQUS 6.10ef2. After installing this according to the following sequence,

1. Microsoft Visual C++

2. Intel Fortran


my verification shows that the status of both visual studio and fortran is failed. Is there anybody who can help me in this regard?


Mesba Bhai.....I am from 03 batch....BUET. I did my Msc and BSc thesis using ABAQUS. where are you now? you are also my FB friend!.

akjon BUETIAN dekhe valo laglo:-)

what are you doing using ABAQUS?


you can personally communicate me

I will  be glad to have a response frome you.


Thanking You
B.K.M.Ashraful islam SAJIB



I need model of femur abaqus and how static loading on cortical bone
What is the relationship between abaqus and umat
thank you
My email:

Hi friends,

I need umat and Vmatsubrotine (GTN model)

please send to me if u have it. Thank u. 

My email:


please i need how to configure abaqus with fortran and I install abaqus 6.11 and FORTRAN 10 and microsoft visual 10 program when run under abaqus there is a error: Problem falling on compilation - ifort.exe not found in PATH.
how will be please
thank you

I need to used the umat Abaqus for simulate the thermomechanical effect in composite plate
please help me for this

I need to used the umat Abaqus for simulate the thermomechanical effect in FGM blade
please tell me if you get any help or guide lines


if i have the equation as below

σ(n+1) = σ (n) +delta(σ)

 if i get the initial stress, σ (n) and the differential stress, delta(σ), ist i need the call the (KEY . EQ. 11) THEN ....

or i need to make the full stiffness matrix and then

do i = 1, ntens

stress (i) = .......


ur help highly appreciated.   


hy every one well iam working abou a project of active vibration control of fgm plate and iam using ansys but i have a probleme in how to do a 3d fgm plate in ansys

so i need a apdl code

my email: 

any help  thank's 


I am starting an academic project and i aim to model cracking of an early age concrete with abaqus..

knowing that I am just starting and this is a classical problem, i need to know if the material behavior already exist in Abaqus ? if not, does any one knows how to create the corresponding UMAT ?



awanibhu's picture

Why nuclear grade graphite is used as a moderator better than other moderator and why it is economical than other.

Hello everyone. My name is Alex, I'm a doctoral student at the University of Salamanca. I am working with an axisymmetric model, steel material. I intend to obtain the residual stresses produced when hits on a major load model (Tera Pascals) for a very short time (picoseconds). The constitutive model is working: Elastic isotropic, platic johnson-cook, eos us-up. First work with Abaqus Explicit Abaqus and then Standart (I import tensions), but Abaqus CAE / Standart (version 6.13.2) does not allow plastic or Johnson-cook or EOS. To get "the spring back and therefore residual stress. 

Please could someone help?. 



Hi all,

I want to learn the steps of the simulation program on the subject of Abacus about (split Hopkinson pressure bar compression)with and without effect temperature .This device gives high strain rate impact where it two long bars between it put the specimen and then strike with short bar. Our data collection system consists of a pair of strain gages midway along each bar. Data is recorded and presented as the initial impulse wave passes the incident bar strain gages, a reflected pulse from when the wave hit the end of the bar and reflected back to the incident bar strain gage and finally a transmitted pulse which is what remains of the wave after it has passed through a sample. In all cases the voltage reading from the strain gages should return to zero once the wave has passed. However with the new bar the signal lingers after the wave has passed a gage.
my version 6.12
Please anyone can help me

<p>I am new to abaqus I am writing VUMAT for neo hookean material. earlier I done UMAT and it run successfully&nbsp;</p>
<p>but while running VUMAT it is giving error that "cannot open include file '' ". I am giving command correctly&nbsp;</p>
<p>is there any other possibility. can input file error cause this kind of error?</p>

Add my 2 cents here:

Can somebody upload a UMAT for Schapery Material Model (1-D or 3-D)of non-linear viscoelastic materials?

ntuecd's picture


can anyone explain me how ABAQUS compute the PEEQ, PENER, and SENER for a elasto- plasticity problem. Second thing i want to know how does ABAQUS compute the equivalent strain for such problem.  It will be really helpful if someone help me out or give me some link.


Thanking you


Subscribe to Comments for "Sharing ABAQUS UMAT and VUMAT subroutines"

Recent comments

More comments


Subscribe to Syndicate