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

Boundary value problems for granular materials subjected to earthquake loading
MOHD SAQIBFinite Elements1 Views 3,481
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 […] 
ABAQUS modeling: Warning about unconnected regions when using hypoplasticity subroutine
Jasmine LokeshFinite Elements0 Views 5,047
I ran a simple 3D model with symmetry boundary conditions first using elastic material and then as hypoplastic using a subroutine. However, in the former case, I don’t get the […] 
Few more Question on calculating lateral displacement of pile using hypoplastic clay model
Muhammad AhmadFinite Elements0 Views 4,073
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 […] 
Problem in GRPS embankment simulation In ABAQUS
Bidur PathakFinite Elements0 Views 3,478
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. […] 
In Abaqus how do we use *DLOAD for defining both userdefined pressure and body force simultaneously?
Ashesh ChoudhuryFinite Elements0 Views 385
In Abaqus how do we use *DLOAD for defining userdefined pressure and body force simultaneously? I want to userdefine both body fore and pressure. 
Geostatic method for Abaqus/Explicit
Joe ChuFinite Elements3 Views 5,174
Dear all, I recently conducted a numerical simulation of pilesoil interaction. I digged a vertical hole and insert the pile into that with general contact. Due to many reasons, I […] 
Modeling of granular soil by using UMAT hypoplastic in Abaqus
Naif AlsanabaniFinite Elements5 Views 4,241
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, […] 
hypoplasticity models in free Tochnog Professional
Silvia ImposimatoFinite Elements0 Views 3,857
In this post we will list builtin hypoplasticity formulations in the free Tochnog Professional program. Since these hypoplasticity options are builtin there is no need to link external umat routines. […]
Search SoilModels Website
Recent posts
 Connecting remote sensors to the Plaxis FEM 3D using python 26.7.2024
 Intergranular Strain Extensions 16.7.2024
 Abaqus Hypoplasticity Primary Variable Is Not Available 1.7.2024
 Vumat to calculate excess pore in Coupled EulerianLagrangian ABAQUS 1.7.2024
 Drucker Prager Cap Hardening Model 20.6.2024
 Soil model for sand under drained cyclic loading 15.6.2024
 In Abaqus how do we use *DLOAD for defining both userdefined pressure and body force simultaneously? 15.6.2024
 Issue in applying SANISAND04 UMAT to BoundaryValueProblem 31.5.2024
 Triaxial compression test using MohrCoulomb model in Abaqus 27.5.2024
 The Mises stress calculated by the hypoplasticity model considering the intergranular strain is all 0 20.5.2024
 Degradation ( disturbance ) of hypoplastic clay 23.4.2024
 Degradation ( disturbance ) of hypoplastic clay 23.4.2024
Recent Comments
 Johan Clausen on Triaxial compression test using MohrCoulomb model in Abaqus
 Ismail Khan on SANISAND for FLAC3D Download
 Tao Zhu on UMAT and oneAPI compiler issues. (SOLVED…i hope :p)
 Lei DAI on SANISAND for FLAC3D Download
 Yu Xiao on “Warning LNK4210” – hypoplastic clay UMAT
 Yu Xiao on “Warning LNK4210” – hypoplastic clay UMAT
 Charbel NOHRA on “Warning LNK4210” – hypoplastic clay UMAT
 Pouya Zahedi on UMAT for hypoplastic subroutine
 Mohsen Misaghian on SANISAND umat issue
 Ashesh Choudhury on Abaqus
 Hu Shihong on Triaxial compression test using MohrCoulomb model in Abaqus
 Mert Eyupgiller on Triaxial compression test using MohrCoulomb model in Abaqus
 Bala Harsha Srusti on ABAQUS job doesn’t make any increment
 PRATIK GOEL on Hardening parameters (yield stress & absolute plastic strain)
 Ashraf AHMED on ExCalibre – SoilModels Automatic Calibration
 Xy Ma on Sand and clay standard datasets
 Xy Ma on Sand and clay standard datasets
 Ramon Varghese on Hypoplastic Interface Model: Does anybody implemented this model in FLAC3D
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