TarletonGroup / CrystalPlasticity

CP UMAT and CZM UEL for Abaqus
152 stars 72 forks source link

question about FS parameter #74

Closed incredible58 closed 2 months ago

incredible58 commented 2 months ago

Dear sir,

I'm using the subroutine to calculate the FS parameter. I've already change the flag value of "statev_outputs(20) " from 0 to 1 ,and that of "statev_outputs(10) " from 1 to 0. But I've got error in abaqus "*** ABAQUS/standard rank 0 encountered a SEGMENTATION FAULT".

I've done the simulation to calculate the total slip per system with the same element type, and materials parameters. What else should I notice to calculate FS parameter?

I would be appreciated if I could get rely from you.

Thank you very much.

Sincerely, David

EralpDemir commented 2 months ago

Hi, One problem might be the number of dependent state variables (DEPVAR) setting in Abaqus. You should increase it to 13 or greater to solve this issue. Best

incredible58 commented 2 months ago

Dear sir,

I'm very glad to hear from you.

I've change DEPVAR to 100, but it doesn't work.

When I calculate FS, is it correct that I only set statev_output(20) to 1 and the rest to 0? My model is a simple single grain with 10000 CPE6 elements. The phase is hcp with 12 slip systems.

Is there any other suggestions I may use?

Thank you very much.

Sincerely, David

EralpDemir commented 2 months ago

You may try to set maxnslip to 30 for hcp.

Please also check your log or dat file. It must write your element number and integration points per element. The integration point shall be 3 for that element with full integration. Also make sure it is not a reduced integration element.

Best

incredible58 commented 2 months ago

Dear sir

I've used the same set(maxnslip : 30, integration point : 3 ) to calculate total slip per system , and the results are perfect with the subroutine. However, the same setting generates an error when calculating FS.

Thank you very much.

Sincerely, David

EralpDemir commented 2 months ago

Hi again,

I have made a minor change in the FS calculation. Please download the latest version in Github.

If it does not work, let me know. I will need to look in detail which may take some time.

Best

incredible58 commented 2 months ago

Dear sir

Thank you for the latest version. I'll try it. Whether it works, i'll tell you. Thank you very much.

Sincerely, David

incredible58 commented 2 months ago

Dear sir

I'm sorry to tell you that it still doesn't work.

And may I ask for references about "plastic dissipation power density" used in the subroutine?

Thank you very much.

Sincerely, David

EralpDemir commented 2 months ago

Ok, I have added a few more checks in case you are hitting into a problem.

Please try the latest version in the Github (v2.23).

I am also sending the input file for testing I used, it works without any issue. Job-1_INPFILE.txt

If you still have an issue, I do not know the problem you are hitting into.

Hope this helps, Eralp

incredible58 commented 2 months ago

Dear sir

Thank you so much for your great patience and huge help. Now I can calculate FS parameter in my model.It's very helpful to my research.

Thanks again, we couldn’t have pulled this off without you

Sincerely, David