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 dsigma11-dsigma33=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 driver-hypoplastic models
Related Articles
-
Hardening soil Model UMAT
isma khabis||Constitutive Modelling|3 |Views 4,409
I have the hardening soil model in the UMAT format. However, when I run the code, it gives me this error as shown in the picture. I tried to make […] -
Getting started with UMAT for Hypoplastic clay model in ABAQUS
Kanika Lamba||Constitutive Modelling|0 |Views 4,378
Hello All, I am getting started with implementing UMAT for soil model in ABAQUS. I am planning to use Hypoplastic clay constitutive soil model to simulate the triaxial tests. I […] -
CEL analysis using hypo plasticity Vumat
abd Elfattah||Constitutive Modelling|8 |Views 4,842
Greetings i used interface Vumat with Umat for simulating CEL problem, after 50% of analysis, i got the error as “Excessive incremental rotation of the elements in element set ” […] -
New features in Excalibre
David Mašín||Constitutive Modelling|0 |Views 7,260
Hi all, We are happy to announce that the Excalibre team received support from the research grant GACR 22-12178S to further develop the calibration software, which is freely available at […] -
-
K0 determination Reloading using HP clay model
muhammad shakeel||Constitutive Modelling|0 |Views 3,928
Hi Prof. Masin , It is related to my previous post. Sorry to bother you again. I simulated Oedometer test in Plaxis as well as centrifuge reloading (using HP clay […] -
Spurious oscillations in calculated accelerations
Piotr Kowalczyk||Constitutive Modelling|2 |Views 4,153
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 […] -
ALERT doctoral school on constitutive modelling – hands on session triax input files
David Mašín||Constitutive Modelling|0 |Views 4,339
ALERT doctoral school on constitutive modelling is happening right now. The school lecture notes are available at the website of ALERT geomaterials (http://alertgeomaterials.eu/publications/) Below you find triax input files for […]
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 Creep-SCLAY model. 1.3.2024
- Cyclic tests with Triax element test driver 14.2.2024
- SUMMER SCHOOL ‘Numerical Modelling in Geotechnical Engineering’, Innsbruck – July 22nd-26th, 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
- 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
- Arie Koot on How to model the settlement in soil due to water drawdown.
- Konstantinos Chatzis on Cyclic tests with Triax element test driver
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, 304-312 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*depsMNR-Nep*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*depsMNR-Nep*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 Newton-Raphson iterations in Matlab driver. It is only calculated from residual stresses, which is insufficient for more complex mixed-controlled 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*depsilon-dy, see Eq. (12).
For you, that would mean either updating Matlab driver formulation or using more advanced element test driver, such as triax