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

Factor of safety in Plaxis 3D using hypoplastic model for sand
Arno CrousFinite Elements1 Views 3,838
I am currently modelling a test conducted in the geotechnical centrifuge facility at HKUST. The test consisted of a sand slope with a rising water table. The slope height in […] 
Plaxis 3D Tunnel
Mesut DemirFinite Elements0 Views 3,636
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 […] 
Tochnog Professional becomes free
Dennis RoddemanFinite Elements0 Views 3,495
After several decades of development we are happy to announce that the Tochnog Professional finite element program has become free. Both for academic and commercial use. We are extremely grateful […] 
Some confusion about the parameters in Abaqus umat of Hypoplastic Modeling
Joe ChuFinite Elements1 Views 3,365
Dear all, When I check the umat file of sand hypoplastic modeling, I find there are two parametesr(p_t(2), bulk_w(15)) without any definitions. Could someone kindly explain the exact meaning for […] 
Few more Question on calculating lateral displacement of pile using hypoplastic clay model
Muhammad AhmadFinite Elements0 Views 3,872
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 […] 
a question on the clay hypoplasticity model with Abaqus umat
Tao ZhuFinite Elements7 Views 3,690
Dear all, I am a PhD student in Hohai University. Recently I have been trying to simulate a triaxial compression test using the hypoplasticity model for clay. The UMAT subroutine […] 
UMAT Implementation in Abaqus (Hypoplastic clay Subroutines)
Lovinder MannFinite Elements4 Views 3,793
Hi, I am working in Abaqus, aiming to model double porosity clay consolidation data by using hypoplastic clay model with meta stable structure. First, i just wanted to verify basic […] 
Calculation of coordinates of climbing formwork using hypoplasticity and pressuremter tests
Christian SchwabFinite Elements0 Views 3,301
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 […]
Who is Online
No one is online right now
Search SoilModels Website
Recent posts
 Degradation ( disturbance ) of hypoplastic clay 23.4.2024
 Degradation ( disturbance ) of hypoplastic clay 23.4.2024
 Paper of Coulomb, C. A. (1773) 15.4.2024
 ABAQUS UMAT of hypoplastic clay model 6.4.2024
 Prague Geotechnical Days 2024 “Geotechnical monitoring” and 30th jubilee Prague Geotechnical Lecture by prof. Eduardo Alonso 5.4.2024
 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
Recent Comments
 Ramon Varghese on Hypoplastic Interface Model: Does anybody implemented this model in FLAC3D
 Jh Xue on UMAT Implementation in Abaqus (Hypoplastic clay Subroutines)
 Santosh Katuwal on PM4Sand for Plaxis Download
 KHA DIDJA on SANISAND Abaqus UMAT and Plaxis implementations
 Debdip Das on SANISAND Abaqus UMAT and Plaxis implementations
 Ramon Varghese on SANISAND for FLAC3D Download
 Akhtyar Gul Shirzoi on ABAQUS UMAT of hypoplastic clay model
 KHA DIDJA on Animating Soil Models
 Jh Xue on Undrained Tests with Clay Hypoplasticity Intergranular
 Zhentao Liu on Problem in simulating CPT using SANISAND04
 Gertraud Medicus on ABAQUS UMAT of hypoplastic clay model
 Chen Zhiming on Download Package of Charles University Implementation of High Cycle Accumulation Model
 Konstantinos Chatzis on Cyclic tests with Triax element test driver
 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
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.