Hi,
I want to change the friction angle with the plastic strain in Abaqus following mohr coulomb model.
Vermeer and de Borst (1984) proposed an equation for frictional hardening and I want to use that. In their equation, the shear hardening parameter obtains from below equation(first picture) and then the friction angle changes with the shear hardening parameter.
So I have to use the usdfld code in my model.
I code that to calculate the shear hardening parameter, and then I put it equal to FIELD1. I attached the usdfld code(second picture) to this post.
can anybody tell me what is wrong in my code?
Do everybody have Vermeer equation code in Fortran for using in Abaqus?
change the friction angle with the plastic strain in Abaqus following mohr coulomb model.
Related Articles

Hypoplastic model for thermal cycles (Abaqus)
Chiara IodiceFinite Elements2 Views 368
Dear all, I’m using the UMAT file in which the hypoplastic model for thermal cycles is implemented. In the Abaqus input file I have added the two lines necessary to […] 
fortran format
sobhan abedinnejadFinite Elements1 Views 60
hi, every one i have problem with fortran format my abaqus is 6.14 and this version accept (.for) format for fortran but some subroutines are in (.f ) format is […] 
Presentations on hypoplasticity applications, D. Mašín, NGI 2017
David MašínFinite Elements0 Views 1,397
Based on some requests, I share here two presentations on hypoplasticity applications. One is related to simulation of various boundary value problems, the other focuses on simulation of cyclic and […] 
UMAT for hypoplastic subroutine
Naif AlsanabaniFinite Elements4 Views 510
Hello all, Could you please anyone provide me umat for Hypoplastic model that used as subroutine in Abqus 6.14. With Best Regards Naif Mohammed 
Hypoplastic model for thermal cycles in clay (Abaqus)
Sherif GomaaFinite Elements2 Views 235
Dear all, I’m using the UMAT file in which the hypoplastic model for thermal cycles is implemented. I have done some runs for an axisymmetrical model for a single energy […] 
Plaxis 3D Tunnel
Mesut DemirFinite Elements0 Views 276
Hi Everybody, I would like to ask an error problem where I am getting while using Plaxis 3D tunnel V 1.2. Previously I was using the software with Windows 7 […] 
new hypoplasticity laws in Tochnog Professional
dennis roddemanFinite Elements0 Views 185
The 2014 Masin hypoplasticity for clays has been implemented in Tochnog Professional under guidance of Dr. David Masin. This implementation includes structure, intergranular strains and visco effects. Extensions to the […] 
How to replicate results for vertical cut in FEM program
Amine AboufirassFinite Elements0 Views 215
The equations for a vertical cut in cohesive material are given in several publications including Heyman (1973): The stability of a vertical cut (see attachment Equations1.jpeg) and Verruijt (2001) Soil […]
Search SoilModels Website
Recent posts
Recent Comments
 Hans Henning Stutz on Some confusion about the parameters in Abaqus umat of Hypoplastic Modeling
 Hans Henning Stutz on abaqus version
 Sobhan Abedinnejad on Hypoplastic rate independent and nonlinearity
 David Mašín on Hypoplastic rate independent and nonlinearity
 David Mašín on SANISAND Abaqus UMAT and Plaxis implementations
 MOHD SAQIB on SANISAND Abaqus UMAT and Plaxis implementations
 Sobhan Abedinnejad on abaqus version
 Hengyi Kang on Question on the Matlab driver – model hypo Cam Clay
 Hashmi Sohawon on Undrained Tests with Basic Clay Hypoplasticity
 David Mašín on Question on the Matlab driver – model hypo Cam Clay
 Hashmi Sohawon on abaqus version
 David Mašín on fortran format
Who is Online
No one is online right now
Use total plastic strain and Cauchy for 3D
strain=sqrt(0.5*(pe11pem)^2+0.5*(pe22pem)^2+0.5*(pe33pem)^2+pe12^2+pe23^2+pe13^2)
pem=(pe11+pe22+pe33)/3
for explixit pe23=array(5) and pe13=array(6)
for implicit pe23=array(6!) and pe13=array(5!)*
for plainstrain: pe13 and pe23 do not exist
*
for explicit vgetvrm gives you e12,e23,e13
for implicit you maybe get y12,y13,y23: you have to check this