Question regarding fluid-soil coupling calculations

[Wangjiongchao]

Honored Professor Shi,
I've recently been learning to use the MPM section in Geotaichi. I've successfully run the DPmaterial2D.py and NewtonianFluid2D.py examples, but I meet issues when attempting fluid-soil coupling calculations. After setting the material_type as TwoPhaseDoubleLayer, I further configured two materials: a Newtonian fluid and a Drucker-Prager soil. However, the simulation failed with the following error: 【File “D:\Geotaichi\GeoTaichi\src\mpm\MaterialManager.py”, line 120. ValueError: Constitutive Model: Newtonian error! Only the following is available: [‘LinearElastic’, ‘HenckyElastic’, ‘NeoHookean’, ‘ElasticPerfectlyPlastic’, ‘MohrCoulomb’, ‘DruckerPrager’, ‘ModifiedCamClay’, ‘CohesiveModifiedCamClay’, ‘UserDefined’].】
I would like to ask how this calculation can be implemented? Or what modifications do I need to make to the current version of the code?
I look forward to your reply!

[wme7]

I had a similar issue
Probably you have to also add a Newtonian material on the TwoPhaseDoubleLayer definition at the MaterialManager

[Wangjiongchao]

Yes, I tried this way and encountered various issues. The first problem was on line 211 of scencemanger.py, where solver_type == ‘’Explicit” does not support TwoPhaseDoubleLayer. However, the root cause seems to point to the enginekernel.py file. My understanding is that the code for the fluid-soil coupling effect might be missing.