-
Hi, This question is first posted in FEniCS discourse and it was suggested that I post it here. After the day of posting, I sort of found a solution, but I would still like to bring this to your attention. Thanks! -- I am interested in doing simulations in FEniCS with Mohr-Coulomb criterion. I followed this tutorial 1 which implements a von Mises criterion with MFront. Here, I successfully compiled the MFront file and ran the sample problem. Now, I moving on to Mohr-Coulomb criterion with this tutorial 2. There are two code blocks that I saved into MFront files. The first one was in the section "Simplification of the MFront file : use of TFEL/Material/MohrCoulombYieldCriterion.hxx file" and other one is in the section "Simplification of the MFront file : use of the StandardElastoViscoPlasticity brick". Both return this error when compiling:
Anyone who had the same problem? Any solution? Thanks! EDIT: I stitched the codes in the “Implementation” section of tutorial 2 to create the most explicit Mfront file without the simplifications in the next sections. I compiled it successfully but I have to include “@parameter pi = 3.14159265359;” and “constexpr auto local_zero_tolerance = real(1e-14);” as they were missing. |
Beta Was this translation helpful? Give feedback.
Replies: 6 comments 12 replies
-
Hi @gitable05, Thanks for you interest to Would you attach you implementations ? Best, BTW, did you used mgis.fenics (https://thelfer.github.io/mgis/web/mgis_fenics.html) |
Beta Was this translation helpful? Give feedback.
-
Please find attached three implementations. They all compile with For All in all, those implementations are not robust as the normal may strongly oscillate between two iterations, so it need to be slightly modified using a small unpublished trick. Just send me an email at tfel-contact@cea.fr, and I'll send you a proper implementation next week. |
Beta Was this translation helpful? Give feedback.
-
@thelfer Dear Thomas,
Best regards and Thanks for your valuable time. Sukumar |
Beta Was this translation helpful? Give feedback.
-
Is the singularity container is an issue? Regards |
Beta Was this translation helpful? Give feedback.
-
Thank you very much Thomas. In fact I am using the container with code-aster 15.4 MPI version available at Would be waiting to hear from you. Sukumar |
Beta Was this translation helpful? Give feedback.
-
Dear Thomas, Thank you very much for your time and advice.
With warm regards and Best wishes Sukumar |
Beta Was this translation helpful? Give feedback.
@gitable05.
Please find attached three implementations. They all compile with
TFEL-3.4.x
although the third one may fail ifC++-17
support is enabled (See Issue #80).For
TFEL-4.x
versions, it seems that PR #5 messed things up and broke theMohrCoulombStressCriterion.ixx
file. We will treat this in Issue #81. Sorry for the inconvenience.All in all, those implementations are not robust as the normal may strongly oscillate between two iterations, so it need to be slightly modified using a small unpublished trick. Just send me an email at tfel-contact@cea.fr, and I'll send you a proper implementation next week.
MohrCoulomb.zip