Dear all,
I’m trying to model multilayer soil in abaqus. To use the SaniSand umat in model I use the following parameters but no plastic strain at all. The unit system used is us(in). Are the input parameters incorrect?
Could anyone kindly help me solve this problem? Thanks a lot!.
problem with sanisand umat
Related Articles

UMAT and oneAPI compiler issues. (SOLVED…i hope :p)
anis kheffacheFinite Elements7 Views 1,578
Dear all, A while ago, a person asked about the hypoplasticity umats, the problem was that the umats were not working when compiled using the oneAPI fortran compiler, an error […] 
Modeling of triaxial test sample by using hypoplastic subroutine
Naif AlsanabaniFinite Elements7 Views 3,929
I tried to model sample of triaxial test in Abaqus program. I used hypoplastic subroutine. The parameters used in hypoplastic model was Hochstetten sand the given in website “https://web.natur.cuni.cz/uhigug/masin/plaxumat/node5.html”. The […] 
DIANA 10.5 new features
jesus mediavillaFinite Elements0 Views 4,146
DIANA 10.5 was released in 2021. It features many powerful capabilities to model soilstructure interection. See more in https://dianafea.com/2021DIANA105Release 
Few more Question on calculating lateral displacement of pile using hypoplastic clay model
Muhammad AhmadFinite Elements0 Views 3,850
Dear Professor David, Thank you very much for your previous clarification. its really helps. In continuation to my previous , there are few more question, i want your kindself to […] 
hypoplasticity model with Abaqus umat on Brno clay
mahdie masoudiFinite Elements5 Views 4,044
Dear all, I am a Ph.D. student at Kharazmi University. Recently I have been trying to simulate a triaxial compression test using the hypoplasticity model for clay using ABAQUS. The […] 
Using pore fluidstress element in SandSand
Kassem DibFinite Elements1 Views 3,457
Hello everyone, Is it possible to use pore fluid stress element in abaqus standard when soil is modeled using SANISAND constitutive model ? (UMAT) 
Problem in GRPS embankment simulation In ABAQUS
Bidur PathakFinite Elements0 Views 3,262
Hi, I am doing Plane strain 2D simulation of GRPS embankment in ABAQUS, i encountered a problem that results of with geogrid and without geogrid both are showing same deflection. […] 
finite element modeling of ISAHP for sand material
arad ahmadi radinFinite Elements1 Views 3,331
On the simulation of the cyclic mobility effect with an ISA hypoplastic hypoplastic hypoplastic model, is it possible to simulate 3D liquefiable sandy soil by available umat file which is […]
Who is Online
No one is online right now
Search SoilModels Website
Recent posts
 UMAT for CreepSCLAY model. 1.3.2024
 Cyclic tests with Triax element test driver 14.2.2024
 SUMMER SCHOOL ‘Numerical Modelling in Geotechnical Engineering’, Innsbruck – July 22nd26th, 2024 1.2.2024
 COURSE IN SOIL MODELING – NTNU, Trondheim – October 14th to 18th, 2024 21.12.2023
 MSE walls design in Plaxis 11.12.2023
 sand liquefaction modelling in Anura3D 8.12.2023
 BCV bentonite experimental and modelling datasets 14.11.2023
 Challenges in simulating geomechanical models using MohrCoulomb with Tension CutOff in Abaqus 20.10.2023
 Analys soil liquefaction with Abaqus under blast load 19.10.2023
 Hardening soil Model UMAT 13.9.2023
 Problem with VUMAT interface 20.8.2023
 UMAT and oneAPI compiler issues. (SOLVED…i hope :p) 5.8.2023
Recent Comments
 Abhay Pratap Singh on Cyclic tests with Triax element test driver
 Giovanni Ciardi on Cyclic tests with Triax element test driver
 Konstantinos Chatzis on Cyclic tests with Triax element test driver
 Ismail Khan on Cyclic tests with Triax element test driver
 Kanika Lamba on Cyclic tests with Triax element test driver
 Arie Koot on How to model the settlement in soil due to water drawdown.
 Konstantinos Chatzis on Cyclic tests with Triax element test driver
 Jose Duque on Cyclic tests with Triax element test driver
 MohamadReza Kamali on How to model the settlement in soil due to water drawdown.
 Leo Alibert on MSE walls design in Plaxis
 Ignacio Zuloaga on MSE walls design in Plaxis
 Giada Orlando on Problem with VUMAT interface
 Michael Spyridis on Problem with VUMAT interface
 Ahmad Moeineddin on Problem with VUMAT interface
 Nikolay Suyo Calla on Hardening soil Model UMAT
 Johan Clausen on Challenges in simulating geomechanical models using MohrCoulomb with Tension CutOff in Abaqus
 Zhenyu Liu on Problem in using Hypoplastic model for sand in Abaqus.
 Dorsa Shadlou on SANISAND for FLAC3D Download
Hey Yousef
The parameters look generally ok. Rather close to the original calibration of Dafalias & Manzari for Toyoura Sand. If you are using US units check if your G0 parameter, which is indeed unitless, is still 125 or is it different value.
Good luck
Piotr
Hello
Thank you for replay
G0 is a dimensionless constant, but by changing its value, PE is STEEL equal to zero.
Hey Yousef
As mentioned by Arie, PE is not calculated by the umat. It is calculated only by the constitutive models implemented in Abaqus. If you want to see if plastic strains occur from the umat, then read umat carefully and try to ask for your own output in log file. I am not sure if you can ask for the actual value of the plastic strain.
Piotr
What do you mean “no plastic strain”?
Sanisand umat does not calculate plastic strain: it is a kinematic hardening model with one backstress.
It only needs and calculates plastic strain rate.
Thank you for replay
PE in visualization always is zero
though with chang in load and b.c
is this correct؟
Or do you mean you miss plastic strain input?
Hardening behavior is described by the hardening parameters
i set depvar=36
and check pe in field output in step
and select umat in job
Something else needs to be done؟
select pe has no use: it is not calculated by sanisand.
You can select SDV output but then you get all 36 statevar. (you cannot select a few)
It will make your odb larger but it gives you more information.
c 1 … alpha_11 back stress, orientation of yield surface cone
c 2 … alpha_22
c 3 … alpha_33
c 4 … alpha_12
c 5 … alpha_13
c 6 … alpha_23
c
c 7 … void void ratio
c
c 8 … Fab_11 fabric tensor z
c 9 … Fab_22
c 10 … Fab_33
c 11 … Fab_12
c 12 … Fab_13
c 13 … Fab_23
c
c 14 … not used
c
c group 2: memory variables for shear reversal (SR) and other purposes
c
c 15 … alpha_sr_11 alpha value at stress reversal points (discrete update)
c 16 … alpha_sr_22
c 17 … alpha_sr_33
c 18 … alpha_sr_12
c 19 … alpha_sr_13
c 20 … alpha_sr_23
c
c 21 … not used
c 22 … not used
c 23 … not used
c 24 … not used
c 25 … not used
c 26 … not used
c 27 … not used
c
c 28 … not used
c
c group 3: variables saved for post processing or other purposes
c
c 29 … pore excess pore pressure (undrained case)
c 30 … p mean effective stress
c 31 … q deviator stress
c 32 … z Lode parameter (cos(3theta))
c 33 … dtsub suggested size of first time substep
c 34 … nfev number of function evaluation
c 35 … not used
c 36 … not used
c