ISSUES RECOMPILING CODE

[Samuel Alfred]

Hello Reza and members, 

I am trying to recompile PRISMS Plasticity to produce outputs for use in PRISMS Fatigue. 

First, I cloned the PRISMS Plasticity and then replaced  the calculateplasticity.cc, updateAfterincrement.cc and the userfunctions.cc from PRISMS Fatigue Github before compiling the code. 

I am trying to run a simple example on Tension on my local computer but it is not working. I have attached the error to this email. 

Thanks.

Samuel

 

[Samuel Alfred]

Okay, I discovered that the  calculateplasticity.cc, updateAfterincrement.cc and the userfunctions.cc from PRISMS Fatigue only work for rate dependent model. 

However, the outputs doesn't give me the Maximum compression and tension that I need to calculate the FIPs. Please, how do I go about this? 

[Mohammadreza Yaghoobi]

Hi Samuel,

My guess is the issue is with your input file.

If you use any of the input files we have in our manuscripts, they should generate the data you need.

You can try the input files for smaller microstructures used here:

https://www.nature.com/articles/s41524-021-00506-8

Link for Input files is presented in data availability.

Thanks

Reza

--
You received this message because you are subscribed to the Google Groups "Prisms-CPFE-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to prisms-cpfe-us...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/prisms-cpfe-users/e16a2dfd-7ed3-4cd6-ba1d-ebdfb30bbc57n%40googlegroups.com.

[Samuel Alfred]

Hello Reza. Thanks for your response. 

I figured out that I have been getting the quadrature outputs all along. I just thought that they will come with a csv file named 'Max_Tension' or 'Max_compression', as seen in the PRISMS fatigue tutorial video.

One more question please. I observed that most of the rate dependent examples and the simulation files used for the PRISMS-Fatigue manuscript have a 'User Defined Material Model Parameters' section enabled. The user manual is not clear on this. I am simulating a single phase FCC for Inconel 718 but I don't understand this lines of code (See below)  and how I can get them for my own IN718 simulation. 

# Number of User Material Constants in a Material model Phase 1
set Number of User Material Constants 1              = 13                                   

# Number of User Material State Variables in a Material model Phase 1
set Number of User Material State Variables 1        = 62

# Material Constants in a Material model Phase 1
set User Material Constants 1             = 1.0e-3, 0.0133333333333333333 , 1.0e-6 , 1.0e-6 ,  100 , 100,  2.0e6, 1.35e5, 2.0e4, 1421.0 , 70.0, 70.0

# Material State Variables in a Material model Phase 1
set User Material State Variables Initial Values 1        = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 , 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0

I just need to know what these numbers represent and how I can get them (i.e. literature or inverse optimization) because the code won't return if they are not defined. 

Thanks for your anticipated response.

Samuel 

[Mohammadreza Yaghoobi]

Please read the supplementary materials of this for the rate-dependent model:

https://www.nature.com/articles/s41524-021-00506-8

This explains the two-term backstress model.

set User Material Constants 1             = 1.0e-3    ( reference shearing rate), 0.04 (inverse strain rate sensitivity exponent), 1.0e-4 (Newton-raphson residual Threshold1), 1.0e-4 (Newton-raphson residual Threshold2),  20000 (Max number of Newton raphson iterations1) , 20000 (Max number of Newton raphson iterations2) ,     (The rest are kinematic hardening parameters; (h,r,m)1,2 Check the paper and try to match those parameters).

You can check the first part of calculatePlasticity,cc for ratedependent model and see what input represents starting from 0. which is this part of the code:

delgam_ref = UserMatConstants(0)*this->userInputs.delT ; // Reference slip increment
    strexp=UserMatConstants(1); // Strain rate sensitivity exponent ; the higher the less sensitive
    locres_tol=UserMatConstants(2); // Tolerance for innter residual of nonlinear constitutive model
    locres_tol2=UserMatConstants(3); // Tolerance for outer slip resistance loop
    nitr1=UserMatConstants(4); // Maximum number of iterations for inner loop convergence
    nitr2=UserMatConstants(5); // Maximum number of iterations for outer loop convergence
    ///////////////////////////////Backstress evolution parameters/////////
    h_back1=UserMatConstants(6);h_back2=UserMatConstants(7);
    r_back1=UserMatConstants(8);r_back2=UserMatConstants(9);
    m_back1=UserMatConstants(10);m_back2=UserMatConstants(11);

Reza

To view this discussion on the web visit https://groups.google.com/d/msgid/prisms-cpfe-users/572e1f65-3478-42ff-9af3-8f3f97dd91e9n%40googlegroups.com.