Dear all
I was simulating a pile in a fully drained soil subjected to horizontal loads using ABAQUS (as shown in Figure 1) and noticed from the monitor that the calculation did not continue for a long time and when I interrupted the job I got a very large “.log” file (100300g).
I have looked at the “.log” file and found that it is a repeated error message (as shown in figure 2), what is the meaning of this message and is there any way to avoid it?
Thank you in advance for your help!
An error in the SANISAND umat runtime
Related Articles

Hypoplastic model for thermal cycles in clay (Abaqus)
Sherif GomaaFinite Elements2 Views 4,829
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 […] 
Tochnog Professional Page moves
Dennis RoddemanFinite Elements1 Views 2,775
Hi All, Starting from Saturday January 26 2019 the Tochnog Professional page moves to http://www.tochnogprofessional.nl Best, Dennis Roddeman 
Boundary value problems for granular materials subjected to earthquake loading
MOHD SAQIBFinite Elements1 Views 2,800
Dear All, Hope you all well, I have been searching for some good publications on the Boundary value problems on Foundation system with granular materials subjected to earthquake loading or […] 
How to replicate results for vertical cut in FEM program
Amine AboufirassFinite Elements0 Views 2,976
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 […] 
Interpreting results of sand hypoplastic model with Abaqus VUMAT
arsalan AhmadFinite Elements0 Views 3,092
Dear concern, Recently i have modelled the soil pile interaction using Abaqus CEL method. The hypoplastic behavior of the soil is model using user subroutine VUMAT. The properties of the […] 
Modeling of granular soil by using UMAT hypoplastic in Abaqus
Naif AlsanabaniFinite Elements5 Views 3,500
Dear all I am using the UMAT of hypoplastic model for simulation soil (3D modeling) under cyclic loading. The analysis was completed at geostatic step. But for consolidation step (soil, […] 
Presentations on hypoplasticity applications, D. Mašín, NGI 2017
David MašínFinite Elements0 Views 4,497
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 […] 
PM4Silt software
Daniel ValenciaFinite Elements5 Views 2,953
Hello everybody, Does somebody know a FEM Software with the implementation of PM4Silt model’s routine? I mean, one different from FLAC. Thanks for your help.
Who is Online
No one is online right now
Search SoilModels Website
Recent posts
 Presentations from Prague Geotechnical Days 2023 including 29th Prague Geotechical Lecture by Yannis Dafalias 25.5.2023
 Problem in using Hypoplastic model for sand in Abaqus. 17.5.2023
 An error in the SANISAND umat runtime 8.5.2023
 Htypoplastic Model in Miidas GTS NX 25.4.2023
 set cavitation limit in ABAQUS/STANDARD when using SANISAND umat 2.4.2023
 FE formulation and software 25.3.2023
 Multiobjective analysis of the Sand Hypoplasticity model calibration 24.3.2023
 Problem with SANISAND UMAT 7.3.2023
 Prague Geotechnical Days 2023 (22nd – 23rd of May) with Prague Geotechnical Lecture by prof. Yannis Dafalias 28.2.2023
 What is an appropriate damping model for dynamic liquefaction analysis of tailings storage facility in FLAC? 10.2.2023
 Abaqus vs Plaxis – Local Convergence Tolerances 2.2.2023
 Zbraslav sand – Cyclic tests by Charles University 17.1.2023
Recent Comments
 Kassem Dib on Problem in using Hypoplastic model for sand in Abaqus.
 Isma Khabis on Problem in using Hypoplastic model for sand in Abaqus.
 Isma Khabis on Problem in using Hypoplastic model for sand in Abaqus.
 Kassem Dib on Problem in using Hypoplastic model for sand in Abaqus.
 Zhentao Liu on An error in the SANISAND umat runtime
 Arie Koot on An error in the SANISAND umat runtime
 Zhentao Liu on An error in the SANISAND umat runtime
 Arie Koot on An error in the SANISAND umat runtime
 Zhentao Liu on An error in the SANISAND umat runtime
 Arie Koot on An error in the SANISAND umat runtime
 Zhentao Liu on An error in the SANISAND umat runtime
 Arie Koot on An error in the SANISAND umat runtime
 Zhentao Liu on An error in the SANISAND umat runtime
 Zhentao Liu on An error in the SANISAND umat runtime
 Zhentao Liu on An error in the SANISAND umat runtime
 Zhentao Liu on An error in the SANISAND umat runtime
 Arie Koot on An error in the SANISAND umat runtime
 Arie Koot on An error in the SANISAND umat runtime
Sorry I got the order of the pictures wrong. Figure 1 is the “.log” file that reported the error, and Figure 2 is the piled soil model
substepsize is too low < 1e18
problem cannot been solved.
perhaps model is wrong
Dear Arie：
Thank you so much for helping again! Does this error mean that the analysis step is too small?
I’m curious why my ABAQUS doesn’t break the calculation but keeps repeating the attempt
It do not know which units you have used, but if the displacement is in meters the displacement is large (20cm), and it lookes there is a gab between pile and soil. IF you have used fully drained and cohesion is zero this is unstable and because you have used an implicit methode it probably not converts. Try to increase cohesion (pshift).
The units in the simulation are m. I am simulating the horizontal displacement of a pile in a sandy soil.
Shouldn’t the cohesion c=0 and parm18 be taken as 0 for sand?
Also, is 20cm a large displacement? I have seen in many papers that the horizontal displacement of the mud surface in the simulation results reaches D/10, for my simulation (pile diameter D=10m) this displacement should reach 1m.
Thanks again for your reply!!!
To add to this, the monitor shows that the incremental step size is not reduced, it stays at the maximum step size I set, but no further calculations are performed, just the machine idles (100% CPU usage) while generating a large LOG file.
If convergence is a problem, ABAQUS should show 1U,2U… Until the calculation breaks isn’t it?
Why does ABAQUS not show 1U.2U etc. after calling SANISAND’s UMAT, but stays at one analysis step and repeats the LOG file writing?
I’m running a full drained simulation, i.e. can’t c take 0 in the case of full drained?
There is tension in the soil: because c=0 it is not possible.
DT_k is reduced to DTmin and DTMin is very low.
Send .inp, .cae, .for in a zipfile to koot@zonneveld.com and I can see.
else
c … if pp_2 lower than zero
DT_k=DT_k/four
if(DT_k.lt.DTmin) then
DT_k= one – T_k
mario_DT_test=1
write(*,*) ‘mario_DT – pp_2<0 – mario9=one, pp_kk=',pp_kk
end if
endif
c
Dear Arie
Thank you very much! I’ve noticed what this code does.
After reading your reply yesterday, I set parm18=0.1 and ran the simulation again, but the simulation did not achieve the target of 1m horizontal displacement of the mud surface either.
I have send my .inp, .cae and .for to your email, Thank you in advance for your help!
use dynamic implicit instead of static general for displacement
take 1 step after gravity.
time for example 100sec and minimal increment 1 sec
Thank you for your advice! I have tried to replace the steps from static to dynamic implicit. However, after the geostress step is completed, ABAQUS begins to idle rather than continue the dynamic implicit step.
The error is still the same as the Fig.1 in my original post.
I am now looking for which step of the UMAT run will cause pp _ kk to become NAN.
geostep is not correct: only soil stresses are calculated: state var are zero (included void)
for geostep you need umat or sigini+sdvini
if you use umat some code have to be added to add stresses to avoid lot of vertical displacement.
But it will not solve your problem that model gives umatrrors and not convert.
Thank you for reply, Arie! I will correct GeoStep with the guide of your advice. As for the errors, I am still looking for the code causing this bug. There is no division operation in the calculation of PP_KK, I don’t think the result of it could get the result of NAN. If you know, can you tell me what other cases would cause PP_KK=(SIG1+SIG2+SIG3)/3 to be calculated to get the answer of NAN?
Just a general question: Why you want to use sanisand umat?
I am study in the direction of constitutive models. And SANISAND model is the first sand model that I have been exposed. Since its code is easily available (posted on this website), I chose this model as the first step in learning the constitutive model.