You are here

Sharing ABAQUS UMAT and VUMAT subroutines

Primary tabs

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

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 http://web.utk.edu/~ygao7/publication.htm

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.

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(baishalang2002@163.com). Thank you very much.

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

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

I am a new learner of CZM. I know that you are willing to share a subroutine about the cohesive zone element. I have searched your webpage, but I couldn't find it. Could you send it to me? My email address is E0256440@u.nus.edu.

I am looking for a UMAT or VUMAT to model Orthotropic Material behaviour. There is by default the Hill's criterion via Abaqus CAE for anisotropic material modeling, but I want to develop something further for my final model.

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.

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?

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) .

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.

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...:

Hello, Re: Sharing Abaqus UMAT and VUMAT subroutines.SIMULIA has an “Abaqus Process Automation Portal”http://www.simulia.com/PAPortal/This 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.

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

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.

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?

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.

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.

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

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.

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.

Hi dear friendsAfter 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!!!).

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.

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.

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.12Please anyone can help meregards

<p>HI&nbsp;</p>
<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 'vaba_param.inc' ". I am giving command correctly&nbsp;</p>
<p>is there any other possibility. can input file error cause this kind of error?</p>

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.

I am going to use two different UMATs for my model. The first one defines super-elasticity and the second one defines bone remodeling proces on my model. My question is, how can I combine them and use them simouloltaniously?

I want to implement the multi-scale model for ceramics under ballistic impact. My reference paper: Dynamic fragmentation of brittle solids: a multi-scale model,Christophe Denoual, François Hild. I have studied the Dassault systems 'Writing user subroutines in ABAQUS'. I have read the examples also. But the model I have has a new parameters which are not similar to the example models. I want to do Abaqus/Explicit analysis.If any one can help me regarding this it would be great. Please read the paper and give me some idea of how to implement in Abaqus using Subroutine VUMAT.If you have the subroutine it would be great help if you please send the subroutine to my email Id:kosuribs99@gmail.com

Dear allI wanna to change an existing material model so that some of it's parameters can be changed as a function of strain-rate. for example, i want to alternate flow stresses in simple plastic material model of abaqus. it should be noted that i have a function for strain rates that has abnormal relation with strain-rates.

I has written a subroutine for the 'probabilistic damage model for the dynamic fragmentation in brittle materials,hild,advances in applied mechanics,2010'..The procedure I followed in UMAT subroutine is as follows:

1.I calculate the principal stresses.2.Randomly generate a failure stress according to equation(84) in 2010 paper.3.For each principal stress I calculate the modified defect growth density given in equation(83) in 2010 paper.4. Corresponding to each principal stress a damage value is found according to the conditions given in equation(81) in 2010 paper(NOTE:Anisotropic damage model is used in three orthogonal directions D1,D2,D3 )5.I update the damage tensor/compliance tensor using (33)(JPEG file is attached)as provided in multi-scale model paper(2002).6.I update the stress and repeat the steps:1-6.

NOTE:1.The material is isotropic elastic needs only E,poisson's ratio,material parameters.The other material parameters like weibull parameters are also given. Refer the papers provided for clarification.

2.I am dealing with Impact problems so I need VUMAT

I have a doubt in VUMAT if I calculate compliance tensor(as show in attached file) at each step and inverse it to obtain the constitutive tensor and then multiply with strainInc to obtain stress increment.I will not update Jacobian.Can I do this in VUMAT to calculate stressNew?.Is my subroutine right?.Please check the VUMAT.txt file.Am I updating the stress correctly?

can any one give me an idea how to simulate a sample specimen while subjected to tensile load and fracture or shear out at the middle of the specimen at 45 degree on Ababqus CAE. please if you have the procedure to simulate it send me throgh my email

I am a Ph.D and also a new user with subroutine and abaqus. I hope some one can help me, I have a hydrodynamic problem I simulate seal (elastomer -hyperelastic material) in contact with fluid wich gouverned by Reynolds equation (I resolve it in subroutine). The deformation of seal is great wich stop simulation don't converge.

do you have idea how to resolve this problem, how to optimase time of calcul.

I want to model the FG Material in VUMAT. For this I need x and y components of each material point. But Abaqus gives me only one x and y per each nblock. While each nblock is composed of multiple material point.