Capillary model

Hi,

This indeed requires modifying the source code, so that this engine measures (for control purposes) the right stress state of the periodic sample.

In the current state, the PeriTriaxController engine only captures the contact force stress contribution, because of these two lines [1] and [2].

For the capillary force stress contribution to be captured as well, and the simulation to work as expected, you basically need to:
- replace the interaction physics "NormShearPhys" type in [1] with "CapillaryPhys"
- and add the capillary force in the interaction force in [2], with such line:
Vector3r f= (-1.)*( nsi->normalForce+nsi->shearForce ) - nsi->fCap;

I think everything should work afterwards, except you would then be unable to use PeriTriaxController engine for classical "dry" simulations (that do not include CapillaryPhys); these would crash.

Jerome

[1] https://github.com/yade/trunk/blob/e4e757f2e98a620e3177b7a36a1d10f69f6a6a28/pkg/dem/PeriIsoCompressor.cpp#L109
[2] https://github.com/yade/trunk/blob/e4e757f2e98a620e3177b7a36a1d10f69f6a6a28/pkg/dem/PeriIsoCompressor.cpp#L112

Hi,
What Jerôme suggests can be done in a script actually, i.e. adapting the deformation rate to the current stress (with user defined stress calculated in the script and including the capillary forces).
B

Thanks for reverting....
I am trying to simulate isotropic compression behaviour of unsaturated granular material. Can you suggest me which engine shall i use such that i am able to capture the yield limit for the packing with capillary bridges.

Thanks

This question was expired because it remained in the 'Open' state without activity for the last 15 days.