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

new hypoplasticity laws in Tochnog Professional
dennis roddemanFinite Elements0 Views 146
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 […] 
Presentations on hypoplasticity applications, D. Mašín, NGI 2017
David MašínFinite Elements0 Views 1,054
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 […] 
How to replicate results for vertical cut in FEM program
Amine AboufirassFinite Elements0 Views 165
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 […] 
Tochnog Professional 100 element free version
dennis roddemanFinite Elements0 Views 192
The Tochnog Professional finite element program is now available as free version with 100 elements. It can be downloaded without registration from http://www.feat.nl. Several commercial and free pre and postprocessors […] 
Interface to study soilstructure interaction in Plaxis 3D (Clay Hypoplasticity)
Leon CortesFinite Elements6 Views 323
Hello, I am trying to simulate the behavior of helical anchors (also called screw anchors) using PLAXIS 3D and a clay hypoplasticity constitutive model. The helical anchor shaft is simulated […] 
Calculation of coordinates of climbing formwork using hypoplasticity and pressuremter tests
Christian SchwabFinite Elements0 Views 89
Our article published in the latest issue of Ground Engineering 10/2017. Great combination of conepressueremetertesting (CPM) by FUGRO and advanced 3D Fem modelling by GEOLINK using hypoplastic constitutive equations to […] 
SANISAND dll
Emilio BilottaFinite Elements2 Views 225
Hi there, I’ve used successfully used SANISAND .dll with Plaxis 2016, but after switching to 2017 and now 2018 the .dll seems not to work (the parameter fields appear without […] 
Hypoplastic model for thermal cycles in clay (Abaqus)
Sherif GomaaFinite Elements2 Views 123
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 […]
Search SoilModels Website
Recent posts
Recent Comments
 David Mašín on SANISAND Abaqus UMAT and Plaxis implementations
 Reza Foroughi on SANISAND Abaqus UMAT and Plaxis implementations
 Frank Richter on Hypoplastic model for thermal cycles in clay (Abaqus)
 Frank Richter on Hypoplastic model for thermal cycles (Abaqus)
 David Mašín on Hypoplastic model for thermal cycles in clay (Abaqus)
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