Welcome to the forums. Please post in English or French.

You are not logged in.

#1 2021-10-06 18:41:31

nicolo.grilli
Member
From: Bristol, United Kingdom
Registered: 2021-07-30
Posts: 10

Assign rotation matrix, link MFront to CodeAster, crystal plasticity

Dear MFront and Code Aster Users,

I have successfully linked MFront (Meric-Cailletaud crystal plasticity model)
h_ttp://tfel.sourceforge.net/MericCailletaudSingleCrystalPlasticity.html
with Code Aster (see attached .mfront and .comm file).

I have recorded the procedure in the following video,
showing a one element, single crystal simulation:
h_ttps://www.youtube.com/watch?v=evNhkEcntaU

The results I obtain agree with MTest
(see attached .test file).

Now the question is: how to change the rotation matrix?

In MTest I know I have to add for instance:

@RotationMatrix {{ 1.0 , 0.0 , 0.0 },
{ 0.0 , 1.0 , 0.0 },
{ 0.0 , 0.0 , 1.0 }};

Ideally, I would like to add additional material properties, like
@MaterialProperty Rij
representing the components of the initial rotation matrix and use them to
initialise the @RotationMatrix object.

This would allow me to pass easily those rotation matrix components
from the Code Aster .comm file to MFront.

Is this possible?

Or how do you envisage a strategy to assign the rotation matrix when running this model in Code Aster?

Next I would like to run polycrystal simulations, but once I figure out how to
pass the rotation matrix to Code Aster, it should be trivial.
I should just define different materials with different parameters in the Code Aster model,
representing the different grains.

Thank you very much in advance for your help.

I am really eager to learn more about the MFront code,
so I ask you about suggestions on how to look myself at the objects in the code,
so I can answer my questions by looking directly at the code.

For instance, if I want to look where the @RotationMatrix object is defined in the MFront code,
where should I look?

I tried to grep 'RotationMatrix' but I can only see that in .mtest files,
so I am a bit confused on where that keyword is defined in the code.
I am familiar with C++ but still a bit confused about the MFront code structure.

Thank you very much,
Best Regards,
Nicolò Grilli


Attachments:
RunCase_1_Stage_1_freebc.comm, Size: 3.17 KiB, Downloads: 46

Offline