Dear all,
After carefully reading the matlab drivers of hypoplastic models, I was wondering why there is not mix control undrained TX compression. I want to achieve the same result as Niemunis and Herle(1997) did. How can I make the value of q fluctuates within a certain range .When I make dsigma11dsigma33=DX , and DX=[30,60,60,60,60,60,60,60,60,60,60,60] , the result is clearly not right. Any suggestions or comments will be appreciated. Thank you!
Question on the Matlab driverhypoplastic models
Related Articles

Velacs Project Data
WaiChing SunConstitutive Modelling1 Views 3,423
Dear colleagues, I want to ask your help — does anyone know where we can find all the test data from the old Velacs Project? If so, that would be […] 
UDSM PLAXIS
nasim aftahiConstitutive Modelling2 Views 3,836
I implemented a unified dynamic soil model, but PLAXIS doesn’t calculate any strains, and it prints zero (or some times NaN amounts) into a text file. I’d like to know […] 
ExCalibre Automatic Calibration
Marco HoferConstitutive Modelling5 Views 3,511
Hello all, i have another question regarding the calibration process. Why do the initial void ratios differ between the input file and the output from excalibre? What is done during […] 
Cyclic simple shear test simulation with hypoplasticity
Piotr KowalczykConstitutive Modelling2 Views 3,948
Dear All, I am using sand hypoplastic model in my research and I have encountered some problem when looking into cyclic simple shear test in the range of small strains. […] 

lade umat file
sobhan abedinnejadConstitutive Modelling0 Views 3,522
Hi everyone, I work on soil constitutive models and I want to compare soil behaviours with various constitutive models but I can’t find Lade umat in the .for format. Can […] 
PhD course “Hypoplasticity for Practical Applications”
Hashmi SohawonConstitutive Modelling2 Views 3,763
Are the slides in regards to Stiffness Anisotropy for Clay Hypoplasticity for the PhD course available online? The available document {Mašín, D. (2015) Part 4: Determination of material parameters, PhD […] 
Spurious oscillations in calculated accelerations
Piotr KowalczykConstitutive Modelling2 Views 3,462
Hello everyone, I have just a general question how people approach the problem of spurious oscillations generated in dynamic analysis when using nonlinear soil constitutive models. I understand sometimes numerical […]
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
Hi C, the mixed control is not there predefined, but you can still add it easily into constraints.m file. Please study Janda, T. and Mašín, D. (2017). General method for simulating laboratory tests with constitutive models for geomechanics. International Journal for Numerical and Analytical Methods in Geomechanics 41, No. 2, 304312 which is available in SoilModels publication database. You need to control axial stress
S(6,3)=1;
and radial strains so that zero volume strain is achieved
E(1,1)=1;
E(1,3)=0.5;
E(2,2)=1;
E(2,3)=0.5;
E(3,4)=1;
E(4,5)=1;
E(5,6)=1;
David
Hi, Prof.David. Thank you for your reply! As you mentioned, I change the constraint. However, there is still something wrong with it. I check the EE after I did what you mentioned, I found that the volume strain isn’t zero in the whole process. I was a little confused about these problems. I want to achieve the same result as Niemunis and Herle(1997) did. They controlled q fluctuates within a certain range. Is it possible for us achieve this through our matlab driver? Looking forward to your response!
Hmm, not sure why, this sets ep22=ep33/2 and ep11=ep33/2.
You may also try
E11=E12=E13=1, which sets ep11+ep22+ep33=0.
E21=1, E22=1, which sets ep11=ep22
E34=E45=E56=1, which sets shear strain components to 0
S63=1, which sets axial (axis direction 3) stress control
This will, however, still not stop the test at predefined q, as radial stress will change during during loading through pore pressure generation. To achieve constant max and min q cyclic loading, you will have to manually set nstep for each load cycle, or add your own stop condition into the code, or better use some of more advanced element test software (such as triax or idriver).
Prof.David. in the f_hyp.m, after I make Laux = Lep and sigpresc = Laux*depsMNRNep*sqrt(depsMNR’*M2*depsMNR), the volume strain is zero, even though other results are not totally right in my views. While your file README.pdf clearly says that sigpresc = Laux*depsMNR rather than sigpresc = Laux*depsMNRNep*sqrt(depsMNR’*M2*depsMNR). I’d appreciate it if you could explain why sigpresc = Laux*depsMNR. Looking forward to your reply!
I see, I think the problem is how residual is calculated in NewtonRaphson iterations in Matlab driver. It is only calculated from residual stresses, which is insufficient for more complex mixedcontrolled tests such as those we discuss here. Correct procedure is defined in Janda & Mašín (2017), where the residual is calculated from S*dsig+E*depsilondy, see Eq. (12).
For you, that would mean either updating Matlab driver formulation or using more advanced element test driver, such as triax