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 in clay (Abaqus)
Sherif GomaaFinite Elements2 Views 205
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 […] 
liquefaction simulation with abaqus
Ruibin HouFinite Elements3 Views 195
Hello, everyone I am an abaqus user and developer. Recently, I want to simulating the liquefaction phenomenon with abaqus/standard or abaqus/explicit. As I know, this cannot be realized using the […] 
EasyTochnog input file generator
Dennis RoddemanFinite Elements0 Views 163
Hi All, Tochnog Professional now has an EasyTochnog extension, see http://www.tochnogprofessional.nl/easy_tochnog/page.php. With simple input data like pile diameter, pile length, slope geometry, sheet pile data, etc. a complete Tochnog Professional […] 
ABAQUS Meshtomesh solution mapping with Modified Cam Clay
Samanthika LiyanapathiranaFinite Elements0 Views 151
I am using the meshtomesh solution mapping with MCC model. The element type used is CPE8RP. My initial analysis runs without any issue and when the deformations reach a certain […] 
Plaxis 3D Tunnel
Mesut DemirFinite Elements0 Views 258
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 […] 
SANISAND dll
Emilio BilottaFinite Elements2 Views 287
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 […] 
Calculation of coordinates of climbing formwork using hypoplasticity and pressuremter tests
Christian SchwabFinite Elements0 Views 108
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 […] 
How to replicate results for vertical cut in FEM program
Amine AboufirassFinite Elements0 Views 205
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
 Mj Saffarian on Error in abaqus when use sanisand
 Mj Saffarian on Error in abaqus when use sanisand
 Hashmi Sohawon on General Questions Hypoplastic Umat on Abaqus
 David Mašín on General Questions Hypoplastic Umat on Abaqus
 Hashmi Sohawon on General Questions Hypoplastic Umat on Abaqus
 David Mašín on General Questions Hypoplastic Umat on Abaqus
 Frank Richter on Nonlocal regularization: Need VUMAT of MohrCoulomb model in abaqus
 Frank Richter on Error in abaqus when use sanisand
 Arie Koot on Error in abaqus when use sanisand
 Zhenhao Shi on Error in abaqus when use sanisand
 Mohamed Abdelmonem on PM4Sand (Boulanger & Ziotopoulou)
 Usman Pervaiz Khan on PM4Sand (Boulanger & Ziotopoulou)
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