Crystal Plasticity

[Cao DucToan]

Hi MOOSE team,

I am testing this example ( modules/tensor_mechanics/test/tests/crystal_plasticity/crysp_cutback.i ), but i cannot see the plasticity field (fp_zz] in Paraview.  Could you please help me for your experiences for this test? If you can estimate a time step  of the input file  for the crystalplasticity model, i really appreciate it.

Thank you very much for any help. 

Kind regards,

Toan D Cao

University of Alabama

Department of Mechanical Engineering

[Wen Jiang]

It is better to use UserObject-based Crystal Plasticity model ( https://mooseframework.inl.gov/source/materials/FiniteStrainUObasedCP.html) and its tests are under /modules/tensor_mechanics/test/tests/cp_user_object/

I assume you can see fp_zz in Paraview, but its field is a constant value, correct? This test has only one element, and we usually plot the the average of Fp for each element. So you would see a constant value for Fp in this case.

The time step depends on many factors, i.e. the loading, material models. etc,  so there is no general way that could give a best estimation. However, if you use a large time step and your solution does not converge, simulation will cut the time step and rerun the step. This is the main purpose of that test file.  In addition, you can use adaptive time stepper, see https://www.mooseframework.org/syntax/Executioner/TimeStepper/index.html 

[Cao DucToan]

Hi Wen Jiang

Thank you very much for your response. Yes, I see it in Paraview but the fp_zz is almost constant. Thanks again for your time. Do you have any examples to see the results and to understand a better program of Moose for beginners like me?

Kind regards,

Toan D Cao

University of Alabama

Department of Mechanical Engineering

--
You received this message because you are subscribed to the Google Groups "moose-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to moose-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/80aeef3e-f4f8-4ea2-9b3b-581849b2474bo%40googlegroups.com.

[Wen Jiang]

The front page of MOOSE website has a training video. The workshop slides can be found at https://www.mooseframework.org/workshop/#/. If you are interested in mechanics problems for instance, the documentation of Mechanics Module can be a good reference, see https://www.mooseframework.org/modules/tensor_mechanics/index.html       

To unsubscribe from this group and stop receiving emails from it, send an email to moose...@googlegroups.com.

[Cao DucToan]

Thanks Wen Jian for a good reference. I will check it and have a good day.

Kind regards,

Toan D Cao

University of Alabama

Department of Mechanical Engineering

Virus-free. www.avg.com

To unsubscribe from this group and stop receiving emails from it, send an email to moose-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/9f26e3cf-af5a-4d40-9a1f-8d4c6d29859ao%40googlegroups.com.

[Nicolò Grilli]

Dear Toan,

I have worked a couple of years on the crystal plasticity model in MOOSE,

I have developed some code that you can find in my github.

https://github.com/ngrilli/Marisol

and you can see the documentation in the papers I published between 2016 and 2019

in Prof. Marisol Koslowski group in Purdue University.

Of course the original reference for the crystal plasticity model

by Chockalingam, Tonks et al. is also a good start:

https://link.springer.com/article/10.1007/s00466-012-0741-7

Best Regards,

Nicolò Grilli

[Cao DucToan]

Dear Nicolò Grilli

I am really happy for your email. Thanks Nicolò Grilli. I appreciate your share.

Kind regards,

Toan D Cao

University of Alabama

Department of Mechanical Engineering

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/2071a804-f02c-467d-b32d-02271d70b392n%40googlegroups.com.

[Cao DucToan]

Dear Nicolò Grilli

I have a difficult problem with the crystalplasticity if I use an EBSD file of mesh. Could you please let me know how I can combine EBSD mesh and the crystalplasticity? I did a combination but unfortunately the program cannot work for the plasticity zone after elasticity zone. If you can help me with your advice, i really appreciate it.

I attached the input file here if you have a quick look for any help.

Kind regards,

Toan D Cao

University of Alabama

Department of Mechanical Engineering

 

On Sat, Oct 17, 2020 at 8:44 PM Nicolò Grilli <nicolog...@gmail.com> wrote:

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/2071a804-f02c-467d-b32d-02271d70b392n%40googlegroups.com.

[Wen Jiang]

EBSD reader only works for phase-field simulation at this time. The capability of reading EBSD, DREAD3D format to CP is still under development. 

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/2071a804-f02c-467d-b32d-02271d70b392n%40googlegroups.com.

[Cao DucToan]

Thanks Wen Jiang.

Kind regards,

Toan D Cao

University of Alabama

Department of Mechanical Engineering

 

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/7c005163-e66a-4641-9b5a-21d84a0a8850o%40googlegroups.com.

[Cao DucToan]

Dear Wen Jiang and Nicolò Grilli

Could you please help me how I can calculate the plastic energy of CP in Moose ? I really need your experience.

Kind regards,

Toan D Cao

University of Alabama

Department of Mechanical Engineering

On Mon, Oct 19, 2020 at 9:25 AM Wen Jiang <jiang...@gmail.com> wrote:

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/7c005163-e66a-4641-9b5a-21d84a0a8850o%40googlegroups.com.

[Wen Jiang]

We do not have any existing object to compute the plastic energy. Could you let us know the expression for your plastic energy? It should not be hard for you to write new codes for that. 

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/7c005163-e66a-4641-9b5a-21d84a0a8850o%40googlegroups.com.

[Cao DucToan]

Hi Wen Jiang

Thanks for your reply. I would like to get plastic work. The plastic energy = the plastic strain multiplies by elastic constant.

Kind regards,

Toan D Cao

University of Alabama

Department of Mechanical Engineering

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/5171656a-dec9-474d-8121-2867742e8be8o%40googlegroups.com.

[Nicolò Grilli]

Dear Toan,

I suggest you look at the EBSD2Abaqus software:

https://github.com/ngrilli/ebsd2abaqusEuler

It is based on the MTEX Matlab package.

You start from a .ctf EBSD file.

You can build a mesh in .inp format (Abaqus)

which is compatible with MOOSE.

I have used this program for Abaqus simulations,

so I am not sure it will work straight away for MOOSE, but it should.

Happy if you test it and let me know.

Best Regards,

Nicolò

[Cao DucToan]

Dear Nicolò

Thanks for your reply. Moose will work with EBSD file. Did you work to compute a stored plastic deformation energy of CP model in Mooose ? Because the stored plastic deformation energy is an important case for Alloy. 

https://www.sciencedirect.com/topics/engineering/plastic-deformation-energy

If you know how I can add a subroutine to compute the stored plastic deformation energy in Moose, I would really appreciate it.

Kind regards,

Toan

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/6dd5dd16-c131-4a86-a5f9-7e60d492d339n%40googlegroups.com.

[Wen Jiang]

You can follow "ElasticEnergyAux" to write a code to compute plastic energy. 

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/6dd5dd16-c131-4a86-a5f9-7e60d492d339n%40googlegroups.com.

[Cao DucToan]

Thanks Wen Jiang for this advice.

Kind regards,

Toan

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/b388bf5e-dc62-4aa2-b88d-6105b101b1abo%40googlegroups.com.

[Cao DucToan]

Hi Wen Jiang and Nicolò Grilli

Thanks for your help during time. My work on the crystal plasticity model using Moose is working well. Now i have my question for coupling between the crystal plasticity and the phase field fracture model. Do you have experience to share  your advice with this work ? How can I do it with the current version of Moose ? If Moose doesn't support it, how can I develop it ? I really appreciate your time.

Kind regards,

Toan D Cao

University of Alabama

Department of Mechanical Engineering

On Wed, Oct 21, 2020 at 9:54 AM Wen Jiang <jiang...@gmail.com> wrote:

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/b388bf5e-dc62-4aa2-b88d-6105b101b1abo%40googlegroups.com.

[Nicolò Grilli]

Dear Toan,

I am not sure about the current version of moose, but coupling CP with phase field damage is relatively easy.

I am attaching an example input file that was used a couple of years ago.

The objects used are part of the repository "Marisol" in my github:

https://github.com/ngrilli/Marisol

and you will see we made many different couplings and several models for phase field damage,

such as anisotropic crack propagation and fracture given by elastic and plastic energy.

Briefly, the phase field model wants as input the elastic strain energy,

or generally an energy that induces damage, which could be a combination of elastic and plastic.

Therefore, you may need to develop a kernel derived from the phase field fracture

that accepts energy quantities that are calculated in the crystal plasticity part

and declare them as material properties in the crystal plasticity part.

In the above repository you will find many examples of this,

but they correspond to a 2018 moose version, so be careful.

This paper:

https://aip.scitation.org/doi/abs/10.1063/1.5109761

and the other papers I published with Prof. Marisol Koslowski in Purdue University

contains the detailed equations of the coupling we did.

Best Regards,

Nicolò

[Cao DucToan]

Thank Nicolò Grilli for your information. That is very useful. I will be following your experience to develop a kernel for coupling between CP and PFF.

Kind regards,

Toan

--
You received this message because you are subscribed to the Google Groups "moose-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to moose-users...@googlegroups.com.

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/9b28531e-aaca-4171-ac6e-391e6984d95en%40googlegroups.com.

[Cao DucToan]

Dear Nicolò

Could you please share the input_slip_sys_HMX_austin.txt and the euler_ang_file_beamorientation0.txt ? I can test how your input file can work with the current moose.

Thank you very much to Nicolò for your help.

Kind regards,

Toan

--
You received this message because you are subscribed to the Google Groups "moose-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to moose-users...@googlegroups.com.

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/9b28531e-aaca-4171-ac6e-391e6984d95en%40googlegroups.com.

[Cao DucToan]

Dear Nicolò and Wen Jian

I got this error when I combined the CP model and the PFF model. Can you let me know what I can do for this error ?

*** ERROR ***
Material property 'L', requested by 'c_ACInterface' is not defined on block 0

Stack frames: 13
0: libMesh::print_trace(std::ostream&)
1: moose::internal::mooseErrorRaw(std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> >, std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> >)
2: callMooseErrorRaw(std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> >&, MooseApp*)
3: SubProblem::checkBlockMatProps()
4: FEProblemBase::checkProblemIntegrity()
5: Action::timedAct()
6: ActionWarehouse::executeActionsWithAction(std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > const&)
7: ActionWarehouse::executeAllActions()
8: MooseApp::runInputFile()
9: MooseApp::run()
10: ./combined-opt(+0x51ec) [0x55c539f311ec]
11: __libc_start_main
12: ./combined-opt(+0x54bf) [0x55c539f314bf]
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=1
:
system msg for write_line failure : Bad file descriptor

kind regards,

Toan

[Nicolò Grilli]

Dear Toan,

About the slip systems: you can just use the FCC example file that is in the moose tests.

When there is a missing material property, you should always check in the examples,

where the material property is defined, see for instance:

https://github.com/idaholab/moose/blob/next/modules/combined/test/tests/phase_field_fracture/crack2d_aniso.i

The missing material, called L, is:

  [./define_mobility]
    type = ParsedMaterial
    material_property_names = 'gc_prop visco'
    f_name = L
    function = '1.0/(gc_prop * visco)'
  [../]

The viscosity controls how fast the damage increases when the critical load is reached.

You can find the meaning of "visco", usually indicated by the eta greek letter in equation 13 of this:

https://aip.scitation.org/doi/abs/10.1063/1.5030656

Best Regards,

Nicolò Grilli

National University of Singapore

[Cao DucToan]

Dear Nicolò

Thanks for your help. I really appreciate it. Have a nice weekend.

Kind regards,

Toan D Cao

University of Alabama

Department of Mechanical Engineering

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/f35d929b-ad1b-4a7d-b680-8c94f7d77c16n%40googlegroups.com.

[Cao DucToan]

Dear Nicolò

I tried understanding the CP+PFF in your input file but ,unfortunately, I received this error. Do you have any ideas for help ?

*** ERROR ***
A 'PFFracBulkRatePlasticFactor' is not a registered object.

Kind regards,

Toan

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/f35d929b-ad1b-4a7d-b680-8c94f7d77c16n%40googlegroups.com.

[Andrew Slaughter]

Check spelling and make sure the phase field module is enabled in your Makefile. Also, be sure to recompile after changing the Makefile.

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/CAGtVYWK3tgg5O-5GHJwin8aOZOPa52ucScWqM6UvaZpufAFS5g%40mail.gmail.com.

[Cao DucToan]

Dear Nicolò

I received this error. Do you have any ideas for help ? thank you very much in advance.

/projects/moose/modules/tensor_mechanics/src/materials/PFFracBulkRateMaterial.C:26:45: error: invalid conversion from 'const Function*' to 'Function*' [-fpermissive]
     _function_prop(isParamValid("function") ? &getFunction("function") : NULL)

The code is written as follows,

InputParameters validParams<PFFracBulkRateMaterial>()
{
  InputParameters params = validParams<Material>();
  params.addClassDescription("Material properties used in phase-field fracture damage evolution kernel");
  params.addParam<FunctionName>("function", "Function describing energy release rate type parameter distribution");
  params.addParam<Real>("gc", 1.0, "Energy release rate type parameter");

  return params;
}

PFFracBulkRateMaterial::PFFracBulkRateMaterial(const InputParameters & parameters) :
    Material(parameters),
    _gc(getParam<Real>("gc")),
    _gc_prop(declareProperty<Real>("gc_prop")),
    _function_prop(isParamValid("function")  &getFunction("function") : NULL)
{
}

Kind regards,

Toan

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/CA%2Bd67H-Z_O8oRZEs2Bj5H%2Bq1nPvdXNXxDUYkx6Q8CBhLJdzvhA%40mail.gmail.com.

[Nicolò Grilli]

Dear Toan,

not fully sure, but that may not be the way Moose currently handle functions.

I suggest you copy from ParsedMaterialBase:

https://github.com/idaholab/moose/blob/next/framework/src/materials/ParsedMaterialBase.C

In the constructor: .C file:

  // get function expression
  _function = parameters.get<std::string>("function");

  // Function expression
  params.addRequiredParam<std::string>("function",
                                       "FParser function expression for the phase free energy");

and in .h file:

  /// function expression
  std::string _function;

Try these expressions.

Also, you may consider defining your function inside a material object,

I mean write a material object that calculates your quantity and assign it to a MaterialProperty

that you can then use in the kernel.

You write functions of the coordinates by accessing the coordinates as:

// Get coordinates of the current IP
x = _q_point[_qp](0);
y = _q_point[_qp](1); z = _q_point[_qp](2);

See for example this one:

https://github.com/ngrilli/c_pfor_am/blob/main/src/auxkernels/LaserScanAux.C

Best Regards,

Nicolò Grilli

National University of Singapore

[Cao DucToan]

Dear Nicolò

Thank you very much for your help. I will try it .

Kind regards ,

Toan

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/122900bf-7521-4bd6-b4a9-3f0ab157ddf3n%40googlegroups.com.

[Cao DucToan]

Dear Nicolò

I did several steps following your suggestion. Unfortunately , it doesn't work.  I don't know what happened with this.

Thanks Nicolò. If you have any idea, please let me know.

Kind regards,

Toan

To view this discussion on the web visit https://groups.google.com/d/msgid/moose-users/122900bf-7521-4bd6-b4a9-3f0ab157ddf3n%40googlegroups.com.