Dear Prof. Masin and Soil Model Community,
I am using HP clay subroutine in Plaxis 3D to simulate soil structure interaction problem. I want void ratio and OCR variation along the depth to back analyse centrifuge test. In Plaxis, either or OCR remains constant depending on input. Moreover, I found HP subroutine in plaxis calculate void ratio without considering kheppa line. Please see details and explanations in the attached images. My objective is to obtained initial void ratio and OCR variation along the model depth. Also, I want to confirm the position of e input in parameters in ln(1+e) – lnp space. Thanks
HP clay subroutine in Plaxis
Related Articles

Consolidation analysis using HP clay model in Plaxis
muhammad shakeelFinite Elements1 Views 3,421
Dear Prof David / SOil model community, I am performing coupled consolidation analysis for deep excavation in clay using HP clay model in Plaxis 3D. Performing the undrained analysis, the […] 
A problem of convergence from using Hypoplastic model in Abaqus
Joe ChuFinite Elements5 Views 3,044
Dear all, The sand hypoplastic model has been used to simulate the interaction between soil and pile. However, there is a bad convergence in the process of calculation. When I […] 
Calibrate hypoplastic model
sahar saabFinite Elements2 Views 2,857
I am trying to calibrate my sand using Excalibre. I tried some triaxial element tests using Abaqus. However the results are different between Excalibre and Abaqus. 
SANISAND umat issue
Philip AlkhouryFinite Elements2 Views 3,008
Hello everyone; I have implemented the Sanisand constitutive model (UMAT) in abaqus in the scope of performing a dynamic analysis on a 3D soilpile interaction problem. I have succeded in […] 
SANISAND dll
Emilio BilottaFinite Elements2 Views 3,131
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 […] 
How to replicate results for vertical cut in FEM program
Amine AboufirassFinite Elements0 Views 2,980
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 […] 
problem with sanisand umat
yousef zandFinite Elements8 Views 3,107
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 […] 
Presentations on hypoplasticity applications, D. Mašín, NGI 2017
David MašínFinite Elements0 Views 4,499
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 […]
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
 Francisco Mendez on Presentations from Prague Geotechnical Days 2023 including 29th Prague Geotechical Lecture by Yannis Dafalias
 Héctor Montenegro on Extended Mohr–Coulomb (EMC)
 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
Dear Muhammad,
Actually, there are these two options only, initialise through constant void ratio or constant OCR. By the way, these only work with “K0 initilisation procedure”. With gravity loading procedure, void ratio will redictribute during gravity loading. If you need more complex distribution of void ratio with depth, you will have to approximate it through using many soil layers which will have the same parameters, just different initial void ratio.
Regards David
Dear Prof. Masin,
Thanks for your reply.
May I know how the model calculate undrained shear strength (cu)? Is it based on initial void ratio like critical state frame work or some other way? My concern is, if void ratio is constant along the model depth how can we obtained increasing cu along the depth with increasing effective stresses having constant void ratio?
Thanks for your help.
Regards,
Muhammad
Hi Muhammad,
Yes, HP model calculates cu based on the critical state soil mechanics framework. If you set constant void ratio with depth, you get constant undrained shear strength with depth. You can either initialise OCR, than cu will be distributed over depth based on that OCR. Or, you need to define layers with different void ratios.
David
Hi both,
Does this mean that the undrained shear strength can therefore be calculated for the hypoplastic model using the same formula for the MCC model:
Su=M/2 exp^(((Γ1e)/λ))
Cheers,
Hashmi
Just with a difference that in the oridnary Modified Cam clay normal compression lines are considered as linear in p vs 1+e plane, whereas hypoplastic model adopts Butterfield’s formulation with normal compression line linear in ln(1+e) vs ln p plane, so the equation needs to be updated accordingly. David