Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

How to couple pore water pressure to simulate triaxial test using Cam-clay model

Please login with a confirmed email address before reporting spam

Here are some questions I want to ask. Thanks!

the figures are in doc file.

Problem description:
We want to simulate the behavior of soil sample from a trialxial test. The cam-clay model here is applied to simulate the soil material. Besides, the pore water pressure is a concerning issue in this simulation. So, we need to couple the pore water pressure in our model.
The simulation is based on an existing example named “Isotropic Compression” in models library.
In order to consider the pore water pressure, we add “Darcy’s law” physics into the model as shown in fig. 1.

In Cam-clay material’s page, the pore water pressure from physics of darcy’s law is set as shown in fig.2 .

From the cam-clay’s “equation view”, we find that the stress contains the term of pore water pressure. Just as the following:
solid.Sl11=-(solid.Sil11+solid.Sil22+solid.Sil33)/3-solid.pf-solid.Pc0*exp(-((1+solid.evoid0)*(solid.eel11+solid.eel22+solid.eel33)+solid.evoid0-solid.Nvoid+solid.lambdaComp*log(solid.Pc0/solid.prefVCL))/solid.kappaSwelling)+2*solid.G*(solid.eel11-(solid.eel11+solid.eel22+solid.eel33)/3)+solid.Sil11 ,
where solid.pf denotes the pore water pressure.

In darcy’s law physics, a “weak contribution” for the domain is added.
-d(solid.evol,TIME)*rhof*test(p)
As fig3. showing. The pore water pressure is coupled into the model with this step.

Parameters setting:
Showed in Fig.4, the initial stress is modified as 400[kPa].
The calculated initial void ration e0 is 0.6556. As fig5 showing:
The above settings are to make sure the state of stress in the yield surface.
The boundary conditions are shown in fig6:
The simulation with the above settings cannot run successfully. And the errors will be tipped as the following.
Failed to find consistent initial values.
Failed to compute elastoplastic strain variables.
Last time step is not converged.
If setting the top surface of the soil sample drained, i.e. p=0, another error is prompted as following:
Failed to find consistent initial values.
Error in elPlastic element.
Lambda non positive for plastic point
Last time step is not converged.

Our questions are:
(1) Is the above procedure to couple the pore water pressure with Cam-clay model to simulate triaxial test correct?
(2)What the meanings of the above error information are?
(3) Is there any similarly example to provide to us?


2 Replies Last Post 24 mag 2013, 10:29 GMT-4
COMSOL Moderator

Hello Liu Ganbin

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.


Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 14 mag 2013, 02:27 GMT-4
Hi,

First question is yes I would say. The prefered way by Comsol for the Cam-Clay model is to couple the Pore-pressure there as you done.

But I´m not sure why you add a weak expression in Darcy´s law.? Why dont you use the mass source in the node instead!

I´m not used to Cam-Clay in Comsol but
Why you get the error message I cant say striaght up if I dont have the mph-file. Sometimes it is due to solver settings e.g. dependent variable could be change from automatic to e.g. 1 and the solver if you use fully coupled change to appropriate settings etc.

I hope this helps. Otherwise post the mph-file so we could have a look at it.

Regards,

Mats
Hi, First question is yes I would say. The prefered way by Comsol for the Cam-Clay model is to couple the Pore-pressure there as you done. But I´m not sure why you add a weak expression in Darcy´s law.? Why dont you use the mass source in the node instead! I´m not used to Cam-Clay in Comsol but Why you get the error message I cant say striaght up if I dont have the mph-file. Sometimes it is due to solver settings e.g. dependent variable could be change from automatic to e.g. 1 and the solver if you use fully coupled change to appropriate settings etc. I hope this helps. Otherwise post the mph-file so we could have a look at it. Regards, Mats

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 24 mag 2013, 10:29 GMT-4
Isn't adding the Darcy's law there enough to guarantee the coupling between the deformations of the soil to the water pressure?

I do not understand why that weak contribution is added. You might want to remove the weak contribution and run it and see if it runs and if yes, if the results make sense.

Please leave a note here if you got any improvements.

I am personally interested in modelling dilation of submerged sand and capturing negative pore pressures in it. Do you think that the Cam clay can be modified to represent fully saturated sands in dilation?

Cheers,
Hossein
Isn't adding the Darcy's law there enough to guarantee the coupling between the deformations of the soil to the water pressure? I do not understand why that weak contribution is added. You might want to remove the weak contribution and run it and see if it runs and if yes, if the results make sense. Please leave a note here if you got any improvements. I am personally interested in modelling dilation of submerged sand and capturing negative pore pressures in it. Do you think that the Cam clay can be modified to represent fully saturated sands in dilation? Cheers, Hossein

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.